CMS 3D CMS Logo
 All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Properties Friends Macros Pages
List of all members | Public Member Functions | Protected Attributes
PhotonIsolationCalculator Class Reference

#include <PhotonIsolationCalculator.h>

Public Member Functions

void calculate (const reco::Photon *, const edm::Event &, const edm::EventSetup &es, reco::Photon::FiducialFlags &phofid, reco::Photon::IsolationVariables &phoisolR03, reco::Photon::IsolationVariables &phoisolR04)
 
double calculateEcalRecHitIso (const reco::Photon *photon, const edm::Event &iEvent, const edm::EventSetup &iSetup, double RCone, double RConeInner, double etaSlice, double eMin, double etMin, bool vetoClusteredHits, bool useNumCrystals)
 
double calculateHcalTowerIso (const reco::Photon *photon, const edm::Event &iEvent, const edm::EventSetup &iSetup, double RCone, double RConeInner, double eMin, signed int depth)
 
void calculateTrackIso (const reco::Photon *photon, const edm::Event &e, double &trkCone, int &ntrkCone, double pTThresh=0, double RCone=.4, double RinnerCone=.1, double etaSlice=0.015, double lip=0.2, double d0=0.1)
 
void classify (const reco::Photon *photon, bool &isEBPho, bool &isEEPho, bool &isEBEtaGap, bool &isEBPhiGap, bool &isEERingGap, bool &isEEDeeGap, bool &isEBEEGap)
 
 PhotonIsolationCalculator ()
 
void setup (const edm::ParameterSet &conf)
 
virtual ~PhotonIsolationCalculator ()
 

Protected Attributes

std::string barrelecalCollection_
 
std::string barrelecalProducer_
 
edm::InputTag beamSpotProducerTag_
 
std::vector< double > ecalIsoBarrelRadiusA_
 
std::vector< double > ecalIsoBarrelRadiusB_
 
std::vector< double > ecalIsoEndcapRadiusA_
 
std::vector< double > ecalIsoEndcapRadiusB_
 
std::string endcapecalCollection_
 
std::string endcapecalProducer_
 
std::string hcalCollection_
 
std::vector< double > hcalIsoBarrelRadiusA_
 
std::vector< double > hcalIsoBarrelRadiusB_
 
std::vector< double > hcalIsoEndcapRadiusA_
 
std::vector< double > hcalIsoEndcapRadiusB_
 
std::string hcalProducer_
 
double isolationtrackEtaSliceA_
 
double isolationtrackEtaSliceB_
 
double isolationtrackThresholdA_
 
double isolationtrackThresholdB_
 
std::vector< double > moduleEtaBoundary_
 
double modulePhiBoundary_
 
double photonEcalRecHitConeInnerRadiusA_
 
double photonEcalRecHitConeInnerRadiusB_
 
double photonEcalRecHitConeOuterRadiusA_
 
double photonEcalRecHitConeOuterRadiusB_
 
double photonEcalRecHitEtaSliceA_
 
double photonEcalRecHitEtaSliceB_
 
double photonEcalRecHitThreshEA_
 
double photonEcalRecHitThreshEB_
 
double photonEcalRecHitThreshEtA_
 
double photonEcalRecHitThreshEtB_
 
double photonHcalDepth1TowerConeInnerRadiusA_
 
double photonHcalDepth1TowerConeInnerRadiusB_
 
double photonHcalDepth1TowerConeOuterRadiusA_
 
double photonHcalDepth1TowerConeOuterRadiusB_
 
double photonHcalDepth1TowerThreshEA_
 
double photonHcalDepth1TowerThreshEB_
 
double photonHcalDepth2TowerConeInnerRadiusA_
 
double photonHcalDepth2TowerConeInnerRadiusB_
 
double photonHcalDepth2TowerConeOuterRadiusA_
 
double photonHcalDepth2TowerConeOuterRadiusB_
 
double photonHcalDepth2TowerThreshEA_
 
double photonHcalDepth2TowerThreshEB_
 
double photonHcalTowerConeInnerRadiusA_
 
double photonHcalTowerConeInnerRadiusB_
 
double photonHcalTowerConeOuterRadiusA_
 
double photonHcalTowerConeOuterRadiusB_
 
double photonHcalTowerThreshEA_
 
double photonHcalTowerThreshEB_
 
int severityLevelCut_
 
float severityRecHitThreshold_
 
EcalSeverityLevelAlgo::SpikeId spId_
 
float spikeIdThreshold_
 
double trackConeInnerRadiusA_
 
double trackConeInnerRadiusB_
 
double trackConeOuterRadiusA_
 
double trackConeOuterRadiusB_
 
double trackD0RadiusA_
 
double trackD0RadiusB_
 
edm::InputTag trackInputTag_
 
double trackLipRadiusA_
 
double trackLipRadiusB_
 
std::vector< double > trkIsoBarrelRadiusA_
 
std::vector< double > trkIsoBarrelRadiusB_
 
std::vector< double > trkIsoEndcapRadiusA_
 
std::vector< double > trkIsoEndcapRadiusB_
 
bool useNumCrystals_
 
std::vector< int > v_chstatus_
 
bool vetoClusteredEcalHits_
 

Detailed Description

Determine and Set quality information on Photon Objects

Author
A. Askew, N. Marinelli, M.B. Anderson

Definition at line 14 of file PhotonIsolationCalculator.h.

Constructor & Destructor Documentation

PhotonIsolationCalculator::PhotonIsolationCalculator ( )
inline

Definition at line 18 of file PhotonIsolationCalculator.h.

18 {};
virtual PhotonIsolationCalculator::~PhotonIsolationCalculator ( )
inlinevirtual

Definition at line 20 of file PhotonIsolationCalculator.h.

20 {};

Member Function Documentation

void PhotonIsolationCalculator::calculate ( const reco::Photon *  pho,
const edm::Event &  e,
const edm::EventSetup &  es,
reco::Photon::FiducialFlags &  phofid,
reco::Photon::IsolationVariables &  phoisolR03,
reco::Photon::IsolationVariables &  phoisolR04 
)

Definition at line 188 of file PhotonIsolationCalculator.cc.

References calculateEcalRecHitIso(), calculateHcalTowerIso(), calculateTrackIso(), classify(), ExpressReco_HICollisions_FallBack::detector, EcalBarrel, EcalEndcap, ecalIsoBarrelRadiusA_, ecalIsoBarrelRadiusB_, ecalIsoEndcapRadiusA_, ecalIsoEndcapRadiusB_, reco::Photon::IsolationVariables::ecalRecHitSumEt, reco::Photon::IsolationVariables::hcalDepth1TowerSumEt, reco::Photon::IsolationVariables::hcalDepth2TowerSumEt, hcalIsoBarrelRadiusA_, hcalIsoBarrelRadiusB_, hcalIsoEndcapRadiusA_, hcalIsoEndcapRadiusB_, reco::Photon::IsolationVariables::hcalTowerSumEt, reco::Photon::FiducialFlags::isEB, reco::Photon::FiducialFlags::isEBEEGap, reco::Photon::FiducialFlags::isEBEtaGap, reco::Photon::FiducialFlags::isEBPhiGap, reco::Photon::FiducialFlags::isEE, reco::Photon::FiducialFlags::isEEDeeGap, reco::Photon::FiducialFlags::isEERingGap, isolationtrackEtaSliceA_, isolationtrackEtaSliceB_, isolationtrackThresholdA_, isolationtrackThresholdB_, reco::Photon::IsolationVariables::nTrkHollowCone, reco::Photon::IsolationVariables::nTrkSolidCone, photonEcalRecHitConeInnerRadiusA_, photonEcalRecHitConeInnerRadiusB_, photonEcalRecHitConeOuterRadiusA_, photonEcalRecHitConeOuterRadiusB_, photonEcalRecHitEtaSliceA_, photonEcalRecHitEtaSliceB_, photonEcalRecHitThreshEA_, photonEcalRecHitThreshEB_, photonEcalRecHitThreshEtA_, photonEcalRecHitThreshEtB_, photonHcalDepth1TowerConeInnerRadiusA_, photonHcalDepth1TowerConeInnerRadiusB_, photonHcalDepth1TowerConeOuterRadiusA_, photonHcalDepth1TowerConeOuterRadiusB_, photonHcalDepth1TowerThreshEA_, photonHcalDepth1TowerThreshEB_, photonHcalDepth2TowerConeInnerRadiusA_, photonHcalDepth2TowerConeInnerRadiusB_, photonHcalDepth2TowerConeOuterRadiusA_, photonHcalDepth2TowerConeOuterRadiusB_, photonHcalDepth2TowerThreshEA_, photonHcalDepth2TowerThreshEB_, photonHcalTowerConeInnerRadiusA_, photonHcalTowerConeInnerRadiusB_, photonHcalTowerConeOuterRadiusA_, photonHcalTowerConeOuterRadiusB_, photonHcalTowerThreshEA_, photonHcalTowerThreshEB_, DetId::subdetId(), reco::Photon::superCluster(), trackConeInnerRadiusA_, trackConeInnerRadiusB_, trackConeOuterRadiusA_, trackConeOuterRadiusB_, trackD0RadiusA_, trackD0RadiusB_, trackLipRadiusA_, trackLipRadiusB_, trkIsoBarrelRadiusA_, trkIsoBarrelRadiusB_, trkIsoEndcapRadiusA_, trkIsoEndcapRadiusB_, reco::Photon::IsolationVariables::trkSumPtHollowCone, reco::Photon::IsolationVariables::trkSumPtSolidCone, useNumCrystals_, and vetoClusteredEcalHits_.

