d4/d30/PhotonIsolationCalculator_8cc_source.html

 #include "RecoEgamma/PhotonIdentification/interface/PhotonIsolationCalculator.h"

 #include "CommonTools/Utils/interface/StringToEnumValue.h"

 #include "DataFormats/EgammaReco/interface/BasicCluster.h"

 #include "DataFormats/EgammaReco/interface/BasicClusterFwd.h"

 #include "DataFormats/EgammaReco/interface/SuperCluster.h"

 #include "DataFormats/TrackReco/interface/Track.h"

 #include "DataFormats/TrackReco/interface/TrackFwd.h"

 #include "DataFormats/EgammaCandidates/interface/GsfElectron.h"

 #include "DataFormats/EgammaCandidates/interface/GsfElectronFwd.h"

 #include "DataFormats/DetId/interface/DetId.h"

 #include "DataFormats/EcalDetId/interface/EEDetId.h"

 #include "DataFormats/EcalDetId/interface/EBDetId.h"

 #include "DataFormats/EcalRecHit/interface/EcalRecHit.h"

 #include "DataFormats/EcalRecHit/interface/EcalRecHitCollections.h"

 #include "DataFormats/HcalRecHit/interface/HcalRecHitCollections.h"

 #include "DataFormats/BeamSpot/interface/BeamSpot.h"


 #include "RecoEgamma/EgammaIsolationAlgos/interface/EgammaHcalIsolation.h"

 #include "RecoEgamma/EgammaIsolationAlgos/interface/EgammaRecHitIsolation.h"

 #include "RecoEgamma/EgammaIsolationAlgos/interface/PhotonTkIsolation.h"

 #include "RecoEgamma/EgammaIsolationAlgos/interface/EgammaEcalIsolation.h"

 #include "RecoEgamma/EgammaIsolationAlgos/interface/EgammaTowerIsolation.h"

 #include "RecoCaloTools/MetaCollections/interface/CaloRecHitMetaCollections.h"

 //#include "RecoEcal/EgammaCoreTools/interface/EcalClusterTools.h"


 #include "Geometry/CaloTopology/interface/CaloTopology.h"

 #include "Geometry/CaloEventSetup/interface/CaloTopologyRecord.h"

 #include "Geometry/CaloGeometry/interface/CaloGeometry.h"

 #include "Geometry/CaloTopology/interface/EcalBarrelTopology.h"

 #include "Geometry/CaloTopology/interface/EcalEndcapTopology.h"

 #include "Geometry/Records/interface/CaloGeometryRecord.h"


 #include "CondFormats/DataRecord/interface/EcalChannelStatusRcd.h"

 #include "RecoLocalCalo/EcalRecAlgos/interface/EcalSeverityLevelAlgoRcd.h"


 #include <string>

 #include <TMath.h>


 void PhotonIsolationCalculator::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) {


   trackInputTag_ = conf.getParameter<edm::InputTag>("trackProducer");

   beamSpotProducerTag_ = conf.getParameter<edm::InputTag>("beamSpotProducer");

   barrelecalCollection_ = conf.getParameter<edm::InputTag>("barrelEcalRecHitCollection");

   endcapecalCollection_ = conf.getParameter<edm::InputTag>("endcapEcalRecHitCollection");

   hcalCollection_ = conf.getParameter<edm::InputTag>("HcalRecHitCollection");


   //  gsfRecoInputTag_ = conf.getParameter<edm::InputTag>("GsfRecoCollection");

   modulePhiBoundary_ = conf.getParameter<double>("modulePhiBoundary");

   moduleEtaBoundary_ = conf.getParameter<std::vector<double> >("moduleEtaBoundary");

   //

   vetoClusteredEcalHits_ = conf.getParameter<bool>("vetoClustered");

   useNumCrystals_ = conf.getParameter<bool>("useNumCrystals");


   int i=0;

   trkIsoBarrelRadiusA_[i++] = (  conf.getParameter<double>("TrackConeOuterRadiusA_Barrel") );

   trkIsoBarrelRadiusA_[i++] = (  conf.getParameter<double>("TrackConeInnerRadiusA_Barrel") ) ;

   trkIsoBarrelRadiusA_[i++] = (  conf.getParameter<double>("isolationtrackThresholdA_Barrel") );

   trkIsoBarrelRadiusA_[i++] = (  conf.getParameter<double>("longImpactParameterA_Barrel") );

   trkIsoBarrelRadiusA_[i++] = (  conf.getParameter<double>("transImpactParameterA_Barrel") );

   trkIsoBarrelRadiusA_[i++] = (  conf.getParameter<double>("isolationtrackEtaSliceA_Barrel") );


   i=0;

   ecalIsoBarrelRadiusA_[i++] = ( conf.getParameter<double>("EcalRecHitInnerRadiusA_Barrel") );

   ecalIsoBarrelRadiusA_[i++] = ( conf.getParameter<double>("EcalRecHitOuterRadiusA_Barrel") );

   ecalIsoBarrelRadiusA_[i++] = ( conf.getParameter<double>("EcalRecHitEtaSliceA_Barrel") );

   ecalIsoBarrelRadiusA_[i++] = ( conf.getParameter<double>("EcalRecHitThreshEA_Barrel") );

   ecalIsoBarrelRadiusA_[i++] = ( conf.getParameter<double>("EcalRecHitThreshEtA_Barrel") );


   i=0;

   hcalIsoBarrelRadiusA_[i++] = ( conf.getParameter<double>("HcalTowerInnerRadiusA_Barrel") );

   hcalIsoBarrelRadiusA_[i++] = ( conf.getParameter<double>("HcalTowerOuterRadiusA_Barrel") );

   hcalIsoBarrelRadiusA_[i++] = ( conf.getParameter<double>("HcalTowerThreshEA_Barrel") );

   hcalIsoBarrelRadiusA_[i++] = ( conf.getParameter<double>("HcalDepth1TowerInnerRadiusA_Barrel") );

   hcalIsoBarrelRadiusA_[i++] = ( conf.getParameter<double>("HcalDepth1TowerOuterRadiusA_Barrel") );

   hcalIsoBarrelRadiusA_[i++] = ( conf.getParameter<double>("HcalDepth1TowerThreshEA_Barrel") );

   hcalIsoBarrelRadiusA_[i++] = ( conf.getParameter<double>("HcalDepth2TowerInnerRadiusA_Barrel") );

   hcalIsoBarrelRadiusA_[i++] = ( conf.getParameter<double>("HcalDepth2TowerOuterRadiusA_Barrel") );

   hcalIsoBarrelRadiusA_[i++] = ( conf.getParameter<double>("HcalDepth2TowerThreshEA_Barrel") );


   i=0;

   trkIsoBarrelRadiusB_[i++] = (  conf.getParameter<double>("TrackConeOuterRadiusB_Barrel") );

   trkIsoBarrelRadiusB_[i++] = (  conf.getParameter<double>("TrackConeInnerRadiusB_Barrel") ) ;

   trkIsoBarrelRadiusB_[i++] = (  conf.getParameter<double>("isolationtrackThresholdB_Barrel") );

   trkIsoBarrelRadiusB_[i++] = (  conf.getParameter<double>("longImpactParameterB_Barrel") );

   trkIsoBarrelRadiusB_[i++] = (  conf.getParameter<double>("transImpactParameterB_Barrel") );

   trkIsoBarrelRadiusB_[i++] = (  conf.getParameter<double>("isolationtrackEtaSliceB_Barrel") );


   i=0;

   ecalIsoBarrelRadiusB_[i++] = ( conf.getParameter<double>("EcalRecHitInnerRadiusB_Barrel") );

   ecalIsoBarrelRadiusB_[i++] = ( conf.getParameter<double>("EcalRecHitOuterRadiusB_Barrel") );

   ecalIsoBarrelRadiusB_[i++] = ( conf.getParameter<double>("EcalRecHitEtaSliceB_Barrel") );

   ecalIsoBarrelRadiusB_[i++] = ( conf.getParameter<double>("EcalRecHitThreshEB_Barrel") );

   ecalIsoBarrelRadiusB_[i++] = ( conf.getParameter<double>("EcalRecHitThreshEtB_Barrel") );


   i=0;

   hcalIsoBarrelRadiusB_[i++] = (  conf.getParameter<double>("HcalTowerInnerRadiusB_Barrel") );

   hcalIsoBarrelRadiusB_[i++] = ( conf.getParameter<double>("HcalTowerOuterRadiusB_Barrel") );

   hcalIsoBarrelRadiusB_[i++] = ( conf.getParameter<double>("HcalTowerThreshEB_Barrel") );

   hcalIsoBarrelRadiusB_[i++] = ( conf.getParameter<double>("HcalDepth1TowerInnerRadiusB_Barrel") );

   hcalIsoBarrelRadiusB_[i++] = ( conf.getParameter<double>("HcalDepth1TowerOuterRadiusB_Barrel") );

   hcalIsoBarrelRadiusB_[i++] = ( conf.getParameter<double>("HcalDepth1TowerThreshEB_Barrel") );

   hcalIsoBarrelRadiusB_[i++] = ( conf.getParameter<double>("HcalDepth2TowerInnerRadiusB_Barrel") );

   hcalIsoBarrelRadiusB_[i++] = ( conf.getParameter<double>("HcalDepth2TowerOuterRadiusB_Barrel") );

   hcalIsoBarrelRadiusB_[i++] = ( conf.getParameter<double>("HcalDepth2TowerThreshEB_Barrel") );


   i=0;

   trkIsoEndcapRadiusA_[i++] = (  conf.getParameter<double>("TrackConeOuterRadiusA_Endcap") );

   trkIsoEndcapRadiusA_[i++] = (  conf.getParameter<double>("TrackConeInnerRadiusA_Endcap") ) ;

   trkIsoEndcapRadiusA_[i++] = (  conf.getParameter<double>("isolationtrackThresholdA_Endcap") );

   trkIsoEndcapRadiusA_[i++] = (  conf.getParameter<double>("longImpactParameterA_Endcap") );

   trkIsoEndcapRadiusA_[i++] = (  conf.getParameter<double>("transImpactParameterA_Endcap") );

   trkIsoEndcapRadiusA_[i++] = (  conf.getParameter<double>("isolationtrackEtaSliceA_Endcap") );


   i=0;

   ecalIsoEndcapRadiusA_[i++] = ( conf.getParameter<double>("EcalRecHitInnerRadiusA_Endcap") );

   ecalIsoEndcapRadiusA_[i++] = ( conf.getParameter<double>("EcalRecHitOuterRadiusA_Endcap") );

   ecalIsoEndcapRadiusA_[i++] = ( conf.getParameter<double>("EcalRecHitEtaSliceA_Endcap") );

   ecalIsoEndcapRadiusA_[i++] = ( conf.getParameter<double>("EcalRecHitThreshEA_Endcap") );

   ecalIsoEndcapRadiusA_[i++] = ( conf.getParameter<double>("EcalRecHitThreshEtA_Endcap") );


   i=0;

   hcalIsoEndcapRadiusA_[i++] = (  conf.getParameter<double>("HcalTowerInnerRadiusA_Endcap") );

   hcalIsoEndcapRadiusA_[i++] = ( conf.getParameter<double>("HcalTowerOuterRadiusA_Endcap") );

   hcalIsoEndcapRadiusA_[i++] = ( conf.getParameter<double>("HcalTowerThreshEA_Endcap") );

   hcalIsoEndcapRadiusA_[i++] = ( conf.getParameter<double>("HcalDepth1TowerInnerRadiusA_Endcap") );

   hcalIsoEndcapRadiusA_[i++] = ( conf.getParameter<double>("HcalDepth1TowerOuterRadiusA_Endcap") );

   hcalIsoEndcapRadiusA_[i++] = ( conf.getParameter<double>("HcalDepth1TowerThreshEA_Endcap") );

   hcalIsoEndcapRadiusA_[i++] = ( conf.getParameter<double>("HcalDepth2TowerInnerRadiusA_Endcap") );

   hcalIsoEndcapRadiusA_[i++] = ( conf.getParameter<double>("HcalDepth2TowerOuterRadiusA_Endcap") );

   hcalIsoEndcapRadiusA_[i++] = ( conf.getParameter<double>("HcalDepth2TowerThreshEA_Endcap") );


   i=0;

   trkIsoEndcapRadiusB_[i++] = (  conf.getParameter<double>("TrackConeOuterRadiusB_Endcap") );

   trkIsoEndcapRadiusB_[i++] = (  conf.getParameter<double>("TrackConeInnerRadiusB_Endcap") ) ;

   trkIsoEndcapRadiusB_[i++] = (  conf.getParameter<double>("isolationtrackThresholdB_Endcap") );

   trkIsoEndcapRadiusB_[i++] = (  conf.getParameter<double>("longImpactParameterB_Endcap") );

   trkIsoEndcapRadiusB_[i++] = (  conf.getParameter<double>("transImpactParameterB_Endcap") );

   trkIsoEndcapRadiusB_[i++] = (  conf.getParameter<double>("isolationtrackEtaSliceB_Endcap") );


   i=0;

   ecalIsoEndcapRadiusB_[i++] = ( conf.getParameter<double>("EcalRecHitInnerRadiusB_Endcap") );

   ecalIsoEndcapRadiusB_[i++] = ( conf.getParameter<double>("EcalRecHitOuterRadiusB_Endcap") );

   ecalIsoEndcapRadiusB_[i++] = ( conf.getParameter<double>("EcalRecHitEtaSliceB_Endcap") );

   ecalIsoEndcapRadiusB_[i++] = ( conf.getParameter<double>("EcalRecHitThreshEB_Endcap") );

   ecalIsoEndcapRadiusB_[i++] = ( conf.getParameter<double>("EcalRecHitThreshEtB_Endcap") );


   i=0;

   hcalIsoEndcapRadiusB_[i++] = (  conf.getParameter<double>("HcalTowerInnerRadiusB_Endcap") );

   hcalIsoEndcapRadiusB_[i++] = ( conf.getParameter<double>("HcalTowerOuterRadiusB_Endcap") );

   hcalIsoEndcapRadiusB_[i++] = ( conf.getParameter<double>("HcalTowerThreshEB_Endcap") );

   hcalIsoEndcapRadiusB_[i++] = ( conf.getParameter<double>("HcalDepth1TowerInnerRadiusB_Endcap") );

   hcalIsoEndcapRadiusB_[i++] = ( conf.getParameter<double>("HcalDepth1TowerOuterRadiusB_Endcap") );

   hcalIsoEndcapRadiusB_[i++] = ( conf.getParameter<double>("HcalDepth1TowerThreshEB_Endcap") );

   hcalIsoEndcapRadiusB_[i++] = ( conf.getParameter<double>("HcalDepth2TowerInnerRadiusB_Endcap") );

   hcalIsoEndcapRadiusB_[i++] = ( conf.getParameter<double>("HcalDepth2TowerOuterRadiusB_Endcap") );

   hcalIsoEndcapRadiusB_[i++] = ( conf.getParameter<double>("HcalDepth2TowerThreshEB_Endcap") );


   //Pick up the variables for the spike removal

   flagsEB_ = flagsEB;

   flagsEE_ = flagsEE;

   severityExclEB_ = severitiesEB;

   severityExclEE_ = severitiesEE;


 }


 void PhotonIsolationCalculator::calculate(const reco::Photon* pho,

                       const edm::Event& e,

                       const edm::EventSetup& es,

                       reco::Photon::FiducialFlags& phofid,

                       reco::Photon::IsolationVariables& phoisolR1,

                       reco::Photon::IsolationVariables& phoisolR2) const {


   //Get fiducial flags. This does not really belong here

   bool isEBPho     = false;

   bool isEEPho     = false;

   bool isEBEtaGap  = false;

   bool isEBPhiGap  = false;

   bool isEERingGap = false;

   bool isEEDeeGap  = false;

   bool isEBEEGap   = false;

   classify(pho, isEBPho, isEEPho, isEBEtaGap, isEBPhiGap, isEERingGap, isEEDeeGap, isEBEEGap);

   phofid.isEB = isEBPho;

   phofid.isEE = isEEPho;

   phofid.isEBEtaGap  = isEBEtaGap;

   phofid.isEBPhiGap  = isEBPhiGap;

   phofid.isEERingGap = isEERingGap;

   phofid.isEEDeeGap  = isEEDeeGap;

   phofid.isEBEEGap   = isEBEEGap;


   // Calculate isolation variables. cone sizes and thresholds

   // are set for Barrel and endcap separately


   reco::SuperClusterRef scRef=pho->superCluster();

   const reco::BasicCluster & seedCluster = *(scRef->seed()) ;

   DetId seedXtalId = seedCluster.hitsAndFractions()[0].first ;

   int detector = seedXtalId.subdetId() ;


   //Isolation parameters variables

   double photonEcalRecHitConeInnerRadiusA_;

   double photonEcalRecHitConeOuterRadiusA_;

   double photonEcalRecHitEtaSliceA_;

   double photonEcalRecHitThreshEA_;

   double photonEcalRecHitThreshEtA_;

   double photonHcalTowerConeInnerRadiusA_;

   double photonHcalTowerConeOuterRadiusA_;

   double photonHcalTowerThreshEA_;

   double photonHcalDepth1TowerConeInnerRadiusA_;

   double photonHcalDepth1TowerConeOuterRadiusA_;

   double photonHcalDepth1TowerThreshEA_;

   double photonHcalDepth2TowerConeInnerRadiusA_;

   double photonHcalDepth2TowerConeOuterRadiusA_;

   double photonHcalDepth2TowerThreshEA_;

   double trackConeOuterRadiusA_;

   double trackConeInnerRadiusA_;

   double isolationtrackThresholdA_;

   double isolationtrackEtaSliceA_;

   double trackLipRadiusA_;

   double trackD0RadiusA_;

   double photonEcalRecHitConeInnerRadiusB_;

   double photonEcalRecHitConeOuterRadiusB_;

   double photonEcalRecHitEtaSliceB_;

   double photonEcalRecHitThreshEB_;

   double photonEcalRecHitThreshEtB_;

   double photonHcalTowerConeInnerRadiusB_;

   double photonHcalTowerConeOuterRadiusB_;

   double photonHcalTowerThreshEB_;

   double photonHcalDepth1TowerConeInnerRadiusB_;

   double photonHcalDepth1TowerConeOuterRadiusB_;

   double photonHcalDepth1TowerThreshEB_;

   double photonHcalDepth2TowerConeInnerRadiusB_;

   double photonHcalDepth2TowerConeOuterRadiusB_;

   double photonHcalDepth2TowerThreshEB_;

   double trackConeOuterRadiusB_;

   double trackConeInnerRadiusB_;

   double isolationtrackThresholdB_;

   double isolationtrackEtaSliceB_;

   double trackLipRadiusB_;

   double trackD0RadiusB_;


   if (detector==EcalBarrel) {


     trackConeOuterRadiusA_    = trkIsoBarrelRadiusA_[0];

     trackConeInnerRadiusA_    = trkIsoBarrelRadiusA_[1];

     isolationtrackThresholdA_ = trkIsoBarrelRadiusA_[2];

     trackLipRadiusA_          = trkIsoBarrelRadiusA_[3];

     trackD0RadiusA_           = trkIsoBarrelRadiusA_[4];

     isolationtrackEtaSliceA_  = trkIsoBarrelRadiusA_[5];


     photonEcalRecHitConeInnerRadiusA_  = ecalIsoBarrelRadiusA_[0];

     photonEcalRecHitConeOuterRadiusA_  = ecalIsoBarrelRadiusA_[1];

     photonEcalRecHitEtaSliceA_         = ecalIsoBarrelRadiusA_[2];

     photonEcalRecHitThreshEA_          = ecalIsoBarrelRadiusA_[3];

     photonEcalRecHitThreshEtA_         = ecalIsoBarrelRadiusA_[4];


     photonHcalTowerConeInnerRadiusA_       = hcalIsoBarrelRadiusA_[0];

     photonHcalTowerConeOuterRadiusA_       = hcalIsoBarrelRadiusA_[1];

     photonHcalTowerThreshEA_               = hcalIsoBarrelRadiusA_[2];

     photonHcalDepth1TowerConeInnerRadiusA_ = hcalIsoBarrelRadiusA_[3];

     photonHcalDepth1TowerConeOuterRadiusA_ = hcalIsoBarrelRadiusA_[4];

     photonHcalDepth1TowerThreshEA_         = hcalIsoBarrelRadiusA_[5];

     photonHcalDepth2TowerConeInnerRadiusA_ = hcalIsoBarrelRadiusA_[6];

     photonHcalDepth2TowerConeOuterRadiusA_ = hcalIsoBarrelRadiusA_[7];

     photonHcalDepth2TowerThreshEA_         = hcalIsoBarrelRadiusA_[8];


     trackConeOuterRadiusB_    = trkIsoBarrelRadiusB_[0];

     trackConeInnerRadiusB_    = trkIsoBarrelRadiusB_[1];

     isolationtrackThresholdB_ = trkIsoBarrelRadiusB_[2];

     trackLipRadiusB_          = trkIsoBarrelRadiusB_[3];

     trackD0RadiusB_           = trkIsoBarrelRadiusB_[4];

     isolationtrackEtaSliceB_  = trkIsoBarrelRadiusB_[5];


     photonEcalRecHitConeInnerRadiusB_  = ecalIsoBarrelRadiusB_[0];

     photonEcalRecHitConeOuterRadiusB_  = ecalIsoBarrelRadiusB_[1];

     photonEcalRecHitEtaSliceB_         = ecalIsoBarrelRadiusB_[2];

     photonEcalRecHitThreshEB_          = ecalIsoBarrelRadiusB_[3];

     photonEcalRecHitThreshEtB_         = ecalIsoBarrelRadiusB_[4];


     photonHcalTowerConeInnerRadiusB_       = hcalIsoBarrelRadiusB_[0];

     photonHcalTowerConeOuterRadiusB_       = hcalIsoBarrelRadiusB_[1];

     photonHcalTowerThreshEB_               = hcalIsoBarrelRadiusB_[2];

     photonHcalDepth1TowerConeInnerRadiusB_ = hcalIsoBarrelRadiusB_[3];

     photonHcalDepth1TowerConeOuterRadiusB_ = hcalIsoBarrelRadiusB_[4];

     photonHcalDepth1TowerThreshEB_         = hcalIsoBarrelRadiusB_[5];

     photonHcalDepth2TowerConeInnerRadiusB_ = hcalIsoBarrelRadiusB_[6];

     photonHcalDepth2TowerConeOuterRadiusB_ = hcalIsoBarrelRadiusB_[7];

     photonHcalDepth2TowerThreshEB_         = hcalIsoBarrelRadiusB_[8];


   } else {

     // detector==EcalEndcap


     trackConeOuterRadiusA_    = trkIsoEndcapRadiusA_[0];

     trackConeInnerRadiusA_    = trkIsoEndcapRadiusA_[1];

     isolationtrackThresholdA_ = trkIsoEndcapRadiusA_[2];

     trackLipRadiusA_          = trkIsoEndcapRadiusA_[3];

     trackD0RadiusA_           = trkIsoEndcapRadiusA_[4];

     isolationtrackEtaSliceA_  = trkIsoEndcapRadiusA_[5];


     photonEcalRecHitConeInnerRadiusA_  = ecalIsoEndcapRadiusA_[0];

     photonEcalRecHitConeOuterRadiusA_  = ecalIsoEndcapRadiusA_[1];

     photonEcalRecHitEtaSliceA_         = ecalIsoEndcapRadiusA_[2];

     photonEcalRecHitThreshEA_          = ecalIsoEndcapRadiusA_[3];

     photonEcalRecHitThreshEtA_         = ecalIsoEndcapRadiusA_[4];


     photonHcalTowerConeInnerRadiusA_       = hcalIsoEndcapRadiusA_[0];

     photonHcalTowerConeOuterRadiusA_       = hcalIsoEndcapRadiusA_[1];

     photonHcalTowerThreshEA_               = hcalIsoEndcapRadiusA_[2];

     photonHcalDepth1TowerConeInnerRadiusA_ = hcalIsoEndcapRadiusA_[3];

     photonHcalDepth1TowerConeOuterRadiusA_ = hcalIsoEndcapRadiusA_[4];

     photonHcalDepth1TowerThreshEA_         = hcalIsoEndcapRadiusA_[5];

     photonHcalDepth2TowerConeInnerRadiusA_ = hcalIsoEndcapRadiusA_[6];

     photonHcalDepth2TowerConeOuterRadiusA_ = hcalIsoEndcapRadiusA_[7];

     photonHcalDepth2TowerThreshEA_         = hcalIsoEndcapRadiusA_[8];


     trackConeOuterRadiusB_    = trkIsoEndcapRadiusB_[0];

     trackConeInnerRadiusB_    = trkIsoEndcapRadiusB_[1];

     isolationtrackThresholdB_ = trkIsoEndcapRadiusB_[2];

     trackLipRadiusB_          = trkIsoEndcapRadiusB_[3];

     trackD0RadiusB_           = trkIsoEndcapRadiusB_[4];

     isolationtrackEtaSliceB_  = trkIsoEndcapRadiusB_[5];


     photonEcalRecHitConeInnerRadiusB_  = ecalIsoEndcapRadiusB_[0];

     photonEcalRecHitConeOuterRadiusB_  = ecalIsoEndcapRadiusB_[1];

     photonEcalRecHitEtaSliceB_         = ecalIsoEndcapRadiusB_[2];

     photonEcalRecHitThreshEB_          = ecalIsoEndcapRadiusB_[3];

     photonEcalRecHitThreshEtB_         = ecalIsoEndcapRadiusB_[4];


     photonHcalTowerConeInnerRadiusB_       = hcalIsoEndcapRadiusB_[0];

     photonHcalTowerConeOuterRadiusB_       = hcalIsoEndcapRadiusB_[1];

     photonHcalTowerThreshEB_               = hcalIsoEndcapRadiusB_[2];

     photonHcalDepth1TowerConeInnerRadiusB_ = hcalIsoEndcapRadiusB_[3];

     photonHcalDepth1TowerConeOuterRadiusB_ = hcalIsoEndcapRadiusB_[4];

     photonHcalDepth1TowerThreshEB_         = hcalIsoEndcapRadiusB_[5];

     photonHcalDepth2TowerConeInnerRadiusB_ = hcalIsoEndcapRadiusB_[6];

     photonHcalDepth2TowerConeOuterRadiusB_ = hcalIsoEndcapRadiusB_[7];

     photonHcalDepth2TowerThreshEB_         = hcalIsoEndcapRadiusB_[8];


   }


   //Calculate hollow cone track isolation, CONE A

   int ntrkA=0;

   double trkisoA=0;

   calculateTrackIso(pho, e, trkisoA, ntrkA, isolationtrackThresholdA_,

             trackConeOuterRadiusA_, trackConeInnerRadiusA_, isolationtrackEtaSliceA_, trackLipRadiusA_, trackD0RadiusA_);


   //Calculate solid cone track isolation, CONE A

   int sntrkA=0;

   double strkisoA=0;

   calculateTrackIso(pho, e, strkisoA, sntrkA, isolationtrackThresholdA_,

             trackConeOuterRadiusA_, 0., isolationtrackEtaSliceA_, trackLipRadiusA_, trackD0RadiusA_ );


   phoisolR1.nTrkHollowCone = ntrkA;

   phoisolR1.trkSumPtHollowCone = trkisoA;

   phoisolR1.nTrkSolidCone = sntrkA;

   phoisolR1.trkSumPtSolidCone = strkisoA;


   //Calculate hollow cone track isolation, CONE B

   int ntrkB=0;

   double trkisoB=0;

   calculateTrackIso(pho, e, trkisoB, ntrkB, isolationtrackThresholdB_,

             trackConeOuterRadiusB_, trackConeInnerRadiusB_,  isolationtrackEtaSliceB_, trackLipRadiusB_, trackD0RadiusB_ );


   //Calculate solid cone track isolation, CONE B

   int sntrkB=0;

   double strkisoB=0;

   calculateTrackIso(pho, e, strkisoB, sntrkB, isolationtrackThresholdB_,

             trackConeOuterRadiusB_, 0., isolationtrackEtaSliceB_, trackLipRadiusB_, trackD0RadiusB_);


   phoisolR2.nTrkHollowCone = ntrkB;

   phoisolR2.trkSumPtHollowCone = trkisoB;

   phoisolR2.nTrkSolidCone = sntrkB;

   phoisolR2.trkSumPtSolidCone = strkisoB;


 //   std::cout << "Output from solid cone track isolation: ";

 //   std::cout << " Sum pT: " << strkiso << " ntrk: " << sntrk << std::endl;


   double EcalRecHitIsoA = calculateEcalRecHitIso(pho, e, es,

                          photonEcalRecHitConeOuterRadiusA_,

                          photonEcalRecHitConeInnerRadiusA_,

                          photonEcalRecHitEtaSliceA_,

                          photonEcalRecHitThreshEA_,

                          photonEcalRecHitThreshEtA_,

                          vetoClusteredEcalHits_,

                          useNumCrystals_);

   phoisolR1.ecalRecHitSumEt = EcalRecHitIsoA;


   double EcalRecHitIsoB = calculateEcalRecHitIso(pho, e, es,

                          photonEcalRecHitConeOuterRadiusB_,

                          photonEcalRecHitConeInnerRadiusB_,

                          photonEcalRecHitEtaSliceB_,

                          photonEcalRecHitThreshEB_,

                          photonEcalRecHitThreshEtB_,

                          vetoClusteredEcalHits_,

                          useNumCrystals_);

   phoisolR2.ecalRecHitSumEt = EcalRecHitIsoB;


   double HcalTowerIsoA = calculateHcalTowerIso(pho, e, es, photonHcalTowerConeOuterRadiusA_,

                           photonHcalTowerConeInnerRadiusA_,

                           photonHcalTowerThreshEA_, -1 );

   phoisolR1.hcalTowerSumEt = HcalTowerIsoA;


   double HcalTowerIsoB = calculateHcalTowerIso(pho, e, es, photonHcalTowerConeOuterRadiusB_,

                           photonHcalTowerConeInnerRadiusB_,

                           photonHcalTowerThreshEB_, -1 );

   phoisolR2.hcalTowerSumEt = HcalTowerIsoB;


   double HcalDepth1TowerIsoA = calculateHcalTowerIso(pho, e, es, photonHcalDepth1TowerConeOuterRadiusA_,

                           photonHcalDepth1TowerConeInnerRadiusA_,

                           photonHcalDepth1TowerThreshEA_, 1 );

   phoisolR1.hcalDepth1TowerSumEt = HcalDepth1TowerIsoA;


   double HcalDepth1TowerIsoB = calculateHcalTowerIso(pho, e, es, photonHcalDepth1TowerConeOuterRadiusB_,

                           photonHcalDepth1TowerConeInnerRadiusB_,

                           photonHcalDepth1TowerThreshEB_, 1 );

   phoisolR2.hcalDepth1TowerSumEt = HcalDepth1TowerIsoB;


   double HcalDepth2TowerIsoA = calculateHcalTowerIso(pho, e, es, photonHcalDepth2TowerConeOuterRadiusA_,

                                photonHcalDepth2TowerConeInnerRadiusA_,

                                photonHcalDepth2TowerThreshEA_, 2 );

   phoisolR1.hcalDepth2TowerSumEt = HcalDepth2TowerIsoA;


   double HcalDepth2TowerIsoB = calculateHcalTowerIso(pho, e, es, photonHcalDepth2TowerConeOuterRadiusB_,

                           photonHcalDepth2TowerConeInnerRadiusB_,

                           photonHcalDepth2TowerThreshEB_, 2 );

   phoisolR2.hcalDepth2TowerSumEt = HcalDepth2TowerIsoB;


   // New Hcal isolation based on the new H/E definition (towers behind the BCs in the SC are used to evaluated H)

   double HcalTowerBcIsoA = calculateHcalTowerIso(pho, e, es, photonHcalTowerConeOuterRadiusA_,

                          photonHcalTowerThreshEA_, -1 );

   phoisolR1.hcalTowerSumEtBc = HcalTowerBcIsoA;


   double HcalTowerBcIsoB = calculateHcalTowerIso(pho, e, es, photonHcalTowerConeOuterRadiusB_,

                          photonHcalTowerThreshEB_, -1 );

   phoisolR2.hcalTowerSumEtBc = HcalTowerBcIsoB;


   double HcalDepth1TowerBcIsoA = calculateHcalTowerIso(pho, e, es, photonHcalDepth1TowerConeOuterRadiusA_,

                                photonHcalDepth1TowerThreshEA_, 1 );

   phoisolR1.hcalDepth1TowerSumEtBc = HcalDepth1TowerBcIsoA;


   double HcalDepth1TowerBcIsoB = calculateHcalTowerIso(pho, e, es, photonHcalDepth1TowerConeOuterRadiusB_,

                                photonHcalDepth1TowerThreshEB_, 1 );

   phoisolR2.hcalDepth1TowerSumEtBc = HcalDepth1TowerBcIsoB;


   double HcalDepth2TowerBcIsoA = calculateHcalTowerIso(pho, e, es, photonHcalDepth2TowerConeOuterRadiusA_,

                                photonHcalDepth2TowerThreshEA_, 2 );

   phoisolR1.hcalDepth2TowerSumEtBc = HcalDepth2TowerBcIsoA;


   double HcalDepth2TowerBcIsoB = calculateHcalTowerIso(pho, e, es, photonHcalDepth2TowerConeOuterRadiusB_,

                                photonHcalDepth2TowerThreshEB_, 2 );

   phoisolR2.hcalDepth2TowerSumEtBc = HcalDepth2TowerBcIsoB;


 }


 void PhotonIsolationCalculator::classify(const reco::Photon* photon,

                      bool &isEBPho,

                      bool &isEEPho,

                      bool &isEBEtaGap,

                      bool &isEBPhiGap,

                      bool &isEERingGap,

                      bool &isEEDeeGap,

                      bool &isEBEEGap){


   const reco::CaloCluster & seedCluster = *(photon->superCluster()->seed()) ;

   // following line is temporary until the D. Evans or F. Ferri submit the new tag for ClusterAlgos

   // DEfinitive will be

   // DetId seedXtalId = scRef->seed()->seed();

   DetId seedXtalId = seedCluster.hitsAndFractions()[0].first ;

   int detector = seedXtalId.subdetId() ;


   //Set fiducial flags for this photon.

   double eta = photon->superCluster()->position().eta();

   double feta = fabs(eta);


   if (detector==EcalBarrel) {


     isEBPho = true;

     if (EBDetId::isNextToEtaBoundary(EBDetId(seedXtalId))) {

       if (fabs(feta-1.479)<.1) isEBEEGap = true ;

       else

     isEBEtaGap = true ;

     }


     if (EBDetId::isNextToPhiBoundary(EBDetId(seedXtalId)))

       isEBPhiGap = true ;


   } else if (detector==EcalEndcap) {


     isEEPho = true;

     if (EEDetId::isNextToRingBoundary(EEDetId(seedXtalId))) {

       if (fabs(feta-1.479)<.1) isEBEEGap = true ;

       else

     isEERingGap = true ;

     }


     if (EEDetId::isNextToDBoundary(EEDetId(seedXtalId)))

       isEEDeeGap = true ;

   }


 }


 void PhotonIsolationCalculator::calculateTrackIso(const reco::Photon* photon,

                                                   const edm::Event& e,

                                                   double &trkCone,

                                                   int &ntrkCone,

                                                   double pTThresh,

                                                   double RCone,

                                                   double RinnerCone,

                                                   double etaSlice,

                                                   double lip,

                                                   double d0)  const {


   ntrkCone =0;trkCone=0;

   //get the tracks

   edm::Handle<reco::TrackCollection> tracks;

   e.getByLabel(trackInputTag_,tracks);

   if(!tracks.isValid()) {

     return;

   }

   const reco::TrackCollection* trackCollection = tracks.product();

   //Photon Eta and Phi.  Hope these are correct.

   reco::BeamSpot vertexBeamSpot;

   edm::Handle<reco::BeamSpot> recoBeamSpotHandle;

   e.getByLabel(beamSpotProducerTag_,recoBeamSpotHandle);

   vertexBeamSpot = *recoBeamSpotHandle;


   PhotonTkIsolation phoIso(RCone,

                RinnerCone,

                etaSlice,

                pTThresh,

                lip ,

                d0,

                trackCollection,

                math::XYZPoint(vertexBeamSpot.x0(),vertexBeamSpot.y0(),vertexBeamSpot.z0()));


   std::pair<int,double> res = phoIso.getIso(photon);

   ntrkCone = res.first;

   trkCone = res.second;

 }


 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 useNumXtals)  const

 {


   edm::Handle<EcalRecHitCollection> ecalhitsCollEB;

   edm::Handle<EcalRecHitCollection> ecalhitsCollEE;

   iEvent.getByLabel(endcapecalCollection_, ecalhitsCollEE);


   iEvent.getByLabel(barrelecalCollection_, ecalhitsCollEB);


   const EcalRecHitCollection* rechitsCollectionEE_ = ecalhitsCollEE.product();

   const EcalRecHitCollection* rechitsCollectionEB_ = ecalhitsCollEB.product();


   edm::ESHandle<EcalSeverityLevelAlgo> sevlv;

   iSetup.get<EcalSeverityLevelAlgoRcd>().get(sevlv);

   const EcalSeverityLevelAlgo* sevLevel = sevlv.product();


   EcalRecHitMetaCollection RecHitsEE(*rechitsCollectionEE_);

   EcalRecHitMetaCollection RecHitsEB(*rechitsCollectionEB_);


   edm::ESHandle<CaloGeometry> geoHandle;

   iSetup.get<CaloGeometryRecord>().get(geoHandle);


   EgammaRecHitIsolation phoIsoEB(RCone,

                                  RConeInner,

                                  etaSlice,

                                  etMin,

                                  eMin,

                                  geoHandle,

                                  &RecHitsEB,

                                  sevLevel,

                                  DetId::Ecal);


   phoIsoEB.setVetoClustered(vetoClusteredHits);

   phoIsoEB.setUseNumCrystals(useNumXtals);

   phoIsoEB.doSeverityChecks(ecalhitsCollEB.product(), severityExclEB_);

   phoIsoEB.doFlagChecks(flagsEB_);

   double ecalIsolEB = phoIsoEB.getEtSum(photon);


   EgammaRecHitIsolation phoIsoEE(RCone,

                                  RConeInner,

                                  etaSlice,

                                  etMin,

                                  eMin,

                                  geoHandle,

                                  &RecHitsEE,

                                  sevLevel,

                                  DetId::Ecal);


   phoIsoEE.setVetoClustered(vetoClusteredHits);

   phoIsoEE.setUseNumCrystals(useNumXtals);

   phoIsoEE.doSeverityChecks(ecalhitsCollEE.product(), severityExclEE_);

   phoIsoEE.doFlagChecks(flagsEE_);


   double ecalIsolEE = phoIsoEE.getEtSum(photon);

   //  delete phoIso;

   double ecalIsol = ecalIsolEB + ecalIsolEE;


   return ecalIsol;

 }


 double PhotonIsolationCalculator::calculateHcalTowerIso(const reco::Photon* photon,

                             const edm::Event& iEvent,

                             const edm::EventSetup& iSetup,

                             double RCone,

                             double RConeInner,

                             double eMin,

                             signed int depth )  const

 {


   edm::Handle<CaloTowerCollection> hcalhitsCollH;


   iEvent.getByLabel(hcalCollection_, hcalhitsCollH);


   const CaloTowerCollection *toww = hcalhitsCollH.product();


   double hcalIsol=0.;


   //std::cout << "before iso call" << std::endl;

   EgammaTowerIsolation phoIso(RCone,

                               RConeInner,

                               eMin,depth,

                               toww);

   hcalIsol = phoIso.getTowerEtSum(photon);

   //  delete phoIso;

   //std::cout << "after call" << std::endl;

   return hcalIsol;


 }


 double PhotonIsolationCalculator::calculateHcalTowerIso(const reco::Photon* photon,

                             const edm::Event& iEvent,

                             const edm::EventSetup& iSetup,

                             double RCone,

                             double eMin,

                             signed int depth )  const

 {


   edm::Handle<CaloTowerCollection> hcalhitsCollH;


   iEvent.getByLabel(hcalCollection_, hcalhitsCollH);


   const CaloTowerCollection *toww = hcalhitsCollH.product();


   double hcalIsol=0.;


   //std::cout << "before iso call" << std::endl;

   EgammaTowerIsolation phoIso(RCone,

                   0.,

                               eMin,depth,

                               toww);

   hcalIsol = phoIso.getTowerEtSum(photon, &(photon->hcalTowersBehindClusters()) );

   //  delete phoIso;

   //std::cout << "after call" << std::endl;

   return hcalIsol;


 }

