PhoMVACut Class Reference

Inheritance diagram for PhoMVACut:

Public Member Functions
CandidateType	candidateType () const final
void	getEventContent (const edm::EventBase &) final
result_type	operator() (const reco::PhotonPtr &) const final
	PhoMVACut (const edm::ParameterSet &c)
void	setConsumes (edm::ConsumesCollector &) final
double	value (const reco::CandidatePtr &cand) const final
Public Member Functions inherited from CutApplicatorWithEventContentBase
	CutApplicatorWithEventContentBase ()
	CutApplicatorWithEventContentBase (const edm::ParameterSet &c)
	CutApplicatorWithEventContentBase (const CutApplicatorWithEventContentBase &)=delete
CutApplicatorWithEventContentBase &	operator= (const CutApplicatorWithEventContentBase &)=delete
	~CutApplicatorWithEventContentBase () override
	Destructor. More...
Public Member Functions inherited from CutApplicatorBase
virtual result_type	asCandidate (const argument_type &) const
	CutApplicatorBase ()
	CutApplicatorBase (const edm::ParameterSet &c)
	CutApplicatorBase (const CutApplicatorBase &)=delete
const std::string &	name () const override
result_type	operator() (const argument_type &) const override final
virtual result_type	operator() (const reco::GsfElectronPtr &) const
virtual result_type	operator() (const pat::ElectronPtr &) const
virtual result_type	operator() (const reco::MuonPtr &) const
virtual result_type	operator() (const pat::MuonPtr &) const
CutApplicatorBase &	operator= (const CutApplicatorBase &)=delete
	~CutApplicatorBase () override
	Destructor. More...
Public Member Functions inherited from candidate_functions::CandidateCut
	CandidateCut ()
virtual	~CandidateCut ()

Private Attributes
edm::Handle< edm::ValueMap< int > >	_mvaCategoriesMap
const std::vector< double >	_mvaCutValues
edm::Handle< edm::ValueMap< float > >	_mvaValueMap

Additional Inherited Members
Public Types inherited from CutApplicatorBase
enum	CandidateType {   NONE, ELECTRON, MUON, PHOTON,   TAU, PATELECTRON, PATMUON, PATPHOTON,   PATTAU }
Public Types inherited from candidate_functions::CandidateCut
using	argument_type = reco::CandidatePtr
using	result_type = bool
Protected Attributes inherited from CutApplicatorWithEventContentBase
std::unordered_map< std::string, edm::InputTag >	contentTags_
std::unordered_map< std::string, edm::EDGetToken >	contentTokens_

Detailed Description

Definition at line 5 of file PhoMVACut.cc.

Constructor & Destructor Documentation

PhoMVACut::PhoMVACut

(

const edm::ParameterSet &

c

)

Definition at line 35 of file PhoMVACut.cc.

References CutApplicatorWithEventContentBase::contentTags_, and edm::ParameterSet::getParameter().

                                             :
  CutApplicatorWithEventContentBase(c),
  _mvaCutValues(c.getParameter<std::vector<double> >("mvaCuts"))
{
  edm::InputTag mvaValTag = c.getParameter<edm::InputTag>("mvaValueMapName");
  contentTags_.emplace("mvaVal",mvaValTag);
  
  edm::InputTag mvaCatTag = c.getParameter<edm::InputTag>("mvaCategoriesMapName");
  contentTags_.emplace("mvaCat",mvaCatTag);
  
}

Member Function Documentation

CandidateType PhoMVACut::candidateType ( ) const

inlinefinalvirtual

Reimplemented from CutApplicatorBase.

Definition at line 16 of file PhoMVACut.cc.

References CutApplicatorBase::PHOTON.

                                            { 
    return PHOTON; 
  }

void PhoMVACut::getEventContent ( const edm::EventBase &  ev )

finalvirtual

Implements CutApplicatorWithEventContentBase.

Definition at line 58 of file PhoMVACut.cc.

References _mvaCategoriesMap, _mvaValueMap, CutApplicatorWithEventContentBase::contentTags_, edm::EventBase::getByLabel(), and operator()().

                                                      {  

  ev.getByLabel(contentTags_["mvaVal"],_mvaValueMap);
  ev.getByLabel(contentTags_["mvaCat"],_mvaCategoriesMap);
}

CutApplicatorBase::result_type PhoMVACut::operator() ( const reco::PhotonPtr &  cand ) const

finalvirtual

Reimplemented from CutApplicatorBase.

Definition at line 66 of file PhoMVACut.cc.

References _mvaCategoriesMap, _mvaCutValues, _mvaValueMap, edm::ValueMap< T >::begin(), eostools::cat(), edm::ValueMap< T >::contains(), CutApplicatorWithEventContentBase::contentTags_, Exception, edm::Ptr< T >::id(), edm::ValueMap< T >::idSize(), edm::HandleBase::isValid(), edm::Ptr< T >::key(), AlCaHLTBitMon_QueryRunRegistry::string, pat::PATObject< ObjectType >::userFloat(), pat::PATObject< ObjectType >::userInt(), and heppy_batch::val.

