SND@LHC Software
Loading...
Searching...
No Matches
Functions | Variables
run_simScript Namespace Reference

Functions

 visualizeMagFields ()
 
 checkOverlapsWithGeant4 ()
 

Variables

 IgnoreCommandLineOptions
 
int debug = 0
 
bool dryrun = False
 
bool DownScaleDiMuon = False
 
float theHNLMass = 1.0*u.GeV
 
 theProductionCouplings = None
 
float theDPmass = 0.2*u.GeV
 
bool motherMode = True
 
str mcEngine = "TGeant4"
 
str simEngine = "Pythia8"
 
str inclusive = "c"
 
bool MCTracksWithHitsOnly = False
 
bool MCTracksWithEnergyCutOnly = True
 
bool MCTracksWithHitsOrEnergyCut = False
 
bool charmonly = False
 
bool HNL = True
 
str inputFile = "/eos/experiment/ship/data/Charm/Cascade-parp16-MSTP82-1-MSEL4-978Bpot.root"
 
bool defaultInputFile = True
 
dict globalDesigns
 
str default = '2018'
 
bool inactivateMuonProcesses = False
 
bool checking4overlaps = False
 
 parser = ArgumentParser()
 
 group = parser.add_mutually_exclusive_group()
 
 dest
 
 help
 
 required
 
 False
 
 action
 
 type
 
 float
 
 None
 
 int
 
 options = parser.parse_args()
 
 Opt_high = int(options.cosmics)
 
 theMass
 
list theCouplings = [float(c) for c in options.thecouplings.split(",")]
 
list theDecayCouplings = [float(c) for c in options.thedeccouplings.split(",")]
 
 ship_geo = ConfigRegistry.loadpy("$FAIRSHIP/geometry/charm-geometry_config.py", Setup = options.CharmdetSetup, cTarget = options.CharmTarget)
 
 Yheight
 
 dy
 
 tankDesign
 
 muShieldDesign
 
 ds
 
 nuTauTargetDesign
 
 nud
 
 CaloDesign
 
 strawDesign
 
 muShieldGeo
 
 muShieldStepGeo
 
 muShieldWithCobaltMagnet
 
 tag
 
str outFile = "%s/ship.%s.root" % (options.outputDir, tag)
 
str parFile = "%s/ship.params.%s.root" % (options.outputDir, tag)
 
 timer = ROOT.TStopwatch()
 
 run = ROOT.FairRunSim()
 
 rtdb = run.GetRuntimeDb()
 
 modules = shipDet_conf.configure(run,ship_geo)
 
 primGen = ROOT.FairPrimaryGenerator()
 
 P8gen = ROOT.HNLPythia8Generator()
 
 passDPconf = pythia8darkphoton_conf.configure(P8gen,options.theMass,options.theDPepsilon,inclusive, motherMode, options.deepCopy)
 
 target = ship_geo.target
 
 startZ = target.z0
 
 lengthZ = target.length
 
 endZ = startZ + lengthZ
 
int SmearBeam = 1*u.cm
 
tuple Lmin = ((ship_geo.Chamber1.z - ship_geo.chambers.Tub1length) - ship_geo.target.z0)/100.
 
tuple Lmax = (ship_geo.TrackStation1.z - ship_geo.target.z0)/100.
 
 ALPACAgen = ROOT.ALPACAGenerator()
 
 test = ROOT.TPythia6()
 
 P6gen = ROOT.tPythia6Generator()
 
 myPgun = ROOT.FairBoxGenerator(options.pID,1)
 
 DISgen = ROOT.MuDISGenerator()
 
 mu_start
 
 mu_end
 
 nEvents = 0
 
 Nuagegen = ROOT.NuageGenerator()
 
int ntt = 6
 
int nXcells = 7
 
int nYcells = 3
 
int nZcells = ntt -1
 
int startx = -ship_geo.NuTauTarget.xdim/2. + nXcells*ship_geo.NuTauTarget.BrX
 
int endx = -ship_geo.NuTauTarget.xdim/2. + (nXcells+1)*ship_geo.NuTauTarget.BrX
 
