PFRecoTauDiscriminationAgainstElectron Class Reference

Inheritance diagram for PFRecoTauDiscriminationAgainstElectron:

Public Member Functions
double	discriminate (const PFTauRef &pfTau)
	PFRecoTauDiscriminationAgainstElectron (const edm::ParameterSet &iConfig)
	~PFRecoTauDiscriminationAgainstElectron ()
Public Member Functions inherited from TauDiscriminationProducerBase< TauType, TauDiscriminator >
virtual void	beginEvent (const edm::Event &evt, const edm::EventSetup &evtSetup)
virtual double	discriminate (const TauRef &tau)=0
virtual void	endEvent (edm::Event &evt)
void	produce (edm::Event &, const edm::EventSetup &)
	TauDiscriminationProducerBase (const edm::ParameterSet &iConfig)
	TauDiscriminationProducerBase ()
virtual	~TauDiscriminationProducerBase ()
Public Member Functions inherited from edm::EDProducer
	EDProducer ()
virtual	~EDProducer ()
Public Member Functions inherited from edm::ProducerBase
	ProducerBase ()
void	registerProducts (ProducerBase *, ProductRegistry *, ModuleDescription const &)
boost::function< void(const BranchDescription &)>	registrationCallback () const
	used by the fwk to register list of products More...
virtual	~ProducerBase ()

Private Member Functions
bool	isInEcalCrack (double) const

Private Attributes
bool	applyCut_bremCombined_
bool	applyCut_bremsRecoveryEOverPLead_
bool	applyCut_ecalCrack_
bool	applyCut_electronPreID_
bool	applyCut_electronPreID_2D_
bool	applyCut_emFraction_
bool	applyCut_EOverPLead_
bool	applyCut_hcal3x3OverPLead_
bool	applyCut_hcalMaxOverPLead_
bool	applyCut_hcalTotOverPLead_
bool	applyCut_PFElectronMVA_
double	bremCombined_fraction_
double	bremCombined_maxHOP_
double	bremCombined_minMass_
double	bremCombined_stripSize_
double	bremsRecoveryEOverPLead_maxValue_
double	bremsRecoveryEOverPLead_minValue_
double	elecPreID0_EOverPLead_maxValue_
double	elecPreID0_HOverPLead_minValue_
double	elecPreID1_EOverPLead_maxValue_
double	elecPreID1_HOverPLead_minValue_
double	emFraction_maxValue_
double	EOverPLead_maxValue_
double	EOverPLead_minValue_
double	hcal3x3OverPLead_minValue_
double	hcalMaxOverPLead_minValue_
double	hcalTotOverPLead_minValue_
double	pfelectronMVA_maxValue_
edm::InputTag	PFTauProducer_

Additional Inherited Members
Public Types inherited from TauDiscriminationProducerBase< TauType, TauDiscriminator >
typedef std::vector< TauType >	TauCollection
typedef edm::Ref< TauCollection >	TauRef
typedef edm::RefProd < TauCollection >	TauRefProd
Public Types inherited from edm::EDProducer
typedef EDProducer	ModuleType
typedef WorkerT< EDProducer >	WorkerType
Public Types inherited from edm::ProducerBase
typedef ProductRegistryHelper::TypeLabelList	TypeLabelList
Static Public Member Functions inherited from edm::EDProducer
static const std::string &	baseType ()
static void	fillDescriptions (ConfigurationDescriptions &descriptions)
static void	prevalidate (ConfigurationDescriptions &descriptions)
Protected Member Functions inherited from edm::EDProducer
CurrentProcessingContext const *	currentContext () const
Protected Member Functions inherited from edm::ProducerBase
template<class TProducer , class TMethod >
void	callWhenNewProductsRegistered (TProducer *iProd, TMethod iMethod)
Protected Attributes inherited from TauDiscriminationProducerBase< TauType, TauDiscriminator >
double	prediscriminantFailValue_
edm::InputTag	TauProducer_

