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