ggPFPhotons Class Reference

#include <ggPFPhotons.h>

Public Member Functions
std::pair< double, double >	CalcRMS (vector< reco::CaloCluster > PFClust, reco::Photon PFPhoton)
bool	EleVeto ()
void	fillPFClusters ()
double	getPFPhoECorr (std::vector< reco::CaloCluster >PFClusters, const GBRForest *ReaderLCEB, const GBRForest *ReaderLCEE)
	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)
bool	hasSLConv ()
bool	isConv ()
bool	isPFEle ()
bool	MatchPFReco ()
float	MustE ()
float	MustEOut ()
float	MustEtOut ()
double	PFClusRMSMust ()
double	PFClusRMSTot ()
std::vector< reco::CaloCluster >	PFClusters ()
std::vector< reco::CaloCluster >	PFClustersSCFP ()
double	PFdEta ()
double	PFdPhi ()
bool	PFElectronVeto (edm::Handle< reco::ConversionCollection > &convH, edm::Handle< reco::GsfElectronCollection > &gsfElectronsHandle)
float	PFLowE ()
float	PFPS1 ()
float	PFPS2 ()
void	PhotonPFCandMatch (reco::SuperCluster sc, std::vector< reco::PFCandidatePtr > &insideBox, edm::Handle< PFCandidateCollection > &pfCandidates, vector< reco::CaloCluster > &PFClust, std::vector< DetId > &MatchedRH)
std::vector< reco::CaloCluster >	recoPhotonClusterLink (reco::Photon phot, edm::Handle< PFCandidateCollection > &pfCandidates)
std::pair< float, float >	SLPoint ()
std::pair< double, double >	SuperClusterSize (reco::Photon phot)
virtual	~ggPFPhotons ()

Static Public Member Functions
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)

Private Attributes
Handle< BeamSpot >	beamSpotHandle_
double	dEtaLowestC_
double	dPhiLowestC_
Handle< EcalRecHitCollection >	EBReducedRecHits_
Handle< EcalRecHitCollection >	EEReducedRecHits_
float	EinMustache_
bool	EleVeto_
Handle< EcalRecHitCollection >	ESRecHits_
const CaloSubdetectorGeometry *	geomBar_
const CaloSubdetectorGeometry *	geomEnd_
bool	hasSLConv_
bool	isConv_
bool	isPFEle_
reco::Photon	matchedPhot_
bool	matchPFReco_
float	MustacheEOut_
float	MustacheEtOut_
double	PFClPhiRMS_
double	PFClPhiRMSMust_
std::vector< reco::CaloCluster >	PFClusters_
reco::GsfElectron	PFElectron_
Handle < reco::GsfElectronCollection >	pfElectrons_
float	PFLowClusE_
double	PFPhoLocallyCorrE_
reco::Photon	PFPhoton_
Handle< reco::PhotonCollection >	pfPhotons_
float	PFPreShower1_
float	PFPreShower2_
std::vector< reco::CaloCluster >	PFSCFootprintClusters_

Detailed Description

Definition at line 17 of file ggPFPhotons.h.

Constructor & Destructor Documentation

ggPFPhotons::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  )

explicit

Definition at line 8 of file ggPFPhotons.cc.

References abs, deltaR(), PFRecoTauDiscriminationAgainstElectronDeadECAL_cfi::dR, isPFEle_, matchedPhot_, matchPFReco_, PFElectron_, PFPhoton_, pfPhotons_, and reco::Photon::superCluster().

              :
  matchedPhot_(phot),
  pfPhotons_(pfPhotons),
  EBReducedRecHits_(EBReducedRecHits),
  EEReducedRecHits_(EEReducedRecHits),
  ESRecHits_(ESRecHits),
  geomBar_(geomBar),
  geomEnd_(geomEnd),
  beamSpotHandle_(beamSpotHandle),
  matchPFReco_(false),
  isPFEle_(false),
  EleVeto_(false),
  isConv_(false),
  hasSLConv_(false)
{
  //do Matching:
  //std::cout<<"HERE IN NEW CLASS "<<std::endl;
  
  reco::PhotonCollection::const_iterator pfphot=pfPhotons_->begin();
  for(;pfphot!=pfPhotons_->end();++pfphot){
    if(pfphot->superCluster()==matchedPhot_.superCluster()){
      PFPhoton_= *(pfphot);
      matchPFReco_=true;
      break;
    }
  }
  reco::GsfElectronCollection::const_iterator pfele=pfElectrons->begin();
  for(;pfele!=pfElectrons->end();++pfele){
    if(pfele->superCluster().isNull())continue;
    
    if(pfele->superCluster()==matchedPhot_.superCluster()){
      if(pfele->pflowSuperCluster().isNull())continue;
      
      PFElectron_= *(pfele);
      
      matchPFReco_=true;
      isPFEle_=true;
      break;
    }
    for(PFCandidateCollection::const_iterator pfParticle =pfCandidates->begin(); pfParticle!=pfCandidates->end(); pfParticle++){
      if(abs(pfParticle->pdgId())!=11)continue;
      if(pfParticle->superClusterRef().isNull())continue;
      double dR=deltaR(pfParticle->superClusterRef()->eta(), 
               pfParticle->superClusterRef()->phi(),
               matchedPhot_.superCluster()->eta(), 
               matchedPhot_.superCluster()->phi());
      if(pfele->superCluster()==pfParticle->superClusterRef() &&dR<0.1){
    PFElectron_= *(pfele);
    matchPFReco_=true;
    isPFEle_=true;
    break;
      } 
    }  
  }

}

ggPFPhotons::~ggPFPhotons ( )

virtual

Definition at line 75 of file ggPFPhotons.cc.

{;}

Member Function Documentation

std::pair< double, double > ggPFPhotons::CalcRMS	(	vector< reco::CaloCluster >	PFClust,
		reco::Photon	PFPhoton
	)

Definition at line 198 of file ggPFPhotons.cc.

References trackerHits::c, funct::cos(), relval_parameters_module::energy, reco::LeafCandidate::eta(), eta(), phi, reco::LeafCandidate::phi(), and mathSSE::sqrt().