PhotonIsolationCalculator::classify
static void classify(const reco::Photon *photon, bool &isEBPho, bool &isEEPho, bool &isEBEtaGap, bool &isEBPhiGap, bool &isEERingGap, bool &isEEDeeGap, bool &isEBEEGap) dso_internal
Definition: PhotonIsolationCalculator.cc:495

edm::ParameterSet::getParameter
T getParameter(std::string const &) const

reco::BeamSpot::z0
double z0() const
z coordinate
Definition: BeamSpot.h:68

reco::Photon::IsolationVariables::hcalDepth1TowerSumEtBc
float hcalDepth1TowerSumEtBc
Definition: Photon.h:303

PhotonIsolationCalculator::ecalIsoBarrelRadiusA_
double ecalIsoBarrelRadiusA_[5]
Definition: PhotonIsolationCalculator.h:109

i
int i
Definition: DBlmapReader.cc:9

reco::Photon::IsolationVariables::hcalDepth1TowerSumEt
float hcalDepth1TowerSumEt
Definition: Photon.h:297

PhotonIsolationCalculator::flagsEE_
std::vector< int > flagsEE_
Definition: PhotonIsolationCalculator.h:124

EBDetId::isNextToEtaBoundary
static bool isNextToEtaBoundary(EBDetId id)
Definition: EBDetId.cc:125

