|
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
|
#ifndef HLTrigger_Muon_Test_HLTMuonRateAnalyzer_H
#define HLTrigger_Muon_Test_HLTMuonRateAnalyzer_H
/** \class HLTMuonRateAnalyzer
* Get L1/HLT efficiency/rate plots
*
* \author J. Alcaraz
*/
// Base Class Headers
#include "FWCore/Framework/interface/one/EDAnalyzer.h"
#include "FWCore/Utilities/interface/InputTag.h"
#include "SimDataFormats/GeneratorProducts/interface/HepMCProduct.h"
#include "DataFormats/HLTReco/interface/TriggerFilterObjectWithRefs.h"
#include "DataFormats/HLTReco/interface/TriggerRefsCollections.h"
#include <vector>
class TFile;
class TH1F;
class HLTMuonRateAnalyzer : public edm::one::EDAnalyzer<> {
public:
/// Constructor
HLTMuonRateAnalyzer(const edm::ParameterSet& pset);
/// Destructor
~HLTMuonRateAnalyzer() override;
// Operations
void analyze(const edm::Event& event, const edm::EventSetup& eventSetup) override;
void beginJob() override;
void endJob() override;
private:
// Input from cfg file
edm::InputTag theGenLabel;
edm::InputTag theL1CollectionLabel;
std::vector<edm::InputTag> theHLTCollectionLabels;
edm::EDGetTokenT<edm::HepMCProduct> theGenToken;
edm::EDGetTokenT<trigger::TriggerFilterObjectWithRefs> theL1CollectionToken;
std::vector<edm::EDGetTokenT<trigger::TriggerFilterObjectWithRefs> > theHLTCollectionTokens;
double theL1ReferenceThreshold;
std::vector<double> theNSigmas;
unsigned int theNumberOfObjects;
double theCrossSection;
double theLuminosity;
double thePtMin;
double thePtMax;
unsigned int theNbins;
std::string theRootFileName;
// The output Root file
TFile* theFile;
// Histograms
TH1F* hL1eff;
TH1F* hL1rate;
std::vector<TH1F*> hHLTeff;
std::vector<TH1F*> hHLTrate;
// Counter of events (weighted in general)
double theNumberOfEvents;
};
#endif
|