Upstream PID simulation: Simulation Results II

MICE PID Review 12 October 2006 1

Upstream PID simulation: Simulation Results II

Hideyuki Sakamoto (Osaka Univ.)

MICE Collaboration Meeting

12th October 2006

MICE PID Review 12 October 2006 2

Contents

• Setup for G4MICE

• Setup for CKOV1– Material– Geometry

• Summary

MICE PID Review 12 October 2006 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

MICE PID Review 12 October 2006 4

CKOV Geometry

Based on the geometry file from Ghislain.

Object per detection unit Shape X,Y [mm] Z [mm] DZ [mm] Material

Entrance window BOX 450.0 2.5 5.0 Polyacetal

Mirror BOX 450.0 6.5 3.0 B270

Aerogel radiator BOX 450.0 20.5 23.0 Aerogel

Glass window BOX 600.0 33.5 3.0 B270

Conical mirror “Conical” Variable Variable Varible Polycarbonate

Exit window BOX 450.0 344.0 5.0 Polyacetal

MICE PID Review 12 October 2006 5

Improvements for G4MICE

• 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

Information from Kenny’s email

MICE PID Review 12 October 2006 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 );

MICE PID Review 12 October 2006 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) se

ems to be easy to build in.

MICE PID Review 12 October 2006 8

Conical mirror to be included

From Ghislain’s 3D CAD file

Beam

Here PMT is located.

From L. Cremaldi From L. Cremaldi

MICE PID Review 12 October 2006 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

Boxes for aerogel or windows are defined bellow

Objects to make Ckov mirror, from which CAD file to be readed.

MICE PID Review 12 October 2006 10

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 ); }

}

Class for CkovMirror

MICE PID Review 12 October 2006 11

Reading CAD file

• Geant4 can make the geometry from CAD file (STEP AP203)– Need to publish in format of STEP with some t

hird-party software, cannot be read with raw CAD format (dwg,dxf)

– STEP/STEPinterface in Geant4 will provide interfaces to read STEP files.

MICE PID Review 12 October 2006 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 ?

MICE PID Review 12 October 2006 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.

MICE PID Review 12 October 2006 14

Beam

Detection unit#1 Detection unit#2

Conical mirror (but box…)

Exit window

Glass window

Radiator (box)

MirrorEntrance window

Side View

MICE PID Review 12 October 2006 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 !!