int starty = -ship_geo.NuTauTarget.ydim/2. + nYcells*ship_geo.NuTauTarget.BrY
 
int endy = - ship_geo.NuTauTarget.ydim/2. + (nYcells+1)*ship_geo.NuTauTarget.BrY
 
int startz = ship_geo.EmuMagnet.zC - ship_geo.NuTauTarget.zdim/2. + ntt *ship_geo.NuTauTT.TTZ + nZcells * ship_geo.NuTauTarget.CellW
 
int endz = ship_geo.EmuMagnet.zC - ship_geo.NuTauTarget.zdim/2. + ntt *ship_geo.NuTauTT.TTZ + nZcells * ship_geo.NuTauTarget.CellW + ship_geo.NuTauTarget.BrZ
 
 Geniegen = ROOT.GenieGenerator()
 
 pdg = ROOT.TDatabasePDG.Instance()
 
 Ntuplegen = ROOT.NtupleGenerator_FLUKA()
 
 fileType = ut.checkFileExists(inputFile)
 
 MuonBackgen = ROOT.MuonBackGenerator()
 
 testf = ROOT.TFile.Open(test)
 
 Cosmicsgen = ROOT.CosmicsGenerator()
 
 n_EVENTS
 
 fastMuon
 
 gMC = ROOT.TVirtualMC.GetMC()
 
 fStack = gMC.GetStack()
 
 gProcessTable = ROOT.G4ProcessTable.GetProcessTable()
 
 procBrems = gProcessTable.FindProcess(ROOT.G4String('muBrems'),ROOT.G4String('mu+'))
 
 trajFilter = ROOT.FairTrajFilter.Instance()
 
 fieldMaker = geomGeant4.addVMCFields(ship_geo, '', True)
 
 onlyWithField
 
 exclude
 
 mygMC = ROOT.TGeant4.GetMC()
 
 procmu = gProcessTable.FindProcess(ROOT.G4String('muIoni'),ROOT.G4String('mu+'))
 
 kParameterMerged = ROOT.kTRUE
 
 parOut = ROOT.FairParRootFileIo(kParameterMerged)
 
 fGeo = ROOT.gGeoManager
 
 rtime = timer.RealTime()
 
 ctime = timer.CpuTime()
 
str tmpFile = outFile+"tmp"
 
str xxx = outFile.split('/')
 
str check = xxx[len(xxx)-1]
 
bool fin = False
 
 nm = ff.GetName().split('/')
 
bool t = fin.cbmsim
 
 fout = ROOT.TFile(tmpFile,'recreate')
 
bool sTree = t.CloneTree(0)
 
list pointContainers = []
 
 name = x.GetName()
 
bool rc = t.GetEvent(n)
 
bool empty = True
 
 rc1 = os.system("rm "+outFile)
 
 rc2 = os.system("mv "+tmpFile+" "+outFile)
 

Function Documentation

◆ checkOverlapsWithGeant4()

run_simScript.checkOverlapsWithGeant4 ( )

Definition at line 736 of file run_simScript.py.

736def checkOverlapsWithGeant4():
737 # after /run/initialize, but prints warning messages, problems with TGeo volume
738 mygMC = ROOT.TGeant4.GetMC()
739 mygMC.ProcessGeantCommand("/geometry/test/recursion_start 0")
740 mygMC.ProcessGeantCommand("/geometry/test/recursion_depth 2")
741 mygMC.ProcessGeantCommand("/geometry/test/run")

◆ visualizeMagFields()

run_simScript.visualizeMagFields ( )

Definition at line 734 of file run_simScript.py.

734def visualizeMagFields():
735 checkMagFields.run()

Variable Documentation

◆ action

run_simScript.action

Definition at line 62 of file run_simScript.py.

◆ ALPACAgen

run_simScript.ALPACAgen = ROOT.ALPACAGenerator()

Definition at line 354 of file run_simScript.py.

◆ CaloDesign

run_simScript.CaloDesign

Definition at line 202 of file run_simScript.py.

◆ charmonly

bool run_simScript.charmonly = False

