|
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
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
|
#ifndef HiggsDQM_H
#define HiggsDQM_H
#include "FWCore/ParameterSet/interface/ParameterSet.h"
#include "FWCore/Framework/interface/ESHandle.h"
#include "FWCore/Framework/interface/LuminosityBlock.h"
#include "FWCore/Framework/interface/Run.h"
#include "FWCore/Framework/interface/Event.h"
#include "FWCore/Framework/interface/EventSetup.h"
#include "FWCore/Utilities/interface/InputTag.h"
#include "DataFormats/Candidate/interface/Candidate.h"
#include "DataFormats/Candidate/interface/CandidateFwd.h"
#include "DataFormats/EgammaCandidates/interface/GsfElectronFwd.h"
#include "DataFormats/MuonReco/interface/MuonFwd.h"
#include "DataFormats/VertexReco/interface/VertexFwd.h"
#include "DataFormats/JetReco/interface/CaloJet.h"
#include "DataFormats/METReco/interface/CaloMETCollection.h"
#include "DataFormats/METReco/interface/PFMETCollection.h"
#include "DataFormats/METReco/interface/METCollection.h"
// Trigger stuff
#include "DataFormats/Common/interface/TriggerResults.h"
#include "DataFormats/Common/interface/Handle.h"
#include "FWCore/Framework/interface/DataKeyTags.h"
#include "Geometry/CaloGeometry/interface/CaloGeometry.h"
#include "Geometry/CaloGeometry/interface/CaloCellGeometry.h"
#include "TrackingTools/TrajectoryState/interface/TrajectoryStateOnSurface.h"
#include "DataFormats/GeometryCommonDetAlgo/interface/Measurement1D.h"
#include "DataFormats/GeometryVector/interface/GlobalVector.h"
#include "DQMServices/Core/interface/DQMStore.h"
#include "DQMServices/Core/interface/DQMEDAnalyzer.h"
#include <iostream>
#include <fstream>
#include <string>
#include <vector>
#include <map>
class HiggsDQM : public DQMEDAnalyzer {
public:
HiggsDQM(const edm::ParameterSet& ps);
~HiggsDQM() override;
protected:
//Book histograms
void bookHistograms(DQMStore::IBooker&, edm::Run const&, edm::EventSetup const&) override;
void analyze(edm::Event const& e, edm::EventSetup const& eSetup) override;
private:
double Distance(const reco::Candidate& c1, const reco::Candidate& c2);
double DistancePhi(const reco::Candidate& c1, const reco::Candidate& c2);
double calcDeltaPhi(double phi1, double phi2);
void bookHistos(DQMStore* bei);
unsigned long long m_cacheID_;
int nEvents_, irun, ievt;
reco::CandidateCollection* leptonscands_;
int leptonflavor;
float pi;
// Variables from config file
std::string theElecTriggerPathToPass;
std::string theMuonTriggerPathToPass;
edm::InputTag theTriggerResultsCollectionTag_;
edm::InputTag theCaloJetCollectionLabel_;
edm::EDGetTokenT<edm::TriggerResults> theTriggerResultsCollection_;
edm::EDGetTokenT<reco::VertexCollection> vertexToken_; // offlinePrimaryVertices
edm::EDGetTokenT<reco::MuonCollection> theMuonCollectionToken_;
edm::EDGetTokenT<reco::GsfElectronCollection> theElectronCollectionToken_;
edm::EDGetTokenT<reco::CaloJetCollection> theCaloJetCollectionToken_;
edm::EDGetTokenT<reco::CaloMETCollection> theCaloMETCollectionToken_;
edm::EDGetTokenT<reco::PFMETCollection> thePfMETCollectionToken_;
double ptThrMu1_; // pt cut on the first muon for the Z^0
double ptThrMu2_; // pt cut on the second muon for the Z^0
// Histograms
MonitorElement* h_vertex_number;
MonitorElement* h_vertex_chi2;
MonitorElement* h_vertex_d0;
MonitorElement* h_vertex_numTrks;
MonitorElement* h_vertex_sumTrks;
MonitorElement* h_jet_et;
MonitorElement* h_jet2_et;
MonitorElement* h_jet_count;
MonitorElement* h_caloMet;
MonitorElement* h_caloMet_phi;
MonitorElement* h_pfMet;
MonitorElement* h_pfMet_phi;
int nfourlept, nElectron, nMuon, nLepton, nZEE, nZMuMu, nHiggs, nLooseIsolEle, nLooseIsolMu;
MonitorElement* h_eMultiplicity;
MonitorElement* h_mMultiplicity;
MonitorElement* h_ePt;
MonitorElement* h_eEta;
MonitorElement* h_ePhi;
MonitorElement* h_mPt_GMTM;
MonitorElement* h_mEta_GMTM;
MonitorElement* h_mPhi_GMTM;
MonitorElement* h_mPt_GMPT;
MonitorElement* h_mEta_GMPT;
MonitorElement* h_mPhi_GMPT;
MonitorElement* h_mPt_GM;
MonitorElement* h_mEta_GM;
MonitorElement* h_mPhi_GM;
MonitorElement* h_mPt_TM;
MonitorElement* h_mEta_TM;
MonitorElement* h_mPhi_TM;
MonitorElement* h_mPt_STAM;
MonitorElement* h_mEta_STAM;
MonitorElement* h_mPhi_STAM;
MonitorElement* h_eCombIso;
MonitorElement* h_mCombIso;
MonitorElement* h_dimumass_GMGM;
MonitorElement* h_dimumass_GMTM;
MonitorElement* h_dimumass_TMTM;
MonitorElement* h_dielemass;
MonitorElement* h_lepcounts;
};
#endif
/* Local Variables: */
/* show-trailing-whitespace: t */
/* truncate-lines: t */
/* End: */
|