Referenced by fillPFClusters(), and recoPhotonClusterLink().

                                                                                                {
  double delPhi2=0;
  double delPhiSum=0;
  double delEta2=0;
  double delEtaSum=0;
  double ClusSum=0;
  float PFPhoPhi=PFPhoton.phi();
  float PFPhoEta=PFPhoton.eta();

  float PFClusPhiRMS=0;
  float PFClusEtaRMS=0;
  std::pair<double, double> RMS;  
  for(unsigned int c=0; c<PFClust.size(); ++c){
    delPhi2=(acos(cos(PFPhoPhi-PFClust[c].phi()))* acos(cos(PFPhoPhi-PFClust[c].phi())) )+delPhi2;
    delPhiSum=delPhiSum+ acos(cos(PFPhoPhi-PFClust[c].phi()))*PFClust[c].energy();
    delEta2=(PFPhoEta-PFClust[c].eta())*(PFPhoEta-PFClust[c].eta()) *PFClust[c].energy() + delEta2;
    delEtaSum=delEtaSum+((PFPhoEta-PFClust[c].eta())*PFClust[c].energy());
    ClusSum=ClusSum+PFClust[c].energy();
  }
  double meandPhi=delPhiSum/ClusSum;
  PFClusPhiRMS=sqrt(fabs(delPhi2/ClusSum - (meandPhi*meandPhi)));
  double meandEta=delEtaSum/ClusSum;
  PFClusEtaRMS=sqrt(fabs(delEta2/ClusSum - (meandEta*meandEta)));
  RMS.first=PFClusEtaRMS;
  RMS.second=PFClusPhiRMS;
  return RMS;
}

bool ggPFPhotons::EleVeto ( )

inline

Definition at line 39 of file ggPFPhotons.h.

{return EleVeto_;}

void ggPFPhotons::fillPFClusters

(

)

Definition at line 100 of file ggPFPhotons.cc.

References CalcRMS(), SiPixelRawToDigiRegional_cfi::deltaPhi, dEtaLowestC_, dPhiLowestC_, EBReducedRecHits_, EEReducedRecHits_, EinMustache_, relval_parameters_module::energy, ESRecHits_, reco::LeafCandidate::eta(), geomBar_, geomEnd_, ggPFClusters::getPFClusters(), ggPFESClusters::getPFESClusters(), ggPFClusters::getPFSuperclusterOverlap(), i, getHLTprescales::index, isPFEle_, matchedPhot_, reco::Mustache::MustacheClust(), MustacheEOut_, MustacheEtOut_, PFClPhiRMS_, PFClPhiRMSMust_, PFClusters_, PFElectron_, PFLowClusE_, reco::GsfElectron::pflowSuperCluster(), PFPhoton_, PFPreShower1_, PFPreShower2_, PFPS1(), PFPS2(), PFSCFootprintClusters_, reco::Photon::pfSuperCluster(), phi, reco::LeafCandidate::phi(), position, funct::sin(), reco::Photon::superCluster(), and theta().

Referenced by ggPFPhotonAnalyzer::analyze().

                                {
  //PFClusterCollection object with appropriate variables:
  
  ggPFClusters PFClusterCollection(EBReducedRecHits_, EEReducedRecHits_, geomBar_,   geomEnd_);
  //if(isPFEle_)cout<<"PFElectron "<<endl;
  //else cout<<"PFPHoton "<<endl;
  //fill PFClusters
  if(isPFEle_)PFClusters_=PFClusterCollection.getPFClusters(*(PFElectron_.pflowSuperCluster()));
  else PFClusters_=PFClusterCollection.getPFClusters(*(PFPhoton_.pfSuperCluster()));
  //fill PFClusters with Cluster energy from Rechits inside SC
  
  PFSCFootprintClusters_.clear();
  for(unsigned int i=0; i<PFClusters_.size(); ++i){
    float SCFootPrintE=PFClusterCollection.getPFSuperclusterOverlap(PFClusters_[i],matchedPhot_);
    CaloCluster calo(SCFootPrintE, PFClusters_[i].position());
    PFSCFootprintClusters_.push_back(calo);
  }
  //Mustache Variables:
  
  
  Mustache Must;
  std::vector<unsigned int> insideMust;
  std::vector<unsigned int> outsideMust;
  Must.MustacheClust(PFClusters_, insideMust, outsideMust);
 
  //Must.MustacheID(PFClusters_, insideMust, outsideMust);
  // cout<<"Inside "<<insideMust.size()<<", Total"<<PFClusters_.size()<<endl;
  //sum MustacheEnergy and order clusters by Energy:
  
  std::multimap<float, unsigned int>OrderedClust;
  float MustacheE=0;
  for(unsigned int i=0; i<insideMust.size(); ++i){
    unsigned int index=insideMust[i];
    OrderedClust.insert(make_pair(PFClusters_[index].energy(), index));
    MustacheE=MustacheE+PFSCFootprintClusters_[index].energy();
  }
  
  float MustacheEOut=0;
  float MustacheEtOut=0;
  for(unsigned int i=0; i<outsideMust.size(); ++i){
    unsigned int index=outsideMust[i];
    MustacheEOut=MustacheEOut+PFClusters_[index].energy();
    MustacheEtOut=MustacheEtOut+PFClusters_[index].energy()*sin(PFClusters_[index].position().theta());
  }
  MustacheEOut_=MustacheEOut;
  MustacheEtOut_=MustacheEtOut;
  EinMustache_=MustacheE;
  
  //find lowest energy Cluster
  std::multimap<float, unsigned int>::iterator it;
  it=OrderedClust.begin();
  unsigned int lowEindex=(*it).second;   
  PFLowClusE_=PFSCFootprintClusters_[lowEindex].energy();
  
  //dEta and dPhi to this Cluster:
  dEtaLowestC_=PFSCFootprintClusters_[lowEindex].eta()-PFPhoton_.eta();
  dPhiLowestC_=deltaPhi(PFSCFootprintClusters_[lowEindex].phi(),PFPhoton_.phi());
  //RMS Of All PFClusters inside SuperCluster:
  
  std::pair<double, double> RMS=CalcRMS(PFSCFootprintClusters_, PFPhoton_);
  PFClPhiRMS_=RMS.second;
  std::vector<CaloCluster>MustacheNLClust;
  for(it=OrderedClust.begin(); it!=OrderedClust.end(); ++it){
    unsigned int index=(*it).second;
    if(index==lowEindex)continue; //skip lowest cluster which could be from PU
    MustacheNLClust.push_back(PFSCFootprintClusters_[index]);
  }
  if(insideMust.size()>3){
  std::pair<double, double> RMSMust=CalcRMS(MustacheNLClust, PFPhoton_);
  PFClPhiRMSMust_=RMSMust.second;
  }
  if(insideMust.size()==2){
    MustacheNLClust.push_back(PFSCFootprintClusters_[lowEindex]);
    std::pair<double, double> RMSMust=CalcRMS(MustacheNLClust, PFPhoton_);
    PFClPhiRMSMust_=RMSMust.second;
  }
  if(insideMust.size()==1){
    PFClPhiRMSMust_=matchedPhot_.superCluster()->phiWidth();
    PFClPhiRMS_=PFClPhiRMSMust_;
  }
  
  //fill ES Clusters
  
  ggPFESClusters PFPSClusterCollection(ESRecHits_, geomEnd_);
  vector<reco::PreshowerCluster>PFPS;
  if(isPFEle_)PFPS=PFPSClusterCollection.getPFESClusters(*((PFElectron_.pflowSuperCluster())));
  else PFPS=PFPSClusterCollection.getPFESClusters(*(PFPhoton_.pfSuperCluster()));
  float PFPS1=0;
  float PFPS2=0;
  for(unsigned int i=0; i<PFPS.size(); ++i){
    if(PFPS[i].plane()==1)PFPS1=PFPS1+PFPS[i].energy();
    if(PFPS[i].plane()==2)PFPS2=PFPS2+PFPS[i].energy();
  }
  PFPreShower1_=PFPS1;
  PFPreShower2_=PFPS2;
  
}

