SND@LHC Software
Loading...
Searching...
No Matches
run_fixedTarget.py File Reference

Go to the source code of this file.

Namespaces

namespace  run_fixedTarget
 

Functions

 run_fixedTarget.get_work_dir (run_number, tag=None)
 
 run_fixedTarget.init ()
 

Variables

str run_fixedTarget.mcEngine = "TGeant4"
 
str run_fixedTarget.simEngine = "Pythia8"
 
int run_fixedTarget.runnr = 1
 
int run_fixedTarget.nev = 1000
 
bool run_fixedTarget.checkOverlap = True
 
bool run_fixedTarget.G4only = False
 
bool run_fixedTarget.storeOnlyMuons = False
 
bool run_fixedTarget.skipNeutrinos = False
 
bool run_fixedTarget.withEvtGen = True
 
int run_fixedTarget.boostDiMuon = 1.
 
int run_fixedTarget.boostFactor = 1.
 
bool run_fixedTarget.charm = False
 
bool run_fixedTarget.beauty = False
 
float run_fixedTarget.chicc = 1.7e-3
 
float run_fixedTarget.chibb = 1.6e-7
 
int run_fixedTarget.npot = 5E13
 
int run_fixedTarget.nStart = 0
 
str run_fixedTarget.charmInputFile = ROOT.gSystem.Getenv("EOSSHIP")+"/eos/experiment/ship/data/Charm/Cascade-parp16-MSTP82-1-MSEL4-76Mpot_1.root"
 
str run_fixedTarget.outputDir = "."
 
str run_fixedTarget.work_dir = "./"
 
float run_fixedTarget.ecut = 0.5
 
int run_fixedTarget.dy = 10.
 
int run_fixedTarget.dv = 6
 
int run_fixedTarget.ds = 9
 
int run_fixedTarget.nud = 3
 
 run_fixedTarget.logger = logging.getLogger(os.path.splitext(os.path.basename(os.sys.argv[0]))[0])
 
 run_fixedTarget.args = init()
 
 run_fixedTarget.ship_geo = ConfigRegistry.loadpy("$FAIRSHIP/geometry/geometry_config.py", Yheight = dy, tankDesign = dv, muShieldDesign = ds, nuTauTargetDesign=nud)
 
str run_fixedTarget.txt = 'pythia8_Geant4_'
 
str run_fixedTarget.outFile = outputDir+'/'+txt+str(runnr)+'_'+str(ecut)+'.root'
 
str run_fixedTarget.parFile = outputDir+'/ship.params.'+txt+str(runnr)+'_'+str(ecut)+'.root'
 
 run_fixedTarget.timer = ROOT.TStopwatch()
 
 run_fixedTarget.run = ROOT.FairRunSim()
 
 run_fixedTarget.rtdb = run.GetRuntimeDb()
 
 run_fixedTarget.cave = ROOT.ShipCave("CAVE")
 
 run_fixedTarget.TargetStation
 
 run_fixedTarget.slices_length = ROOT.std.vector('float')()
 
 run_fixedTarget.slices_material = ROOT.std.vector('std::string')()
 
 run_fixedTarget.MuonShield
 
 run_fixedTarget.sensPlane = ROOT.exitHadronAbsorber()
 
 run_fixedTarget.primGen = ROOT.FairPrimaryGenerator()
 
 run_fixedTarget.P8gen = ROOT.FixedTargetGenerator()
 
 run_fixedTarget.gMC = ROOT.TVirtualMC.GetMC()
 
 run_fixedTarget.fStack = gMC.GetStack()
 
 run_fixedTarget.gProcessTable = ROOT.G4ProcessTable.GetProcessTable()
 
 run_fixedTarget.procAnnihil = gProcessTable.FindProcess(ROOT.G4String('AnnihiToMuPair'),ROOT.G4String('e+'))
 
 run_fixedTarget.procGMuPair = gProcessTable.FindProcess(ROOT.G4String('GammaToMuPair'),ROOT.G4String('gamma'))
 
 run_fixedTarget.rtime = timer.RealTime()
 
 run_fixedTarget.ctime = timer.CpuTime()
 
str run_fixedTarget.tmpFile = outFile+"tmp"
 
 run_fixedTarget.fin = ROOT.gROOT.GetListOfFiles()[0]
 
 run_fixedTarget.fHeader = fin.FileHeader
 
 run_fixedTarget.poteq = P8gen.GetPotForCharm()
 
str run_fixedTarget.info = "POT equivalent = %7.3G"%(poteq)
 
str run_fixedTarget.conditions = " with ecut="+str(ecut)
 
 run_fixedTarget.nt = fin.Get('4DP')
 
 run_fixedTarget.tf = ROOT.TFile('FourDP.root','recreate')
 
 run_fixedTarget.tnt = nt.CloneTree(0)
 
 run_fixedTarget.rc = nt.GetEvent(i)
 
 run_fixedTarget.t = fin.cbmsim
 
 run_fixedTarget.fout = ROOT.TFile(tmpFile,'recreate' )
 
 run_fixedTarget.sTree = t.CloneTree(0)
 
int run_fixedTarget.nEvents = 0
 
 run_fixedTarget.x = fin.Get(k.GetName())
 
 run_fixedTarget.className = x.Class().GetName()
 
 run_fixedTarget.xcopy = x.Clone()
 
 run_fixedTarget.ff = fin.FileHeader.Clone(fout.GetName())
 
 run_fixedTarget.rc1 = os.system("rm "+outFile)
 
 run_fixedTarget.rc2 = os.system("mv "+tmpFile+" "+outFile)
 
 run_fixedTarget.sGeo = ROOT.gGeoManager