CMS 3D CMS Logo
List of all members | Classes | Public Member Functions | Private Member Functions | Private Attributes
GEDPhotonProducer Class Reference

#include <GEDPhotonProducer.h>

Inheritance diagram for GEDPhotonProducer:
edm::stream::EDProducer<>

Classes

class  RecoStepInfo
 

Public Member Functions

void beginRun (edm::Run const &r, edm::EventSetup const &es) final
 
void endRun (edm::Run const &, edm::EventSetup const &) final
 
 GEDPhotonProducer (const edm::ParameterSet &ps)
 
void produce (edm::Event &evt, const edm::EventSetup &es) override
 
 ~GEDPhotonProducer () override
 
- Public Member Functions inherited from edm::stream::EDProducer<>
 EDProducer ()=default
 
bool hasAbilityToProduceInLumis () const final
 
bool hasAbilityToProduceInRuns () const final
 

Private Member Functions

void fillPhotonCollection (edm::Event &evt, edm::EventSetup const &es, const edm::Handle< reco::PhotonCoreCollection > &photonCoreHandle, const CaloTopology *topology, const EcalRecHitCollection *ecalBarrelHits, const EcalRecHitCollection *ecalEndcapHits, const EcalRecHitCollection *preshowerHits, const edm::Handle< CaloTowerCollection > &hcalTowersHandle, const reco::VertexCollection &pvVertices, reco::PhotonCollection &outputCollection, int &iSC)
 
void fillPhotonCollection (edm::Event &evt, edm::EventSetup const &es, const edm::Handle< reco::PhotonCollection > &photonHandle, const edm::Handle< reco::PFCandidateCollection > pfCandidateHandle, const edm::Handle< reco::PFCandidateCollection > pfEGCandidateHandle, edm::ValueMap< reco::PhotonRef > pfEGCandToPhotonMap, edm::Handle< reco::VertexCollection > &pvVertices, reco::PhotonCollection &outputCollection, int &iSC, const edm::Handle< edm::ValueMap< float >> &chargedHadrons, const edm::Handle< edm::ValueMap< float >> &neutralHadrons, const edm::Handle< edm::ValueMap< float >> &photons, const edm::Handle< edm::ValueMap< float >> &chargedHadronsWorstVtx, const edm::Handle< edm::ValueMap< float >> &chargedHadronsWorstVtxGeomVeto, const edm::Handle< edm::ValueMap< float >> &chargedHadronsPFPV, const edm::Handle< edm::ValueMap< float >> &pfEcalClusters, const edm::Handle< edm::ValueMap< float >> &pfHcalClusters)
 

Private Attributes

edm::EDGetTokenT< EcalRecHitCollection > barrelEcalHits_
 
std::string candidateP4type_
 
bool checkHcalStatus_
 
edm::ParameterSet conf_
 
std::string conversionCollection_
 
std::string conversionProducer_
 
edm::EDGetTokenT< EcalRecHitCollection > endcapEcalHits_
 
EcalClusterFunctionBaseClass * energyCorrectionF
 
std::vector< int > flagsexclEB_
 
std::vector< int > flagsexclEE_
 
edm::EDGetTokenT< CaloTowerCollection > hcalTowers_
 
double highEt_
 
double hOverEConeSize_
 
double maxHOverE_
 
double minR9Barrel_
 
double minR9Endcap_
 
double minSCEt_
 
edm::EDGetTokenT< edm::ValueMap< std::vector< reco::PFCandidateRef > > > particleBasedIsolationToken
 
edm::EDGetTokenT< reco::PFCandidateCollection > pfCandidates_
 
edm::EDGetTokenT< reco::PFCandidateCollection > pfEgammaCandidates_
 
edm::EDGetTokenT< edm::ValueMap< float > > phoChargedIsolationToken_
 
edm::EDGetTokenT< edm::ValueMap< float > > phoChargedPFPVIsoToken_
 
edm::EDGetTokenT< edm::ValueMap< float > > phoChargedWorstVtxGeomVetoIsoToken_
 
edm::EDGetTokenT< edm::ValueMap< float > > phoChargedWorstVtxIsoToken_
 
edm::EDGetTokenT< edm::ValueMap< float > > phoNeutralHadronIsolationToken_
 
edm::EDGetTokenT< edm::ValueMap< float > > phoPFECALClusIsolationToken_
 
edm::EDGetTokenT< edm::ValueMap< float > > phoPFHCALClusIsolationToken_
 
edm::EDGetTokenT< edm::ValueMap< float > > phoPhotonIsolationToken_
 
std::string photonCollection_
 
edm::EDGetTokenT< reco::PhotonCoreCollection > photonCoreProducerT_
 
edm::InputTag photonProducer_
 
edm::EDGetTokenT< reco::PhotonCollection > photonProducerT_
 
PositionCalc posCalculator_
 
std::vector< double > preselCutValuesBarrel_
 
std::vector< double > preselCutValuesEndcap_
 
edm::EDGetTokenT< EcalRecHitCollection > preshowerHits_
 
RecoStepInfo recoStep_
 
bool runMIPTagger_
 
std::vector< int > severitiesexclEB_
 
std::vector< int > severitiesexclEE_
 
edm::ESHandle< CaloGeometry > theCaloGeom_
 
edm::ESHandle< CaloTopology > theCaloTopo_
 
PhotonEnergyCorrector * thePhotonEnergyCorrector_
 
PhotonIsolationCalculator * thePhotonIsolationCalculator_
 
PhotonMIPHaloTagger * thePhotonMIPHaloTagger_
 
bool usePrimaryVertex_
 
bool validConversions_
 
bool validPixelSeeds_
 
std::string valueMapPFCandPhoton_
 
edm::EDGetTokenT< reco::VertexCollection > vertexProducer_
 

Additional Inherited Members

- Public Types inherited from edm::stream::EDProducer<>
typedef CacheContexts< T... > CacheTypes
 
typedef CacheTypes::GlobalCache GlobalCache
 
typedef AbilityChecker< T... > HasAbility
 
typedef CacheTypes::LuminosityBlockCache LuminosityBlockCache
 
typedef LuminosityBlockContextT< LuminosityBlockCache, RunCache, GlobalCache > LuminosityBlockContext
 
typedef CacheTypes::LuminosityBlockSummaryCache LuminosityBlockSummaryCache
 
typedef CacheTypes::RunCache RunCache
 
typedef RunContextT< RunCache, GlobalCache > RunContext
 
typedef CacheTypes::RunSummaryCache RunSummaryCache
 

Detailed Description

Author
Nancy Marinelli, U. of Notre Dame, US

Definition at line 42 of file GEDPhotonProducer.h.

Constructor & Destructor Documentation

GEDPhotonProducer::GEDPhotonProducer ( const edm::ParameterSet &  ps)

Definition at line 68 of file GEDPhotonProducer.cc.

References barrelEcalHits_, candidateP4type_, checkHcalStatus_, conf_, endcapEcalHits_, edm::ParameterSet::exists(), flagsexclEB_, flagsexclEE_, edm::ParameterSet::getParameter(), hcalTowers_, highEt_, hOverEConeSize_, GEDPhotonProducer::RecoStepInfo::isFinal(), gedPhotons_cfi::isolationSumsCalculatorSet, edm::EDGetTokenT< T >::isUninitialized(), minR9Barrel_, minR9Endcap_, gedPhotons_cfi::mipVariableSet, dataset::name, pfCandidates_, pfEgammaCandidates_, gedPhotonSequence_cff::pfIsolCfg, phoChargedIsolationToken_, phoChargedPFPVIsoToken_, phoChargedWorstVtxGeomVetoIsoToken_, phoChargedWorstVtxIsoToken_, phoNeutralHadronIsolationToken_, phoPFECALClusIsolationToken_, phoPFHCALClusIsolationToken_, phoPhotonIsolationToken_, photonCollection_, photonCoreProducerT_, photonProducer_, photonProducerT_, cleanAndMergeSuperClusters_cfi::posCalcParameters, posCalculator_, preselCutValuesBarrel_, preselCutValuesEndcap_, preshowerHits_, recoStep_, runMIPTagger_, PhotonIsolationCalculator::setup(), PhotonMIPHaloTagger::setup(), severitiesexclEB_, severitiesexclEE_, AlCaHLTBitMon_QueryRunRegistry::string, thePhotonEnergyCorrector_, thePhotonIsolationCalculator_, thePhotonMIPHaloTagger_, usePrimaryVertex_, valueMapPFCandPhoton_, and vertexProducer_.

