P8gen.SetParameters("ProcessLevel:all = off")
if ds==7: # short muon shield
P8gen.SetLmin(44*u.m)
P8gen.SetLmax(107*u.m)
if inputFile:
ut.checkFileExists(inputFile)
# read from external file
P8gen.UseExternalFile(inputFile, firstEvent)
if DarkPhoton:
P8gen = ROOT.DPPythia8Generator()
if inclusive=='qcd':
P8gen.SetDPId(4900023)
else:
P8gen.SetDPId(9900015)
import pythia8darkphoton_conf
passDPconf = pythia8darkphoton_conf.configure(P8gen,theDPmass,theDPepsilon,inclusive,deepCopy)
if (passDPconf!=1): sys.exit()
P8gen.SetSmearBeam(1*u.cm) # finite beam size
if ds==7: # short muon shield
P8gen.SetLmin(44*u.m)
P8gen.SetLmax(107*u.m)
if charmonly:
ut.checkFileExists(inputFile)
primGen.SetBeam(0.,0., ship_geo.Box.TX-2., ship_geo.Box.TY-2.) #Uniform distribution in x/y on the target (1 cm of margin at both sides)
primGen.SmearVertexXY(True)
P8gen = ROOT.Pythia8Generator()
P8gen.UseExternalFile(inputFile, firstEvent)
if ship_geo.MufluxSpectrometer.muflux == False :
P8gen.SetTarget("volTarget_1",0.,0.) # will distribute PV inside target, beam offset x=y=0.
else:
print "ERROR: charmonly option should not be used for the muonflux measurement"
theCouplings[2], RPVSUSYbench,
inclusive, deepCopy)
P8gen.SetParameters("ProcessLevel:all = off")
if inputFile:
ut.checkFileExists(inputFile)
# read from external file
P8gen.UseExternalFile(inputFile, firstEvent)
if DarkPhoton:
P8gen = ROOT.DPPythia8Generator()
if inclusive == 'qcd':
P8gen.SetDPId(4900023)
else:
P8gen.SetDPId(9900015)
import pythia8darkphoton_conf
passDPconf = pythia8darkphoton_conf.configure(P8gen, theDPmass,
theDPepsilon, inclusive,
deepCopy)
if (passDPconf != 1): sys.exit()
if HNL or RPVSUSY or DarkPhoton:
P8gen.SetSmearBeam(1 * u.cm) # finite beam size
P8gen.SetLmin((ship_geo.Chamber1.z - ship_geo.chambers.Tub1length) -
ship_geo.target.z0)
P8gen.SetLmax(ship_geo.TrackStation1.z - ship_geo.target.z0)
if charmonly:
primGen.SetTarget(0., 0.) #vertex is setted in pythia8Generator
ut.checkFileExists(inputFile)
if ship_geo.Box.gausbeam:
primGen.SetBeam(
0., 0., 0.5,
0.5) #more central beam, for hits in downstream detectors
primGen.SmearGausVertexXY(True) #sigma = x
options.RPVSUSYbench, inclusive,
options.deepCopy)
P8gen.SetParameters("ProcessLevel:all = off")
if inputFile:
ut.checkFileExists(inputFile)
# read from external file
P8gen.UseExternalFile(inputFile, options.firstEvent)
if options.DarkPhoton:
P8gen = ROOT.DPPythia8Generator()
if inclusive == 'qcd':
P8gen.SetDPId(4900023)
else:
P8gen.SetDPId(9900015)
import pythia8darkphoton_conf
passDPconf = pythia8darkphoton_conf.configure(P8gen, options.theMass,
options.theDPepsilon,
inclusive, motherMode,
options.deepCopy)
if (passDPconf != 1): sys.exit()
if HNL or options.RPVSUSY or options.DarkPhoton:
P8gen.SetSmearBeam(1 * u.cm) # finite beam size
P8gen.SetLmin((ship_geo.Chamber1.z - ship_geo.chambers.Tub1length) -
ship_geo.target.z0)
P8gen.SetLmax(ship_geo.TrackStation1.z - ship_geo.target.z0)
if charmonly:
primGen.SetTarget(0., 0.) #vertex is setted in pythia8Generator
ut.checkFileExists(inputFile)
if ship_geo.Box.gausbeam:
primGen.SetBeam(
0., 0., 0.5,
0.5) #more central beam, for hits in downstream detectors
primGen.SmearGausVertexXY(True) #sigma = x