EBDetId
Definition: EBDetId.h:19

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:180

EcalSeverityLevelAlgoRcd.h

EcalChannelStatusRcd.h

PhotonIsolationCalculator::modulePhiBoundary_
double modulePhiBoundary_
Definition: PhotonIsolationCalculator.h:103

edm::SortedCollection< EcalRecHit >

PhotonIsolationCalculator::hcalIsoBarrelRadiusB_
double hcalIsoBarrelRadiusB_[9]
Definition: PhotonIsolationCalculator.h:113

SimpleCaloRecHitMetaCollection
Definition: CaloRecHitMetaCollections.h:9

reco::Photon::FiducialFlags::isEEDeeGap
bool isEEDeeGap
Definition: Photon.h:99

EgammaHcalIsolation.h

reco::Photon::FiducialFlags::isEBEEGap
bool isEBEEGap
Definition: Photon.h:100

PhotonIsolationCalculator.h

EgammaRecHitIsolation::getEtSum
double getEtSum(const reco::Candidate *emObject) const
Definition: EgammaRecHitIsolation.h:41

EcalSeverityLevelAlgoRcd
Definition: EcalSeverityLevelAlgoRcd.h:13

PhotonIsolationCalculator::useNumCrystals_
bool useNumCrystals_
Definition: PhotonIsolationCalculator.h:106

