/afs/cern.ch/work/a/aaltunda/public/www/CMSSW_5_3_13_patch3/src/RecoEgamma/EgammaTools/interface/ggPFPhotons.h

Go to the documentation of this file.

#ifndef ggPFPhotons_h
#define ggPFPhotons_h
#include "DataFormats/PatCandidates/interface/Photon.h"

#include "DataFormats/VertexReco/interface/Vertex.h"
#include "DataFormats/VertexReco/interface/VertexFwd.h"
#include "DataFormats/ParticleFlowCandidate/interface/PFCandidateFwd.h"
#include "RecoEgamma/EgammaTools/interface/ggPFClusters.h"
#include "RecoEgamma/EgammaTools/interface/ggPFESClusters.h"
#include "RecoEgamma/EgammaTools/interface/ggPFTracks.h"
#include "RecoEcal/EgammaCoreTools/interface/Mustache.h"
#include <memory>
using namespace edm;
using namespace std;
using namespace reco;
//PFPHOTON CLASS BY RISHI PATEL FOR HiggsToGammaGamma rpatel@cern.ch
class ggPFPhotons  {

 public:
  
  explicit ggPFPhotons(
                       reco::Photon phot,
                       edm::Handle<reco::PhotonCollection>& pfPhotons,
                       edm::Handle<reco::GsfElectronCollection>& pfElectrons,
                       edm::Handle<PFCandidateCollection>& pfCandidates,
                       edm::Handle<EcalRecHitCollection>& EBReducedRecHits,
                       edm::Handle<EcalRecHitCollection>& EEReducedRecHits,
                       edm::Handle<EcalRecHitCollection>& ESRecHits,
                       const CaloSubdetectorGeometry* geomBar,
                       const CaloSubdetectorGeometry* geomEnd,
                       edm::Handle<BeamSpot>& beamSpotHandle
                       );
  //add Geometry for Cluster Shape Calculation??
  virtual ~ggPFPhotons();
  bool MatchPFReco(){ return matchPFReco_;}
  bool isConv(){ return isConv_;}
  bool hasSLConv(){return hasSLConv_;}
  bool isPFEle(){return isPFEle_;}
  bool EleVeto(){return EleVeto_;}
  float PFPS1(){return PFPreShower1_;}
  float PFPS2(){return PFPreShower2_;}
  float MustE(){return EinMustache_;}
  float MustEOut(){return MustacheEOut_;}
  float MustEtOut(){return MustacheEtOut_;}
  float PFLowE(){return PFLowClusE_;}
  double PFdEta(){return dEtaLowestC_;}
  double PFdPhi(){return dPhiLowestC_;}
  double PFClusRMSTot(){return PFClPhiRMS_;}
  double PFClusRMSMust(){return PFClPhiRMSMust_;}
  std::vector<reco::CaloCluster>PFClusters(){return PFClusters_;}
  std::vector<reco::CaloCluster>PFClustersSCFP(){return PFSCFootprintClusters_;}//Electron Veto:
  bool PFElectronVeto(edm::Handle<reco::ConversionCollection>& convH,
                      edm::Handle<reco::GsfElectronCollection>& gsfElectronsHandle);
  //for Vertex
  std::pair<float, float> SLPoint();
  //fill PFCLuster Variables
  void fillPFClusters();
  //Local Containment Correction
  double getPFPhoECorr( std::vector<reco::CaloCluster>PFClusters, const GBRForest *ReaderLCEB, const GBRForest *ReaderLCEE);
  std::pair<double, double>CalcRMS(vector<reco::CaloCluster> PFClust, reco::Photon PFPhoton);
  std::vector<reco::CaloCluster>recoPhotonClusterLink(          
                                                      reco::Photon phot, 
                                                      edm::Handle<PFCandidateCollection>& pfCandidates
                                                      );
  std::pair<double, double>SuperClusterSize(reco::Photon phot
                                              );
  static void recoPhotonClusterLink(
                                                             reco::SuperCluster sc, 
                                                             std::vector<reco::PFCandidatePtr>&insideMust, 
                                                             std::vector<reco::PFCandidatePtr>&outsideMust,
                                                             edm::Handle<PFCandidateCollection>& pfCandidates,
                                                             double etabound,
                                                             double phibound
                                                             );
  static std::pair<double, double>SuperClusterSize(reco::SuperCluster sc,
                                                   Handle<EcalRecHitCollection>&   EBReducedRecHits,
                                                   Handle<EcalRecHitCollection>&   EEReducedRecHits,
                                                   const CaloSubdetectorGeometry* geomBar,
                                                   const CaloSubdetectorGeometry* geomEnd
                                            );
  void PhotonPFCandMatch(  
                          reco::SuperCluster sc,                
                          std::vector<reco::PFCandidatePtr>&insideBox, 
                          edm::Handle<PFCandidateCollection>& pfCandidates,
                          vector<reco::CaloCluster> &PFClust,
                          std::vector<DetId>& MatchedRH
                          );
  //for filling PFCluster Variables
 private:
  reco::Photon matchedPhot_;
  Handle<reco::PhotonCollection> pfPhotons_;
  Handle<reco::GsfElectronCollection> pfElectrons_;
  Handle<EcalRecHitCollection>  EBReducedRecHits_;
  Handle<EcalRecHitCollection>  EEReducedRecHits_;
  Handle<EcalRecHitCollection> ESRecHits_;
  const CaloSubdetectorGeometry* geomBar_;
  const CaloSubdetectorGeometry* geomEnd_;
  Handle<BeamSpot> beamSpotHandle_;
  reco::Photon PFPhoton_;
  reco::GsfElectron PFElectron_;
  bool matchPFReco_;
  bool isPFEle_;
  bool EleVeto_;
  bool isConv_;
  bool hasSLConv_;
  std::vector<reco::CaloCluster>PFClusters_;
  std::vector<reco::CaloCluster>PFSCFootprintClusters_;
  float EinMustache_;
  float MustacheEOut_;
  float MustacheEtOut_;
  float PFPreShower1_;
  float PFPreShower2_;
  float PFLowClusE_;
  double dEtaLowestC_;
  double dPhiLowestC_;
  double PFClPhiRMS_;
  double PFClPhiRMSMust_;  
  double PFPhoLocallyCorrE_;
};
#endif