'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,
deepCopy)
if (passDPconf != 1): sys.exit()
if HNL or RPVSUSY or DarkPhoton:
P8gen.SetSmearBeam(1 * u.cm) # finite beam size
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)
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)
theProductionCouplings = theDecayCouplings = theCouplings
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,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 )
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)
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 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)
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
