|
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
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
|
#ifndef MuonPFAnalyzer_H
#define MuonPFAnalyzer_H
/** \class MuRecoAnalyzer
*
* DQM monitoring source for PF muons
*
* \author C. Battilana - CIEMAT
*/
//Base class
#include "DQMServices/Core/interface/DQMEDAnalyzer.h"
#include "FWCore/Utilities/interface/InputTag.h"
#include "DQMServices/Core/interface/DQMStore.h"
#include "FWCore/ServiceRegistry/interface/Service.h"
#include "DataFormats/MuonReco/interface/Muon.h"
#include "DataFormats/MuonReco/interface/MuonFwd.h"
#include "SimDataFormats/GeneratorProducts/interface/HepMCProduct.h"
#include "DataFormats/HepMCCandidate/interface/GenParticle.h"
#include "DataFormats/HepMCCandidate/interface/GenParticleFwd.h"
#include "DataFormats/VertexReco/interface/Vertex.h"
#include "DataFormats/VertexReco/interface/VertexFwd.h"
#include "DataFormats/BeamSpot/interface/BeamSpot.h"
#include <map>
#include <string>
class MuonPFAnalyzer : public DQMEDAnalyzer {
public:
typedef std::pair<const reco::Muon *, const reco::GenParticle *> RecoGenPair;
typedef std::vector<RecoGenPair> RecoGenCollection;
/// Constructor
explicit MuonPFAnalyzer(const edm::ParameterSet &);
/// Destructor
~MuonPFAnalyzer() override;
void analyze(const edm::Event &, const edm::EventSetup &) override;
void bookHistograms(DQMStore::IBooker &, edm::Run const &, edm::EventSetup const &) override;
private:
// Book histos for a given group of plots (e.g. for Tight TUNEP muons)
void bookHistos(DQMStore::IBooker &, const std::string &);
// Get a specific plot for a given group
MonitorElement *getPlot(const std::string &group, const std::string &type);
// Algorithm to identify muon pt track type
int muonTrackType(const reco::Muon *muon, bool usePF);
// Compute comb. rel. iso. (RECO based) for a given muon
inline float combRelIso(const reco::Muon *muon);
// Compute delta phi taking into account overflows
inline float fDeltaPhi(float phi1, float phi2);
// Set labels for code plots
void setCodeLabels(MonitorElement *plot, int nAxis);
// Fill plot within its range limits
void fillInRange(MonitorElement *plot, int nAxis, double x, double y = 0);
// Perform reco-gen geometrical matching on a best effort basis
// (if runOnMC == false or no matched gen particles are available gen is set to 0 in theRecoGen)
void recoToGenMatch(edm::Handle<reco::MuonCollection> &reco, edm::Handle<reco::GenParticleCollection> &gen);
const reco::Vertex getPrimaryVertex(edm::Handle<reco::VertexCollection> &vertex,
edm::Handle<reco::BeamSpot> &beamSpot);
edm::EDGetTokenT<reco::GenParticleCollection> theGenLabel_;
edm::EDGetTokenT<reco::MuonCollection> theRecoLabel_;
edm::EDGetTokenT<reco::VertexCollection> theVertexLabel_;
edm::EDGetTokenT<reco::BeamSpot> theBeamSpotLabel_;
std::vector<std::string> theMuonKinds;
std::map<std::string, std::map<std::string, MonitorElement *> > thePlots;
RecoGenCollection theRecoGen;
double theHighPtTh;
double theRecoGenR;
double theIsoCut;
bool theRunOnMC;
std::string theFolder;
};
#endif
|