In This Package:

diffuser_ball.py

Go to the documentation of this file.

#!/usr/bin/env python

'''
usage example:

  nuwa.py -n 10 diffuser_ball.py 

'''

### GenTools
import GenTools
from GenTools.Helpers import DiffuserBall
import GaudiKernel.SystemOfUnits as units
ball = DiffuserBall()
ball.setPosition( [0.0*units.cm, 0.0*units.cm, 0.0*units.cm] )
#ball.ball.Anisotropy = False
ball.ball.BrightSpotTheta = 180*units.degree
ball.ball.BrightSpotPhi = 0
GenTools.Configure().register(ball)

#dumper = GenTools.Dumper()

### DetSim
import DetSim
detsim = DetSim.Configure(physlist=DetSim.physics_list_basic)
params = {
    'start' :"(start > 0)",
    'track1':"(id==1 and ProcessType==1)",
    'track2':"(id==2 and ProcessType==1)",
    'GD':    "MaterialName == '/dd/Materials/GdDopedLS'",
    'LS':    "MaterialName == '/dd/Materials/LiquidScintillator'",
    'oil':   "MaterialName == '/dd/Materials/MineralOil'",
    'iAV':   "MaterialName == '/dd/Materials/Acrylic'",
    'oAV':   "MaterialName == '/dd/Materials/Acrylic'"
    }
detsim.unobserver(stats=[
    # for track 1
    ["pdgId_Trk1","pdg","%(track1)s and %(start)s"%params],
    ["t_Trk1",    "t" , "%(track1)s and %(start)s"%params],
    ["x_Trk1",    "x", "%(track1)s and %(start)s"%params],
    ["y_Trk1",    "y", "%(track1)s and %(start)s"%params],
    ["z_Trk1",    "z", "%(track1)s and %(start)s"%params],
    ["e_Trk1",    "E",  "%(track1)s and %(start)s"%params],
    ["p_Trk1",    "p",  "%(track1)s and %(start)s"%params],
    ["ke_Trk1",   "KE", "%(track1)s and %(start)s"%params],
    ["vx_Trk1",   "lvx","%(track1)s and %(start)s"%params],
    ["vy_Trk1",   "lvy","%(track1)s and %(start)s"%params],
    ["vz_Trk1",   "lvz","%(track1)s and %(start)s"%params],
    ["TrkLength_GD_Trk1",  "dx","%(track1)s and %(GD)s"%params],
    ["TrkLength_iAV_Trk1", "dx","%(track1)s and %(iAV)s"%params],
    ["TrkLength_LS_Trk1",  "dx","%(track1)s and %(LS)s"%params],
    ["TrkLength_oAV_Trk1", "dx","%(track1)s and %(oAV)s"%params],
    ["TrkLength_Oil_Trk1", "dx","%(track1)s and %(oil)s"%params],
    # for track 2 now
    ["pdgId_Trk2","pdg","%(track2)s and %(start)s"%params],
    ["t_Trk2",    "t" , "%(track2)s and %(start)s"%params],
    ["x_Trk2",    "x", "%(track2)s and %(start)s"%params],
    ["y_Trk2",    "y", "%(track2)s and %(start)s"%params],
    ["z_Trk2",    "z", "%(track2)s and %(start)s"%params],
    ["e_Trk2",    "E",  "%(track2)s and %(start)s"%params],
    ["p_Trk2",    "p",  "%(track2)s and %(start)s"%params],
    ["ke_Trk2",   "KE", "%(track2)s and %(start)s"%params],
    ["vx_Trk2",   "lvx","%(track2)s and %(start)s"%params],
    ["vy_Trk2",   "lvy","%(track2)s and %(start)s"%params],
    ["vz_Trk2",   "lvz","%(track2)s and %(start)s"%params],
    ["TrkLength_GD_Trk2",  "dx","%(track2)s and %(GD)s"%params],
    ["TrkLength_iAV_Trk2", "dx","%(track2)s and %(iAV)s"%params],
    ["TrkLength_LS_Trk2",  "dx","%(track2)s and %(LS)s"%params],
    ["TrkLength_oAV_Trk2", "dx","%(track2)s and %(oAV)s"%params],
    ["TrkLength_Oil_Trk2", "dx","%(track2)s and %(oil)s"%params]
    ])

## DetSimVali
from DetSimValidation.DetSimValidationConf import DetSimVali
dsv = DetSimVali()
dsv.Volume = '/dd/Structure/AD/far-oil1'

from Gaudi.Configuration import ApplicationMgr
app = ApplicationMgr()
app.TopAlg.append(dsv)

from GaudiSvc.GaudiSvcConf import THistSvc
histsvc = THistSvc()
histsvc.Output =["file1 DATAFILE='../mc_ntuple/diffuser.root' OPT='RECREATE' TYP='ROOT' "]

### ElecSim
import ElecSim
elecsim = ElecSim.Configure()

### TrigSim
import TrigSim
trigsim = TrigSim.Configure()

### ReadoutSim
import ReadoutSim
rosim = ReadoutSim.Configure()

---