EgammaEcalIsolation.h

reco::TrackCollection
std::vector< Track > TrackCollection
collection of Tracks
Definition: TrackFwd.h:10

reco::Photon::IsolationVariables::nTrkSolidCone
int nTrkSolidCone
Definition: Photon.h:311

EcalRecHitCollections.h

PhotonIsolationCalculator::trkIsoEndcapRadiusB_
double trkIsoEndcapRadiusB_[6]
Definition: PhotonIsolationCalculator.h:118

reco::Photon::superCluster
reco::SuperClusterRef superCluster() const
Ref to SuperCluster.
Definition: Photon.cc:58

EcalRecHit.h

PhotonIsolationCalculator::barrelecalCollection_
edm::InputTag barrelecalCollection_
Definition: PhotonIsolationCalculator.h:97

reco::CaloCluster::hitsAndFractions
const std::vector< std::pair< DetId, float > > & hitsAndFractions() const
Definition: CaloCluster.h:190

TrackFwd.h

PhotonIsolationCalculator::severityExclEE_
std::vector< int > severityExclEE_
Definition: PhotonIsolationCalculator.h:126

PhotonIsolationCalculator::trkIsoBarrelRadiusB_
double trkIsoBarrelRadiusB_[6]
Definition: PhotonIsolationCalculator.h:111

