Inheritance diagram for PATTauDiscriminantCutMultiplexer:

Classes
struct	DiscriminantCutEntry

Public Member Functions
void	beginEvent (const edm::Event &event, const edm::EventSetup &eventSetup) override

double	discriminate (const pat::TauRef &) const override

	PATTauDiscriminantCutMultiplexer (const edm::ParameterSet &pset)

	~PATTauDiscriminantCutMultiplexer () override

Public Member Functions inherited from TauDiscriminationProducerBase< TauType, TauDiscriminator >
virtual void	endEvent (edm::Event &)

void	produce (edm::Event &, const edm::EventSetup &) override

	TauDiscriminationProducerBase (const edm::ParameterSet &iConfig)

	TauDiscriminationProducerBase ()

	~TauDiscriminationProducerBase () override

Public Member Functions inherited from edm::stream::EDProducer<>
	EDProducer ()=default

bool	hasAbilityToProduceInBeginLumis () const final

bool	hasAbilityToProduceInBeginRuns () const final

bool	hasAbilityToProduceInEndLumis () const final

bool	hasAbilityToProduceInEndRuns () const final

Static Public Member Functions
static void	fillDescriptions (edm::ConfigurationDescriptions &descriptions)

Static Public Member Functions inherited from TauDiscriminationProducerBase< TauType, TauDiscriminator >
static void	fillProducerDescriptions (edm::ParameterSetDescription &desc)

Private Types
typedef std::map< int, std::unique_ptr< DiscriminantCutEntry > >	DiscriminantCutMap

Private Attributes
DiscriminantCutMap	cuts_

edm::FileInPath	inputFileName_

bool	isInitialized_

edm::InputTag	key_

edm::EDGetTokenT< pat::PATTauDiscriminator >	key_token

edm::Handle< pat::PATTauDiscriminator >	keyHandle_

bool	loadMVAfromDB_

std::string	moduleLabel_

std::unique_ptr< const TFormula >	mvaOutput_normalization_

std::string	mvaOutputNormalizationName_

edm::InputTag	toMultiplex_

edm::EDGetTokenT< pat::PATTauDiscriminator >	toMultiplex_token

edm::Handle< pat::PATTauDiscriminator >	toMultiplexHandle_

int	verbosity_

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::stream::EDProducer<>
typedef CacheContexts< T... >	CacheTypes

typedef CacheTypes::GlobalCache	GlobalCache

typedef AbilityChecker< T... >	HasAbility

typedef CacheTypes::LuminosityBlockCache	LuminosityBlockCache

typedef LuminosityBlockContextT< LuminosityBlockCache, RunCache, GlobalCache >	LuminosityBlockContext

typedef CacheTypes::LuminosityBlockSummaryCache	LuminosityBlockSummaryCache

typedef CacheTypes::RunCache	RunCache

typedef RunContextT< RunCache, GlobalCache >	RunContext

typedef CacheTypes::RunSummaryCache	RunSummaryCache

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

double	prediscriminantFailValue_

edm::EDGetTokenT< TauCollection >	Tau_token

size_t	tauIndex_

edm::InputTag	TauProducer_

Detailed Description

Definition at line 39 of file PATTauDiscriminantCutMultiplexer.cc.

Member Typedef Documentation

typedef std::map<int, std::unique_ptr<DiscriminantCutEntry> > PATTauDiscriminantCutMultiplexer::DiscriminantCutMap

private

Definition at line 65 of file PATTauDiscriminantCutMultiplexer.cc.

Constructor & Destructor Documentation

PATTauDiscriminantCutMultiplexer::PATTauDiscriminantCutMultiplexer ( const edm::ParameterSet & pset )

explicit

Definition at line 135 of file PATTauDiscriminantCutMultiplexer.cc.

