CMS 3D CMS Logo

/afs/cern.ch/work/a/aaltunda/public/www/CMSSW_5_3_14/src/Configuration/Skimming/python/PDWG_TauSkim_cff.py

Go to the documentation of this file.
00001 import FWCore.ParameterSet.Config as cms
00002 
00003 from HLTrigger.HLTfilters.hltHighLevel_cfi import hltHighLevel as hlt_selector
00004 hlt_selector.throw = cms.bool(False)
00005 TauSkimMuTauMETHLT = hlt_selector.clone()
00006 TauSkimMuTauMETHLT.TriggerResultsTag = cms.InputTag('TriggerResults', '', 'HLT')
00007 TauSkimMuTauMETHLT.andOr=cms.bool(True)
00008 TauSkimMuTauMETHLT.HLTPaths = cms.vstring("HLT_IsoMu15_eta2p1_L1ETM20_v*", "HLT_IsoMu15Rho_eta2p1_L1ETM20_v*")
00009 
00010 TauSkimPFTausSelected = cms.EDFilter("PFTauSelector",
00011    src = cms.InputTag("hpsPFTauProducer"),
00012    discriminators = cms.VPSet(
00013     cms.PSet( discriminator=cms.InputTag("hpsPFTauDiscriminationByDecayModeFinding"),
00014                    selectionCut=cms.double(0.5)           
00015         ),
00016     cms.PSet( discriminator=cms.InputTag("hpsPFTauDiscriminationByLooseCombinedIsolationDBSumPtCorr"),
00017                    selectionCut=cms.double(0.5)           
00018         ),
00019    ),
00020    cut = cms.string('pt > 22. && abs(eta) < 2.3') 
00021 )
00022 
00023 TauSkimPFTauSkimmedBy1 = cms.EDFilter("CandViewCountFilter",
00024  src = cms.InputTag('TauSkimPFTausSelected'),
00025  minNumber = cms.uint32(1)
00026 )
00027 
00028 TauSkimPFTauSkimmedBy2 = cms.EDFilter("CandViewCountFilter",
00029  src = cms.InputTag('TauSkimPFTausSelected'),
00030  minNumber = cms.uint32(2)
00031 )
00032 
00033 TauSkimDiTauPairs  = cms.EDProducer("CandViewShallowCloneCombiner",
00034                                     decay = cms.string("TauSkimPFTausSelected TauSkimPFTausSelected"),
00035                                     checkCharge = cms.bool(False),
00036                                     cut         = cms.string("sqrt((daughter(0).eta-daughter(1).eta)*(daughter(0).eta-daughter(1).eta)+  min( abs(daughter(0).phi-daughter(1).phi), 2*3.1415926 - abs(daughter(0).phi-daughter(1).phi)  ) *  min( abs(daughter(0).phi-daughter(1).phi), 2*3.1415926 - abs(daughter(0).phi-daughter(1).phi)  )  )>0.5"),
00037                                     )
00038 
00039 TauSkimDiTauPairFilter = cms.EDFilter("CandViewCountFilter",
00040                                       src = cms.InputTag("TauSkimDiTauPairs"),
00041                                       minNumber = cms.uint32(1)
00042                                       )
00043 
00044 TauSkimPFTausSelectedForMuTau = TauSkimPFTausSelected.clone()
00045 TauSkimPFTausSelectedForMuTau.discriminators = cms.VPSet(
00046     cms.PSet( discriminator=cms.InputTag("hpsPFTauDiscriminationByDecayModeFinding"),
00047               selectionCut=cms.double(0.5)
00048               ),
00049     cms.PSet( discriminator=cms.InputTag("hpsPFTauDiscriminationByLooseCombinedIsolationDBSumPtCorr"),
00050               selectionCut=cms.double(0.5)
00051               ),
00052     cms.PSet( discriminator=cms.InputTag("hpsPFTauDiscriminationByLooseMuonRejection"),
00053               selectionCut=cms.double(0.5)
00054               ),
00055     )
00056 TauSkimPFTausSelectedForMuTau.cut = cms.string('pt > 18. && abs(eta) < 2.3')
00057 TauSkimPFTauSkimmedForMuTau = cms.EDFilter("CandViewCountFilter",
00058                                            src = cms.InputTag('TauSkimPFTausSelectedForMuTau'),
00059                                            minNumber = cms.uint32(1)
00060                                            )
00061 
00062 TauSkimPFTausSelectedForMuTauMET = TauSkimPFTausSelected.clone()
00063 TauSkimPFTausSelectedForMuTauMET.cut = cms.string('pt > 18. && abs(eta) < 2.3')
00064 TauSkimPFTauSkimmedForMuTauMET = cms.EDFilter("CandViewCountFilter",
00065                                               src = cms.InputTag('TauSkimPFTausSelectedForMuTauMET'),
00066                                               minNumber = cms.uint32(1)
00067                                               )
00068 
00069 TauSkimMuonSelected = cms.EDFilter("MuonRefSelector",
00070                                    src = cms.InputTag("muons"),
00071                                    cut = cms.string("pt > 15 && " + "abs(eta) < 2.4 && " +
00072                                                     "isGlobalMuon && isTrackerMuon" +
00073                                                     " && globalTrack.isNonnull "+
00074                                                     " && globalTrack.hitPattern.numberOfValidTrackerHits>=5"+
00075                                                     " && globalTrack.normalizedChi2<20"+
00076                                                     " && (pfIsolationR03.sumChargedHadronPt/pt) < 0.3"
00077                                                     ),
00078                                    )
00079 
00080 TauSkimMuonSkimmedBy1 = cms.EDFilter("CandViewCountFilter",
00081                                      src = cms.InputTag('TauSkimMuonSelected'),
00082                                      minNumber = cms.uint32(1)
00083                                      )
00084 
00085 TauSkimMuTauPairs  = cms.EDProducer("CandViewShallowCloneCombiner",
00086                                     decay = cms.string("TauSkimMuonSelected TauSkimPFTausSelectedForMuTau"),
00087                                     checkCharge = cms.bool(False),
00088                                     cut         = cms.string("sqrt((daughter(0).eta-daughter(1).eta)*(daughter(0).eta-daughter(1).eta)+  min( abs(daughter(0).phi-daughter(1).phi), 2*3.1415926 - abs(daughter(0).phi-daughter(1).phi)  ) *  min( abs(daughter(0).phi-daughter(1).phi), 2*3.1415926 - abs(daughter(0).phi-daughter(1).phi)  )  )>0.5"),
00089                                     )
00090 
00091 TauSkimMuTauMETPairs  = cms.EDProducer("CandViewShallowCloneCombiner",
00092                                        decay = cms.string("TauSkimMuonSelected TauSkimPFTausSelectedForMuTauMET"),
00093                                        checkCharge = cms.bool(False),
00094                                        cut         = cms.string("sqrt((daughter(0).eta-daughter(1).eta)*(daughter(0).eta-daughter(1).eta)+  min( abs(daughter(0).phi-daughter(1).phi), 2*3.1415926 - abs(daughter(0).phi-daughter(1).phi)  ) *  min( abs(daughter(0).phi-daughter(1).phi), 2*3.1415926 - abs(daughter(0).phi-daughter(1).phi)  )  )>0.5"),
00095                                        )
00096 
00097 TauSkimMuTauPairFilter = cms.EDFilter("CandViewCountFilter",
00098                                       src = cms.InputTag("TauSkimMuTauPairs"),
00099                                       minNumber = cms.uint32(1)
00100                                       )
00101 TauSkimMuTauMETPairFilter = cms.EDFilter("CandViewCountFilter",
00102                                          src = cms.InputTag("TauSkimMuTauMETPairs"),
00103                                          minNumber = cms.uint32(1)
00104                                          )
00105 
00106 
00107 tauSkim1Sequence = cms.Sequence(
00108     TauSkimPFTausSelected *
00109     TauSkimPFTauSkimmedBy1
00110     )
00111 
00112 tauSkim2Sequence = cms.Sequence(
00113     TauSkimPFTausSelected *
00114     TauSkimPFTauSkimmedBy2 *
00115     TauSkimDiTauPairs *
00116     TauSkimDiTauPairFilter
00117     )
00118 
00119 mutauSkimSequence = cms.Sequence(
00120     TauSkimMuonSelected *
00121     TauSkimMuonSkimmedBy1 *
00122     TauSkimPFTausSelectedForMuTau *
00123     TauSkimPFTauSkimmedForMuTau *
00124     TauSkimMuTauPairs *
00125     TauSkimMuTauPairFilter
00126     )
00127 
00128 mutauMETSkimSequence = cms.Sequence(
00129     TauSkimMuTauMETHLT *
00130     TauSkimMuonSelected *
00131     TauSkimMuonSkimmedBy1 *
00132     TauSkimPFTausSelectedForMuTauMET *
00133     TauSkimPFTauSkimmedForMuTauMET *
00134     TauSkimMuTauMETPairs *
00135     TauSkimMuTauMETPairFilter
00136     )