Detailed Description

Definition at line 12 of file PFRecoTauDiscriminationAgainstElectron.cc.

Constructor & Destructor Documentation

PFRecoTauDiscriminationAgainstElectron::PFRecoTauDiscriminationAgainstElectron ( const edm::ParameterSet &  iConfig )

inlineexplicit

Definition at line 14 of file PFRecoTauDiscriminationAgainstElectron.cc.

References edm::ParameterSet::getParameter().

                                                                                     :PFTauDiscriminationProducerBase(iConfig) {

         emFraction_maxValue_              = iConfig.getParameter<double>("EmFraction_maxValue");
         applyCut_emFraction_              = iConfig.getParameter<bool>("ApplyCut_EmFraction");
         hcalTotOverPLead_minValue_        = iConfig.getParameter<double>("HcalTotOverPLead_minValue");
         applyCut_hcalTotOverPLead_        = iConfig.getParameter<bool>("ApplyCut_HcalTotOverPLead");
         hcalMaxOverPLead_minValue_        = iConfig.getParameter<double>("HcalMaxOverPLead_minValue");
         applyCut_hcalMaxOverPLead_        = iConfig.getParameter<bool>("ApplyCut_HcalMaxOverPLead");
         hcal3x3OverPLead_minValue_        = iConfig.getParameter<double>("Hcal3x3OverPLead_minValue");

         applyCut_hcal3x3OverPLead_        = iConfig.getParameter<bool>("ApplyCut_Hcal3x3OverPLead");
         EOverPLead_minValue_              = iConfig.getParameter<double>("EOverPLead_minValue");
         EOverPLead_maxValue_              = iConfig.getParameter<double>("EOverPLead_maxValue");
         applyCut_EOverPLead_              = iConfig.getParameter<bool>("ApplyCut_EOverPLead");
         bremsRecoveryEOverPLead_minValue_ = iConfig.getParameter<double>("BremsRecoveryEOverPLead_minValue");
         bremsRecoveryEOverPLead_maxValue_ = iConfig.getParameter<double>("BremsRecoveryEOverPLead_maxValue");

         applyCut_bremsRecoveryEOverPLead_ = iConfig.getParameter<bool>("ApplyCut_BremsRecoveryEOverPLead");

         applyCut_electronPreID_           = iConfig.getParameter<bool>("ApplyCut_ElectronPreID");

         applyCut_electronPreID_2D_        = iConfig.getParameter<bool>("ApplyCut_ElectronPreID_2D");

         elecPreID0_EOverPLead_maxValue_   = iConfig.getParameter<double>("ElecPreID0_EOverPLead_maxValue");
         elecPreID0_HOverPLead_minValue_   = iConfig.getParameter<double>("ElecPreID0_HOverPLead_minValue");
         elecPreID1_EOverPLead_maxValue_   = iConfig.getParameter<double>("ElecPreID1_EOverPLead_maxValue");
         elecPreID1_HOverPLead_minValue_   = iConfig.getParameter<double>("ElecPreID1_HOverPLead_minValue");


         applyCut_PFElectronMVA_           = iConfig.getParameter<bool>("ApplyCut_PFElectronMVA");
         pfelectronMVA_maxValue_           = iConfig.getParameter<double>("PFElectronMVA_maxValue");



         applyCut_ecalCrack_               = iConfig.getParameter<bool>("ApplyCut_EcalCrackCut");

         applyCut_bremCombined_            = iConfig.getParameter<bool>("ApplyCut_BremCombined");
         bremCombined_fraction_            = iConfig.getParameter<double>("BremCombined_Fraction");
         bremCombined_maxHOP_              = iConfig.getParameter<double>("BremCombined_HOP");
         bremCombined_minMass_             = iConfig.getParameter<double>("BremCombined_Mass");
         bremCombined_stripSize_           = iConfig.getParameter<double>("BremCombined_StripSize");

      }