68  :
69  recoStep_(config.getParameter<std::string>("reconstructionStep")),
70  conf_(config)
71 {
72 
73  // use configuration file to setup input/output collection names
74  //
75  photonProducer_ = conf_.getParameter<edm::InputTag>("photonProducer");
76 
77  if ( recoStep_.isFinal() ) {
78  photonProducerT_ =
79  consumes<reco::PhotonCollection>(photonProducer_);
80  pfCandidates_ =
81  consumes<reco::PFCandidateCollection>(conf_.getParameter<edm::InputTag>("pfCandidates"));
82 
83  const edm::ParameterSet& pfIsolCfg = conf_.getParameter<edm::ParameterSet>("pfIsolCfg");
84  auto getVMToken = [&pfIsolCfg,this](const std::string& name){
85  return consumes<edm::ValueMap<float> >(pfIsolCfg.getParameter<edm::InputTag>(name));
86  };
87  phoChargedIsolationToken_ = getVMToken("chargedHadronIso");
88  phoNeutralHadronIsolationToken_ = getVMToken("neutralHadronIso");
89  phoPhotonIsolationToken_ = getVMToken("photonIso");
90  phoChargedWorstVtxIsoToken_ = getVMToken("chargedHadronWorstVtxIso");
91  phoChargedWorstVtxGeomVetoIsoToken_ = getVMToken("chargedHadronWorstVtxGeomVetoIso");
92  phoChargedPFPVIsoToken_ = getVMToken("chargedHadronPFPVIso");
93 
94  //OOT photons in legacy 80X re-miniAOD do not have PF cluster embeded into the reco object
95  //to preserve 80X behaviour
96  if(conf_.exists("pfECALClusIsolation")){
97  phoPFECALClusIsolationToken_ =
98  consumes<edm::ValueMap<float>>(conf_.getParameter<edm::InputTag>("pfECALClusIsolation"));
99  }
100  if(conf_.exists("pfHCALClusIsolation")){
101  phoPFHCALClusIsolationToken_ =
102  consumes<edm::ValueMap<float>>(conf_.getParameter<edm::InputTag>("pfHCALClusIsolation"));
103  }
104  } else {
105 
106  photonCoreProducerT_ =
107  consumes<reco::PhotonCoreCollection>(photonProducer_);
108 
109  }
110 
111  auto pfEg = conf_.getParameter<edm::InputTag>("pfEgammaCandidates");
112  if (not pfEg.label().empty())
113  pfEgammaCandidates_ =
114  consumes<reco::PFCandidateCollection>(pfEg);
115  barrelEcalHits_ =
116  consumes<EcalRecHitCollection>(conf_.getParameter<edm::InputTag>("barrelEcalHits"));
117  endcapEcalHits_ =
118  consumes<EcalRecHitCollection>(conf_.getParameter<edm::InputTag>("endcapEcalHits"));
119  preshowerHits_ =
120  consumes<EcalRecHitCollection>(conf_.getParameter<edm::InputTag>("preshowerHits"));
121  vertexProducer_ =
122  consumes<reco::VertexCollection>(conf_.getParameter<edm::InputTag>("primaryVertexProducer"));
123 
124  auto hcTow = conf_.getParameter<edm::InputTag>("hcalTowers");
125  if (not hcTow.label().empty())
126  hcalTowers_ =
127  consumes<CaloTowerCollection>(hcTow);
128  //
129  photonCollection_ = conf_.getParameter<std::string>("outputPhotonCollection");
130  hOverEConeSize_ = conf_.getParameter<double>("hOverEConeSize");
131  highEt_ = conf_.getParameter<double>("highEt");
132  // R9 value to decide converted/unconverted
133  minR9Barrel_ = conf_.getParameter<double>("minR9Barrel");
134  minR9Endcap_ = conf_.getParameter<double>("minR9Endcap");
135  usePrimaryVertex_ = conf_.getParameter<bool>("usePrimaryVertex");
136  runMIPTagger_ = conf_.getParameter<bool>("runMIPTagger");
137 
138  candidateP4type_ = config.getParameter<std::string>("candidateP4type") ;
139  valueMapPFCandPhoton_ = config.getParameter<std::string>("valueMapPhotons");
140 
141 
142  edm::ParameterSet posCalcParameters =
143  config.getParameter<edm::ParameterSet>("posCalcParameters");
144  posCalculator_ = PositionCalc(posCalcParameters);
145 
146 
147  //AA
148  //Flags and Severities to be excluded from photon calculations
149  const std::vector<std::string> flagnamesEB =
150  config.getParameter<std::vector<std::string> >("RecHitFlagToBeExcludedEB");
151 
152  const std::vector<std::string> flagnamesEE =
153  config.getParameter<std::vector<std::string> >("RecHitFlagToBeExcludedEE");
154 
155  flagsexclEB_=
156  StringToEnumValue<EcalRecHit::Flags>(flagnamesEB);
157 
158  flagsexclEE_=
159  StringToEnumValue<EcalRecHit::Flags>(flagnamesEE);
160 
161  const std::vector<std::string> severitynamesEB =
162  config.getParameter<std::vector<std::string> >("RecHitSeverityToBeExcludedEB");
163 
164  severitiesexclEB_=
165  StringToEnumValue<EcalSeverityLevel::SeverityLevel>(severitynamesEB);
166 
167  const std::vector<std::string> severitynamesEE =
168  config.getParameter<std::vector<std::string> >("RecHitSeverityToBeExcludedEE");
169 
170  severitiesexclEE_=
171  StringToEnumValue<EcalSeverityLevel::SeverityLevel>(severitynamesEE);
172 
173  thePhotonEnergyCorrector_ =
174  new PhotonEnergyCorrector(conf_, consumesCollector());
175 
176  //AA
177 
178  //
179 
180  // Parameters for the position calculation:
181  // std::map<std::string,double> providedParameters;
182  // providedParameters.insert(std::make_pair("LogWeighted",conf_.getParameter<bool>("posCalc_logweight")));
183  //providedParameters.insert(std::make_pair("T0_barl",conf_.getParameter<double>("posCalc_t0_barl")));
184  //providedParameters.insert(std::make_pair("T0_endc",conf_.getParameter<double>("posCalc_t0_endc")));
185  //providedParameters.insert(std::make_pair("T0_endcPresh",conf_.getParameter<double>("posCalc_t0_endcPresh")));
186  //providedParameters.insert(std::make_pair("W0",conf_.getParameter<double>("posCalc_w0")));
187  //providedParameters.insert(std::make_pair("X0",conf_.getParameter<double>("posCalc_x0")));
188  //posCalculator_ = PositionCalc(providedParameters);
189  // cut values for pre-selection
190  preselCutValuesBarrel_.push_back(conf_.getParameter<double>("minSCEtBarrel"));
191  preselCutValuesBarrel_.push_back(conf_.getParameter<double>("maxHoverEBarrel"));
192  preselCutValuesBarrel_.push_back(conf_.getParameter<double>("ecalRecHitSumEtOffsetBarrel"));
193  preselCutValuesBarrel_.push_back(conf_.getParameter<double>("ecalRecHitSumEtSlopeBarrel"));
194  preselCutValuesBarrel_.push_back(conf_.getParameter<double>("hcalTowerSumEtOffsetBarrel"));
195  preselCutValuesBarrel_.push_back(conf_.getParameter<double>("hcalTowerSumEtSlopeBarrel"));
196  preselCutValuesBarrel_.push_back(conf_.getParameter<double>("nTrackSolidConeBarrel"));
197  preselCutValuesBarrel_.push_back(conf_.getParameter<double>("nTrackHollowConeBarrel"));
198  preselCutValuesBarrel_.push_back(conf_.getParameter<double>("trackPtSumSolidConeBarrel"));
199  preselCutValuesBarrel_.push_back(conf_.getParameter<double>("trackPtSumHollowConeBarrel"));
200  preselCutValuesBarrel_.push_back(conf_.getParameter<double>("sigmaIetaIetaCutBarrel"));
201  //
202  preselCutValuesEndcap_.push_back(conf_.getParameter<double>("minSCEtEndcap"));
203  preselCutValuesEndcap_.push_back(conf_.getParameter<double>("maxHoverEEndcap"));
204  preselCutValuesEndcap_.push_back(conf_.getParameter<double>("ecalRecHitSumEtOffsetEndcap"));
205  preselCutValuesEndcap_.push_back(conf_.getParameter<double>("ecalRecHitSumEtSlopeEndcap"));
206  preselCutValuesEndcap_.push_back(conf_.getParameter<double>("hcalTowerSumEtOffsetEndcap"));
207  preselCutValuesEndcap_.push_back(conf_.getParameter<double>("hcalTowerSumEtSlopeEndcap"));
208  preselCutValuesEndcap_.push_back(conf_.getParameter<double>("nTrackSolidConeEndcap"));
209  preselCutValuesEndcap_.push_back(conf_.getParameter<double>("nTrackHollowConeEndcap"));
210  preselCutValuesEndcap_.push_back(conf_.getParameter<double>("trackPtSumSolidConeEndcap"));
211  preselCutValuesEndcap_.push_back(conf_.getParameter<double>("trackPtSumHollowConeEndcap"));
212  preselCutValuesEndcap_.push_back(conf_.getParameter<double>("sigmaIetaIetaCutEndcap"));
213  //
214 
215  //moved from beginRun to here, I dont see how this could cause harm as its just reading in the exactly same parameters each run
216  if ( !recoStep_.isFinal()){
217  thePhotonIsolationCalculator_ = new PhotonIsolationCalculator();
218  edm::ParameterSet isolationSumsCalculatorSet = conf_.getParameter<edm::ParameterSet>("isolationSumsCalculatorSet");
219  thePhotonIsolationCalculator_->setup(isolationSumsCalculatorSet, flagsexclEB_, flagsexclEE_, severitiesexclEB_, severitiesexclEE_,consumesCollector());
220  thePhotonMIPHaloTagger_ = new PhotonMIPHaloTagger();
221  edm::ParameterSet mipVariableSet = conf_.getParameter<edm::ParameterSet>("mipVariableSet");
222  thePhotonMIPHaloTagger_->setup(mipVariableSet,consumesCollector());
223 
224  }else{
225  thePhotonIsolationCalculator_=nullptr;
226  thePhotonMIPHaloTagger_=nullptr;
227  }
228 
229  checkHcalStatus_ = conf_.getParameter<bool>("checkHcalStatus");
230 
231  // Register the product
232  produces< reco::PhotonCollection >(photonCollection_);
233  if (not pfEgammaCandidates_.isUninitialized())
234  produces< edm::ValueMap<reco::PhotonRef> > (valueMapPFCandPhoton_);
235 
236 
237 }
GEDPhotonProducer::photonProducer_
edm::InputTag photonProducer_
Definition: GEDPhotonProducer.h:102
edm::ParameterSet::getParameter
T getParameter(std::string const &) const
GEDPhotonProducer::thePhotonEnergyCorrector_
PhotonEnergyCorrector * thePhotonEnergyCorrector_
Definition: GEDPhotonProducer.h:168
GEDPhotonProducer::thePhotonMIPHaloTagger_
PhotonMIPHaloTagger * thePhotonMIPHaloTagger_
Definition: GEDPhotonProducer.h:162
GEDPhotonProducer::pfCandidates_
edm::EDGetTokenT< reco::PFCandidateCollection > pfCandidates_
Definition: GEDPhotonProducer.h:110
GEDPhotonProducer::hcalTowers_
edm::EDGetTokenT< CaloTowerCollection > hcalTowers_
Definition: GEDPhotonProducer.h:111
gedPhotons_cfi.isolationSumsCalculatorSet
isolationSumsCalculatorSet
Definition: gedPhotons_cfi.py:34
AlCaHLTBitMon_QueryRunRegistry.string
string
Definition: AlCaHLTBitMon_QueryRunRegistry.py:256
gedPhotonSequence_cff.pfIsolCfg
pfIsolCfg
Definition: gedPhotonSequence_cff.py:27
GEDPhotonProducer::recoStep_
RecoStepInfo recoStep_
Definition: GEDPhotonProducer.h:150
edm::ParameterSet::exists
bool exists(std::string const &parameterName) const
checks if a parameter exists
Definition: ParameterSet.cc:674
GEDPhotonProducer::flagsexclEB_
std::vector< int > flagsexclEB_
Definition: GEDPhotonProducer.h:135
config
Definition: config.py:1
GEDPhotonProducer::endcapEcalHits_
edm::EDGetTokenT< EcalRecHitCollection > endcapEcalHits_
Definition: GEDPhotonProducer.h:107
GEDPhotonProducer::thePhotonIsolationCalculator_
PhotonIsolationCalculator * thePhotonIsolationCalculator_
Definition: GEDPhotonProducer.h:130
GEDPhotonProducer::severitiesexclEE_
std::vector< int > severitiesexclEE_
Definition: GEDPhotonProducer.h:138
GEDPhotonProducer::photonCollection_
std::string photonCollection_
Definition: GEDPhotonProducer.h:101
gedPhotons_cfi.mipVariableSet
mipVariableSet
Definition: gedPhotons_cfi.py:36
GEDPhotonProducer::phoChargedPFPVIsoToken_
edm::EDGetTokenT< edm::ValueMap< float > > phoChargedPFPVIsoToken_
Definition: GEDPhotonProducer.h:121
GEDPhotonProducer::barrelEcalHits_
edm::EDGetTokenT< EcalRecHitCollection > barrelEcalHits_
Definition: GEDPhotonProducer.h:106
GEDPhotonProducer::posCalculator_
PositionCalc posCalculator_
Definition: GEDPhotonProducer.h:154
GEDPhotonProducer::phoChargedWorstVtxGeomVetoIsoToken_
edm::EDGetTokenT< edm::ValueMap< float > > phoChargedWorstVtxGeomVetoIsoToken_
Definition: GEDPhotonProducer.h:120
GEDPhotonProducer::phoChargedIsolationToken_
edm::EDGetTokenT< edm::ValueMap< float > > phoChargedIsolationToken_
Definition: GEDPhotonProducer.h:116
PhotonIsolationCalculator::setup
void setup(const edm::ParameterSet &conf, std::vector< int > const &flagsEB_, std::vector< int > const &flagsEE_, std::vector< int > const &severitiesEB_, std::vector< int > const &severitiesEE_, edm::ConsumesCollector &&iC)
Definition: PhotonIsolationCalculator.cc:45
GEDPhotonProducer::flagsexclEE_
std::vector< int > flagsexclEE_
Definition: GEDPhotonProducer.h:136
GEDPhotonProducer::photonCoreProducerT_
edm::EDGetTokenT< reco::PhotonCoreCollection > photonCoreProducerT_
Definition: GEDPhotonProducer.h:104
GEDPhotonProducer::phoPFHCALClusIsolationToken_
edm::EDGetTokenT< edm::ValueMap< float > > phoPFHCALClusIsolationToken_
Definition: GEDPhotonProducer.h:124
PhotonMIPHaloTagger::setup
void setup(const edm::ParameterSet &conf, edm::ConsumesCollector &&iC)
Definition: PhotonMIPHaloTagger.cc:48
GEDPhotonProducer::pfEgammaCandidates_
edm::EDGetTokenT< reco::PFCandidateCollection > pfEgammaCandidates_
Definition: GEDPhotonProducer.h:109
GEDPhotonProducer::severitiesexclEB_
std::vector< int > severitiesexclEB_
Definition: GEDPhotonProducer.h:137
cleanAndMergeSuperClusters_cfi.posCalcParameters
posCalcParameters
Definition: cleanAndMergeSuperClusters_cfi.py:17
GEDPhotonProducer::phoNeutralHadronIsolationToken_
edm::EDGetTokenT< edm::ValueMap< float > > phoNeutralHadronIsolationToken_
Definition: GEDPhotonProducer.h:117
GEDPhotonProducer::phoChargedWorstVtxIsoToken_
edm::EDGetTokenT< edm::ValueMap< float > > phoChargedWorstVtxIsoToken_
Definition: GEDPhotonProducer.h:119
GEDPhotonProducer::preselCutValuesBarrel_
std::vector< double > preselCutValuesBarrel_
Definition: GEDPhotonProducer.h:164
GEDPhotonProducer::photonProducerT_
edm::EDGetTokenT< reco::PhotonCollection > photonProducerT_
Definition: GEDPhotonProducer.h:105
GEDPhotonProducer::preselCutValuesEndcap_
std::vector< double > preselCutValuesEndcap_
Definition: GEDPhotonProducer.h:165
GEDPhotonProducer::candidateP4type_
std::string candidateP4type_
Definition: GEDPhotonProducer.h:169
GEDPhotonProducer::phoPhotonIsolationToken_
edm::EDGetTokenT< edm::ValueMap< float > > phoPhotonIsolationToken_
Definition: GEDPhotonProducer.h:118
GEDPhotonProducer::phoPFECALClusIsolationToken_
edm::EDGetTokenT< edm::ValueMap< float > > phoPFECALClusIsolationToken_
Definition: GEDPhotonProducer.h:123
GEDPhotonProducer::conf_
edm::ParameterSet conf_
Definition: GEDPhotonProducer.h:153
edm::EDGetTokenT::isUninitialized
bool isUninitialized() const
Definition: EDGetToken.h:70
GEDPhotonProducer::vertexProducer_
edm::EDGetTokenT< reco::VertexCollection > vertexProducer_
Definition: GEDPhotonProducer.h:112
GEDPhotonProducer::preshowerHits_
edm::EDGetTokenT< EcalRecHitCollection > preshowerHits_
Definition: GEDPhotonProducer.h:108
GEDPhotonProducer::valueMapPFCandPhoton_
std::string valueMapPFCandPhoton_
Definition: GEDPhotonProducer.h:128
GEDPhotonProducer::~GEDPhotonProducer ( )
override

Definition at line 239 of file GEDPhotonProducer.cc.

References thePhotonEnergyCorrector_, thePhotonIsolationCalculator_, and thePhotonMIPHaloTagger_.

