#include <DeepTauBase.h>

Public Member Functions
double	operator() (const reco::BaseTau &tau, bool isPFTau) const

	TauWPThreshold (const std::string &cut_str)

Private Attributes
std::unique_ptr< TF1 >	fn_

double	value_

Detailed Description

Definition at line 38 of file DeepTauBase.h.

Constructor & Destructor Documentation

◆ TauWPThreshold()

deep_tau::TauWPThreshold::TauWPThreshold ( const std::string & cut_str )

explicit

Definition at line 17 of file DeepTauBase.cc.

                                                          {
     bool simple_value = false;
     try {
       size_t pos = 0;
       value_ = std::stod(cut_str, &pos);
       simple_value = (pos == cut_str.size());
     } catch (std::invalid_argument&) {
     } catch (std::out_of_range&) {
     }
     if (!simple_value) {
       static const std::string prefix =
           "[&](double *x, double *p) { const int decayMode = p[0];"
           "const double pt = p[1]; const double eta = p[2];";
       static const int n_params = 3;
       static const auto handler = [](int, Bool_t, const char*, const char*) -> void {};
  
       const std::string fn_str = prefix + cut_str + "}";
       auto old_handler = SetErrorHandler(handler);
       fn_ = std::make_unique<TF1>("fn_", fn_str.c_str(), 0, 1, n_params);
       SetErrorHandler(old_handler);
       if (!fn_->IsValid())
         throw cms::Exception("TauWPThreshold: invalid formula") << "Invalid WP cut formula = '" << cut_str << "'.";
     }
   }

References fn_, createfilelist::int, hcallasereventfilter2012_cfi::prefix, AlCaHLTBitMon_QueryRunRegistry::string, and value_.

Member Function Documentation

◆ operator()()

double deep_tau::TauWPThreshold::operator()	(	const reco::BaseTau &	tau,
		bool	isPFTau
	)		const

Definition at line 42 of file DeepTauBase.cc.

                                                                               {
     if (!fn_)
       return value_;
  
     if (isPFTau)
       fn_->SetParameter(0, dynamic_cast<const reco::PFTau&>(tau).decayMode());
     else
       fn_->SetParameter(0, dynamic_cast<const pat::Tau&>(tau).decayMode());
     fn_->SetParameter(1, tau.pt());
     fn_->SetParameter(2, tau.eta());
     return fn_->Eval(0);
   }