double ggPFPhotons::getPFPhoECorr	(	std::vector< reco::CaloCluster >	PFClusters,
		const GBRForest *	ReaderLCEB,
		const GBRForest *	ReaderLCEE
	)

Definition at line 226 of file ggPFPhotons.cc.

References beamSpotHandle_, EBReducedRecHits_, EEReducedRecHits_, geomBar_, geomEnd_, ggPFClusters::getPFClusters(), i, getHLTprescales::index, isPFEle_, ggPFClusters::LocalEnergyCorrection(), reco::Mustache::MustacheClust(), PFClusters_, PFElectron_, reco::GsfElectron::pflowSuperCluster(), PFPhoLocallyCorrE_, PFPhoton_, and reco::Photon::pfSuperCluster().

Referenced by ggPFPhotonAnalyzer::analyze().

                                                                                                                                {
  TVector3 bs(beamSpotHandle_->position().x(),beamSpotHandle_->position().y(),
          beamSpotHandle_->position().z());
  float beamspotZ=bs.Z();
  //PFClusterCollection object with appropriate variables:
  ggPFClusters PFClusterCollection(EBReducedRecHits_, EEReducedRecHits_, geomBar_,   geomEnd_);
  //fill PFClusters
  PFClusters_.clear();
  if(isPFEle_)PFClusters_=PFClusterCollection.getPFClusters(*(PFElectron_.pflowSuperCluster()));
  else PFClusters_=PFClusterCollection.getPFClusters(*(PFPhoton_.pfSuperCluster()));
  Mustache Must;
  std::vector<unsigned int> insideMust;
  std::vector<unsigned int> outsideMust;
  Must.MustacheClust(PFClusters_, insideMust, outsideMust);
  float ECorr=0;
  for(unsigned int i=0; i<insideMust.size();++i){
    unsigned int index=insideMust[i];
    ECorr=ECorr+PFClusterCollection.LocalEnergyCorrection(ReaderLCEB, ReaderLCEE, PFClusters_[index], beamspotZ);
  }
  PFPhoLocallyCorrE_=ECorr;
  
  return PFPhoLocallyCorrE_;
}

bool ggPFPhotons::hasSLConv ( )

inline

Definition at line 37 of file ggPFPhotons.h.

Referenced by ggPFPhotonAnalyzer::analyze().

{return hasSLConv_;}

bool ggPFPhotons::isConv ( )

inline

Definition at line 36 of file ggPFPhotons.h.

Referenced by ggPFPhotonAnalyzer::analyze().

{ return isConv_;}

bool ggPFPhotons::isPFEle ( )

inline

Definition at line 38 of file ggPFPhotons.h.

{return isPFEle_;}

bool ggPFPhotons::MatchPFReco ( )

inline

Definition at line 35 of file ggPFPhotons.h.

Referenced by ggPFPhotonAnalyzer::analyze().

{ return matchPFReco_;}

float ggPFPhotons::MustE ( )

inline

Definition at line 42 of file ggPFPhotons.h.

Referenced by ggPFPhotonAnalyzer::analyze().

{return EinMustache_;}

float ggPFPhotons::MustEOut ( )

inline

Definition at line 43 of file ggPFPhotons.h.

Referenced by ggPFPhotonAnalyzer::analyze().

{return MustacheEOut_;}

float ggPFPhotons::MustEtOut ( )

inline

Definition at line 44 of file ggPFPhotons.h.

{return MustacheEtOut_;}

double ggPFPhotons::PFClusRMSMust ( )

inline

Definition at line 49 of file ggPFPhotons.h.

Referenced by ggPFPhotonAnalyzer::analyze().

{return PFClPhiRMSMust_;}

double ggPFPhotons::PFClusRMSTot ( )

inline

Definition at line 48 of file ggPFPhotons.h.

Referenced by ggPFPhotonAnalyzer::analyze().

{return PFClPhiRMS_;}

std::vector<reco::CaloCluster> ggPFPhotons::PFClusters ( )

inline

Definition at line 50 of file ggPFPhotons.h.

Referenced by ggPFPhotonAnalyzer::analyze(), and recoPhotonClusterLink().

{return PFClusters_;}

