Functions | |
| def | configure |
| def | run |
| def NeutronSim::configure | ( | argv = [] |
) |
Configure this module with neutron energy
Definition at line 13 of file NeutronSim.py.
00013 : 00014 """Configure this module with neutron energy""" 00015 00016 import sys, getopt 00017 opts,args = getopt.getopt(argv,"k:") 00018 neutronE = 100.0 00019 for opt,arg in opts: 00020 if opt == "-k": 00021 neutronE = float(arg) 00022 print "======================================================" 00023 print "Neutron energy = ", neutronE, " eV" 00024 print "======================================================" 00025 00026 # Set run info 00027 from RunDataSvc.RunDataSvcConf import RunDataSvc 00028 runDataSvc = RunDataSvc() 00029 runDataSvc.SimRunType = "Physics" 00030 00031 # Generator 00032 import GaudiKernel.SystemOfUnits as units 00033 from GenTools.Helpers import Gun 00034 from GenTools.GenToolsConf import GtGunGenTool, GtPositionerTool 00035 gun = Gun(name = "neutronGun", 00036 volume = "/dd/Structure/AD/db-oil1", 00037 gun = GtGunGenTool("gun", 00038 ParticlesPerEvent = 1, 00039 ParticleName = "n0", 00040 Momentum = neutronE*units.eV, 00041 MomentumMode = "Fixed", 00042 MomentumSpread = 0.*units.MeV, 00043 DirectionMode = "Uniform", 00044 Direction = [ 1, 0, 0 ], 00045 DirectionSpread = 3), 00046 positioner = GtPositionerTool("pos", 00047 Strategy = "FullVolume", 00048 Mode = "Uniform", 00049 Spread = 2.5*units.meter, 00050 Position = [0,0,0*units.meter]) 00051 ) 00052 gun.timerator.LifeTime = 0.1*units.second 00053 import GenTools 00054 GenTools.Configure().register(gun) 00055 00056 # Detector Simulation 00057 import DetSim 00058 detsim = DetSim.Configure() 00059 detsim.historian(trackSelection="(pdg == 2112)",\ 00060 vertexSelection="(pdg == 2112)") 00061 params = { 00062 'start' :"(start > 0)", 00063 'track1':"(id==1)", 00064 'track2':"(id==2)", 00065 'GD': "MaterialName == '/dd/Materials/GdDopedLS'", 00066 'LS': "MaterialName == '/dd/Materials/LiquidScintillator'", 00067 'MO': "MaterialName == '/dd/Materials/MineralOil'", 00068 'IAV': "DetectorElementName == 'db-iav1'", 00069 'OAV': "DetectorElementName == 'db-oav1'", 00070 'IWS': "MaterialName == '/dd/Materials/IwsWater'", 00071 'OWS': "MaterialName == '/dd/Materials/OwsWater'", 00072 'lastvtx': "IsStopping == 1", 00073 'firstvtx': "IsStarting == 1", 00074 'Neutron': "pdg == 2112", 00075 'NeutronMom': "creator == 2112", 00076 'Gamma': "pdg == 22", 00077 'Muon': "(pdg == 13 or pdg == -13)" 00078 } 00079 00080 detsim.unobserver(stats=[ 00081 ["EDepInGdLS", "dE", "%(GD)s"%params], 00082 ["EDepInLS", "dE", "%(LS)s"%params], 00083 ["EDepInIAV", "dE", "%(IAV)s"%params], 00084 ["EDepInOAV", "dE", "%(OAV)s"%params], 00085 ["EDepInOIL", "dE", "%(MO)s"%params], 00086 00087 ["QEDepInGdLS", "qdE", "%(GD)s"%params], 00088 ["QEDepInLS", "qdE", "%(LS)s"%params], 00089 ["QEDepInIAV", "qdE", "%(IAV)s"%params], 00090 ["QEDepInOAV", "qdE", "%(OAV)s"%params], 00091 ["QEDepInOIL", "qdE", "%(MO)s"%params], 00092 00093 ["tQESumGdLS", "qEt", "%(GD)s"%params], 00094 ["xQESumGdLS", "qEx", "%(GD)s"%params], 00095 ["yQESumGdLS", "qEy", "%(GD)s"%params], 00096 ["zQESumGdLS", "qEz", "%(GD)s"%params], 00097 00098 ["tQESumLS", "qEt", "%(LS)s"%params], 00099 ["xQESumLS", "qEx", "%(LS)s"%params], 00100 ["yQESumLS", "qEy", "%(LS)s"%params], 00101 ["zQESumLS", "qEz", "%(LS)s"%params], 00102 00103 ["tQESumMO", "qEt", "%(MO)s"%params], 00104 ["xQESumMO", "qEx", "%(MO)s"%params], 00105 ["yQESumMO", "qEy", "%(MO)s"%params], 00106 ["zQESumMO", "qEz", "%(MO)s"%params], 00107 00108 # track 1 00109 ["pdgId_Trk1","pdg","%(track1)s and %(start)s"%params], 00110 ["t_Trk1", "t" , "%(track1)s and %(start)s"%params], 00111 ["lx_Trk1", "lx", "%(track1)s and %(start)s"%params], 00112 ["ly_Trk1", "ly", "%(track1)s and %(start)s"%params], 00113 ["lz_Trk1", "lz", "%(track1)s and %(start)s"%params], 00114 ["tEnd_Trk1", "t" , "%(track1)s and %(lastvtx)s"%params], 00115 ["lxEnd_Trk1", "lx", "%(track1)s and %(lastvtx)s"%params], 00116 ["lyEnd_Trk1", "ly", "%(track1)s and %(lastvtx)s"%params], 00117 ["lzEnd_Trk1", "lz", "%(track1)s and %(lastvtx)s"%params], 00118 ["e_Trk1", "E", "%(track1)s and %(start)s"%params], 00119 ["p_Trk1", "p", "%(track1)s and %(start)s"%params], 00120 ["ke_Trk1", "KE", "%(track1)s and %(start)s"%params], 00121 ["vx_Trk1", "lvx","%(track1)s and %(start)s"%params], 00122 ["vy_Trk1", "lvy","%(track1)s and %(start)s"%params], 00123 ["vz_Trk1", "lvz","%(track1)s and %(start)s"%params], 00124 ["TrkLength_GD_Trk1", "dx","%(track1)s and %(GD)s"%params], 00125 ["TrkLength_iAV_Trk1", "dx","%(track1)s and %(IAV)s"%params], 00126 ["TrkLength_LS_Trk1", "dx","%(track1)s and %(LS)s"%params], 00127 ["TrkLength_oAV_Trk1", "dx","%(track1)s and %(OAV)s"%params], 00128 ["TrkLength_Oil_Trk1", "dx","%(track1)s and %(MO)s"%params], 00129 00130 ["capTarget", "capTargetZ","%(track1)s and %(lastvtx)s"%params] 00131 00132 ]) 00133 00134 import ElecSim 00135 elecsim = ElecSim.Configure() 00136 00137 import TrigSim 00138 trigsim = TrigSim.Configure() 00139 00140 import ReadoutSim 00141 rosim = ReadoutSim.Configure() 00142 def run(app):
| def NeutronSim::run | ( | app | ) |
1.4.7