elif theProductionCouplings is not None and theDecayCouplings is not None:
print 'and with couplings', theProductionCouplings, 'at production'
print 'and', theDecayCouplings, 'at decay'
else:
raise ValueError(
'Either both production and decay couplings must be specified, or neither.'
)
pythia8_conf.configure(P8gen, theMass, theProductionCouplings,
theDecayCouplings, inclusive, deepCopy)
if RPVSUSY:
print 'Generating RPVSUSY events of mass %.3f GeV' % theMass
print 'and with couplings=[%.3f,%.3f]' % (theCouplings[0],
theCouplings[1])
print 'and with stop mass=%.3f GeV\n' % theCouplings[2]
pythia8_conf.configurerpvsusy(P8gen, theMass,
[theCouplings[0], theCouplings[1]],
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,
primGen = ROOT.FairPrimaryGenerator()
if simEngine == "Pythia8":
primGen.SetTarget(ship_geo.target.z0, 0.)
# -----Pythia8--------------------------------------
if HNL or RPVSUSY:
P8gen = ROOT.HNLPythia8Generator()
import pythia8_conf
if HNL:
print 'Generating HNL events of mass %.3f GeV\n'%theMass
print 'and with couplings=',theCouplings
pythia8_conf.configure(P8gen,theMass,theCouplings,inclusive,deepCopy)
if RPVSUSY:
print 'Generating RPVSUSY events of mass %.3f GeV\n'%theMass
print 'and with couplings=[%.3f,%.3f]\n'%(theCouplings[0],theCouplings[1])
print 'and with stop mass=\%.3f GeV\n',theCouplings[2]
pythia8_conf.configurerpvsusy(P8gen,theMass,[theCouplings[0],theCouplings[1]],
theCouplings[2],RPVSUSYbench,'c',deepCopy)
P8gen.SetSmearBeam(1*u.cm) # finite beam size
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)
