Updating jobOptions files and adding initial state protons to EventFiller

Generators/Superchic_i/python/EventFiller.py

@@ -49,6 +49,34 @@ class LheEVNTFiller(EvgenAlg):
                 if event and firstLine:
                     firstLine = False
                     evt.weights().push_back(float(line.split()[2]))
+
+                    #Add the initial state protons
+                    pos = HepMC.FourVector(0.0, 0.0, 0.0, 0.0)
+                    gv = HepMC.GenVertex(pos)
+                    ROOT.SetOwnership(gv, False)
+                    evt.add_vertex(gv)
+                    mom = HepMC.FourVector( 0. , 0. , 6499999.9323 , 6500000.00 )
+                    gp = HepMC.GenParticle()
+                    gp.set_status( 2 )
+                    gp.set_pdg_id( 2212 )
+                    gp.set_momentum(mom)
+                    gp.set_generated_mass(938.272046)
+                    ROOT.SetOwnership(gp, False)
+                    gv.add_particle_out(gp)
+
+                    pos = HepMC.FourVector(0.0, 0.0, 0.0, 0.0)
+                    gv = HepMC.GenVertex(pos)
+                    ROOT.SetOwnership(gv, False)
+                    evt.add_vertex(gv)
+                    mom = HepMC.FourVector( 0. , 0. , -6499999.9323 , 6500000.00 )
+                    gp = HepMC.GenParticle()
+                    gp.set_status( 2 )
+                    gp.set_pdg_id( 2212 )
+                    gp.set_momentum(mom)
+                    gp.set_generated_mass(938.272046)
+                    ROOT.SetOwnership(gp, False)
+                    gv.add_particle_out(gp)
+
                     continue
                 if event:
                     pos = HepMC.FourVector(0.0, 0.0, 0.0, 0.0)
@@ -60,7 +88,6 @@ class LheEVNTFiller(EvgenAlg):
                     gp.set_status(int(line.split()[1]) )
                     gp.set_pdg_id(int(line.split()[0]) )
                     gp.set_momentum(mom)
-                    print(float(line.split()[10]) * 1000.)
                     gp.set_generated_mass(float(line.split()[10]) * 1000.)
                     ROOT.SetOwnership(gp, False)
                     gv.add_particle_out(gp)

Generators/Superchic_i/share/jobOptions.SuperchicPy8EG_yyMuMu_DD.py → Generators/Superchic_i/share/jobOptions.SuperChicPy8_yyMuMu_DD.py

@@ -35,54 +35,48 @@ Init.genunw  = ".true."
 hasNev = hasattr(runArgs,"maxEvents")
 if hasNev:
     Init.nev = str(runArgs.maxEvents)
-Init.ymin  = "-2.4d0"              # Minimum object rapidity Y_X
-Init.ymax  = "2.4d0"               # Maximum object rapidity Y_X
-Init.mmin  = "6d0"                 # Minimum object mass M_X
-Init.mmax  = "500d0"               # Maximum object mass M_X
+Init.ymin  = "-5.0d0"              # Minimum object rapidity Y_X
+Init.ymax  = "5.0d0"               # Maximum object rapidity Y_X
+Init.mmin  = "20d0"                 # Minimum object mass M_X
+Init.mmax  = "2000d0"               # Maximum object mass M_X
 Init.gencuts  = ".true."           # Generate cuts below
 Init.ptxmax  = "100d0"         #cut on proton pt
-Init.ptamin  = "3.0d0"             # Minimum pT of outgoing object a 
-Init.ptbmin  = "3.0d0"             # Minimum pT of outgoing object b 
-Init.etaamin  = "-2.4d0"           # Minimum eta of outgoing object a
-Init.etaamax   = "2.4d0"           # Maximum eta of outgoing object a
-Init.etabmin  = "-2.4d0"           # Minimum eta of outgoing object b
-Init.etabmax   = "2.4d0"           # Maximum eta of outgoing object b
+Init.ptamin  = "10.0d0"             # Minimum pT of outgoing object a 
+Init.ptbmin  = "10.0d0"             # Minimum pT of outgoing object b 
+Init.etaamin  = "-2.5d0"           # Minimum eta of outgoing object a
+Init.etaamax   = "2.5d0"           # Maximum eta of outgoing object a
+Init.etabmin  = "-2.5d0"           # Minimum eta of outgoing object b
+Init.etabmax   = "2.5d0"           # Maximum eta of outgoing object b
 Init.acoabmax  = "100d0"
 
 
 SuperChicRun(Init)
 
 