240 {
241  delete thePhotonEnergyCorrector_;
242  delete thePhotonIsolationCalculator_;
243  delete thePhotonMIPHaloTagger_;
244  //delete energyCorrectionF;
245 }
GEDPhotonProducer::thePhotonEnergyCorrector_
PhotonEnergyCorrector * thePhotonEnergyCorrector_
Definition: GEDPhotonProducer.h:168
GEDPhotonProducer::thePhotonMIPHaloTagger_
PhotonMIPHaloTagger * thePhotonMIPHaloTagger_
Definition: GEDPhotonProducer.h:162
GEDPhotonProducer::thePhotonIsolationCalculator_
PhotonIsolationCalculator * thePhotonIsolationCalculator_
Definition: GEDPhotonProducer.h:130

Member Function Documentation

void GEDPhotonProducer::beginRun ( edm::Run const &  r,
edm::EventSetup const &  es 
)
final

Definition at line 249 of file GEDPhotonProducer.cc.

References init, GEDPhotonProducer::RecoStepInfo::isFinal(), recoStep_, and thePhotonEnergyCorrector_.

249  {
250 
251  if ( !recoStep_.isFinal() ) {
252  thePhotonEnergyCorrector_ -> init(theEventSetup);
253  }
254 
255 }
GEDPhotonProducer::thePhotonEnergyCorrector_
PhotonEnergyCorrector * thePhotonEnergyCorrector_
Definition: GEDPhotonProducer.h:168
init
int init
Definition: HydjetWrapper.h:67
GEDPhotonProducer::recoStep_
RecoStepInfo recoStep_
Definition: GEDPhotonProducer.h:150
void GEDPhotonProducer::endRun ( edm::Run const &  r,
edm::EventSetup const &  theEventSetup 
)
final

Definition at line 257 of file GEDPhotonProducer.cc.

257  {
258 }
void GEDPhotonProducer::fillPhotonCollection ( edm::Event &  evt,
edm::EventSetup const &  es,
const edm::Handle< reco::PhotonCoreCollection > &  photonCoreHandle,
const CaloTopology *  topology,
const EcalRecHitCollection *  ecalBarrelHits,
const EcalRecHitCollection *  ecalEndcapHits,
const EcalRecHitCollection *  preshowerHits,
const edm::Handle< CaloTowerCollection > &  hcalTowersHandle,
const reco::VertexCollection &  pvVertices,
reco::PhotonCollection &  outputCollection,
int &  iSC 
)
private

fill shower shape block

fill extra shower shapes

fill full5x5 shower shape block

fill extra full5x5 shower shapes

get ecal photon specific corrected energy plus values from regressions and store them in the Photon

Pre-selection loose isolation cuts

Definition at line 501 of file GEDPhotonProducer.cc.

References Cluster2ndMoments::alpha, barrelEcalHits_, PhotonIsolationCalculator::calculate(), PhotonEnergyCorrector::calculate(), candidateP4type_, checkHcalStatus_, EcalClusterToolsT< noZS >::cluster2ndMoments(), EcalClusterToolsT< noZS >::covariances(), EcalClusterToolsT< noZS >::e1x3(), reco::Photon::ShowerShape::e1x3, EcalClusterToolsT< noZS >::e1x5(), reco::Photon::ShowerShape::e1x5, EcalClusterToolsT< noZS >::e2nd(), reco::Photon::ShowerShape::e2nd, EcalClusterToolsT< noZS >::e2x2(), reco::Photon::ShowerShape::e2x2, reco::Photon::ShowerShape::e2x5, EcalClusterToolsT< noZS >::e2x5Bottom(), reco::Photon::ShowerShape::e2x5Bottom, EcalClusterToolsT< noZS >::e2x5Left(), reco::Photon::ShowerShape::e2x5Left, EcalClusterToolsT< noZS >::e2x5Max(), reco::Photon::ShowerShape::e2x5Max, EcalClusterToolsT< noZS >::e2x5Right(), reco::Photon::ShowerShape::e2x5Right, EcalClusterToolsT< noZS >::e2x5Top(), reco::Photon::ShowerShape::e2x5Top, EcalClusterToolsT< noZS >::e3x3(), reco::Photon::ShowerShape::e3x3, EcalClusterToolsT< noZS >::e5x5(), reco::Photon::ShowerShape::e5x5, EcalClusterToolsT< noZS >::eBottom(), reco::Photon::ShowerShape::eBottom, reco::Photon::ecal_photons, EcalBarrel, EcalEndcap, ecalRecalibRecHit_cfi::ecalRecHit, reco::Photon::ShowerShape::effSigmaRR, EcalClusterToolsT< noZS >::eLeft(), reco::Photon::ShowerShape::eLeft, photonPostprocessing_cfi::eMax, EcalClusterToolsT< noZS >::eMax(), edm::SortedCollection< T, SORT >::end(), endcapEcalHits_, EcalClusterToolsT< noZS >::eRight(), reco::Photon::ShowerShape::eRight, EcalClusterToolsT< noZS >::eTop(), reco::Photon::ShowerShape::eTop, f, edm::SortedCollection< T, SORT >::find(), flagsexclEB_, flagsexclEE_, particleFlow_cfi::full5x5_sigmaIetaIeta, geometry, EgammaHadTower::getDepth1HcalESum(), EgammaHadTower::getDepth2HcalESum(), EgammaTowerIsolation::getTowerESum(), EgammaHadTower::hasActiveHcal(), reco::Photon::ShowerShape::hcalDepth1OverEcal, reco::Photon::ShowerShape::hcalDepth1OverEcalBc, reco::Photon::ShowerShape::hcalDepth2OverEcal, reco::Photon::ShowerShape::hcalDepth2OverEcalBc, hcalTowers_, reco::Photon::ShowerShape::hcalTowersBehindClusters, highEt_, hfClusterShapes_cfi::hits, hOverEConeSize_, createfilelist::int, reco::Photon::ShowerShape::invalidHcal, reco::Photon::FiducialFlags::isEE, edm::isFinite(), EcalTools::isHGCalDet(), edm::Ref< C, T, F >::isNonnull(), reco::Photon::SaturationInfo::isSeedSaturated, EcalClusterToolsT< noZS >::localCovariances(), reco::Photon::ShowerShape::maxEnergyXtal, reco::Photon::SaturationInfo::nSaturatedXtals, nullptr, p4, preselCutValuesBarrel_, preselCutValuesEndcap_, preshowerHits_, edm::Handle< T >::product(), edm::ESHandle< T >::product(), ptFast(), reco::Photon::regression1, reco::Photon::regression2, runMIPTagger_, reco::CaloCluster::seed(), EgammaHadTower::setTowerCollection(), severitiesexclEB_, severitiesexclEE_, electronIdCutBased_cfi::sigmaEtaEta, reco::Photon::ShowerShape::sigmaEtaEta, reco::Photon::ShowerShape::sigmaIetaIeta, reco::Photon::ShowerShape::sigmaIetaIphi, reco::Photon::ShowerShape::sigmaIphiIphi, Cluster2ndMoments::sMaj, reco::Photon::ShowerShape::smAlpha, Cluster2ndMoments::sMin, reco::Photon::ShowerShape::smMajor, reco::Photon::ShowerShape::smMinor, mathSSE::sqrt(), theCaloGeom_, thePhotonEnergyCorrector_, thePhotonIsolationCalculator_, thePhotonMIPHaloTagger_, EgammaHadTower::towersOf(), and badGlobalMuonTaggersAOD_cff::vtx.

Referenced by produce().

