NICADD Test Beam Simulato r 10/28/03 1 Jeremy I. McCormick Project Description Classes Geometry Hits Processing Command Interface Particle Generation Text Output Format Histograms Visualization Additional Information NICADD Test Beam Simulator
Jan 14, 2016
NICADD Test Beam Simulator10/28/03 1
Jeremy I. McCormick
Project DescriptionClasses
Geometry
Hits Processing
Command Interface
Particle Generation
Text Output Format
Histograms
Visualization
Additional Information
NICADD Test Beam Simulator
NICADD Test Beam Simulator10/28/03 2
•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
Project Description and Goals
NICADD Test Beam Simulator10/28/03 3
• Geant4.5.2.p01
• GPS
• g4phys SO (Mokka v01-04 Physics List)
• popt for CL
• Composite Materials def from PEGS4 (CERN)
Software Libraries
NICADD Test Beam Simulator10/28/03 4
G4 Utility• Cellulizer• BoxCellulizer• MaterialMap
User• TBEventAction• TBPrimaryGeneratorAction• TBRunAction• TBSteppingAction• TBTracking Action• TBVisManager
Geometry• GeometryMap• PhysVolInfo• TBDetectorConstruction• TBCal• TBTailCatcher• TBECal• TBHCal
Hits and Tracks• TBCellHit
• TBSensitiveDetector
• TBTrackInformation
• TBEdepMap
Classes
NICADD Test Beam Simulator10/28/03 5
Class Diagram
Directed Assoc.
Association
Inheritance
NICADD Test Beam Simulator10/28/03 6
GeometryMap ComputeWorldVectors() AddDaughters(vol, layno)
PhysVolInfo GetLayNo() GetVolType() GetGlobId() GetHit() GetThreeVector() ComputeWorldVector()
TBCellHit
TBDetectorConstructionTBSensitiveDetector
GeometryMap
NICADD Test Beam Simulator10/28/03 7
ECal HCal
TC
~2 m
1.5 m1 m
Detector Geometry
NICADD Test Beam Simulator10/28/03 8
Same as HCal1.5 x 1.56 *
10.5 cm =
63 cm
Tail Catcher
-max 1.5 x 1.5203 cmTotal
MaterialsX x Y dim. (m)No. Layers * thickness =
total thickness
Component
Stainless Steel, Polystyrene
1 x 150 *
2.5 cm =
125 cm
HCal
W, G10, Si, Cu, Air
1 x 130 *
5 mm =
15 cm
ECal
Geometry Overview
NICADD Test Beam Simulator10/28/03 9
323,290-86Total
90156Tail Catcher
300,00010,00030ECal
23,200100 inner +
364 outer = 464
50HCal
Total CellsCells / Layer
LayersComponent
1.5 m
Cells
NICADD Test Beam Simulator10/28/03 10
Front Side (Angled)
1 m 1.25 m
30 cm
HCal
NICADD Test Beam Simulator10/28/03 11
Front Side (partial)
1 m ~6 cm
ECal
NICADD Test Beam Simulator10/28/03 12
Front (Horizontal Layer)
Front (Vertical Layer)
SideFront (transparent)
1.5 m1.5 m
63 cm
Tail Catcher
NICADD Test Beam Simulator10/28/03 13
• 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
Command Interface
NICADD Test Beam Simulator10/28/03 14
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
Particle Generation
NICADD Test Beam Simulator10/28/03 15
10 pi+pos = (0, 0 ,1990)theta = 0-10phi = allE = 2-10 GeVcharged traj only
GPS Conical Example
NICADD Test Beam Simulator10/28/03 16
10 pi+pos = (rX,rY,1990)square = 15 cm / sideE = 10-20 GeVcharged traj only
GPS Square Example
NICADD Test Beam Simulator10/28/03 17
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); }}
Hits Processing
NICADD Test Beam Simulator10/28/03 18
TBSensitiveDetector::EndOfEvent()
{ Copy Hits to HC }
TBEventAction::EndOfEventAction()
WriteEvent(evt)
WriteHitsCollection(HC)
aHit->Write(ostr)
TrackEdep->Write(ostr)
aHit->Draw()
VisualizeTrajectories(evt)
End of Event
NICADD Test Beam Simulator10/28/03 19
EVENT evt_no
ECAL no_hits
layer_no … x y z cell_edep abs_edep
track_no track_edep
…
HCAL …
…
TAILCATCHER …
…
Text Output Format
NICADD Test Beam Simulator10/28/03 20
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
Output Example
NICADD Test Beam Simulator10/28/03 21
HCal Live Edep: 20 GeV Pi+
NICADD Test Beam Simulator10/28/03 22
HCal Abs. Edep: 20 GeV Pi+
NICADD Test Beam Simulator10/28/03 23
Edep in Inner Cells: 20 GeV Pi+
NICADD Test Beam Simulator10/28/03 24
ECal Live Edep: 20 GeV Pi+
NICADD Test Beam Simulator10/28/03 25
Cell Hits (all layers)Tracks
Blue: positiveRed: negativeGreen: neutral
Trajectory Charge
Color
(front view of HCal & ECal hits)
Visualization Scheme
(Geant4 defaults)
NICADD Test Beam Simulator10/28/03 26
pi+pos = (0, 0, 1990)mom = (0, 0, -1)E = 20 GeVall traj & lit cells
All Trajectories and Lit Cells: Side View
NICADD Test Beam Simulator10/28/03 27
100 pi+small squarecharged traj and hits
GPS Square: Side View
NICADD Test Beam Simulator10/28/03 28
20 pi+theta = 0-10phi = allmom = (0, 0, -1)E = 5-10 GeVcharged traj & lit cells
GPS Cone: Top View
NICADD Test Beam Simulator10/28/03 29
100 pi+random circlecharged traj
GPS Cone: Front View
NICADD Test Beam Simulator10/28/03 30
Hits and Tracks in ECal
e-pos = (0, 0, 1990)mom = (0, 0, -1)E = 50 GeVall traj & hits
NICADD Test Beam Simulator10/28/03 31
•Homepage: http://nicadd.niu.edu/~jeremy/g4proj
•Source Browser: http://nicadd.niu.edu/~jeremy/tbeam/html/
•CVS (bash):
export CVSROOT=:ext:[email protected]:/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
Additional Information