CaloRecHitMetaCollections.h

CaloTopologyRecord.h

eta
T eta() const
Definition: Basic3DVectorLD.h:179

edm::Handle< reco::TrackCollection >

reco::Photon::IsolationVariables::hcalDepth2TowerSumEt
float hcalDepth2TowerSumEt
Definition: Photon.h:299

reco::Photon
Definition: Photon.h:22

reco::BasicCluster

CaloGeometryRecord
Definition: CaloGeometryRecord.h:26

PhotonIsolationCalculator::hcalIsoBarrelRadiusA_
double hcalIsoBarrelRadiusA_[9]
Definition: PhotonIsolationCalculator.h:110

EgammaRecHitIsolation::doFlagChecks
void doFlagChecks(const std::vector< int > &v)
Definition: EgammaRecHitIsolation.h:57

funct::true
true
Definition: Factorize.h:183

PhotonIsolationCalculator::endcapecalCollection_
edm::InputTag endcapecalCollection_
Definition: PhotonIsolationCalculator.h:98

EgammaRecHitIsolation::setUseNumCrystals
void setUseNumCrystals(bool b=true)
Definition: EgammaRecHitIsolation.h:47

PhotonIsolationCalculator::hcalIsoEndcapRadiusA_
double hcalIsoEndcapRadiusA_[9]
Definition: PhotonIsolationCalculator.h:117