Definition at line 46 of file run_simScript.py.

◆ check

str run_simScript.check = xxx[len(xxx)-1]

Definition at line 704 of file run_simScript.py.

◆ checking4overlaps

run_simScript.checking4overlaps = False

Definition at line 57 of file run_simScript.py.

◆ Cosmicsgen

run_simScript.Cosmicsgen = ROOT.CosmicsGenerator()

Definition at line 557 of file run_simScript.py.

◆ ctime

run_simScript.ctime = timer.CpuTime()

Definition at line 687 of file run_simScript.py.

◆ debug

int run_simScript.debug = 0

Definition at line 21 of file run_simScript.py.

◆ default

run_simScript.default = '2018'

Definition at line 54 of file run_simScript.py.

◆ defaultInputFile

bool run_simScript.defaultInputFile = True

Definition at line 50 of file run_simScript.py.

◆ dest

run_simScript.dest

Definition at line 62 of file run_simScript.py.

◆ DISgen

run_simScript.DISgen = ROOT.MuDISGenerator()

Definition at line 393 of file run_simScript.py.

◆ DownScaleDiMuon

bool run_simScript.DownScaleDiMuon = False

Definition at line 25 of file run_simScript.py.

◆ dryrun

bool run_simScript.dryrun = False

Definition at line 23 of file run_simScript.py.

◆ ds

run_simScript.ds

Definition at line 202 of file run_simScript.py.

◆ dy

run_simScript.dy

Definition at line 201 of file run_simScript.py.

◆ empty

run_simScript.empty = True

Definition at line 720 of file run_simScript.py.

◆ endx

int run_simScript.endx = -ship_geo.NuTauTarget.xdim/2. + (nXcells+1)*ship_geo.NuTauTarget.BrX

Definition at line 428 of file run_simScript.py.

◆ endy

int run_simScript.endy = - ship_geo.NuTauTarget.ydim/2. + (nYcells+1)*ship_geo.NuTauTarget.BrY

Definition at line 430 of file run_simScript.py.

◆ endZ

run_simScript.endZ = startZ + lengthZ

Definition at line 346 of file run_simScript.py.

◆ endz

int run_simScript.endz = ship_geo.EmuMagnet.zC - ship_geo.NuTauTarget.zdim/2. + ntt *ship_geo.NuTauTT.TTZ + nZcells * ship_geo.NuTauTarget.CellW + ship_geo.NuTauTarget.BrZ

Definition at line 432 of file run_simScript.py.

◆ exclude

run_simScript.exclude

Definition at line 639 of file run_simScript.py.

◆ False

run_simScript.False

Definition at line 62 of file run_simScript.py.

◆ fastMuon

run_simScript.fastMuon

Definition at line 571 of file run_simScript.py.

◆ fGeo

run_simScript.fGeo = ROOT.gGeoManager

Definition at line 676 of file run_simScript.py.

◆ fieldMaker

run_simScript.fieldMaker = geomGeant4.addVMCFields(ship_geo, '', True)

Definition at line 633 of file run_simScript.py.

◆ fileType

run_simScript.fileType = ut.checkFileExists(inputFile)

Definition at line 527 of file run_simScript.py.

◆ fin

run_simScript.fin = False

Definition at line 705 of file run_simScript.py.

◆ float

run_simScript.float

Definition at line 79 of file run_simScript.py.

◆ fout

run_simScript.fout = ROOT.TFile(tmpFile,'recreate')

Definition at line 711 of file run_simScript.py.

◆ fStack

run_simScript.fStack = gMC.GetStack()

Definition at line 591 of file run_simScript.py.

◆ Geniegen

run_simScript.Geniegen = ROOT.GenieGenerator()

Definition at line 463 of file run_simScript.py.

◆ globalDesigns

dict run_simScript.globalDesigns
Initial value:
1= {'2016':{'dy':10.,'dv':5,'ds':7,'nud':1,'caloDesign':0,'strawDesign':4},\
2 '2018':{'dy':10.,'dv':6,'ds':9,'nud':3,'caloDesign':3,'strawDesign':10}}

