Skip to content

Commit

Permalink
Create new RelVal workflows for ExternalGeneratorFilter in relval_ext…
Browse files Browse the repository at this point in the history
…endedgen
  • Loading branch information
colizz committed Aug 2, 2021
1 parent 2c7b50e commit 337cbfb
Show file tree
Hide file tree
Showing 10 changed files with 441 additions and 0 deletions.
Original file line number Diff line number Diff line change
@@ -0,0 +1,57 @@
import FWCore.ParameterSet.Config as cms

_generator = cms.EDFilter("AMPTGeneratorFilter",
diquarky = cms.double(0.0),
diquarkx = cms.double(0.0),
diquarkpx = cms.double(7.0),
ntmax = cms.int32(1000),
dpcoal = cms.double(1000000.0),
diquarkembedding = cms.int32(0),
maxmiss = cms.int32(1000),
ktkick = cms.int32(1),
mu = cms.double(3.2264),
quenchingpar = cms.double(2.0),
popcornpar = cms.double(1.0),
drcoal = cms.double(1000000.0),
amptmode = cms.int32(1),
izpc = cms.int32(0),
popcornmode = cms.bool(True),
minijetpt = cms.double(-7.0),
ks0decay = cms.bool(False),
alpha = cms.double(0.47140452),
dt = cms.double(0.2),
rotateEventPlane = cms.bool(True),
shadowingmode = cms.bool(True),
diquarkpy = cms.double(0.0),
deuteronfactor = cms.int32(5),
stringFragB = cms.double(0.9),#default value in Hijing. Good for pA
quenchingmode = cms.bool(False),
stringFragA = cms.double(0.5),
deuteronmode = cms.int32(0),
doInitialAndFinalRadiation = cms.int32(3),
phidecay = cms.bool(True),
deuteronxsec = cms.int32(1),
pthard = cms.double(2.0),
firstRun = cms.untracked.uint32(1),
frame = cms.string('CMS'),
targ = cms.string('P'),
izp = cms.int32(82),
bMin = cms.double(0),
firstEvent = cms.untracked.uint32(1),
izt = cms.int32(1),
proj = cms.string('A'),
comEnergy = cms.double(5020.0),
iat = cms.int32(1),
bMax = cms.double(15),
iap = cms.int32(208)
)

from GeneratorInterface.Core.ExternalGeneratorFilter import ExternalGeneratorFilter
generator = ExternalGeneratorFilter(_generator)

configurationMetadata = cms.untracked.PSet(
version = cms.untracked.string('$Revision$'),
name = cms.untracked.string('$Source$'),
annotation = cms.untracked.string('AMPT PPb 5020 GeV Minimum Bias')
)

Original file line number Diff line number Diff line change
@@ -0,0 +1,68 @@
import FWCore.ParameterSet.Config as cms
from Configuration.Generator.Pythia8CommonSettings_cfi import *
from Configuration.Generator.MCTunes2017.PythiaCP5Settings_cfi import *
from GeneratorInterface.EvtGenInterface.EvtGenSetting_cff import *

_generator = cms.EDFilter("Pythia8GeneratorFilter",
comEnergy = cms.double(13000.0),
crossSection = cms.untracked.double(54000000000),
filterEfficiency = cms.untracked.double(3.0e-4),
pythiaHepMCVerbosity = cms.untracked.bool(False),
maxEventsToPrint = cms.untracked.int32(0),
pythiaPylistVerbosity = cms.untracked.int32(0),
ExternalDecays = cms.PSet(
#using alternative name for decayer
EvtGen1 = cms.untracked.PSet(
#uses latest evt and decay tables from evtgen
decay_table = cms.string('GeneratorInterface/EvtGenInterface/data/DECAY_NOLONGLIFE.DEC'),
particle_property_file = cms.FileInPath('GeneratorInterface/EvtGenInterface/data/evt.pdl'),
convertPythiaCodes = cms.untracked.bool(False),
#here we will use the user.dec store in the release
user_decay_file = cms.vstring('GeneratorInterface/ExternalDecays/data/Bs_mumu.dec'),
list_forced_decays = cms.vstring('MyB_s0','Myanti-B_s0'),
operates_on_particles = cms.vint32()
),
parameterSets = cms.vstring('EvtGen1')
),
PythiaParameters = cms.PSet(
pythia8CommonSettingsBlock,
pythia8CP5SettingsBlock,
processParameters = cms.vstring(
#filter of a b-quark before hadronizing, and use a better data-like process
'PTFilter:filter = on',
'PTFilter:quarkToFilter = 5',
'PTFilter:scaleToFilter = 1.0',
'SoftQCD:nonDiffractive = on',
),
parameterSets = cms.vstring('pythia8CommonSettings',
'pythia8CP5Settings',
'processParameters',
)
)
)