PFRecoTauDiscriminationAgainstElectron::~PFRecoTauDiscriminationAgainstElectron ( )

inline

Definition at line 60 of file PFRecoTauDiscriminationAgainstElectron.cc.

{}

Member Function Documentation

double PFRecoTauDiscriminationAgainstElectron::discriminate

(

const PFTauRef &

pfTau

)

Definition at line 103 of file PFRecoTauDiscriminationAgainstElectron.cc.

References edm::Ref< C, T, F >::isNonnull(), and edm::Ref< C, T, F >::isNull().

{



    // ensure tau has at least one charged object

    if( (*thePFTauRef).leadPFChargedHadrCand().isNull() )
    {
       return 0.;
    } else
    {
       // Check if track goes to Ecal crack
       TrackRef myleadTk;
       myleadTk=(*thePFTauRef).leadPFChargedHadrCand()->trackRef();
       math::XYZPointF myleadTkEcalPos = (*thePFTauRef).leadPFChargedHadrCand()->positionAtECALEntrance();
       if(myleadTk.isNonnull())
       {
          if (applyCut_ecalCrack_ && isInEcalCrack(myleadTkEcalPos.eta()))
          {
             return 0.;
          }
       }
    }

    bool decision = false;
    bool emfPass = true, htotPass = true, hmaxPass = true;
    bool h3x3Pass = true, estripPass = true, erecovPass = true;
    bool epreidPass = true, epreid2DPass = true;
    bool mvaPass = true, bremCombinedPass = true;

    if (applyCut_emFraction_) {
      if ((*thePFTauRef).emFraction() > emFraction_maxValue_) {
    emfPass = false;
      }
    }
    if (applyCut_hcalTotOverPLead_) {
      if ((*thePFTauRef).hcalTotOverPLead() < hcalTotOverPLead_minValue_) {
    htotPass = false;
      }
    }
    if (applyCut_hcalMaxOverPLead_) {
      if ((*thePFTauRef).hcalMaxOverPLead() < hcalMaxOverPLead_minValue_) {
    hmaxPass = false;
      }
    }
    if (applyCut_hcal3x3OverPLead_) {
      if ((*thePFTauRef).hcal3x3OverPLead() < hcal3x3OverPLead_minValue_) {
    h3x3Pass = false;
      }
    }
    if (applyCut_EOverPLead_) {
      if ((*thePFTauRef).ecalStripSumEOverPLead() > EOverPLead_minValue_ &&
      (*thePFTauRef).ecalStripSumEOverPLead() < EOverPLead_maxValue_) {
    estripPass = false;
      } else {
    estripPass = true;
      }
    }
    if (applyCut_bremsRecoveryEOverPLead_) {
      if ((*thePFTauRef).bremsRecoveryEOverPLead() > bremsRecoveryEOverPLead_minValue_ &&
      (*thePFTauRef).bremsRecoveryEOverPLead() < bremsRecoveryEOverPLead_maxValue_) {
    erecovPass = false;
      } else {
    erecovPass = true;
      }
    }
    if (applyCut_electronPreID_) {
      if ((*thePFTauRef).electronPreIDDecision()) {
    epreidPass = false;
      }  else {
    epreidPass = true;
      }
    }

    if (applyCut_electronPreID_2D_) {
      if (
      ((*thePFTauRef).electronPreIDDecision() &&
       ((*thePFTauRef).ecalStripSumEOverPLead() < elecPreID1_EOverPLead_maxValue_ ||
        (*thePFTauRef).hcal3x3OverPLead() > elecPreID1_HOverPLead_minValue_))
      ||
      (!(*thePFTauRef).electronPreIDDecision() &&
       ((*thePFTauRef).ecalStripSumEOverPLead() < elecPreID0_EOverPLead_maxValue_ ||
        (*thePFTauRef).hcal3x3OverPLead() > elecPreID0_HOverPLead_minValue_))
      ){
    epreid2DPass = true;
      }  else {
    epreid2DPass = false;
      }
    }

    if (applyCut_PFElectronMVA_) {
      if ((*thePFTauRef).electronPreIDOutput()>pfelectronMVA_maxValue_) {
    mvaPass = false;
      }
    }
    if (applyCut_bremCombined_) {
      if (thePFTauRef->leadPFChargedHadrCand()->trackRef().isNull()) {
        // No KF track found
        return 0;
      }
      if(thePFTauRef->signalPFChargedHadrCands().size()==1 && thePFTauRef->signalPFGammaCands().size()==0) {
    if(thePFTauRef->leadPFChargedHadrCand()->hcalEnergy()/thePFTauRef->leadPFChargedHadrCand()->trackRef()->p()<bremCombined_maxHOP_)
      bremCombinedPass = false;
      }
      else if(thePFTauRef->signalPFChargedHadrCands().size()==1 && thePFTauRef->signalPFGammaCands().size()>0) {
    //calculate the brem ratio energy
    float bremEnergy=0.;
    float emEnergy=0.;
    for(unsigned int Nc = 0 ;Nc < thePFTauRef->signalPFGammaCands().size();++Nc)
      {
        PFCandidateRef cand = thePFTauRef->signalPFGammaCands().at(Nc);
        if(fabs(thePFTauRef->leadPFChargedHadrCand()->trackRef()->eta()-cand->eta())<bremCombined_stripSize_)
          bremEnergy+=cand->energy();
        emEnergy+=cand->energy();
      }
    if(bremEnergy/emEnergy>bremCombined_fraction_&&thePFTauRef->mass()<bremCombined_minMass_)
      bremCombinedPass = false;

      }
    }

    decision = emfPass && htotPass && hmaxPass &&
      h3x3Pass && estripPass && erecovPass && epreidPass && epreid2DPass && mvaPass &&bremCombinedPass;

    return (decision ? 1. : 0.);
}

