Download
nicadd test beam simulator n.
Skip this Video
Loading SlideShow in 5 Seconds..
NICADD Test Beam Simulator PowerPoint Presentation
Download Presentation
NICADD Test Beam Simulator

NICADD Test Beam Simulator

122 Views Download Presentation
Download Presentation

NICADD Test Beam Simulator

- - - - - - - - - - - - - - - - - - - - - - - - - - - E N D - - - - - - - - - - - - - - - - - - - - - - - - - - -
Presentation Transcript

  1. NICADD Test Beam Simulator Project Description Classes Geometry Hits Processing Command Interface Particle Generation Text Output Format Histograms Visualization Additional Information Jeremy I. McCormick NICADD Test Beam Simulator

  2. Project Description and Goals • Test Beam Prototype: HCal, ECal, Tail Catcher • Sensitive Material: Polystyrene, Silicon • HCal & ECal layering similar to SDJan03 • Single volume per cell geometry • Edep by cell with track contributions • Investigate E-flow & cluster separability NICADD Test Beam Simulator

  3. Software Libraries • Geant4.5.2.p01 • GPS • g4phys SO (Mokka v01-04 Physics List) • popt for CL • Composite Materials def from PEGS4 (CERN) NICADD Test Beam Simulator

  4. Classes User • TBEventAction • TBPrimaryGeneratorAction • TBRunAction • TBSteppingAction • TBTracking Action • TBVisManager G4 Utility • Cellulizer • BoxCellulizer • MaterialMap Geometry • GeometryMap • PhysVolInfo • TBDetectorConstruction • TBCal • TBTailCatcher • TBECal • TBHCal Hits and Tracks • TBCellHit • TBSensitiveDetector • TBTrackInformation • TBEdepMap NICADD Test Beam Simulator

  5. Class Diagram Association Directed Assoc. Inheritance NICADD Test Beam Simulator

  6. GeometryMap GeometryMap ComputeWorldVectors() AddDaughters(vol, layno) TBDetectorConstruction TBSensitiveDetector PhysVolInfo GetLayNo() GetVolType() GetGlobId() GetHit() GetThreeVector() ComputeWorldVector() TBCellHit NICADD Test Beam Simulator

  7. Detector Geometry TC ECal HCal 1.5 m 1 m ~2 m NICADD Test Beam Simulator

  8. 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 - NICADD Test Beam Simulator

  9. 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 NICADD Test Beam Simulator

  10. HCal Front Side (Angled) 30 cm 1.25 m 1 m NICADD Test Beam Simulator

  11. ECal Front Side (partial) 1 m ~6 cm NICADD Test Beam Simulator

  12. Tail Catcher Side Front (Horizontal Layer) Front (transparent) Front (Vertical Layer) 63 cm 1.5 m 1.5 m NICADD Test Beam Simulator

  13. Command Interface • Command-line Usage: -? : print usage -c : visualize cells -t : visualize tracks -h : visualize hits -o outfile_path : output filepath -m macro_file : macro -v : start with vis.mac • All Options Example: testbeam -vhtc -o outfile_txt -m run.mac NICADD Test Beam Simulator

  14. Particle Generation General Particle Source • 1 to n MC particles • fixed or random E, pos & mom • theta, phi within min, max • energy spectra • /tbeam/doRun no_events for multiple events NICADD Test Beam Simulator

  15. GPS Conical Example 10 pi+ pos = (0, 0 ,1990) theta = 0-10 phi = all E = 2-10 GeV charged traj only NICADD Test Beam Simulator

  16. GPS Square Example 10 pi+ pos = (rX,rY,1990) square = 15 cm / side E = 10-20 GeV charged traj only NICADD Test Beam Simulator

  17. 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); } } NICADD Test Beam Simulator

  18. End of Event TBSensitiveDetector::EndOfEvent() { Copy Hits to HC } TBEventAction::EndOfEventAction() WriteEvent(evt) WriteHitsCollection(HC) aHit->Write(ostr) TrackEdep->Write(ostr) aHit->Draw() VisualizeTrajectories(evt) NICADD Test Beam Simulator

  19. Text Output Format EVENT evt_no ECAL no_hits layer_no … x y z cell_edep abs_edep track_no track_edep … HCAL … … TAILCATCHER … … NICADD Test Beam Simulator

  20. Output Example EVENT 0 HCAL 2 1 0 0 0 0 -352.5 0 407.334 0 1 0 0 65 -85 435 1.22144 0 1 2 1.22144 ECAL 1 1 0 0 85 95 584.3 0.113163 0 1 2 0.113163 TAILCATCHER 1 1 0 0 0 300 -690 0.00100138 0 1 1 0.00100138 NICADD Test Beam Simulator

  21. HCal Live Edep: 20 GeV Pi+ NICADD Test Beam Simulator

  22. HCal Abs. Edep: 20 GeV Pi+ NICADD Test Beam Simulator

  23. Edep in Inner Cells: 20 GeV Pi+ NICADD Test Beam Simulator

  24. ECal Live Edep: 20 GeV Pi+ NICADD Test Beam Simulator

  25. Visualization Scheme Tracks Cell Hits (all layers) Trajectory Charge Color Blue: positive Red: negative Green: neutral (Geant4 defaults) (front view of HCal & ECal hits) NICADD Test Beam Simulator

  26. All Trajectories and Lit Cells: Side View pi+ pos = (0, 0, 1990) mom = (0, 0, -1) E = 20 GeV all traj & lit cells NICADD Test Beam Simulator

  27. GPS Square: Side View 100 pi+ small square charged traj and hits NICADD Test Beam Simulator

  28. GPS Cone: Top View 20 pi+ theta = 0-10 phi = all mom = (0, 0, -1) E = 5-10 GeV charged traj & lit cells NICADD Test Beam Simulator

  29. GPS Cone: Front View 100 pi+ random circle charged traj NICADD Test Beam Simulator

  30. Hits and Tracks in ECal e- pos = (0, 0, 1990) mom = (0, 0, -1) E = 50 GeV all traj & hits NICADD Test Beam Simulator

  31. 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 NICADD Test Beam Simulator