+## Base config for Pythia8
+from Pythia8_i.Pythia8_iConf import Pythia8_i
+genSeq += Pythia8_i("Pythia8")
 evgenConfig.generators += ["Pythia8"]
 
-include("Pythia8_i/Pythia8_Base_Fragment.py")
+## Control storing LHE in the HepMC record
+if "StoreLHE" in genSeq.Pythia8.__slots__.keys():
+    print "Pythia8_Base_Fragment.py: DISABLING storage of LHE record in HepMC by default. Please re-enable storage if desired"
+    genSeq.Pythia8.StoreLHE = False
 
 genSeq.Pythia8.Commands += [
-  "PDF:pSet= LHAPDF6:NNPDF23_nnlo_as_0118_qed",
-    "SpaceShower:rapidityOrder = on",
-    "SigmaProcess:alphaSvalue = 0.140",
-    "SpaceShower:pT0Ref = 1.56",
-    "SpaceShower:pTmaxFudge = 0.91",
-    "SpaceShower:pTdampFudge = 1.05",
-    "SpaceShower:alphaSvalue = 0.127",
-    "TimeShower:alphaSvalue = 0.127",
-    "SpaceShower:pTdampMatch = 1",
-    "BeamRemnants:primordialKThard = 1.88",
     "PartonLevel:MPI = off",
-    "PartonLevel:FSR = on",
     "SpaceShower:dipoleRecoil = on",
     "SpaceShower:pTmaxMatch = 2",
-    "SpaceShower:QEDshowerByQ = off"
+    "BeamRemnants:primordialKT = off",
+    "SpaceShower:QEDshowerByQ = off",
+    "BeamRemnants:unresolvedHadron = 0",
+    "SpaceShower:pTdampMatch=1"
 ]
 
-evgenConfig.tune = "NNPDF23_QED"
-
-#EvtGen for b fragmentation as default.  No EvtGen is available in "nonStandard"
-include("Pythia8_i/Pythia8_EvtGen.py")
-
-include("Pythia8_i/Pythia8_ShowerWeights.py")
-
 assert hasattr(genSeq, "Pythia8")
 fileName = "evrecs/evrec"+Init.outtg[1:-1]+".dat"
 genSeq.Pythia8.LHEFile = fileName
 genSeq.Pythia8.CollisionEnergy = int(runArgs.ecmEnergy)
+
+testSeq.TestHepMC.MaxTransVtxDisp = 1000000
+testSeq.TestHepMC.MaxVtxDisp      = 1000000000

Generators/Superchic_i/share/jobOptions.SuperChicPy8_yyMuMu_El.py

