270 likes | 446 Views
NICADD Test Beam Simulations. Project Geometry Standalone TBMokka Comparisons Plans. Project. NICADD collaboration with DESY/Calice Test Beam Prototype: ECal, HCal, Tail Catcher ECal/HCal similar to SDJan03 layering Sensitive Materials: Polystyrene, Silicon
E N D
NICADD Test Beam Simulations Project Geometry Standalone TBMokka Comparisons Plans Test Beam Simulations
Project • NICADD collaboration with DESY/Calice • Test Beam Prototype: ECal, HCal, Tail Catcher • ECal/HCal similar to SDJan03 layering • Sensitive Materials: Polystyrene, Silicon • NICADD concentration on Tail Catcher • TB/Mokka, TB standalone • Other projects at DESY (Geant4/Geant3) • concept generalizable to any “box-like” calorimeter Test Beam Simulations
Geometry Overview Component No. Layers * thickness = total thickness X x Y dim. (m) Materials HCal 50 * 2.5 cm = 125 cm 1 x 1 Stainless Steel, Polystyrene ECal 30 * 5 mm = 15 cm 1 x 1 W, G10, Si, Cu, Air Tail Catcher 6 * 10.5 cm = 63 cm 1.5 x 1.5 Same as HCal Total 203 cm max 1.5 x 1.5 - Test Beam Simulations
ECal Front Side 1 m ~6 cm Test Beam Simulations
HCal Front Side (Angled) 1.25 m 1 m Test Beam Simulations
Tail Catcher Side Front (Hor) Front (transparent) Front (Vert) 63 cm 1.5 m 1.5 m Test Beam Simulations
Alternate TC Configurations 12 layers 10-layer Fine Test Beam Simulations
Cells Component Layers Cells / Layer Total Cells HCal 50 100 inner + 364 outer = 464 23,200 ECal 30 10,000 300,000 Tail Catcher 6 15 90 Total 86 - 323,290 1.5 m Test Beam Simulations
Test Beam Standalone • G4 project • g4phys library (Mokka v01-05), popt • GPS • placement for each cell • custom geometry map • text output format • available via CVS Test Beam Simulations
TB Hits Processing TBSensitiveDetector::ProcessHits() PhysVolInfo *pinfo = PhysVolMap->GetPhysVolInfo(physVol); TBCellHit *theHit = CellHit(physVol); if (!theHit) { theHit = new TBCellHit(); if (pinfo) { theHit->SetPhysVolInfo(pinfo); mapHits.insert(HitMapPair(physVol,theHit)); PhysVolInfo::PhysVolType vtype = pinfo->GetVolType(); if (vtype == PhysVolInfo::CELL) { theHit->IncrEdep(edep); TBTrackInformation *tinfo = (TBTrackInformation*)(aStep->GetTrack()->GetUserInformation()); theHit->IncrTrkEdep(tinfo->GetTrackID(),edep); } else if (vtype == PhysVolInfo::ABSORBER) { theHit->IncrAbsEdep(edep); } } Test Beam Simulations
Text Output Format EVENT evt_no TRACK n_tracks PDGcode x y z px py pz ECAL n_hits layer_n … x y z cell_edep abs_edep n_tracks track_n track_edep … HCAL … … TAILCATCHER … … Test Beam Simulations
Event Display Test Beam Simulations
TBMokka • Mokka v02-01 • Geometry Driver • Sensitive Detector • Cell Replicas • X/Z cell dimensions • Indexing scheme started • LCIO output • “proof of concept” for Mokka driver authoring • includes additions to Mokka code, such as working indexing scheme in Geometry/CGA/CellHit.cc Test Beam Simulations
Mokka CGA Model TBecal00 Databases TBhcal00 hcal layer_inner layer_outer layer_thickness TBecal00 ecal layer layer_thickness TBcatcher00 catcher layer layer_thickness Drivers TBhcal00 TBecal00 TBcatcher00 Test Beam Simulations
Histograms • initial TB/TBMokka comparisons • all samples 20 GeV pions unless otherwise noted • performed with JAS3: AIDA, LCIO reader, txt FileMgr • ... Test Beam Simulations
TB: HCal Live Edep Test Beam Simulations
TB: HCal Live Edep Test Beam Simulations
TB: % Edep in Inner Cells Test Beam Simulations
TBMokka: Component Hit Number Test Beam Simulations
Number of Hits per Event Test Beam Simulations
Edep by Event 200 piplus @ 20 GeV (all comparison histos) Test Beam Simulations
Hit Edep Test Beam Simulations
ECal: Number of Hits per Event Test Beam Simulations
ECal: Number of Hits by Layer Test Beam Simulations
HCal: Number of Hits Test Beam Simulations
Plans • improve TBMokka performance • TB runtime geometry specification • find source of discrepancies in # hits/edep and eliminate • common physics list • introduce TBMokka to European sim. teams (Calice/DESY) • TBMokka into Mokka CVS Test Beam Simulations
Additional Information • Homepage: http://nicadd.niu.edu/~jeremy/g4proj • Source Browser: http://nicadd.niu.edu/~jeremy/tbeam/html/ • CVS (bash): • export CVSROOT=:ext:lcdguest@nicadd.niu.edu:/cvs • export CVS_RSH=ssh • cvs co tbeam • password: lcd_2004 • g4phys: http://nicadd.niu.edu/~jeremy/download/g4phys.tgz • popt: http://www.gnu.org/directory/libs/popt.html Test Beam Simulations