PhotonIsolationCalculator::flagsEB_
std::vector< int > flagsEB_
Definition: PhotonIsolationCalculator.h:123

iEvent
int iEvent
Definition: GenABIO.cc:243

DetId::Ecal
Definition: DetId.h:24

EEDetId
Definition: EEDetId.h:15

CaloGeometryRecord.h

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) const dso_internal
Definition: PhotonIsolationCalculator.cc:658

PhotonIsolationCalculator::hcalCollection_
edm::InputTag hcalCollection_
Definition: PhotonIsolationCalculator.h:99

reco::Photon::FiducialFlags::isEBPhiGap
bool isEBPhiGap
Definition: Photon.h:97

EcalBarrelTopology.h

PhotonIsolationCalculator::severityExclEB_
std::vector< int > severityExclEB_
Definition: PhotonIsolationCalculator.h:125

PhotonIsolationCalculator::trkIsoBarrelRadiusA_
double trkIsoBarrelRadiusA_[6]
Definition: PhotonIsolationCalculator.h:108

EcalEndcapTopology.h

edm::ESHandle< EcalSeverityLevelAlgo >

EBDetId::isNextToPhiBoundary
static bool isNextToPhiBoundary(EBDetId id)
Definition: EBDetId.cc:130

PhotonIsolationCalculator::trackInputTag_
edm::InputTag trackInputTag_
Definition: PhotonIsolationCalculator.h:101