Definition at line 52 of file run_simScript.py.

◆ gMC

run_simScript.gMC = ROOT.TVirtualMC.GetMC()

Definition at line 590 of file run_simScript.py.

◆ gProcessTable

run_simScript.gProcessTable = ROOT.G4ProcessTable.GetProcessTable()

Definition at line 609 of file run_simScript.py.

◆ group

run_simScript.group = parser.add_mutually_exclusive_group()

Definition at line 61 of file run_simScript.py.

◆ help

run_simScript.help

Definition at line 62 of file run_simScript.py.

◆ HNL

bool run_simScript.HNL = True

Definition at line 47 of file run_simScript.py.

◆ IgnoreCommandLineOptions

run_simScript.IgnoreCommandLineOptions

Definition at line 13 of file run_simScript.py.

◆ inactivateMuonProcesses

bool run_simScript.inactivateMuonProcesses = False

Definition at line 56 of file run_simScript.py.

◆ inclusive

run_simScript.inclusive = "c"

Definition at line 40 of file run_simScript.py.

◆ inputFile

str run_simScript.inputFile = "/eos/experiment/ship/data/Charm/Cascade-parp16-MSTP82-1-MSEL4-978Bpot.root"

Definition at line 49 of file run_simScript.py.

◆ int

run_simScript.int

Definition at line 122 of file run_simScript.py.

◆ kParameterMerged

run_simScript.kParameterMerged = ROOT.kTRUE

Definition at line 661 of file run_simScript.py.

◆ lengthZ

run_simScript.lengthZ = target.length

Definition at line 345 of file run_simScript.py.

◆ Lmax

tuple run_simScript.Lmax = (ship_geo.TrackStation1.z - ship_geo.target.z0)/100.

Definition at line 349 of file run_simScript.py.

◆ Lmin

tuple run_simScript.Lmin = ((ship_geo.Chamber1.z - ship_geo.chambers.Tub1length) - ship_geo.target.z0)/100.

Definition at line 348 of file run_simScript.py.

◆ mcEngine

str run_simScript.mcEngine = "TGeant4"

Definition at line 37 of file run_simScript.py.

◆ MCTracksWithEnergyCutOnly

bool run_simScript.MCTracksWithEnergyCutOnly = True

Definition at line 43 of file run_simScript.py.

◆ MCTracksWithHitsOnly

bool run_simScript.MCTracksWithHitsOnly = False

Definition at line 42 of file run_simScript.py.

◆ MCTracksWithHitsOrEnergyCut

bool run_simScript.MCTracksWithHitsOrEnergyCut = False

Definition at line 44 of file run_simScript.py.

◆ modules

run_simScript.modules = shipDet_conf.configure(run,ship_geo)

Definition at line 263 of file run_simScript.py.

◆ motherMode

run_simScript.motherMode = True

Definition at line 35 of file run_simScript.py.

◆ mu_end

run_simScript.mu_end

Definition at line 395 of file run_simScript.py.

◆ mu_start

run_simScript.mu_start

Definition at line 395 of file run_simScript.py.

◆ MuonBackgen

run_simScript.MuonBackgen = ROOT.MuonBackGenerator()

Definition at line 534 of file run_simScript.py.

◆ muShieldDesign

run_simScript.muShieldDesign

Definition at line 202 of file run_simScript.py.

◆ muShieldGeo

run_simScript.muShieldGeo

Definition at line 203 of file run_simScript.py.

◆ muShieldStepGeo

run_simScript.muShieldStepGeo

Definition at line 204 of file run_simScript.py.

◆ muShieldWithCobaltMagnet

run_simScript.muShieldWithCobaltMagnet

Definition at line 204 of file run_simScript.py.

◆ mygMC

run_simScript.mygMC = ROOT.TGeant4.GetMC()

Definition at line 646 of file run_simScript.py.

◆ myPgun

run_simScript.myPgun = ROOT.FairBoxGenerator(options.pID,1)

Definition at line 378 of file run_simScript.py.

◆ n_EVENTS

run_simScript.n_EVENTS