bool PFRecoTauDiscriminationAgainstElectron::isInEcalCrack ( double  eta ) const

private

Definition at line 232 of file PFRecoTauDiscriminationAgainstElectron.cc.

{
  eta = fabs(eta);
  return (eta < 0.018 ||
      (eta>0.423 && eta<0.461) ||
      (eta>0.770 && eta<0.806) ||
      (eta>1.127 && eta<1.163) ||
      (eta>1.460 && eta<1.558));
}

Member Data Documentation

bool PFRecoTauDiscriminationAgainstElectron::applyCut_bremCombined_

private

Definition at line 93 of file PFRecoTauDiscriminationAgainstElectron.cc.

bool PFRecoTauDiscriminationAgainstElectron::applyCut_bremsRecoveryEOverPLead_

private

Definition at line 77 of file PFRecoTauDiscriminationAgainstElectron.cc.

bool PFRecoTauDiscriminationAgainstElectron::applyCut_ecalCrack_

private

Definition at line 91 of file PFRecoTauDiscriminationAgainstElectron.cc.

bool PFRecoTauDiscriminationAgainstElectron::applyCut_electronPreID_

private

Definition at line 81 of file PFRecoTauDiscriminationAgainstElectron.cc.

bool PFRecoTauDiscriminationAgainstElectron::applyCut_electronPreID_2D_

private

Definition at line 83 of file PFRecoTauDiscriminationAgainstElectron.cc.

bool PFRecoTauDiscriminationAgainstElectron::applyCut_emFraction_

private

Definition at line 65 of file PFRecoTauDiscriminationAgainstElectron.cc.

bool PFRecoTauDiscriminationAgainstElectron::applyCut_EOverPLead_

private

Definition at line 74 of file PFRecoTauDiscriminationAgainstElectron.cc.

bool PFRecoTauDiscriminationAgainstElectron::applyCut_hcal3x3OverPLead_

private

Definition at line 71 of file PFRecoTauDiscriminationAgainstElectron.cc.