std::vector<reco::CaloCluster> ggPFPhotons::PFClustersSCFP ( )

inline

Definition at line 51 of file ggPFPhotons.h.

{return PFSCFootprintClusters_;}//Electron Veto:

double ggPFPhotons::PFdEta ( )

inline

Definition at line 46 of file ggPFPhotons.h.

Referenced by ggPFPhotonAnalyzer::analyze().

{return dEtaLowestC_;}

double ggPFPhotons::PFdPhi ( )

inline

Definition at line 47 of file ggPFPhotons.h.

Referenced by ggPFPhotonAnalyzer::analyze().

{return dPhiLowestC_;}

bool ggPFPhotons::PFElectronVeto	(	edm::Handle< reco::ConversionCollection > &	convH,
		edm::Handle< reco::GsfElectronCollection > &	gsfElectronsHandle
	)

Definition at line 77 of file ggPFPhotons.cc.

References beamSpotHandle_, ConversionTools::hasMatchedPromptElectron(), isPFEle_, PFElectron_, and reco::GsfElectron::superCluster().

                                                                                                                                    {
  bool isprompt=false;
  if(!isPFEle_)return isprompt;
  isprompt= ConversionTools::hasMatchedPromptElectron(PFElectron_.superCluster(), gsfElectronsHandle, convH, beamSpotHandle_->position());
  return isprompt;
}

float ggPFPhotons::PFLowE ( )

inline

Definition at line 45 of file ggPFPhotons.h.

Referenced by ggPFPhotonAnalyzer::analyze().

{return PFLowClusE_;}

float ggPFPhotons::PFPS1 ( )

inline

Definition at line 40 of file ggPFPhotons.h.

Referenced by ggPFPhotonAnalyzer::analyze(), and fillPFClusters().

{return PFPreShower1_;}

float ggPFPhotons::PFPS2 ( )

inline

Definition at line 41 of file ggPFPhotons.h.

Referenced by ggPFPhotonAnalyzer::analyze(), and fillPFClusters().

{return PFPreShower2_;}

void ggPFPhotons::PhotonPFCandMatch	(	reco::SuperCluster	sc,
		std::vector< reco::PFCandidatePtr > &	insideBox,
		edm::Handle< PFCandidateCollection > &	pfCandidates,
		vector< reco::CaloCluster > &	PFClust,
		std::vector< DetId > &	MatchedRH
	)

Definition at line 530 of file ggPFPhotons.cc.

References abs, ggPFClusters::BasicClusterPFCandLink(), trackerHits::c, reco::SuperCluster::clustersBegin(), reco::SuperCluster::clustersEnd(), filterCSVwithJSON::copy, funct::cos(), deltaR(), PFRecoTauDiscriminationAgainstElectronDeadECAL_cfi::dR, EBReducedRecHits_, EEReducedRecHits_, reco::tau::disc::Eta(), reco::CaloCluster::eta(), cropTnPTrees::frac, geomBar_, geomEnd_, i, getHLTprescales::index, reco::Photon::isEB(), reco::Photon::isEE(), matchedPhot_, benchmark_cfg::pdgId, colinearityKinematic::Phi, reco::CaloCluster::phi(), reco::CaloCluster::setEnergy(), reco::CaloCluster::setPosition(), mathSSE::sqrt(), SuperClusterSize(), x, reco::CaloCluster::x(), detailsBasic3DVector::y, reco::CaloCluster::y(), detailsBasic3DVector::z, and reco::CaloCluster::z().

Referenced by ggPFPhotonAnalyzer::analyze().