510  {
511 
512 
513  const CaloGeometry* geometry = theCaloGeom_.product();
514  const EcalRecHitCollection* hits = nullptr ;
515  std::vector<double> preselCutValues;
516  std::vector<int> flags_, severitiesexcl_;
517 
518  for(unsigned int lSC=0; lSC < photonCoreHandle->size(); lSC++) {
519 
520  reco::PhotonCoreRef coreRef(reco::PhotonCoreRef(photonCoreHandle, lSC));
521  reco::SuperClusterRef parentSCRef = coreRef->parentSuperCluster();
522  reco::SuperClusterRef scRef=coreRef->superCluster();
523 
524 
525 
526  // const reco::SuperCluster* pClus=&(*scRef);
527  iSC++;
528 
529  DetId::Detector thedet = scRef->seed()->hitsAndFractions()[0].first.det();
530  int subdet = scRef->seed()->hitsAndFractions()[0].first.subdetId();
531  if (subdet==EcalBarrel) {
532  preselCutValues = preselCutValuesBarrel_;
533  hits = ecalBarrelHits;
534  flags_ = flagsexclEB_;
535  severitiesexcl_ = severitiesexclEB_;
536  } else if (subdet==EcalEndcap) {
537  preselCutValues = preselCutValuesEndcap_;
538  hits = ecalEndcapHits;
539  flags_ = flagsexclEE_;
540  severitiesexcl_ = severitiesexclEE_;
541  } else if ( EcalTools::isHGCalDet(thedet) ) {
542  preselCutValues = preselCutValuesEndcap_;
543  hits = nullptr;
544  flags_ = flagsexclEE_;
545  severitiesexcl_ = severitiesexclEE_;
546  } else {
547  edm::LogWarning("")<<"GEDPhotonProducer: do not know if it is a barrel or endcap SuperCluster: " << thedet << ' ' << subdet;
548  }
549 
550 
551 
552 
553  // SC energy preselection
554  if (parentSCRef.isNonnull() &&
555  ptFast(parentSCRef->energy(),parentSCRef->position(),math::XYZPoint(0,0,0)) <= preselCutValues[0] ) continue;
556  // calculate HoE
557 
558  double HoE1,HoE2;
559  HoE1=HoE2=0.;
560 
561  std::vector<CaloTowerDetId> TowersBehindClus;
562  float hcalDepth1OverEcalBc,hcalDepth2OverEcalBc;
563  hcalDepth1OverEcalBc=hcalDepth2OverEcalBc=0.f;
564  bool invalidHcal = false;
565 
566  if (not hcalTowers_.isUninitialized()) {
567  const CaloTowerCollection* hcalTowersColl = hcalTowersHandle.product();
568  EgammaTowerIsolation towerIso1(hOverEConeSize_,0.,0.,1,hcalTowersColl) ;
569  EgammaTowerIsolation towerIso2(hOverEConeSize_,0.,0.,2,hcalTowersColl) ;
570  HoE1=towerIso1.getTowerESum(&(*scRef))/scRef->energy();
571  HoE2=towerIso2.getTowerESum(&(*scRef))/scRef->energy();
572 
573  EgammaHadTower towerIsoBehindClus(es);
574  towerIsoBehindClus.setTowerCollection(hcalTowersHandle.product());
575  TowersBehindClus = towerIsoBehindClus.towersOf(*scRef);
576  hcalDepth1OverEcalBc = towerIsoBehindClus.getDepth1HcalESum(TowersBehindClus)/scRef->energy();
577  hcalDepth2OverEcalBc = towerIsoBehindClus.getDepth2HcalESum(TowersBehindClus)/scRef->energy();
578 
579  if (checkHcalStatus_ && hcalDepth1OverEcalBc == 0 && hcalDepth2OverEcalBc == 0) {
580  invalidHcal = !towerIsoBehindClus.hasActiveHcal(TowersBehindClus);
581  }
582  }
583 
584  // std::cout << " GEDPhotonProducer calculation of HoE with towers in a cone " << HoE1 << " " << HoE2 << std::endl;
585  //std::cout << " GEDPhotonProducer calcualtion of HoE with towers behind the BCs " << hcalDepth1OverEcalBc << " " << hcalDepth2OverEcalBc << std::endl;
586 
587  float maxXtal = ( hits != nullptr ? EcalClusterTools::eMax( *(scRef->seed()), &(*hits) ) : 0.f );
588  //AA
589  //Change these to consider severity level of hits
590  float e1x5 = ( hits != nullptr ? EcalClusterTools::e1x5( *(scRef->seed()), &(*hits), &(*topology)) : 0.f );
591  float e2x5 = ( hits != nullptr ? EcalClusterTools::e2x5Max( *(scRef->seed()), &(*hits), &(*topology)) : 0.f );
592  float e3x3 = ( hits != nullptr ? EcalClusterTools::e3x3( *(scRef->seed()), &(*hits), &(*topology)) : 0.f );
593  float e5x5 = ( hits != nullptr ? EcalClusterTools::e5x5( *(scRef->seed()), &(*hits), &(*topology)) : 0.f );
594  std::vector<float> cov = ( hits != nullptr ? EcalClusterTools::covariances( *(scRef->seed()), &(*hits), &(*topology), geometry) : std::vector<float>( {0.f,0.f,0.f} ) );
595  std::vector<float> locCov = ( hits != nullptr ? EcalClusterTools::localCovariances( *(scRef->seed()), &(*hits), &(*topology)) : std::vector<float>( {0.f,0.f,0.f} ) );
596 
597  float sigmaEtaEta = sqrt(cov[0]);
598  float sigmaIetaIeta = sqrt(locCov[0]);
599 
600  float full5x5_maxXtal = ( hits != nullptr ? noZS::EcalClusterTools::eMax( *(scRef->seed()), &(*hits) ) : 0.f );
601  //AA
602  //Change these to consider severity level of hits
603  float full5x5_e1x5 = ( hits != nullptr ? noZS::EcalClusterTools::e1x5( *(scRef->seed()), &(*hits), &(*topology)) : 0.f );
604  float full5x5_e2x5 = ( hits != nullptr ? noZS::EcalClusterTools::e2x5Max( *(scRef->seed()), &(*hits), &(*topology)) : 0.f );
605  float full5x5_e3x3 = ( hits != nullptr ? noZS::EcalClusterTools::e3x3( *(scRef->seed()), &(*hits), &(*topology)) : 0.f );
606  float full5x5_e5x5 = ( hits != nullptr ? noZS::EcalClusterTools::e5x5( *(scRef->seed()), &(*hits), &(*topology)) : 0.f );
607  std::vector<float> full5x5_cov = ( hits != nullptr ? noZS::EcalClusterTools::covariances( *(scRef->seed()), &(*hits), &(*topology), geometry) : std::vector<float>( {0.f,0.f,0.f} ) );
608  std::vector<float> full5x5_locCov = ( hits != nullptr ? noZS::EcalClusterTools::localCovariances( *(scRef->seed()), &(*hits), &(*topology)) : std::vector<float>( {0.f,0.f,0.f} ) );
609 
610  float full5x5_sigmaEtaEta = sqrt(full5x5_cov[0]);
611  float full5x5_sigmaIetaIeta = sqrt(full5x5_locCov[0]);
612 
613  // compute position of ECAL shower
614  math::XYZPoint caloPosition = scRef->position();
615 
616 
618  double photonEnergy=1.;
619  math::XYZPoint vtx(0.,0.,0.);
620  if (!vertexCollection.empty()) vtx = vertexCollection.begin()->position();
621  // compute momentum vector of photon from primary vertex and cluster position
622  math::XYZVector direction = caloPosition - vtx;
623  //math::XYZVector momentum = direction.unit() * photonEnergy ;
624  math::XYZVector momentum = direction.unit() ;
625 
626  // Create dummy candidate with unit momentum and zero energy to allow setting of all variables. The energy is set for last.
627  math::XYZTLorentzVectorD p4(momentum.x(), momentum.y(), momentum.z(), photonEnergy );
628  reco::Photon newCandidate(p4, caloPosition, coreRef, vtx);
629 
630  //std::cout << " standard p4 before " << newCandidate.p4() << " energy " << newCandidate.energy() << std::endl;
631  //std::cout << " type " <<newCandidate.getCandidateP4type() << " standard p4 after " << newCandidate.p4() << " energy " << newCandidate.energy() << std::endl;
632 
633  // Calculate fiducial flags and isolation variable. Blocked are filled from the isolationCalculator
634  reco::Photon::FiducialFlags fiducialFlags;
635  reco::Photon::IsolationVariables isolVarR03, isolVarR04;
636  if( !EcalTools::isHGCalDet(thedet) ) {
637  thePhotonIsolationCalculator_->calculate( &newCandidate,evt,es,fiducialFlags,isolVarR04, isolVarR03);
638  }
639  newCandidate.setFiducialVolumeFlags( fiducialFlags );
640  newCandidate.setIsolationVariables(isolVarR04, isolVarR03 );
641 
642 
644  reco::Photon::ShowerShape showerShape;
645  showerShape.e1x5= e1x5;
646  showerShape.e2x5= e2x5;
647  showerShape.e3x3= e3x3;
648  showerShape.e5x5= e5x5;
649  showerShape.maxEnergyXtal = maxXtal;
650  showerShape.sigmaEtaEta = sigmaEtaEta;
651  showerShape.sigmaIetaIeta = sigmaIetaIeta;
652  showerShape.hcalDepth1OverEcal = HoE1;
653  showerShape.hcalDepth2OverEcal = HoE2;
654  showerShape.hcalDepth1OverEcalBc = hcalDepth1OverEcalBc;
655  showerShape.hcalDepth2OverEcalBc = hcalDepth2OverEcalBc;
656  showerShape.hcalTowersBehindClusters = TowersBehindClus;
657  showerShape.invalidHcal = invalidHcal;
659  const float spp = (!edm::isFinite(locCov[2]) ? 0. : sqrt(locCov[2]));
660  const float sep = locCov[1];
661  showerShape.sigmaIetaIphi = sep;
662  showerShape.sigmaIphiIphi = spp;
663  showerShape.e2nd = ( hits != nullptr ? EcalClusterTools::e2nd(*(scRef->seed()),&(*hits)) : 0.f );
664  showerShape.eTop = ( hits != nullptr ? EcalClusterTools::eTop(*(scRef->seed()), &(*hits), &(*topology)) : 0.f );
665  showerShape.eLeft = ( hits != nullptr ? EcalClusterTools::eLeft(*(scRef->seed()), &(*hits), &(*topology)) : 0.f );
666  showerShape.eRight = ( hits != nullptr ? EcalClusterTools::eRight(*(scRef->seed()), &(*hits), &(*topology)) : 0.f );
667  showerShape.eBottom = ( hits != nullptr ? EcalClusterTools::eBottom(*(scRef->seed()), &(*hits), &(*topology)) : 0.f );
668  showerShape.e1x3 = ( hits != nullptr ? EcalClusterTools::e1x3(*(scRef->seed()), &(*hits), &(*topology)) : 0.f );
669  showerShape.e2x2 = ( hits != nullptr ? EcalClusterTools::e2x2(*(scRef->seed()), &(*hits), &(*topology)) : 0.f );
670  showerShape.e2x5Max = ( hits != nullptr ? EcalClusterTools::e2x5Max(*(scRef->seed()), &(*hits), &(*topology)) : 0.f );
671  showerShape.e2x5Left = ( hits != nullptr ? EcalClusterTools::e2x5Left(*(scRef->seed()), &(*hits), &(*topology)) : 0.f );
672  showerShape.e2x5Right = ( hits != nullptr ? EcalClusterTools::e2x5Right(*(scRef->seed()), &(*hits), &(*topology)) : 0.f );
673  showerShape.e2x5Top = ( hits != nullptr ? EcalClusterTools::e2x5Top(*(scRef->seed()), &(*hits), &(*topology)) : 0.f );
674  showerShape.e2x5Bottom = ( hits != nullptr ? EcalClusterTools::e2x5Bottom(*(scRef->seed()), &(*hits), &(*topology)) : 0.f );
675  if(hits){
676  Cluster2ndMoments clus2ndMoments = EcalClusterTools::cluster2ndMoments(*(scRef->seed()),*hits);
677  showerShape.smMajor = clus2ndMoments.sMaj;
678  showerShape.smMinor = clus2ndMoments.sMin;
679  showerShape.smAlpha = clus2ndMoments.alpha;
680  }else{
681  showerShape.smMajor = 0.f;
682  showerShape.smMinor = 0.f;
683  showerShape.smAlpha = 0.f;
684  }
685 
686 
687  // fill preshower shapes
688  EcalClusterLazyTools toolsforES(evt, es, barrelEcalHits_, endcapEcalHits_, preshowerHits_);
689  const float sigmaRR = toolsforES.eseffsirir( *scRef );
690  showerShape.effSigmaRR = sigmaRR;
691  newCandidate.setShowerShapeVariables ( showerShape );
692 
693  reco::Photon::SaturationInfo saturationInfo;
694  const reco::CaloCluster& seedCluster = *(scRef->seed()) ;
695  DetId seedXtalId = seedCluster.seed();
696  int nSaturatedXtals = 0;
697  bool isSeedSaturated = false;
698  if (hits != nullptr) {
699  const auto hitsAndFractions = scRef->hitsAndFractions();
700  for (auto&& hitFractionPair : hitsAndFractions) {
701  auto&& ecalRecHit = hits->find(hitFractionPair.first);
702  if (ecalRecHit == hits->end()) continue;
703  if (ecalRecHit->checkFlag(EcalRecHit::Flags::kSaturated)) {
704  nSaturatedXtals++;
705  if (seedXtalId == ecalRecHit->detid())
706  isSeedSaturated = true;
707  }
708  }
709  }
710  saturationInfo.nSaturatedXtals = nSaturatedXtals;
711  saturationInfo.isSeedSaturated = isSeedSaturated;
712  newCandidate.setSaturationInfo(saturationInfo);
713 
715  reco::Photon::ShowerShape full5x5_showerShape;
716  full5x5_showerShape.e1x5= full5x5_e1x5;
717  full5x5_showerShape.e2x5= full5x5_e2x5;
718  full5x5_showerShape.e3x3= full5x5_e3x3;
719  full5x5_showerShape.e5x5= full5x5_e5x5;
720  full5x5_showerShape.maxEnergyXtal = full5x5_maxXtal;
721  full5x5_showerShape.sigmaEtaEta = full5x5_sigmaEtaEta;
722  full5x5_showerShape.sigmaIetaIeta = full5x5_sigmaIetaIeta;
724  const float full5x5_spp = (!edm::isFinite(full5x5_locCov[2]) ? 0. : sqrt(full5x5_locCov[2]));
725  const float full5x5_sep = full5x5_locCov[1];
726  full5x5_showerShape.sigmaIetaIphi = full5x5_sep;
727  full5x5_showerShape.sigmaIphiIphi = full5x5_spp;
728  full5x5_showerShape.e2nd = ( hits != nullptr ? noZS::EcalClusterTools::e2nd(*(scRef->seed()),&(*hits)) : 0.f );
729  full5x5_showerShape.eTop = ( hits != nullptr ? noZS::EcalClusterTools::eTop(*(scRef->seed()), &(*hits), &(*topology)) : 0.f );
730  full5x5_showerShape.eLeft = ( hits != nullptr ? noZS::EcalClusterTools::eLeft(*(scRef->seed()), &(*hits), &(*topology)) : 0.f );
731  full5x5_showerShape.eRight = ( hits != nullptr ? noZS::EcalClusterTools::eRight(*(scRef->seed()), &(*hits), &(*topology)) : 0.f );
732  full5x5_showerShape.eBottom = ( hits != nullptr ? noZS::EcalClusterTools::eBottom(*(scRef->seed()), &(*hits), &(*topology)) : 0.f );
733  full5x5_showerShape.e1x3 = ( hits != nullptr ? noZS::EcalClusterTools::e1x3(*(scRef->seed()), &(*hits), &(*topology)) : 0.f );
734  full5x5_showerShape.e2x2 = ( hits != nullptr ? noZS::EcalClusterTools::e2x2(*(scRef->seed()), &(*hits), &(*topology)) : 0.f );
735  full5x5_showerShape.e2x5Max = ( hits != nullptr ? noZS::EcalClusterTools::e2x5Max(*(scRef->seed()), &(*hits), &(*topology)) : 0.f );
736  full5x5_showerShape.e2x5Left = ( hits != nullptr ? noZS::EcalClusterTools::e2x5Left(*(scRef->seed()), &(*hits), &(*topology)) : 0.f );
737  full5x5_showerShape.e2x5Right = ( hits != nullptr ? noZS::EcalClusterTools::e2x5Right(*(scRef->seed()), &(*hits), &(*topology)) : 0.f );
738  full5x5_showerShape.e2x5Top = ( hits != nullptr ? noZS::EcalClusterTools::e2x5Top(*(scRef->seed()), &(*hits), &(*topology)) : 0.f );
739  full5x5_showerShape.e2x5Bottom = ( hits != nullptr ? noZS::EcalClusterTools::e2x5Bottom(*(scRef->seed()), &(*hits), &(*topology)) : 0.f );
740  if(hits){
741  Cluster2ndMoments clus2ndMoments = noZS::EcalClusterTools::cluster2ndMoments(*(scRef->seed()),*hits);
742  full5x5_showerShape.smMajor = clus2ndMoments.sMaj;
743  full5x5_showerShape.smMinor = clus2ndMoments.sMin;
744  full5x5_showerShape.smAlpha = clus2ndMoments.alpha;
745  }else{
746  full5x5_showerShape.smMajor = 0.f;
747  full5x5_showerShape.smMinor = 0.f;
748  full5x5_showerShape.smAlpha = 0.f;
749  }
750  // fill preshower shapes
751  full5x5_showerShape.effSigmaRR = sigmaRR;
752  newCandidate.full5x5_setShowerShapeVariables ( full5x5_showerShape );
753 
754 
755 
758  // Photon candidate takes by default (set in photons_cfi.py)
759  // a 4-momentum derived from the ecal photon-specific corrections.
760  if( !EcalTools::isHGCalDet(thedet) ) {
761  thePhotonEnergyCorrector_->calculate(evt, newCandidate, subdet, vertexCollection, es);
762  if ( candidateP4type_ == "fromEcalEnergy") {
763  newCandidate.setP4( newCandidate.p4(reco::Photon::ecal_photons) );
764  newCandidate.setCandidateP4type(reco::Photon::ecal_photons);
765  } else if ( candidateP4type_ == "fromRegression1") {
766  newCandidate.setP4( newCandidate.p4(reco::Photon::regression1) );
767  newCandidate.setCandidateP4type(reco::Photon::regression1);
768  } else if ( candidateP4type_ == "fromRegression2") {
769  newCandidate.setP4( newCandidate.p4(reco::Photon::regression2) );
770  newCandidate.setCandidateP4type(reco::Photon::regression2);
771  } else if ( candidateP4type_ == "fromRefinedSCRegression" ) {
772  newCandidate.setP4( newCandidate.p4(reco::Photon::regression2) );
773  newCandidate.setCandidateP4type(reco::Photon::regression2);
774  }
775  } else {
776  math::XYZVector gamma_momentum = direction.unit() * scRef->energy();
777  math::XYZTLorentzVectorD p4(gamma_momentum.x(),
778  gamma_momentum.y(),
779  gamma_momentum.z(),
780  scRef->energy());
781  newCandidate.setP4(p4);
782  newCandidate.setCandidateP4type(reco::Photon::ecal_photons);
783  // Make it an EE photon
784  reco::Photon::FiducialFlags fiducialFlags;
785  fiducialFlags.isEE = true;
786  newCandidate.setFiducialVolumeFlags(fiducialFlags);
787  }
788 
789  // std::cout << " final p4 " << newCandidate.p4() << " energy " << newCandidate.energy() << std::endl;
790 
791 
792  // std::cout << " GEDPhotonProducer from candidate HoE with towers in a cone " << newCandidate.hadronicOverEm() << " " << newCandidate.hadronicDepth1OverEm() << " " << newCandidate.hadronicDepth2OverEm() << std::endl;
793  // std::cout << " GEDPhotonProducer from candidate of HoE with towers behind the BCs " << newCandidate.hadTowOverEm() << " " << newCandidate.hadTowDepth1OverEm() << " " << newCandidate.hadTowDepth2OverEm() << std::endl;
794 
795 
796  // fill MIP Vairables for Halo: Block for MIP are filled from PhotonMIPHaloTagger
797  reco::Photon::MIPVariables mipVar ;
798  if(subdet==EcalBarrel && runMIPTagger_ )
799  {
800 
801  thePhotonMIPHaloTagger_-> MIPcalculate( &newCandidate,evt,es,mipVar);
802  newCandidate.setMIPVariables(mipVar);
803  }
804 
805 
806 
808  bool isLooseEM=true;
809  if ( newCandidate.pt() < highEt_) {
810  if ( newCandidate.hadronicOverEm() >= preselCutValues[1] ) isLooseEM=false;
811  if ( newCandidate.ecalRecHitSumEtConeDR04() > preselCutValues[2]+ preselCutValues[3]*newCandidate.pt() ) isLooseEM=false;
812  if ( newCandidate.hcalTowerSumEtConeDR04() > preselCutValues[4]+ preselCutValues[5]*newCandidate.pt() ) isLooseEM=false;
813  if ( newCandidate.nTrkSolidConeDR04() > int(preselCutValues[6]) ) isLooseEM=false;
814  if ( newCandidate.nTrkHollowConeDR04() > int(preselCutValues[7]) ) isLooseEM=false;
815  if ( newCandidate.trkSumPtSolidConeDR04() > preselCutValues[8] ) isLooseEM=false;
816  if ( newCandidate.trkSumPtHollowConeDR04() > preselCutValues[9] ) isLooseEM=false;
817  if ( newCandidate.sigmaIetaIeta() > preselCutValues[10] ) isLooseEM=false;
818  }
819 
820 
821 
822  if ( isLooseEM)
823  outputPhotonCollection.push_back(newCandidate);
824 
825 
826  }
827 }
EcalTools::isHGCalDet
static bool isHGCalDet(DetId::Detector thedet)
identify HGCal cells
Definition: EcalTools.h:54
GEDPhotonProducer::thePhotonEnergyCorrector_
PhotonEnergyCorrector * thePhotonEnergyCorrector_
Definition: GEDPhotonProducer.h:168
photonPostprocessing_cfi.eMax
eMax
Definition: photonPostprocessing_cfi.py:22
GEDPhotonProducer::thePhotonMIPHaloTagger_
PhotonMIPHaloTagger * thePhotonMIPHaloTagger_
Definition: GEDPhotonProducer.h:162
PhotonIsolationCalculator::calculate
void calculate(const reco::Photon *, const edm::Event &, const edm::EventSetup &es, reco::Photon::FiducialFlags &phofid, reco::Photon::IsolationVariables &phoisolR03, reco::Photon::IsolationVariables &phoisolR04) const
Definition: PhotonIsolationCalculator.cc:182
edm::Ref::isNonnull
bool isNonnull() const
Checks for non-null.
Definition: Ref.h:251
GEDPhotonProducer::theCaloGeom_
edm::ESHandle< CaloGeometry > theCaloGeom_
Definition: GEDPhotonProducer.h:156
EcalClusterToolsT::e2x5Bottom
static float e2x5Bottom(const reco::BasicCluster &cluster, const EcalRecHitCollection *recHits, const CaloTopology *topology)
Definition: EcalClusterTools.h:464
reco::Photon::ShowerShape::hcalTowersBehindClusters
std::vector< CaloTowerDetId > hcalTowersBehindClusters
Definition: Photon.h:152
math::XYZTLorentzVectorD
ROOT::Math::LorentzVector< ROOT::Math::PxPyPzE4D< double > > XYZTLorentzVectorD
Lorentz vector with cylindrical internal representation using pseudorapidity.
Definition: LorentzVector.h:14
EcalClusterToolsT::covariances
static std::vector< float > covariances(const reco::BasicCluster &cluster, const EcalRecHitCollection *recHits, const CaloTopology *topology, const CaloGeometry *geometry, float w0=4.7)
Definition: EcalClusterTools.h:789
GEDPhotonProducer::hcalTowers_
edm::EDGetTokenT< CaloTowerCollection > hcalTowers_
Definition: GEDPhotonProducer.h:111
EcalClusterToolsT::eMax
static float eMax(const reco::BasicCluster &cluster, const EcalRecHitCollection *recHits)
Definition: EcalClusterTools.h:423
ecalRecalibRecHit_cfi.ecalRecHit
ecalRecHit
Definition: ecalRecalibRecHit_cfi.py:4
hfClusterShapes_cfi.hits
hits
Definition: hfClusterShapes_cfi.py:5
nullptr
#define nullptr
Definition: GCC11Compatibility.h:37
GEDPhotonProducer::flagsexclEB_
std::vector< int > flagsexclEB_
Definition: GEDPhotonProducer.h:135
particleFlow_cfi.full5x5_sigmaIetaIeta
full5x5_sigmaIetaIeta
Definition: particleFlow_cfi.py:77
GEDPhotonProducer::endcapEcalHits_
edm::EDGetTokenT< EcalRecHitCollection > endcapEcalHits_
Definition: GEDPhotonProducer.h:107
edm::isFinite
constexpr bool isFinite(T x)
GEDPhotonProducer::thePhotonIsolationCalculator_
PhotonIsolationCalculator * thePhotonIsolationCalculator_
Definition: GEDPhotonProducer.h:130
GEDPhotonProducer::severitiesexclEE_
std::vector< int > severitiesexclEE_
Definition: GEDPhotonProducer.h:138
ptFast
double ptFast(const double energy, const math::XYZPoint &position, const math::XYZPoint &origin)
Definition: SuperClusterImporter.cc:14
EcalClusterToolsT::e2x5Top
static float e2x5Top(const reco::BasicCluster &cluster, const EcalRecHitCollection *recHits, const CaloTopology *topology)
Definition: EcalClusterTools.h:457
gedPhotons_cfi.outputPhotonCollection
outputPhotonCollection
Definition: gedPhotons_cfi.py:30
GEDPhotonProducer::barrelEcalHits_
edm::EDGetTokenT< EcalRecHitCollection > barrelEcalHits_
Definition: GEDPhotonProducer.h:106
EcalClusterToolsT::e2x2
static float e2x2(const reco::BasicCluster &cluster, const EcalRecHitCollection *recHits, const CaloTopology *topology)
Definition: EcalClusterTools.h:368
mathSSE::sqrt
T sqrt(T t)
Definition: SSEVec.h:18
p4
double p4[4]
Definition: TauolaWrapper.h:92
EcalClusterToolsT::e2nd
static float e2nd(const reco::BasicCluster &cluster, const EcalRecHitCollection *recHits)
Definition: EcalClusterTools.h:429
GEDPhotonProducer::flagsexclEE_
std::vector< int > flagsexclEE_
Definition: GEDPhotonProducer.h:136
EcalClusterToolsT::e2x5Max
static float e2x5Max(const reco::BasicCluster &cluster, const EcalRecHitCollection *recHits, const CaloTopology *topology)
Definition: EcalClusterTools.h:473
f
double f[11][100]
Definition: MuScleFitUtils.cc:78
EcalClusterToolsT::eBottom
static float eBottom(const reco::BasicCluster &cluster, const EcalRecHitCollection *recHits, const CaloTopology *topology)
Definition: EcalClusterTools.h:538
EcalClusterToolsT::e2x5Right
static float e2x5Right(const reco::BasicCluster &cluster, const EcalRecHitCollection *recHits, const CaloTopology *topology)
Definition: EcalClusterTools.h:443
GEDPhotonProducer::severitiesexclEB_
std::vector< int > severitiesexclEB_
Definition: GEDPhotonProducer.h:137
edm::SortedCollection::end
const_iterator end() const
Definition: SortedCollection.h:310
EcalClusterToolsT::cluster2ndMoments
static Cluster2ndMoments cluster2ndMoments(const reco::BasicCluster &basicCluster, const EcalRecHitCollection &recHits, double phiCorrectionFactor=0.8, double w0=4.7, bool useLogWeights=true)
Definition: EcalClusterTools.h:1231
electronIdCutBased_cfi.sigmaEtaEta
sigmaEtaEta
Definition: electronIdCutBased_cfi.py:125
EcalClusterToolsT::eTop
static float eTop(const reco::BasicCluster &cluster, const EcalRecHitCollection *recHits, const CaloTopology *topology)
Definition: EcalClusterTools.h:531
DetId
Definition: DetId.h:18
PhotonEnergyCorrector::calculate
void calculate(edm::Event &evt, reco::Photon &, int subdet, const reco::VertexCollection &vtxcol, const edm::EventSetup &iSetup)
Definition: PhotonEnergyCorrector.cc:92
reco::CaloCluster::seed
DetId seed() const
return DetId of seed
Definition: CaloCluster.h:207
edm::Handle::product
T const * product() const
Definition: Handle.h:74
GEDPhotonProducer::preselCutValuesBarrel_
std::vector< double > preselCutValuesBarrel_
Definition: GEDPhotonProducer.h:164
math::XYZVector
XYZVectorD XYZVector
spatial vector with cartesian internal representation
Definition: Vector3D.h:30
math::XYZPoint
XYZPointD XYZPoint
point in space with cartesian internal representation
Definition: Point3D.h:12
DetId::Detector
Detector
Definition: DetId.h:26
GEDPhotonProducer::preselCutValuesEndcap_
std::vector< double > preselCutValuesEndcap_
Definition: GEDPhotonProducer.h:165
EcalClusterToolsT::e2x5Left
static float e2x5Left(const reco::BasicCluster &cluster, const EcalRecHitCollection *recHits, const CaloTopology *topology)
Definition: EcalClusterTools.h:450
GEDPhotonProducer::candidateP4type_
std::string candidateP4type_
Definition: GEDPhotonProducer.h:169
particleFlowSuperClusterECAL_cfi.vertexCollection
vertexCollection
Definition: particleFlowSuperClusterECAL_cfi.py:41
geometry
ESHandle< TrackerGeometry > geometry
Definition: TkLasBeamFitter.cc:176
edm::SortedCollection::find
iterator find(key_type k)
Definition: SortedCollection.h:271
EcalClusterToolsT::eRight
static float eRight(const reco::BasicCluster &cluster, const EcalRecHitCollection *recHits, const CaloTopology *topology)
Definition: EcalClusterTools.h:524
EcalClusterToolsT::e3x3
static float e3x3(const reco::BasicCluster &cluster, const EcalRecHitCollection *recHits, const CaloTopology *topology)
Definition: EcalClusterTools.h:391
badGlobalMuonTaggersAOD_cff.vtx
vtx
Definition: badGlobalMuonTaggersAOD_cff.py:5
EcalClusterToolsT::e1x3
static float e1x3(const reco::BasicCluster &cluster, const EcalRecHitCollection *recHits, const CaloTopology *topology)
Definition: EcalClusterTools.h:503
edm::ESHandle::product
T const * product() const
Definition: ESHandle.h:86
EcalClusterToolsT::eLeft
static float eLeft(const reco::BasicCluster &cluster, const EcalRecHitCollection *recHits, const CaloTopology *topology)
Definition: EcalClusterTools.h:517
EcalClusterToolsT::localCovariances
static std::vector< float > localCovariances(const reco::BasicCluster &cluster, const EcalRecHitCollection *recHits, const CaloTopology *topology, float w0=4.7)
Definition: EcalClusterTools.h:858
GEDPhotonProducer::preshowerHits_
edm::EDGetTokenT< EcalRecHitCollection > preshowerHits_
Definition: GEDPhotonProducer.h:108
EcalClusterToolsT::e1x5
static float e1x5(const reco::BasicCluster &cluster, const EcalRecHitCollection *recHits, const CaloTopology *topology)
Definition: EcalClusterTools.h:489
EcalClusterToolsT::e5x5
static float e5x5(const reco::BasicCluster &cluster, const EcalRecHitCollection *recHits, const CaloTopology *topology)
Definition: EcalClusterTools.h:409
void GEDPhotonProducer::fillPhotonCollection ( edm::Event &  evt,
edm::EventSetup const &  es,
const edm::Handle< reco::PhotonCollection > &  photonHandle,
const edm::Handle< reco::PFCandidateCollection >  pfCandidateHandle,
const edm::Handle< reco::PFCandidateCollection >  pfEGCandidateHandle,
edm::ValueMap< reco::PhotonRef >  pfEGCandToPhotonMap,
edm::Handle< reco::VertexCollection > &  pvVertices,
reco::PhotonCollection &  outputCollection,
int &  iSC,
const edm::Handle< edm::ValueMap< float >> &  chargedHadrons,
const edm::Handle< edm::ValueMap< float >> &  neutralHadrons,
const edm::Handle< edm::ValueMap< float >> &  photons,
const edm::Handle< edm::ValueMap< float >> &  chargedHadronsWorstVtx,
const edm::Handle< edm::ValueMap< float >> &  chargedHadronsWorstVtxGeomVeto,
const edm::Handle< edm::ValueMap< float >> &  chargedHadronsPFPV,
const edm::Handle< edm::ValueMap< float >> &  pfEcalClusters,
const edm::Handle< edm::ValueMap< float >> &  pfHcalClusters 
)
private

