1 / 15

Upstream PID simulation: Simulation Results II

Upstream PID simulation: Simulation Results II. Hideyuki Sakamoto (Osaka Univ.) MICE Collaboration Meeting 12 th October 2006. Contents. Setup for G4MICE Setup for CKOV1 Material Geometry Summary. Setup G4MICE. Setup the latest version into the OSAKA local machine. Compiler gcc3.2.3

gudrun
Download Presentation

Upstream PID simulation: Simulation Results II

An Image/Link below is provided (as is) to download presentation Download Policy: Content on the Website is provided to you AS IS for your information and personal use and may not be sold / licensed / shared on other websites without getting consent from its author. Content is provided to you AS IS for your information and personal use only. Download presentation by click this link. While downloading, if for some reason you are not able to download a presentation, the publisher may have deleted the file from their server. During download, if you can't get a presentation, the file might be deleted by the publisher.

E N D

Presentation Transcript

  1. Upstream PID simulation: Simulation Results II Hideyuki Sakamoto (Osaka Univ.) MICE Collaboration Meeting 12th October 2006

  2. Contents • Setup for G4MICE • Setup for CKOV1 • Material • Geometry • Summary

  3. Setup G4MICE • Setup the latest version into the OSAKA local machine. • Compiler gcc3.2.3 • GEANT4 version 4.8.1 • CLHEP version 1.9.2.2 • CERNLIB version 2003 • ROOT version 5.10.0 • GSL version 1.8

  4. CKOV Geometry Based on the geometry file from Ghislain.

  5. Improvements for G4MICE Information from Kenny’s email • Adding Materials • aerogel1  :Density=0.251 (g/cm3) Si:O = 0.467 : 0.533 • aerogel2  :Density=0.430 (g/cm3) Si:O = 0.467 : 0.533 • polyacetal :Density=1.41 (g/cm3) C:H:O=0.375 : 0.125 : 0.5 • B270 :Density=2.44 (g/cm3) Si:O:B:Na:K:Ca = 0.319 :0.486 : 0.033 : 0.057 :0.099 :0.007 • Adding Geometry for Mirror • Conical shapes: will be imported from CAD file

  6. Materials in G4MICE • Materials are easy to be added into G4MICE • Defined ${MICE_PATH}/MICE/Simulation/src/FillMaterials.cc e.g. aerogel1 ****************************************************** // Aerogel for CKOV1 in Vessel1 G4double density = 0.251*g/cm3; std::string name = "AEROGEL_1"; G4Material* aero1 = new G4Material(name,density, nComp=2, kStateSolid); aero1->AddElement(elSi, 46.7*perCent); aero1->AddElement(elO, 53.3*perCent); run.miceMaterials->addMaterial( aero1, name );

  7. Geometry in G4MICE • Easy to include geometries with simple or pre-installed in Geant4 shapes(e.g. BOX,TUBE..) • BOX : for Entrace/Exit Window, Aerogel radiator • Not so easy to include complex shapes. • “Conical”: for Mirror • Installing by reading from a CAD file(STEPAP203) seems to be easy to build in.

  8. Beam Here PMT is located. Conical mirror to be included From L. Cremaldi From Ghislain’s 3D CAD file

  9. Source of geometry in G4MICE • ${MICE_PATH}/MICE/Simulation/src/MICEDetectorConstruction.cc ***************************************************************************** else if ( detector == "EmCalKLFiber" ) { EmCalKLFiber* fiber = new EmCalKLFiber( mod, mat, moth); logic = fiber->logicalFiber(); place = fiber->placementFiber(); } else if ( detector == "CkovMirror" ) { CkovMirror* mirror = new CkovMirror( mod, mat, moth); logic = mirror->logicalMirror(); place = mirror->placementMirror(); } Code for EmCal Objects to make Ckov mirror, from which CAD file to be readed. Boxes for aerogel or windows are defined bellow

  10. Class for CkovMirror CkovMirror::CkovMirror( MiceModule* mod, G4Material* mater, G4VPhysicalVolume *mlv ) { G4AssemblyCreator MyAC("Ckov_mirror.stp"); // Associate a creator to a given STEP file. MyAC.ReadStepFile(); // Reads the STEP file. STEPentity* ent=0; // No predefined STEP entity in this example. A dummy pointer is used. MyAC.CreateG4Geometry(*ent); // Generates GEANT4 geometry objects. void *pl = MyAC.GetCreatedObject(); // Retrieve vector of placed entities. G4Assembly* assembly = new G4Assembly(); // An assembly is an aggregation of placed entities. assembly->SetPlacedVector(*(G4PlacedVector*)pl); // Initialise the assembly. G4int solids = assembly->GetNumberOfSolids(); // Get the total number of solids among all entities. for(G4int c=0; c<solids; c++) // Generate logical volumes and placements for each defined solid. { G4PlacedSolid* ps = assembly->GetPlacedSolid(c); G4LogicalVolume* logicMirror = new G4LogicalVolume(ps->GetSolid(), mater, "STEPlog"); G4RotationMatrix* hr = ps->GetRotation(); G4ThreeVector* tr = ps->GetTranslation(); // G4VPhysicalVolume* pv = new G4PVPlacement(hr, *tr, ps->GetSolid()->GetName(), lv, // experimentalHall_phys, false, c); G4VPhysicalVolume* physiMirror = new G4PVPlacement( hr, *tr, logicMirror, ps->GetSolid()->GetName(), mlv->GetLogicalVolume(), false, c ); } }

  11. Reading CAD file • Geant4 can make the geometry from CAD file (STEP AP203) • Need to publish in format of STEP with some third-party software, cannot be read with raw CAD format (dwg,dxf) • STEP/STEPinterface in Geant4 will provide interfaces to read STEP files.

  12. !!..BUT, Unfortunately, those interfaces were dropped off from Geant4 when version4.6 was released. Not supported any more…>< • Is there any interface to read CAD file? • Anyway, I tried using Geant4.5.p02 source code for STEP/STEPinterface, but not success to compile so far. • [My Question]: Is this no problem to use such dropped off codes even if compile will succeed ?

  13. Present Status • Material: AIR for all • Geometry: BOX for all • Conical mirror is replaced with BOX shapes. • I did not commit the source code yet, but will be done soon. • So far, the only problem is geometry for Conical mirror.

  14. Side View Detection unit#1 Detection unit#2 Entrance window Mirror Radiator (box) Glass window Beam Exit window Conical mirror (but box…)

  15. Summary • Setup for CKOV1 was almost done except for Conical shapes for mirror. • There is a problem to read CAD file for building the conical mirror. • If the above problems will be solved, setup for CKOV1 can be finished => GOTO Simulation !!

More Related