CMS 3D CMS Logo
 All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Properties Friends Macros Pages
PFMETMonitor.cc
Go to the documentation of this file.
1 #include "DataFormats/METReco/interface/MET.h"
2 #include "FWCore/ParameterSet/interface/ParameterSet.h"
3 #include "DQMOffline/PFTau/interface/Matchers.h"
4 
5 #include "DQMOffline/PFTau/interface/PFMETMonitor.h"
6 
7 #include <TROOT.h>
8 #include <TFile.h>
9 #include <TH1.h>
10 #include <TH2.h>
11 
12 #include <TProfile.h>
13 
14 //
15 // -- Constructor
16 //
17 PFMETMonitor::PFMETMonitor( Benchmark::Mode mode) :
18  Benchmark(mode),
19  candBench_(mode),
20  matchCandBench_(mode) {
21 
22  setRange( 0.0, 10e10, -10.0, 10.0, -3.14, 3.14);
23 
24  px_ = 0;
25  sumEt_ = 0;
26  delta_ex_ = 0;
27  delta_ex_VS_set_ = 0;
28  delta_set_VS_set_ = 0;
29  delta_set_Over_set_VS_set_ = 0;
30 
31  profile_delta_ex_VS_set_ = 0;
32  profile_delta_set_VS_set_ = 0;
33  profile_delta_set_Over_set_VS_set_ = 0;
34  profileRMS_delta_ex_VS_set_ = 0;
35  profileRMS_delta_set_VS_set_ = 0;
36  profileRMS_delta_set_Over_set_VS_set_ = 0;
37 
38  createMETSpecificHistos_ = false;
39  histogramBooked_ = false;
40 }
41 //
42 // -- Destructor
43 //
44 PFMETMonitor::~PFMETMonitor() {}
45 
46 //
47 // -- Set Parameters accessing them from ParameterSet
48 //
49 void PFMETMonitor::setParameters( const edm::ParameterSet & parameterSet) {
50 
51  mode_ = (Benchmark::Mode) parameterSet.getParameter<int>( "mode" );
52  createMETSpecificHistos_ = parameterSet.getParameter<bool>( "CreateMETSpecificHistos" );
53  setRange( parameterSet.getParameter<double>("ptMin"),
54  parameterSet.getParameter<double>("ptMax"),
55  parameterSet.getParameter<double>("etaMin"),
56  parameterSet.getParameter<double>("etaMax"),
57  parameterSet.getParameter<double>("phiMin"),
58  parameterSet.getParameter<double>("phiMax") );
59 
60 
61  candBench_.setParameters(mode_);
62  matchCandBench_.setParameters(mode_);
63 }
64 //
65 // -- Set Parameters
66 //
67 void PFMETMonitor::setParameters(Benchmark::Mode mode, float ptmin, float ptmax,
68  float etamin, float etamax, float phimin,
69  float phimax, bool metSpHistos) {
70  mode_ = mode;
71  createMETSpecificHistos_ = metSpHistos;
72 
73  setRange( ptmin, ptmax, etamin, etamax, phimin, phimax);
74 
75  candBench_.setParameters(mode_);
76  matchCandBench_.setParameters(mode_);
77 }
78 //
79 // -- Create histograms accessing parameters from ParameterSet
80 //
81 void PFMETMonitor::setup(const edm::ParameterSet & parameterSet) {
82  candBench_.setup(parameterSet);
83  matchCandBench_.setup(parameterSet);
84 
85  if (createMETSpecificHistos_ && !histogramBooked_) {
86 
87  edm::ParameterSet pxPS = parameterSet.getParameter<edm::ParameterSet>("DeltaPxHistoParameter");
88  edm::ParameterSet dpxPS = parameterSet.getParameter<edm::ParameterSet>("DeltaPxHistoParameter");
89  edm::ParameterSet dptPS = parameterSet.getParameter<edm::ParameterSet>("DeltaPtHistoParameter");
90  edm::ParameterSet setPS = parameterSet.getParameter<edm::ParameterSet>("SumEtHistoParameter");
91  edm::ParameterSet dsetPS = parameterSet.getParameter<edm::ParameterSet>("DeltaSumEtHistoParameter");
92  edm::ParameterSet setOvsetPS = parameterSet.getParameter<edm::ParameterSet>("DeltaSumEtOvSumEtHistoParameter");
93 
94  if (pxPS.getParameter<bool>("switchOn")) {
95  px_ = book1D("px_", "px_;p_{X} (GeV)",
96  pxPS.getParameter<int32_t>("nBin"),
97  pxPS.getParameter<double>("xMin"),
98  pxPS.getParameter<double>("xMax"));
99  }
100  if (setPS.getParameter<bool>("switchOn")) {
101  sumEt_ = book1D("sumEt_", "sumEt_;#sumE_{T}",
102  setPS.getParameter<int32_t>("nBin"),
103  setPS.getParameter<double>("xMin"),
104  setPS.getParameter<double>("xMax"));
105  }
106  if (dpxPS.getParameter<bool>("switchOn")) {
107  delta_ex_ = book1D("delta_ex_", "#DeltaME_{X}",
108  dpxPS.getParameter<int32_t>("nBin"),
109  dpxPS.getParameter<double>("xMin"),
110  dpxPS.getParameter<double>("xMax"));
111  }
112 
113  if (dpxPS.getParameter<bool>("switchOn")) {
114  delta_ex_VS_set_ = book2D("delta_ex_VS_set_", ";SE_{T, true} (GeV);#DeltaE_{X}",
115  setPS.getParameter<int32_t>("nBin"),
116  setPS.getParameter<double>("xMin"),
117  setPS.getParameter<double>("xMax"),
118  dptPS.getParameter<int32_t>("nBin"),
119  dptPS.getParameter<double>("xMin"),
120  dptPS.getParameter<double>("xMax"));
121  }
122  if (dsetPS.getParameter<bool>("switchOn")) {
123  delta_set_VS_set_ = book2D("delta_set_VS_set_", ";SE_{T, true} (GeV);#DeltaSE_{T}",
124  setPS.getParameter<int32_t>("nBin"),
125  setPS.getParameter<double>("xMin"),
126  setPS.getParameter<double>("xMax"),
127  dsetPS.getParameter<int32_t>("nBin"),
128  dsetPS.getParameter<double>("xMin"),
129  dsetPS.getParameter<double>("xMax"));
130  }
131  if (setOvsetPS.getParameter<bool>("switchOn")) {
132  delta_set_Over_set_VS_set_ = book2D("delta_set_Over_set_VS_set_", ";SE_{T, true} (GeV);#DeltaSE_{T}/SE_{T}",
133  setPS.getParameter<int32_t>("nBin"),
134  setPS.getParameter<double>("xMin"),
135  setPS.getParameter<double>("xMax"),
136  setOvsetPS.getParameter<int32_t>("nBin"),
137  setOvsetPS.getParameter<double>("xMin"),
138  setOvsetPS.getParameter<double>("xMax"));
139  }
140  // TProfile
141  if (dpxPS.getParameter<bool>("switchOn")) {
142  profile_delta_ex_VS_set_ = bookProfile("profile_delta_ex_VS_set_", ";SE_{T, true} (GeV);#DeltaE_{X}",
143  setPS.getParameter<int32_t>("nBin"),
144  setPS.getParameter<double>("xMin"),
145  setPS.getParameter<double>("xMax"),
146  dptPS.getParameter<double>("xMin"),
147  dptPS.getParameter<double>("xMax"), "" );
148  profileRMS_delta_ex_VS_set_ = bookProfile("profileRMS_delta_ex_VS_set_", ";SE_{T, true} (GeV);#DeltaE_{X}",
149  setPS.getParameter<int32_t>("nBin"),
150  setPS.getParameter<double>("xMin"),
151  setPS.getParameter<double>("xMax"),
152  dptPS.getParameter<double>("xMin"),
153  dptPS.getParameter<double>("xMax"), "s" );
154  }
155  if (dsetPS.getParameter<bool>("switchOn")) {
156  profile_delta_set_VS_set_ = bookProfile("profile_delta_set_VS_set_", ";SE_{T, true} (GeV);#DeltaSE_{T}",
157  setPS.getParameter<int32_t>("nBin"),
158  setPS.getParameter<double>("xMin"),
159  setPS.getParameter<double>("xMax"),
160  dsetPS.getParameter<double>("xMin"),
161  dsetPS.getParameter<double>("xMax"), "" );
162  profileRMS_delta_set_VS_set_ = bookProfile("profileRMS_delta_set_VS_set_", ";SE_{T, true} (GeV);#DeltaSE_{T}",
163  setPS.getParameter<int32_t>("nBin"),
164  setPS.getParameter<double>("xMin"),
165  setPS.getParameter<double>("xMax"),
166  dsetPS.getParameter<double>("xMin"),
167  dsetPS.getParameter<double>("xMax"), "s" );
168  }
169  if (setOvsetPS.getParameter<bool>("switchOn")) {
170  profile_delta_set_Over_set_VS_set_ = bookProfile("profile_delta_set_Over_set_VS_set_", ";SE_{T, true} (GeV);#DeltaSE_{T}/SE_{T}",
171  setPS.getParameter<int32_t>("nBin"),
172  setPS.getParameter<double>("xMin"),
173  setPS.getParameter<double>("xMax"),
174  setOvsetPS.getParameter<double>("xMin"),
175  setOvsetPS.getParameter<double>("xMax"), "" );
176  profileRMS_delta_set_Over_set_VS_set_ = bookProfile("profileRMS_delta_set_Over_set_VS_set_", ";SE_{T, true} (GeV);#DeltaSE_{T}/SE_{T}",
177  setPS.getParameter<int32_t>("nBin"),
178  setPS.getParameter<double>("xMin"),
179  setPS.getParameter<double>("xMax"),
180  setOvsetPS.getParameter<double>("xMin"),
181  setOvsetPS.getParameter<double>("xMax"), "s" );
182  }
183  histogramBooked_ = true;
184  }
185 }
186 //
187 // -- Create histograms using local parameters
188 //
189 void PFMETMonitor::setup() {
190  candBench_.setup();
191  matchCandBench_.setup();
192 
193  if (createMETSpecificHistos_ && !histogramBooked_) {
194 
195  PhaseSpace pxPS = PhaseSpace( 50, 0, 200);
196  PhaseSpace dpxPS = PhaseSpace( 50, -500, 500);
197  PhaseSpace setPS = PhaseSpace( 50, 0.0, 3000);
198  PhaseSpace dsetPS = PhaseSpace( 50, -1000.0, 1000);
199  PhaseSpace setOvsetPS = PhaseSpace( 100,0., 2.);
200 
201  px_ = book1D("px_", "px_;p_{X} (GeV)", pxPS.n, pxPS.m, pxPS.M);
202 
203  sumEt_ = book1D("sumEt_", "sumEt_;#sumE_{T}", setPS.n, setPS.m, setPS.M);
204 
205  delta_ex_ = book1D("delta_ex_", "#DeltaME_{X}", dpxPS.n, dpxPS.m, dpxPS.M);
206 
207  delta_ex_VS_set_ = book2D("delta_ex_VS_set_", ";SE_{T, true} (GeV);#DeltaE_{X}",
208  setPS.n, setPS.m, setPS.M,
209  dpxPS.n, dpxPS.m, dpxPS.M );
210 
211  delta_set_VS_set_ = book2D("delta_set_VS_set_", ";SE_{T, true} (GeV);#DeltaSE_{T}",
212  setPS.n, setPS.m, setPS.M,
213  dsetPS.n, dsetPS.m, dsetPS.M );
214 
215  delta_set_Over_set_VS_set_ = book2D("delta_set_Over_set_VS_set_", ";SE_{T, true} (GeV);#DeltaSE_{T}/SE_{T}",
216  setPS.n, setPS.m, setPS.M,
217  setOvsetPS.n, setOvsetPS.m, setOvsetPS.M );
218  // TProfile
219  profile_delta_ex_VS_set_ = bookProfile("profile_delta_ex_VS_set_", ";SE_{T, true} (GeV);#DeltaE_{X}",
220  setPS.n, setPS.m, setPS.M,
221  setOvsetPS.m, setOvsetPS.M, "" );
222 
223  profile_delta_set_VS_set_ = bookProfile("profile_delta_set_VS_set_", ";SE_{T, true} (GeV);#DeltaSE_{T}",
224  setPS.n, setPS.m, setPS.M,
225  setOvsetPS.m, setOvsetPS.M, "" );
226 
227 
228  profile_delta_set_Over_set_VS_set_ = bookProfile("profile_delta_set_Over_set_VS_set_", ";SE_{T, true} (GeV);#DeltaSE_{T}/SE_{T}",
229  setPS.n, setPS.m, setPS.M,
230  setOvsetPS.m, setOvsetPS.M, "" );
231  // TProfile RMS
232  profileRMS_delta_ex_VS_set_ = bookProfile("profileRMS_delta_ex_VS_set_", ";SE_{T, true} (GeV);#DeltaE_{X}",
233  setPS.n, setPS.m, setPS.M,
234  setOvsetPS.m, setOvsetPS.M, "s" );
235 
236  profileRMS_delta_set_VS_set_ = bookProfile("profileRMS_delta_set_VS_set_", ";SE_{T, true} (GeV);#DeltaSE_{T}",
237  setPS.n, setPS.m, setPS.M,
238  setOvsetPS.m, setOvsetPS.M, "s" );
239 
240 
241  profileRMS_delta_set_Over_set_VS_set_ = bookProfile("profileRMS_delta_set_Over_set_VS_set_", ";SE_{T, true} (GeV);#DeltaSE_{T}/SE_{T}",
242  setPS.n, setPS.m, setPS.M,
243  setOvsetPS.m, setOvsetPS.M, "s" );
244  histogramBooked_ = true;
245  }
246 }
247 
248 void PFMETMonitor::setDirectory(TDirectory* dir) {
249  Benchmark::setDirectory(dir);
250 
251  candBench_.setDirectory(dir);
252  matchCandBench_.setDirectory(dir);
253 }
254 
255 void PFMETMonitor::fillOne(const reco::MET& met,
256  const reco::MET& matchedMet, float& minVal, float& maxVal) {
257  /*void PFMETMonitor::fillOne(const reco::MET& met,
258  const reco::MET& matchedMet, float& minVal, float& maxVal,
259  const edm::ParameterSet & parameterSet) {*/
260 
261  candBench_.fillOne(met); //std::cout <<"\nfillone MET candBench" <<std::endl;
262  matchCandBench_.fillOne(met, matchedMet); //std::cout <<"\nfillone MET MatchCandBench done" <<std::endl;
263  //matchCandBench_.fillOne(met, matchedMet, parameterSet); std::cout <<"\nfillone MET MatchCandBench done" <<std::endl;
264 
265  if (createMETSpecificHistos_ && histogramBooked_) {
266  if( !isInRange(met.pt(), met.eta(), met.phi() ) ) return;
267 
268  if (px_) px_->Fill(met.px());
269  if (delta_ex_) {
270  delta_ex_->Fill( met.px() - matchedMet.px() );
271  delta_ex_->Fill( met.py() - matchedMet.py() );
272  }
273  if (sumEt_) sumEt_->Fill( met.sumEt() );
274 
275  if (delta_ex_VS_set_) {
276  delta_ex_VS_set_->Fill( matchedMet.sumEt(), met.px() - matchedMet.px() );
277  delta_ex_VS_set_->Fill( matchedMet.sumEt(), met.py() - matchedMet.py() );
278  profile_delta_ex_VS_set_->Fill( matchedMet.sumEt(), met.px() - matchedMet.px() );
279  profile_delta_ex_VS_set_->Fill( matchedMet.sumEt(), met.py() - matchedMet.py() );
280  profileRMS_delta_ex_VS_set_->Fill( matchedMet.sumEt(), met.px() - matchedMet.px() );
281  profileRMS_delta_ex_VS_set_->Fill( matchedMet.sumEt(), met.py() - matchedMet.py() );
282  }
283  if (delta_set_VS_set_) {
284  delta_set_VS_set_->Fill( matchedMet.sumEt(), met.sumEt() - matchedMet.sumEt() );
285  profile_delta_set_VS_set_->Fill( matchedMet.sumEt(), met.sumEt() - matchedMet.sumEt() );
286  profileRMS_delta_set_VS_set_->Fill( matchedMet.sumEt(), met.sumEt() - matchedMet.sumEt() );
287  }
288  if (delta_set_Over_set_VS_set_ && matchedMet.sumEt()>0.001 ) {
289  float setRes = (met.sumEt() - matchedMet.sumEt()) / matchedMet.sumEt();
290  if (setRes > maxVal) maxVal = setRes;
291  if (setRes < minVal) minVal = setRes;
292  delta_set_Over_set_VS_set_->Fill( matchedMet.sumEt(), setRes );
293  profile_delta_set_Over_set_VS_set_->Fill( matchedMet.sumEt(), setRes );
294  profileRMS_delta_set_Over_set_VS_set_->Fill( matchedMet.sumEt(), setRes );
295  }
296  }
297 }
PFMETMonitor::sumEt_
TH1F * sumEt_
Definition: PFMETMonitor.h:44
PFMETMonitor::profile_delta_ex_VS_set_
TProfile * profile_delta_ex_VS_set_
Definition: PFMETMonitor.h:50
edm::ParameterSet::getParameter
T getParameter(std::string const &) const
PFMETMonitor::setParameters
void setParameters(Benchmark::Mode mode, float ptmin, float ptmax, float etamin, float etamax, float phimin, float phimax, bool metSpHistos)
set the parameters locally
Definition: PFMETMonitor.cc:67
PFMETMonitor::setDirectory
void setDirectory(TDirectory *dir)
set directory (to use in ROOT)
Definition: PFMETMonitor.cc:248
PFMETMonitor::profileRMS_delta_set_VS_set_
TProfile * profileRMS_delta_set_VS_set_
Definition: PFMETMonitor.h:55
PFMETMonitor::profile_delta_set_VS_set_
TProfile * profile_delta_set_VS_set_
Definition: PFMETMonitor.h:51
CaloMET_cfi.met
tuple met
____________________________________________________________________________||
Definition: CaloMET_cfi.py:4
MatchCandidateBenchmark::fillOne
void fillOne(const reco::Candidate &candidate, const reco::Candidate &matchedCandidate)
fill histograms with a given particle
Definition: MatchCandidateBenchmark.cc:265
PFMETMonitor::delta_ex_VS_set_
TH2F * delta_ex_VS_set_
Definition: PFMETMonitor.h:46
PFMETMonitor::delta_set_VS_set_
TH2F * delta_set_VS_set_
Definition: PFMETMonitor.h:47
CandidateBenchmark::fillOne
void fillOne(const reco::Candidate &candidate)
fill histograms with a given particle
Definition: CandidateBenchmark.cc:87
MatchCandidateBenchmark::setup
void setup()
book histograms
Definition: MatchCandidateBenchmark.cc:35
CandidateBenchmark::setup
void setup()
book histograms
Definition: CandidateBenchmark.cc:29
Benchmark
abstract base class
Definition: Benchmark.h:21
PFMETMonitor::histogramBooked_
bool histogramBooked_
Definition: PFMETMonitor.h:63
PFMETMonitor::delta_ex_
TH1F * delta_ex_
Definition: PFMETMonitor.h:45
Benchmark::setDirectory
virtual void setDirectory(TDirectory *dir)
Definition: Benchmark.cc:19
alignBH_cfg.mode
list mode
Definition: alignBH_cfg.py:7
PFMETMonitor::profileRMS_delta_set_Over_set_VS_set_
TProfile * profileRMS_delta_set_Over_set_VS_set_
Definition: PFMETMonitor.h:56
reco::LeafCandidate::py
virtual double py() const GCC11_FINAL
y coordinate of momentum vector
Definition: LeafCandidate.h:174
PFMETMonitor::profileRMS_delta_ex_VS_set_
TProfile * profileRMS_delta_ex_VS_set_
Definition: PFMETMonitor.h:54
reco::LeafCandidate::phi
virtual float phi() const GCC11_FINAL
momentum azimuthal angle
Definition: LeafCandidate.h:180
reco::MET::sumEt
double sumEt() const
Definition: MET.h:53
Benchmark::book2D
TH2F * book2D(const char *histname, const char *title, int nbinsx, float xmin, float xmax, int nbinsy, float ymin, float ymax)
book a 2D histogram, either with DQM or plain root.
Definition: Benchmark.cc:43
reco::MET
Definition: MET.h:39
reco::LeafCandidate::px
virtual double px() const GCC11_FINAL
x coordinate of momentum vector
Definition: LeafCandidate.h:172
PFMETMonitor::PFMETMonitor
PFMETMonitor(Benchmark::Mode mode=Benchmark::DEFAULT)
Definition: PFMETMonitor.cc:17
Benchmark::setParameters
void setParameters(Mode mode)
Definition: Benchmark.h:50
reco::LeafCandidate::eta
virtual float eta() const GCC11_FINAL
momentum pseudorapidity
Definition: LeafCandidate.h:184
PFMETMonitor::candBench_
CandidateBenchmark candBench_
Definition: PFMETMonitor.h:59
PFMETMonitor::fillOne
void fillOne(const reco::MET &met, const reco::MET &matchedMet, float &minVal, float &maxVal)
Definition: PFMETMonitor.cc:255
PFMETMonitor::~PFMETMonitor
virtual ~PFMETMonitor()
Definition: PFMETMonitor.cc:44
PFMETMonitor::delta_set_Over_set_VS_set_
TH2F * delta_set_Over_set_VS_set_
Definition: PFMETMonitor.h:48
PFMETMonitor::createMETSpecificHistos_
bool createMETSpecificHistos_
Definition: PFMETMonitor.h:62
Benchmark::setRange
void setRange(float ptMin, float ptMax, float etaMin, float etaMax, float phiMin, float phiMax)
Definition: Benchmark.h:52
PFMETMonitor::profile_delta_set_Over_set_VS_set_
TProfile * profile_delta_set_Over_set_VS_set_
Definition: PFMETMonitor.h:52
ptmin
double ptmin
Definition: HydjetWrapper.h:85
Benchmark::mode_
Mode mode_
Definition: Benchmark.h:96
dir
dbl *** dir
Definition: mlp_gen.cc:35
PFMETMonitor::setup
void setup()
book histograms
Definition: PFMETMonitor.cc:189
reco::LeafCandidate::pt
virtual float pt() const GCC11_FINAL
transverse momentum
Definition: LeafCandidate.h:178
Benchmark::isInRange
bool isInRange(float pt, float eta, float phi) const
Definition: Benchmark.h:59
PFMETMonitor::matchCandBench_
MatchCandidateBenchmark matchCandBench_
Definition: PFMETMonitor.h:60
edm::parameterSet
ParameterSet const & parameterSet(Provenance const &provenance)
Definition: Provenance.cc:11
Benchmark::bookProfile
TProfile * bookProfile(const char *histname, const char *title, int nbinsx, float xmin, float xmax, float ymin, float ymax, const char *option)
book a TProfile histogram, either with DQM or plain root.
Definition: Benchmark.cc:95
Benchmark::book1D
TH1F * book1D(const char *histname, const char *title, int nbins, float xmin, float xmax)
book a 1D histogram, either with DQM or plain root.
Definition: Benchmark.cc:25
Generated for CMSSW Reference Manual by doxygen 1.8.5