"PoolSource",
# replace 'myfile.root' with the source file you want to use
#fileNames = cms.untracked.vstring(
# 'file:/afs/cern.ch/cms/Tutorials/TWIKI_DATA/TTJets_8TeV_53X.root'
#)
fileNames=cms.untracked.vstring(ivars.inputFiles))
# the fragmentation systematics calculator
process.load('SimGeneral.HepPDTESSource.pythiapdt_cfi')
process.mergedGenParticles = cms.EDProducer(
"MergedGenParticleProducer",
inputPruned=cms.InputTag("prunedGenParticles"),
inputPacked=cms.InputTag("packedGenParticles"),
)
from GeneratorInterface.RivetInterface.genParticles2HepMC_cfi import genParticles2HepMC
process.genParticles2HepMC = genParticles2HepMC.clone(
genParticles=cms.InputTag("mergedGenParticles"))
process.load("GeneratorInterface.RivetInterface.particleLevel_cfi")
process.particleLevel.excludeNeutrinosFromJetClustering = False
process.load('TopQuarkAnalysis.BFragmentationAnalyzer.bfragWgtProducer_cfi')
#process.dump=cms.EDAnalyzer('EventContentAnalyzer')
#process.Tracer = cms.Service("Tracer")
### Set output
process.TFileService = cms.Service("TFileService",
fileName=cms.string(ivars.outputFile))
# this is loaded in ../../RecoMET/METFilters/test/test_badChargedCandidateFilter.py
# before met filters
process.load('Configuration.StandardSequences.Services_cff')
# .. and I still get MissingParameter: Parameter 'BadChargedCandidateFilter' not found.
)
)
)
if options.frag=='P':
process.generator.PythiaParameters.processParameters.append('StringZ:usePetersonB = on')
print "Running Peterson with epsilonB = %s" % options.param
process.generator.PythiaParameters.processParameters.append('StringZ:epsilonB = %f'%options.param)
if options.frag=='BL':
process.generator.PythiaParameters.processParameters.append('StringZ:rFactB = %f'%options.param)
print "Running Bowler-Lund with rFactB = %s" % options.param
print "max events = %s" % options.maxEvents
#pseudo-top config
from GeneratorInterface.RivetInterface.genParticles2HepMC_cfi import genParticles2HepMC
process.genParticles2HepMC = genParticles2HepMC.clone( genParticles = cms.InputTag("genParticles") )
process.load("TopQuarkAnalysis.TopEventProducers.producers.pseudoTop_cfi")
process.pseudoTop.jetMaxEta=cms.double(5.0)
#analysis config
process.TFileService = cms.Service("TFileService",
fileName = cms.string(options.outputFile)
)
process.load('TopLJets2015.TopAnalysis.bfragAnalysis_cfi')
# Path and EndPath definitions
process.ProductionFilterSequence = cms.Sequence(process.generator)
process.generation_step = cms.Path(process.pgen)
process.AnalysisSequence = cms.Path(process.genParticles2HepMC*process.pseudoTop*process.bfragAnalysis)
process.genfiltersummary_step = cms.EndPath(process.genFilterSummary)
process.endjob_step = cms.EndPath(process.endOfProcess)