/afs/cern.ch/work/a/aaltunda/public/www/CMSSW_5_3_14/src/DPGAnalysis/Skims/python/singleElectronSkim_cff.py

Go to the documentation of this file.

import FWCore.ParameterSet.Config as cms

#from HLTrigger.HLTfilters.hltHighLevel_cfi import *
#exoticaMuHLT = hltHighLevel
#Define the HLT path to be used.
#exoticaMuHLT.HLTPaths =['HLT_L1MuOpen']
#exoticaMuHLT.TriggerResultsTag = cms.InputTag("TriggerResults","","HLT8E29")

#Define the HLT quality cut 
#exoticaHLTMuonFilter = cms.EDFilter("HLTSummaryFilter",
#    summary = cms.InputTag("hltTriggerSummaryAOD","","HLT8E29"), # trigger summary
#    member  = cms.InputTag("hltL3MuonCandidates","","HLT8E29"),      # filter or collection                                                                    
#    cut     = cms.string("pt>0"),                     # cut on trigger object
#    minN    = cms.int32(0)                  # min. # of passing objects needed
# )
                               

#Define the Reco quality cut

singleRecoElectronPt20Filter = cms.EDFilter("GsfElectronRefSelector",
                                        src = cms.InputTag("gsfElectrons"),
                                        cut = cms.string('pt > 20.0 && abs(eta) < 2.5 ' ),
                                        filter = cms.bool(True),
                                        minN    = cms.int32(1) 
                                        )

singleRecoElectronPt15Filter = cms.EDFilter("GsfElectronRefSelector",
                                        src = cms.InputTag("gsfElectrons"),
                                        cut = cms.string('pt > 15.0 && abs(eta) < 2.5 ' ),
                                        filter = cms.bool(True),
                                        minN    = cms.int32(1)
                                        )

singleRecoElectronPt10Filter = cms.EDFilter("GsfElectronRefSelector",
                                        src = cms.InputTag("gsfElectrons"),
                                        cut = cms.string('pt > 10.0 && abs(eta) < 2.5 ' ),
                                        filter = cms.bool(True),
                                        minN    = cms.int32(1)
                                        )

singleRecoElectronPt5Filter = cms.EDFilter("GsfElectronRefSelector",
                                        src = cms.InputTag("gsfElectrons"),
                                        cut = cms.string('pt > 5.0 && abs(eta) < 2.5 ' ),
                                        filter = cms.bool(True),
                                        minN    = cms.int32(1)
                                        )

singleRecoElectronPt1Filter = cms.EDFilter("GsfElectronRefSelector",
                                           src = cms.InputTag("gsfElectrons"),
                                           cut = cms.string('pt > 1.0 && abs(eta) < 2.5 ' ),
                                           filter = cms.bool(True),
                                           minN    = cms.int32(1)
                                           )


#superClusterMerger =  cms.EDFilter("EgammaSuperClusterMerger",
#                                   src = cms.VInputTag(cms.InputTag('correctedHybridSuperClusters'),
#                                                       cms.InputTag('correctedMulti5x5SuperClustersWithPreshower'))
#                                   )
#superClusterCands = cms.EDProducer("ConcreteEcalCandidateProducer",
#                                   src = cms.InputTag("superClusterMerger"),
#                                   particleType = cms.string('e-')
#                                   )

#goodSuperClusters = cms.EDFilter("CandViewRefSelector",
#                                 src = cms.InputTag("superClusterCands"),
#                                 cut = cms.string('et > 1.0')
#                                 )

#superClusterPt5Filter = cms.EDFilter("CandViewCountFilter",
#                                      src = cms.InputTag("goodSuperClusters"),
#                                      minNumber = cms.uint32(1)
#                                      )

#electronCounter = cms.EDFilter("CandViewCountFilter",
#                               src = cms.InputTag("gsfElectrons"),
#                               minNumber = cms.uint32(1)
#                               )

#from RecoEgamma.EgammaElectronProducers.ecalDrivenElectronSeeds_cfi import *

#ecalDrivenElectronSeeds.SCEtCut = cms.double(1.0)
#ecalDrivenElectronSeeds.applyHOverECut = cms.bool(False)
#ecalDrivenElectronSeeds.SeedConfiguration.SeedConfigurationnSigmasDeltaZ1 = cms.double(5.)
#ecalDrivenElectronSeeds.SeedConfiguration.z2MinB = cms.double(-0.9)
#ecalDrivenElectronSeeds.SeedConfiguration.z2MaxB = cms.double(0.9)
#ecalDrivenElectronSeeds.SeedConfiguration.r2MinF = cms.double(-1.5)
#ecalDrivenElectronSeeds.SeedConfiguration.r2MaxF = cms.double(1.5)
#ecalDrivenElectronSeeds.SeedConfiguration.rMinI = cms.double(-2.)
#ecalDrivenElectronSeeds.SeedConfiguration.rMaxI = cms.double(2.)
#ecalDrivenElectronSeeds.SeedConfiguration.DeltaPhi1Low = cms.double(0.3)
#ecalDrivenElectronSeeds.SeedConfiguration.DeltaPhi1High = cms.double(0.3)
#ecalDrivenElectronSeeds.SeedConfiguration.DeltaPhi2 = cms.double(0.3)
#gsfElectrons.applyPreselection = cms.bool(False)
#gsfElectrons.applyAmbResolution = cms.bool(True)

#Define group sequence, using HLT/Reco quality cut. 
#singleMuHLTQualitySeq = cms.Sequence()
singleElectronPt20RecoQualitySeq = cms.Sequence(
    #singleElectronHLT+
    singleRecoElectronPt20Filter
)

singleElectronPt15RecoQualitySeq = cms.Sequence(
    #singleElectronHLT+
    singleRecoElectronPt15Filter
)

singleElectronPt10RecoQualitySeq = cms.Sequence(
    #singleElectronHLT+
    singleRecoElectronPt10Filter
)

singleElectronPt5RecoQualitySeq = cms.Sequence(
        #singleElectronHLT+
        singleRecoElectronPt5Filter
        )

singleElectronPt1RecoQualitySeq = cms.Sequence(
            #singleElectronHLT+
            singleRecoElectronPt1Filter
                    )

#singleElectronSCRecoQualitySeq = cms.Sequence(
#    singleRecoElectronPt5Filter
#    )