Skims/python/singlePhotonSkim_cff.py

0001 import FWCore.ParameterSet.Config as cms
0002
0003 #from HLTrigger.HLTfilters.hltHighLevel_cfi import *

0004 #exoticaMuHLT = hltHighLevel

0005 #Define the HLT path to be used.

0006 #exoticaMuHLT.HLTPaths =['HLT_L1MuOpen']

0007 #exoticaMuHLT.TriggerResultsTag = cms.InputTag("TriggerResults","","HLT8E29")

0008
0009 #Define the HLT quality cut

0010 #exoticaHLTMuonFilter = cms.EDFilter("HLTSummaryFilter",

0011 #    summary = cms.InputTag("hltTriggerSummaryAOD","","HLT8E29"), # trigger summary

0012 #    member  = cms.InputTag("hltL3MuonCandidates","","HLT8E29"),      # filter or collection

0013 #    cut     = cms.string("pt>0"),                     # cut on trigger object

0014 #    minN    = cms.int32(0)                  # min. # of passing objects needed

0015 # )

0016
0017 #Define the Reco quality cut

0018 singlePhotonPt20Filter = cms.EDFilter("PhotonSelector",
0019                                      src = cms.InputTag("photons"),
0020                                      cut = cms.string('pt > 20 && abs(eta) < 2.7' ),
0021                                      filter = cms.bool(True)
0022 )
0023 singlePhotonPt15Filter = cms.EDFilter("PhotonSelector",
0024                                       src = cms.InputTag("photons"),
0025                                       cut = cms.string('pt > 15 && abs(eta) < 2.7' ),
0026                                       filter = cms.bool(True)
0027                                       )
0028 singlePhotonPt10Filter = cms.EDFilter("PhotonSelector",
0029                                       src = cms.InputTag("photons"),
0030                                       cut = cms.string('pt > 10 && abs(eta) < 2.7' ),
0031                                       filter = cms.bool(True)
0032                                       )
0033 singlePhotonPt5Filter = cms.EDFilter("PhotonSelector",
0034                                      src = cms.InputTag("photons"),
0035                                      cut = cms.string('pt > 5 && abs(eta) < 2.7' ),
0036                                      filter = cms.bool(True)
0037                                      )
0038
0039
0040 egSuperClusterMerger =  cms.EDProducer("EgammaSuperClusterMerger",
0041                                    src = cms.VInputTag(cms.InputTag('correctedHybridSuperClusters'),
0042                                                        cms.InputTag('correctedMulti5x5SuperClustersWithPreshower'))
0043                                    )
0044 egSuperClusterCands = cms.EDProducer("ConcreteEcalCandidateProducer",
0045                                    src = cms.InputTag("egSuperClusterMerger"),
0046                                    particleType = cms.string('gamma')
0047                                    )
0048
0049 eggoodSuperClusters = cms.EDFilter("CandViewRefSelector",
0050                                  src = cms.InputTag("egSuperClusterCands"),
0051                                  cut = cms.string('et > 3.0')
0052                                  )
0053
0054 egSuperClusterPt5Filter = cms.EDFilter("CandViewCountFilter",
0055                                      src = cms.InputTag("eggoodSuperClusters"),
0056                                      minNumber = cms.uint32(2)
0057                                      )
0058
0059 oneEmCluster = cms.Sequence(
0060         egSuperClusterMerger+egSuperClusterCands+eggoodSuperClusters+egSuperClusterPt5Filter
0061             )
0062
0063 gammaJet = cms.EDFilter("EtaPtMinCandViewSelector",
0064                           src = cms.InputTag("iterativeCone5CaloJets"),
0065                           ptMin   = cms.double(5),
0066                           etaMin = cms.double(-2),
0067                           etaMax = cms.double(2)
0068                           )
0069
0070 gammaJetFilter = cms.EDFilter("CandViewCountFilter",
0071                               src = cms.InputTag("gammaJet"),
0072                               minNumber = cms.uint32(2)
0073                               )
0074
0075 #Define group sequence, using HLT/Reco quality cut.

0076 #exoticaMuHLTQualitySeq = cms.Sequence()

0077 singlePhotonPt20QualitySeq = cms.Sequence(
0078     #exoticaMuHLT+

0079     singlePhotonPt20Filter
0080 )
0081 singlePhotonPt15QualitySeq = cms.Sequence(
0082         #exoticaMuHLT+

0083         singlePhotonPt15Filter
0084         )
0085 singlePhotonPt10QualitySeq = cms.Sequence(
0086         #exoticaMuHLT+

0087         singlePhotonPt10Filter
0088         )
0089 singlePhotonPt5QualitySeq = cms.Sequence(
0090         #exoticaMuHLT+

0091         oneEmCluster+gammaJet+gammaJetFilter+singlePhotonPt5Filter
0092         )
0093