Commit 0b2b20c3 authored by Andrzej Olszewski's avatar Andrzej Olszewski Committed by Jan Kretzschmar
Browse files

Py8 gamgam->jets (wmin20) new DSIDs

parent 34030242
Pipeline #2949471 passed with stages
in 1 minute and 13 seconds
evgenConfig.description = "Pythia8 jet gamma gamma -> qqbar, resolved x resolved, using nuclear photon flux and NCTEQ15 PDFs with A14 tune"
evgenConfig.keywords = ["QED","QCD","coherent","jets"]
evgenConfig.contact = ["angerami@cern.ch"]
evgenConfig.nEventsPerJob = 10000
include("Pythia8_i/Pythia8_A14_NNPDF23LO_EvtGen_Common.py")
from Pythia8_i.Pythia8_iConf import UserPhotonFlux
photonFluxTool=UserPhotonFlux()
photonFluxTool.Process=3
photonFluxTool.NuclearCharge=82;
photonFluxTool.MinimumB=6.62;
from AthenaCommon.AppMgr import ToolSvc
ToolSvc += UserPhotonFlux()
genSeq.Pythia8.CustomInterface = photonFluxTool
#options to set up gamma* p collisions
genSeq.Pythia8.Beam1='ELECTRON'
genSeq.Pythia8.Beam2='ELECTRON'
genSeq.Pythia8.Commands += ["Photon:sampleQ2 = off",
"PDF:lepton2gamma = on",
"PDF:lepton2gammaSet = 2",
"PDF:lepton2gammaApprox = 2",
"PhaseSpace:pTHatMin = 5.0",
"Photon:Wmin = 20.0"] #Wmin = 2 x pThatmin
#resolved x resolved
genSeq.Pythia8.Commands +=["Photon:ProcessType = 1",
"HardQCD:all = on"]
if hasattr(testSeq, "TestHepMC"):
testSeq.TestHepMC.TachyonsTest=False
include("GeneratorFilters/FindJets.py")
CreateJets(prefiltSeq, 0.4)
AddJetsFilter(filtSeq,runArgs.ecmEnergy, 0.4)
from AthenaCommon.SystemOfUnits import GeV
filtSeq.QCDTruthJetFilter.MinPt = 10.0*GeV
evgenConfig.description = "Pythia8 jet gamma gamma -> qqbar, resolved x direct, using nuclear photon flux and NCTEQ15 PDFs with A14 tune"
evgenConfig.keywords = ["QED","QCD","coherent","jets"]
evgenConfig.contact = ["angerami@cern.ch"]
evgenConfig.nEventsPerJob = 10000
include("Pythia8_i/Pythia8_A14_NNPDF23LO_EvtGen_Common.py")
from Pythia8_i.Pythia8_iConf import UserPhotonFlux
photonFluxTool=UserPhotonFlux()
photonFluxTool.Process=3
photonFluxTool.NuclearCharge=82;
photonFluxTool.MinimumB=6.62;
from AthenaCommon.AppMgr import ToolSvc
ToolSvc += UserPhotonFlux()
genSeq.Pythia8.CustomInterface = photonFluxTool
#options to set up gamma* p collisions
genSeq.Pythia8.Beam1='ELECTRON'
genSeq.Pythia8.Beam2='ELECTRON'
genSeq.Pythia8.Commands += ["Photon:sampleQ2 = off",
"PDF:lepton2gamma = on",
"PDF:lepton2gammaSet = 2",
"PDF:lepton2gammaApprox = 2",
"PhaseSpace:pTHatMin = 5.0",
"Photon:Wmin = 20.0"] #Wmin = 2 x pThatmin
#resolved x direct
genSeq.Pythia8.Commands +=["Photon:ProcessType = 2",
"PhotonParton:all = on"]
if hasattr(testSeq, "TestHepMC"):
testSeq.TestHepMC.TachyonsTest=False
include("GeneratorFilters/FindJets.py")
CreateJets(prefiltSeq, 0.4)
AddJetsFilter(filtSeq,runArgs.ecmEnergy, 0.4)
from AthenaCommon.SystemOfUnits import GeV
filtSeq.QCDTruthJetFilter.MinPt = 10.0*GeV
evgenConfig.description = "Pythia8 jet gamma gamma -> qqbar, direct x resolved, using nuclear photon flux and NCTEQ15 PDFs with A14 tune"
evgenConfig.keywords = ["QED","QCD","coherent","jets"]
evgenConfig.contact = ["angerami@cern.ch"]
evgenConfig.nEventsPerJob = 10000
include("Pythia8_i/Pythia8_A14_NNPDF23LO_EvtGen_Common.py")
from Pythia8_i.Pythia8_iConf import UserPhotonFlux
photonFluxTool=UserPhotonFlux()
photonFluxTool.Process=3
photonFluxTool.NuclearCharge=82;
photonFluxTool.MinimumB=6.62;
from AthenaCommon.AppMgr import ToolSvc
ToolSvc += UserPhotonFlux()
genSeq.Pythia8.CustomInterface = photonFluxTool
#options to set up gamma* p collisions
genSeq.Pythia8.Beam1='ELECTRON'
genSeq.Pythia8.Beam2='ELECTRON'
genSeq.Pythia8.Commands += ["Photon:sampleQ2 = off",
"PDF:lepton2gamma = on",
"PDF:lepton2gammaSet = 2",
"PDF:lepton2gammaApprox = 2",
"PhaseSpace:pTHatMin = 5.0",
"Photon:Wmin = 20.0"] #Wmin = 2 x pThatmin
#direct x resolved
genSeq.Pythia8.Commands +=["Photon:ProcessType = 3",
"PhotonParton:all = on"]
if hasattr(testSeq, "TestHepMC"):
testSeq.TestHepMC.TachyonsTest=False
include("GeneratorFilters/FindJets.py")
CreateJets(prefiltSeq, 0.4)
AddJetsFilter(filtSeq,runArgs.ecmEnergy, 0.4)
from AthenaCommon.SystemOfUnits import GeV
filtSeq.QCDTruthJetFilter.MinPt = 10.0*GeV
evgenConfig.description = "Pythia8 jet gamma gamma -> qqbar, direct x direct, using nuclear photon flux and NCTEQ15 PDFs with A14 tune"
evgenConfig.keywords = ["QED","QCD","coherent","jets"]
evgenConfig.contact = ["angerami@cern.ch"]
evgenConfig.nEventsPerJob = 10000
include("Pythia8_i/Pythia8_A14_NNPDF23LO_EvtGen_Common.py")
from Pythia8_i.Pythia8_iConf import UserPhotonFlux
photonFluxTool=UserPhotonFlux()
photonFluxTool.Process=3
photonFluxTool.NuclearCharge=82;
photonFluxTool.MinimumB=6.62;
from AthenaCommon.AppMgr import ToolSvc
ToolSvc += UserPhotonFlux()
genSeq.Pythia8.CustomInterface = photonFluxTool
#options to set up gamma* p collisions
genSeq.Pythia8.Beam1='ELECTRON'
genSeq.Pythia8.Beam2='ELECTRON'
genSeq.Pythia8.Commands += ["Photon:sampleQ2 = off",
"PDF:lepton2gamma = on",
"PDF:lepton2gammaSet = 2",
"PDF:lepton2gammaApprox = 2",
"PhaseSpace:pTHatMin = 5.0",
"Photon:Wmin = 20.0"] #Wmin = 2 x pThatmin
#direct x direct
genSeq.Pythia8.Commands +=["Photon:ProcessType = 4",
"PhotonCollision:gmgm2qqbar = on",
"PhotonCollision:gmgm2ccbar = on",
"PhotonCollision:gmgm2bbbar = on"]
if hasattr(testSeq, "TestHepMC"):
testSeq.TestHepMC.TachyonsTest=False
include("GeneratorFilters/FindJets.py")
CreateJets(prefiltSeq, 0.4)
AddJetsFilter(filtSeq,runArgs.ecmEnergy, 0.4)
from AthenaCommon.SystemOfUnits import GeV
filtSeq.QCDTruthJetFilter.MinPt = 10.0*GeV
Markdown is supported
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment