CaloRecoTauDiscriminationByIsolation.cc

Go to the documentation of this file.

/* class CaloRecoTauDiscriminationByIsolation
 * created : Jul 23 2007,
 * revised : Sep 5 2007,
 * contributors : Ludovic Houchu (Ludovic.Houchu@cern.ch ; IPHC, Strasbourg), Christian Veelken (veelken@fnal.gov ; UC Davis), Evan Friis (UC Davis)
 */

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

#include <FWCore/ParameterSet/interface/ConfigurationDescriptions.h>
#include <FWCore/ParameterSet/interface/ParameterSetDescription.h>

namespace {

using namespace reco;

class CaloRecoTauDiscriminationByIsolation final : public CaloTauDiscriminationProducerBase {
 public:
  explicit CaloRecoTauDiscriminationByIsolation(const edm::ParameterSet& iConfig):CaloTauDiscriminationProducerBase(iConfig){   
    applyDiscriminationByTrackerIsolation_ = iConfig.getParameter<bool>("ApplyDiscriminationByTrackerIsolation");
    TrackerIsolAnnulus_maximumOccupancy_   = iConfig.getParameter<unsigned>("TrackerIsolAnnulus_maximumOccupancy");   

    applyDiscriminationByECALIsolation_    = iConfig.getParameter<bool>("ApplyDiscriminationByECALIsolation");
    EcalIsolAnnulus_maximumSumEtCut_       = iConfig.getParameter<double>("ECALisolAnnulus_maximumSumEtCut");   
  }
  ~CaloRecoTauDiscriminationByIsolation() override{} 
  double discriminate(const CaloTauRef&) const override;
  static void fillDescriptions(edm::ConfigurationDescriptions & descriptions);

 private:  
  bool applyDiscriminationByTrackerIsolation_;
  unsigned TrackerIsolAnnulus_maximumOccupancy_;   
  bool applyDiscriminationByECALIsolation_;
  double EcalIsolAnnulus_maximumSumEtCut_;
};

double CaloRecoTauDiscriminationByIsolation::discriminate(const CaloTauRef& caloTau) const 
{
  if ( applyDiscriminationByTrackerIsolation_ ){  
    if ( caloTau->isolationTracks().size() > TrackerIsolAnnulus_maximumOccupancy_ ) return 0.;
  }
  
  if ( applyDiscriminationByECALIsolation_ ) {
    if ( caloTau->isolationECALhitsEtSum() > EcalIsolAnnulus_maximumSumEtCut_ ) return 0.;
  }
  
  // N.B. the lead track requirement must be included in the discriminants
  return 1.;
}

}

void
CaloRecoTauDiscriminationByIsolation::fillDescriptions(edm::ConfigurationDescriptions& descriptions) {
  // caloRecoTauDiscriminationByIsolation
  edm::ParameterSetDescription desc;
  desc.add<bool>("ApplyDiscriminationByTrackerIsolation");
  desc.add<unsigned>("TrackerIsolAnnulus_maximumOccupancy"); // unsigned means unsigned int, in my test it resulted in uint32
  desc.add<bool>("ApplyDiscriminationByECALIsolation");
  desc.add<double>("ECALisolAnnulus_maximumSumEtCut");

  fillProducerDescriptions(desc); // inherited from the base

  descriptions.add("caloRecoTauDiscriminationByIsolation", desc);
}

DEFINE_FWK_MODULE(CaloRecoTauDiscriminationByIsolation);