References taus_updatedMVAIds_cff::category, gather_cfg::cout, TkAlMuonSelectors_cfi::cut, cuts_, Exception, edm::ParameterSet::getParameter(), inputFileName_, key_, key_token, PATTauDiscriminantCutMultiplexer::DiscriminantCutEntry::kFixedCut, PATTauDiscriminantCutMultiplexer::DiscriminantCutEntry::kVariableCut, loadMVAfromDB_, taus_updatedMVAIds_cff::mapping, moduleLabel_, eostools::move(), mvaOutputNormalizationName_, AlCaHLTBitMon_QueryRunRegistry::string, toMultiplex_, toMultiplex_token, verbosity_, and HistogramManager_cfi::VPSet().

     : PATTauDiscriminationProducerBase(cfg),
       moduleLabel_(cfg.getParameter<std::string>("@module_label")),
       mvaOutput_normalization_(),
       isInitialized_(false) {
   toMultiplex_ = cfg.getParameter<edm::InputTag>("toMultiplex");
   toMultiplex_token = consumes<pat::PATTauDiscriminator>(toMultiplex_);
   key_ = cfg.getParameter<edm::InputTag>("key");
   key_token = consumes<pat::PATTauDiscriminator>(key_);
 
   verbosity_ = cfg.getParameter<int>("verbosity");
   loadMVAfromDB_ = cfg.getParameter<bool>("loadMVAfromDB");
   if (!loadMVAfromDB_) {
     inputFileName_ = cfg.getParameter<edm::FileInPath>("inputFileName");
   }
   if (verbosity_)
     std::cout << moduleLabel_ << " loadMVA = " << loadMVAfromDB_ << std::endl;
   mvaOutputNormalizationName_ = cfg.getParameter<std::string>(
       "mvaOutput_normalization");  // default value is "" which should just overwrite existing value with same empty string
 
   // Setup our cut map
   typedef std::vector<edm::ParameterSet> VPSet;
   VPSet mapping = cfg.getParameter<VPSet>("mapping");
   for (VPSet::const_iterator mappingEntry = mapping.begin(); mappingEntry != mapping.end(); ++mappingEntry) {
     unsigned category = mappingEntry->getParameter<uint32_t>("category");
     std::unique_ptr<DiscriminantCutEntry> cut{new DiscriminantCutEntry()};
     if (mappingEntry->existsAs<double>("cut")) {
       cut->cutValue_ = mappingEntry->getParameter<double>("cut");
       cut->mode_ = DiscriminantCutEntry::kFixedCut;
     } else if (mappingEntry->existsAs<std::string>("cut")) {
       cut->cutName_ = mappingEntry->getParameter<std::string>("cut");
       std::string cutVariable_string = mappingEntry->getParameter<std::string>("variable");
       cut->cutVariable_.reset(new StringObjectFunction<pat::Tau>(cutVariable_string));
       cut->mode_ = DiscriminantCutEntry::kVariableCut;
     } else {
       throw cms::Exception("PATTauDiscriminantCutMultiplexer") << " Undefined Configuration Parameter 'cut' !!\n";
     }
     cuts_[category] = std::move(cut);
   }
   verbosity_ = cfg.getParameter<int>("verbosity");
   if (verbosity_)
     std::cout << "constructed " << moduleLabel_ << std::endl;
 }

PATTauDiscriminantCutMultiplexer::~PATTauDiscriminantCutMultiplexer ( )

override

Definition at line 179 of file PATTauDiscriminantCutMultiplexer.cc.

179 {}

Member Function Documentation

void PATTauDiscriminantCutMultiplexer::beginEvent	(	const edm::Event &	event,
		const edm::EventSetup &	eventSetup
	)

overridevirtual

Reimplemented from TauDiscriminationProducerBase< TauType, TauDiscriminator >.

Definition at line 181 of file PATTauDiscriminantCutMultiplexer.cc.

