/data/doxygen/doxygen-1.7.3/gen/CMSSW_4_2_8/src/RecoTauTag/RecoTau/plugins/CaloRecoTauDiscriminationAgainstHadronicJets.cc

Go to the documentation of this file.

/*
 * class CaloRecoTauDiscriminationAgainstHadronicJets
 * created : April 21 2010,
 * revised : ,
 * Authors : Sami Lehti (HIP)
 */

#include "RecoTauTag/RecoTau/interface/TauDiscriminationProducerBase.h"
#include "RecoTauTag/RecoTau/interface/TCTauAlgorithm.h"

using namespace reco;
using namespace edm;

class CaloRecoTauDiscriminationAgainstHadronicJets : public CaloTauDiscriminationProducerBase {
  public:
    explicit CaloRecoTauDiscriminationAgainstHadronicJets(
        const edm::ParameterSet& iConfig)
        :CaloTauDiscriminationProducerBase(iConfig){
          tcTauAlgorithm = new TCTauAlgorithm(iConfig);
        }
    ~CaloRecoTauDiscriminationAgainstHadronicJets(){}
    double discriminate(const CaloTauRef& theCaloTauRef);
    void beginEvent(const edm::Event&, const edm::EventSetup&);

  private:
    TCTauAlgorithm*  tcTauAlgorithm;
};

void CaloRecoTauDiscriminationAgainstHadronicJets::beginEvent(
    const edm::Event& iEvent, const edm::EventSetup& iSetup){
  tcTauAlgorithm->eventSetup(iEvent,iSetup);
}


double CaloRecoTauDiscriminationAgainstHadronicJets::discriminate(
    const CaloTauRef& theCaloTauRef){
  math::XYZTLorentzVector p4 =
      tcTauAlgorithm->recalculateEnergy(*theCaloTauRef);
  return ((tcTauAlgorithm->algoComponent() !=
           TCTauAlgorithm::TCAlgoHadronicJet) ? 1. : 0.);
}

DEFINE_FWK_MODULE(CaloRecoTauDiscriminationAgainstHadronicJets);