+evgenConfig.description = "SuperChic4 MC gamma + gamma pp collisions at 13000 GeV to 2 muons with dissociation on both side"
+evgenConfig.keywords = ["2photon","2muon","dissociation"]
+#evgenConfig.weighting = 0
+evgenConfig.contact = ["william.patrick.mccormack.iii@cern.ch"]
+
+if not os.path.exists('inputs'):
+    os.makedirs('inputs')
+if not os.path.exists('evrecs'):
+    os.makedirs('evrecs')
+if not os.path.exists('outputs'):
+    os.makedirs('outputs')
+
+from Superchic_i.SuperChicUtils import *
+
+#class with the superchic initialization parameters.  Please see SuperChicUtils for a complete list of tunable parameters.
+Init = SuperChicConfig()
+
+hasECM = hasattr(runArgs,"ecmEnergy")
+if hasECM:
+    ecm = int(runArgs.ecmEnergy)
+    Init.rts = str(ecm)+"d0" # set the COM collision energy (in fortran syntax)
+Init.isurv = "4"                  # Model of soft survival (from 1 -> 4, corresponding to arXiv:1306.2149)
+Init.intag = "'in5'"              # for input files
+Init.PDFname = "'MMHT2014lo68cl'" # PDF set name
+Init.PDFmember = "0"              # PDF member
+Init.proc = "57"                  # Process number (59 = gg->gg, 56 = gg->ee, 68 = gg->a->gg ); Please consult Superchic Manual https://superchic.hepforge.org/
+Init.beam = "'prot'"              # Beam type ('prot', 'ion')
+Init.outtg = "'out'"              # for output file name
+Init.sfaci = ".true."             # Include soft survival effects
+Init.diff = "'el'"                # interaction: elastic ('el'), single ('sd','sda','sdb') and double ('dd') dissociation.
+hasSeed = hasattr(runArgs,"randomSeed")
+if hasSeed:
+    Init.iseed  = str(int(runArgs.randomSeed))
+Init.genunw  = ".true."
+hasNev = hasattr(runArgs,"maxEvents")
+if hasNev:
+    Init.nev = str(runArgs.maxEvents)
+Init.ymin  = "-5.0d0"              # Minimum object rapidity Y_X
+Init.ymax  = "5.0d0"               # Maximum object rapidity Y_X
+Init.mmin  = "20d0"                 # Minimum object mass M_X
+Init.mmax  = "2000d0"               # Maximum object mass M_X
+Init.gencuts  = ".true."           # Generate cuts below
+Init.ptxmax  = "100d0"         #cut on proton pt
+Init.ptamin  = "10.0d0"             # Minimum pT of outgoing object a 
+Init.ptbmin  = "10.0d0"             # Minimum pT of outgoing object b 
+Init.etaamin  = "-2.5d0"           # Minimum eta of outgoing object a
+Init.etaamax   = "2.5d0"           # Maximum eta of outgoing object a
+Init.etabmin  = "-2.5d0"           # Minimum eta of outgoing object b
+Init.etabmax   = "2.5d0"           # Maximum eta of outgoing object b
+Init.acoabmax  = "100d0"
+
+
+SuperChicRun(Init)
+
+
+## Base config for Pythia8
+from Pythia8_i.Pythia8_iConf import Pythia8_i
+genSeq += Pythia8_i("Pythia8")
+evgenConfig.generators += ["Pythia8"]
+
+## Control storing LHE in the HepMC record
+if "StoreLHE" in genSeq.Pythia8.__slots__.keys():
+    print "Pythia8_Base_Fragment.py: DISABLING storage of LHE record in HepMC by default. Please re-enable storage if desired"
+    genSeq.Pythia8.StoreLHE = False
+
+genSeq.Pythia8.Commands += [
+    "PartonLevel:MPI = off",
+    "SpaceShower:pTmaxMatch = 2",
+    "BeamRemnants:primordialKT = off",
+    "BeamRemnants:unresolvedHadron = 3",
+    "SpaceShower:pTdampMatch=1",
+    "PartonLevel:ISR = off",
+    "LesHouches:matchInOut = off"
+]
+
+assert hasattr(genSeq, "Pythia8")
+fileName = "evrecs/evrec"+Init.outtg[1:-1]+".dat"
+genSeq.Pythia8.LHEFile = fileName
+genSeq.Pythia8.CollisionEnergy = int(runArgs.ecmEnergy)
+
+testSeq.TestHepMC.MaxTransVtxDisp = 1000000
+testSeq.TestHepMC.MaxVtxDisp      = 1000000000

Generators/Superchic_i/share/jobOptions.SuperChicPy8_yyMuMu_SDA.py

