Inheritance diagram for PFRecoTauDiscriminationAgainstMuon:

Public Member Functions
double	discriminate (const PFTauRef &pfTau) override

	PFRecoTauDiscriminationAgainstMuon (const edm::ParameterSet &iConfig)

	~PFRecoTauDiscriminationAgainstMuon ()

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 ()

ModuleDescription const &	moduleDescription () const

virtual	~EDProducer ()

Public Member Functions inherited from edm::ProducerBase
	ProducerBase ()

void	registerProducts (ProducerBase , ProductRegistry , ModuleDescription const &)

std::function< void(BranchDescription const &)>	registrationCallback () const
	used by the fwk to register list of products More...

virtual	~ProducerBase ()

Public Member Functions inherited from edm::EDConsumerBase
	EDConsumerBase ()

ProductHolderIndexAndSkipBit	indexFrom (EDGetToken, BranchType, TypeID const &) const

void	itemsMayGet (BranchType, std::vector< ProductHolderIndexAndSkipBit > &) const

void	itemsToGet (BranchType, std::vector< ProductHolderIndexAndSkipBit > &) const

std::vector < ProductHolderIndexAndSkipBit > const &	itemsToGetFromEvent () const

void	labelsForToken (EDGetToken iToken, Labels &oLabels) const

bool	registeredToConsume (ProductHolderIndex, bool, BranchType) const

bool	registeredToConsumeMany (TypeID const &, BranchType) const

void	updateLookup (BranchType iBranchType, ProductHolderIndexHelper const &)

virtual	~EDConsumerBase ()

Private Attributes
double	a

double	b

double	c

bool	checkNumMatches_

std::string	discriminatorOption_

double	hop_

int	maxNumberOfMatches_

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

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::ProducerBase
void	callWhenNewProductsRegistered (std::function< void(BranchDescription const &)> const &func)

Protected Member Functions inherited from edm::EDConsumerBase
template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType >	consumes (edm::InputTag const &tag)

EDGetToken	consumes (const TypeToGet &id, edm::InputTag const &tag)

template<BranchType B>
EDGetToken	consumes (TypeToGet const &id, edm::InputTag const &tag)

ConsumesCollector	consumesCollector ()
	Use a ConsumesCollector to gather consumes information from helper functions. More...

template<typename ProductType , BranchType B = InEvent>
void	consumesMany ()

void	consumesMany (const TypeToGet &id)

template<BranchType B>
void	consumesMany (const TypeToGet &id)

template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType >	mayConsume (edm::InputTag const &tag)

EDGetToken	mayConsume (const TypeToGet &id, edm::InputTag const &tag)

template<BranchType B>
EDGetToken	mayConsume (const TypeToGet &id, edm::InputTag const &tag)

Protected Attributes inherited from TauDiscriminationProducerBase< TauType, TauDiscriminator >
std::string	moduleLabel_

double	prediscriminantFailValue_

edm::EDGetTokenT< TauCollection >	Tau_token

edm::InputTag	TauProducer_

Detailed Description

Definition at line 20 of file PFRecoTauDiscriminationAgainstMuon.cc.

Constructor & Destructor Documentation

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

inlineexplicit

Definition at line 23 of file PFRecoTauDiscriminationAgainstMuon.cc.

References a, b, trackerHits::c, edm::ParameterSet::exists(), edm::ParameterSet::getParameter(), and AlCaHLTBitMon_QueryRunRegistry::string.

     : PFTauDiscriminationProducerBase(iConfig) 
   {   
     discriminatorOption_  = iConfig.getParameter<std::string>("discriminatorOption");  
     hop_  = iConfig.getParameter<double>("HoPMin");  
     a  = iConfig.getParameter<double>("a");  
     b  = iConfig.getParameter<double>("b");  
     c  = iConfig.getParameter<double>("c");      
     maxNumberOfMatches_ = iConfig.exists("maxNumberOfMatches") ? iConfig.getParameter<int>("maxNumberOfMatches") : 0;
     checkNumMatches_ = iConfig.exists("checkNumMatches") ? iConfig.getParameter<bool>("checkNumMatches") : false;
   }

