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

Loading in 2 Seconds...

play fullscreen
1 / 31

NICADD Test Beam Simulator - PowerPoint PPT Presentation


  • 116 Views
  • Uploaded on

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

loader
I am the owner, or an agent authorized to act on behalf of the owner, of the copyrighted work described.
capcha
Download Presentation

PowerPoint Slideshow about 'NICADD Test Beam Simulator' - nora


An Image/Link below is provided (as is) to download presentation

Download Policy: Content on the Website is provided to you AS IS for your information and personal use and may not be sold / licensed / shared on other websites without getting consent from its author.While downloading, if for some reason you are not able to download a presentation, the publisher may have deleted the file from their server.


- - - - - - - - - - - - - - - - - - - - - - - - - - E N D - - - - - - - - - - - - - - - - - - - - - - - - - -
Presentation Transcript
nicadd test beam simulator
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

project description and goals
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

software libraries
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

classes
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

class diagram
Class Diagram

Association

Directed Assoc.

Inheritance

NICADD Test Beam Simulator

geometrymap
GeometryMap

GeometryMap

ComputeWorldVectors()

AddDaughters(vol, layno)

TBDetectorConstruction

TBSensitiveDetector

PhysVolInfo

GetLayNo()

GetVolType()

GetGlobId()

GetHit()

GetThreeVector()

ComputeWorldVector()

TBCellHit

NICADD Test Beam Simulator

detector geometry
Detector Geometry

TC

ECal

HCal

1.5 m

1 m

~2 m

NICADD Test Beam Simulator

geometry overview
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

cells
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

slide10
HCal

Front

Side (Angled)

30 cm

1.25 m

1 m

NICADD Test Beam Simulator

slide11
ECal

Front

Side (partial)

1 m

~6 cm

NICADD Test Beam Simulator

tail catcher
Tail Catcher

Side

Front (Horizontal Layer)

Front (transparent)

Front (Vertical Layer)

63 cm

1.5 m

1.5 m

NICADD Test Beam Simulator

command interface
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

particle generation
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

gps conical example
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

gps square example
GPS Square Example

10 pi+

pos = (rX,rY,1990)

square = 15 cm / side

E = 10-20 GeV

charged traj only

NICADD Test Beam Simulator

hits processing
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

end of event
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

text output format
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

output example
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

hcal live edep 20 gev pi
HCal Live Edep: 20 GeV Pi+

NICADD Test Beam Simulator

hcal abs edep 20 gev pi
HCal Abs. Edep: 20 GeV Pi+

NICADD Test Beam Simulator

edep in inner cells 20 gev pi
Edep in Inner Cells: 20 GeV Pi+

NICADD Test Beam Simulator

ecal live edep 20 gev pi
ECal Live Edep: 20 GeV Pi+

NICADD Test Beam Simulator

visualization scheme
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

all trajectories and lit cells side view
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

gps square side view
GPS Square: Side View

100 pi+

small square

charged traj

and hits

NICADD Test Beam Simulator

gps cone top view
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

gps cone front view
GPS Cone: Front View

100 pi+

random circle

charged traj

NICADD Test Beam Simulator

hits and tracks in ecal
Hits and Tracks in ECal

e-

pos = (0, 0, 1990)

mom = (0, 0, -1)

E = 50 GeV

all traj & hits

NICADD Test Beam Simulator

additional information
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