Definition at line 832 of file GEDPhotonProducer.cc.

References PhotonEnergyCorrector::calculate(), candidateP4type_, reco::Photon::PflowIsolationVariables::chargedHadronIso, reco::Photon::PflowIsolationVariables::chargedHadronPFPVIso, reco::Photon::PflowIsolationVariables::chargedHadronWorstVtxGeomVetoIso, reco::Photon::PflowIsolationVariables::chargedHadronWorstVtxIso, reco::Photon::ecal_photons, EcalBarrel, EcalEndcap, EcalTools::isHGCalDet(), edm::Ref< C, T, F >::isNonnull(), GEDPhotonProducer::RecoStepInfo::isOOT(), edm::EDGetTokenT< T >::isUninitialized(), reco::Photon::PflowIsolationVariables::neutralHadronIso, reco::Photon::p4(), phoPFECALClusIsolationToken_, phoPFHCALClusIsolationToken_, reco::Photon::PflowIsolationVariables::photonIso, preselCutValuesBarrel_, preselCutValuesEndcap_, ptFast(), recoStep_, reco::Photon::regression1, reco::Photon::regression2, reco::Photon::setCandidateP4type(), reco::Photon::setP4(), reco::Photon::setPflowIDVariables(), reco::Photon::setPflowIsolationVariables(), reco::Photon::PflowIsolationVariables::sumEcalClusterEt, reco::Photon::PflowIsolationVariables::sumHcalClusterEt, and thePhotonEnergyCorrector_.

