---

CaloRecoTauDiscriminationAgainstElectron Class Reference

#include <RecoTauTag/RecoTau/interface/CaloRecoTauDiscriminationAgainstElectron.h>

Inheritance diagram for CaloRecoTauDiscriminationAgainstElectron:

List of all members.

Public Member Functions
	CaloRecoTauDiscriminationAgainstElectron (const ParameterSet &iConfig)
virtual void	produce (Event &, const EventSetup &)
	~CaloRecoTauDiscriminationAgainstElectron ()
Private Attributes
bool	ApplyCut_leadTrackavoidsECALcrack_
bool	ApplyCut_maxleadTrackHCAL3x3hottesthitDEta_
string	CaloTauProducer_
double	leadTrack_HCAL3x3hitsEtSumOverPt_minvalue_
double	maxleadTrackHCAL3x3hottesthitDEta_

---

Detailed Description

Definition at line 29 of file CaloRecoTauDiscriminationAgainstElectron.h.

---

Constructor & Destructor Documentation

CaloRecoTauDiscriminationAgainstElectron::CaloRecoTauDiscriminationAgainstElectron

(

const ParameterSet &

iConfig

)

[inline, explicit]

Definition at line 31 of file CaloRecoTauDiscriminationAgainstElectron.h.

References ApplyCut_leadTrackavoidsECALcrack_, ApplyCut_maxleadTrackHCAL3x3hottesthitDEta_, CaloTauProducer_, edm::ParameterSet::getParameter(), leadTrack_HCAL3x3hitsEtSumOverPt_minvalue_, and maxleadTrackHCAL3x3hottesthitDEta_.

                                                                                {   
    CaloTauProducer_                            = iConfig.getParameter<string>("CaloTauProducer");
    leadTrack_HCAL3x3hitsEtSumOverPt_minvalue_  = iConfig.getParameter<double>("leadTrack_HCAL3x3hitsEtSumOverPt_minvalue");  
    ApplyCut_maxleadTrackHCAL3x3hottesthitDEta_ = iConfig.getParameter<bool>("ApplyCut_maxleadTrackHCAL3x3hottesthitDEta");
    maxleadTrackHCAL3x3hottesthitDEta_          = iConfig.getParameter<double>("maxleadTrackHCAL3x3hottesthitDEta");
    ApplyCut_leadTrackavoidsECALcrack_          = iConfig.getParameter<bool>("ApplyCut_leadTrackavoidsECALcrack");
    
    produces<CaloTauDiscriminator>();
  }

CaloRecoTauDiscriminationAgainstElectron::~CaloRecoTauDiscriminationAgainstElectron

(

)

[inline]

Definition at line 40 of file CaloRecoTauDiscriminationAgainstElectron.h.

{} 

---

Member Function Documentation

void CaloRecoTauDiscriminationAgainstElectron::produce	(	Event &	iEvent,
		const EventSetup &	iEventSetup
	)			[virtual]

Implements edm::EDProducer.

Definition at line 3 of file CaloRecoTauDiscriminationAgainstElectron.cc.