bool PFRecoTauDiscriminationAgainstElectron::applyCut_hcalMaxOverPLead_

private

Definition at line 69 of file PFRecoTauDiscriminationAgainstElectron.cc.

bool PFRecoTauDiscriminationAgainstElectron::applyCut_hcalTotOverPLead_

private

Definition at line 67 of file PFRecoTauDiscriminationAgainstElectron.cc.

bool PFRecoTauDiscriminationAgainstElectron::applyCut_PFElectronMVA_

private

Definition at line 89 of file PFRecoTauDiscriminationAgainstElectron.cc.

double PFRecoTauDiscriminationAgainstElectron::bremCombined_fraction_

private

Definition at line 94 of file PFRecoTauDiscriminationAgainstElectron.cc.

double PFRecoTauDiscriminationAgainstElectron::bremCombined_maxHOP_

private

Definition at line 95 of file PFRecoTauDiscriminationAgainstElectron.cc.

double PFRecoTauDiscriminationAgainstElectron::bremCombined_minMass_

private

Definition at line 96 of file PFRecoTauDiscriminationAgainstElectron.cc.

double PFRecoTauDiscriminationAgainstElectron::bremCombined_stripSize_

private

Definition at line 97 of file PFRecoTauDiscriminationAgainstElectron.cc.

double PFRecoTauDiscriminationAgainstElectron::bremsRecoveryEOverPLead_maxValue_

private

Definition at line 79 of file PFRecoTauDiscriminationAgainstElectron.cc.

double PFRecoTauDiscriminationAgainstElectron::bremsRecoveryEOverPLead_minValue_

private

Definition at line 78 of file PFRecoTauDiscriminationAgainstElectron.cc.

double PFRecoTauDiscriminationAgainstElectron::elecPreID0_EOverPLead_maxValue_

private

Definition at line 84 of file PFRecoTauDiscriminationAgainstElectron.cc.

double PFRecoTauDiscriminationAgainstElectron::elecPreID0_HOverPLead_minValue_

private

Definition at line 85 of file PFRecoTauDiscriminationAgainstElectron.cc.

double PFRecoTauDiscriminationAgainstElectron::elecPreID1_EOverPLead_maxValue_

private

Definition at line 86 of file PFRecoTauDiscriminationAgainstElectron.cc.

double PFRecoTauDiscriminationAgainstElectron::elecPreID1_HOverPLead_minValue_

private

Definition at line 87 of file PFRecoTauDiscriminationAgainstElectron.cc.

double PFRecoTauDiscriminationAgainstElectron::emFraction_maxValue_

private

Definition at line 66 of file PFRecoTauDiscriminationAgainstElectron.cc.

double PFRecoTauDiscriminationAgainstElectron::EOverPLead_maxValue_

private

Definition at line 76 of file PFRecoTauDiscriminationAgainstElectron.cc.

double PFRecoTauDiscriminationAgainstElectron::EOverPLead_minValue_

private

Definition at line 75 of file PFRecoTauDiscriminationAgainstElectron.cc.

double PFRecoTauDiscriminationAgainstElectron::hcal3x3OverPLead_minValue_

private

Definition at line 72 of file PFRecoTauDiscriminationAgainstElectron.cc.

double PFRecoTauDiscriminationAgainstElectron::hcalMaxOverPLead_minValue_

private

Definition at line 70 of file PFRecoTauDiscriminationAgainstElectron.cc.

double PFRecoTauDiscriminationAgainstElectron::hcalTotOverPLead_minValue_

private

Definition at line 68 of file PFRecoTauDiscriminationAgainstElectron.cc.

double PFRecoTauDiscriminationAgainstElectron::pfelectronMVA_maxValue_

private

Definition at line 90 of file PFRecoTauDiscriminationAgainstElectron.cc.

edm::InputTag PFRecoTauDiscriminationAgainstElectron::PFTauProducer_

private

Definition at line 64 of file PFRecoTauDiscriminationAgainstElectron.cc.