Definition at line 563 of file run_simScript.py.

◆ name

run_simScript.name = x.GetName()

Definition at line 716 of file run_simScript.py.

◆ nEvents

int run_simScript.nEvents = 0

Definition at line 409 of file run_simScript.py.

◆ nm

run_simScript.nm = ff.GetName().split('/')

Definition at line 707 of file run_simScript.py.

◆ None

run_simScript.None

Definition at line 87 of file run_simScript.py.

◆ ntt

int run_simScript.ntt = 6

Definition at line 423 of file run_simScript.py.

◆ Ntuplegen

run_simScript.Ntuplegen = ROOT.NtupleGenerator_FLUKA()

Definition at line 509 of file run_simScript.py.

◆ Nuagegen

run_simScript.Nuagegen = ROOT.NuageGenerator()

Definition at line 415 of file run_simScript.py.

◆ nud

run_simScript.nud

Definition at line 202 of file run_simScript.py.

◆ nuTauTargetDesign

run_simScript.nuTauTargetDesign

Definition at line 202 of file run_simScript.py.

◆ nXcells

int run_simScript.nXcells = 7

Definition at line 424 of file run_simScript.py.

◆ nYcells

int run_simScript.nYcells = 3

Definition at line 425 of file run_simScript.py.

◆ nZcells

int run_simScript.nZcells = ntt -1

Definition at line 426 of file run_simScript.py.

◆ onlyWithField

run_simScript.onlyWithField

Definition at line 638 of file run_simScript.py.

◆ Opt_high

run_simScript.Opt_high = int(options.cosmics)

Definition at line 153 of file run_simScript.py.

◆ options

run_simScript.options = parser.parse_args()

Definition at line 127 of file run_simScript.py.

◆ outFile

str run_simScript.outFile = "%s/ship.%s.root" % (options.outputDir, tag)

Definition at line 235 of file run_simScript.py.

◆ P6gen

run_simScript.P6gen = ROOT.tPythia6Generator()

Definition at line 372 of file run_simScript.py.

◆ P8gen

run_simScript.P8gen = ROOT.HNLPythia8Generator()

Definition at line 271 of file run_simScript.py.

◆ parFile

str run_simScript.parFile = "%s/ship.params.%s.root" % (options.outputDir, tag)

Definition at line 241 of file run_simScript.py.

◆ parOut

run_simScript.parOut = ROOT.FairParRootFileIo(kParameterMerged)

Definition at line 662 of file run_simScript.py.

◆ parser

run_simScript.parser = ArgumentParser()

Definition at line 60 of file run_simScript.py.

◆ passDPconf

run_simScript.passDPconf = pythia8darkphoton_conf.configure(P8gen,options.theMass,options.theDPepsilon,inclusive, motherMode, options.deepCopy)

Definition at line 302 of file run_simScript.py.

◆ pdg

run_simScript.pdg = ROOT.TDatabasePDG.Instance()

Definition at line 499 of file run_simScript.py.

◆ pointContainers

list run_simScript.pointContainers = []

Definition at line 714 of file run_simScript.py.

◆ primGen

run_simScript.primGen = ROOT.FairPrimaryGenerator()

Definition at line 265 of file run_simScript.py.

◆ procBrems

run_simScript.procBrems = gProcessTable.FindProcess(ROOT.G4String('muBrems'),ROOT.G4String('mu+'))

Definition at line 610 of file run_simScript.py.

◆ procmu

run_simScript.procmu = gProcessTable.FindProcess(ROOT.G4String('muIoni'),ROOT.G4String('mu+'))

Definition at line 654 of file run_simScript.py.

◆ rc

bool run_simScript.rc = t.GetEvent(n)

Definition at line 719 of file run_simScript.py.

◆ rc1

run_simScript.rc1 = os.system("rm "+outFile)

Definition at line 729 of file run_simScript.py.

◆ rc2

run_simScript.rc2 = os.system("mv "+tmpFile+" "+outFile)

Definition at line 730 of file run_simScript.py.

◆ required

run_simScript.required

Definition at line 62 of file run_simScript.py.

