CMS 3D CMS Logo
 All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Properties Friends Macros Pages
PhotonProducer.cc
Go to the documentation of this file.
1 #include <iostream>
2 #include <vector>
3 #include <memory>
4 
5 // Framework
6 #include "FWCore/MessageLogger/interface/MessageLogger.h"
7 #include "FWCore/Utilities/interface/Exception.h"
8 
9 #include "CommonTools/Utils/interface/StringToEnumValue.h"
10 
11 #include "Geometry/Records/interface/CaloGeometryRecord.h"
12 #include "Geometry/CaloGeometry/interface/CaloSubdetectorGeometry.h"
13 #include "Geometry/CaloTopology/interface/CaloTopology.h"
14 #include "Geometry/CaloEventSetup/interface/CaloTopologyRecord.h"
15 
16 
17 #include "DataFormats/VertexReco/interface/Vertex.h"
18 #include "DataFormats/VertexReco/interface/VertexFwd.h"
19 #include "DataFormats/EgammaReco/interface/ClusterShape.h"
20 #include "DataFormats/EgammaCandidates/interface/PhotonCore.h"
21 #include "DataFormats/EgammaCandidates/interface/Photon.h"
22 #include "DataFormats/EgammaCandidates/interface/PhotonFwd.h"
23 #include "DataFormats/EgammaCandidates/interface/Conversion.h"
24 
25 #include "DataFormats/EgammaReco/interface/ElectronSeed.h"
26 #include "RecoCaloTools/Selectors/interface/CaloConeSelector.h"
27 
28 #include "RecoLocalCalo/EcalRecAlgos/interface/EcalSeverityLevelAlgoRcd.h"
29 #include "RecoLocalCalo/EcalRecAlgos/interface/EcalSeverityLevelAlgo.h"
30 
31 #include "RecoEgamma/EgammaPhotonProducers/interface/PhotonProducer.h"
32 #include "RecoEgamma/EgammaIsolationAlgos/interface/EgammaTowerIsolation.h"
33 #include "RecoEcal/EgammaCoreTools/interface/EcalClusterFunctionBaseClass.h"
34 #include "RecoEcal/EgammaCoreTools/interface/EcalClusterFunctionFactory.h"
35 #include "RecoEcal/EgammaCoreTools/plugins/EcalClusterCrackCorrection.h"
36 #include "RecoEgamma/EgammaIsolationAlgos/interface/EgammaHadTower.h"
37 
38 PhotonProducer::PhotonProducer(const edm::ParameterSet& config) :
39 
40  conf_(config)
41 {
42 
43  // use onfiguration file to setup input/output collection names
44 
45  photonCoreProducer_ =
46  consumes<reco::PhotonCoreCollection>(conf_.getParameter<edm::InputTag>("photonCoreProducer"));
47  barrelEcalHits_ =
48  consumes<EcalRecHitCollection>(conf_.getParameter<edm::InputTag>("barrelEcalHits"));
49  endcapEcalHits_ =
50  consumes<EcalRecHitCollection>(conf_.getParameter<edm::InputTag>("endcapEcalHits"));
51  vertexProducer_ =
52  consumes<reco::VertexCollection>(conf_.getParameter<edm::InputTag>("primaryVertexProducer"));
53  hcalTowers_ =
54  consumes<CaloTowerCollection>(conf_.getParameter<edm::InputTag>("hcalTowers"));
55  hOverEConeSize_ = conf_.getParameter<double>("hOverEConeSize");
56  highEt_ = conf_.getParameter<double>("highEt");
57  // R9 value to decide converted/unconverted
58  minR9Barrel_ = conf_.getParameter<double>("minR9Barrel");
59  minR9Endcap_ = conf_.getParameter<double>("minR9Endcap");
60  usePrimaryVertex_ = conf_.getParameter<bool>("usePrimaryVertex");
61  runMIPTagger_ = conf_.getParameter<bool>("runMIPTagger");
62 
63  candidateP4type_ = config.getParameter<std::string>("candidateP4type") ;
64 
65  edm::ParameterSet posCalcParameters =
66  config.getParameter<edm::ParameterSet>("posCalcParameters");
67  posCalculator_ = PositionCalc(posCalcParameters);
68 
69 
70  //AA
71  //Flags and Severities to be excluded from photon calculations
72  const std::vector<std::string> flagnamesEB =
73  config.getParameter<std::vector<std::string> >("RecHitFlagToBeExcludedEB");
74 
75  const std::vector<std::string> flagnamesEE =
76  config.getParameter<std::vector<std::string> >("RecHitFlagToBeExcludedEE");
77 
78  flagsexclEB_=
79  StringToEnumValue<EcalRecHit::Flags>(flagnamesEB);
80 
81  flagsexclEE_=
82  StringToEnumValue<EcalRecHit::Flags>(flagnamesEE);
83 
84  const std::vector<std::string> severitynamesEB =
85  config.getParameter<std::vector<std::string> >("RecHitSeverityToBeExcludedEB");
86 
87  severitiesexclEB_=
88  StringToEnumValue<EcalSeverityLevel::SeverityLevel>(severitynamesEB);
89 
90  const std::vector<std::string> severitynamesEE =
91  config.getParameter<std::vector<std::string> >("RecHitSeverityToBeExcludedEE");
92 
93  severitiesexclEE_=
94  StringToEnumValue<EcalSeverityLevel::SeverityLevel>(severitynamesEE);
95 
96  //AA
97 
98  //
99 
100  // Parameters for the position calculation:
101  // std::map<std::string,double> providedParameters;
102  // providedParameters.insert(std::make_pair("LogWeighted",conf_.getParameter<bool>("posCalc_logweight")));
103  //providedParameters.insert(std::make_pair("T0_barl",conf_.getParameter<double>("posCalc_t0_barl")));
104  //providedParameters.insert(std::make_pair("T0_endc",conf_.getParameter<double>("posCalc_t0_endc")));
105  //providedParameters.insert(std::make_pair("T0_endcPresh",conf_.getParameter<double>("posCalc_t0_endcPresh")));
106  //providedParameters.insert(std::make_pair("W0",conf_.getParameter<double>("posCalc_w0")));
107  //providedParameters.insert(std::make_pair("X0",conf_.getParameter<double>("posCalc_x0")));
108  //posCalculator_ = PositionCalc(providedParameters);
109  // cut values for pre-selection
110  preselCutValuesBarrel_.push_back(conf_.getParameter<double>("minSCEtBarrel"));
111  preselCutValuesBarrel_.push_back(conf_.getParameter<double>("maxHoverEBarrel"));
112  preselCutValuesBarrel_.push_back(conf_.getParameter<double>("ecalRecHitSumEtOffsetBarrel"));
113  preselCutValuesBarrel_.push_back(conf_.getParameter<double>("ecalRecHitSumEtSlopeBarrel"));
114  preselCutValuesBarrel_.push_back(conf_.getParameter<double>("hcalTowerSumEtOffsetBarrel"));
115  preselCutValuesBarrel_.push_back(conf_.getParameter<double>("hcalTowerSumEtSlopeBarrel"));
116  preselCutValuesBarrel_.push_back(conf_.getParameter<double>("nTrackSolidConeBarrel"));
117  preselCutValuesBarrel_.push_back(conf_.getParameter<double>("nTrackHollowConeBarrel"));
118  preselCutValuesBarrel_.push_back(conf_.getParameter<double>("trackPtSumSolidConeBarrel"));
119  preselCutValuesBarrel_.push_back(conf_.getParameter<double>("trackPtSumHollowConeBarrel"));
120  preselCutValuesBarrel_.push_back(conf_.getParameter<double>("sigmaIetaIetaCutBarrel"));
121  //
122  preselCutValuesEndcap_.push_back(conf_.getParameter<double>("minSCEtEndcap"));
123  preselCutValuesEndcap_.push_back(conf_.getParameter<double>("maxHoverEEndcap"));
124  preselCutValuesEndcap_.push_back(conf_.getParameter<double>("ecalRecHitSumEtOffsetEndcap"));
125  preselCutValuesEndcap_.push_back(conf_.getParameter<double>("ecalRecHitSumEtSlopeEndcap"));
126  preselCutValuesEndcap_.push_back(conf_.getParameter<double>("hcalTowerSumEtOffsetEndcap"));
127  preselCutValuesEndcap_.push_back(conf_.getParameter<double>("hcalTowerSumEtSlopeEndcap"));
128  preselCutValuesEndcap_.push_back(conf_.getParameter<double>("nTrackSolidConeEndcap"));
129  preselCutValuesEndcap_.push_back(conf_.getParameter<double>("nTrackHollowConeEndcap"));
130  preselCutValuesEndcap_.push_back(conf_.getParameter<double>("trackPtSumSolidConeEndcap"));
131  preselCutValuesEndcap_.push_back(conf_.getParameter<double>("trackPtSumHollowConeEndcap"));
132  preselCutValuesEndcap_.push_back(conf_.getParameter<double>("sigmaIetaIetaCutEndcap"));
133  //
134 
135  thePhotonEnergyCorrector_ = new PhotonEnergyCorrector(conf_, consumesCollector());
136  thePhotonIsolationCalculator_ = new PhotonIsolationCalculator();
137  edm::ParameterSet isolationSumsCalculatorSet = conf_.getParameter<edm::ParameterSet>("isolationSumsCalculatorSet");
138  thePhotonIsolationCalculator_->setup(isolationSumsCalculatorSet, flagsexclEB_, flagsexclEE_, severitiesexclEB_, severitiesexclEE_,consumesCollector());
139 
140 
141  thePhotonMIPHaloTagger_ = new PhotonMIPHaloTagger();
142  edm::ParameterSet mipVariableSet = conf_.getParameter<edm::ParameterSet>("mipVariableSet");
143  thePhotonMIPHaloTagger_->setup(mipVariableSet,consumesCollector());
144 
145  // Register the product
146  produces< reco::PhotonCollection >(PhotonCollection_);
147 
148 }
149 
150 PhotonProducer::~PhotonProducer()
151 {
152  delete thePhotonEnergyCorrector_;
153  delete thePhotonIsolationCalculator_;
154  delete thePhotonMIPHaloTagger_;
155  //delete energyCorrectionF;
156 }
157 
158 
159 
160 void PhotonProducer::beginRun (edm::Run const& r, edm::EventSetup const & theEventSetup) {
161 
162 
163 
164  thePhotonEnergyCorrector_ -> init(theEventSetup);
165 }
166 
167 void PhotonProducer::endRun (edm::Run const& r, edm::EventSetup const & theEventSetup) {
168 }
169 
170 
171 void PhotonProducer::produce(edm::Event& theEvent, const edm::EventSetup& theEventSetup) {
172 
173  using namespace edm;
174  // nEvt_++;
175 
176  reco::PhotonCollection outputPhotonCollection;
177  std::auto_ptr< reco::PhotonCollection > outputPhotonCollection_p(new reco::PhotonCollection);
178 
179 
180  // Get the PhotonCore collection
181  bool validPhotonCoreHandle=true;
182  Handle<reco::PhotonCoreCollection> photonCoreHandle;
183  theEvent.getByToken(photonCoreProducer_,photonCoreHandle);
184  if (!photonCoreHandle.isValid()) {
185  edm::LogError("PhotonProducer")
186  << "Error! Can't get the photonCoreProducer";
187  validPhotonCoreHandle=false;
188  }
189 
190  // Get EcalRecHits
191  bool validEcalRecHits=true;
192  Handle<EcalRecHitCollection> barrelHitHandle;
193  EcalRecHitCollection barrelRecHits;
194  theEvent.getByToken(barrelEcalHits_, barrelHitHandle);
195  if (!barrelHitHandle.isValid()) {
196  edm::LogError("PhotonProducer")
197  << "Error! Can't get the barrelEcalHits";
198  validEcalRecHits=false;
199  }
200  if ( validEcalRecHits) barrelRecHits = *(barrelHitHandle.product());
201 
202 
203  Handle<EcalRecHitCollection> endcapHitHandle;
204  theEvent.getByToken(endcapEcalHits_, endcapHitHandle);
205  EcalRecHitCollection endcapRecHits;
206  if (!endcapHitHandle.isValid()) {
207  edm::LogError("PhotonProducer")
208  << "Error! Can't get the endcapEcalHits";
209  validEcalRecHits=false;
210  }
211  if( validEcalRecHits) endcapRecHits = *(endcapHitHandle.product());
212 
213  //AA
214  //Get the severity level object
215  edm::ESHandle<EcalSeverityLevelAlgo> sevLv;
216  theEventSetup.get<EcalSeverityLevelAlgoRcd>().get(sevLv);
217  //
218 
219 
220 // get Hcal towers collection
221  Handle<CaloTowerCollection> hcalTowersHandle;
222  theEvent.getByToken(hcalTowers_, hcalTowersHandle);
223 
224 
225  // get the geometry from the event setup:
226  theEventSetup.get<CaloGeometryRecord>().get(theCaloGeom_);
227 
228  //
229  // update energy correction function
230  // energyCorrectionF->init(theEventSetup);
231 
232  edm::ESHandle<CaloTopology> pTopology;
233  theEventSetup.get<CaloTopologyRecord>().get(theCaloTopo_);
234  const CaloTopology *topology = theCaloTopo_.product();
235 
236  // Get the primary event vertex
237  Handle<reco::VertexCollection> vertexHandle;
238  reco::VertexCollection vertexCollection;
239  bool validVertex=true;
240  if ( usePrimaryVertex_ ) {
241  theEvent.getByToken(vertexProducer_, vertexHandle);
242  if (!vertexHandle.isValid()) {
243  edm::LogWarning("PhotonProducer") << "Error! Can't get the product primary Vertex Collection "<< "\n";
244  validVertex=false;
245  }
246  if (validVertex) vertexCollection = *(vertexHandle.product());
247  }
248  // math::XYZPoint vtx(0.,0.,0.);
249  //if (vertexCollection.size()>0) vtx = vertexCollection.begin()->position();
250 
251 
252  int iSC=0; // index in photon collection
253  // Loop over barrel and endcap SC collections and fill the photon collection
254  if ( validPhotonCoreHandle)
255  fillPhotonCollection(theEvent,
256  theEventSetup,
257  photonCoreHandle,
258  topology,
259  &barrelRecHits,
260  &endcapRecHits,
261  hcalTowersHandle,
262  //vtx,
263  vertexCollection,
264  outputPhotonCollection,
265  iSC,
266  sevLv.product());
267 
268 
269  // put the product in the event
270  edm::LogInfo("PhotonProducer") << " Put in the event " << iSC << " Photon Candidates \n";
271  outputPhotonCollection_p->assign(outputPhotonCollection.begin(),outputPhotonCollection.end());
272  theEvent.put( outputPhotonCollection_p, PhotonCollection_);
273 
274 }
275 
276 void PhotonProducer::fillPhotonCollection(edm::Event& evt,
277  edm::EventSetup const & es,
278  const edm::Handle<reco::PhotonCoreCollection> & photonCoreHandle,
279  const CaloTopology* topology,
280  const EcalRecHitCollection* ecalBarrelHits,
281  const EcalRecHitCollection* ecalEndcapHits,
282  const edm::Handle<CaloTowerCollection> & hcalTowersHandle,
283  // math::XYZPoint & vtx,
284  reco::VertexCollection & vertexCollection,
285  reco::PhotonCollection & outputPhotonCollection, int& iSC,
286  const EcalSeverityLevelAlgo * sevLv) {
287 
288  const CaloGeometry* geometry = theCaloGeom_.product();
289  const CaloSubdetectorGeometry* subDetGeometry =0 ;
290  const CaloSubdetectorGeometry* geometryES = theCaloGeom_->getSubdetectorGeometry(DetId::Ecal, EcalPreshower);
291  const EcalRecHitCollection* hits = 0 ;
292  std::vector<double> preselCutValues;
293  float minR9=0;
294 
295 
296  std::vector<int> flags_, severitiesexcl_;
297 
298  for(unsigned int lSC=0; lSC < photonCoreHandle->size(); lSC++) {
299 
300  reco::PhotonCoreRef coreRef(reco::PhotonCoreRef(photonCoreHandle, lSC));
301  reco::SuperClusterRef scRef=coreRef->superCluster();
302  // const reco::SuperCluster* pClus=&(*scRef);
303  iSC++;
304 
305  int subdet = scRef->seed()->hitsAndFractions()[0].first.subdetId();
306  subDetGeometry = theCaloGeom_->getSubdetectorGeometry(DetId::Ecal, subdet);
307 
308  if (subdet==EcalBarrel) {
309  preselCutValues = preselCutValuesBarrel_;
310  minR9 = minR9Barrel_;
311  hits = ecalBarrelHits;
312  flags_ = flagsexclEB_;
313  severitiesexcl_ = severitiesexclEB_;
314  } else if (subdet==EcalEndcap) {
315  preselCutValues = preselCutValuesEndcap_;
316  minR9 = minR9Endcap_;
317  hits = ecalEndcapHits;
318  flags_ = flagsexclEE_;
319  severitiesexcl_ = severitiesexclEE_;
320  } else {
321  edm::LogWarning("")<<"PhotonProducer: do not know if it is a barrel or endcap SuperCluster";
322  }
323 
324 
325  // SC energy preselection
326  if (scRef->energy()/cosh(scRef->eta()) <= preselCutValues[0] ) continue;
327  // calculate HoE
328 
329  const CaloTowerCollection* hcalTowersColl = hcalTowersHandle.product();
330  EgammaTowerIsolation towerIso1(hOverEConeSize_,0.,0.,1,hcalTowersColl) ;
331  EgammaTowerIsolation towerIso2(hOverEConeSize_,0.,0.,2,hcalTowersColl) ;
332  double HoE1=towerIso1.getTowerESum(&(*scRef))/scRef->energy();
333  double HoE2=towerIso2.getTowerESum(&(*scRef))/scRef->energy();
334 
335  EgammaHadTower towerIsoBehindClus(es);
336  towerIsoBehindClus.setTowerCollection(hcalTowersHandle.product());
337  std::vector<CaloTowerDetId> TowersBehindClus = towerIsoBehindClus.towersOf(*scRef);
338  float hcalDepth1OverEcalBc = towerIsoBehindClus.getDepth1HcalESum(TowersBehindClus)/scRef->energy();
339  float hcalDepth2OverEcalBc = towerIsoBehindClus.getDepth2HcalESum(TowersBehindClus)/scRef->energy();
340  // std::cout << " PhotonProducer calculation of HoE with towers in a cone " << HoE1 << " " << HoE2 << std::endl;
341  //std::cout << " PhotonProducer calcualtion of HoE with towers behind the BCs " << hcalDepth1OverEcalBc << " " << hcalDepth2OverEcalBc << std::endl;
342 
343 
344  // recalculate position of seed BasicCluster taking shower depth for unconverted photon
345  math::XYZPoint unconvPos = posCalculator_.Calculate_Location(scRef->seed()->hitsAndFractions(),hits,subDetGeometry,geometryES);
346 
347 
348  float maxXtal = EcalClusterTools::eMax( *(scRef->seed()), &(*hits) );
349  //AA
350  //Change these to consider severity level of hits
351  float e1x5 = EcalClusterTools::e1x5( *(scRef->seed()), &(*hits), &(*topology));
352  float e2x5 = EcalClusterTools::e2x5Max( *(scRef->seed()), &(*hits), &(*topology));
353  float e3x3 = EcalClusterTools::e3x3( *(scRef->seed()), &(*hits), &(*topology));
354  float e5x5 = EcalClusterTools::e5x5( *(scRef->seed()), &(*hits), &(*topology));
355  std::vector<float> cov = EcalClusterTools::covariances( *(scRef->seed()), &(*hits), &(*topology), geometry);
356  std::vector<float> locCov = EcalClusterTools::localCovariances( *(scRef->seed()), &(*hits), &(*topology));
357 
358  float sigmaEtaEta = sqrt(cov[0]);
359  float sigmaIetaIeta = sqrt(locCov[0]);
360  float r9 =e3x3/(scRef->rawEnergy());
361 
362  // compute position of ECAL shower
363  math::XYZPoint caloPosition;
364  if (r9>minR9) {
365  caloPosition = unconvPos;
366  } else {
367  caloPosition = scRef->position();
368  }
369 
371  double photonEnergy=1.;
372  math::XYZPoint vtx(0.,0.,0.);
373  if (vertexCollection.size()>0) vtx = vertexCollection.begin()->position();
374  // compute momentum vector of photon from primary vertex and cluster position
375  math::XYZVector direction = caloPosition - vtx;
376  //math::XYZVector momentum = direction.unit() * photonEnergy ;
377  math::XYZVector momentum = direction.unit() ;
378 
379  // Create dummy candidate with unit momentum and zero energy to allow setting of all variables. The energy is set for last.
380  math::XYZTLorentzVectorD p4(momentum.x(), momentum.y(), momentum.z(), photonEnergy );
381  reco::Photon newCandidate(p4, caloPosition, coreRef, vtx);
382  //std::cout << " standard p4 before " << newCandidate.p4() << " energy " << newCandidate.energy() << std::endl;
383  //std::cout << " type " <<newCandidate.getCandidateP4type() << " standard p4 after " << newCandidate.p4() << " energy " << newCandidate.energy() << std::endl;
384 
385  // Calculate fiducial flags and isolation variable. Blocked are filled from the isolationCalculator
386  reco::Photon::FiducialFlags fiducialFlags;
387  reco::Photon::IsolationVariables isolVarR03, isolVarR04;
388  thePhotonIsolationCalculator_-> calculate ( &newCandidate,evt,es,fiducialFlags,isolVarR04, isolVarR03);
389  newCandidate.setFiducialVolumeFlags( fiducialFlags );
390  newCandidate.setIsolationVariables(isolVarR04, isolVarR03 );
391 
392 
394  reco::Photon::ShowerShape showerShape;
395  showerShape.e1x5= e1x5;
396  showerShape.e2x5= e2x5;
397  showerShape.e3x3= e3x3;
398  showerShape.e5x5= e5x5;
399  showerShape.maxEnergyXtal = maxXtal;
400  showerShape.sigmaEtaEta = sigmaEtaEta;
401  showerShape.sigmaIetaIeta = sigmaIetaIeta;
402  showerShape.hcalDepth1OverEcal = HoE1;
403  showerShape.hcalDepth2OverEcal = HoE2;
404  showerShape.hcalDepth1OverEcalBc = hcalDepth1OverEcalBc;
405  showerShape.hcalDepth2OverEcalBc = hcalDepth2OverEcalBc;
406  showerShape.hcalTowersBehindClusters = TowersBehindClus;
407  newCandidate.setShowerShapeVariables ( showerShape );
408 
411  // Photon candidate takes by default (set in photons_cfi.py) a 4-momentum derived from the ecal photon-specific corrections.
412  thePhotonEnergyCorrector_->calculate(evt, newCandidate, subdet, vertexCollection,es);
413  if ( candidateP4type_ == "fromEcalEnergy") {
414  newCandidate.setP4( newCandidate.p4(reco::Photon::ecal_photons) );
415  newCandidate.setCandidateP4type(reco::Photon::ecal_photons);
416  } else if ( candidateP4type_ == "fromRegression") {
417  newCandidate.setP4( newCandidate.p4(reco::Photon::regression1) );
418  newCandidate.setCandidateP4type(reco::Photon::regression1);
419  }
420 
421  // std::cout << " final p4 " << newCandidate.p4() << " energy " << newCandidate.energy() << std::endl;
422 
423 
424  // std::cout << " PhotonProducer from candidate HoE with towers in a cone " << newCandidate.hadronicOverEm() << " " << newCandidate.hadronicDepth1OverEm() << " " << newCandidate.hadronicDepth2OverEm() << std::endl;
425  // std::cout << " PhotonProducer from candidate of HoE with towers behind the BCs " << newCandidate.hadTowOverEm() << " " << newCandidate.hadTowDepth1OverEm() << " " << newCandidate.hadTowDepth2OverEm() << std::endl;
426 
427 
428  // fill MIP Vairables for Halo: Block for MIP are filled from PhotonMIPHaloTagger
429  reco::Photon::MIPVariables mipVar ;
430  if(subdet==EcalBarrel && runMIPTagger_ )
431  {
432 
433  thePhotonMIPHaloTagger_-> MIPcalculate( &newCandidate,evt,es,mipVar);
434  newCandidate.setMIPVariables(mipVar);
435  }
436 
437 
438 
440  bool isLooseEM=true;
441  if ( newCandidate.pt() < highEt_) {
442  if ( newCandidate.hadronicOverEm() >= preselCutValues[1] ) isLooseEM=false;
443  if ( newCandidate.ecalRecHitSumEtConeDR04() > preselCutValues[2]+ preselCutValues[3]*newCandidate.pt() ) isLooseEM=false;
444  if ( newCandidate.hcalTowerSumEtConeDR04() > preselCutValues[4]+ preselCutValues[5]*newCandidate.pt() ) isLooseEM=false;
445  if ( newCandidate.nTrkSolidConeDR04() > int(preselCutValues[6]) ) isLooseEM=false;
446  if ( newCandidate.nTrkHollowConeDR04() > int(preselCutValues[7]) ) isLooseEM=false;
447  if ( newCandidate.trkSumPtSolidConeDR04() > preselCutValues[8] ) isLooseEM=false;
448  if ( newCandidate.trkSumPtHollowConeDR04() > preselCutValues[9] ) isLooseEM=false;
449  if ( newCandidate.sigmaIetaIeta() > preselCutValues[10] ) isLooseEM=false;
450  }
451 
452 
453 
454  if ( isLooseEM)
455  outputPhotonCollection.push_back(newCandidate);
456 
457 
458  }
459 }
460 
edm::ParameterSet::getParameter
T getParameter(std::string const &) const
PhotonProducer::vertexProducer_
edm::EDGetTokenT< reco::VertexCollection > vertexProducer_
Definition: PhotonProducer.h:69
PhotonProducer::hcalTowers_
edm::EDGetTokenT< CaloTowerCollection > hcalTowers_
Definition: PhotonProducer.h:68
PhotonProducer::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 edm::Handle< CaloTowerCollection > &hcalTowersHandle, reco::VertexCollection &pvVertices, reco::PhotonCollection &outputCollection, int &iSC, const EcalSeverityLevelAlgo *sevLv)
Definition: PhotonProducer.cc:276
PhotonProducer::preselCutValuesEndcap_
std::vector< double > preselCutValuesEndcap_
Definition: PhotonProducer.h:109
reco::Photon::ShowerShape::hcalTowersBehindClusters
std::vector< CaloTowerDetId > hcalTowersBehindClusters
Definition: Photon.h:151
math::XYZTLorentzVectorD
ROOT::Math::LorentzVector< ROOT::Math::PxPyPzE4D< double > > XYZTLorentzVectorD
Lorentz vector with cylindrical internal representation using pseudorapidity.
Definition: LorentzVector.h:14
ecaldqm::topology
CaloTopology const * topology(0)
edm::Event::getByToken
bool getByToken(EDGetToken token, Handle< PROD > &result) const
Definition: Event.h:449
PhotonProducer::conf_
edm::ParameterSet conf_
Definition: PhotonProducer.h:95
init
int init
Definition: HydjetWrapper.h:62
reco::VertexCollection
std::vector< Vertex > VertexCollection
collection of Vertex objects
Definition: VertexFwd.h:9
PhotonProducer::posCalculator_
PositionCalc posCalculator_
Definition: PhotonProducer.h:97
GoodVertex_cfg.vertexCollection
tuple vertexCollection
Definition: GoodVertex_cfg.py:61
AlCaHLTBitMon_QueryRunRegistry.string
string string
Definition: AlCaHLTBitMon_QueryRunRegistry.py:255
PhotonProducer::photonCoreProducer_
edm::EDGetTokenT< reco::PhotonCoreCollection > photonCoreProducer_
Definition: PhotonProducer.h:65
PhotonProducer::hOverEConeSize_
double hOverEConeSize_
Definition: PhotonProducer.h:83
EgammaHadTower::setTowerCollection
void setTowerCollection(const CaloTowerCollection *towercollection)
Definition: EgammaHadTower.cc:110
PhotonProducer::endcapEcalHits_
edm::EDGetTokenT< EcalRecHitCollection > endcapEcalHits_
Definition: PhotonProducer.h:67
PhotonProducer::beginRun
virtual void beginRun(edm::Run const &r, edm::EventSetup const &es) overridefinal
Definition: PhotonProducer.cc:160
PhotonProducer::produce
virtual void produce(edm::Event &evt, const edm::EventSetup &es)
Definition: PhotonProducer.cc:171
HLT_25ns14e33_v3_cff.posCalcParameters
tuple posCalcParameters
Definition: HLT_25ns14e33_v3_cff.py:52748
PhotonProducer::flagsexclEB_
std::vector< int > flagsexclEB_
Definition: PhotonProducer.h:77
EgammaHadTower::getDepth1HcalESum
double getDepth1HcalESum(const reco::SuperCluster &sc) const
Definition: EgammaHadTower.cc:102
edm::Event::put
OrphanHandle< PROD > put(std::auto_ptr< PROD > product)
Put a new product.
Definition: Event.h:113
edm::EDConsumerBase::consumesCollector
ConsumesCollector consumesCollector()
Use a ConsumesCollector to gather consumes information from helper functions.
Definition: EDConsumerBase.cc:71
mathSSE::sqrt
T sqrt(T t)
Definition: SSEVec.h:48
p4
double p4[4]
Definition: TauolaWrapper.h:92
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
PhotonProducer::theCaloTopo_
edm::ESHandle< CaloTopology > theCaloTopo_
Definition: PhotonProducer.h:100
PhotonProducer::flagsexclEE_
std::vector< int > flagsexclEE_
Definition: PhotonProducer.h:78
PhotonProducer::candidateP4type_
std::string candidateP4type_
Definition: PhotonProducer.h:113
PhotonProducer::severitiesexclEE_
std::vector< int > severitiesexclEE_
Definition: PhotonProducer.h:80
PhotonMIPHaloTagger::setup
void setup(const edm::ParameterSet &conf, edm::ConsumesCollector &&iC)
Definition: PhotonMIPHaloTagger.cc:48
EgammaTowerIsolation::getTowerESum
double getTowerESum(const reco::Candidate *cand, const std::vector< CaloTowerDetId > *detIdToExclude=0) const
Definition: EgammaTowerIsolation.h:228
PhotonProducer::thePhotonIsolationCalculator_
PhotonIsolationCalculator * thePhotonIsolationCalculator_
Definition: PhotonProducer.h:103
PhotonProducer::severitiesexclEB_
std::vector< int > severitiesexclEB_
Definition: PhotonProducer.h:79
PhotonProducer::PhotonCollection_
std::string PhotonCollection_
Definition: PhotonProducer.h:64
PhotonProducer::theCaloGeom_
edm::ESHandle< CaloGeometry > theCaloGeom_
Definition: PhotonProducer.h:99
PhotonProducer::thePhotonEnergyCorrector_
PhotonEnergyCorrector * thePhotonEnergyCorrector_
Definition: PhotonProducer.h:112
PhotonEnergyCorrector::calculate
void calculate(edm::Event &evt, reco::Photon &, int subdet, const reco::VertexCollection &vtxcol, const edm::EventSetup &iSetup)
Definition: PhotonEnergyCorrector.cc:102
edm::Handle::product
T const * product() const
Definition: Handle.h:81
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
reco::PhotonCollection
std::vector< Photon > PhotonCollection
collectin of Photon objects
Definition: PhotonFwd.h:9
edm::EventSetup::get
const T & get() const
Definition: EventSetup.h:55
PhotonProducer::preselCutValuesBarrel_
std::vector< double > preselCutValuesBarrel_
Definition: PhotonProducer.h:108
edm::ESHandle::product
T const * product() const
Definition: ESHandle.h:86
EgammaHadTower::towersOf
std::vector< CaloTowerDetId > towersOf(const reco::SuperCluster &sc) const
Definition: EgammaHadTower.cc:31
EgammaHadTower::getDepth2HcalESum
double getDepth2HcalESum(const reco::SuperCluster &sc) const
Definition: EgammaHadTower.cc:106
geometry
ESHandle< TrackerGeometry > geometry
Definition: TkLasBeamFitter.cc:176
PhotonProducer::PhotonProducer
PhotonProducer(const edm::ParameterSet &ps)
Definition: PhotonProducer.cc:38
PhotonProducer::barrelEcalHits_
edm::EDGetTokenT< EcalRecHitCollection > barrelEcalHits_
Definition: PhotonProducer.h:66
PositionCalc::Calculate_Location
math::XYZPoint Calculate_Location(const HitsAndFractions &iDetIds, const edm::SortedCollection< HitType > *iRecHits, const CaloSubdetectorGeometry *iSubGeom, const CaloSubdetectorGeometry *iESGeom=0)
Definition: PositionCalc.h:68
PhotonProducer::thePhotonMIPHaloTagger_
PhotonMIPHaloTagger * thePhotonMIPHaloTagger_
Definition: PhotonProducer.h:106
PhotonProducer::endRun
virtual void endRun(edm::Run const &, edm::EventSetup const &) overridefinal
Definition: PhotonProducer.cc:167
edm::Run
Definition: Run.h:41
Generated for CMSSW Reference Manual by doxygen 1.8.5