+evgenConfig.description = "SuperChic4 MC gamma + gamma pp collisions at 13000 GeV to 2 muons with dissociation on both side"
+evgenConfig.keywords = ["2photon","2muon","dissociation"]
+#evgenConfig.weighting = 0
+evgenConfig.contact = ["william.patrick.mccormack.iii@cern.ch"]
+
+if not os.path.exists('inputs'):
+    os.makedirs('inputs')
+if not os.path.exists('evrecs'):
+    os.makedirs('evrecs')
+if not os.path.exists('outputs'):
+    os.makedirs('outputs')
+
+from Superchic_i.SuperChicUtils import *
+
+#class with the superchic initialization parameters.  Please see SuperChicUtils for a complete list of tunable parameters.
+Init = SuperChicConfig()
+
+hasECM = hasattr(runArgs,"ecmEnergy")
+if hasECM:
+    ecm = int(runArgs.ecmEnergy)
+    Init.rts = str(ecm)+"d0" # set the COM collision energy (in fortran syntax)
+Init.isurv = "4"                  # Model of soft survival (from 1 -> 4, corresponding to arXiv:1306.2149)
+Init.intag = "'in5'"              # for input files
+Init.PDFname = "'MMHT2014lo68cl'" # PDF set name
+Init.PDFmember = "0"              # PDF member
+Init.proc = "57"                  # Process number (59 = gg->gg, 56 = gg->ee, 68 = gg->a->gg ); Please consult Superchic Manual https://superchic.hepforge.org/
+Init.beam = "'prot'"              # Beam type ('prot', 'ion')
+Init.outtg = "'out'"              # for output file name
+Init.sfaci = ".true."             # Include soft survival effects
+Init.diff = "'sda'"                # interaction: elastic ('el'), single ('sd','sda','sdb') and double ('dd') dissociation.
+hasSeed = hasattr(runArgs,"randomSeed")
+if hasSeed:
+    Init.iseed  = str(int(runArgs.randomSeed))
+Init.genunw  = ".true."
+hasNev = hasattr(runArgs,"maxEvents")
+if hasNev:
+    Init.nev = str(runArgs.maxEvents)
+Init.ymin  = "-5.0d0"              # Minimum object rapidity Y_X
+Init.ymax  = "5.0d0"               # Maximum object rapidity Y_X
+Init.mmin  = "20d0"                 # Minimum object mass M_X
+Init.mmax  = "2000d0"               # Maximum object mass M_X
+Init.gencuts  = ".true."           # Generate cuts below
+Init.ptxmax  = "100d0"         #cut on proton pt
+Init.ptamin  = "10.0d0"             # Minimum pT of outgoing object a 
+Init.ptbmin  = "10.0d0"             # Minimum pT of outgoing object b 
+Init.etaamin  = "-2.5d0"           # Minimum eta of outgoing object a
+Init.etaamax   = "2.5d0"           # Maximum eta of outgoing object a
+Init.etabmin  = "-2.5d0"           # Minimum eta of outgoing object b
+Init.etabmax   = "2.5d0"           # Maximum eta of outgoing object b
+Init.acoabmax  = "100d0"
+
+
+SuperChicRun(Init)
+
+
+## Base config for Pythia8
+from Pythia8_i.Pythia8_iConf import Pythia8_i
+genSeq += Pythia8_i("Pythia8")
+evgenConfig.generators += ["Pythia8"]
+
+## Control storing LHE in the HepMC record
+if "StoreLHE" in genSeq.Pythia8.__slots__.keys():
+    print "Pythia8_Base_Fragment.py: DISABLING storage of LHE record in HepMC by default. Please re-enable storage if desired"
+    genSeq.Pythia8.StoreLHE = False
+
+genSeq.Pythia8.Commands += [
+    "PartonLevel:MPI = off",
+    "SpaceShower:dipoleRecoil = on",
+    "SpaceShower:pTmaxMatch = 2",
+    "BeamRemnants:primordialKT = off",
+    "SpaceShower:QEDshowerByQ = off",
+    "BeamRemnants:unresolvedHadron = 2",
+    "SpaceShower:pTdampMatch=1"
+]
+
+assert hasattr(genSeq, "Pythia8")
+fileName = "evrecs/evrec"+Init.outtg[1:-1]+".dat"
+genSeq.Pythia8.LHEFile = fileName
+genSeq.Pythia8.CollisionEnergy = int(runArgs.ecmEnergy)
+
+testSeq.TestHepMC.MaxTransVtxDisp = 1000000
+testSeq.TestHepMC.MaxVtxDisp      = 1000000000

Generators/Superchic_i/share/jobOptions.SuperChicPy8_yyMuMu_SDB.py