847  {
848 
849 
850 
851  std::vector<double> preselCutValues;
852 
853 
854  for(unsigned int lSC=0; lSC < photonHandle->size(); lSC++) {
855  reco::PhotonRef phoRef(reco::PhotonRef(photonHandle, lSC));
856  reco::SuperClusterRef parentSCRef = phoRef->parentSuperCluster();
857  reco::SuperClusterRef scRef=phoRef->superCluster();
858  DetId::Detector thedet = scRef->seed()->hitsAndFractions()[0].first.det();
859  int subdet = scRef->seed()->hitsAndFractions()[0].first.subdetId();
860  if (subdet==EcalBarrel) {
861  preselCutValues = preselCutValuesBarrel_;
862  } else if (subdet==EcalEndcap) {
863  preselCutValues = preselCutValuesEndcap_;
864  } else if (EcalTools::isHGCalDet(thedet)) {
865  preselCutValues = preselCutValuesEndcap_;
866  } else {
867  edm::LogWarning("")<<"GEDPhotonProducer: do not know if it is a barrel or endcap SuperCluster" << thedet << ' ' << subdet;
868  }
869 
870 
871 
872  // SC energy preselection
873  if (parentSCRef.isNonnull() &&
874  ptFast(parentSCRef->energy(),parentSCRef->position(),math::XYZPoint(0,0,0)) <= preselCutValues[0] ) continue;
875  reco::Photon newCandidate(*phoRef);
876  iSC++;
877 
878 
879  // Calculate the PF isolation and ID - for the time being there is no calculation. Only the setting
880  reco::Photon::PflowIsolationVariables pfIso;
881  reco::Photon::PflowIDVariables pfID;
882 
883  //get the pointer for the photon object
884  edm::Ptr<reco::Photon> photonPtr(photonHandle, lSC);
885 
886  if(!recoStep_.isOOT()){ //out of time photons do not have PF info so skip in this case
887  pfIso.chargedHadronIso = (*chargedHadrons)[photonPtr] ;
888  pfIso.neutralHadronIso = (*neutralHadrons)[photonPtr];
889  pfIso.photonIso = (*photons)[photonPtr];
890  pfIso.chargedHadronWorstVtxIso = (*chargedHadronsWorstVtx)[photonPtr];
891  pfIso.chargedHadronWorstVtxGeomVetoIso = (*chargedHadronsWorstVtxGeomVeto)[photonPtr];
892  pfIso.chargedHadronPFPVIso = (*chargedHadronsPFPV)[photonPtr];
893  }
894 
895  //OOT photons in legacy 80X reminiAOD workflow dont have pf cluster isolation embeded into them at this stage
896  if(!phoPFECALClusIsolationToken_.isUninitialized()){
897  pfIso.sumEcalClusterEt = (*pfEcalClusters)[photonPtr];
898  }else pfIso.sumEcalClusterEt = 0.;
899 
900  if(!phoPFHCALClusIsolationToken_.isUninitialized()){
901  pfIso.sumHcalClusterEt = (*pfHcalClusters)[photonPtr];
902  }else pfIso.sumHcalClusterEt = 0.;
903 
904  newCandidate.setPflowIsolationVariables(pfIso);
905  newCandidate.setPflowIDVariables(pfID);
906 
907  // do the regression
908  thePhotonEnergyCorrector_->calculate(evt, newCandidate, subdet, *vertexHandle, es);
909  if ( candidateP4type_ == "fromEcalEnergy") {
910  newCandidate.setP4( newCandidate.p4(reco::Photon::ecal_photons) );
911  newCandidate.setCandidateP4type(reco::Photon::ecal_photons);
912  } else if ( candidateP4type_ == "fromRegression1") {
913  newCandidate.setP4( newCandidate.p4(reco::Photon::regression1) );
914  newCandidate.setCandidateP4type(reco::Photon::regression1);
915  } else if ( candidateP4type_ == "fromRegression2") {
916  newCandidate.setP4( newCandidate.p4(reco::Photon::regression2) );
917  newCandidate.setCandidateP4type(reco::Photon::regression2);
918  } else if ( candidateP4type_ == "fromRefinedSCRegression" ) {
919  newCandidate.setP4( newCandidate.p4(reco::Photon::regression2) );
920  newCandidate.setCandidateP4type(reco::Photon::regression2);
921  }
922 
923  // std::cout << " GEDPhotonProducer pf based isolation chargedHadron " << newCandidate.chargedHadronIso() << " neutralHadron " << newCandidate.neutralHadronIso() << " Photon " << newCandidate.photonIso() << std::endl;
924  //std::cout << " GEDPhotonProducer from candidate HoE with towers in a cone " << newCandidate.hadronicOverEm() << " " << newCandidate.hadronicDepth1OverEm() << " " << newCandidate.hadronicDepth2OverEm() << std::endl;
925  //std::cout << " GEDPhotonProducer from candidate of HoE with towers behind the BCs " << newCandidate.hadTowOverEm() << " " << newCandidate.hadTowDepth1OverEm() << " " << newCandidate.hadTowDepth2OverEm() << std::endl;
926  //std::cout << " standard p4 before " << newCandidate.p4() << " energy " << newCandidate.energy() << std::endl;
927  //std::cout << " type " <<newCandidate.getCandidateP4type() << " standard p4 after " << newCandidate.p4() << " energy " << newCandidate.energy() << std::endl;
928  //std::cout << " final p4 " << newCandidate.p4() << " energy " << newCandidate.energy() << std::endl;
929 
930  outputPhotonCollection.push_back(newCandidate);
931 
932  }
933 
934 }
EcalTools::isHGCalDet
static bool isHGCalDet(DetId::Detector thedet)
identify HGCal cells
Definition: EcalTools.h:54
GEDPhotonProducer::thePhotonEnergyCorrector_
PhotonEnergyCorrector * thePhotonEnergyCorrector_
Definition: GEDPhotonProducer.h:168
edm::Ref::isNonnull
bool isNonnull() const
Checks for non-null.
Definition: Ref.h:251
GEDPhotonProducer::recoStep_
RecoStepInfo recoStep_
Definition: GEDPhotonProducer.h:150
ptFast
double ptFast(const double energy, const math::XYZPoint &position, const math::XYZPoint &origin)
Definition: SuperClusterImporter.cc:14
gedPhotons_cfi.outputPhotonCollection
outputPhotonCollection
Definition: gedPhotons_cfi.py:30
GEDPhotonProducer::phoPFHCALClusIsolationToken_
edm::EDGetTokenT< edm::ValueMap< float > > phoPFHCALClusIsolationToken_
Definition: GEDPhotonProducer.h:124
PhotonEnergyCorrector::calculate
void calculate(edm::Event &evt, reco::Photon &, int subdet, const reco::VertexCollection &vtxcol, const edm::EventSetup &iSetup)
Definition: PhotonEnergyCorrector.cc:92
GEDPhotonProducer::preselCutValuesBarrel_
std::vector< double > preselCutValuesBarrel_
Definition: GEDPhotonProducer.h:164
math::XYZPoint
XYZPointD XYZPoint
point in space with cartesian internal representation
Definition: Point3D.h:12
DetId::Detector
Detector
Definition: DetId.h:26
GEDPhotonProducer::preselCutValuesEndcap_
std::vector< double > preselCutValuesEndcap_
Definition: GEDPhotonProducer.h:165
GEDPhotonProducer::candidateP4type_
std::string candidateP4type_
Definition: GEDPhotonProducer.h:169
GEDPhotonProducer::phoPFECALClusIsolationToken_
edm::EDGetTokenT< edm::ValueMap< float > > phoPFECALClusIsolationToken_
Definition: GEDPhotonProducer.h:123
edm::EDGetTokenT::isUninitialized
bool isUninitialized() const
Definition: EDGetToken.h:70
void GEDPhotonProducer::produce ( edm::Event &  evt,
const edm::EventSetup &  es 
)
override

Definition at line 261 of file GEDPhotonProducer.cc.

References barrelEcalHits_, particleFlowSuperClusterOOTECAL_cff::barrelRecHits, endcapEcalHits_, particleFlowSuperClusterOOTECAL_cff::endcapRecHits, Exception, objects.autophobj::filler, fillPhotonCollection(), PhotonEnergyCorrector::gedRegression(), edm::EventSetup::get(), edm::Event::getByToken(), hcalTowers_, PhotonEnergyCorrector::init(), GEDPhotonProducer::RecoStepInfo::isFinal(), GEDPhotonProducer::RecoStepInfo::isOOT(), edm::EDGetTokenT< T >::isUninitialized(), edm::HandleBase::isValid(), edm::InputTag::label(), eostools::move(), gedPhotons_cfi::outputPhotonCollection, pfCandidates_, pfEgammaCandidates_, phoChargedIsolationToken_, phoChargedPFPVIsoToken_, phoChargedWorstVtxGeomVetoIsoToken_, phoChargedWorstVtxIsoToken_, phoNeutralHadronIsolationToken_, phoPFECALClusIsolationToken_, phoPFHCALClusIsolationToken_, phoPhotonIsolationToken_, photonCollection_, photonCoreProducerT_, photonProducer_, photonProducerT_, preshowerHits_, edm::Handle< T >::product(), edm::ESHandle< T >::product(), edm::Event::put(), recoStep_, theCaloGeom_, theCaloTopo_, thePhotonEnergyCorrector_, ecaldqm::topology(), usePrimaryVertex_, valueMapPFCandPhoton_, MuonErrorMatrixValues_cff::values, particleFlowSuperClusterECAL_cfi::vertexCollection, and vertexProducer_.

