CMS 3D CMS Logo
 All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Properties Friends Macros Pages
VertexHistogramMaker.cc
Go to the documentation of this file.
1 #include "Validation/RecoVertex/interface/VertexHistogramMaker.h"
2 #include "FWCore/MessageLogger/interface/MessageLogger.h"
3 #include "FWCore/ParameterSet/interface/ParameterSet.h"
4 #include "FWCore/Framework/interface/Event.h"
5 #include "FWCore/Framework/interface/Run.h"
6 #include "FWCore/Framework/interface/LuminosityBlock.h"
7 #include "FWCore/ServiceRegistry/interface/Service.h"
8 #include "CommonTools/UtilAlgos/interface/TFileService.h"
9 #include "DataFormats/VertexReco/interface/Vertex.h"
10 #include "DataFormats/Luminosity/interface/LumiDetails.h"
11 #include "TH2F.h"
12 #include "TH1F.h"
13 #include "TProfile.h"
14 
15 
16 VertexHistogramMaker::VertexHistogramMaker():
17  m_currdir(0), m_maxLS(100), m_weightThreshold(0.5), m_trueOnly(true),
18  m_runHisto(true), m_runHistoProfile(true), m_runHistoBXProfile(true), m_runHistoBXProfile2D(false), m_runHisto2D(false),
19  m_bsConstrained(false),
20  m_histoParameters() { }
21 
22 VertexHistogramMaker::VertexHistogramMaker(const edm::ParameterSet& iConfig):
23  m_currdir(0),
24  m_maxLS(iConfig.getParameter<unsigned int>("maxLSBeforeRebin")),
25  m_weightThreshold(iConfig.getUntrackedParameter<double>("weightThreshold",0.5)),
26  m_trueOnly(iConfig.getUntrackedParameter<bool>("trueOnly",true)),
27  m_runHisto(iConfig.getUntrackedParameter<bool>("runHisto",true)),
28  m_runHistoProfile(iConfig.getUntrackedParameter<bool>("runHistoProfile",true)),
29  m_runHistoBXProfile(iConfig.getUntrackedParameter<bool>("runHistoBXProfile",true)),
30  m_runHistoBXProfile2D(iConfig.getUntrackedParameter<bool>("runHistoBXProfile2D",false)),
31  m_runHisto2D(iConfig.getUntrackedParameter<bool>("runHisto2D",false)),
32  m_bsConstrained(iConfig.getParameter<bool>("bsConstrained")),
33  m_histoParameters(iConfig.getUntrackedParameter<edm::ParameterSet>("histoParameters",edm::ParameterSet())),
34  m_rhm(false),m_fhm(true)
35 { }
36 
37 
38 VertexHistogramMaker::~VertexHistogramMaker() {
39 
40  delete m_currdir;
41 
42 }
43 
44 
45 void VertexHistogramMaker::book(const std::string dirname) {
46 
47  edm::Service<TFileService> tfserv;
48  TFileDirectory* currdir = &(*tfserv);
49 
50  if(dirname!="") {
51  currdir = new TFileDirectory(tfserv->mkdir(dirname));
52  m_currdir = currdir;
53  }
54 
55  edm::LogInfo("HistogramBooking") << "Vertex histogram booking in directory " << dirname;
56 
57  m_hnvtx = currdir->make<TH1F>("nvtx","Number of Vertices",60,-0.5,59.5);
58  m_hnvtx->GetXaxis()->SetTitle("vertices"); m_hnvtx->GetYaxis()->SetTitle("Events");
59 
60  m_hntruevtx = currdir->make<TH1F>("ntruevtx","Number of True Vertices",60,-0.5,59.5);
61  m_hntruevtx->GetXaxis()->SetTitle("vertices"); m_hntruevtx->GetYaxis()->SetTitle("Events");
62 
63  m_hntruevtxvslumi = currdir->make<TProfile>("ntruevtxvslumi","Number of True Vertices vs BX lumi",250,0.,10.);
64  m_hntruevtxvslumi->GetXaxis()->SetTitle("BX lumi [10^{30}cm^{-2}s^{-1}]"); m_hntruevtxvslumi->GetYaxis()->SetTitle("Vertices");
65 
66  m_hntruevtxvslumi2D = currdir->make<TH2D>("ntruevtxvslumi2D","Number of True Vertices vs BX lumi",250,0.,10.,100,-0.5,99.5);
67  m_hntruevtxvslumi2D->GetXaxis()->SetTitle("BX lumi [10^{30}cm^{-2}s^{-1}]"); m_hntruevtxvslumi2D->GetYaxis()->SetTitle("Vertices");
68 
69  m_hntracks = currdir->make<TH1F>("ntracks","Number of Tracks",300,-0.5,299.5);
70  m_hntracks->GetXaxis()->SetTitle("tracks"); m_hntracks->GetYaxis()->SetTitle("Vertices");
71 
72  m_hsqsumptsq = currdir->make<TH1F>("sqsumptsq","sqrt(sum pt**2)",1000,0.,1000.);
73  m_hsqsumptsq->GetXaxis()->SetTitle("sqrt(#Sigma pt^{2}) (GeV)"); m_hsqsumptsq->GetYaxis()->SetTitle("Vertices");
74 
75  char htitle[300];
76  sprintf(htitle,"sqrt(sum pt**2) of Tracks weight > %f",m_weightThreshold);
77  m_hsqsumptsqheavy = currdir->make<TH1F>("sqsumptsqheavy",htitle,1000,0.,1000.);
78  m_hsqsumptsqheavy->GetXaxis()->SetTitle("sqrt(#Sigma pt^{2}) (GeV)"); m_hsqsumptsqheavy->GetYaxis()->SetTitle("Vertices");
79 
80  sprintf(htitle,"Number of Tracks weight > %f",m_weightThreshold);
81  m_hnheavytracks = currdir->make<TH1F>("nheavytracks",htitle,200,-0.5,199.5);
82  m_hnheavytracks->GetXaxis()->SetTitle("tracks"); m_hnheavytracks->GetYaxis()->SetTitle("Vertices");
83 
84  m_hndof = currdir->make<TH1F>("ndof","Number of degree of freedom",250,-0.5,499.5);
85  m_hndof->GetXaxis()->SetTitle("ndof"); m_hndof->GetYaxis()->SetTitle("Vertices");
86 
87  m_hndofvstracks = currdir->make<TH2F>("ndofvstracks","Ndof vs Ntracks",300,-0.5,299.5,250,-0.5,499.5);
88  m_hndofvstracks->GetXaxis()->SetTitle("tracks"); m_hndofvstracks->GetYaxis()->SetTitle("ndof");
89 
90  m_hndofvsvtxz = currdir->make<TProfile>("ndofvsvtxz","Ndof vs Vertex Z position",200,
91  m_histoParameters.getUntrackedParameter<double>("zMin",-20.),
92  m_histoParameters.getUntrackedParameter<double>("zMax",20.));
93  m_hndofvsvtxz->GetXaxis()->SetTitle("Z [cm]"); m_hndofvsvtxz->GetYaxis()->SetTitle("ndof");
94 
95  m_hntracksvsvtxz = currdir->make<TProfile>("ntracksvsvtxz","Ntracks vs Vertex Z position",200,
96  m_histoParameters.getUntrackedParameter<double>("zMin",-20.),
97  m_histoParameters.getUntrackedParameter<double>("zMax",20.));
98  m_hntracksvsvtxz->GetXaxis()->SetTitle("Z [cm]"); m_hntracksvsvtxz->GetYaxis()->SetTitle("tracks");
99 
100  m_haveweightvsvtxz = currdir->make<TProfile>("aveweightvsvtxz","Average weight vs Vertex Z position",200,
101  m_histoParameters.getUntrackedParameter<double>("zMin",-20.),
102  m_histoParameters.getUntrackedParameter<double>("zMax",20.));
103  m_haveweightvsvtxz->GetXaxis()->SetTitle("Z [cm]"); m_haveweightvsvtxz->GetYaxis()->SetTitle("Average weight");
104 
105  m_haveweightvsvtxzchk = currdir->make<TProfile>("aveweightvsvtxzchk","Average weight vs Vertex Z position (check)",200,
106  m_histoParameters.getUntrackedParameter<double>("zMin",-20.),
107  m_histoParameters.getUntrackedParameter<double>("zMax",20.));
108  m_haveweightvsvtxzchk->GetXaxis()->SetTitle("Z [cm]"); m_haveweightvsvtxzchk->GetYaxis()->SetTitle("Average weight");
109 
110  m_hweights = currdir->make<TH1F>("weights","Tracks weights",51,0.,1.02);
111  m_hweights->GetXaxis()->SetTitle("weights"); m_hweights->GetYaxis()->SetTitle("Tracks");
112 
113  m_haveweight = currdir->make<TH1F>("aveweight","Tracks average weights sum",51,0.,1.02);
114  m_haveweight->GetXaxis()->SetTitle("Average weight"); m_haveweight->GetYaxis()->SetTitle("Vertices");
115 
116 
117  m_hvtxx = currdir->make<TH1F>("vtxx","Vertex X position",
118  m_histoParameters.getUntrackedParameter<unsigned int>("nBinX",200),
119  m_histoParameters.getUntrackedParameter<double>("xMin",-1.),
120  m_histoParameters.getUntrackedParameter<double>("xMax",1.)
121  );
122  m_hvtxx->GetXaxis()->SetTitle("X [cm]"); m_hvtxx->GetYaxis()->SetTitle("Vertices");
123 
124  m_hvtxy = currdir->make<TH1F>("vtxy","Vertex Y position",
125  m_histoParameters.getUntrackedParameter<unsigned int>("nBinY",200),
126  m_histoParameters.getUntrackedParameter<double>("yMin",-1.),
127  m_histoParameters.getUntrackedParameter<double>("yMax",1.)
128  );
129  m_hvtxy->GetXaxis()->SetTitle("Y [cm]"); m_hvtxy->GetYaxis()->SetTitle("Vertices");
130 
131  m_hvtxz = currdir->make<TH1F>("vtxz","Vertex Z position",
132  m_histoParameters.getUntrackedParameter<unsigned int>("nBinZ",200),
133  m_histoParameters.getUntrackedParameter<double>("zMin",-20.),
134  m_histoParameters.getUntrackedParameter<double>("zMax",20.)
135  );
136  m_hvtxz->GetXaxis()->SetTitle("Z [cm]"); m_hvtxz->GetYaxis()->SetTitle("Vertices");
137 
138  if(m_runHisto) {
139  m_hvtxxrun = m_rhm.makeTH1F("vtxxrun","Vertex X position",
140  m_histoParameters.getUntrackedParameter<unsigned int>("nBinX",200),
141  m_histoParameters.getUntrackedParameter<double>("xMin",-1.),
142  m_histoParameters.getUntrackedParameter<double>("xMax",1.));
143 
144  m_hvtxyrun = m_rhm.makeTH1F("vtxyrun","Vertex Y position",
145  m_histoParameters.getUntrackedParameter<unsigned int>("nBinY",200),
146  m_histoParameters.getUntrackedParameter<double>("yMin",-1.),
147  m_histoParameters.getUntrackedParameter<double>("yMax",1.));
148 
149  m_hvtxzrun = m_rhm.makeTH1F("vtxzrun","Vertex Z position",
150  m_histoParameters.getUntrackedParameter<unsigned int>("nBinZ",200),
151  m_histoParameters.getUntrackedParameter<double>("zMin",-20.),
152  m_histoParameters.getUntrackedParameter<double>("zMax",20.));
153 
154  if(m_runHistoProfile) {
155  m_hvtxxvsorbrun = m_rhm.makeTProfile("vtxxvsorbrun","Vertex X position vs orbit number",4*m_maxLS,0.5,m_maxLS*262144+0.5);
156  m_hvtxyvsorbrun = m_rhm.makeTProfile("vtxyvsorbrun","Vertex Y position vs orbit number",4*m_maxLS,0.5,m_maxLS*262144+0.5);
157  m_hvtxzvsorbrun = m_rhm.makeTProfile("vtxzvsorbrun","Vertex Z position vs orbit number",4*m_maxLS,0.5,m_maxLS*262144+0.5);
158  m_hnvtxvsorbrun = m_rhm.makeTProfile("nvtxvsorbrun","Number of true vertices vs orbit number",m_maxLS,0.5,m_maxLS*262144+0.5);
159  }
160 
161  if(m_runHisto2D) {
162  m_hnvtxvsbxvsorbrun = m_rhm.makeTProfile2D("nvtxvsbxvsorbrun","Number of true vertices vs BX vs orbit number",
163  3564,-0.5,3563.5,m_maxLS,0.5,m_maxLS*262144+0.5);
164  m_hnvtxvsorbrun2D = m_rhm.makeTH2F("nvtxvsorbrun2D","Number of true vertices vs orbit number",
165  m_maxLS,0.5,m_maxLS*262144+0.5,60,-0.5,59.5);
166  }
167 
168  if(m_runHistoBXProfile) {
169  m_hvtxxvsbxrun = m_fhm.makeTProfile("vtxxvsbxrun","Vertex X position vs BX number",3564,-0.5,3563.5);
170  m_hvtxyvsbxrun = m_fhm.makeTProfile("vtxyvsbxrun","Vertex Y position vs BX number",3564,-0.5,3563.5);
171  m_hvtxzvsbxrun = m_fhm.makeTProfile("vtxzvsbxrun","Vertex Z position vs BX number",3564,-0.5,3563.5);
172 
173  m_hnvtxvsbxrun = m_rhm.makeTProfile("nvtxvsbxrun","Number of true vertices vs BX number",3564,-0.5,3563.5);
174 
175  if(m_runHistoBXProfile2D) {
176  m_hnvtxvsbxvslumirun = m_fhm.makeTProfile2D("nvtxvsbxvslumirun","Number of vertices vs BX and BX lumi",3564,-0.5,3563.5,250,0.,10.);
177  }
178  if(m_runHisto2D) {
179  m_hvtxxvsbx2drun = m_fhm.makeTH2F("vtxxvsbx2drun","Vertex X position vs BX number",3564,-0.5,3563.5,
180  m_histoParameters.getUntrackedParameter<unsigned int>("nBinX",200),
181  m_histoParameters.getUntrackedParameter<double>("xMin",-1.),
182  m_histoParameters.getUntrackedParameter<double>("xMax",1.));
183  m_hvtxyvsbx2drun = m_fhm.makeTH2F("vtxyvsbx2drun","Vertex Y position vs BX number",3564,-0.5,3563.5,
184  m_histoParameters.getUntrackedParameter<unsigned int>("nBinY",200),
185  m_histoParameters.getUntrackedParameter<double>("yMin",-1.),
186  m_histoParameters.getUntrackedParameter<double>("yMax",1.));
187  m_hvtxzvsbx2drun = m_fhm.makeTH2F("vtxzvsbx2drun","Vertex Z position vs BX number",3564,-0.5,3563.5,
188  m_histoParameters.getUntrackedParameter<unsigned int>("nBinZ",200),
189  m_histoParameters.getUntrackedParameter<double>("zMin",-20.),
190  m_histoParameters.getUntrackedParameter<double>("zMax",20.));
191  }
192  }
193 
194 
195  }
196 }
197 
198 void VertexHistogramMaker::beginRun(const edm::Run& iRun) {
199 
200  TFileDirectory* currdir = m_currdir;
201  if(currdir==0) {
202  edm::Service<TFileService> tfserv;
203  currdir = &(*tfserv);
204  }
205 
206  m_rhm.beginRun(iRun,*currdir);
207  m_fhm.beginRun(iRun,*currdir);
208 
209 
210  if(m_runHisto) {
211  (*m_hvtxxrun)->GetXaxis()->SetTitle("X [cm]"); (*m_hvtxxrun)->GetYaxis()->SetTitle("Vertices");
212  (*m_hvtxyrun)->GetXaxis()->SetTitle("Y [cm]"); (*m_hvtxyrun)->GetYaxis()->SetTitle("Vertices");
213  (*m_hvtxzrun)->GetXaxis()->SetTitle("Z [cm]"); (*m_hvtxzrun)->GetYaxis()->SetTitle("Vertices");
214 
215  if(m_runHistoProfile) {
216  (*m_hvtxxvsorbrun)->GetXaxis()->SetTitle("time [orbit#]"); (*m_hvtxxvsorbrun)->GetYaxis()->SetTitle("X [cm]");
217  (*m_hvtxxvsorbrun)->SetBit(TH1::kCanRebin);
218  (*m_hvtxyvsorbrun)->GetXaxis()->SetTitle("time [orbit#]"); (*m_hvtxyvsorbrun)->GetYaxis()->SetTitle("Y [cm]");
219  (*m_hvtxyvsorbrun)->SetBit(TH1::kCanRebin);
220  (*m_hvtxzvsorbrun)->GetXaxis()->SetTitle("time [orbit#]"); (*m_hvtxzvsorbrun)->GetYaxis()->SetTitle("Z [cm]");
221  (*m_hvtxzvsorbrun)->SetBit(TH1::kCanRebin);
222  (*m_hnvtxvsorbrun)->GetXaxis()->SetTitle("time [orbit#]"); (*m_hnvtxvsorbrun)->GetYaxis()->SetTitle("Nvertices");
223  (*m_hnvtxvsorbrun)->SetBit(TH1::kCanRebin);
224  }
225 
226  if(m_runHistoBXProfile) {
227  (*m_hvtxxvsbxrun)->GetXaxis()->SetTitle("BX"); (*m_hvtxxvsbxrun)->GetYaxis()->SetTitle("X [cm]");
228  (*m_hvtxyvsbxrun)->GetXaxis()->SetTitle("BX"); (*m_hvtxyvsbxrun)->GetYaxis()->SetTitle("Y [cm]");
229  (*m_hvtxzvsbxrun)->GetXaxis()->SetTitle("BX"); (*m_hvtxzvsbxrun)->GetYaxis()->SetTitle("Z [cm]");
230  (*m_hnvtxvsbxrun)->GetXaxis()->SetTitle("BX"); (*m_hnvtxvsbxrun)->GetYaxis()->SetTitle("Nvertices");
231  if(m_runHistoBXProfile2D) {
232  (*m_hnvtxvsbxvslumirun)->GetXaxis()->SetTitle("BX"); (*m_hnvtxvsbxvslumirun)->GetYaxis()->SetTitle("BX lumi [10^{30}cm^{-2}s^{-1}]");
233  }
234  if(m_runHisto2D) {
235  (*m_hvtxxvsbx2drun)->GetXaxis()->SetTitle("BX"); (*m_hvtxxvsbx2drun)->GetYaxis()->SetTitle("X [cm]");
236  (*m_hvtxyvsbx2drun)->GetXaxis()->SetTitle("BX"); (*m_hvtxyvsbx2drun)->GetYaxis()->SetTitle("Y [cm]");
237  (*m_hvtxzvsbx2drun)->GetXaxis()->SetTitle("BX"); (*m_hvtxzvsbx2drun)->GetYaxis()->SetTitle("Z [cm]");
238  }
239  }
240 
241  if(m_runHisto2D) {
242  (*m_hnvtxvsbxvsorbrun)->GetXaxis()->SetTitle("BX#"); (*m_hnvtxvsbxvsorbrun)->GetYaxis()->SetTitle("time [orbit#]");
243  (*m_hnvtxvsbxvsorbrun)->SetBit(TH1::kCanRebin);
244  (*m_hnvtxvsorbrun2D)->GetXaxis()->SetTitle("time [orbit#]"); (*m_hnvtxvsorbrun2D)->GetYaxis()->SetTitle("Nvertices");
245  (*m_hnvtxvsorbrun2D)->SetBit(TH1::kCanRebin);
246  }
247  }
248 }
249 
250 void VertexHistogramMaker::fill(const unsigned int orbit, const int bx, const float bxlumi, const reco::VertexCollection& vertices, const double weight) {
251 
252  m_hnvtx->Fill(vertices.size(),weight);
253 
254  int ntruevtx = 0;
255  for(reco::VertexCollection::const_iterator vtx=vertices.begin();vtx!=vertices.end();++vtx) {
256  if(!vtx->isFake()) ntruevtx++;
257 
258  if(!(m_trueOnly && vtx->isFake())) {
259 
260  double aveweight = m_bsConstrained ? vtx->ndof()/(2.*vtx->tracksSize()) : (vtx->ndof()+3)/(2.*vtx->tracksSize());
261 
262  m_hntracks->Fill(vtx->tracksSize(),weight);
263  m_hndof->Fill(vtx->ndof(),weight);
264  m_haveweight->Fill(aveweight,weight);
265  m_hndofvstracks->Fill(vtx->tracksSize(),vtx->ndof(),weight);
266  m_hndofvsvtxz->Fill(vtx->z(),vtx->ndof(),weight);
267  m_hntracksvsvtxz->Fill(vtx->z(),vtx->tracksSize(),weight);
268  m_haveweightvsvtxz->Fill(vtx->z(),aveweight,weight);
269 
270  m_hvtxx->Fill(vtx->x(),weight);
271  m_hvtxy->Fill(vtx->y(),weight);
272  m_hvtxz->Fill(vtx->z(),weight);
273 
274  if(m_runHisto) {
275  if(m_hvtxxrun && *m_hvtxxrun ) (*m_hvtxxrun)->Fill(vtx->x(),weight);
276  if(m_hvtxyrun && *m_hvtxyrun ) (*m_hvtxyrun)->Fill(vtx->y(),weight);
277  if(m_hvtxzrun && *m_hvtxzrun ) (*m_hvtxzrun)->Fill(vtx->z(),weight);
278  if(m_runHistoProfile) {
279  if(m_hvtxxvsorbrun && *m_hvtxxvsorbrun ) (*m_hvtxxvsorbrun)->Fill(orbit,vtx->x(),weight);
280  if(m_hvtxyvsorbrun && *m_hvtxyvsorbrun ) (*m_hvtxyvsorbrun)->Fill(orbit,vtx->y(),weight);
281  if(m_hvtxzvsorbrun && *m_hvtxzvsorbrun ) (*m_hvtxzvsorbrun)->Fill(orbit,vtx->z(),weight);
282  }
283  if(m_runHistoBXProfile) {
284  if(m_hvtxxvsbxrun && *m_hvtxxvsbxrun ) (*m_hvtxxvsbxrun)->Fill(bx,vtx->x(),weight);
285  if(m_hvtxyvsbxrun && *m_hvtxyvsbxrun ) (*m_hvtxyvsbxrun)->Fill(bx,vtx->y(),weight);
286  if(m_hvtxzvsbxrun && *m_hvtxzvsbxrun ) (*m_hvtxzvsbxrun)->Fill(bx,vtx->z(),weight);
287  if(m_runHisto2D) {
288  if(m_hvtxxvsbx2drun && *m_hvtxxvsbx2drun ) (*m_hvtxxvsbx2drun)->Fill(bx,vtx->x(),weight);
289  if(m_hvtxyvsbx2drun && *m_hvtxyvsbx2drun ) (*m_hvtxyvsbx2drun)->Fill(bx,vtx->y(),weight);
290  if(m_hvtxzvsbx2drun && *m_hvtxzvsbx2drun ) (*m_hvtxzvsbx2drun)->Fill(bx,vtx->z(),weight);
291  }
292  }
293  }
294 
295  int nheavytracks = 0;
296  double sumpt2 = 0.;
297  double sumpt2heavy = 0.;
298 
299  for(reco::Vertex::trackRef_iterator trk=vtx->tracks_begin();trk!=vtx->tracks_end();++trk) {
300 
301  sumpt2 += (*trk)->pt()*(*trk)->pt();
302 
303  if(vtx->trackWeight(*trk) > m_weightThreshold) {
304  nheavytracks++;
305  sumpt2heavy += (*trk)->pt()*(*trk)->pt();
306  }
307 
308  m_hweights->Fill(vtx->trackWeight(*trk),weight);
309  m_haveweightvsvtxzchk->Fill(vtx->z(),vtx->trackWeight(*trk),weight);
310 
311  }
312 
313  m_hnheavytracks->Fill(nheavytracks,weight);
314  m_hsqsumptsq->Fill(sqrt(sumpt2),weight);
315  m_hsqsumptsqheavy->Fill(sqrt(sumpt2heavy),weight);
316 
317 
318  }
319 
320 
321  }
322 
323  m_hntruevtx->Fill(ntruevtx,weight);
324 
325  if(bxlumi >= 0.) {
326  m_hntruevtxvslumi->Fill(bxlumi,ntruevtx,weight);
327  m_hntruevtxvslumi2D->Fill(bxlumi,ntruevtx,weight);
328  }
329 
330  if(m_runHisto) {
331  if(m_runHistoProfile) {
332  if(m_hnvtxvsorbrun && *m_hnvtxvsorbrun ) (*m_hnvtxvsorbrun)->Fill(orbit,ntruevtx,weight);
333  }
334  if(m_runHistoBXProfile) {
335  if(m_hnvtxvsbxrun && *m_hnvtxvsbxrun ) (*m_hnvtxvsbxrun)->Fill(bx,ntruevtx,weight);
336  if(m_runHistoBXProfile2D) {
337  if(m_hnvtxvsbxvslumirun && *m_hnvtxvsbxvslumirun && bxlumi >= 0.) (*m_hnvtxvsbxvslumirun)->Fill(bx,bxlumi,ntruevtx,weight);
338  }
339  }
340  if(m_runHisto2D) {
341  if(m_hnvtxvsbxvsorbrun && *m_hnvtxvsbxvsorbrun ) (*m_hnvtxvsbxvsorbrun)->Fill(bx,orbit,ntruevtx,weight);
342  if(m_hnvtxvsorbrun2D && *m_hnvtxvsorbrun2D ) {
343  if(ntruevtx < (*m_hnvtxvsorbrun2D)->GetYaxis()->GetXmax() && ntruevtx > (*m_hnvtxvsorbrun2D)->GetYaxis()->GetXmin()) {
344  (*m_hnvtxvsorbrun2D)->Fill(orbit,ntruevtx,weight);
345  }
346  }
347  }
348  }
349 
350 
351 }
352 
353 void VertexHistogramMaker::fill(const edm::Event& iEvent, const reco::VertexCollection& vertices, const double weight) {
354 
355  TFileDirectory* currdir = m_currdir;
356  if(currdir==0) {
357  edm::Service<TFileService> tfserv;
358  currdir = &(*tfserv);
359  }
360 
361 
362  // get luminosity
363 
364  edm::Handle<LumiDetails> ld;
365  iEvent.getLuminosityBlock().getByLabel("lumiProducer",ld);
366 
367  float bxlumi = -1.;
368 
369  if(ld.isValid()) {
370  if(ld->isValid()) {
371  bxlumi = ld->lumiValue(LumiDetails::kOCC1,iEvent.bunchCrossing())*6.37;
372  }
373  }
374 
375  fill(iEvent.orbitNumber(),iEvent.bunchCrossing(),bxlumi,vertices,weight);
376 
377 }
edm::ParameterSet::getUntrackedParameter
T getUntrackedParameter(std::string const &, T const &) const
VertexHistogramMaker::m_currdir
TFileDirectory * m_currdir
Definition: VertexHistogramMaker.h:36
VertexHistogramMaker::m_histoParameters
const edm::ParameterSet m_histoParameters
Definition: VertexHistogramMaker.h:46
VertexHistogramMaker::m_hnvtxvsbxvsorbrun
TProfile2D ** m_hnvtxvsbxvsorbrun
Definition: VertexHistogramMaker.h:76
VertexHistogramMaker::m_rhm
RunHistogramManager m_rhm
Definition: VertexHistogramMaker.h:48
edm::EventBase::bunchCrossing
int bunchCrossing() const
Definition: EventBase.h:62
VertexHistogramMaker::m_hnvtxvsbxvslumirun
TProfile2D ** m_hnvtxvsbxvslumirun
Definition: VertexHistogramMaker.h:84
reco::VertexCollection
std::vector< Vertex > VertexCollection
collection of Vertex objects
Definition: VertexFwd.h:9
VertexHistogramMaker::fill
void fill(const edm::Event &iEvent, const reco::VertexCollection &vertices, const double weight=1.)
Definition: VertexHistogramMaker.cc:353
edm::LuminosityBlock::getByLabel
bool getByLabel(std::string const &label, Handle< PROD > &result) const
Definition: LuminosityBlock.h:168
AlCaHLTBitMon_QueryRunRegistry.string
string string
Definition: AlCaHLTBitMon_QueryRunRegistry.py:255
RunHistogramManager::makeTH1F
TH1F ** makeTH1F(const char *name, const char *title, const unsigned int nbinx, const double xmin, const double xmax)
Definition: RunHistogramManager.cc:29
iEvent
int iEvent
Definition: GenABIO.cc:243
mathSSE::sqrt
T sqrt(T t)
Definition: SSEVec.h:48
RunHistogramManager::makeTProfile
TProfile ** makeTProfile(const char *name, const char *title, const unsigned int nbinx, const double xmin, const double xmax)
Definition: RunHistogramManager.cc:53
edm::Event::getLuminosityBlock
LuminosityBlock const & getLuminosityBlock() const
Definition: Event.h:68
edm::EventBase::orbitNumber
int orbitNumber() const
Definition: EventBase.h:63
edm::HandleBase::isValid
bool isValid() const
Definition: HandleBase.h:76
VertexHistogramMaker::book
void book(const std::string dirname="")
Definition: VertexHistogramMaker.cc:45
VertexHistogramMaker::beginRun
void beginRun(const edm::Run &iRun)
Definition: VertexHistogramMaker.cc:198
TFileDirectory::mkdir
TFileDirectory mkdir(const std::string &dir, const std::string &descr="")
create a new subdirectory
Definition: TFileDirectory.cc:162
RunHistogramManager::makeTProfile2D
TProfile2D ** makeTProfile2D(const char *name, const char *title, const unsigned int nbinx, const double xmin, const double xmax, const unsigned int nbiny, const double ymin, const double ymax)
Definition: RunHistogramManager.cc:78
RunHistogramManager::beginRun
void beginRun(const edm::Run &iRun)
Definition: RunHistogramManager.cc:90
reco::Vertex::trackRef_iterator
std::vector< TrackBaseRef >::const_iterator trackRef_iterator
The iteratator for the vector&lt;TrackRef&gt;
Definition: Vertex.h:38
TFileDirectory::make
T * make() const
make new ROOT object
Definition: TFileDirectory.h:55
RunHistogramManager::makeTH2F
TH2F ** makeTH2F(const char *name, const char *title, const unsigned int nbinx, const double xmin, const double xmax, const unsigned int nbiny, const double ymin, const double ymax)
Definition: RunHistogramManager.cc:66
histoStyle.weight
int weight
Definition: histoStyle.py:50
VertexHistogramMaker::m_maxLS
const unsigned int m_maxLS
Definition: VertexHistogramMaker.h:37
edm::Run
Definition: Run.h:36
VertexHistogramMaker::m_fhm
RunHistogramManager m_fhm
Definition: VertexHistogramMaker.h:49
Generated for CMSSW Reference Manual by doxygen 1.8.5