EEDetId::isNextToRingBoundary
static bool isNextToRingBoundary(EEDetId id)
Definition: EEDetId.cc:375

reco::Photon::FiducialFlags::isEE
bool isEE
Definition: Photon.h:95

PhotonTkIsolation
Definition: PhotonTkIsolation.h:25

reco::CaloCluster
Definition: CaloCluster.h:29

edm::EventSetup
Definition: EventSetup.h:44

EgammaTowerIsolation.h

PhotonIsolationCalculator::ecalIsoBarrelRadiusB_
double ecalIsoBarrelRadiusB_[5]
Definition: PhotonIsolationCalculator.h:112

PhotonTkIsolation.h

PhotonIsolationCalculator::ecalIsoEndcapRadiusB_
double ecalIsoEndcapRadiusB_[5]
Definition: PhotonIsolationCalculator.h:119

EEDetId::isNextToDBoundary
static bool isNextToDBoundary(EEDetId id)
Definition: EEDetId.cc:367

edm::HandleBase::isValid
bool isValid() const
Definition: HandleBase.h:76

DetId::subdetId
int subdetId() const
get the contents of the subdetector field (not cast into any detector&#39;s numbering enum) ...
Definition: DetId.h:37

edm::Event::getByLabel
bool getByLabel(InputTag const &tag, Handle< PROD > &result) const
Definition: Event.h:390

PhotonTkIsolation::getIso
std::pair< int, float > getIso(const reco::Candidate *) const
Definition: PhotonTkIsolation.cc:72

EgammaRecHitIsolation::doSeverityChecks
void doSeverityChecks(const EcalRecHitCollection *const recHits, const std::vector< int > &v)
Definition: EgammaRecHitIsolation.h:49

reco::Photon::IsolationVariables::hcalDepth2TowerSumEtBc
float hcalDepth2TowerSumEtBc
Definition: Photon.h:305

dbtoconf.conf
tuple conf
Definition: dbtoconf.py:185

BasicClusterFwd.h

reco::Photon::IsolationVariables
Definition: Photon.h:288

StringToEnumValue.h

GsfElectron.h

reco::Photon::IsolationVariables::hcalTowerSumEt
float hcalTowerSumEt
Definition: Photon.h:295

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) const dso_internal
Definition: PhotonIsolationCalculator.cc:588