References gather_cfg::cout, TkAlMuonSelectors_cfi::cut, cuts_, edm::Event::getByToken(), dtResolutionTest_cfi::inputFile, inputFileName_, isInitialized_, key_token, keyHandle_, PATTauDiscriminantCutMultiplexer::DiscriminantCutEntry::kVariableCut, loadMVAfromDB_, moduleLabel_, eostools::move(), mvaOutput_normalization_, mvaOutputNormalizationName_, groupFilesInBlocks::temp, toMultiplex_token, toMultiplexHandle_, and verbosity_.

                                                                                               {
   if (verbosity_)
     std::cout << " begin! " << moduleLabel_ << " " << isInitialized_ << std::endl;
   if (!isInitialized_) {
     //Only open the file once and we can close it when this routine is done
     // since all objects gotten from the file will have been copied
     std::unique_ptr<TFile> inputFile;
     if (!mvaOutputNormalizationName_.empty()) {
       if (!loadMVAfromDB_) {
         inputFile = openInputFile(inputFileName_);
         mvaOutput_normalization_ = loadObjectFromFile<TFormula>(*inputFile, mvaOutputNormalizationName_);
       } else {
         auto temp = loadTFormulaFromDB(
             es, mvaOutputNormalizationName_, Form("%s_mvaOutput_normalization", moduleLabel_.data()), verbosity_);
         mvaOutput_normalization_ = std::move(temp);
       }
     }
     for (DiscriminantCutMap::iterator cut = cuts_.begin(); cut != cuts_.end(); ++cut) {
       if (cut->second->mode_ == DiscriminantCutEntry::kVariableCut) {
         if (!loadMVAfromDB_) {
           if (not inputFile) {
             inputFile = openInputFile(inputFileName_);
           }
           if (verbosity_)
             std::cout << "Loading from file" << inputFileName_ << std::endl;
           cut->second->cutFunction_ = loadObjectFromFile<TGraph>(*inputFile, cut->second->cutName_);
         } else {
           if (verbosity_)
             std::cout << "Loading from DB" << std::endl;
           cut->second->cutFunction_ = loadTGraphFromDB(es, cut->second->cutName_, verbosity_);
         }
       }
     }
     isInitialized_ = true;
   }
 
   evt.getByToken(toMultiplex_token, toMultiplexHandle_);
   evt.getByToken(key_token, keyHandle_);
 }

double PATTauDiscriminantCutMultiplexer::discriminate ( const pat::TauRef & tau ) const

overridevirtual

Implements TauDiscriminationProducerBase< TauType, TauDiscriminator >.

Definition at line 221 of file PATTauDiscriminantCutMultiplexer.cc.

