Skip to content

Commit

Permalink
Create new RelVal workflows for ExternalGeneratorFilter (corrected)
Browse files Browse the repository at this point in the history
  • Loading branch information
colizz committed Jul 31, 2021
1 parent 2c7b50e commit 932a8f2
Show file tree
Hide file tree
Showing 29 changed files with 1,370 additions and 1 deletion.
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,66 @@
import FWCore.ParameterSet.Config as cms

from Configuration.Generator.Pythia8CommonSettings_cfi import *
from Configuration.Generator.Pythia8CUEP8M1Settings_cfi import *
from GeneratorInterface.EvtGenInterface.EvtGenSetting_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.38e-3),
crossSection = cms.untracked.double(540000000.),
comEnergy = cms.double(13000.0),

ExternalDecays = cms.PSet(
EvtGen130 = cms.untracked.PSet(
decay_table = cms.string('GeneratorInterface/EvtGenInterface/data/DECAY_2010.DEC'),
particle_property_file = cms.FileInPath('GeneratorInterface/EvtGenInterface/data/evt_Bmm.pdl'),
user_decay_file = cms.vstring('GeneratorInterface/EvtGenInterface/data/Bd_mumu.dec'),
list_forced_decays = cms.vstring('MyB0',
'Myanti-B0'),
operates_on_particles = cms.vint32(),
),
parameterSets = cms.vstring('EvtGen130')
),

PythiaParameters = cms.PSet(pythia8CommonSettingsBlock,
pythia8CUEP8M1SettingsBlock,
processParameters = cms.vstring("SoftQCD:nonDiffractive = on"),
parameterSets = cms.vstring('pythia8CommonSettings',
'pythia8CUEP8M1Settings',
'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/PYTHIA8_Bd2MuMu_EtaPtFilter_CUEP8M1_13TeV_cff.py $'),
annotation = cms.untracked.string('Spring 2015: Pythia8+EvtGen130 generation of B0 --> Mu+Mu-, 13TeV, Tune CUETP8M1')
)

bfilter = cms.EDFilter(
"PythiaFilter",
MaxEta = cms.untracked.double(9999.),
MinEta = cms.untracked.double(-9999.),
ParticleID = cms.untracked.int32(511)
)

decayfilter = cms.EDFilter(
"PythiaDauVFilter",
verbose = cms.untracked.int32(1),
NumberDaughters = cms.untracked.int32(2),
ParticleID = cms.untracked.int32(511),
DaughterIDs = cms.untracked.vint32(13, -13),
MinPt = cms.untracked.vdouble(2.5, 2.5),
MinEta = cms.untracked.vdouble(-2.5, -2.5),
MaxEta = cms.untracked.vdouble( 2.5, 2.5)
)

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

from Configuration.Generator.Pythia8CommonSettings_cfi import *
from Configuration.Generator.Pythia8CUEP8M1Settings_cfi import *
from GeneratorInterface.EvtGenInterface.EvtGenSetting_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.38e-3),
crossSection = cms.untracked.double(540000000.),
comEnergy = cms.double(13000.0),

ExternalDecays = cms.PSet(
EvtGen130 = cms.untracked.PSet(
decay_table = cms.string('GeneratorInterface/EvtGenInterface/data/DECAY_2010.DEC'),
particle_property_file = cms.FileInPath('GeneratorInterface/EvtGenInterface/data/evt_Bmm.pdl'),
user_decay_file = cms.vstring('GeneratorInterface/EvtGenInterface/data/Bs_JpsiPhi_V3.dec'),
list_forced_decays = cms.vstring('MyB_s0',
'Myanti-B_s0'),
),
operates_on_particles = cms.vint32(),
parameterSets = cms.vstring('EvtGen130')
),

PythiaParameters = cms.PSet(pythia8CommonSettingsBlock,
pythia8CUEP8M1SettingsBlock,
processParameters = cms.vstring("SoftQCD:nonDiffractive = on"),
parameterSets = cms.vstring('pythia8CommonSettings',
'pythia8CUEP8M1Settings',
'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/PYTHIA8_BsJpsiPhi_EtaPtFilter_CUEP8M1_13TeV_cff.py $'),
annotation = cms.untracked.string('Spring 2015: Pythia8+EvtGen130 generation of Bs --> J/psi phi, 13TeV, Tune CUETP8M1')
)

