Strategy

Description

Different strategies are possible to determine the trigger efficiency of jets systems, e.g. should we trigger on

• highest pt?
• most central?
• most forward?
• best match?
• any?
• all jets?
• ...

This class includes the selection and the corresponding histograms to fill to get the response and the efficiency.

Collaboration diagram for Strategy:

Public Member Functions
	Strategy (string name, const vector< double > &thresholds)
tuple< float, float, float >	GetFillValues (const auto &hlt_it, const auto &rec_it, const double threshold)
void	Fill (const RecJet &recjet, const FourVector &hltjet, const vector< bool > &Bits, double w, size_t shift=1)
void	Write (TDirectory *d)

Private Attributes
const char *	axistitle = ";p_{T}^{HLT};|#eta^{rec}|;#frac{p_{T}^{rec}}{p_{T}^{HLT bin low edge}}"
string	name
vector< double >	thresholds
unique_ptr< TEfficiency >	efficiency
unique_ptr< TH3 >	response

Constructor & Destructor Documentation

Strategy()

Strategy ( string  name, const vector< double > &  thresholds  )

inline

Extracts the appropriate values from the input, fills the efficiency and the response for each trigger and stores the histograms in the output file.

        : name(name)
        , thresholds(thresholds)
        , efficiency(makeRespHist<TEfficiency>(name + "_efficiency", thresholds, axistitle))
        , response  (makeRespHist<TH3D>       (name + "_response"  , thresholds, axistitle))
    { }

Member Function Documentation

Fill()

void Fill ( const RecJet &  recjet, const FourVector &  hltjet, const vector< bool > &  Bits, double  w, size_t  shift = 1  )

inline

    {
        for (size_t i = shift; i < Bits.size(); ++i) {
 
            // the reference trigger must have fired
            if (!Bits.at(i-1)) continue;
 
            double test_threshold = thresholds[i  ],
                   ref_threshold  = thresholds[i-shift];
 
            float recpt = recjet.CorrPt(),
                  norm_recpt = recpt/test_threshold,
                  hltpt = hltjet.Pt(),
                  receta = abs( recjet.p4.Eta() );
 
            bool pass = hltpt > test_threshold;
 
            double jw = recjet.Weight();
            efficiency->FillWeighted(pass, jw*w  , test_threshold, receta, norm_recpt);
            response  ->Fill(test_threshold, receta, norm_recpt, jw*w);
        }
    }

GetFillValues()

tuple<float, float, float> GetFillValues ( const auto &  hlt_it, const auto &  rec_it, const double  threshold  )

inline

    {
        auto ptrec = rec_it->CorrPt();
        auto pthlt = hlt_it->Pt();
 
        auto delta  = ptrec/threshold;
        auto etarec = abs(rec_it.p4.Eta());
 
        return make_tuple(pthlt, etarec, ptrec);
    }

Write()

void Write ( TDirectory *  d )

inline

    {
        d->cd();
        auto dd = d->mkdir(name.c_str());
        dd->cd();
        efficiency->SetDirectory(dd);
        efficiency->Write("efficiency");
        response->SetDirectory(dd);
        response->Write("response");
    }

Member Data Documentation

axistitle

const char* axistitle = ";p_{T}^{HLT};|#eta^{rec}|;#frac{p_{T}^{rec}}{p_{T}^{HLT bin low edge}}"

private

efficiency

unique_ptr<TEfficiency> efficiency

private

Trigger efficiency.

name

string name

private

criterion (see global description of the class)

response

unique_ptr<TH3> response

private

Actual response (different from that of Core/JEC/bin/getHLTJetResponse)

thresholds

vector<double> thresholds

private

HLT single jet trigger thresholds.

---

The documentation for this class was generated from the following file:

• /builds/cms-analysis/general/DasAnalysisSystem/Core/Installer/Core/Trigger/bin/getTriggerEfficiencyEmulation.cc