References gather_cfg::cout, cuts_, geometryDiff::epsilon, PATTauDiscriminantCutMultiplexer::DiscriminantCutEntry::kFixedCut, PATTauDiscriminantCutMultiplexer::DiscriminantCutEntry::kVariableCut, moduleLabel_, mvaOutput_normalization_, TauDiscriminationProducerBase< TauType, TauDiscriminator >::prediscriminantFailValue_, metsig::tau, verbosity_, multiplicitycorr_cfi::xMax, and photonAnalyzer_cfi::xMin.

                                                                                 {
   if (verbosity_) {
     std::cout << "<PATTauDiscriminantCutMultiplexer::discriminate>:" << std::endl;
     std::cout << " moduleLabel = " << moduleLabel_ << std::endl;
   }
 
   double disc_result = (*toMultiplexHandle_)[tau];
   if (verbosity_) {
     std::cout << "disc_result = " << disc_result << std::endl;
   }
   if (mvaOutput_normalization_) {
     disc_result = mvaOutput_normalization_->Eval(disc_result);
     //if ( disc_result > 1. ) disc_result = 1.;
     //if ( disc_result < 0. ) disc_result = 0.;
     if (verbosity_) {
       std::cout << "disc_result (normalized) = " << disc_result << std::endl;
     }
   }
   double key_result = (*keyHandle_)[tau];
   DiscriminantCutMap::const_iterator cutIter = cuts_.find(TMath::Nint(key_result));
 
   // Return null if it doesn't exist
   if (cutIter == cuts_.end()) {
     return prediscriminantFailValue_;
   }
   // See if the discriminator passes our cuts
   bool passesCuts = false;
   if (cutIter->second->mode_ == DiscriminantCutEntry::kFixedCut) {
     passesCuts = (disc_result > cutIter->second->cutValue_);
     if (verbosity_) {
       std::cout << "cutValue (fixed) = " << cutIter->second->cutValue_ << " --> passesCuts = " << passesCuts
                 << std::endl;
     }
   } else if (cutIter->second->mode_ == DiscriminantCutEntry::kVariableCut) {
     double cutVariable = (*cutIter->second->cutVariable_)(*tau);
     double xMin, xMax, dummy;
     cutIter->second->cutFunction_->GetPoint(0, xMin, dummy);
     cutIter->second->cutFunction_->GetPoint(cutIter->second->cutFunction_->GetN() - 1, xMax, dummy);
     const double epsilon = 1.e-3;
     if (cutVariable < (xMin + epsilon))
       cutVariable = xMin + epsilon;
     else if (cutVariable > (xMax - epsilon))
       cutVariable = xMax - epsilon;
     double cutValue = cutIter->second->cutFunction_->Eval(cutVariable);
     passesCuts = (disc_result > cutValue);
     if (verbosity_) {
       std::cout << "cutValue (@" << cutVariable << ") = " << cutValue << " --> passesCuts = " << passesCuts
                 << std::endl;
     }
   } else
     assert(0);
 
   return passesCuts;
 }

void PATTauDiscriminantCutMultiplexer::fillDescriptions ( edm::ConfigurationDescriptions & descriptions )

static

Definition at line 276 of file PATTauDiscriminantCutMultiplexer.cc.

References edm::ConfigurationDescriptions::add(), edm::ParameterSetDescription::add(), edm::ParameterSet::addParameter(), edm::ParameterSetDescription::addVPSet(), DEFINE_FWK_MODULE, TauDiscriminationProducerBase< TauType, TauDiscriminator >::fillProducerDescriptions(), HLT_2018_cff::InputTag, edm::ParameterSetDescription::setAllowAnything(), and AlCaHLTBitMon_QueryRunRegistry::string.

                                                                                                   {
   // patTauDiscriminantCutMultiplexer
   edm::ParameterSetDescription desc;
   desc.add<edm::InputTag>("toMultiplex", edm::InputTag("fixme"));
   desc.add<int>("verbosity", 0);
   {
     edm::ParameterSetDescription desc_mapping;
     desc_mapping.add<unsigned int>("category");
     desc_mapping.setAllowAnything();  //option cut can be double or (string plus additional option variable)
 
     std::vector<edm::ParameterSet> vpsd_mapping;
     edm::ParameterSet psd1;
     psd1.addParameter<unsigned int>("category", 0);
     psd1.addParameter<double>("cut", 0.5);
     vpsd_mapping.push_back(psd1);
     edm::ParameterSet psd2;
     psd2.addParameter<unsigned int>("category", 1);
     psd2.addParameter<double>("cut", 0.2);
     vpsd_mapping.push_back(psd2);
     desc.addVPSet("mapping", desc_mapping, vpsd_mapping);
   }
   desc.add<edm::FileInPath>("inputFileName", edm::FileInPath("RecoTauTag/RecoTau/data/emptyMVAinputFile"));
   desc.add<bool>("loadMVAfromDB", true);
   fillProducerDescriptions(desc);  // inherited from the base
   desc.add<std::string>("mvaOutput_normalization", "");
   desc.add<edm::InputTag>("key", edm::InputTag("fixme"));
   descriptions.add("patTauDiscriminantCutMultiplexer", desc);
 }