bfilter = cms.EDFilter(
"PythiaFilter",
MaxEta = cms.untracked.double(9999.),
MinEta = cms.untracked.double(-9999.),
ParticleID = cms.untracked.int32(531)
)

jpsifilter = cms.EDFilter(
"PythiaDauVFilter",
verbose = cms.untracked.int32(1),
NumberDaughters = cms.untracked.int32(2),
MotherID = cms.untracked.int32(531),
ParticleID = cms.untracked.int32(443),
DaughterIDs = cms.untracked.vint32(13, -13),
MinPt = cms.untracked.vdouble(2.5, 2.5),
MinEta = cms.untracked.vdouble(-2.5, -2.5),
MaxEta = cms.untracked.vdouble( 2.5, 2.5)
)

phifilter = cms.EDFilter(
"PythiaDauVFilter",
verbose = cms.untracked.int32(1),
NumberDaughters = cms.untracked.int32(2),
MotherID = cms.untracked.int32(531),
ParticleID = cms.untracked.int32(333),
DaughterIDs = cms.untracked.vint32(321, -321),
MinPt = cms.untracked.vdouble(0.4, 0.4),
MinEta = cms.untracked.vdouble(-2.5, -2.5),
MaxEta = cms.untracked.vdouble( 2.5, 2.5)
)

ProductionFilterSequence = cms.Sequence(generator*bfilter*jpsifilter*phifilter)

Original file line number Diff line number Diff line change
@@ -0,0 +1,66 @@
import FWCore.ParameterSet.Config as cms

from Configuration.Generator.Pythia8CommonSettings_cfi import *
from Configuration.Generator.Pythia8CUEP8M1Settings_cfi import *
from GeneratorInterface.EvtGenInterface.EvtGenSetting_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.38e-3),
crossSection = cms.untracked.double(540000000.),
comEnergy = cms.double(13000.0),

ExternalDecays = cms.PSet(
EvtGen130 = cms.untracked.PSet(
decay_table = cms.string('GeneratorInterface/EvtGenInterface/data/DECAY_2010.DEC'),
particle_property_file = cms.FileInPath('GeneratorInterface/EvtGenInterface/data/evt_Bsmm.pdl'),
user_decay_file = cms.vstring('GeneratorInterface/EvtGenInterface/data/Bs_mumu.dec'),
list_forced_decays = cms.vstring('MyB_s0',
'Myanti-B_s0'),
operates_on_particles = cms.vint32(),
),
parameterSets = cms.vstring('EvtGen130')
),

PythiaParameters = cms.PSet(pythia8CommonSettingsBlock,
pythia8CUEP8M1SettingsBlock,
processParameters = cms.vstring("SoftQCD:nonDiffractive = on"),
parameterSets = cms.vstring('pythia8CommonSettings',
'pythia8CUEP8M1Settings',
'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/PYTHIA8_Bs2MuMu_EtaPtFilter_CUEP8M1_13TeV_cff.py $'),
annotation = cms.untracked.string('Spring 2015: Pythia8+EvtGen130 generation of Bs --> Mu+Mu-, 13TeV, Tune CUETP8M1')
)

bfilter = cms.EDFilter(
"PythiaFilter",
MaxEta = cms.untracked.double(9999.),
MinEta = cms.untracked.double(-9999.),
ParticleID = cms.untracked.int32(531)
)

decayfilter = cms.EDFilter(
"PythiaDauVFilter",
verbose = cms.untracked.int32(1),
NumberDaughters = cms.untracked.int32(2),
ParticleID = cms.untracked.int32(531),
DaughterIDs = cms.untracked.vint32(13, -13),
MinPt = cms.untracked.vdouble(2.5, 2.5),
MinEta = cms.untracked.vdouble(-2.5, -2.5),
MaxEta = cms.untracked.vdouble( 2.5, 2.5)
)

ProductionFilterSequence = cms.Sequence(generator*bfilter*decayfilter)
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)
Loading

0 comments on commit 932a8f2

Please sign in to comment.