'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