from GeneratorInterface.Core.ExternalGeneratorFilter import ExternalGeneratorFilter
generator = ExternalGeneratorFilter(_generator)

generator.PythiaParameters.processParameters.extend(EvtGenExtraParticles)

MuMuFilter = cms.EDFilter("MCParticlePairFilter",
Status = cms.untracked.vint32(1, 1),
MinPt = cms.untracked.vdouble(3., 3.),
MaxEta = cms.untracked.vdouble(2.5, 2.5),
MinEta = cms.untracked.vdouble(-2.5, -2.5),
ParticleCharge = cms.untracked.int32(-1),
ParticleID1 = cms.untracked.vint32(13,-13),
)

# -- Require Muon from Bs
MuFilter = cms.EDFilter("PythiaFilter",
Status = cms.untracked.int32(1),
MotherID = cms.untracked.int32(531),
MinPt = cms.untracked.double(3.),
ParticleID = cms.untracked.int32(13),
MaxEta = cms.untracked.double(2.5),
MinEta = cms.untracked.double(-2.5)
)

ProductionFilterSequence = cms.Sequence(generator*MuMuFilter*MuFilter)
Original file line number Diff line number Diff line change
@@ -0,0 +1,106 @@
import FWCore.ParameterSet.Config as cms
from Configuration.Generator.Pythia8CommonSettings_cfi import *
from Configuration.Generator.MCTunes2017.PythiaCP5Settings_cfi import *
from GeneratorInterface.EvtGenInterface.EvtGenSetting_cff import *

_generator = cms.EDFilter("Pythia8GeneratorFilter",
pythiaPylistVerbosity = cms.untracked.int32(0),
pythiaHepMCVerbosity = cms.untracked.bool(False),
comEnergy = cms.double(13000.0),
maxEventsToPrint = cms.untracked.int32(0),
ExternalDecays = cms.PSet(
EvtGen130 = cms.untracked.PSet(
#uses latest evt and decay tables from evtgen
decay_table = cms.string('GeneratorInterface/EvtGenInterface/data/DECAY_NOLONGLIFE.DEC'),
particle_property_file = cms.FileInPath('GeneratorInterface/EvtGenInterface/data/evt.pdl'),
convertPythiaCodes = cms.untracked.bool(False),
#user_decay_file = cms.vstring('GeneratorInterface/ExternalDecays/data/Bu_Kstarmumu_Kspi.dec'),
#content was dump in the embed string below. This should test this feature.
list_forced_decays = cms.vstring('MyB+','MyB-'),
operates_on_particles = cms.vint32(),
user_decay_embedded= cms.vstring(
'# This is the decay file for the decay B+ -> MU+ MU- K*+(-> Ks pi+)',
'#',
'# Descriptor: [B+ -> mu+ mu- {,gamma} {,gamma} (K*+ -> Ks pi+)]cc',
'#',
'# NickName:',
'#',
'# Physics: Includes radiative mode',
'#',
'# Tested: Yes',
'# By: K. Ulmer',
'# Date: 2-26-08',
'#',
'Alias MyB+ B+',
'Alias MyB- B-',
'ChargeConj MyB+ MyB-',
'Alias MyK*+ K*+',
'Alias MyK*- K*-',
'ChargeConj MyK*+ MyK*-',
'Alias MyK_S0 K_S0',
'ChargeConj MyK_S0 MyK_S0',
'#',
'Decay MyB+',
' 1.000 MyK*+ mu+ mu- BTOSLLBALL;',
'Enddecay',
'CDecay MyB-',
'#',
'Decay MyK*+',
' 1.000 MyK_S0 pi+ VSS;',
'Enddecay',
'CDecay MyK*-',
'#',
'Decay MyK_S0',
' 1.000 pi+ pi- PHSP;',
'Enddecay',
'End'
),
),
parameterSets = cms.vstring('EvtGen130')
),
PythiaParameters = cms.PSet(
pythia8CommonSettingsBlock,
pythia8CP5SettingsBlock,
processParameters = cms.vstring(
'HardQCD:gg2bbbar = on ',
'HardQCD:qqbar2bbbar = on ',
'HardQCD:hardbbbar = on',
'PhaseSpace:pTHatMin = 20.',
),
parameterSets = cms.vstring('pythia8CommonSettings',
'pythia8CP5Settings',
'processParameters',
)
)
)

