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     maxHcalRecHitSeverity = egammaHBHERecHit.maxHcalRecHitSeverity,
0067     hcalRun2EffDepth = cms.bool(False),
0068     posCalc_x0 = cms.double(0.89),
0069     posCalc_t0_barl = cms.double(7.7),
0070     minSCEtBarrel = cms.double(10.0),
0071     minSCEtEndcap = cms.double(10.0),
0072     maxHoverEBarrel = cms.double(0.5),
0073     maxHoverEEndcap = cms.double(0.5),
0074     ecalRecHitSumEtOffsetBarrel = cms.double(999999999),
0075     ecalRecHitSumEtSlopeBarrel = cms.double(0.),
0076     ecalRecHitSumEtOffsetEndcap = cms.double(999999999),
0077     ecalRecHitSumEtSlopeEndcap = cms.double(0.),
0078     hcalRecHitSumEtOffsetBarrel = cms.double(999999999),
0079     hcalRecHitSumEtSlopeBarrel = cms.double(0.),
0080     hcalRecHitSumEtOffsetEndcap = cms.double(999999999),
0081     hcalRecHitSumEtSlopeEndcap = cms.double(0.),
0082     nTrackSolidConeBarrel =cms.double(999999999),
0083     nTrackSolidConeEndcap =cms.double(999999999),
0084     nTrackHollowConeBarrel =cms.double(999999999),
0085     nTrackHollowConeEndcap =cms.double(999999999),
0086     trackPtSumSolidConeBarrel =cms.double(999999999),
0087     trackPtSumSolidConeEndcap =cms.double(999999999),
0088     trackPtSumHollowConeBarrel =cms.double(999999999),
0089     trackPtSumHollowConeEndcap =cms.double(999999999),
0090     sigmaIetaIetaCutBarrel=cms.double(999999999),
0091     sigmaIetaIetaCutEndcap=cms.double(999999999),
0092     posCalcParameters = cms.PSet( T0_barl      = cms.double(7.4),
0093                                   T0_endc      = cms.double(6.3),
0094                                   T0_endcPresh = cms.double(3.6),
0095                                   LogWeighted  = cms.bool(True),
0096                                   W0           = cms.double(4.2),
0097                                   X0           = cms.double(0.89)
0098                                  ),
0099
0100     RecHitFlagToBeExcludedEB = cleanedHybridSuperClusters.RecHitFlagToBeExcluded,
0101     RecHitSeverityToBeExcludedEB = cleanedHybridSuperClusters.RecHitSeverityToBeExcluded,
0102     RecHitFlagToBeExcludedEE = multi5x5BasicClustersCleaned.RecHitFlagToBeExcluded,
0103     RecHitSeverityToBeExcludedEE = cleanedHybridSuperClusters.RecHitSeverityToBeExcluded,
0104     checkHcalStatus = cms.bool(True),
0105     PhotonDNNPFid = cms.PSet(
0106         enabled = cms.bool(False),
0107         inputTensorName = cms.string("FirstLayer_input"),
0108         outputTensorName = cms.string("sequential/FinalLayer/Sigmoid"),
0109         modelsFiles = cms.vstring(
0110                                 'RecoEgamma/PhotonIdentification/data/Photon_PFID_dnn/v1/EB/EB_modelDNN.pb',
0111                                 'RecoEgamma/PhotonIdentification/data/Photon_PFID_dnn/v1/EE/EE_modelDNN.pb'),
0112         scalersFiles = cms.vstring(
0113                     'RecoEgamma/PhotonIdentification/data/Photon_PFID_dnn/v1/EB/EB_scaler.txt',
0114                     'RecoEgamma/PhotonIdentification/data/Photon_PFID_dnn/v1/EE/EE_scaler.txt'
0115         ),
0116         outputDim = cms.vuint32(1,1),
0117         useEBModelInGap = cms.bool(True)
0118     ),
0119
0120     pfECALClusIsolCfg = cms.PSet(
0121         pfClusterProducer = egammaEcalPFClusterIsolationProducerRecoPhoton.pfClusterProducer,
0122         drMax = egammaEcalPFClusterIsolationProducerRecoPhoton.drMax,
0123         drVetoBarrel = egammaEcalPFClusterIsolationProducerRecoPhoton.drVetoBarrel,
0124         drVetoEndcap = egammaEcalPFClusterIsolationProducerRecoPhoton.drVetoEndcap,
0125         etaStripBarrel = egammaEcalPFClusterIsolationProducerRecoPhoton.etaStripBarrel,
0126         etaStripEndcap = egammaEcalPFClusterIsolationProducerRecoPhoton.etaStripEndcap,
0127         energyBarrel = egammaEcalPFClusterIsolationProducerRecoPhoton.energyBarrel,
0128         energyEndcap = egammaEcalPFClusterIsolationProducerRecoPhoton.energyEndcap
0129     ),
0130
0131     pfHCALClusIsolCfg = cms.PSet(
0132
0133         pfClusterProducerHCAL = egammaHcalPFClusterIsolationProducerRecoPhoton.pfClusterProducerHCAL,
0134         useHF = egammaHcalPFClusterIsolationProducerRecoPhoton.useHF,
0135         pfClusterProducerHFEM = egammaHcalPFClusterIsolationProducerRecoPhoton.pfClusterProducerHFEM,
0136         pfClusterProducerHFHAD = egammaHcalPFClusterIsolationProducerRecoPhoton.pfClusterProducerHFHAD,
0137         drMax = egammaHcalPFClusterIsolationProducerRecoPhoton.drMax,
0138         drVetoBarrel = egammaHcalPFClusterIsolationProducerRecoPhoton.drVetoBarrel,
0139         drVetoEndcap = egammaHcalPFClusterIsolationProducerRecoPhoton.drVetoEndcap,
0140         etaStripBarrel = egammaHcalPFClusterIsolationProducerRecoPhoton.etaStripBarrel,
0141         etaStripEndcap = egammaHcalPFClusterIsolationProducerRecoPhoton.etaStripEndcap,
0142         energyBarrel = egammaHcalPFClusterIsolationProducerRecoPhoton.energyBarrel,
0143         energyEndcap = egammaHcalPFClusterIsolationProducerRecoPhoton.energyEndcap,
0144         useEt = egammaHcalPFClusterIsolationProducerRecoPhoton.useEt,
0145
0146     )
0147
0148 )
0149
0150