Referenced by getEventContent().

                                           {  

  // in case we are by-value
  const std::string& val_name = contentTags_.find("mvaVal")->second.instance();
  const std::string& cat_name = contentTags_.find("mvaCat")->second.instance();
  edm::Ptr<pat::Photon> pat(cand);
  float val = -1.0;
  int   cat = -1;
  if( _mvaCategoriesMap.isValid() && _mvaCategoriesMap->contains( cand.id() ) &&
      _mvaValueMap.isValid() && _mvaValueMap->contains( cand.id() ) ) {
    cat = (*_mvaCategoriesMap)[cand];
    val = (*_mvaValueMap)[cand];
  } else if ( _mvaCategoriesMap.isValid() && _mvaValueMap.isValid() &&
              _mvaCategoriesMap->idSize() == 1 && _mvaValueMap->idSize() == 1 &&
              cand.id() == edm::ProductID() ) {
    // in case we have spoofed a ptr
    //note this must be a 1:1 valuemap (only one product input)
    cat = _mvaCategoriesMap->begin()[cand.key()];
    val = _mvaValueMap->begin()[cand.key()];
  } else if ( _mvaCategoriesMap.isValid() && _mvaValueMap.isValid() ){ // throw an exception
    cat = (*_mvaCategoriesMap)[cand];
    val = (*_mvaValueMap)[cand];
  }

  // Find the cut value
  const int iCategory = _mvaCategoriesMap.isValid() ? cat : pat->userInt( cat_name );
  if( iCategory >= (int)(_mvaCutValues.size()) )
    throw cms::Exception(" Error in MVA categories: ")
      << " found a particle with a category larger than max configured " << std::endl;
  const float cutValue = _mvaCutValues[iCategory];

  // Look up the MVA value for this particle
  const float mvaValue = _mvaValueMap.isValid() ? val : pat->userFloat( val_name );

  // Apply the cut and return the result
  return mvaValue > cutValue;
}

void PhoMVACut::setConsumes ( edm::ConsumesCollector &  cc )

finalvirtual

Implements CutApplicatorWithEventContentBase.

Definition at line 47 of file PhoMVACut.cc.

References edm::ConsumesCollector::consumes(), CutApplicatorWithEventContentBase::contentTags_, and CutApplicatorWithEventContentBase::contentTokens_.

                                                    {

  auto mvaVal = 
    cc.consumes<edm::ValueMap<float> >(contentTags_["mvaVal"]);
  contentTokens_.emplace("mvaVal",mvaVal);

  auto mvaCat = 
    cc.consumes<edm::ValueMap<int> >(contentTags_["mvaCat"]);
  contentTokens_.emplace("mvaCat",mvaCat);
}

double PhoMVACut::value ( const reco::CandidatePtr &  cand ) const

finalvirtual

Implements candidate_functions::CandidateCut.

Definition at line 104 of file PhoMVACut.cc.

References _mvaCategoriesMap, _mvaValueMap, edm::ValueMap< T >::begin(), edm::ValueMap< T >::contains(), CutApplicatorWithEventContentBase::contentTags_, edm::Ptr< T >::id(), edm::ValueMap< T >::idSize(), edm::HandleBase::isValid(), edm::Ptr< T >::key(), AlCaHLTBitMon_QueryRunRegistry::string, pat::PATObject< ObjectType >::userFloat(), and heppy_batch::val.

                                                          {
  
  // in case we are by-value
  const std::string& val_name =contentTags_.find("mvaVal")->second.instance();
  edm::Ptr<pat::Photon> pat(cand);
  float val = 0.0;
  if( _mvaCategoriesMap.isValid() && _mvaCategoriesMap->contains( cand.id() ) &&
      _mvaValueMap.isValid() && _mvaValueMap->contains( cand.id() ) ) {
    val = (*_mvaValueMap)[cand];
  } else if ( _mvaCategoriesMap.isValid() && _mvaValueMap.isValid() &&
              _mvaCategoriesMap->idSize() == 1 && _mvaValueMap->idSize() == 1 &&
              cand.id() == edm::ProductID() ) {
    // in case we have spoofed a ptr
    //note this must be a 1:1 valuemap (only one product input)
    val = _mvaValueMap->begin()[cand.key()];
  } else if ( _mvaCategoriesMap.isValid() && _mvaValueMap.isValid() ) {
    val = (*_mvaValueMap)[cand];
  }
  const float mvaValue = _mvaValueMap.isValid() ? val : pat->userFloat( val_name );
  return mvaValue;
}

Member Data Documentation

edm::Handle<edm::ValueMap<int> > PhoMVACut::_mvaCategoriesMap

private

Definition at line 27 of file PhoMVACut.cc.

Referenced by getEventContent(), operator()(), and value().

const std::vector<double> PhoMVACut::_mvaCutValues

private

Definition at line 23 of file PhoMVACut.cc.

Referenced by operator()().

edm::Handle<edm::ValueMap<float> > PhoMVACut::_mvaValueMap

private

Definition at line 26 of file PhoMVACut.cc.

Referenced by getEventContent(), operator()(), and value().