261  {
262 
263  using namespace edm;
264  // nEvt_++;
265 
266  reco::PhotonCollection outputPhotonCollection;
267  auto outputPhotonCollection_p = std::make_unique<reco::PhotonCollection>();
268  edm::ValueMap<reco::PhotonRef> pfEGCandToPhotonMap;
269 
270 
271  // Get the PhotonCore collection
272  bool validPhotonCoreHandle=false;
273  Handle<reco::PhotonCoreCollection> photonCoreHandle;
274  bool validPhotonHandle= false;
275  Handle<reco::PhotonCollection> photonHandle;
276  //value maps for isolation
277  edm::Handle<edm::ValueMap<float> > phoChargedIsolationMap;
278  edm::Handle<edm::ValueMap<float> > phoNeutralHadronIsolationMap;
279  edm::Handle<edm::ValueMap<float> > phoPhotonIsolationMap;
280  edm::Handle<edm::ValueMap<float> > phoChargedWorstVtxIsoMap;
281  edm::Handle<edm::ValueMap<float> > phoChargedWorstVtxGeomVetoIsoMap;
282  edm::Handle<edm::ValueMap<float> > phoChargedPFPVIsoMap;
283 
284  edm::Handle<edm::ValueMap<float> > phoPFECALClusIsolationMap;
285  edm::Handle<edm::ValueMap<float> > phoPFHCALClusIsolationMap;
286 
287  if ( recoStep_.isFinal() ) {
288  theEvent.getByToken(photonProducerT_,photonHandle);
289  //get isolation objects
290  theEvent.getByToken(phoChargedIsolationToken_,phoChargedIsolationMap);
291  theEvent.getByToken(phoNeutralHadronIsolationToken_,phoNeutralHadronIsolationMap);
292  theEvent.getByToken(phoPhotonIsolationToken_,phoPhotonIsolationMap);
293  theEvent.getByToken(phoChargedWorstVtxIsoToken_,phoChargedWorstVtxIsoMap);
294  theEvent.getByToken(phoChargedWorstVtxGeomVetoIsoToken_,phoChargedWorstVtxGeomVetoIsoMap);
295  theEvent.getByToken(phoChargedPFPVIsoToken_,phoChargedPFPVIsoMap);
296 
297  //OOT photons in legacy 80X re-miniAOD workflow dont have cluster isolation embed in them
298  if(!phoPFECALClusIsolationToken_.isUninitialized()) {
299  theEvent.getByToken(phoPFECALClusIsolationToken_,phoPFECALClusIsolationMap);
300  }
301  if(!phoPFHCALClusIsolationToken_.isUninitialized()){
302  theEvent.getByToken(phoPFHCALClusIsolationToken_,phoPFHCALClusIsolationMap);
303  }
304 
305  if ( photonHandle.isValid()) {
306  validPhotonHandle=true;
307  } else {
308  throw cms::Exception("GEDPhotonProducer") << "Error! Can't get the product " << photonProducer_.label() << "\n";
309  }
310  } else {
311 
312  theEvent.getByToken(photonCoreProducerT_,photonCoreHandle);
313  if (photonCoreHandle.isValid()) {
314  validPhotonCoreHandle=true;
315  } else {
316  throw cms::Exception("GEDPhotonProducer")
317  << "Error! Can't get the photonCoreProducer " << photonProducer_.label() << "\n";
318  }
319  }
320 
321  // Get EcalRecHits
322  bool validEcalRecHits=true;
323  Handle<EcalRecHitCollection> barrelHitHandle;
324  const EcalRecHitCollection dummyEB;
325  theEvent.getByToken(barrelEcalHits_, barrelHitHandle);
326  if (!barrelHitHandle.isValid()) {
327  throw cms::Exception("GEDPhotonProducer")
328  << "Error! Can't get the barrelEcalHits";
329  }
330  const EcalRecHitCollection& barrelRecHits(validEcalRecHits ? *(barrelHitHandle.product()) : dummyEB);
331 
332  Handle<EcalRecHitCollection> endcapHitHandle;
333  theEvent.getByToken(endcapEcalHits_, endcapHitHandle);
334  const EcalRecHitCollection dummyEE;
335  if (!endcapHitHandle.isValid()) {
336  throw cms::Exception("GEDPhotonProducer")
337  << "Error! Can't get the endcapEcalHits";
338  }
339  const EcalRecHitCollection& endcapRecHits(validEcalRecHits ? *(endcapHitHandle.product()) : dummyEE);
340 
341  bool validPreshowerRecHits=true;
342  Handle<EcalRecHitCollection> preshowerHitHandle;
343  theEvent.getByToken(preshowerHits_, preshowerHitHandle);
344  EcalRecHitCollection preshowerRecHits;
345  if (!preshowerHitHandle.isValid()) {
346  throw cms::Exception("GEDPhotonProducer")
347  << "Error! Can't get the preshowerEcalHits";
348  }
349  if( validPreshowerRecHits ) preshowerRecHits = *(preshowerHitHandle.product());
350 
351 
352 
353  Handle<reco::PFCandidateCollection> pfEGCandidateHandle;
354  // Get the PF refined cluster collection
355  if (not pfEgammaCandidates_.isUninitialized()){
356  theEvent.getByToken(pfEgammaCandidates_,pfEGCandidateHandle);
357  if (!pfEGCandidateHandle.isValid()) {
358  throw cms::Exception("GEDPhotonProducer")
359  << "Error! Can't get the pfEgammaCandidates";
360  }
361  }
362 
363  Handle<reco::PFCandidateCollection> pfCandidateHandle;
364 
365  if ( recoStep_.isFinal() ) {
366  // Get the PF candidates collection
367  theEvent.getByToken(pfCandidates_,pfCandidateHandle);
368  //OOT photons have no PF candidates so its not an error in this case
369  if (!pfCandidateHandle.isValid() && !recoStep_.isOOT()) {
370  throw cms::Exception("GEDPhotonProducer")
371  << "Error! Can't get the pfCandidates";
372  }
373  }
374 
375  //AA
376  //Get the severity level object
377  edm::ESHandle<EcalSeverityLevelAlgo> sevLv;
378  theEventSetup.get<EcalSeverityLevelAlgoRcd>().get(sevLv);
379  //
380 
381 
382 // get Hcal towers collection
383  Handle<CaloTowerCollection> hcalTowersHandle;
384  if (not hcalTowers_.isUninitialized()){
385  theEvent.getByToken(hcalTowers_, hcalTowersHandle);
386  }
387 
388  // get the geometry from the event setup:
389  theEventSetup.get<CaloGeometryRecord>().get(theCaloGeom_);
390 
391  //
392  // update energy correction function
393  // energyCorrectionF->init(theEventSetup);
394 
395  edm::ESHandle<CaloTopology> pTopology;
396  theEventSetup.get<CaloTopologyRecord>().get(theCaloTopo_);
397  const CaloTopology *topology = theCaloTopo_.product();
398 
399  // Get the primary event vertex
400  Handle<reco::VertexCollection> vertexHandle;
401  const reco::VertexCollection dummyVC;
402  bool validVertex=true;
403  if ( usePrimaryVertex_ ) {
404  theEvent.getByToken(vertexProducer_, vertexHandle);
405  if (!vertexHandle.isValid()) {
406  throw cms::Exception("GEDPhotonProducer")
407  << "Error! Can't get the product primary Vertex Collection";
408  }
409  }
410  const reco::VertexCollection& vertexCollection(usePrimaryVertex_ && validVertex ? *(vertexHandle.product()) : dummyVC);
411 
412  // math::XYZPoint vtx(0.,0.,0.);
413  //if (vertexCollection.size()>0) vtx = vertexCollection.begin()->position();
414 
415  // get the regression calculator ready
416  thePhotonEnergyCorrector_->init(theEventSetup);
417  if( thePhotonEnergyCorrector_->gedRegression() ) {
418  thePhotonEnergyCorrector_->gedRegression()->setEvent(theEvent);
419  thePhotonEnergyCorrector_->gedRegression()->setEventContent(theEventSetup);
420  }
421 
422 
423  int iSC=0; // index in photon collection
424  // Loop over barrel and endcap SC collections and fill the photon collection
425  if ( validPhotonCoreHandle)
426  fillPhotonCollection(theEvent,
427  theEventSetup,
428  photonCoreHandle,
429  topology,
430  &barrelRecHits,
431  &endcapRecHits,
432  &preshowerRecHits,
433  hcalTowersHandle,
434  //vtx,
435  vertexCollection,
436  outputPhotonCollection,
437  iSC);
438 
439  iSC=0;
440  if ( validPhotonHandle && recoStep_.isFinal() )
441  fillPhotonCollection(theEvent,
442  theEventSetup,
443  photonHandle,
444  pfCandidateHandle,
445  pfEGCandidateHandle,
446  pfEGCandToPhotonMap,
447  vertexHandle,
448  outputPhotonCollection,
449  iSC,
450  phoChargedIsolationMap,
451  phoNeutralHadronIsolationMap,
452  phoPhotonIsolationMap,
453  phoChargedWorstVtxIsoMap,
454  phoChargedWorstVtxGeomVetoIsoMap,
455  phoChargedPFPVIsoMap,
456  phoPFECALClusIsolationMap,
457  phoPFHCALClusIsolationMap);
458 
459 
460 
461  // put the product in the event
462  edm::LogInfo("GEDPhotonProducer") << " Put in the event " << iSC << " Photon Candidates \n";
463  outputPhotonCollection_p->assign(outputPhotonCollection.begin(),outputPhotonCollection.end());
464  const edm::OrphanHandle<reco::PhotonCollection> photonOrphHandle = theEvent.put(std::move(outputPhotonCollection_p), photonCollection_);
465 
466 
467  if ( !recoStep_.isFinal() && not pfEgammaCandidates_.isUninitialized()) {
469  auto pfEGCandToPhotonMap_p = std::make_unique<edm::ValueMap<reco::PhotonRef>>();
470  edm::ValueMap<reco::PhotonRef>::Filler filler(*pfEGCandToPhotonMap_p);
471  unsigned nObj = pfEGCandidateHandle->size();
472  std::vector<reco::PhotonRef> values(nObj);
474  for(unsigned int lCand=0; lCand < nObj; lCand++) {
475  reco::PFCandidateRef pfCandRef (reco::PFCandidateRef(pfEGCandidateHandle,lCand));
476  reco::SuperClusterRef pfScRef = pfCandRef -> superClusterRef();
477 
478  for(unsigned int lSC=0; lSC < photonOrphHandle->size(); lSC++) {
479  reco::PhotonRef photonRef(reco::PhotonRef(photonOrphHandle, lSC));
480  reco::SuperClusterRef scRef=photonRef->superCluster();
481  if ( pfScRef != scRef ) continue;
482  values[lCand] = photonRef;
483  }
484  }
485 
486 
487  filler.insert(pfEGCandidateHandle,values.begin(),values.end());
488  filler.fill();
489  theEvent.put(std::move(pfEGCandToPhotonMap_p),valueMapPFCandPhoton_);
490 
491 
492  }
493 
494 
495 
496 
497 
498 
499 }
GEDPhotonProducer::photonProducer_
edm::InputTag photonProducer_
Definition: GEDPhotonProducer.h:102
GEDPhotonProducer::thePhotonEnergyCorrector_
PhotonEnergyCorrector * thePhotonEnergyCorrector_
Definition: GEDPhotonProducer.h:168
GEDPhotonProducer::pfCandidates_
edm::EDGetTokenT< reco::PFCandidateCollection > pfCandidates_
Definition: GEDPhotonProducer.h:110
GEDPhotonProducer::theCaloGeom_
edm::ESHandle< CaloGeometry > theCaloGeom_
Definition: GEDPhotonProducer.h:156
GEDPhotonProducer::hcalTowers_
edm::EDGetTokenT< CaloTowerCollection > hcalTowers_
Definition: GEDPhotonProducer.h:111
ecaldqm::topology
CaloTopology const * topology(0)
GEDPhotonProducer::recoStep_
RecoStepInfo recoStep_
Definition: GEDPhotonProducer.h:150
reco::VertexCollection
std::vector< Vertex > VertexCollection
collection of Vertex objects
Definition: VertexFwd.h:9
GEDPhotonProducer::endcapEcalHits_
edm::EDGetTokenT< EcalRecHitCollection > endcapEcalHits_
Definition: GEDPhotonProducer.h:107
PhotonEnergyCorrector::gedRegression
std::unique_ptr< ModifyObjectValueBase > & gedRegression()
Definition: PhotonEnergyCorrector.h:33
GEDPhotonProducer::photonCollection_
std::string photonCollection_
Definition: GEDPhotonProducer.h:101
gedPhotons_cfi.outputPhotonCollection
outputPhotonCollection
Definition: gedPhotons_cfi.py:30
GEDPhotonProducer::theCaloTopo_
edm::ESHandle< CaloTopology > theCaloTopo_
Definition: GEDPhotonProducer.h:157
GEDPhotonProducer::phoChargedPFPVIsoToken_
edm::EDGetTokenT< edm::ValueMap< float > > phoChargedPFPVIsoToken_
Definition: GEDPhotonProducer.h:121
GEDPhotonProducer::barrelEcalHits_
edm::EDGetTokenT< EcalRecHitCollection > barrelEcalHits_
Definition: GEDPhotonProducer.h:106
GEDPhotonProducer::phoChargedWorstVtxGeomVetoIsoToken_
edm::EDGetTokenT< edm::ValueMap< float > > phoChargedWorstVtxGeomVetoIsoToken_
Definition: GEDPhotonProducer.h:120
GEDPhotonProducer::phoChargedIsolationToken_
edm::EDGetTokenT< edm::ValueMap< float > > phoChargedIsolationToken_
Definition: GEDPhotonProducer.h:116
GEDPhotonProducer::photonCoreProducerT_
edm::EDGetTokenT< reco::PhotonCoreCollection > photonCoreProducerT_
Definition: GEDPhotonProducer.h:104
GEDPhotonProducer::phoPFHCALClusIsolationToken_
edm::EDGetTokenT< edm::ValueMap< float > > phoPFHCALClusIsolationToken_
Definition: GEDPhotonProducer.h:124
PhotonEnergyCorrector::init
void init(const edm::EventSetup &theEventSetup)
Definition: PhotonEnergyCorrector.cc:70
edm::HandleBase::isValid
bool isValid() const
Definition: HandleBase.h:74
GEDPhotonProducer::pfEgammaCandidates_
edm::EDGetTokenT< reco::PFCandidateCollection > pfEgammaCandidates_
Definition: GEDPhotonProducer.h:109
GEDPhotonProducer::phoNeutralHadronIsolationToken_
edm::EDGetTokenT< edm::ValueMap< float > > phoNeutralHadronIsolationToken_
Definition: GEDPhotonProducer.h:117
GEDPhotonProducer::phoChargedWorstVtxIsoToken_
edm::EDGetTokenT< edm::ValueMap< float > > phoChargedWorstVtxIsoToken_
Definition: GEDPhotonProducer.h:119
edm::Handle::product
T const * product() const
Definition: Handle.h:74
GEDPhotonProducer::photonProducerT_
edm::EDGetTokenT< reco::PhotonCollection > photonProducerT_
Definition: GEDPhotonProducer.h:105
particleFlowSuperClusterOOTECAL_cff.barrelRecHits
barrelRecHits
Definition: particleFlowSuperClusterOOTECAL_cff.py:16
reco::PhotonCollection
std::vector< Photon > PhotonCollection
collectin of Photon objects
Definition: PhotonFwd.h:9
GEDPhotonProducer::fillPhotonCollection
void fillPhotonCollection(edm::Event &evt, edm::EventSetup const &es, const edm::Handle< reco::PhotonCoreCollection > &photonCoreHandle, const CaloTopology *topology, const EcalRecHitCollection *ecalBarrelHits, const EcalRecHitCollection *ecalEndcapHits, const EcalRecHitCollection *preshowerHits, const edm::Handle< CaloTowerCollection > &hcalTowersHandle, const reco::VertexCollection &pvVertices, reco::PhotonCollection &outputCollection, int &iSC)
Definition: GEDPhotonProducer.cc:501
GEDPhotonProducer::phoPhotonIsolationToken_
edm::EDGetTokenT< edm::ValueMap< float > > phoPhotonIsolationToken_
Definition: GEDPhotonProducer.h:118
particleFlowSuperClusterECAL_cfi.vertexCollection
vertexCollection
Definition: particleFlowSuperClusterECAL_cfi.py:41
edm::InputTag::label
std::string const & label() const
Definition: InputTag.h:36
MuonErrorMatrixValues_cff.values
values
Definition: MuonErrorMatrixValues_cff.py:6
edm
HLT enums.
Definition: AlignableModifier.h:17
particleFlowSuperClusterOOTECAL_cff.endcapRecHits
endcapRecHits
Definition: particleFlowSuperClusterOOTECAL_cff.py:17
GEDPhotonProducer::phoPFECALClusIsolationToken_
edm::EDGetTokenT< edm::ValueMap< float > > phoPFECALClusIsolationToken_
Definition: GEDPhotonProducer.h:123
edm::EDGetTokenT::isUninitialized
bool isUninitialized() const
Definition: EDGetToken.h:70
GEDPhotonProducer::vertexProducer_
edm::EDGetTokenT< reco::VertexCollection > vertexProducer_
Definition: GEDPhotonProducer.h:112
edm::ESHandle::product
T const * product() const
Definition: ESHandle.h:86
GEDPhotonProducer::preshowerHits_
edm::EDGetTokenT< EcalRecHitCollection > preshowerHits_
Definition: GEDPhotonProducer.h:108
eostools.move
def move(src, dest)
Definition: eostools.py:511
GEDPhotonProducer::valueMapPFCandPhoton_
std::string valueMapPFCandPhoton_
Definition: GEDPhotonProducer.h:128