{
  std::vector<reco::PFCandidatePtr>ChgHad;
  std::vector<reco::PFCandidatePtr>PFPho;
  std::pair<double, double>scSize=SuperClusterSize(matchedPhot_);
  double etabound=scSize.first;
  double phibound=scSize.second;
  double dRbound=sqrt(etabound*etabound+phibound*phibound);
  double seedEta=sc.eta();
  double seedPhi=sc.phi();
  for(PFCandidateCollection::const_iterator pfParticle =pfCandidates->begin(); pfParticle!=pfCandidates->end(); pfParticle++){
    if(pfParticle->pdgId()==130)continue;
    unsigned int index=pfParticle - pfCandidates->begin();
    math::XYZVector photon_directionWrtVtx(sc.x() - pfParticle->vx(),
                       sc.y() - pfParticle->vy(),
                       sc.z() - pfParticle->vz());
    seedEta= photon_directionWrtVtx.Eta();
    seedPhi=photon_directionWrtVtx.Phi();
    double dphi=acos(cos(seedPhi-pfParticle->momentum().Phi()));
    double deta=fabs(seedEta-pfParticle->momentum().Eta());
    if(deta<etabound && dphi<phibound && matchedPhot_.isEB()){
      reco::PFCandidatePtr pfRef(pfCandidates,index);
      insideBox.push_back(pfRef); //Fill PFCandidates in a box around SC
      if(pfParticle->pdgId()==22)PFPho.push_back(pfRef);
      if(abs(pfParticle->pdgId())==211)ChgHad.push_back(pfRef);
    }
    double dR=deltaR(seedEta,seedPhi, pfParticle->momentum().Eta(), pfParticle->momentum().Phi());
    if(matchedPhot_.isEE() && dR<dRbound){
      //if(dR<0.4){
      reco::PFCandidatePtr pfRef(pfCandidates,index);
      insideBox.push_back(pfRef); //Fill PFCandidates in a box around SC
      if(pfParticle->pdgId()==22)PFPho.push_back(pfRef);
      if(abs(pfParticle->pdgId())==211)ChgHad.push_back(pfRef);
    }
  }
  
  ggPFClusters PFClusterCollection(EBReducedRecHits_, EEReducedRecHits_, geomBar_,   geomEnd_);
  //Link PFCandidates to Basic Clusters
  std::vector<reco::PFCandidatePtr>copy=insideBox;
  
  PFClusterCollection.BasicClusterPFCandLink(sc, insideBox, MatchedRH);
  PFPho.clear();
  ChgHad.clear();
  for(unsigned int i=0; i<insideBox.size(); ++i){
    if(insideBox[i]->pdgId()==22)PFPho.push_back(insideBox[i]);
    if(abs(insideBox[i]->pdgId())==211)ChgHad.push_back(insideBox[i]);
  }
  //finally compute frac of shared ECal Energy between PFPhotons and PFChgHad
  std::vector<unsigned int > orbitPho;
  std::vector<unsigned int > overlapChgHad;
  std::vector<float > fractionShared;
  for(unsigned int c=0; c<ChgHad.size(); ++c){
    bool satPho=false;
    float frac=1.0;
    CaloCluster calo;
    float drmin=999;
    reco::CaloCluster_iterator cit=sc.clustersBegin();
    math::XYZPointF positionvtx(ChgHad[c]->positionAtECALEntrance().x()-ChgHad[c]->vx(), 
                ChgHad[c]->positionAtECALEntrance().y()-ChgHad[c]->vy(),
                ChgHad[c]->positionAtECALEntrance().z()-ChgHad[c]->vz());
    for(; cit!=sc.clustersEnd(); ++cit){
      float dR=deltaR((*cit)->eta(), (*cit)->phi(), positionvtx.eta(), positionvtx.phi());
      if(dR<drmin){
    calo.setPosition((*cit)->position());
    drmin=dR;
      }
    }
    
    for(unsigned int ipho=0; ipho<PFPho.size(); ++ipho){
      double deta=fabs(ChgHad[c]->positionAtECALEntrance().Eta()- PFPho[ipho]->positionAtECALEntrance().Eta());
      double dphi=acos(cos(ChgHad[c]->positionAtECALEntrance().Phi()-PFPho[ipho]->positionAtECALEntrance().Phi()));
      
      if(deta<0.05 && dphi<0.07){
    orbitPho.push_back(ipho);
    frac=ChgHad[c]->ecalEnergy()/(PFPho[ipho]->ecalEnergy()+ChgHad[c]->ecalEnergy());
    fractionShared.push_back(frac);
    satPho=true;
      }
    }
    if(satPho){
      calo.setEnergy(frac*ChgHad[c]->rawEcalEnergy());
    }
    else{
      calo.setEnergy(ChgHad[c]->rawEcalEnergy());
    }
    PFClust.push_back(calo); //pushBack PFClusters from Charged Hadrons
  }
  for(unsigned int ipho=0; ipho<PFPho.size(); ++ipho){
    CaloCluster calo;
    float frac=1.0;
    math::XYZPoint caloPos(PFPho[ipho]->positionAtECALEntrance().x(), 
               PFPho[ipho]->positionAtECALEntrance().y(),
               PFPho[ipho]->positionAtECALEntrance().z());
    calo.setPosition(caloPos);
    for(unsigned int i=0; i<orbitPho.size(); ++i){
      if(ipho==orbitPho[i]){
    frac=1-fractionShared[i];
    break;
      }
    }
    calo.setEnergy(frac*PFPho[ipho]->rawEcalEnergy());
    PFClust.push_back(calo);//PFClusters from PFPhotons
  }
}

std::vector< reco::CaloCluster > ggPFPhotons::recoPhotonClusterLink	(	reco::Photon	phot,
		edm::Handle< PFCandidateCollection > &	pfCandidates
	)

Definition at line 250 of file ggPFPhotons.cc.

References CalcRMS(), funct::cos(), SiPixelRawToDigiRegional_cfi::deltaPhi, dEtaLowestC_, dPhiLowestC_, EinMustache_, relval_parameters_module::energy, reco::LeafCandidate::eta(), i, getHLTprescales::index, matchedPhot_, reco::Mustache::MustacheClust(), MustacheEOut_, MustacheEtOut_, PFClPhiRMS_, PFClPhiRMSMust_, PFClusters_, PFLowClusE_, PFPhoton_, PFPreShower1_, PFPreShower2_, phi, reco::LeafCandidate::phi(), position, funct::sin(), reco::Photon::superCluster(), SuperClusterSize(), and theta().

                                  {
  PFClusters_.clear();
  //initialize variables:
  EinMustache_=0;
  MustacheEOut_=0;
  MustacheEtOut_=0;
  PFPreShower1_=0;
  PFPreShower2_=0;
  PFLowClusE_=0;
  dEtaLowestC_=0;
  dPhiLowestC_=0;
  PFClPhiRMS_=0;
  PFClPhiRMSMust_=0;
  //find SC boundaries from RecHits:
  std::pair<double,double>Boundary=SuperClusterSize(phot);      
  double etabound=Boundary.first;
  double phibound=Boundary.second;
  double seedEta=phot.superCluster()->seed()->eta();
  double seedPhi=phot.superCluster()->seed()->phi();
  for (PFCandidateCollection::const_iterator pfParticle =pfCandidates->begin(); pfParticle!=pfCandidates->end(); pfParticle++){
    if(pfParticle->pdgId()!=22) continue; //if photon
    //use a box about superCluster Width:
    float dphi=acos(cos(seedPhi-pfParticle->phi()));
    float deta=fabs(seedEta-pfParticle->eta());
    PFPreShower1_=0;
    PFPreShower2_=0;
    if(deta<etabound && dphi<phibound){
      //hard coded size for now but will make it more dynamic
      //based on SuperCluster Shape
      math::XYZPoint position(pfParticle->positionAtECALEntrance().X(), pfParticle->positionAtECALEntrance().Y(), pfParticle->positionAtECALEntrance().Z()) ;
      CaloCluster calo(pfParticle->rawEcalEnergy() ,position );
      PFClusters_.push_back(calo);
      PFPreShower1_=PFPreShower1_+pfParticle->pS1Energy();
      PFPreShower2_=PFPreShower2_+pfParticle->pS2Energy();
    }
  }
 
    //Mustache Variables:
  
  if(PFClusters_.size()>0){
  Mustache Must;
  std::vector<unsigned int> insideMust;
  std::vector<unsigned int> outsideMust;
  Must.MustacheClust(PFClusters_, insideMust, outsideMust);
  //sum MustacheEnergy and order clusters by Energy:
  
  std::multimap<float, unsigned int>OrderedClust;
  float MustacheE=0;
  for(unsigned int i=0; i<insideMust.size(); ++i){
    unsigned int index=insideMust[i];
    OrderedClust.insert(make_pair(PFClusters_[index].energy(), index));
    MustacheE=MustacheE+PFClusters_[index].energy();
  }
  
  float MustacheEOut=0;
  float MustacheEtOut=0;
  for(unsigned int i=0; i<outsideMust.size(); ++i){
    unsigned int index=outsideMust[i];
    MustacheEOut=MustacheEOut+PFClusters_[index].energy();
    MustacheEtOut=MustacheEtOut+PFClusters_[index].energy()*sin(PFClusters_[index].position().theta());
  }
  MustacheEOut_=MustacheEOut;
  MustacheEtOut_=MustacheEtOut;
  EinMustache_=MustacheE;
  
  //find lowest energy Cluster
  std::multimap<float, unsigned int>::iterator it;
  it=OrderedClust.begin();
  unsigned int lowEindex=(*it).second;   
  PFLowClusE_=PFClusters_[lowEindex].energy();
  //dEta and dPhi to this Cluster:
  dEtaLowestC_=PFClusters_[lowEindex].eta()-PFPhoton_.eta();
  dPhiLowestC_=deltaPhi(PFClusters_[lowEindex].phi(),PFPhoton_.phi());
  //RMS Of All PFClusters inside SuperCluster:
  if(insideMust.size()==1){
    PFLowClusE_=0;
    dEtaLowestC_=0;
    dPhiLowestC_=0;
  }
  std::pair<double, double> RMS=CalcRMS(PFClusters_, PFPhoton_);
  PFClPhiRMS_=RMS.second;
  std::vector<CaloCluster>MustacheNLClust;
  for(it=OrderedClust.begin(); it!=OrderedClust.end(); ++it){
    unsigned int index=(*it).second;
    if(index==lowEindex)continue; //skip lowest cluster which could be from PU
    MustacheNLClust.push_back(PFClusters_[index]);
  }
  if(insideMust.size()>3){
  std::pair<double, double> RMSMust=CalcRMS(MustacheNLClust, PFPhoton_);
  PFClPhiRMSMust_=RMSMust.second;
  }
  if(insideMust.size()==2){
    MustacheNLClust.push_back(PFClusters_[lowEindex]);
    std::pair<double, double> RMSMust=CalcRMS(MustacheNLClust, PFPhoton_);
    PFClPhiRMSMust_=RMSMust.second;
  }
  if(insideMust.size()==1){
    PFClPhiRMSMust_=matchedPhot_.superCluster()->phiWidth();
    PFClPhiRMS_=PFClPhiRMSMust_;
  }
  
  }
  return PFClusters_;
  
}

