Skimming/python/PDWG_DiPhoton_SD_cff.py

0001 import FWCore.ParameterSet.Config as cms
0002
0003 #                                __                 __
0004 #    ____  ____   ____   _______/  |______    _____/  |_  ______
0005 #  _/ ___\/  _ \ /    \ /  ___/\   __\__  \  /    \   __\/  ___/
0006 #  \  \__(  <_> )   |  \\___ \  |  |  / __ \|   |  \  |  \___ \
0007 #   \___  >____/|___|  /____  > |__| (____  /___|  /__| /____  >
0008 #        \/           \/     \/            \/     \/          \/
0009
0010 PHOTON_CALOIDISO_ET_HIGH_CUT_MIN = 26.
0011 PHOTON_CALOIDISO_ET_LOW_CUT_MIN = 18.
0012 MASS_DIPHOTON_CALOIDISO_CUT_MIN = 70.
0013
0014 PHOTON_R9ID_ET_HIGH_CUT_MIN = 26.
0015 PHOTON_R9ID_ET_LOW_CUT_MIN = 18.
0016 MASS_DIPHOTON_R9ID_CUT_MIN = 70.
0017
0018 MASS_DIPHOTON_MIXEDID_CUT_MIN = 70.
0019
0020 #  _____ __  _____            _         _   _
0021 # |  |  |  ||_   _|   ___ ___| |___ ___| |_|_|___ ___
0022 # |     |  |__| |    |_ -| -_| | -_|  _|  _| | . |   |
0023 # |__|__|_____|_|    |___|___|_|___|___|_| |_|___|_|_|
0024
0025 import copy
0026 from HLTrigger.HLTfilters.hltHighLevel_cfi import *
0027 DiPhotonHltFilter = copy.deepcopy(hltHighLevel)
0028 DiPhotonHltFilter.throw = cms.bool(False)
0029 DiPhotonHltFilter.HLTPaths = ["HLT_Photon*_Photon*"]
0030
0031 #  _____     _     _____   _ _____
0032 # |     |___| |___|     |_| |     |___ ___
0033 # |   --| .'| | . |-   -| . |-   -|_ -| . |
0034 # |_____|__,|_|___|_____|___|_____|___|___|
0035
0036 hltDiPhotonCaloIdIsoObjectProducer = cms.EDProducer("CandidateTriggerObjectProducer",
0037                                              triggerName = cms.string("HLT_Photon.*_CaloId.*_Iso.*_Photon.*_CaloId.*_Iso.*_.*"),
0038                                              triggerResults = cms.InputTag("TriggerResults","","HLT"),
0039                                              triggerEvent   = cms.InputTag("hltTriggerSummaryAOD","","HLT"),
0040                                              stageL1Trigger = cms.uint32(1)
0041                                              )
0042 from Configuration.Eras.Modifier_stage2L1Trigger_cff import stage2L1Trigger
0043 stage2L1Trigger.toModify(hltDiPhotonCaloIdIsoObjectProducer, stageL1Trigger = 2)
0044
0045 TrailingPtCaloIdIsoPhotons = cms.EDFilter("CandViewRefSelector",
0046     filter = cms.bool(True),
0047     src = cms.InputTag("hltDiPhotonCaloIdIsoObjectProducer"),
0048     cut = cms.string('pt > '+str(PHOTON_CALOIDISO_ET_LOW_CUT_MIN))
0049 )
0050
0051 LeadingPtCaloIdIsoPhotons = cms.EDFilter("CandViewRefSelector",
0052     filter = cms.bool(True),
0053     src = cms.InputTag("hltDiPhotonCaloIdIsoObjectProducer"),
0054     cut = cms.string('pt > '+str(PHOTON_CALOIDISO_ET_HIGH_CUT_MIN))
0055 )
0056
0057 CaloIdIsoPhotonPairs = cms.EDProducer("CandViewShallowCloneCombiner",
0058     decay = cms.string("LeadingPtCaloIdIsoPhotons TrailingPtCaloIdIsoPhotons"), # charge coniugate states are implied
0059     checkCharge = cms.bool(False),
0060     cut   = cms.string("mass > " + str(MASS_DIPHOTON_CALOIDISO_CUT_MIN))
0061 )
0062
0063 CaloIdIsoPhotonPairsCounter = cms.EDFilter("CandViewCountFilter",
0064                                     src = cms.InputTag("CaloIdIsoPhotonPairs"),
0065                                     minNumber = cms.uint32(1)
0066                                     )
0067
0068 CaloIdIsoPhotonPairsFilter = cms.Sequence(DiPhotonHltFilter*hltDiPhotonCaloIdIsoObjectProducer*TrailingPtCaloIdIsoPhotons*LeadingPtCaloIdIsoPhotons*CaloIdIsoPhotonPairs*CaloIdIsoPhotonPairsCounter)
0069
0070 #  _____ ___ _____   _
0071 # | __  | . |     |_| |
0072 # |    -|_  |-   -| . |
0073 # |__|__|___|_____|___|
0074
0075 hltDiPhotonR9IdObjectProducer = hltDiPhotonCaloIdIsoObjectProducer.clone(
0076                                              triggerName = cms.string("HLT_Photon.*_R9Id.*_Photon.*_R9Id.*_.*"),
0077                                              )
0078
0079 TrailingPtR9IdPhotons = TrailingPtCaloIdIsoPhotons.clone(
0080     src = cms.InputTag("hltDiPhotonR9IdObjectProducer"),
0081     cut = cms.string('pt > '+str(PHOTON_R9ID_ET_LOW_CUT_MIN))
0082 )
0083
0084 LeadingPtR9IdPhotons = LeadingPtCaloIdIsoPhotons.clone(
0085     src = cms.InputTag("hltDiPhotonR9IdObjectProducer"),
0086     cut = cms.string('pt > '+str(PHOTON_R9ID_ET_HIGH_CUT_MIN))
0087 )
0088
0089 R9IdPhotonPairs = CaloIdIsoPhotonPairs.clone(
0090     decay = cms.string("LeadingPtR9IdPhotons TrailingPtR9IdPhotons"), # charge coniugate states are implied
0091     cut   = cms.string("mass > " + str(MASS_DIPHOTON_R9ID_CUT_MIN))
0092 )
0093
0094 R9IdPhotonPairsCounter = CaloIdIsoPhotonPairsCounter.clone(
0095                                     src = cms.InputTag("R9IdPhotonPairs"),
0096                                     )
0097
0098 R9IdPhotonPairsFilter = cms.Sequence(DiPhotonHltFilter*hltDiPhotonR9IdObjectProducer*TrailingPtR9IdPhotons*LeadingPtR9IdPhotons*R9IdPhotonPairs*R9IdPhotonPairsCounter)
0099
0100 #  __   __ _ _   _ ____   _    _____     _     _____   _ _____           _____ ___ _____   _
0101 # |  \ /  | | \_/ |  __|_| |  |     |___| |___|     |_| |     |___ ___  | __  | . |     |_| |
0102 # | | v | | |\ _ /|  __| . |  |   --| .'| | . |-   -| . |-   -|_ -| . | |    -|_  |-   -| . |
0103 # |_|   |_|_|_/ \_|____|___|  |_____|__,|_|___|_____|___|_____|___|___| |__|__|___|_____|___|
0104
0105
0106 hltDiPhotonMixedCaloR9IdObjectProducer = hltDiPhotonCaloIdIsoObjectProducer.clone(
0107                                                     triggerName = cms.string("HLT_Photon.*_CaloId.*_Iso.*_Photon.*_R9Id.*_.*"),
0108                                              )
0109
0110 TrailingPtMixedCaloR9IdPhotons = TrailingPtCaloIdIsoPhotons.clone(
0111    src = cms.InputTag("hltDiPhotonMixedCaloR9IdObjectProducer"),
0112    cut = cms.string('pt > '+str(PHOTON_R9ID_ET_LOW_CUT_MIN))
0113 )
0114
0115 LeadingPtMixedCaloR9IdPhotons = LeadingPtCaloIdIsoPhotons.clone(
0116     src = cms.InputTag("hltDiPhotonMixedCaloR9IdObjectProducer"),
0117     cut = cms.string('pt > '+str(PHOTON_CALOIDISO_ET_HIGH_CUT_MIN))
0118 )
0119
0120 MixedCaloR9IdPhotonPairs = CaloIdIsoPhotonPairs.clone(
0121     decay = cms.string("LeadingPtMixedCaloR9IdPhotons TrailingPtMixedCaloR9IdPhotons"), # charge coniugate states are implied
0122     cut   = cms.string("mass > " + str(MASS_DIPHOTON_MIXEDID_CUT_MIN))
0123 )
0124
0125 MixedCaloR9IdPhotonPairsCounter = CaloIdIsoPhotonPairsCounter.clone(
0126                                     src = cms.InputTag("MixedCaloR9IdPhotonPairs"),
0127                                     )
0128
0129 MixedCaloR9IdPhotonPairsFilter = cms.Sequence(DiPhotonHltFilter*hltDiPhotonMixedCaloR9IdObjectProducer*TrailingPtMixedCaloR9IdPhotons*LeadingPtMixedCaloR9IdPhotons*MixedCaloR9IdPhotonPairs*MixedCaloR9IdPhotonPairsCounter)
0130
0131 #  __   __ _ _   _ ____   _    _____ ___ _____   _   _____     _     _____   _ _____
0132 # |  \ /  | | \_/ |  __|_| |  | __  | . |     |_| | |     |___| |___|     |_| |     |___ ___
0133 # | | v | | |\ _ /|  __| . |  |    -|_  |-   -| . | |   --| .'| | . |-   -| . |-   -|_ -| . |
0134 # |_|   |_|_|_/ \_|____|___|  |__|__|___|_____|___| |_____|__,|_|___|_____|___|_____|___|___|
0135
0136 hltDiPhotonMixedR9CaloIdObjectProducer = hltDiPhotonCaloIdIsoObjectProducer.clone(
0137                                              triggerName = cms.string("HLT_Photon.*_R9Id.*_Photon.*_CaloId.*_Iso.*_.*"),
0138                                              )
0139
0140 TrailingPtMixedR9CaloIdPhotons = TrailingPtCaloIdIsoPhotons.clone(
0141     src = cms.InputTag("hltDiPhotonMixedR9CaloIdObjectProducer"),
0142     cut = cms.string('pt > '+str(PHOTON_CALOIDISO_ET_LOW_CUT_MIN))
0143 )
0144
0145 LeadingPtMixedR9CaloIdPhotons = LeadingPtCaloIdIsoPhotons.clone(
0146     src = cms.InputTag("hltDiPhotonMixedR9CaloIdObjectProducer"),
0147     cut = cms.string('pt > '+str(PHOTON_R9ID_ET_HIGH_CUT_MIN))
0148 )
0149
0150 MixedR9CaloIdPhotonPairs = CaloIdIsoPhotonPairs.clone(
0151     decay = cms.string("LeadingPtMixedR9CaloIdPhotons TrailingPtMixedR9CaloIdPhotons"), # charge coniugate states are implied
0152     cut   = cms.string("mass > " + str(MASS_DIPHOTON_MIXEDID_CUT_MIN))
0153 )
0154
0155 MixedR9CaloIdPhotonPairsCounter = CaloIdIsoPhotonPairsCounter.clone(
0156                                     src = cms.InputTag("MixedR9CaloIdPhotonPairs"),
0157                                     )
0158
0159 MixedR9CaloIdPhotonPairsFilter = cms.Sequence(DiPhotonHltFilter*hltDiPhotonMixedR9CaloIdObjectProducer*TrailingPtMixedR9CaloIdPhotons*LeadingPtMixedR9CaloIdPhotons*MixedR9CaloIdPhotonPairs*MixedR9CaloIdPhotonPairsCounter)