EgammaPhotonProducers/python/gedPhotons_cfi.py

0001 import FWCore.ParameterSet.Config as cms
0002
0003 from RecoEgamma.PhotonIdentification.pfIsolationCalculator_cfi import *
0004 from RecoEgamma.PhotonIdentification.isolationCalculator_cfi import *
0005 from RecoEgamma.PhotonIdentification.mipVariable_cfi import *
0006 from RecoEgamma.PhotonIdentification.mvaHaloVariable_cfi import *
0007 from RecoEcal.EgammaClusterProducers.hybridSuperClusters_cfi import *
0008 from RecoEcal.EgammaClusterProducers.multi5x5BasicClusters_cfi import *
0009
0010 from RecoEgamma.EgammaTools.regressionModifier_cfi import *
0011
0012 from RecoEgamma.EgammaIsolationAlgos.egammaHBHERecHitThreshold_cff import egammaHBHERecHit
0013
0014 from RecoEgamma.EgammaIsolationAlgos.egammaEcalPFClusterIsolationProducerRecoPhoton_cfi import egammaEcalPFClusterIsolationProducerRecoPhoton
0015 from RecoEgamma.EgammaIsolationAlgos.egammaHcalPFClusterIsolationProducerRecoPhoton_cfi import egammaHcalPFClusterIsolationProducerRecoPhoton
0016
0017 #
0018 # producer for photons
0019 #
0020 gedPhotons = cms.EDProducer("GEDPhotonProducer",
0021     photonProducer = cms.InputTag("gedPhotonsTmp"),
0022     reconstructionStep = cms.string("tmp"),
0023     #old regression <<<<<< >>>>>> do not use
0024     regressionWeightsFromDB =   cms.bool(True),
0025     energyRegressionWeightsFileLocation = cms.string('/afs/cern.ch/user/b/bendavid/cmspublic/regweights/gbrph.root'),
0026     energyRegressionWeightsDBLocation = cms.string('wgbrph'),
0027     # refined SC regression setup
0028     useRegression = cms.bool(True),
0029     regressionConfig = regressionModifier.clone(),
0030     superClusterEnergyCorrFunction =  cms.string("EcalClusterEnergyCorrection"),
0031     superClusterEnergyErrorFunction = cms.string("EcalClusterEnergyUncertainty"),
0032     superClusterCrackEnergyCorrFunction =  cms.string("EcalClusterCrackCorrection"),
0033     photonEcalEnergyCorrFunction = cms.string("EcalClusterEnergyCorrectionObjectSpecific"),
0034     pfEgammaCandidates = cms.InputTag("particleFlowEGamma"),
0035     pfCandidates = cms.InputTag("particleFlowTmp"),
0036     outputPhotonCollection = cms.string(""),
0037     valueMapPhotons = cms.string("valMapPFEgammaCandToPhoton"),
0038     #candidateP4type = cms.string("fromRegression"),
0039     candidateP4type = cms.string("fromRefinedSCRegression"),
0040     isolationSumsCalculatorSet = cms.PSet(isolationSumsCalculator),
0041     PFIsolationCalculatorSet = cms.PSet(pfIsolationCalculator),
0042     mipVariableSet = cms.PSet(mipVariable),
0043     mvaBasedHaloVariableSet = cms.PSet(mvaHaloVariable),
0044     usePrimaryVertex = cms.bool(True),
0045     primaryVertexProducer = cms.InputTag('offlinePrimaryVerticesWithBS'),
0046     posCalc_t0_endcPresh = cms.double(3.6),
0047     posCalc_logweight = cms.bool(True),
0048     posCalc_w0 = cms.double(4.2),
0049     hbheInstance = cms.string(''),
0050     posCalc_t0_endc = cms.double(6.3),
0051     barrelEcalHits = cms.InputTag("ecalRecHit","EcalRecHitsEB"),
0052     hbheModule = cms.string('hbhereco'),
0053     endcapEcalHits = cms.InputTag("ecalRecHit","EcalRecHitsEE"),
0054     preshowerHits = cms.InputTag("ecalPreshowerRecHit","EcalRecHitsES"),
0055     runMIPTagger = cms.bool(True),
0056     runMVABasedHaloTagger = cms.bool(True),
0057     highEt  = cms.double(100.),
0058     minR9Barrel = cms.double(0.94),
0059     minR9Endcap = cms.double(0.95),
0060     multThresEB = cms.double(1.00),
0061     multThresEE = cms.double(1.25),
0062     hOverEConeSize = cms.double(0.15),
0063     hbheRecHits = egammaHBHERecHit.hbheRecHits,
0064     recHitEThresholdHB = egammaHBHERecHit.recHitEThresholdHB,
0065     recHitEThresholdHE = egammaHBHERecHit.recHitEThresholdHE,
0066     usePFThresholdsFromDB = egammaHBHERecHit.usePFThresholdsFromDB,
0067     maxHcalRecHitSeverity = egammaHBHERecHit.maxHcalRecHitSeverity,
0068     hcalRun2EffDepth = cms.bool(False),
0069     posCalc_x0 = cms.double(0.89),
0070     posCalc_t0_barl = cms.double(7.7),
0071     minSCEtBarrel = cms.double(10.0),
0072     minSCEtEndcap = cms.double(10.0),
0073     maxHoverEBarrel = cms.double(0.5),
0074     maxHoverEEndcap = cms.double(0.5),
0075     ecalRecHitSumEtOffsetBarrel = cms.double(999999999),
0076     ecalRecHitSumEtSlopeBarrel = cms.double(0.),
0077     ecalRecHitSumEtOffsetEndcap = cms.double(999999999),
0078     ecalRecHitSumEtSlopeEndcap = cms.double(0.),
0079     hcalRecHitSumEtOffsetBarrel = cms.double(999999999),
0080     hcalRecHitSumEtSlopeBarrel = cms.double(0.),
0081     hcalRecHitSumEtOffsetEndcap = cms.double(999999999),
0082     hcalRecHitSumEtSlopeEndcap = cms.double(0.),
0083     nTrackSolidConeBarrel =cms.double(999999999),
0084     nTrackSolidConeEndcap =cms.double(999999999),
0085     nTrackHollowConeBarrel =cms.double(999999999),
0086     nTrackHollowConeEndcap =cms.double(999999999),
0087     trackPtSumSolidConeBarrel =cms.double(999999999),
0088     trackPtSumSolidConeEndcap =cms.double(999999999),
0089     trackPtSumHollowConeBarrel =cms.double(999999999),
0090     trackPtSumHollowConeEndcap =cms.double(999999999),
0091     sigmaIetaIetaCutBarrel=cms.double(999999999),
0092     sigmaIetaIetaCutEndcap=cms.double(999999999),
0093     posCalcParameters = cms.PSet( T0_barl      = cms.double(7.4),
0094                                   T0_endc      = cms.double(6.3),
0095                                   T0_endcPresh = cms.double(3.6),
0096                                   LogWeighted  = cms.bool(True),
0097                                   W0           = cms.double(4.2),
0098                                   X0           = cms.double(0.89)
0099                                  ),
0100
0101     RecHitFlagToBeExcludedEB = cleanedHybridSuperClusters.RecHitFlagToBeExcluded,
0102     RecHitSeverityToBeExcludedEB = cleanedHybridSuperClusters.RecHitSeverityToBeExcluded,
0103     RecHitFlagToBeExcludedEE = multi5x5BasicClustersCleaned.RecHitFlagToBeExcluded,
0104     RecHitSeverityToBeExcludedEE = cleanedHybridSuperClusters.RecHitSeverityToBeExcluded,
0105     checkHcalStatus = cms.bool(True),
0106     PhotonDNNPFid = cms.PSet(
0107         enabled = cms.bool(False),
0108         inputTensorName = cms.string("FirstLayer_input"),
0109         outputTensorName = cms.string("sequential/FinalLayer/Sigmoid"),
0110         modelsFiles = cms.vstring(
0111                                 'RecoEgamma/PhotonIdentification/data/Photon_PFID_dnn/v1/EB/EB_modelDNN.pb',
0112                                 'RecoEgamma/PhotonIdentification/data/Photon_PFID_dnn/v1/EE/EE_modelDNN.pb'),
0113         scalersFiles = cms.vstring(
0114                     'RecoEgamma/PhotonIdentification/data/Photon_PFID_dnn/v1/EB/EB_scaler.txt',
0115                     'RecoEgamma/PhotonIdentification/data/Photon_PFID_dnn/v1/EE/EE_scaler.txt'
0116         ),
0117         outputDim = cms.vuint32(1,1),
0118         useEBModelInGap = cms.bool(True)
0119     ),
0120
0121     pfECALClusIsolCfg = cms.PSet(
0122         pfClusterProducer = egammaEcalPFClusterIsolationProducerRecoPhoton.pfClusterProducer,
0123         drMax = egammaEcalPFClusterIsolationProducerRecoPhoton.drMax,
0124         drVetoBarrel = egammaEcalPFClusterIsolationProducerRecoPhoton.drVetoBarrel,
0125         drVetoEndcap = egammaEcalPFClusterIsolationProducerRecoPhoton.drVetoEndcap,
0126         etaStripBarrel = egammaEcalPFClusterIsolationProducerRecoPhoton.etaStripBarrel,
0127         etaStripEndcap = egammaEcalPFClusterIsolationProducerRecoPhoton.etaStripEndcap,
0128         energyBarrel = egammaEcalPFClusterIsolationProducerRecoPhoton.energyBarrel,
0129         energyEndcap = egammaEcalPFClusterIsolationProducerRecoPhoton.energyEndcap
0130     ),
0131
0132     pfHCALClusIsolCfg = cms.PSet(
0133
0134         pfClusterProducerHCAL = egammaHcalPFClusterIsolationProducerRecoPhoton.pfClusterProducerHCAL,
0135         useHF = egammaHcalPFClusterIsolationProducerRecoPhoton.useHF,
0136         pfClusterProducerHFEM = egammaHcalPFClusterIsolationProducerRecoPhoton.pfClusterProducerHFEM,
0137         pfClusterProducerHFHAD = egammaHcalPFClusterIsolationProducerRecoPhoton.pfClusterProducerHFHAD,
0138         drMax = egammaHcalPFClusterIsolationProducerRecoPhoton.drMax,
0139         drVetoBarrel = egammaHcalPFClusterIsolationProducerRecoPhoton.drVetoBarrel,
0140         drVetoEndcap = egammaHcalPFClusterIsolationProducerRecoPhoton.drVetoEndcap,
0141         etaStripBarrel = egammaHcalPFClusterIsolationProducerRecoPhoton.etaStripBarrel,
0142         etaStripEndcap = egammaHcalPFClusterIsolationProducerRecoPhoton.etaStripEndcap,
0143         energyBarrel = egammaHcalPFClusterIsolationProducerRecoPhoton.energyBarrel,
0144         energyEndcap = egammaHcalPFClusterIsolationProducerRecoPhoton.energyEndcap,
0145         useEt = egammaHcalPFClusterIsolationProducerRecoPhoton.useEt,
0146
0147     )
0148
0149 )
0150
0151