#include <TauWPThreshold.h>

Public Member Functions
double	operator() (int dm, double pt, double eta) const

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

double	operator() (const reco::Candidate &tau) const

	TauWPThreshold (const std::string &cut_str)

Private Attributes
std::unique_ptr< TF1 >	fn_

double	value_

Detailed Description

Definition at line 13 of file TauWPThreshold.h.

Constructor & Destructor Documentation

◆ TauWPThreshold()

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

inlineexplicit

Definition at line 15 of file TauWPThreshold.h.

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

                                                       {
       bool simple_value = false;
       const char* cut_cstr = cut_str.c_str();
       char* end_cstr{};
       value_ = std::strtod(cut_cstr, &end_cstr);
       // simple number if end_cstr is empty and not equal to initial string (cut_cstr), and returned value is finite
       simple_value = !*end_cstr && cut_cstr != end_cstr && std::isfinite(value_);
       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 {};
 
         std::string fn_str = prefix;
         if (cut_str.find("return") == std::string::npos)
           fn_str += " return " + cut_str + ";}";
         else
           fn_str += 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 << "'.";
       }
     }

Member Function Documentation

◆ operator()() [1/3]

double tau::TauWPThreshold::operator()	(	int	dm,
		double	pt,
		double	eta
	)		const

inline

Definition at line 42 of file TauWPThreshold.h.

References symbols::dm, PVValHelper::eta, fn_, DiDispStaMuonMonitor_cfi::pt, and value_.

                                                            {
       if (!fn_)
         return value_;
 
       fn_->SetParameter(0, dm);
       fn_->SetParameter(1, pt);
       fn_->SetParameter(2, eta);
       return fn_->Eval(0);
     }

◆ operator()() [2/3]

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

inline

Definition at line 52 of file TauWPThreshold.h.

References boostedTaus_cff::decayMode, symbols::dm, and metsig::tau.

                                                                   {
       const int dm =
           isPFTau ? dynamic_cast<const reco::PFTau&>(tau).decayMode() : dynamic_cast<const pat::Tau&>(tau).decayMode();
       return (*this)(dm, tau.pt(), tau.eta());
     }