EgammaValidationAutoConf_cff.py

Go to the documentation of this file.

import FWCore.ParameterSet.Config as cms

class dummy:
    pass

samples=dummy()

##########################################################
# Define which preselections to run                      #
##########################################################

samples.names = ['Wenu',
                 'Zee',
                 #'TripleEle',
                 'GammaJet',
                 'DiGamma']
samples.pdgid = [ 11,
                  11,
                  #11,
                  22,
                  22]
samples.num   = [1,
                 2,
                 #3,
                 1,
                 2]

##########################################################
# produce generated paricles in acceptance               #
##########################################################

genp = cms.EDFilter("PdgIdAndStatusCandViewSelector",
    status = cms.vint32(3),
    src = cms.InputTag("genParticles"),
    pdgId = cms.vint32(11)  # replaced in loop
)

fiducial = cms.EDFilter("EtaPtMinCandViewSelector",
    src = cms.InputTag("genp"),
    etaMin = cms.double(-2.5),  # to be replaced in loop ?
    etaMax = cms.double(2.5),   # to be replaced in loop ?
    ptMin = cms.double(2.0)     # to be replaced in loop ?
)

##########################################################
# loop over samples to create modules and sequence       #
##########################################################

tmp = cms.SequencePlaceholder("tmp")
egammaSelectors = cms.Sequence(tmp) # no empty sequences allowed, start with dummy

#loop over samples
for samplenum in range(len(samples.names)):

    # clone genparticles and select correct type
    genpartname = "genpart"+samples.names[samplenum]
    globals()[genpartname] = genp.clone()
    setattr(globals()[genpartname],"pdgId",cms.vint32(samples.pdgid[samplenum]) ) # set pdgId
    egammaSelectors *= globals()[genpartname]                            # add to sequence

    # clone generator fiducial region
    fiducialname = "fiducial"+samples.names[samplenum]
    globals()[fiducialname] = fiducial.clone()
    setattr(globals()[fiducialname],"src",cms.InputTag(genpartname) ) # set input collection
    egammaSelectors *= globals()[fiducialname]               # add to sequence

egammaSelectors.remove(tmp)  # remove the initial dummy

dqmFeeder = cms.EDAnalyzer('EmDQMFeeder',
                           #processname = cms.string("HLT"), # can be obtained from triggerobject
                           triggerobject = cms.InputTag("hltTriggerSummaryRAW","","HLT"),
                           genEtaAcc = cms.double(2.5),
                           genEtAcc = cms.double(2.0),
                           PtMax = cms.untracked.double(100.0),
                           isData = cms.bool(False),
                           verbosity = cms.untracked.uint32(0),
                           mcMatchedOnly = cms.untracked.bool(True),
                           noPhiPlots = cms.untracked.bool(True),
                           noIsolationPlots = cms.untracked.bool(True)
                          )

# selectors go into separate "prevalidation" sequence
egammaValidationSequence   = cms.Sequence(dqmFeeder)
egammaValidationSequenceFS = cms.Sequence(dqmFeeder)