193  {
194 
195 
196  //Get fiducial flags. This does not really belong here
197  bool isEBPho = false;
198  bool isEEPho = false;
199  bool isEBEtaGap = false;
200  bool isEBPhiGap = false;
201  bool isEERingGap = false;
202  bool isEEDeeGap = false;
203  bool isEBEEGap = false;
204  classify(pho, isEBPho, isEEPho, isEBEtaGap, isEBPhiGap, isEERingGap, isEEDeeGap, isEBEEGap);
205  phofid.isEB = isEBPho;
206  phofid.isEE = isEEPho;
207  phofid.isEBEtaGap = isEBEtaGap;
208  phofid.isEBPhiGap = isEBPhiGap;
209  phofid.isEERingGap = isEERingGap;
210  phofid.isEEDeeGap = isEEDeeGap;
211  phofid.isEBEEGap = isEBEEGap;
212 
213  // Calculate isolation variables. cone sizes and thresholds
214  // are set for Barrel and endcap separately
215 
216  reco::SuperClusterRef scRef=pho->superCluster();
217  const reco::BasicCluster & seedCluster = *(scRef->seed()) ;
218  DetId seedXtalId = seedCluster.hitsAndFractions()[0].first ;
219  int detector = seedXtalId.subdetId() ;
220 
221  if (detector==EcalBarrel) {
222 
223  trackConeOuterRadiusA_ = trkIsoBarrelRadiusA_[0];
224  trackConeInnerRadiusA_ = trkIsoBarrelRadiusA_[1];
225  isolationtrackThresholdA_ = trkIsoBarrelRadiusA_[2];
226  trackLipRadiusA_ = trkIsoBarrelRadiusA_[3];
227  trackD0RadiusA_ = trkIsoBarrelRadiusA_[4];
228  isolationtrackEtaSliceA_ = trkIsoBarrelRadiusA_[5];
229 
230  photonEcalRecHitConeInnerRadiusA_ = ecalIsoBarrelRadiusA_[0];
231  photonEcalRecHitConeOuterRadiusA_ = ecalIsoBarrelRadiusA_[1];
232  photonEcalRecHitEtaSliceA_ = ecalIsoBarrelRadiusA_[2];
233  photonEcalRecHitThreshEA_ = ecalIsoBarrelRadiusA_[3];
234  photonEcalRecHitThreshEtA_ = ecalIsoBarrelRadiusA_[4];
235 
236  photonHcalTowerConeInnerRadiusA_ = hcalIsoBarrelRadiusA_[0];
237  photonHcalTowerConeOuterRadiusA_ = hcalIsoBarrelRadiusA_[1];
238  photonHcalTowerThreshEA_ = hcalIsoBarrelRadiusA_[2];
239  photonHcalDepth1TowerConeInnerRadiusA_ = hcalIsoBarrelRadiusA_[3];
240  photonHcalDepth1TowerConeOuterRadiusA_ = hcalIsoBarrelRadiusA_[4];
241  photonHcalDepth1TowerThreshEA_ = hcalIsoBarrelRadiusA_[5];
242  photonHcalDepth2TowerConeInnerRadiusA_ = hcalIsoBarrelRadiusA_[6];
243  photonHcalDepth2TowerConeOuterRadiusA_ = hcalIsoBarrelRadiusA_[7];
244  photonHcalDepth2TowerThreshEA_ = hcalIsoBarrelRadiusA_[8];
245 
246 
247  trackConeOuterRadiusB_ = trkIsoBarrelRadiusB_[0];
248  trackConeInnerRadiusB_ = trkIsoBarrelRadiusB_[1];
249  isolationtrackThresholdB_ = trkIsoBarrelRadiusB_[2];
250  trackLipRadiusB_ = trkIsoBarrelRadiusB_[3];
251  trackD0RadiusB_ = trkIsoBarrelRadiusB_[4];
252  isolationtrackEtaSliceB_ = trkIsoBarrelRadiusB_[5];
253 
254  photonEcalRecHitConeInnerRadiusB_ = ecalIsoBarrelRadiusB_[0];
255  photonEcalRecHitConeOuterRadiusB_ = ecalIsoBarrelRadiusB_[1];
256  photonEcalRecHitEtaSliceB_ = ecalIsoBarrelRadiusB_[2];
257  photonEcalRecHitThreshEB_ = ecalIsoBarrelRadiusB_[3];
258  photonEcalRecHitThreshEtB_ = ecalIsoBarrelRadiusB_[4];
259 
260  photonHcalTowerConeInnerRadiusB_ = hcalIsoBarrelRadiusB_[0];
261  photonHcalTowerConeOuterRadiusB_ = hcalIsoBarrelRadiusB_[1];
262  photonHcalTowerThreshEB_ = hcalIsoBarrelRadiusB_[2];
263  photonHcalDepth1TowerConeInnerRadiusB_ = hcalIsoBarrelRadiusB_[3];
264  photonHcalDepth1TowerConeOuterRadiusB_ = hcalIsoBarrelRadiusB_[4];
265  photonHcalDepth1TowerThreshEB_ = hcalIsoBarrelRadiusB_[5];
266  photonHcalDepth2TowerConeInnerRadiusB_ = hcalIsoBarrelRadiusB_[6];
267  photonHcalDepth2TowerConeOuterRadiusB_ = hcalIsoBarrelRadiusB_[7];
268  photonHcalDepth2TowerThreshEB_ = hcalIsoBarrelRadiusB_[8];
269 
270 
271 
272 
273  } else if
274  (detector==EcalEndcap) {
275 
276  trackConeOuterRadiusA_ = trkIsoEndcapRadiusA_[0];
277  trackConeInnerRadiusA_ = trkIsoEndcapRadiusA_[1];
278  isolationtrackThresholdA_ = trkIsoEndcapRadiusA_[2];
279  trackLipRadiusA_ = trkIsoEndcapRadiusA_[3];
280  trackD0RadiusA_ = trkIsoEndcapRadiusA_[4];
281  isolationtrackEtaSliceA_ = trkIsoEndcapRadiusA_[5];
282 
283  photonEcalRecHitConeInnerRadiusA_ = ecalIsoEndcapRadiusA_[0];
284  photonEcalRecHitConeOuterRadiusA_ = ecalIsoEndcapRadiusA_[1];
285  photonEcalRecHitEtaSliceA_ = ecalIsoEndcapRadiusA_[2];
286  photonEcalRecHitThreshEA_ = ecalIsoEndcapRadiusA_[3];
287  photonEcalRecHitThreshEtA_ = ecalIsoEndcapRadiusA_[4];
288 
289  photonHcalTowerConeInnerRadiusA_ = hcalIsoEndcapRadiusA_[0];
290  photonHcalTowerConeOuterRadiusA_ = hcalIsoEndcapRadiusA_[1];
291  photonHcalTowerThreshEA_ = hcalIsoEndcapRadiusA_[2];
292  photonHcalDepth1TowerConeInnerRadiusA_ = hcalIsoEndcapRadiusA_[3];
293  photonHcalDepth1TowerConeOuterRadiusA_ = hcalIsoEndcapRadiusA_[4];
294  photonHcalDepth1TowerThreshEA_ = hcalIsoEndcapRadiusA_[5];
295  photonHcalDepth2TowerConeInnerRadiusA_ = hcalIsoEndcapRadiusA_[6];
296  photonHcalDepth2TowerConeOuterRadiusA_ = hcalIsoEndcapRadiusA_[7];
297  photonHcalDepth2TowerThreshEA_ = hcalIsoEndcapRadiusA_[8];
298 
299 
300  trackConeOuterRadiusB_ = trkIsoEndcapRadiusB_[0];
301  trackConeInnerRadiusB_ = trkIsoEndcapRadiusB_[1];
302  isolationtrackThresholdB_ = trkIsoEndcapRadiusB_[2];
303  trackLipRadiusB_ = trkIsoEndcapRadiusB_[3];
304  trackD0RadiusB_ = trkIsoEndcapRadiusB_[4];
305  isolationtrackEtaSliceB_ = trkIsoEndcapRadiusB_[5];
306 
307 
308  photonEcalRecHitConeInnerRadiusB_ = ecalIsoEndcapRadiusB_[0];
309  photonEcalRecHitConeOuterRadiusB_ = ecalIsoEndcapRadiusB_[1];
310  photonEcalRecHitEtaSliceB_ = ecalIsoEndcapRadiusB_[2];
311  photonEcalRecHitThreshEB_ = ecalIsoEndcapRadiusB_[3];
312  photonEcalRecHitThreshEtB_ = ecalIsoEndcapRadiusB_[4];
313 
314  photonHcalTowerConeInnerRadiusB_ = hcalIsoEndcapRadiusB_[0];
315  photonHcalTowerConeOuterRadiusB_ = hcalIsoEndcapRadiusB_[1];
316  photonHcalTowerThreshEB_ = hcalIsoEndcapRadiusB_[2];
317  photonHcalDepth1TowerConeInnerRadiusB_ = hcalIsoEndcapRadiusB_[3];
318  photonHcalDepth1TowerConeOuterRadiusB_ = hcalIsoEndcapRadiusB_[4];
319  photonHcalDepth1TowerThreshEB_ = hcalIsoEndcapRadiusB_[5];
320  photonHcalDepth2TowerConeInnerRadiusB_ = hcalIsoEndcapRadiusB_[6];
321  photonHcalDepth2TowerConeOuterRadiusB_ = hcalIsoEndcapRadiusB_[7];
322  photonHcalDepth2TowerThreshEB_ = hcalIsoEndcapRadiusB_[8];
323 
324  }
325 
326 
327  //Calculate hollow cone track isolation, CONE A
328  int ntrkA=0;
329  double trkisoA=0;
330  calculateTrackIso(pho, e, trkisoA, ntrkA, isolationtrackThresholdA_,
331  trackConeOuterRadiusA_, trackConeInnerRadiusA_, isolationtrackEtaSliceA_, trackLipRadiusA_, trackD0RadiusA_);
332 
333  //Calculate solid cone track isolation, CONE A
334  int sntrkA=0;
335  double strkisoA=0;
336  calculateTrackIso(pho, e, strkisoA, sntrkA, isolationtrackThresholdA_,
337  trackConeOuterRadiusA_, 0., isolationtrackEtaSliceA_, trackLipRadiusA_, trackD0RadiusA_ );
338 
339  phoisolR1.nTrkHollowCone = ntrkA;
340  phoisolR1.trkSumPtHollowCone = trkisoA;
341  phoisolR1.nTrkSolidCone = sntrkA;
342  phoisolR1.trkSumPtSolidCone = strkisoA;
343 
344  //Calculate hollow cone track isolation, CONE B
345  int ntrkB=0;
346  double trkisoB=0;
347  calculateTrackIso(pho, e, trkisoB, ntrkB, isolationtrackThresholdB_,
348  trackConeOuterRadiusB_, trackConeInnerRadiusB_, isolationtrackEtaSliceB_, trackLipRadiusB_, trackD0RadiusB_ );
349 
350  //Calculate solid cone track isolation, CONE B
351  int sntrkB=0;
352  double strkisoB=0;
353  calculateTrackIso(pho, e, strkisoB, sntrkB, isolationtrackThresholdB_,
354  trackConeOuterRadiusB_, 0., isolationtrackEtaSliceB_, trackLipRadiusB_, trackD0RadiusB_);
355 
356  phoisolR2.nTrkHollowCone = ntrkB;
357  phoisolR2.trkSumPtHollowCone = trkisoB;
358  phoisolR2.nTrkSolidCone = sntrkB;
359  phoisolR2.trkSumPtSolidCone = strkisoB;
360 
361 // std::cout << "Output from solid cone track isolation: ";
362 // std::cout << " Sum pT: " << strkiso << " ntrk: " << sntrk << std::endl;
363 
364  double EcalRecHitIsoA = calculateEcalRecHitIso(pho, e, es,
365  photonEcalRecHitConeOuterRadiusA_,
366  photonEcalRecHitConeInnerRadiusA_,
367  photonEcalRecHitEtaSliceA_,
368  photonEcalRecHitThreshEA_,
369  photonEcalRecHitThreshEtA_,
370  vetoClusteredEcalHits_,
371  useNumCrystals_);
372  phoisolR1.ecalRecHitSumEt = EcalRecHitIsoA;
373 
374  double EcalRecHitIsoB = calculateEcalRecHitIso(pho, e, es,
375  photonEcalRecHitConeOuterRadiusB_,
376  photonEcalRecHitConeInnerRadiusB_,
377  photonEcalRecHitEtaSliceB_,
378  photonEcalRecHitThreshEB_,
379  photonEcalRecHitThreshEtB_,
380  vetoClusteredEcalHits_,
381  useNumCrystals_);
382  phoisolR2.ecalRecHitSumEt = EcalRecHitIsoB;
383 
384  double HcalTowerIsoA = calculateHcalTowerIso(pho, e, es, photonHcalTowerConeOuterRadiusA_,
385  photonHcalTowerConeInnerRadiusA_,
386  photonHcalTowerThreshEA_, -1 );
387  phoisolR1.hcalTowerSumEt = HcalTowerIsoA;
388 
389 
390  double HcalTowerIsoB = calculateHcalTowerIso(pho, e, es, photonHcalTowerConeOuterRadiusB_,
391  photonHcalTowerConeInnerRadiusB_,
392  photonHcalTowerThreshEB_, -1 );
393  phoisolR2.hcalTowerSumEt = HcalTowerIsoB;
394 
396 
397  double HcalDepth1TowerIsoA = calculateHcalTowerIso(pho, e, es, photonHcalDepth1TowerConeOuterRadiusA_,
398  photonHcalDepth1TowerConeInnerRadiusA_,
399  photonHcalDepth1TowerThreshEA_, 1 );
400  phoisolR1.hcalDepth1TowerSumEt = HcalDepth1TowerIsoA;
401 
402 
403  double HcalDepth1TowerIsoB = calculateHcalTowerIso(pho, e, es, photonHcalDepth1TowerConeOuterRadiusB_,
404  photonHcalDepth1TowerConeInnerRadiusB_,
405  photonHcalDepth1TowerThreshEB_, 1 );
406  phoisolR2.hcalDepth1TowerSumEt = HcalDepth1TowerIsoB;
407 
408 
409 
411 
412  double HcalDepth2TowerIsoA = calculateHcalTowerIso(pho, e, es, photonHcalDepth2TowerConeOuterRadiusA_,
413  photonHcalDepth2TowerConeInnerRadiusA_,
414  photonHcalDepth2TowerThreshEA_, 2 );
415  phoisolR1.hcalDepth2TowerSumEt = HcalDepth2TowerIsoA;
416 
417 
418  double HcalDepth2TowerIsoB = calculateHcalTowerIso(pho, e, es, photonHcalDepth2TowerConeOuterRadiusB_,
419  photonHcalDepth2TowerConeInnerRadiusB_,
420  photonHcalDepth2TowerThreshEB_, 2 );
421  phoisolR2.hcalDepth2TowerSumEt = HcalDepth2TowerIsoB;
422 
423 
424 
425 
426 
427 
428 }
PhotonIsolationCalculator::ecalIsoBarrelRadiusA_
std::vector< double > ecalIsoBarrelRadiusA_
Definition: PhotonIsolationCalculator.h:91
reco::Photon::superCluster
reco::SuperClusterRef superCluster() const
Retrieve photonCore attributes.
Definition: Photon.h:54
PhotonIsolationCalculator::hcalIsoBarrelRadiusB_
std::vector< double > hcalIsoBarrelRadiusB_
Definition: PhotonIsolationCalculator.h:95
PhotonIsolationCalculator::ecalIsoEndcapRadiusB_
std::vector< double > ecalIsoEndcapRadiusB_
Definition: PhotonIsolationCalculator.h:101
PhotonIsolationCalculator::trkIsoEndcapRadiusB_
std::vector< double > trkIsoEndcapRadiusB_
Definition: PhotonIsolationCalculator.h:100
PhotonIsolationCalculator::calculateTrackIso
void calculateTrackIso(const reco::Photon *photon, const edm::Event &e, double &trkCone, int &ntrkCone, double pTThresh=0, double RCone=.4, double RinnerCone=.1, double etaSlice=0.015, double lip=0.2, double d0=0.1)
Definition: PhotonIsolationCalculator.cc:483
PhotonIsolationCalculator::calculateHcalTowerIso
double calculateHcalTowerIso(const reco::Photon *photon, const edm::Event &iEvent, const edm::EventSetup &iSetup, double RCone, double RConeInner, double eMin, signed int depth)
Definition: PhotonIsolationCalculator.cc:598
PhotonIsolationCalculator::hcalIsoEndcapRadiusA_
std::vector< double > hcalIsoEndcapRadiusA_
Definition: PhotonIsolationCalculator.h:99
PhotonIsolationCalculator::trkIsoBarrelRadiusB_
std::vector< double > trkIsoBarrelRadiusB_
Definition: PhotonIsolationCalculator.h:93
PhotonIsolationCalculator::hcalIsoBarrelRadiusA_
std::vector< double > hcalIsoBarrelRadiusA_
Definition: PhotonIsolationCalculator.h:92
DetId::subdetId
int subdetId() const
get the contents of the subdetector field (not cast into any detector&#39;s numbering enum) ...
Definition: DetId.h:39
DetId
Definition: DetId.h:20
ExpressReco_HICollisions_FallBack.detector
tuple detector
Definition: ExpressReco_HICollisions_FallBack.py:25265
PhotonIsolationCalculator::trkIsoEndcapRadiusA_
std::vector< double > trkIsoEndcapRadiusA_
Definition: PhotonIsolationCalculator.h:97
PhotonIsolationCalculator::hcalIsoEndcapRadiusB_
std::vector< double > hcalIsoEndcapRadiusB_
Definition: PhotonIsolationCalculator.h:102
PhotonIsolationCalculator::ecalIsoEndcapRadiusA_
std::vector< double > ecalIsoEndcapRadiusA_
Definition: PhotonIsolationCalculator.h:98
PhotonIsolationCalculator::ecalIsoBarrelRadiusB_
std::vector< double > ecalIsoBarrelRadiusB_
Definition: PhotonIsolationCalculator.h:94
PhotonIsolationCalculator::calculateEcalRecHitIso
double calculateEcalRecHitIso(const reco::Photon *photon, const edm::Event &iEvent, const edm::EventSetup &iSetup, double RCone, double RConeInner, double etaSlice, double eMin, double etMin, bool vetoClusteredHits, bool useNumCrystals)
Definition: PhotonIsolationCalculator.cc:527
PhotonIsolationCalculator::classify
void classify(const reco::Photon *photon, bool &isEBPho, bool &isEEPho, bool &isEBEtaGap, bool &isEBPhiGap, bool &isEERingGap, bool &isEEDeeGap, bool &isEBEEGap)
Definition: PhotonIsolationCalculator.cc:431
PhotonIsolationCalculator::trkIsoBarrelRadiusA_
std::vector< double > trkIsoBarrelRadiusA_
Definition: PhotonIsolationCalculator.h:90
double PhotonIsolationCalculator::calculateEcalRecHitIso ( const reco::Photon *  photon,
const edm::Event &  iEvent,
const edm::EventSetup &  iSetup,
double  RCone,
double  RConeInner,
double  etaSlice,
double  eMin,
double  etMin,
bool  vetoClusteredHits,
bool  useNumCrystals 
)

Definition at line 527 of file PhotonIsolationCalculator.cc.

References barrelecalCollection_, barrelecalProducer_, EgammaRecHitIsolation::doSpikeRemoval(), DetId::Ecal, endcapecalCollection_, endcapecalProducer_, edm::EventSetup::get(), edm::Event::getByLabel(), EgammaRecHitIsolation::getEtSum(), edm::Handle< T >::product(), edm::ESHandle< class >::product(), EgammaRecHitIsolation::setUseNumCrystals(), EgammaRecHitIsolation::setVetoClustered(), severityLevelCut_, severityRecHitThreshold_, spId_, and spikeIdThreshold_.

Referenced by calculate().

537 {
538 
539  edm::Handle<EcalRecHitCollection> ecalhitsCollEB;
540  edm::Handle<EcalRecHitCollection> ecalhitsCollEE;
541 
542  iEvent.getByLabel(endcapecalProducer_,endcapecalCollection_, ecalhitsCollEE);
543 
544  iEvent.getByLabel(barrelecalProducer_,barrelecalCollection_, ecalhitsCollEB);
545 
546  const EcalRecHitCollection* rechitsCollectionEE_ = ecalhitsCollEE.product();
547  const EcalRecHitCollection* rechitsCollectionEB_ = ecalhitsCollEB.product();
548 
549  //Get the channel status from the db
550  edm::ESHandle<EcalChannelStatus> chStatus;
551  iSetup.get<EcalChannelStatusRcd>().get(chStatus);
552 
553  std::auto_ptr<CaloRecHitMetaCollectionV> RecHitsEE(0);
554  RecHitsEE = std::auto_ptr<CaloRecHitMetaCollectionV>(new EcalRecHitMetaCollection(*rechitsCollectionEE_));
555 
556  std::auto_ptr<CaloRecHitMetaCollectionV> RecHitsEB(0);
557  RecHitsEB = std::auto_ptr<CaloRecHitMetaCollectionV>(new EcalRecHitMetaCollection(*rechitsCollectionEB_));
558 
559  edm::ESHandle<CaloGeometry> geoHandle;
560  iSetup.get<CaloGeometryRecord>().get(geoHandle);
561 
562  EgammaRecHitIsolation phoIsoEB(RCone,
563  RConeInner,
564  etaSlice,
565  etMin,
566  eMin,
567  geoHandle,
568  &(*RecHitsEB),
569  DetId::Ecal);
570 
571  phoIsoEB.setVetoClustered(vetoClusteredHits);
572  phoIsoEB.setUseNumCrystals(useNumXtals);
573  phoIsoEB.doSpikeRemoval(ecalhitsCollEB.product(),chStatus.product(),severityLevelCut_,severityRecHitThreshold_,spId_,spikeIdThreshold_);
574 
575  double ecalIsolEB = phoIsoEB.getEtSum(photon);
576 
577  EgammaRecHitIsolation phoIsoEE(RCone,
578  RConeInner,
579  etaSlice,
580  etMin,
581  eMin,
582  geoHandle,
583  &(*RecHitsEE),
584  DetId::Ecal);
585 
586  phoIsoEE.setVetoClustered(vetoClusteredHits);
587  phoIsoEE.setUseNumCrystals(useNumXtals);
588 
589  double ecalIsolEE = phoIsoEE.getEtSum(photon);
590  // delete phoIso;
591  double ecalIsol = ecalIsolEB + ecalIsolEE;
592 
593  return ecalIsol;
594 
595 
596 }
ExpressReco_HICollisions_FallBack.etMin
tuple etMin
Definition: ExpressReco_HICollisions_FallBack.py:25165
PhotonIsolationCalculator::spId_
EcalSeverityLevelAlgo::SpikeId spId_
Definition: PhotonIsolationCalculator.h:129
edm::Event::getByLabel
bool getByLabel(InputTag const &tag, Handle< PROD > &result) const
Definition: Event.h:359
edm::EventSetup::get
const T & get() const
Definition: EventSetup.h:55
edm::ESHandle::product
T const * product() const
Definition: ESHandle.h:62
edm::Handle::product
T const * product() const
Definition: Handle.h:74
EcalRecHitMetaCollection
SimpleCaloRecHitMetaCollection< EcalRecHitCollection > EcalRecHitMetaCollection
Definition: CaloRecHitMetaCollections.h:34
double PhotonIsolationCalculator::calculateHcalTowerIso ( const reco::Photon *  photon,
const edm::Event &  iEvent,
const edm::EventSetup &  iSetup,
double  RCone,
double  RConeInner,
double  eMin,
signed int  depth 
)

Definition at line 598 of file PhotonIsolationCalculator.cc.

References edm::Event::getByLabel(), EgammaTowerIsolation::getTowerEtSum(), hcalCollection_, hcalProducer_, and edm::Handle< T >::product().

Referenced by calculate().

605 {
606 
607  edm::Handle<CaloTowerCollection> hcalhitsCollH;
608 
609  iEvent.getByLabel(hcalProducer_,hcalCollection_, hcalhitsCollH);
610 
611  const CaloTowerCollection *toww = hcalhitsCollH.product();
612 
613  double ecalIsol=0.;
614 
615  //std::cout << "before iso call" << std::endl;
616  EgammaTowerIsolation phoIso(RCone,
617  RConeInner,
618  eMin,depth,
619  toww);
620  ecalIsol = phoIso.getTowerEtSum(photon);
621  // delete phoIso;
622  //std::cout << "after call" << std::endl;
623  return ecalIsol;
624 
625 
626 }
edm::Event::getByLabel
bool getByLabel(InputTag const &tag, Handle< PROD > &result) const
Definition: Event.h:359
edm::Handle::product
T const * product() const
Definition: Handle.h:74
void PhotonIsolationCalculator::calculateTrackIso ( const reco::Photon *  photon,
const edm::Event &  e,
double &  trkCone,
int &  ntrkCone,
double  pTThresh = 0,
double  RCone = .4,
double  RinnerCone = .1,
double  etaSlice = 0.015,
double  lip = 0.2,
double  d0 = 0.1 
)

Definition at line 483 of file PhotonIsolationCalculator.cc.

References beamSpotProducerTag_, cmsDriverOptions::counter, edm::Event::getByLabel(), PhotonTkIsolation::getNumberTracks(), PhotonTkIsolation::getPtTracks(), edm::HandleBase::isValid(), edm::Handle< T >::product(), ExpressReco_HICollisions_FallBack::trackCollection, trackInputTag_, testEve_cfg::tracks, reco::BeamSpot::x0(), reco::BeamSpot::y0(), and reco::BeamSpot::z0().

Referenced by calculate().

492  {
493  int counter =0;
494  double ptSum =0.;
495 
496 
497  //get the tracks
498  edm::Handle<reco::TrackCollection> tracks;
499  e.getByLabel(trackInputTag_,tracks);
500  if(!tracks.isValid()) {
501  return;
502  }
503  const reco::TrackCollection* trackCollection = tracks.product();
504  //Photon Eta and Phi. Hope these are correct.
505  reco::BeamSpot vertexBeamSpot;
506  edm::Handle<reco::BeamSpot> recoBeamSpotHandle;
507  e.getByLabel(beamSpotProducerTag_,recoBeamSpotHandle);
508  vertexBeamSpot = *recoBeamSpotHandle;
509 
510  PhotonTkIsolation phoIso(RCone,
511  RinnerCone,
512  etaSlice,
513  pTThresh,
514  lip ,
515  d0,
516  trackCollection,
517  math::XYZPoint(vertexBeamSpot.x0(),vertexBeamSpot.y0(),vertexBeamSpot.z0()));
518 
519  counter = phoIso.getNumberTracks(photon);
520  ptSum = phoIso.getPtTracks(photon);
521  ntrkCone = counter;
522  trkCone = ptSum;
523 }
reco::BeamSpot::z0
double z0() const
z coordinate
Definition: BeamSpot.h:69
ExpressReco_HICollisions_FallBack.trackCollection
tuple trackCollection
Definition: ExpressReco_HICollisions_FallBack.py:7877
debug_cff.d0
tuple d0
Definition: debug_cff.py:3
reco::TrackCollection
std::vector< Track > TrackCollection
collection of Tracks
Definition: TrackFwd.h:10
edm::HandleBase::isValid
bool isValid() const
Definition: HandleBase.h:76
edm::Event::getByLabel
bool getByLabel(InputTag const &tag, Handle< PROD > &result) const
Definition: Event.h:359
testEve_cfg.tracks
tuple tracks
Definition: testEve_cfg.py:39
math::XYZPoint
XYZPointD XYZPoint
point in space with cartesian internal representation
Definition: Point3D.h:13
edm::Handle::product
T const * product() const
Definition: Handle.h:74
reco::BeamSpot::y0
double y0() const
y coordinate
Definition: BeamSpot.h:67
reco::BeamSpot::x0
double x0() const
x coordinate
Definition: BeamSpot.h:65
void PhotonIsolationCalculator::classify ( const reco::Photon *  photon,
bool &  isEBPho,
bool &  isEEPho,
bool &  isEBEtaGap,
bool &  isEBPhiGap,
bool &  isEERingGap,
bool &  isEEDeeGap,
bool &  isEBEEGap 
)

Definition at line 431 of file PhotonIsolationCalculator.cc.

References ExpressReco_HICollisions_FallBack::detector, EcalBarrel, EcalEndcap, eta(), reco::CaloCluster::hitsAndFractions(), EEDetId::isNextToDBoundary(), EBDetId::isNextToEtaBoundary(), EBDetId::isNextToPhiBoundary(), EEDetId::isNextToRingBoundary(), DetId::subdetId(), reco::Photon::superCluster(), and funct::true.

Referenced by calculate().

438  {
439 
440 
441  const reco::CaloCluster & seedCluster = *(photon->superCluster()->seed()) ;
442  // following line is temporary until the D. Evans or F. Ferri submit the new tag for ClusterAlgos
443  // DEfinitive will be
444  // DetId seedXtalId = scRef->seed()->seed();
445  DetId seedXtalId = seedCluster.hitsAndFractions()[0].first ;
446  int detector = seedXtalId.subdetId() ;
447 
448  //Set fiducial flags for this photon.
449  double eta = photon->superCluster()->position().eta();
450  double feta = fabs(eta);
451 
452  if (detector==EcalBarrel) {
453 
454 
455  isEBPho = true;
456  if (EBDetId::isNextToEtaBoundary(EBDetId(seedXtalId))) {
457  if (fabs(feta-1.479)<.1) isEBEEGap = true ;
458  else
459  isEBEtaGap = true ;
460  }
461 
462  if (EBDetId::isNextToPhiBoundary(EBDetId(seedXtalId)))
463  isEBPhiGap = true ;
464 
465 
466 
467  } else if (detector==EcalEndcap) {
468 
469  isEEPho = true;
470  if (EEDetId::isNextToRingBoundary(EEDetId(seedXtalId))) {
471  if (fabs(feta-1.479)<.1) isEBEEGap = true ;
472  else
473  isEERingGap = true ;
474  }
475 
476  if (EEDetId::isNextToDBoundary(EEDetId(seedXtalId)))
477  isEEDeeGap = true ;
478  }
479 
480 
481 }
EBDetId::isNextToEtaBoundary
static bool isNextToEtaBoundary(EBDetId id)
Definition: EBDetId.cc:175
reco::Photon::superCluster
reco::SuperClusterRef superCluster() const
Retrieve photonCore attributes.
Definition: Photon.h:54
reco::CaloCluster::hitsAndFractions
const std::vector< std::pair< DetId, float > > & hitsAndFractions() const
Definition: CaloCluster.h:178
eta
T eta() const
Definition: newBasic3DVector.h:122
EBDetId::isNextToPhiBoundary
static bool isNextToPhiBoundary(EBDetId id)
Definition: EBDetId.cc:180
EEDetId::isNextToRingBoundary
static bool isNextToRingBoundary(EEDetId id)
Definition: EEDetId.cc:391
EEDetId::isNextToDBoundary
static bool isNextToDBoundary(EEDetId id)
Definition: EEDetId.cc:383
DetId::subdetId
int subdetId() const
get the contents of the subdetector field (not cast into any detector&#39;s numbering enum) ...
Definition: DetId.h:39
DetId
Definition: DetId.h:20
ExpressReco_HICollisions_FallBack.detector
tuple detector
Definition: ExpressReco_HICollisions_FallBack.py:25265
void PhotonIsolationCalculator::setup ( const edm::ParameterSet &  conf)

Isolation parameters for barrel and for two different cone sizes

Isolation parameters for Endcap and for two different cone sizes

Definition at line 46 of file PhotonIsolationCalculator.cc.

References barrelecalCollection_, barrelecalProducer_, beamSpotProducerTag_, ecalIsoBarrelRadiusA_, ecalIsoBarrelRadiusB_, ecalIsoEndcapRadiusA_, ecalIsoEndcapRadiusB_, endcapecalCollection_, endcapecalProducer_, edm::ParameterSet::getParameter(), hcalCollection_, hcalIsoBarrelRadiusA_, hcalIsoBarrelRadiusB_, hcalIsoEndcapRadiusA_, hcalIsoEndcapRadiusB_, hcalProducer_, EcalSeverityLevelAlgo::kE1OverE9, EcalSeverityLevelAlgo::kSwissCross, EcalSeverityLevelAlgo::kSwissCrossBordersIncluded, moduleEtaBoundary_, modulePhiBoundary_, severityLevelCut_, severityRecHitThreshold_, python.multivaluedict::sort(), spId_, spikeIdThreshold_, trackInputTag_, trkIsoBarrelRadiusA_, trkIsoBarrelRadiusB_, trkIsoEndcapRadiusA_, trkIsoEndcapRadiusB_, useNumCrystals_, v_chstatus_, and vetoClusteredEcalHits_.

Referenced by PhotonProducer::beginRun().

46  {
47 
48 
49  trackInputTag_ = conf.getParameter<edm::InputTag>("trackProducer");
50  beamSpotProducerTag_ = conf.getParameter<edm::InputTag>("beamSpotProducer");
51  barrelecalCollection_ = conf.getParameter<std::string>("barrelEcalRecHitCollection");
52  barrelecalProducer_ = conf.getParameter<std::string>("barrelEcalRecHitProducer");
53  endcapecalCollection_ = conf.getParameter<std::string>("endcapEcalRecHitCollection");
54  endcapecalProducer_ = conf.getParameter<std::string>("endcapEcalRecHitProducer");
55  hcalCollection_ = conf.getParameter<std::string>("HcalRecHitCollection");
56  hcalProducer_ = conf.getParameter<std::string>("HcalRecHitProducer");
57 
58  // gsfRecoInputTag_ = conf.getParameter<edm::InputTag>("GsfRecoCollection");
59  modulePhiBoundary_ = conf.getParameter<double>("modulePhiBoundary");
60  moduleEtaBoundary_ = conf.getParameter<std::vector<double> >("moduleEtaBoundary");
61  //
62  vetoClusteredEcalHits_ = conf.getParameter<bool>("vetoClustered");
63  useNumCrystals_ = conf.getParameter<bool>("useNumCrystals");
65  trkIsoBarrelRadiusA_.push_back( conf.getParameter<double>("TrackConeOuterRadiusA_Barrel") );
66  trkIsoBarrelRadiusA_.push_back( conf.getParameter<double>("TrackConeInnerRadiusA_Barrel") ) ;
67  trkIsoBarrelRadiusA_.push_back( conf.getParameter<double>("isolationtrackThresholdA_Barrel") );
68  trkIsoBarrelRadiusA_.push_back( conf.getParameter<double>("longImpactParameterA_Barrel") );
69  trkIsoBarrelRadiusA_.push_back( conf.getParameter<double>("transImpactParameterA_Barrel") );
70  trkIsoBarrelRadiusA_.push_back( conf.getParameter<double>("isolationtrackEtaSliceA_Barrel") );
71 
72  ecalIsoBarrelRadiusA_.push_back( conf.getParameter<double>("EcalRecHitInnerRadiusA_Barrel") );
73  ecalIsoBarrelRadiusA_.push_back( conf.getParameter<double>("EcalRecHitOuterRadiusA_Barrel") );
74  ecalIsoBarrelRadiusA_.push_back( conf.getParameter<double>("EcalRecHitEtaSliceA_Barrel") );
75  ecalIsoBarrelRadiusA_.push_back( conf.getParameter<double>("EcalRecHitThreshEA_Barrel") );
76  ecalIsoBarrelRadiusA_.push_back( conf.getParameter<double>("EcalRecHitThreshEtA_Barrel") );
77 
78  hcalIsoBarrelRadiusA_.push_back( conf.getParameter<double>("HcalTowerInnerRadiusA_Barrel") );
79  hcalIsoBarrelRadiusA_.push_back( conf.getParameter<double>("HcalTowerOuterRadiusA_Barrel") );
80  hcalIsoBarrelRadiusA_.push_back( conf.getParameter<double>("HcalTowerThreshEA_Barrel") );
81  hcalIsoBarrelRadiusA_.push_back( conf.getParameter<double>("HcalDepth1TowerInnerRadiusA_Barrel") );
82  hcalIsoBarrelRadiusA_.push_back( conf.getParameter<double>("HcalDepth1TowerOuterRadiusA_Barrel") );
83  hcalIsoBarrelRadiusA_.push_back( conf.getParameter<double>("HcalDepth1TowerThreshEA_Barrel") );
84  hcalIsoBarrelRadiusA_.push_back( conf.getParameter<double>("HcalDepth2TowerInnerRadiusA_Barrel") );
85  hcalIsoBarrelRadiusA_.push_back( conf.getParameter<double>("HcalDepth2TowerOuterRadiusA_Barrel") );
86  hcalIsoBarrelRadiusA_.push_back( conf.getParameter<double>("HcalDepth2TowerThreshEA_Barrel") );
87 
88 
89  trkIsoBarrelRadiusB_.push_back( conf.getParameter<double>("TrackConeOuterRadiusB_Barrel") );
90  trkIsoBarrelRadiusB_.push_back( conf.getParameter<double>("TrackConeInnerRadiusB_Barrel") ) ;
91  trkIsoBarrelRadiusB_.push_back( conf.getParameter<double>("isolationtrackThresholdB_Barrel") );
92  trkIsoBarrelRadiusB_.push_back( conf.getParameter<double>("longImpactParameterB_Barrel") );
93  trkIsoBarrelRadiusB_.push_back( conf.getParameter<double>("transImpactParameterB_Barrel") );
94  trkIsoBarrelRadiusB_.push_back( conf.getParameter<double>("isolationtrackEtaSliceB_Barrel") );
95 
96  ecalIsoBarrelRadiusB_.push_back( conf.getParameter<double>("EcalRecHitInnerRadiusB_Barrel") );
97  ecalIsoBarrelRadiusB_.push_back( conf.getParameter<double>("EcalRecHitOuterRadiusB_Barrel") );
98  ecalIsoBarrelRadiusB_.push_back( conf.getParameter<double>("EcalRecHitEtaSliceB_Barrel") );
99  ecalIsoBarrelRadiusB_.push_back( conf.getParameter<double>("EcalRecHitThreshEB_Barrel") );
100  ecalIsoBarrelRadiusB_.push_back( conf.getParameter<double>("EcalRecHitThreshEtB_Barrel") );
101 
102  hcalIsoBarrelRadiusB_.push_back( conf.getParameter<double>("HcalTowerInnerRadiusB_Barrel") );
103  hcalIsoBarrelRadiusB_.push_back( conf.getParameter<double>("HcalTowerOuterRadiusB_Barrel") );
104  hcalIsoBarrelRadiusB_.push_back( conf.getParameter<double>("HcalTowerThreshEB_Barrel") );
105  hcalIsoBarrelRadiusB_.push_back( conf.getParameter<double>("HcalDepth1TowerInnerRadiusB_Barrel") );
106  hcalIsoBarrelRadiusB_.push_back( conf.getParameter<double>("HcalDepth1TowerOuterRadiusB_Barrel") );
107  hcalIsoBarrelRadiusB_.push_back( conf.getParameter<double>("HcalDepth1TowerThreshEB_Barrel") );
108  hcalIsoBarrelRadiusB_.push_back( conf.getParameter<double>("HcalDepth2TowerInnerRadiusB_Barrel") );
109  hcalIsoBarrelRadiusB_.push_back( conf.getParameter<double>("HcalDepth2TowerOuterRadiusB_Barrel") );
110  hcalIsoBarrelRadiusB_.push_back( conf.getParameter<double>("HcalDepth2TowerThreshEB_Barrel") );
111 
113  trkIsoEndcapRadiusA_.push_back( conf.getParameter<double>("TrackConeOuterRadiusA_Endcap") );
114  trkIsoEndcapRadiusA_.push_back( conf.getParameter<double>("TrackConeInnerRadiusA_Endcap") ) ;
115  trkIsoEndcapRadiusA_.push_back( conf.getParameter<double>("isolationtrackThresholdA_Endcap") );
116  trkIsoEndcapRadiusA_.push_back( conf.getParameter<double>("longImpactParameterA_Endcap") );
117  trkIsoEndcapRadiusA_.push_back( conf.getParameter<double>("transImpactParameterA_Endcap") );
118  trkIsoEndcapRadiusA_.push_back( conf.getParameter<double>("isolationtrackEtaSliceA_Endcap") );
119 
120 
121  ecalIsoEndcapRadiusA_.push_back( conf.getParameter<double>("EcalRecHitInnerRadiusA_Endcap") );
122  ecalIsoEndcapRadiusA_.push_back( conf.getParameter<double>("EcalRecHitOuterRadiusA_Endcap") );
123  ecalIsoEndcapRadiusA_.push_back( conf.getParameter<double>("EcalRecHitEtaSliceA_Endcap") );
124  ecalIsoEndcapRadiusA_.push_back( conf.getParameter<double>("EcalRecHitThreshEA_Endcap") );
125  ecalIsoEndcapRadiusA_.push_back( conf.getParameter<double>("EcalRecHitThreshEtA_Endcap") );
126 
127 
128  hcalIsoEndcapRadiusA_.push_back( conf.getParameter<double>("HcalTowerInnerRadiusA_Endcap") );
129  hcalIsoEndcapRadiusA_.push_back( conf.getParameter<double>("HcalTowerOuterRadiusA_Endcap") );
130  hcalIsoEndcapRadiusA_.push_back( conf.getParameter<double>("HcalTowerThreshEA_Endcap") );
131  hcalIsoEndcapRadiusA_.push_back( conf.getParameter<double>("HcalDepth1TowerInnerRadiusA_Endcap") );
132  hcalIsoEndcapRadiusA_.push_back( conf.getParameter<double>("HcalDepth1TowerOuterRadiusA_Endcap") );
133  hcalIsoEndcapRadiusA_.push_back( conf.getParameter<double>("HcalDepth1TowerThreshEA_Endcap") );
134  hcalIsoEndcapRadiusA_.push_back( conf.getParameter<double>("HcalDepth2TowerInnerRadiusA_Endcap") );
135  hcalIsoEndcapRadiusA_.push_back( conf.getParameter<double>("HcalDepth2TowerOuterRadiusA_Endcap") );
136  hcalIsoEndcapRadiusA_.push_back( conf.getParameter<double>("HcalDepth2TowerThreshEA_Endcap") );
137 
138 
139  trkIsoEndcapRadiusB_.push_back( conf.getParameter<double>("TrackConeOuterRadiusB_Endcap") );
140  trkIsoEndcapRadiusB_.push_back( conf.getParameter<double>("TrackConeInnerRadiusB_Endcap") ) ;
141  trkIsoEndcapRadiusB_.push_back( conf.getParameter<double>("isolationtrackThresholdB_Endcap") );
142  trkIsoEndcapRadiusB_.push_back( conf.getParameter<double>("longImpactParameterB_Endcap") );
143  trkIsoEndcapRadiusB_.push_back( conf.getParameter<double>("transImpactParameterB_Endcap") );
144  trkIsoEndcapRadiusB_.push_back( conf.getParameter<double>("isolationtrackEtaSliceB_Endcap") );
145 
146  ecalIsoEndcapRadiusB_.push_back( conf.getParameter<double>("EcalRecHitInnerRadiusB_Endcap") );
147  ecalIsoEndcapRadiusB_.push_back( conf.getParameter<double>("EcalRecHitOuterRadiusB_Endcap") );
148  ecalIsoEndcapRadiusB_.push_back( conf.getParameter<double>("EcalRecHitEtaSliceB_Endcap") );
149  ecalIsoEndcapRadiusB_.push_back( conf.getParameter<double>("EcalRecHitThreshEB_Endcap") );
150  ecalIsoEndcapRadiusB_.push_back( conf.getParameter<double>("EcalRecHitThreshEtB_Endcap") );
151 
152  hcalIsoEndcapRadiusB_.push_back( conf.getParameter<double>("HcalTowerInnerRadiusB_Endcap") );
153  hcalIsoEndcapRadiusB_.push_back( conf.getParameter<double>("HcalTowerOuterRadiusB_Endcap") );
154  hcalIsoEndcapRadiusB_.push_back( conf.getParameter<double>("HcalTowerThreshEB_Endcap") );
155  hcalIsoEndcapRadiusB_.push_back( conf.getParameter<double>("HcalDepth1TowerInnerRadiusB_Endcap") );
156  hcalIsoEndcapRadiusB_.push_back( conf.getParameter<double>("HcalDepth1TowerOuterRadiusB_Endcap") );
157  hcalIsoEndcapRadiusB_.push_back( conf.getParameter<double>("HcalDepth1TowerThreshEB_Endcap") );
158  hcalIsoEndcapRadiusB_.push_back( conf.getParameter<double>("HcalDepth2TowerInnerRadiusB_Endcap") );
159  hcalIsoEndcapRadiusB_.push_back( conf.getParameter<double>("HcalDepth2TowerOuterRadiusB_Endcap") );
160  hcalIsoEndcapRadiusB_.push_back( conf.getParameter<double>("HcalDepth2TowerThreshEB_Endcap") );
161 
162  //Pick up the variables for the spike removal
163  severityLevelCut_ = conf.getParameter<int>("severityLevelCut");
164  severityRecHitThreshold_ = conf.getParameter<double>("severityRecHitThreshold");
165  spikeIdThreshold_ = conf.getParameter<double>("spikeIdThreshold");
166 
167  const std::vector<int> &tempV(conf.getParameter<std::vector<int> >("recHitFlagsToBeExcluded"));
168  v_chstatus_.clear();
169  v_chstatus_.insert(v_chstatus_.begin(),tempV.begin(),tempV.end());
170  std::sort( v_chstatus_.begin(), v_chstatus_.end() );
171 
172 
173  //Need to figure out which algo to use
174  if(!conf.getParameter<std::string>("spikeIdString").compare("kE1OverE9") ) {
175  spId_ = EcalSeverityLevelAlgo::kE1OverE9;
176  } else if(!conf.getParameter<std::string>("spikeIdString").compare("kSwissCross") ) {
177  spId_ = EcalSeverityLevelAlgo::kSwissCross;
178  } else if(!conf.getParameter<std::string>("spikeIdString").compare("kSwissCrossBordersIncluded") ) {
179  spId_ = EcalSeverityLevelAlgo::kSwissCrossBordersIncluded;
180  } else {
181  spId_ = EcalSeverityLevelAlgo::kSwissCrossBordersIncluded;
182  edm::LogWarning("PhotonIsolationCalculator|SpikeRemovalForIsolation")
183  << "Cannot find the requested method. kSwissCross set instead.";
184  }
185 }
edm::ParameterSet::getParameter
T getParameter(std::string const &) const
PhotonIsolationCalculator::ecalIsoBarrelRadiusA_
std::vector< double > ecalIsoBarrelRadiusA_
Definition: PhotonIsolationCalculator.h:91
PhotonIsolationCalculator::hcalIsoBarrelRadiusB_
std::vector< double > hcalIsoBarrelRadiusB_
Definition: PhotonIsolationCalculator.h:95
PhotonIsolationCalculator::ecalIsoEndcapRadiusB_
std::vector< double > ecalIsoEndcapRadiusB_
Definition: PhotonIsolationCalculator.h:101
PhotonIsolationCalculator::trkIsoEndcapRadiusB_
std::vector< double > trkIsoEndcapRadiusB_
Definition: PhotonIsolationCalculator.h:100
PhotonIsolationCalculator::hcalIsoEndcapRadiusA_
std::vector< double > hcalIsoEndcapRadiusA_
Definition: PhotonIsolationCalculator.h:99
PhotonIsolationCalculator::trkIsoBarrelRadiusB_
std::vector< double > trkIsoBarrelRadiusB_
Definition: PhotonIsolationCalculator.h:93
PhotonIsolationCalculator::hcalIsoBarrelRadiusA_
std::vector< double > hcalIsoBarrelRadiusA_
Definition: PhotonIsolationCalculator.h:92
PhotonIsolationCalculator::spId_
EcalSeverityLevelAlgo::SpikeId spId_
Definition: PhotonIsolationCalculator.h:129
PhotonIsolationCalculator::trkIsoEndcapRadiusA_
std::vector< double > trkIsoEndcapRadiusA_
Definition: PhotonIsolationCalculator.h:97
PhotonIsolationCalculator::moduleEtaBoundary_
std::vector< double > moduleEtaBoundary_
Definition: PhotonIsolationCalculator.h:86
PhotonIsolationCalculator::hcalIsoEndcapRadiusB_
std::vector< double > hcalIsoEndcapRadiusB_
Definition: PhotonIsolationCalculator.h:102
PhotonIsolationCalculator::ecalIsoEndcapRadiusA_
std::vector< double > ecalIsoEndcapRadiusA_
Definition: PhotonIsolationCalculator.h:98
PhotonIsolationCalculator::ecalIsoBarrelRadiusB_
std::vector< double > ecalIsoBarrelRadiusB_
Definition: PhotonIsolationCalculator.h:94
PhotonIsolationCalculator::trkIsoBarrelRadiusA_
std::vector< double > trkIsoBarrelRadiusA_
Definition: PhotonIsolationCalculator.h:90

Member Data Documentation

std::string PhotonIsolationCalculator::barrelecalCollection_
protected

Definition at line 75 of file PhotonIsolationCalculator.h.

Referenced by calculateEcalRecHitIso(), and setup().

std::string PhotonIsolationCalculator::barrelecalProducer_
protected

Definition at line 76 of file PhotonIsolationCalculator.h.

Referenced by calculateEcalRecHitIso(), and setup().

edm::InputTag PhotonIsolationCalculator::beamSpotProducerTag_
protected

Definition at line 84 of file PhotonIsolationCalculator.h.

Referenced by calculateTrackIso(), and setup().

std::vector<double> PhotonIsolationCalculator::ecalIsoBarrelRadiusA_
protected

Definition at line 91 of file PhotonIsolationCalculator.h.

Referenced by calculate(), and setup().

std::vector<double> PhotonIsolationCalculator::ecalIsoBarrelRadiusB_
protected

Definition at line 94 of file PhotonIsolationCalculator.h.

Referenced by calculate(), and setup().

std::vector<double> PhotonIsolationCalculator::ecalIsoEndcapRadiusA_
protected

Definition at line 98 of file PhotonIsolationCalculator.h.

Referenced by calculate(), and setup().

std::vector<double> PhotonIsolationCalculator::ecalIsoEndcapRadiusB_
protected

Definition at line 101 of file PhotonIsolationCalculator.h.

Referenced by calculate(), and setup().

std::string PhotonIsolationCalculator::endcapecalCollection_
protected

Definition at line 77 of file PhotonIsolationCalculator.h.

Referenced by calculateEcalRecHitIso(), and setup().

std::string PhotonIsolationCalculator::endcapecalProducer_
protected

Definition at line 78 of file PhotonIsolationCalculator.h.

Referenced by calculateEcalRecHitIso(), and setup().

std::string PhotonIsolationCalculator::hcalCollection_
protected

Definition at line 79 of file PhotonIsolationCalculator.h.

Referenced by calculateHcalTowerIso(), and setup().

std::vector<double> PhotonIsolationCalculator::hcalIsoBarrelRadiusA_
protected

Definition at line 92 of file PhotonIsolationCalculator.h.

Referenced by calculate(), and setup().

std::vector<double> PhotonIsolationCalculator::hcalIsoBarrelRadiusB_
protected

Definition at line 95 of file PhotonIsolationCalculator.h.

Referenced by calculate(), and setup().

std::vector<double> PhotonIsolationCalculator::hcalIsoEndcapRadiusA_
protected

Definition at line 99 of file PhotonIsolationCalculator.h.

Referenced by calculate(), and setup().

std::vector<double> PhotonIsolationCalculator::hcalIsoEndcapRadiusB_
protected

Definition at line 102 of file PhotonIsolationCalculator.h.

Referenced by calculate(), and setup().

std::string PhotonIsolationCalculator::hcalProducer_
protected

Definition at line 80 of file PhotonIsolationCalculator.h.

Referenced by calculateHcalTowerIso(), and setup().

double PhotonIsolationCalculator::isolationtrackEtaSliceA_
protected

Definition at line 122 of file PhotonIsolationCalculator.h.

Referenced by calculate().

double PhotonIsolationCalculator::isolationtrackEtaSliceB_
protected

Definition at line 149 of file PhotonIsolationCalculator.h.

Referenced by calculate().

double PhotonIsolationCalculator::isolationtrackThresholdA_
protected

Definition at line 121 of file PhotonIsolationCalculator.h.

Referenced by calculate().

double PhotonIsolationCalculator::isolationtrackThresholdB_
protected

Definition at line 148 of file PhotonIsolationCalculator.h.

Referenced by calculate().

std::vector<double> PhotonIsolationCalculator::moduleEtaBoundary_
protected

Definition at line 86 of file PhotonIsolationCalculator.h.

Referenced by setup().

double PhotonIsolationCalculator::modulePhiBoundary_
protected

Definition at line 85 of file PhotonIsolationCalculator.h.

Referenced by setup().

double PhotonIsolationCalculator::photonEcalRecHitConeInnerRadiusA_
protected

Definition at line 105 of file PhotonIsolationCalculator.h.

Referenced by calculate().

double PhotonIsolationCalculator::photonEcalRecHitConeInnerRadiusB_
protected

Definition at line 132 of file PhotonIsolationCalculator.h.

Referenced by calculate().

double PhotonIsolationCalculator::photonEcalRecHitConeOuterRadiusA_
protected

Definition at line 106 of file PhotonIsolationCalculator.h.

Referenced by calculate().

double PhotonIsolationCalculator::photonEcalRecHitConeOuterRadiusB_
protected

Definition at line 133 of file PhotonIsolationCalculator.h.

Referenced by calculate().

double PhotonIsolationCalculator::photonEcalRecHitEtaSliceA_
protected

Definition at line 107 of file PhotonIsolationCalculator.h.

Referenced by calculate().

double PhotonIsolationCalculator::photonEcalRecHitEtaSliceB_
protected

Definition at line 134 of file PhotonIsolationCalculator.h.

Referenced by calculate().

double PhotonIsolationCalculator::photonEcalRecHitThreshEA_
protected

Definition at line 108 of file PhotonIsolationCalculator.h.

Referenced by calculate().

double PhotonIsolationCalculator::photonEcalRecHitThreshEB_
protected

Definition at line 135 of file PhotonIsolationCalculator.h.

Referenced by calculate().

double PhotonIsolationCalculator::photonEcalRecHitThreshEtA_
protected

Definition at line 109 of file PhotonIsolationCalculator.h.

Referenced by calculate().

double PhotonIsolationCalculator::photonEcalRecHitThreshEtB_
protected

Definition at line 136 of file PhotonIsolationCalculator.h.

Referenced by calculate().

double PhotonIsolationCalculator::photonHcalDepth1TowerConeInnerRadiusA_
protected

Definition at line 113 of file PhotonIsolationCalculator.h.

Referenced by calculate().

double PhotonIsolationCalculator::photonHcalDepth1TowerConeInnerRadiusB_
protected

Definition at line 140 of file PhotonIsolationCalculator.h.

Referenced by calculate().

double PhotonIsolationCalculator::photonHcalDepth1TowerConeOuterRadiusA_
protected

Definition at line 114 of file PhotonIsolationCalculator.h.

Referenced by calculate().

double PhotonIsolationCalculator::photonHcalDepth1TowerConeOuterRadiusB_
protected

Definition at line 141 of file PhotonIsolationCalculator.h.

Referenced by calculate().

double PhotonIsolationCalculator::photonHcalDepth1TowerThreshEA_
protected

Definition at line 115 of file PhotonIsolationCalculator.h.

Referenced by calculate().

double PhotonIsolationCalculator::photonHcalDepth1TowerThreshEB_
protected

Definition at line 142 of file PhotonIsolationCalculator.h.

Referenced by calculate().

double PhotonIsolationCalculator::photonHcalDepth2TowerConeInnerRadiusA_
protected

Definition at line 116 of file PhotonIsolationCalculator.h.

Referenced by calculate().

double PhotonIsolationCalculator::photonHcalDepth2TowerConeInnerRadiusB_
protected

Definition at line 143 of file PhotonIsolationCalculator.h.

Referenced by calculate().

double PhotonIsolationCalculator::photonHcalDepth2TowerConeOuterRadiusA_
protected

Definition at line 117 of file PhotonIsolationCalculator.h.

Referenced by calculate().

double PhotonIsolationCalculator::photonHcalDepth2TowerConeOuterRadiusB_
protected

Definition at line 144 of file PhotonIsolationCalculator.h.

Referenced by calculate().

double PhotonIsolationCalculator::photonHcalDepth2TowerThreshEA_
protected

Definition at line 118 of file PhotonIsolationCalculator.h.

Referenced by calculate().

double PhotonIsolationCalculator::photonHcalDepth2TowerThreshEB_
protected

Definition at line 145 of file PhotonIsolationCalculator.h.

Referenced by calculate().

double PhotonIsolationCalculator::photonHcalTowerConeInnerRadiusA_
protected

Definition at line 110 of file PhotonIsolationCalculator.h.

Referenced by calculate().

double PhotonIsolationCalculator::photonHcalTowerConeInnerRadiusB_
protected

Definition at line 137 of file PhotonIsolationCalculator.h.

Referenced by calculate().

double PhotonIsolationCalculator::photonHcalTowerConeOuterRadiusA_
protected

Definition at line 111 of file PhotonIsolationCalculator.h.

Referenced by calculate().

double PhotonIsolationCalculator::photonHcalTowerConeOuterRadiusB_
protected

Definition at line 138 of file PhotonIsolationCalculator.h.

Referenced by calculate().

double PhotonIsolationCalculator::photonHcalTowerThreshEA_
protected

Definition at line 112 of file PhotonIsolationCalculator.h.

Referenced by calculate().

double PhotonIsolationCalculator::photonHcalTowerThreshEB_
protected

Definition at line 139 of file PhotonIsolationCalculator.h.

Referenced by calculate().

int PhotonIsolationCalculator::severityLevelCut_
protected

Definition at line 126 of file PhotonIsolationCalculator.h.

Referenced by calculateEcalRecHitIso(), and setup().

float PhotonIsolationCalculator::severityRecHitThreshold_
protected

Definition at line 127 of file PhotonIsolationCalculator.h.

Referenced by calculateEcalRecHitIso(), and setup().

EcalSeverityLevelAlgo::SpikeId PhotonIsolationCalculator::spId_
protected

Definition at line 129 of file PhotonIsolationCalculator.h.

Referenced by calculateEcalRecHitIso(), and setup().

float PhotonIsolationCalculator::spikeIdThreshold_
protected

Definition at line 128 of file PhotonIsolationCalculator.h.

Referenced by calculateEcalRecHitIso(), and setup().

double PhotonIsolationCalculator::trackConeInnerRadiusA_
protected

Definition at line 120 of file PhotonIsolationCalculator.h.

Referenced by calculate().

double PhotonIsolationCalculator::trackConeInnerRadiusB_
protected

Definition at line 147 of file PhotonIsolationCalculator.h.

Referenced by calculate().

double PhotonIsolationCalculator::trackConeOuterRadiusA_
protected

Definition at line 119 of file PhotonIsolationCalculator.h.

Referenced by calculate().

double PhotonIsolationCalculator::trackConeOuterRadiusB_
protected

Definition at line 146 of file PhotonIsolationCalculator.h.

Referenced by calculate().

double PhotonIsolationCalculator::trackD0RadiusA_
protected

Definition at line 124 of file PhotonIsolationCalculator.h.

Referenced by calculate().

double PhotonIsolationCalculator::trackD0RadiusB_
protected

Definition at line 151 of file PhotonIsolationCalculator.h.

Referenced by calculate().

edm::InputTag PhotonIsolationCalculator::trackInputTag_
protected

Definition at line 82 of file PhotonIsolationCalculator.h.

Referenced by calculateTrackIso(), and setup().

double PhotonIsolationCalculator::trackLipRadiusA_
protected

Definition at line 123 of file PhotonIsolationCalculator.h.

Referenced by calculate().

double PhotonIsolationCalculator::trackLipRadiusB_
protected

Definition at line 150 of file PhotonIsolationCalculator.h.

Referenced by calculate().

std::vector<double> PhotonIsolationCalculator::trkIsoBarrelRadiusA_
protected

Definition at line 90 of file PhotonIsolationCalculator.h.

Referenced by calculate(), and setup().

std::vector<double> PhotonIsolationCalculator::trkIsoBarrelRadiusB_
protected

Definition at line 93 of file PhotonIsolationCalculator.h.

Referenced by calculate(), and setup().

std::vector<double> PhotonIsolationCalculator::trkIsoEndcapRadiusA_
protected

Definition at line 97 of file PhotonIsolationCalculator.h.

Referenced by calculate(), and setup().

std::vector<double> PhotonIsolationCalculator::trkIsoEndcapRadiusB_
protected

Definition at line 100 of file PhotonIsolationCalculator.h.

Referenced by calculate(), and setup().

bool PhotonIsolationCalculator::useNumCrystals_
protected

Definition at line 88 of file PhotonIsolationCalculator.h.

Referenced by calculate(), and setup().

std::vector<int> PhotonIsolationCalculator::v_chstatus_
protected

Definition at line 130 of file PhotonIsolationCalculator.h.

Referenced by setup().

bool PhotonIsolationCalculator::vetoClusteredEcalHits_
protected

Definition at line 87 of file PhotonIsolationCalculator.h.

Referenced by calculate(), and setup().

Generated for CMSSW Reference Manual by doxygen 1.8.5