References ApplyCut_leadTrackavoidsECALcrack_, ApplyCut_maxleadTrackHCAL3x3hottesthitDEta_, CaloTauProducer_, ECALBounds::crack_absEtaIntervalA(), ECALBounds::crack_absEtaIntervalB(), ECALBounds::crack_absEtaIntervalC(), ECALBounds::crack_absEtaIntervalD(), ECALBounds::crack_absEtaIntervalE(), edm::EventSetup::get(), edm::Event::getByLabel(), edm::isnan(), leadTrack_HCAL3x3hitsEtSumOverPt_minvalue_, maxleadTrackHCAL3x3hottesthitDEta_, edm::ESHandle< T >::product(), TauTagTools::propagTrackECALSurfContactPoint(), and edm::Event::put().

                                                                                                 {
  Handle<CaloTauCollection> theCaloTauCollection;
  iEvent.getByLabel(CaloTauProducer_,theCaloTauCollection);

  // fill the AssociationVector object
  auto_ptr<CaloTauDiscriminator> theCaloTauDiscriminatorAgainstElectron(new CaloTauDiscriminator(CaloTauRefProd(theCaloTauCollection)));

  for(size_t iCaloTau=0;iCaloTau<theCaloTauCollection->size();++iCaloTau) {
    CaloTauRef theCaloTauRef(theCaloTauCollection,iCaloTau);
    if (!(*theCaloTauRef).leadTrack()){
      theCaloTauDiscriminatorAgainstElectron->setValue(iCaloTau,0);
      continue;
    }
    if (ApplyCut_maxleadTrackHCAL3x3hottesthitDEta_){
      // optional selection : ask for small |deta| between direction of propag. leading Track - ECAL inner surf. contact point and direction of highest Et hit among HCAL hits inside a 3x3 calo. tower matrix centered on direction of propag. leading Track - ECAL inner surf. contact point
      if (isnan((*theCaloTauRef).leadTrackHCAL3x3hottesthitDEta()) || (*theCaloTauRef).leadTrackHCAL3x3hottesthitDEta()>maxleadTrackHCAL3x3hottesthitDEta_){
        theCaloTauDiscriminatorAgainstElectron->setValue(iCaloTau,0);
        continue;
      }
    }
    if (ApplyCut_leadTrackavoidsECALcrack_){
      // optional selection : ask that leading track - ECAL inner surface contact point does not fall inside any ECAL eta crack 
      ESHandle<MagneticField> theMagneticField;
      iEventSetup.get<IdealMagneticFieldRecord>().get(theMagneticField);
      math::XYZPoint thepropagleadTrackECALSurfContactPoint=TauTagTools::propagTrackECALSurfContactPoint(theMagneticField.product(),(*theCaloTauRef).leadTrack());
      if(thepropagleadTrackECALSurfContactPoint.R()==0. ||
         fabs(thepropagleadTrackECALSurfContactPoint.eta())<ECALBounds::crack_absEtaIntervalA().second || 
         (fabs(thepropagleadTrackECALSurfContactPoint.eta())>ECALBounds::crack_absEtaIntervalB().first && fabs(thepropagleadTrackECALSurfContactPoint.eta())<ECALBounds::crack_absEtaIntervalB().second) ||
         (fabs(thepropagleadTrackECALSurfContactPoint.eta())>ECALBounds::crack_absEtaIntervalC().first && fabs(thepropagleadTrackECALSurfContactPoint.eta())<ECALBounds::crack_absEtaIntervalC().second) ||
         (fabs(thepropagleadTrackECALSurfContactPoint.eta())>ECALBounds::crack_absEtaIntervalD().first && fabs(thepropagleadTrackECALSurfContactPoint.eta())<ECALBounds::crack_absEtaIntervalD().second) ||
         (fabs(thepropagleadTrackECALSurfContactPoint.eta())>ECALBounds::crack_absEtaIntervalE().first && fabs(thepropagleadTrackECALSurfContactPoint.eta())<ECALBounds::crack_absEtaIntervalE().second)
         ){
        theCaloTauDiscriminatorAgainstElectron->setValue(iCaloTau,0);
        continue;
      }     
    }
    if (isnan((*theCaloTauRef).leadTrackHCAL3x3hitsEtSum())){
      theCaloTauDiscriminatorAgainstElectron->setValue(iCaloTau,0);
    }else{
      if ((*theCaloTauRef).leadTrackHCAL3x3hitsEtSum()/(*theCaloTauRef).leadTrack()->pt()<=leadTrack_HCAL3x3hitsEtSumOverPt_minvalue_) theCaloTauDiscriminatorAgainstElectron->setValue(iCaloTau,0);
      else theCaloTauDiscriminatorAgainstElectron->setValue(iCaloTau,1);
    }
  }
   
  iEvent.put(theCaloTauDiscriminatorAgainstElectron);
}

---

Member Data Documentation

bool CaloRecoTauDiscriminationAgainstElectron::ApplyCut_leadTrackavoidsECALcrack_ [private]

Definition at line 47 of file CaloRecoTauDiscriminationAgainstElectron.h.

Referenced by CaloRecoTauDiscriminationAgainstElectron(), and produce().

bool CaloRecoTauDiscriminationAgainstElectron::ApplyCut_maxleadTrackHCAL3x3hottesthitDEta_ [private]

Definition at line 45 of file CaloRecoTauDiscriminationAgainstElectron.h.

Referenced by CaloRecoTauDiscriminationAgainstElectron(), and produce().

string CaloRecoTauDiscriminationAgainstElectron::CaloTauProducer_ [private]

Definition at line 43 of file CaloRecoTauDiscriminationAgainstElectron.h.

Referenced by CaloRecoTauDiscriminationAgainstElectron(), and produce().

double CaloRecoTauDiscriminationAgainstElectron::leadTrack_HCAL3x3hitsEtSumOverPt_minvalue_ [private]

Definition at line 44 of file CaloRecoTauDiscriminationAgainstElectron.h.

Referenced by CaloRecoTauDiscriminationAgainstElectron(), and produce().

double CaloRecoTauDiscriminationAgainstElectron::maxleadTrackHCAL3x3hottesthitDEta_ [private]

Definition at line 46 of file CaloRecoTauDiscriminationAgainstElectron.h.

Referenced by CaloRecoTauDiscriminationAgainstElectron(), and produce().

---

The documentation for this class was generated from the following files:

• RecoTauTag/RecoTau/interface/CaloRecoTauDiscriminationAgainstElectron.h
• RecoTauTag/RecoTau/src/CaloRecoTauDiscriminationAgainstElectron.cc

---