from GeneratorInterface.Core.ExternalGeneratorFilter import ExternalGeneratorFilter
generator = ExternalGeneratorFilter(_generator)

generator.PythiaParameters.processParameters.extend(EvtGenExtraParticles)

configurationMetadata = cms.untracked.PSet(
version = cms.untracked.string('$Revision: 1.1 $'),
name = cms.untracked.string('$Source: Configuration/Generator/python/BuToKstarMuMu_forSTEAM_13TeV_cfi.py $'),
annotation = cms.untracked.string('Summer14: Pythia8+EvtGen130 generation of Bu --> K* Mu+Mu-, 13TeV, Tune CP5')
)

###########
# Filters #
###########
# Filter only pp events which produce a B+:
bufilter = cms.EDFilter("PythiaFilter", ParticleID = cms.untracked.int32(521))

# Filter on final state muons
mumugenfilter = cms.EDFilter("MCParticlePairFilter",
Status = cms.untracked.vint32(1, 1),
MinPt = cms.untracked.vdouble(2.8, 2.8),
MinP = cms.untracked.vdouble(2.8, 2.8),
MaxEta = cms.untracked.vdouble(2.3, 2.3),
MinEta = cms.untracked.vdouble(-2.3, -2.3),
ParticleID1 = cms.untracked.vint32(13,-13),
ParticleID2 = cms.untracked.vint32(13,-13)
)


ProductionFilterSequence = cms.Sequence(generator*bufilter*mumugenfilter)
Original file line number Diff line number Diff line change
@@ -0,0 +1,31 @@
import FWCore.ParameterSet.Config as cms

from Configuration.Generator.Pyquen2015Settings_cff import *

_generator = cms.EDFilter("HydjetGeneratorFilter",
collisionParameters5020GeV,
qgpParameters,
hydjetParameters,
hydjetMode = cms.string('kHydroQJets'),
PythiaParameters = cms.PSet(pyquenPythiaDefaultBlock,
# Quarkonia and Weak Bosons added back upon dilepton group's request.
parameterSets = cms.vstring('pythiaUESettings',
'hydjetPythiaDefault',
'myParameters',
'pythiaJets',
'pythiaPromptPhotons',
'pythiaZjets',
'pythiaBottomoniumNRQCD',
'pythiaCharmoniumNRQCD',
'pythiaQuarkoniaSettings',
'pythiaWeakBosons'
)
),
cFlag = cms.int32(0),
bMin = cms.double(0),
bMax = cms.double(30),
bFixed = cms.double(12)
)

from GeneratorInterface.Core.ExternalGeneratorFilter import ExternalGeneratorFilter
generator = ExternalGeneratorFilter(_generator)
Original file line number Diff line number Diff line change
@@ -0,0 +1,36 @@
import FWCore.ParameterSet.Config as cms

from Configuration.Generator.PyquenDefaultSettings_cff import *

