NICADD Test Beam Simulator - PowerPoint PPT Presentation

Presentation Transcript

NICADD Test Beam Simulator

Project Description

Classes

Geometry

Hits Processing

Command Interface

Particle Generation

Text Output Format

Histograms

Visualization

Additional Information

Jeremy I. McCormick

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

Software Libraries

Geant4.5.2.p01
GPS
g4phys SO (Mokka v01-04 Physics List)
popt for CL
Composite Materials def from PEGS4 (CERN)

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

Class Diagram

Association

Directed Assoc.

Inheritance

GeometryMap

ComputeWorldVectors()

AddDaughters(vol, layno)

TBDetectorConstruction

TBSensitiveDetector

PhysVolInfo

GetLayNo()

GetVolType()

GetGlobId()

GetHit()

GetThreeVector()

ComputeWorldVector()

TBCellHit

Detector Geometry

TC

ECal

HCal

1.5 m

1 m

~2 m

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

-

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

HCal

Front

Side (Angled)

30 cm

1.25 m

1 m

ECal

Front

Side (partial)

1 m

~6 cm

Tail Catcher

Side

Front (Horizontal Layer)

Front (transparent)

Front (Vertical Layer)

63 cm

1.5 m

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

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

GPS Conical Example

10 pi+

pos = (0, 0 ,1990)

theta = 0-10

phi = all

E = 2-10 GeV

charged traj only

GPS Square Example

10 pi+

pos = (rX,rY,1990)

square = 15 cm / side

E = 10-20 GeV

charged traj only

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

}

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)

Text Output Format

EVENT evt_no

ECAL no_hits

layer_no … x y z cell_edep abs_edep

track_no track_edep

…

HCAL …

…

TAILCATCHER …

…

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

HCal Live Edep: 20 GeV Pi+

HCal Abs. Edep: 20 GeV Pi+

Edep in Inner Cells: 20 GeV Pi+

ECal Live Edep: 20 GeV Pi+

Visualization Scheme

Tracks

Cell Hits (all layers)

Trajectory Charge

Color

Blue: positive

Red: negative

Green: neutral

(Geant4 defaults)

(front view of HCal & ECal hits)

All Trajectories and Lit Cells: Side View

pi+

pos = (0, 0, 1990)

mom = (0, 0, -1)

E = 20 GeV

all traj & lit cells

GPS Square: Side View

100 pi+

small square

charged traj

and hits

GPS Cone: Top View

20 pi+

theta = 0-10

phi = all

mom = (0, 0, -1)

E = 5-10 GeV

charged traj & lit cells

GPS Cone: Front View

100 pi+

random circle

charged traj

Hits and Tracks in ECal

e-

pos = (0, 0, 1990)

mom = (0, 0, -1)

E = 50 GeV

all traj & hits

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