void ggPFPhotons::recoPhotonClusterLink ( reco::SuperCluster  sc, std::vector< reco::PFCandidatePtr > &  insideMust, std::vector< reco::PFCandidatePtr > &  outsideMust, edm::Handle< PFCandidateCollection > &  pfCandidates, double  etabound, double  phibound  )

static

Definition at line 416 of file ggPFPhotons.cc.

References funct::cos(), i, getHLTprescales::index, reco::Mustache::MustacheClust(), PFClusters(), position, and reco::SuperCluster::seed().

                     {
  std::vector<reco::CaloCluster>PFClusters;
  std::vector<reco::PFCandidatePtr>PFCand;
  double seedEta=sc.seed()->eta();
  double seedPhi=sc.seed()->phi();
  for(PFCandidateCollection::const_iterator pfParticle =pfCandidates->begin(); pfParticle!=pfCandidates->end(); pfParticle++){
    unsigned int index=pfParticle - pfCandidates->begin();
    if(pfParticle->pdgId()!=22) continue; //if photon
    //use a box about superCluster Width:
    float dphi=acos(cos(seedPhi-pfParticle->phi()));
    float deta=fabs(seedEta-pfParticle->eta());
    if(deta<etabound && dphi<phibound){
      //hard coded size for now but will make it more dynamic
      //based on SuperCluster Shape
      math::XYZPoint position(pfParticle->positionAtECALEntrance().X(), pfParticle->positionAtECALEntrance().Y(), pfParticle->positionAtECALEntrance().Z()) ;
      CaloCluster calo(pfParticle->rawEcalEnergy() ,position );
      PFClusters.push_back(calo);
      reco::PFCandidatePtr pfRef(pfCandidates,index);
      PFCand.push_back(pfRef);
    }
    
  }
  Mustache Must;
  std::vector<unsigned int> insideMustindex;
  std::vector<unsigned int> outsideMustindex;
  Must.MustacheClust(PFClusters, insideMustindex, outsideMustindex);
  for(unsigned int i=0; i<insideMustindex.size(); ++i){
    unsigned int index=insideMustindex[i];
    insideMust.push_back(PFCand[index]);
  }
  for(unsigned int i=0; i<outsideMustindex.size(); ++i){
    unsigned int index=outsideMustindex[i];
    outsideMust.push_back(PFCand[index]);
  }
    
}

std::pair< float, float > ggPFPhotons::SLPoint

(

)

Definition at line 85 of file ggPFPhotons.cc.

References beamSpotHandle_, ggPFTracks::gsfElectronProj(), hasSLConv_, ggPFTracks::isConv(), isConv_, isPFEle_, PFElectron_, PFPhoton_, and ggPFTracks::SLCombZVtx().

Referenced by ggPFPhotonAnalyzer::analyze().

                                          {
  ggPFTracks pfTks(beamSpotHandle_);
  std::pair<float, float> SLPoint(0,0);
  TVector3 bs(beamSpotHandle_->position().x(),beamSpotHandle_->position().y(),
          beamSpotHandle_->position().z());
  if(isPFEle_){
    isConv_=true;
    SLPoint=pfTks.gsfElectronProj(PFElectron_);
    return SLPoint;
  }
  SLPoint=pfTks.SLCombZVtx(PFPhoton_, hasSLConv_);
  isConv_=pfTks.isConv();//bool to flag if there are conv tracks
  return SLPoint; 
}