_generator = cms.EDFilter("PyquenGeneratorFilter",
collisionParameters,
qgpParameters,
pyquenParameters,
doQuench = cms.bool(True),
bFixed = cms.double(0.0), ## fixed impact param (fm); valid only if cflag_=0
PythiaParameters = cms.PSet(pyquenPythiaDefaultBlock,
parameterSets = cms.vstring('pythiaUESettings','customProcesses','pythiaZjets','pythiaZtoMuonsAndElectrons','kinematics'),
kinematics = cms.vstring ("CKIN(3)=10", #min pthat
"CKIN(4)=9999", #max pthat
"CKIN(7)=-2.", #min rapidity
"CKIN(8)=2." #max rapidity
)

),
cFlag = cms.int32(0), ## centrality flag
bMin = cms.double(0.0), ## min impact param (fm); valid only if cflag_!=0
bMax = cms.double(0.0) ## max impact param (fm); valid only if cflag_!=0
)

from GeneratorInterface.Core.ExternalGeneratorFilter import ExternalGeneratorFilter
generator = ExternalGeneratorFilter(_generator)

_generator.embeddingMode = False

configurationMetadata = cms.untracked.PSet(
version = cms.untracked.string('$Revision: 1.3 $'),
name = cms.untracked.string('$Source: /local/projects/CMSSW/rep/CMSSW/Configuration/Generator/python/Pyquen_ZeemumuJets_pt10_2760GeV_cfi.py,v $'),
annotation = cms.untracked.string('PYQUEN Z->mumu and Z->ee (pt-hat > 10 GeV) at sqrt(s) = 2.76TeV')
)

ProductionFilterSequence = cms.Sequence(generator)
Original file line number Diff line number Diff line change
@@ -0,0 +1,17 @@
import FWCore.ParameterSet.Config as cms

_generator = cms.EDFilter("ReggeGribovPartonMCGeneratorFilter",

bmin = cms.double(0), #impact parameter min in fm
bmax = cms.double(10000),#impact parameter max in fm
paramFileName = cms.untracked.string("Configuration/Generator/data/ReggeGribovPartonMC.param"), #file with more parameters specific to crmc interface
skipNuclFrag = cms.bool(True), #in HI collisions nuclear fragments with pt=0 can be in the hep event. to skip those activate this option
beammomentum = cms.double(4080),
targetmomentum = cms.double(-4080),
beamid = cms.int32(208),
targetid = cms.int32(1),
model = cms.int32(0),
)

from GeneratorInterface.Core.ExternalGeneratorFilter import ExternalGeneratorFilter
generator = ExternalGeneratorFilter(_generator)
Original file line number Diff line number Diff line change
@@ -0,0 +1,40 @@
import FWCore.ParameterSet.Config as cms
from Configuration.Generator.Pythia8CommonSettings_cfi import *
from Configuration.Generator.Pythia8CUEP8M1Settings_cfi import *
from GeneratorInterface.ExternalDecays.TauolaSettings_cff import *

_generator = cms.EDFilter("Pythia8GeneratorFilter",
pythiaHepMCVerbosity = cms.untracked.bool(False),
maxEventsToPrint = cms.untracked.int32(0),
pythiaPylistVerbosity = cms.untracked.int32(0),
filterEfficiency = cms.untracked.double(1.0),
comEnergy = cms.double(8000.0),
ExternalDecays = cms.PSet(
Tauola = cms.untracked.PSet(
TauolaPolar,
InputCards = cms.PSet(
pjak1 = cms.int32(0),
pjak2 = cms.int32(0),
mdtau = cms.int32(214)
),
),
parameterSets = cms.vstring('Tauola')
),
PythiaParameters = cms.PSet(
pythia8CommonSettingsBlock,
pythia8CUEP8M1SettingsBlock,
processParameters = cms.vstring(
'WeakSingleBoson:ffbar2gmZ = on',
'23:onMode = off',
'23:onIfAny = 15',
'PhaseSpace:mHatMin = 75.',
),
parameterSets = cms.vstring('pythia8CommonSettings',
'pythia8CUEP8M1Settings',
'processParameters',
)
)
)

from GeneratorInterface.Core.ExternalGeneratorFilter import ExternalGeneratorFilter
generator = ExternalGeneratorFilter(_generator, _external_process_components_=cms.vstring("HepPDTESSource"))
Loading

0 comments on commit 337cbfb

Please sign in to comment.