Member Data Documentation

edm::EDGetTokenT<EcalRecHitCollection> GEDPhotonProducer::barrelEcalHits_
private

Definition at line 106 of file GEDPhotonProducer.h.

Referenced by fillPhotonCollection(), GEDPhotonProducer(), and produce().

std::string GEDPhotonProducer::candidateP4type_
private

Definition at line 169 of file GEDPhotonProducer.h.

Referenced by fillPhotonCollection(), and GEDPhotonProducer().

bool GEDPhotonProducer::checkHcalStatus_
private

Definition at line 171 of file GEDPhotonProducer.h.

Referenced by fillPhotonCollection(), and GEDPhotonProducer().

edm::ParameterSet GEDPhotonProducer::conf_
private

Definition at line 153 of file GEDPhotonProducer.h.

Referenced by GEDPhotonProducer().

std::string GEDPhotonProducer::conversionCollection_
private

Definition at line 127 of file GEDPhotonProducer.h.

std::string GEDPhotonProducer::conversionProducer_
private

Definition at line 126 of file GEDPhotonProducer.h.

edm::EDGetTokenT<EcalRecHitCollection> GEDPhotonProducer::endcapEcalHits_
private

Definition at line 107 of file GEDPhotonProducer.h.

Referenced by fillPhotonCollection(), GEDPhotonProducer(), and produce().

EcalClusterFunctionBaseClass* GEDPhotonProducer::energyCorrectionF
private

Definition at line 167 of file GEDPhotonProducer.h.

std::vector<int> GEDPhotonProducer::flagsexclEB_
private

Definition at line 135 of file GEDPhotonProducer.h.

Referenced by fillPhotonCollection(), and GEDPhotonProducer().

std::vector<int> GEDPhotonProducer::flagsexclEE_
private

Definition at line 136 of file GEDPhotonProducer.h.

Referenced by fillPhotonCollection(), and GEDPhotonProducer().

edm::EDGetTokenT<CaloTowerCollection> GEDPhotonProducer::hcalTowers_
private

Definition at line 111 of file GEDPhotonProducer.h.

Referenced by fillPhotonCollection(), GEDPhotonProducer(), and produce().

double GEDPhotonProducer::highEt_
private

Definition at line 144 of file GEDPhotonProducer.h.

Referenced by fillPhotonCollection(), and GEDPhotonProducer().

double GEDPhotonProducer::hOverEConeSize_
private

Definition at line 141 of file GEDPhotonProducer.h.

Referenced by fillPhotonCollection(), and GEDPhotonProducer().

double GEDPhotonProducer::maxHOverE_
private

Definition at line 142 of file GEDPhotonProducer.h.

double GEDPhotonProducer::minR9Barrel_
private

Definition at line 145 of file GEDPhotonProducer.h.

Referenced by GEDPhotonProducer().

double GEDPhotonProducer::minR9Endcap_
private

Definition at line 146 of file GEDPhotonProducer.h.

Referenced by GEDPhotonProducer().

double GEDPhotonProducer::minSCEt_
private

Definition at line 143 of file GEDPhotonProducer.h.

edm::EDGetTokenT<edm::ValueMap<std::vector<reco::PFCandidateRef > > > GEDPhotonProducer::particleBasedIsolationToken
private

Definition at line 114 of file GEDPhotonProducer.h.

edm::EDGetTokenT<reco::PFCandidateCollection> GEDPhotonProducer::pfCandidates_
private

Definition at line 110 of file GEDPhotonProducer.h.

Referenced by GEDPhotonProducer(), and produce().

edm::EDGetTokenT<reco::PFCandidateCollection> GEDPhotonProducer::pfEgammaCandidates_
private

Definition at line 109 of file GEDPhotonProducer.h.

Referenced by GEDPhotonProducer(), and produce().

edm::EDGetTokenT<edm::ValueMap<float> > GEDPhotonProducer::phoChargedIsolationToken_
private

Definition at line 116 of file GEDPhotonProducer.h.

Referenced by GEDPhotonProducer(), and produce().

edm::EDGetTokenT<edm::ValueMap<float> > GEDPhotonProducer::phoChargedPFPVIsoToken_
private

Definition at line 121 of file GEDPhotonProducer.h.

Referenced by GEDPhotonProducer(), and produce().

edm::EDGetTokenT<edm::ValueMap<float> > GEDPhotonProducer::phoChargedWorstVtxGeomVetoIsoToken_
private

Definition at line 120 of file GEDPhotonProducer.h.

Referenced by GEDPhotonProducer(), and produce().

edm::EDGetTokenT<edm::ValueMap<float> > GEDPhotonProducer::phoChargedWorstVtxIsoToken_
private

Definition at line 119 of file GEDPhotonProducer.h.

Referenced by GEDPhotonProducer(), and produce().

edm::EDGetTokenT<edm::ValueMap<float> > GEDPhotonProducer::phoNeutralHadronIsolationToken_
private

Definition at line 117 of file GEDPhotonProducer.h.

Referenced by GEDPhotonProducer(), and produce().

edm::EDGetTokenT<edm::ValueMap<float> > GEDPhotonProducer::phoPFECALClusIsolationToken_
private

Definition at line 123 of file GEDPhotonProducer.h.

Referenced by fillPhotonCollection(), GEDPhotonProducer(), and produce().

edm::EDGetTokenT<edm::ValueMap<float> > GEDPhotonProducer::phoPFHCALClusIsolationToken_
private

Definition at line 124 of file GEDPhotonProducer.h.

Referenced by fillPhotonCollection(), GEDPhotonProducer(), and produce().

edm::EDGetTokenT<edm::ValueMap<float> > GEDPhotonProducer::phoPhotonIsolationToken_
private

Definition at line 118 of file GEDPhotonProducer.h.

Referenced by GEDPhotonProducer(), and produce().

std::string GEDPhotonProducer::photonCollection_
private

Definition at line 101 of file GEDPhotonProducer.h.

Referenced by GEDPhotonProducer(), and produce().

edm::EDGetTokenT<reco::PhotonCoreCollection> GEDPhotonProducer::photonCoreProducerT_
private

Definition at line 104 of file GEDPhotonProducer.h.

Referenced by GEDPhotonProducer(), and produce().

edm::InputTag GEDPhotonProducer::photonProducer_
private

Definition at line 102 of file GEDPhotonProducer.h.

Referenced by GEDPhotonProducer(), and produce().

edm::EDGetTokenT<reco::PhotonCollection> GEDPhotonProducer::photonProducerT_
private

Definition at line 105 of file GEDPhotonProducer.h.

Referenced by GEDPhotonProducer(), and produce().

PositionCalc GEDPhotonProducer::posCalculator_
private

Definition at line 154 of file GEDPhotonProducer.h.

Referenced by GEDPhotonProducer().

std::vector<double> GEDPhotonProducer::preselCutValuesBarrel_
private

Definition at line 164 of file GEDPhotonProducer.h.

Referenced by fillPhotonCollection(), and GEDPhotonProducer().

std::vector<double> GEDPhotonProducer::preselCutValuesEndcap_
private

Definition at line 165 of file GEDPhotonProducer.h.

Referenced by fillPhotonCollection(), and GEDPhotonProducer().

edm::EDGetTokenT<EcalRecHitCollection> GEDPhotonProducer::preshowerHits_
private

Definition at line 108 of file GEDPhotonProducer.h.

Referenced by fillPhotonCollection(), GEDPhotonProducer(), and produce().

RecoStepInfo GEDPhotonProducer::recoStep_
private

Definition at line 150 of file GEDPhotonProducer.h.

Referenced by beginRun(), fillPhotonCollection(), GEDPhotonProducer(), and produce().

bool GEDPhotonProducer::runMIPTagger_
private

Definition at line 147 of file GEDPhotonProducer.h.

Referenced by fillPhotonCollection(), and GEDPhotonProducer().

std::vector<int> GEDPhotonProducer::severitiesexclEB_
private

Definition at line 137 of file GEDPhotonProducer.h.

Referenced by fillPhotonCollection(), and GEDPhotonProducer().

std::vector<int> GEDPhotonProducer::severitiesexclEE_
private

Definition at line 138 of file GEDPhotonProducer.h.

Referenced by fillPhotonCollection(), and GEDPhotonProducer().

edm::ESHandle<CaloGeometry> GEDPhotonProducer::theCaloGeom_
private

Definition at line 156 of file GEDPhotonProducer.h.

Referenced by fillPhotonCollection(), and produce().

edm::ESHandle<CaloTopology> GEDPhotonProducer::theCaloTopo_
private

Definition at line 157 of file GEDPhotonProducer.h.

Referenced by produce().

PhotonEnergyCorrector* GEDPhotonProducer::thePhotonEnergyCorrector_
private

Definition at line 168 of file GEDPhotonProducer.h.

Referenced by beginRun(), fillPhotonCollection(), GEDPhotonProducer(), produce(), and ~GEDPhotonProducer().

PhotonIsolationCalculator* GEDPhotonProducer::thePhotonIsolationCalculator_
private

Definition at line 130 of file GEDPhotonProducer.h.

Referenced by fillPhotonCollection(), GEDPhotonProducer(), and ~GEDPhotonProducer().

PhotonMIPHaloTagger* GEDPhotonProducer::thePhotonMIPHaloTagger_
private

Definition at line 162 of file GEDPhotonProducer.h.

Referenced by fillPhotonCollection(), GEDPhotonProducer(), and ~GEDPhotonProducer().

bool GEDPhotonProducer::usePrimaryVertex_
private

Definition at line 152 of file GEDPhotonProducer.h.

Referenced by GEDPhotonProducer(), and produce().

bool GEDPhotonProducer::validConversions_
private

Definition at line 149 of file GEDPhotonProducer.h.

bool GEDPhotonProducer::validPixelSeeds_
private

Definition at line 159 of file GEDPhotonProducer.h.

std::string GEDPhotonProducer::valueMapPFCandPhoton_
private

Definition at line 128 of file GEDPhotonProducer.h.

Referenced by GEDPhotonProducer(), and produce().

edm::EDGetTokenT<reco::VertexCollection> GEDPhotonProducer::vertexProducer_
private

Definition at line 112 of file GEDPhotonProducer.h.

Referenced by GEDPhotonProducer(), and produce().

Generated for CMSSW Reference Manual by doxygen 1.8.11