◆ rtdb

run_simScript.rtdb = run.GetRuntimeDb()

Definition at line 255 of file run_simScript.py.

◆ rtime

run_simScript.rtime = timer.RealTime()

Definition at line 686 of file run_simScript.py.

◆ run

run_simScript.run = ROOT.FairRunSim()

Definition at line 251 of file run_simScript.py.

◆ ship_geo

run_simScript.ship_geo = ConfigRegistry.loadpy("$FAIRSHIP/geometry/charm-geometry_config.py", Setup = options.CharmdetSetup, cTarget = options.CharmTarget)

Definition at line 201 of file run_simScript.py.

◆ simEngine

str run_simScript.simEngine = "Pythia8"

Definition at line 38 of file run_simScript.py.

◆ SmearBeam

int run_simScript.SmearBeam = 1*u.cm

Definition at line 347 of file run_simScript.py.

◆ startx

int run_simScript.startx = -ship_geo.NuTauTarget.xdim/2. + nXcells*ship_geo.NuTauTarget.BrX

Definition at line 427 of file run_simScript.py.

◆ starty

int run_simScript.starty = -ship_geo.NuTauTarget.ydim/2. + nYcells*ship_geo.NuTauTarget.BrY

Definition at line 429 of file run_simScript.py.

◆ startZ

run_simScript.startZ = target.z0

Definition at line 344 of file run_simScript.py.

◆ startz

int run_simScript.startz = ship_geo.EmuMagnet.zC - ship_geo.NuTauTarget.zdim/2. + ntt *ship_geo.NuTauTT.TTZ + nZcells * ship_geo.NuTauTarget.CellW

Definition at line 431 of file run_simScript.py.

◆ strawDesign

run_simScript.strawDesign

Definition at line 203 of file run_simScript.py.

◆ sTree

bool run_simScript.sTree = t.CloneTree(0)

Definition at line 712 of file run_simScript.py.

◆ t

bool run_simScript.t = fin.cbmsim

Definition at line 710 of file run_simScript.py.

◆ tag

run_simScript.tag

Definition at line 225 of file run_simScript.py.

◆ tankDesign

run_simScript.tankDesign

Definition at line 201 of file run_simScript.py.

◆ target

run_simScript.target = ship_geo.target

Definition at line 343 of file run_simScript.py.

◆ test

run_simScript.test = ROOT.TPythia6()

Definition at line 371 of file run_simScript.py.

◆ testf

run_simScript.testf = ROOT.TFile.Open(test)

Definition at line 541 of file run_simScript.py.

◆ theCouplings

list run_simScript.theCouplings = [float(c) for c in options.thecouplings.split(",")]

Definition at line 164 of file run_simScript.py.

◆ theDecayCouplings

list run_simScript.theDecayCouplings = [float(c) for c in options.thedeccouplings.split(",")]

Definition at line 168 of file run_simScript.py.

◆ theDPmass

float run_simScript.theDPmass = 0.2*u.GeV

Definition at line 32 of file run_simScript.py.

◆ theHNLMass

float run_simScript.theHNLMass = 1.0*u.GeV

Definition at line 28 of file run_simScript.py.

◆ theMass

run_simScript.theMass

Definition at line 161 of file run_simScript.py.

◆ theProductionCouplings

list run_simScript.theProductionCouplings = None

Definition at line 29 of file run_simScript.py.

◆ timer

run_simScript.timer = ROOT.TStopwatch()

Definition at line 247 of file run_simScript.py.

◆ tmpFile

str run_simScript.tmpFile = outFile+"tmp"

Definition at line 702 of file run_simScript.py.

◆ trajFilter

run_simScript.trajFilter = ROOT.FairTrajFilter.Instance()

Definition at line 616 of file run_simScript.py.

◆ type

run_simScript.type

Definition at line 65 of file run_simScript.py.

◆ xxx

str run_simScript.xxx = outFile.split('/')

Definition at line 703 of file run_simScript.py.

◆ Yheight

run_simScript.Yheight

Definition at line 201 of file run_simScript.py.