std::pair< double, double > ggPFPhotons::SuperClusterSize

(

reco::Photon

phot

)

Definition at line 358 of file ggPFPhotons.cc.

References funct::cos(), deltaR(), dPhi(), PFRecoTauDiscriminationAgainstElectronDeadECAL_cfi::dR, EBReducedRecHits_, EEReducedRecHits_, PV3DBase< T, PVType, FrameType >::eta(), eta(), geomBar_, geomEnd_, CaloSubdetectorGeometry::getGeometry(), CaloCellGeometry::getPosition(), i, reco::Photon::isEB(), step1_ZMM_7Tev::MaxEta, ecalTB2006H4_GenSimDigiReco_cfg::MaxPhi, phi, PV3DBase< T, PVType, FrameType >::phi(), and reco::Photon::superCluster().

Referenced by PhotonPFCandMatch(), and recoPhotonClusterLink().

                                {
  std::pair<double, double>SCsize(0.1,0.4);//Eta, Phi
  //find maximum distance between SuperCluster Seed and Rec Hits
  SuperClusterRef recoSC=phot.superCluster();
  reco::CaloCluster_iterator cit=recoSC->clustersBegin();
  double seedeta=recoSC->seed()->eta();
  double seedphi=recoSC->seed()->phi();
  EcalRecHitCollection::const_iterator eb;
  EcalRecHitCollection::const_iterator ee;
  double MaxEta=-99;
  double MaxPhi=-99;
  double MaxR=-99;
  for(;cit!=recoSC->clustersEnd();++cit){
    std::vector< std::pair<DetId, float> >bcCells=(*cit)->hitsAndFractions();
    if(phot.isEB()){
      for(unsigned int i=0; i<bcCells.size(); ++i){
    for(eb=EBReducedRecHits_->begin();eb!= EBReducedRecHits_->end();++eb){
      if(eb->id().rawId()==bcCells[i].first.rawId()){
        DetId id=bcCells[i].first;
        float eta=geomBar_->getGeometry(id)->getPosition().eta();
        float dEta = fabs(seedeta-eta);
        if(dEta>MaxEta){
          MaxEta=dEta;
          float phi=geomBar_->getGeometry(id)->getPosition().phi();
          float dPhi = acos(cos( seedphi-phi));
          if(dPhi>MaxPhi){
        MaxPhi=dPhi;
          }
        }
      }
    }
    
      }
    }
    else{
      for(unsigned int i=0; i<bcCells.size(); ++i){
        for(ee=EEReducedRecHits_->begin();ee!= EEReducedRecHits_->end();++ee){
      if(ee->id().rawId()==bcCells[i].first.rawId()){
        DetId id=bcCells[i].first;
        float eta=geomEnd_->getGeometry(id)->getPosition().eta();
        float phi=geomEnd_->getGeometry(id)->getPosition().phi();
        float dR=deltaR(eta, phi, seedeta, seedphi);
        if(dR>MaxR){
          MaxR=dR;
          MaxEta=fabs(seedeta-eta);
          MaxPhi= acos(cos(seedphi-phi));
        }
      }
    }
      }
    }
  }
  SCsize.first=MaxEta; SCsize.second=MaxPhi;
  return SCsize;
}

std::pair< double, double > ggPFPhotons::SuperClusterSize ( reco::SuperCluster  sc, Handle< EcalRecHitCollection > &  EBReducedRecHits, Handle< EcalRecHitCollection > &  EEReducedRecHits, const CaloSubdetectorGeometry *  geomBar, const CaloSubdetectorGeometry *  geomEnd  )

static

Definition at line 461 of file ggPFPhotons.cc.

References reco::SuperCluster::clustersBegin(), reco::SuperCluster::clustersEnd(), funct::cos(), dPhi(), PV3DBase< T, PVType, FrameType >::eta(), reco::CaloCluster::eta(), eta(), CaloSubdetectorGeometry::getGeometry(), CaloCellGeometry::getPosition(), i, step1_ZMM_7Tev::MaxEta, ecalTB2006H4_GenSimDigiReco_cfg::MaxPhi, phi, PV3DBase< T, PVType, FrameType >::phi(), reco::CaloCluster::phi(), and DetId::subdetId().

                          {
  std::pair<double, double>SCsize(0.1,0.4);//Eta, Phi
  //find maximum distance between SuperCluster Seed and Rec Hits
  reco::CaloCluster_iterator cit=sc.clustersBegin();
  double seedeta=sc.eta();
  double seedphi=sc.phi();
  EcalRecHitCollection::const_iterator eb;
  EcalRecHitCollection::const_iterator ee;
  double MaxEta=-99;
  double MaxPhi=-99;
  for(;cit!=sc.clustersEnd();++cit){
    std::vector< std::pair<DetId, float> >bcCells=(*cit)->hitsAndFractions();
    DetId seedXtalId = bcCells[0].first ;
    int detector = seedXtalId.subdetId(); //use Seed to check if EB or EE
    bool isEb;
    if(detector==1)isEb=true;
    else isEb=false;
    if(isEb){
      for(unsigned int i=0; i<bcCells.size(); ++i){
    for(eb=EBReducedRecHits->begin();eb!= EBReducedRecHits->end();++eb){
      if(eb->id().rawId()==bcCells[i].first.rawId()){
        DetId id=bcCells[i].first;
        float eta=geomBar->getGeometry(id)->getPosition().eta();
        float dEta = fabs(seedeta-eta);
        if(dEta>MaxEta){
          MaxEta=dEta;
          float phi=geomBar->getGeometry(id)->getPosition().phi();
          float dPhi = acos(cos( seedphi-phi));
          if(dPhi>MaxPhi){
        MaxPhi=dPhi;
          }
        }
      }
    }
    
      }
    }
    else{
      for(unsigned int i=0; i<bcCells.size(); ++i){
    for(ee=EEReducedRecHits->begin();ee!= EEReducedRecHits->end();++ee){
      if(ee->id().rawId()==bcCells[i].first.rawId()){
        DetId id=bcCells[i].first;
        float eta=geomEnd->getGeometry(id)->getPosition().eta();
        float dEta = fabs(seedeta-eta);
        if(dEta>MaxEta){
          MaxEta=dEta;
        }
        float phi=geomEnd->getGeometry(id)->getPosition().phi();
        float dPhi = acos(cos(seedphi-phi));
        if(dPhi>MaxPhi){
          MaxPhi=dPhi;
        }
      }
    }
    
      }
    }
  }
  SCsize.first=MaxEta; SCsize.second=MaxPhi;
  return SCsize;
  
  
}