DetId
Definition: DetId.h:18

reco::Photon::IsolationVariables::trkSumPtSolidCone
float trkSumPtSolidCone
Definition: Photon.h:307

GsfElectronFwd.h

EgammaRecHitIsolation.h

CaloGeometry.h

reco::Photon::FiducialFlags::isEB
bool isEB
Definition: Photon.h:94

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_)
Definition: PhotonIsolationCalculator.cc:46

testEve_cfg.tracks
tuple tracks
Definition: testEve_cfg.py:39

PhotonIsolationCalculator::beamSpotProducerTag_
edm::InputTag beamSpotProducerTag_
Definition: PhotonIsolationCalculator.h:102

math::XYZPoint
XYZPointD XYZPoint
point in space with cartesian internal representation
Definition: Point3D.h:12

edm::EventSetup::get
const T & get() const
Definition: EventSetup.h:55

edm::ESHandle::product
T const * product() const
Definition: ESHandle.h:62

reco::Photon::FiducialFlags::isEERingGap
bool isEERingGap
Definition: Photon.h:98

edm::Handle::product
T const * product() const
Definition: Handle.h:81

reco::Photon::IsolationVariables::ecalRecHitSumEt
float ecalRecHitSumEt
Definition: Photon.h:293

PhotonIsolationCalculator::trkIsoEndcapRadiusA_
double trkIsoEndcapRadiusA_[6]
Definition: PhotonIsolationCalculator.h:115

alignCSCRings.e
list e
Definition: alignCSCRings.py:90

HcalRecHitCollections.h

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) const dso_internal
Definition: PhotonIsolationCalculator.cc:547

edm::InputTag
Definition: InputTag.h:17

PhotonIsolationCalculator::moduleEtaBoundary_
std::vector< double > moduleEtaBoundary_
Definition: PhotonIsolationCalculator.h:104

EgammaRecHitIsolation::setVetoClustered
void setVetoClustered(bool b=true)
Definition: EgammaRecHitIsolation.h:48

EgammaTowerIsolation::getTowerEtSum
double getTowerEtSum(const reco::Candidate *cand, const std::vector< CaloTowerDetId > *detIdToExclude=0) const
Definition: EgammaTowerIsolation.h:212

BeamSpot.h

EcalBarrel
Definition: EcalSubdetector.h:11

reco::BeamSpot::y0
double y0() const
y coordinate
Definition: BeamSpot.h:66

reco::Photon::hcalTowersBehindClusters
const std::vector< CaloTowerDetId > & hcalTowersBehindClusters() const
Definition: Photon.h:181

edm::ParameterSet
Definition: ParameterSet.h:35

EBDetId.h

EEDetId.h

EcalEndcap
Definition: EcalSubdetector.h:11

reco::Photon::FiducialFlags::isEBEtaGap
bool isEBEtaGap
Definition: Photon.h:96

reco::Photon::IsolationVariables::nTrkHollowCone
int nTrkHollowCone
Definition: Photon.h:313

reco::Photon::IsolationVariables::hcalTowerSumEtBc
float hcalTowerSumEtBc
Definition: Photon.h:301

edm::Event
Definition: Event.h:62

EgammaRecHitIsolation
Definition: EgammaRecHitIsolation.h:27

Track.h

BasicCluster.h

PhotonIsolationCalculator::ecalIsoEndcapRadiusA_
double ecalIsoEndcapRadiusA_[5]
Definition: PhotonIsolationCalculator.h:116

DetId.h

reco::BeamSpot
Definition: BeamSpot.h:22

reco::Photon::FiducialFlags
Definition: Photon.h:90

edm::Ref< SuperClusterCollection >

EgammaTowerIsolation
Definition: EgammaTowerIsolation.h:200

EcalSeverityLevelAlgo
Definition: EcalSeverityLevelAlgo.h:35

CaloTopology.h

reco::Photon::IsolationVariables::trkSumPtHollowCone
float trkSumPtHollowCone
Definition: Photon.h:309

reco::BeamSpot::x0
double x0() const
x coordinate
Definition: BeamSpot.h:64

PhotonIsolationCalculator::hcalIsoEndcapRadiusB_
double hcalIsoEndcapRadiusB_[9]
Definition: PhotonIsolationCalculator.h:120

PhotonIsolationCalculator::vetoClusteredEcalHits_
bool vetoClusteredEcalHits_
Definition: PhotonIsolationCalculator.h:105

SuperCluster.h