| Classes | Job Modules | Data Objects | Services | Algorithms | Tools | Packages | Directories | Tracs |

In This Package:

TrigSim::gammaESum::__init__

TrigSim::gammaESum::__init__ Namespace Reference


Functions

def configure
def run

Variables

string DATETIME_FORMAT = '%Y-%m-%dT%H:%M:%S'

Function Documentation

def TrigSim::gammaESum::__init__::configure (   argv = []  )  

Configure this module with gamma energy

Definition at line 26 of file __init__.py.

00026                         :
00027     """Configure this module with gamma energy"""
00028     
00029     import sys, getopt
00030     from time import localtime, gmtime, mktime, strftime, strptime, timezone
00031     opts,args = getopt.getopt(argv,"p:w:k:v:")
00032     wallTime = 0
00033     gammaE = 1.0
00034     pmtDataPath = None
00035     tagVolume = "/dd/Structure/AD/db-oil1"
00036     for opt,arg in opts:
00037         if opt == "-p":
00038             pmtDataPath = arg
00039         if opt == "-w":
00040             if -1 != arg.find('T'):
00041                 wallTime = int(mktime(strptime(arg,
00042                                                DATETIME_FORMAT)) - timezone)
00043             else:
00044                 wallTime = int(arg)
00045  
00046         if opt == "-k":
00047             gammaE = float(arg)
00048             print "======================================================"
00049             print "Photon energy = ", gammaE, " MeV"
00050             print "======================================================"
00051 
00052         if opt == "-v":
00053             tagVolume = arg
00054             print "======================================================"
00055             print "Photon energy = ", gammaE, " MeV"
00056             print "======================================================"
00057  
00058     print "======================================================"
00059     print "Begin JOB TIME = ", strftime(DATETIME_FORMAT,
00060                                               gmtime())
00061     print "======================================================"
00062 
00063     # Generator
00064     import GaudiKernel.SystemOfUnits as units
00065     from GenTools.Helpers import Gun
00066     from GenTools.GenToolsConf import GtGunGenTool, GtPositionerTool
00067     gun = Gun(name = "gammaGun",
00068               volume = tagVolume,
00069               gun = GtGunGenTool("gun",
00070                                  ParticlesPerEvent = 1,
00071                                  ParticleName = "gamma",
00072                                  Momentum = gammaE*units.MeV,
00073                                  MomentumMode = "Fixed",
00074                                  MomentumSpread = 0.*units.MeV,
00075                                  DirectionMode = "Uniform",
00076                                  Direction = [ 1, 0, 0 ],
00077                                  DirectionSpread = 3),
00078               positioner = GtPositionerTool("pos",
00079                                             Strategy = "FullVolume",
00080                                             Mode = "Uniform",
00081                                             Spread = 2.5*units.meter,
00082                                             Position = [0,0,0*units.meter])
00083               )
00084     gun.timerator.LifeTime = 0.020*units.second
00085     import GenTools
00086     GenTools.Configure().register(gun)
00087     
00088     # Detector Simulation
00089     import DetSim
00090     detsim = DetSim.Configure(site="dayabay", \
00091                               physlist = DetSim.physics_list_basic)
00092     detsim.historian(trackSelection="(pdg == 22)",\
00093                      vertexSelection="(pdg == 22)")
00094     params = {
00095         'start' :"(start > 0)",
00096         'track1':"(id==1)",
00097         'track2':"(id==2)",
00098         'GD':    "MaterialName == '/dd/Materials/GdDopedLS'",
00099         'LS':    "MaterialName == '/dd/Materials/LiquidScintillator'",
00100         'MO':   "MaterialName == '/dd/Materials/MineralOil'",
00101         'IAV':   "DetectorElementName == 'db-iav1'",
00102         'OAV':   "DetectorElementName == 'db-oav1'",
00103         'IWS': "MaterialName == '/dd/Materials/IwsWater'",
00104         'OWS': "MaterialName == '/dd/Materials/OwsWater'",
00105         'lastvtx': "IsStopping == 1",
00106         'firstvtx': "IsStarting == 1",
00107         'Neutron': "pdg == 2112",
00108         'NeutronMom': "creator == 2112",
00109         'Gamma': "pdg == 22",
00110         'Muon': "(pdg == 13 or pdg == -13)"
00111         }
00112 
00113     detsim.unobserver(stats=[
00114             ["EDepInGdLS", "dE", "%(GD)s"%params],
00115             ["EDepInLS", "dE", "%(LS)s"%params],
00116             ["EDepInIAV", "dE", "%(IAV)s"%params],
00117             ["EDepInOAV", "dE", "%(OAV)s"%params],
00118             ["EDepInOIL", "dE", "%(MO)s"%params],
00119 
00120             ["QEDepInGdLS", "qdE", "%(GD)s"%params],
00121             ["QEDepInLS", "qdE", "%(LS)s"%params],
00122             ["QEDepInIAV", "qdE", "%(IAV)s"%params],
00123             ["QEDepInOAV", "qdE", "%(OAV)s"%params],
00124             ["QEDepInOIL", "qdE", "%(MO)s"%params],
00125 
00126             ["tQESumGdLS", "qEt", "%(GD)s"%params],
00127             ["xQESumGdLS", "qEx", "%(GD)s"%params],
00128             ["yQESumGdLS", "qEy", "%(GD)s"%params],
00129             ["zQESumGdLS", "qEz", "%(GD)s"%params],
00130 
00131             ["tQESumLS", "qEt", "%(LS)s"%params],
00132             ["xQESumLS", "qEx", "%(LS)s"%params],
00133             ["yQESumLS", "qEy", "%(LS)s"%params],
00134             ["zQESumLS", "qEz", "%(LS)s"%params],
00135 
00136             ["tQESumMO", "qEt", "%(MO)s"%params],
00137             ["xQESumMO", "qEx", "%(MO)s"%params],
00138             ["yQESumMO", "qEy", "%(MO)s"%params],
00139             ["zQESumMO", "qEz", "%(MO)s"%params],
00140             
00141             # track 1
00142             ["pdgId_Trk1","pdg","%(track1)s and %(start)s"%params],
00143             ["t_Trk1",    "t" , "%(track1)s and %(start)s"%params],
00144             ["x_Trk1",    "x", "%(track1)s and %(start)s"%params],
00145             ["y_Trk1",    "y", "%(track1)s and %(start)s"%params],
00146             ["z_Trk1",    "z", "%(track1)s and %(start)s"%params],
00147             ["tEnd_Trk1",    "t" , "%(track1)s and %(lastvtx)s"%params],
00148             ["xEnd_Trk1",    "x", "%(track1)s and %(lastvtx)s"%params],
00149             ["yEnd_Trk1",    "y", "%(track1)s and %(lastvtx)s"%params],
00150             ["zEnd_Trk1",    "z", "%(track1)s and %(lastvtx)s"%params],
00151             ["e_Trk1",    "E",  "%(track1)s and %(start)s"%params],
00152             ["p_Trk1",    "p",  "%(track1)s and %(start)s"%params],
00153             ["ke_Trk1",   "KE", "%(track1)s and %(start)s"%params],
00154             ["vx_Trk1",   "lvx","%(track1)s and %(start)s"%params],
00155             ["vy_Trk1",   "lvy","%(track1)s and %(start)s"%params],
00156             ["vz_Trk1",   "lvz","%(track1)s and %(start)s"%params],
00157             ["TrkLength_GD_Trk1",  "dx","%(track1)s and %(GD)s"%params],
00158             ["TrkLength_iAV_Trk1", "dx","%(track1)s and %(IAV)s"%params],
00159             ["TrkLength_LS_Trk1",  "dx","%(track1)s and %(LS)s"%params],
00160             ["TrkLength_oAV_Trk1", "dx","%(track1)s and %(OAV)s"%params],
00161             ["TrkLength_Oil_Trk1", "dx","%(track1)s and %(MO)s"%params]
00162             ])
00163             
00164     import ElecSim
00165     elecsim = ElecSim.Configure()
00166     if pmtDataPath != None:
00167         # change PMT properties
00168         elecsim.dataSvc.setPmtSimData( pmtDataPath )
00169     
00170     import TrigSim
00171     trigsim = TrigSim.Configure()
00172 
00173 
00174     import TrigSim.TrigSimConf as TsConf
00175     TsConf.TsTriggerAlg().TrigTools=["TsESumTriggerTool"]
00176     TsConf.TsESumTriggerTool().OutputLevel=VERBOSE
00177     TsConf.TsESumTriggerTool().ADThreshold=0.07*SystemOfUnits.volt
00178 
00179     import ReadoutSim
00180     rosim = ReadoutSim.Configure()
00181 
def run(app):

def TrigSim::gammaESum::__init__::run (   app  ) 

Definition at line 182 of file __init__.py.

00182             :
00183     pass
00184

Variable Documentation

string TrigSim::gammaESum::__init__::DATETIME_FORMAT = '%Y-%m-%dT%H:%M:%S' [static]

Definition at line 20 of file __init__.py.

| Classes | Job Modules | Data Objects | Services | Algorithms | Tools | Packages | Directories | Tracs |
Generated on Mon Apr 11 20:50:05 2011 for TrigSim by doxygen 1.4.7