Member Data Documentation

Handle<BeamSpot> ggPFPhotons::beamSpotHandle_

private

Definition at line 98 of file ggPFPhotons.h.

Referenced by getPFPhoECorr(), PFElectronVeto(), and SLPoint().

double ggPFPhotons::dEtaLowestC_

private

Definition at line 114 of file ggPFPhotons.h.

Referenced by fillPFClusters(), and recoPhotonClusterLink().

double ggPFPhotons::dPhiLowestC_

private

Definition at line 115 of file ggPFPhotons.h.

Referenced by fillPFClusters(), and recoPhotonClusterLink().

Handle<EcalRecHitCollection> ggPFPhotons::EBReducedRecHits_

private

Definition at line 93 of file ggPFPhotons.h.

Referenced by fillPFClusters(), getPFPhoECorr(), PhotonPFCandMatch(), and SuperClusterSize().

Handle<EcalRecHitCollection> ggPFPhotons::EEReducedRecHits_

private

Definition at line 94 of file ggPFPhotons.h.

Referenced by fillPFClusters(), getPFPhoECorr(), PhotonPFCandMatch(), and SuperClusterSize().

float ggPFPhotons::EinMustache_

private

Definition at line 108 of file ggPFPhotons.h.

Referenced by fillPFClusters(), and recoPhotonClusterLink().

bool ggPFPhotons::EleVeto_

private

Definition at line 103 of file ggPFPhotons.h.

Handle<EcalRecHitCollection> ggPFPhotons::ESRecHits_

private

Definition at line 95 of file ggPFPhotons.h.

Referenced by fillPFClusters().

const CaloSubdetectorGeometry* ggPFPhotons::geomBar_

private

Definition at line 96 of file ggPFPhotons.h.

Referenced by fillPFClusters(), getPFPhoECorr(), PhotonPFCandMatch(), and SuperClusterSize().

const CaloSubdetectorGeometry* ggPFPhotons::geomEnd_

private

Definition at line 97 of file ggPFPhotons.h.

Referenced by fillPFClusters(), getPFPhoECorr(), PhotonPFCandMatch(), and SuperClusterSize().

bool ggPFPhotons::hasSLConv_

private

Definition at line 105 of file ggPFPhotons.h.

Referenced by SLPoint().

bool ggPFPhotons::isConv_

private

Definition at line 104 of file ggPFPhotons.h.

Referenced by SLPoint().

bool ggPFPhotons::isPFEle_

private

Definition at line 102 of file ggPFPhotons.h.

Referenced by fillPFClusters(), getPFPhoECorr(), ggPFPhotons(), PFElectronVeto(), and SLPoint().

reco::Photon ggPFPhotons::matchedPhot_

private

Definition at line 90 of file ggPFPhotons.h.

Referenced by fillPFClusters(), ggPFPhotons(), PhotonPFCandMatch(), and recoPhotonClusterLink().

bool ggPFPhotons::matchPFReco_

private

Definition at line 101 of file ggPFPhotons.h.

Referenced by ggPFPhotons().

float ggPFPhotons::MustacheEOut_

private

Definition at line 109 of file ggPFPhotons.h.

Referenced by fillPFClusters(), and recoPhotonClusterLink().

float ggPFPhotons::MustacheEtOut_

private

Definition at line 110 of file ggPFPhotons.h.

Referenced by fillPFClusters(), and recoPhotonClusterLink().

double ggPFPhotons::PFClPhiRMS_

private

Definition at line 116 of file ggPFPhotons.h.

Referenced by fillPFClusters(), and recoPhotonClusterLink().

double ggPFPhotons::PFClPhiRMSMust_

private

Definition at line 117 of file ggPFPhotons.h.

Referenced by fillPFClusters(), and recoPhotonClusterLink().

std::vector<reco::CaloCluster> ggPFPhotons::PFClusters_

private

Definition at line 106 of file ggPFPhotons.h.

Referenced by fillPFClusters(), getPFPhoECorr(), and recoPhotonClusterLink().

reco::GsfElectron ggPFPhotons::PFElectron_

private

Definition at line 100 of file ggPFPhotons.h.

Referenced by fillPFClusters(), getPFPhoECorr(), ggPFPhotons(), PFElectronVeto(), and SLPoint().

Handle<reco::GsfElectronCollection> ggPFPhotons::pfElectrons_

private

Definition at line 92 of file ggPFPhotons.h.

float ggPFPhotons::PFLowClusE_

private

Definition at line 113 of file ggPFPhotons.h.

Referenced by fillPFClusters(), and recoPhotonClusterLink().

double ggPFPhotons::PFPhoLocallyCorrE_

private

Definition at line 118 of file ggPFPhotons.h.

Referenced by getPFPhoECorr().

reco::Photon ggPFPhotons::PFPhoton_

private

Definition at line 99 of file ggPFPhotons.h.

Referenced by fillPFClusters(), getPFPhoECorr(), ggPFPhotons(), recoPhotonClusterLink(), and SLPoint().

Handle<reco::PhotonCollection> ggPFPhotons::pfPhotons_

private

Definition at line 91 of file ggPFPhotons.h.

Referenced by ggPFPhotons().

float ggPFPhotons::PFPreShower1_

private

Definition at line 111 of file ggPFPhotons.h.

Referenced by fillPFClusters(), and recoPhotonClusterLink().

float ggPFPhotons::PFPreShower2_

private

Definition at line 112 of file ggPFPhotons.h.

Referenced by fillPFClusters(), and recoPhotonClusterLink().

std::vector<reco::CaloCluster> ggPFPhotons::PFSCFootprintClusters_

private

Definition at line 107 of file ggPFPhotons.h.

Referenced by fillPFClusters().