PFRecoTauDiscriminationAgainstMuon::~PFRecoTauDiscriminationAgainstMuon ( )

inline

Definition at line 35 of file PFRecoTauDiscriminationAgainstMuon.cc.

35 {}

Member Function Documentation

double PFRecoTauDiscriminationAgainstMuon::discriminate ( const PFTauRef & pfTau )

override

Definition at line 49 of file PFRecoTauDiscriminationAgainstMuon.cc.

References a, muon::AllArbitrated, b, trackerHits::c, edm::hlt::Exception, muon::isGoodMuon(), edm::Ref< C, T, F >::isNonnull(), muon::segmentCompatibility(), and edm::errors::UnimplementedFeature.

 {
   bool decision = true;
 
   if ( thePFTauRef->hasMuonReference() ) {
 
     MuonRef muonref = thePFTauRef->leadPFChargedHadrCand()->muonRef();
     if ( discriminatorOption_ == "noSegMatch" ) {
       if ( muonref ->numberOfMatches() > maxNumberOfMatches_ ) decision = false;
     } else if (discriminatorOption_ == "twoDCut") {
       double seg = muon::segmentCompatibility(*muonref);
       double calo= muonref->caloCompatibility(); 
       double border = calo * a + seg * b +c;
       if ( border > 0 ) decision = false; 
     } else if ( discriminatorOption_ == "merePresence" ) {
       decision = false;
     } else if (discriminatorOption_ == "combined" ) { // testing purpose only
       unsigned int muType = 0;
       if      ( muonref->isGlobalMuon()  ) muType = 1;
       else if ( muonref->isCaloMuon()    ) muType = 2;
       else if ( muonref->isTrackerMuon() ) muType = 3;
       double muonEnergyFraction = thePFTauRef->pfTauTagInfoRef()->pfjetRef()->chargedMuEnergyFraction();
       bool eta_veto = false;
       bool phi_veto = false;
       if ( fabs(muonref->eta()) > 2.3 || (fabs(muonref->eta()) > 1.4 && fabs(muonref->eta()) < 1.6)) eta_veto = true;
       if ( muonref->phi() < 0.1 && muonref->phi() > -0.1) phi_veto = true;
       if ( muType != 1 || muonref ->numberOfMatches() > 0 || eta_veto || phi_veto || muonEnergyFraction > 0.9 ) decision = false; // as place holder
     } else if ( discriminatorOption_ == "noAllArbitrated" || discriminatorOption_ == "noAllArbitratedWithHOP" ) {
       if(checkNumMatches_ && muonref ->numberOfMatches() > maxNumberOfMatches_) decision = false;
       if ( muon::isGoodMuon(*muonref, muon::AllArbitrated) ) decision = false;      
     } else if ( discriminatorOption_ == "HOP" ) { 
       decision = true; // only calo. muon cut requested: keep all tau candidates, regardless of signals in muon system
     } else {
       throw edm::Exception(edm::errors::UnimplementedFeature) 
     << " Invalid Discriminator option = " << discriminatorOption_ << " --> please check cfi file !!\n";
     }
   } // valid muon ref
 
   // Additional calo. muon cut: veto one prongs compatible with MIP signature
   if ( discriminatorOption_ == "HOP" || discriminatorOption_ == "noAllArbitratedWithHOP" ) {
     if ( thePFTauRef->leadPFChargedHadrCand().isNonnull() ) {
       double muonCaloEn = thePFTauRef->leadPFChargedHadrCand()->hcalEnergy() + thePFTauRef->leadPFChargedHadrCand()->ecalEnergy();
       if ( thePFTauRef->decayMode() == 0 && muonCaloEn < (hop_*thePFTauRef->leadPFChargedHadrCand()->p()) ) decision = false;
     }
   }
 
   return (decision ? 1. : 0.);
 } 