/data/refman/pasoursint/CMSSW_4_4_5_patch3/src/RecoParticleFlow/PFProducer/src/PhotonSelectorAlgo.cc

Go to the documentation of this file.

// 
// Original Authors: Nicholas Wardle, Florian Beaudette
//
#include "RecoParticleFlow/PFProducer/interface/PhotonSelectorAlgo.h"

PhotonSelectorAlgo::PhotonSelectorAlgo(
                                       float choice,
                                       float c_Et,
                                       float c_iso_track_a, float c_iso_track_b,
                                       float c_iso_ecal_a, float c_iso_ecal_b,
                                       float c_iso_hcal_a, float c_iso_hcal_b,
                                       float c_hoe,
                                       float comb_iso,
                                       float loose_hoe
                                       ):
  choice_(choice),
  c_Et_(c_Et),
  c_iso_track_a_(c_iso_track_a),  c_iso_track_b_(c_iso_track_b),
  c_iso_ecal_a_(c_iso_ecal_a),  c_iso_ecal_b_(c_iso_ecal_b),
  c_iso_hcal_a_(c_iso_hcal_a), c_iso_hcal_b_(c_iso_hcal_b),
  c_hoe_(c_hoe),
  comb_iso_(comb_iso),
  loose_hoe_(loose_hoe)
{
  ;
}

bool PhotonSelectorAlgo::passPhotonSelection(const reco::Photon & photon) const {

  // Photon ET
  float photonPt=photon.pt();
  if( photonPt < c_Et_ ) return false;
  if(choice_<0.1) //EGM Loose
    {
      //std::cout<<"Cuts:"<<c_Et_<<" H/E "<<c_hoe_<<"ECal Iso "<<c_iso_ecal_a_<<"HCal Iso "<<c_iso_hcal_a_<<"Track Iso "<<c_iso_track_a_<<std::endl;
      // HoE
      if (photon.hadronicOverEm() > c_hoe_ ) return false;
      
      // Track iso
      if( photon.trkSumPtHollowConeDR04() > c_iso_track_a_ + c_iso_track_b_*photonPt) return false;
      
      // ECAL iso
      if (photon.ecalRecHitSumEtConeDR04() > c_iso_ecal_a_ + c_iso_ecal_b_*photonPt) return false;
      
      // HCAL iso
      if (photon.hcalTowerSumEtConeDR04() > c_iso_hcal_a_ + c_iso_hcal_b_*photonPt) return false ;
    }
  if(choice_>0.99)
    {
      
      //std::cout<<"Cuts "<<comb_iso_<<" H/E "<<loose_hoe_<<std::endl;
      if (photon.hadronicOverEm() >loose_hoe_ ) return false;
      //Isolation variables in 0.3 cone combined
        if(photon.trkSumPtHollowConeDR03()+photon.ecalRecHitSumEtConeDR03()+photon.hcalTowerSumEtConeDR03()>comb_iso_)return false;             
    }

  return true;
}