+evgenConfig.description = "SuperChic4 MC gamma + gamma pp collisions at 13000 GeV to 2 muons with dissociation on both side"
+evgenConfig.keywords = ["2photon","2muon","dissociation"]
+#evgenConfig.weighting = 0
+evgenConfig.contact = ["william.patrick.mccormack.iii@cern.ch"]
+
+if not os.path.exists('inputs'):
+    os.makedirs('inputs')
+if not os.path.exists('evrecs'):
+    os.makedirs('evrecs')
+if not os.path.exists('outputs'):
+    os.makedirs('outputs')
+
+from Superchic_i.SuperChicUtils import *
+
+#class with the superchic initialization parameters.  Please see SuperChicUtils for a complete list of tunable parameters.
+Init = SuperChicConfig()
+
+hasECM = hasattr(runArgs,"ecmEnergy")
+if hasECM:
+    ecm = int(runArgs.ecmEnergy)
+    Init.rts = str(ecm)+"d0" # set the COM collision energy (in fortran syntax)
+Init.isurv = "4"                  # Model of soft survival (from 1 -> 4, corresponding to arXiv:1306.2149)
+Init.intag = "'in5'"              # for input files
+Init.PDFname = "'MMHT2014lo68cl'" # PDF set name
+Init.PDFmember = "0"              # PDF member
+Init.proc = "57"                  # Process number (59 = gg->gg, 56 = gg->ee, 68 = gg->a->gg ); Please consult Superchic Manual https://superchic.hepforge.org/
+Init.beam = "'prot'"              # Beam type ('prot', 'ion')
+Init.outtg = "'out'"              # for output file name
+Init.sfaci = ".true."             # Include soft survival effects
+Init.diff = "'sdb'"                # interaction: elastic ('el'), single ('sd','sda','sdb') and double ('dd') dissociation.
+hasSeed = hasattr(runArgs,"randomSeed")
+if hasSeed:
+    Init.iseed  = str(int(runArgs.randomSeed))
+Init.genunw  = ".true."
+hasNev = hasattr(runArgs,"maxEvents")
+if hasNev:
+    Init.nev = str(runArgs.maxEvents)
+Init.ymin  = "-5.0d0"              # Minimum object rapidity Y_X
+Init.ymax  = "5.0d0"               # Maximum object rapidity Y_X
+Init.mmin  = "20d0"                 # Minimum object mass M_X
+Init.mmax  = "2000d0"               # Maximum object mass M_X
+Init.gencuts  = ".true."           # Generate cuts below
+Init.ptxmax  = "100d0"         #cut on proton pt
+Init.ptamin  = "10.0d0"             # Minimum pT of outgoing object a 
+Init.ptbmin  = "10.0d0"             # Minimum pT of outgoing object b 
+Init.etaamin  = "-2.5d0"           # Minimum eta of outgoing object a
+Init.etaamax   = "2.5d0"           # Maximum eta of outgoing object a
+Init.etabmin  = "-2.5d0"           # Minimum eta of outgoing object b
+Init.etabmax   = "2.5d0"           # Maximum eta of outgoing object b
+Init.acoabmax  = "100d0"
+
+
+SuperChicRun(Init)
+
+
+## Base config for Pythia8
+from Pythia8_i.Pythia8_iConf import Pythia8_i
+genSeq += Pythia8_i("Pythia8")
+evgenConfig.generators += ["Pythia8"]
+
+## Control storing LHE in the HepMC record
+if "StoreLHE" in genSeq.Pythia8.__slots__.keys():
+    print "Pythia8_Base_Fragment.py: DISABLING storage of LHE record in HepMC by default. Please re-enable storage if desired"
+    genSeq.Pythia8.StoreLHE = False
+
+genSeq.Pythia8.Commands += [
+    "PartonLevel:MPI = off",
+    "SpaceShower:dipoleRecoil = on",
+    "SpaceShower:pTmaxMatch = 2",
+    "BeamRemnants:primordialKT = off",
+    "SpaceShower:QEDshowerByQ = off",
+    "BeamRemnants:unresolvedHadron = 1",
+    "SpaceShower:pTdampMatch=1"
+]
+
+assert hasattr(genSeq, "Pythia8")
+fileName = "evrecs/evrec"+Init.outtg[1:-1]+".dat"
+genSeq.Pythia8.LHEFile = fileName
+genSeq.Pythia8.CollisionEnergy = int(runArgs.ecmEnergy)
+
+testSeq.TestHepMC.MaxTransVtxDisp = 1000000
+testSeq.TestHepMC.MaxVtxDisp      = 1000000000