Jeremy%20I.%20McCormick

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

• Test Beam Prototype
• Hadron EM Calorimeters, Tail Catcher
• HCal ECal layering similar to SDJan03
• Polystyrene and Silicon sensitive regions
• one volume per cell geometries
• aggregate step hits by cell with track
  contributions

3
Libraries

• Geant4.5.2.p01
• GPS
• g4phys SO
• popt for CL
• Composite Materials def from PEGS4 (CERN)

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

5
Class Diagram Associations
6
GeometryMap
GeometryMap ComputeWorldVectors()
AddDaughters(vol, layno)
TBDetectorConstruction
TBSensitiveDetector
PhysVolInfo GetLayNo() GetVolType()
GetGlobId() GetHit() GetThreeVector()
ComputeWorldVector()
TBCellHit
7
Detector Geometry
TC
ECal
HCal
1.5 m
1 m
2 m
8
Geometry Overview
9
Cells
1.5 m
10
HCal
Front
Side (Angled)
30 cm
1.25 m
1 m
11
ECal
Front
Side
1 m
6 cm
12
Tail Catcher
Front (Hor)
Side
Front (transparent)
Front (Vert)
63 cm
1.5 m
1.5 m
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
• /tbeam/doRun no_events

14
Particle Generation

• General Particle Source
• 1 to n MC particles
• random starting regions
• random theta / phi within max / min
• random energy
• /tbeam/doRun no_events for multiple events

15
GPS Conical Example
Pi pos (0, 0 ,1990) theta 0-10 phi all e
2-10 GeV charged traj only
16
GPS Square Example
Pi pos (rX,rY,1990) square 15 cm / side e
10-20 GeV charged traj only
17
Hits Processing
TBSensitiveDetectorProcessHits() PhysVolInfo
pinfo PhysVolMap-gtGetPhysVolInfo(physVol)
TBCellHit theHit CellHit(physVol) if
(!theHit) theHit new TBCellHit() if
(pinfo) theHit-gtSetPhysVolInfo(pinfo)
mapHits.insert(HitMapPair(physVol,theHit))
PhysVolInfoPhysVolType vtype
pinfo-gtGetVolType() if (vtype
PhysVolInfoCELL) theHit-gtIncrEdep(edep)
TBTrackInformation tinfo (TBTrackInformation
)(aStep-gtGetTrack()-gtGetUserInformation())
theHit-gtIncrTrkEdep(tinfo-gtGetTrackID(),edep)
else if (vtype PhysVolInfoABSORBER)
theHit-gtIncrAbsEdep(edep)
18
End of Event
TBSensitiveDetectorEndOfEvent() Copy Hits
to HC TBEventActionEndOfEventAction()
WriteEvent(evt) WriteHitsCollection(HC)
aHit-gtWrite(ostr)
TrackEdep-gtWrite(ostr) aHit-gtDraw()
VisualizeTrajectories(evt)
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
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
21
Hist 1
22
Hist 2
23
Visualization Scheme
Cell Hits (all layers)
Tracks
Trajectory Charge
Color
Blue positive Red negative Green
neutral
(front view of HCal ECal hits)
24
All Trajectories and Lit Cells Side View
pi pos (0, 0, 1990) mom (0, 0, -1) E 20
GeV all traj lit cells
25
GPS Square Side View
26
GPS Cone Top View
20 pi theta 0-10 phi all mom (0, 0, -1) E
5-10 GeV charged traj lit cells
27
GPS Cone Front View
100 pi charged traj
28
Hits and Tracks in ECal
e- pos (0, 0, 1990) mom (0, 0, -1) E 50
GeV all traj hits
29
Additional Information

• Homepage http//nicadd.niu.edu/jeremy/g4proj
• Source Browser http//nicadd.niu.edu/jeremy/tbe
  am/html/
• CVSROOTextlcdguest_at_nicadd.niu.edu/cvs
• CVS_RSHssh
• Password lcd_2004
• g4phys http//nicadd.niu.edu/jeremy/download/g4
  phys.tgz
• popt http//www.gnu.org/directory/libs/popt.html