test
CMS 3D CMS Logo
 All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Properties Friends Macros Pages
PrimaryVertexMonitor.cc
Go to the documentation of this file.
1 #include "DQMOffline/RecoB/plugins/PrimaryVertexMonitor.h"
2 #include "FWCore/ServiceRegistry/interface/Service.h"
3 #include "FWCore/Framework/interface/MakerMacros.h"
4 #include "FWCore/Utilities/interface/isFinite.h"
5 
6 
7 #include "CommonTools/Statistics/interface/ChiSquaredProbability.h"
8 
9 #include "FWCore/MessageLogger/interface/MessageLogger.h"
10 
11 #include "TMath.h"
12 
13 using namespace reco;
14 using namespace edm;
15 
16 PrimaryVertexMonitor::PrimaryVertexMonitor(const edm::ParameterSet& pSet)
17  : conf_ ( pSet )
18  , TopFolderName_ ( pSet.getParameter<std::string>("TopFolderName") )
19  , AlignmentLabel_( pSet.getParameter<std::string>("AlignmentLabel"))
20  , nbvtx(NULL)
21  , bsX(NULL)
22  , bsY(NULL)
23  , bsZ(NULL)
24  , bsSigmaZ(NULL)
25  , bsDxdz(NULL)
26  , bsDydz(NULL)
27  , bsBeamWidthX(NULL)
28  , bsBeamWidthY(NULL)
29  , bsType(NULL)
30  , sumpt(NULL)
31  , ntracks(NULL)
32  , weight(NULL)
33  , chi2ndf(NULL)
34  , chi2prob(NULL)
35  , dxy(NULL)
36  , dxy2(NULL)
37  , dz(NULL)
38  , dxyErr(NULL)
39  , dzErr(NULL)
40  , dxyVsPhi_pt1(NULL)
41  , dzVsPhi_pt1(NULL)
42  , dxyVsEta_pt1(NULL)
43  , dzVsEta_pt1(NULL)
44  , dxyVsPhi_pt10(NULL)
45  , dzVsPhi_pt10(NULL)
46  , dxyVsEta_pt10(NULL)
47  , dzVsEta_pt10(NULL)
48 {
49  // dqmStore_ = edm::Service<DQMStore>().operator->();
50 
51 
52  vertexInputTag_ = pSet.getParameter<InputTag>("vertexLabel");
53  beamSpotInputTag_ = pSet.getParameter<InputTag>("beamSpotLabel");
54  vertexToken_ = consumes<reco::VertexCollection>(vertexInputTag_);
55  scoreToken_ = consumes<VertexScore>(vertexInputTag_);
56  beamspotToken_ = consumes<reco::BeamSpot> (beamSpotInputTag_);
57 
58 }
59 
60 // -- BeginRun
61 //---------------------------------------------------------------------------------//
62 void
63 PrimaryVertexMonitor::bookHistograms(DQMStore::IBooker &iBooker,
64  edm::Run const &, edm::EventSetup const &) {
65 
66  std::string dqmLabel = "";
67 
68  //
69  // Book all histograms.
70  //
71 
72  // get the store
73  dqmLabel = TopFolderName_+"/"+vertexInputTag_.label();
74  iBooker.setCurrentFolder(dqmLabel);
75 
76 // xPos = iBooker.book1D ("xPos","x Coordinate" ,100, -0.1, 0.1);
77 
78  nbvtx = iBooker.book1D("vtxNbr","Reconstructed Vertices in Event",50,-0.5,49.5);
79  nbgvtx = iBooker.book1D("goodvtxNbr","Reconstructed Good Vertices in Event",50,-0.5,49.5);
80 
81  // to be configured each year...
82  auto vposx = conf_.getParameter<double>("Xpos");
83  auto vposy = conf_.getParameter<double>("Ypos");
84 
85  nbtksinvtx[0] = iBooker.book1D("otherVtxTrksNbr","Reconstructed Tracks in Vertex (other Vtx)",40,-0.5,99.5);
86  ntracksVsZ[0] = iBooker.bookProfile("otherVtxTrksVsZ","Reconstructed Tracks in Vertex (other Vtx) vs Z",80,-20.,20.,50,0,100,"");
87  ntracksVsZ[0]->setAxisTitle("z-bs",1);
88  ntracksVsZ[0]->setAxisTitle("#tracks",2);
89 
90  score[0] = iBooker.book1D("otherVtxScore","sqrt(score) (other Vtx)",100,0.,400.);
91  trksWeight[0] = iBooker.book1D("otherVtxTrksWeight","Total weight of Tracks in Vertex (other Vtx)",40,0,100.);
92  vtxchi2[0] = iBooker.book1D("otherVtxChi2","#chi^{2} (other Vtx)",100,0.,200.);
93  vtxndf[0] = iBooker.book1D("otherVtxNdf","ndof (other Vtx)",100,0.,200.);
94  vtxprob[0] = iBooker.book1D("otherVtxProb","#chi^{2} probability (other Vtx)",100,0.,1.);
95  nans[0] = iBooker.book1D("otherVtxNans","Illegal values for x,y,z,xx,xy,xz,yy,yz,zz (other Vtx)",9,0.5,9.5);
96 
97  nbtksinvtx[1] = iBooker.book1D("tagVtxTrksNbr","Reconstructed Tracks in Vertex (tagged Vtx)",100,-0.5,99.5);
98  ntracksVsZ[1] = iBooker.bookProfile("tagVtxTrksVsZ","Reconstructed Tracks in Vertex (tagged Vtx) vs Z",80,-20.,20.,50,0,100,"");
99  ntracksVsZ[1]->setAxisTitle("z-bs",1);
100  ntracksVsZ[1]->setAxisTitle("#tracks",2);
101 
102  score[1] = iBooker.book1D("tagVtxScore","sqrt(score) (tagged Vtx)",100,0.,400.);
103  trksWeight[1] = iBooker.book1D("tagVtxTrksWeight","Total weight of Tracks in Vertex (tagged Vtx)",100,0,100.);
104  vtxchi2[1] = iBooker.book1D("tagVtxChi2","#chi^{2} (tagged Vtx)",100,0.,200.);
105  vtxndf[1] = iBooker.book1D("tagVtxNdf","ndof (tagged Vtx)",100,0.,200.);
106  vtxprob[1] = iBooker.book1D("tagVtxProb","#chi^{2} probability (tagged Vtx)",100,0.,1.);
107  nans[1] = iBooker.book1D("tagVtxNans","Illegal values for x,y,z,xx,xy,xz,yy,yz,zz (tagged Vtx)",9,0.5,9.5);
108 
109  xrec[0] = iBooker.book1D("otherPosX","Position x Coordinate (other Vtx)",100,vposx-0.1,vposx+0.1);
110  yrec[0] = iBooker.book1D("otherPosY","Position y Coordinate (other Vtx)",100,vposy-0.1,vposy+0.1);
111  zrec[0] = iBooker.book1D("otherPosZ","Position z Coordinate (other Vtx)",100,-20.,20.);
112  xDiff[0] = iBooker.book1D("otherDiffX","X distance from BeamSpot (other Vtx)",100,-500,500);
113  yDiff[0] = iBooker.book1D("otherDiffY","Y distance from BeamSpot (other Vtx)",100,-500,500);
114  xerr[0] = iBooker.book1D("otherErrX","Uncertainty x Coordinate (other Vtx)",100,0.,100);
115  yerr[0] = iBooker.book1D("otherErrY","Uncertainty y Coordinate (other Vtx)",100,0.,100);
116  zerr[0] = iBooker.book1D("otherErrZ","Uncertainty z Coordinate (other Vtx)",100,0.,100);
117  xerrVsTrks[0] = iBooker.book2D("otherErrVsWeightX","Uncertainty x Coordinate vs. track weight (other Vtx)",100,0,100.,100,0.,100);
118  yerrVsTrks[0] = iBooker.book2D("otherErrVsWeightY","Uncertainty y Coordinate vs. track weight (other Vtx)",100,0,100.,100,0.,100);
119  zerrVsTrks[0] = iBooker.book2D("otherErrVsWeightZ","Uncertainty z Coordinate vs. track weight (other Vtx)",100,0,100.,100,0.,100);
120 
121 
122  xrec[1] = iBooker.book1D("tagPosX","Position x Coordinate (tagged Vtx)",100,vposx-0.1,vposx+0.1);
123  yrec[1] = iBooker.book1D("tagPosY","Position y Coordinate (tagged Vtx)",100,vposx-0.1,vposy+0.1);
124  zrec[1] = iBooker.book1D("tagPosZ","Position z Coordinate (tagged Vtx)",100,-20.,20.);
125  xDiff[1] = iBooker.book1D("tagDiffX","X distance from BeamSpot (tagged Vtx)",100,-500, 500);
126  yDiff[1] = iBooker.book1D("tagDiffY","Y distance from BeamSpot (tagged Vtx)",100,-500, 500);
127  xerr[1] = iBooker.book1D("tagErrX","Uncertainty x Coordinate (tagged Vtx)",100,0.,100);
128  yerr[1] = iBooker.book1D("tagErrY","Uncertainty y Coordinate (tagged Vtx)",100,0.,100);
129  zerr[1] = iBooker.book1D("tagErrZ","Uncertainty z Coordinate (tagged Vtx)",100,0.,100);
130  xerrVsTrks[1] = iBooker.book2D("tagErrVsWeightX","Uncertainty x Coordinate vs. track weight (tagged Vtx)",100,0,100.,100,0.,100);
131  yerrVsTrks[1] = iBooker.book2D("tagErrVsWeightY","Uncertainty y Coordinate vs. track weight (tagged Vtx)",100,0,100.,100,0.,100);
132  zerrVsTrks[1] = iBooker.book2D("tagErrVsWeightZ","Uncertainty z Coordinate vs. track weight (tagged Vtx)",100,0,100.,100,0.,100);
133 
134  type[0] = iBooker.book1D("otherType","Vertex type (other Vtx)",3,-0.5,2.5);
135  type[1] = iBooker.book1D("tagType","Vertex type (tagged Vtx)",3,-0.5,2.5);
136  for (int i=0;i<2;++i){
137  type[i]->getTH1F()->GetXaxis()->SetBinLabel(1,"Valid, real");
138  type[i]->getTH1F()->GetXaxis()->SetBinLabel(2,"Valid, fake");
139  type[i]->getTH1F()->GetXaxis()->SetBinLabel(3,"Invalid");
140  }
141 
142 
143  // get the store
144  dqmLabel = TopFolderName_+"/"+beamSpotInputTag_.label();
145  iBooker.setCurrentFolder(dqmLabel);
146 
147  bsX = iBooker.book1D("bsX", "BeamSpot x0", 100,-0.1,0.1);
148  bsY = iBooker.book1D("bsY", "BeamSpot y0", 100,-0.1,0.1);
149  bsZ = iBooker.book1D("bsZ", "BeamSpot z0", 100,-2.,2.);
150  bsSigmaZ = iBooker.book1D("bsSigmaZ", "BeamSpot sigmaZ", 100, 0., 10. );
151  bsDxdz = iBooker.book1D("bsDxdz", "BeamSpot dxdz", 100, -0.0003, 0.0003);
152  bsDydz = iBooker.book1D("bsDydz", "BeamSpot dydz", 100, -0.0003, 0.0003);
153  bsBeamWidthX = iBooker.book1D("bsBeamWidthX", "BeamSpot BeamWidthX", 100, 0., 100.);
154  bsBeamWidthY = iBooker.book1D("bsBeamWidthY", "BeamSpot BeamWidthY", 100, 0., 100.);
155  bsType = iBooker.book1D("bsType", "BeamSpot type", 4, -1.5, 2.5);
156  bsType->getTH1F()->GetXaxis()->SetBinLabel(1,"Unknown");
157  bsType->getTH1F()->GetXaxis()->SetBinLabel(2,"Fake");
158  bsType->getTH1F()->GetXaxis()->SetBinLabel(3,"LHC");
159  bsType->getTH1F()->GetXaxis()->SetBinLabel(4,"Tracker");
160 
161 
162  // get the store
163  dqmLabel = TopFolderName_+"/"+AlignmentLabel_;
164  iBooker.setCurrentFolder(dqmLabel);
165 
166  int TKNoBin = conf_.getParameter<int>( "TkSizeBin");
167  double TKNoMin = conf_.getParameter<double>("TkSizeMin");
168  double TKNoMax = conf_.getParameter<double>("TkSizeMax");
169 
170  int DxyBin = conf_.getParameter<int>( "DxyBin");
171  double DxyMin = conf_.getParameter<double>("DxyMin");
172  double DxyMax = conf_.getParameter<double>("DxyMax");
173 
174  int DzBin = conf_.getParameter<int>( "DzBin");
175  double DzMin = conf_.getParameter<double>("DzMin");
176  double DzMax = conf_.getParameter<double>("DzMax");
177 
178  int PhiBin = conf_.getParameter<int>( "PhiBin");
179  double PhiMin = conf_.getParameter<double>("PhiMin");
180  double PhiMax = conf_.getParameter<double>("PhiMax");
181 
182  int EtaBin = conf_.getParameter<int>( "EtaBin");
183  double EtaMin = conf_.getParameter<double>("EtaMin");
184  double EtaMax = conf_.getParameter<double>("EtaMax");
185 
186 
187  ntracks = iBooker.book1D("ntracks","number of PV tracks (p_{T} > 1 GeV)", TKNoBin, TKNoMin, TKNoMax);
188  ntracks->setAxisTitle("Number of PV Tracks (p_{T} > 1 GeV) per Event", 1);
189  ntracks->setAxisTitle("Number of Event", 2);
190 
191  weight = iBooker.book1D("weight","weight of PV tracks (p_{T} > 1 GeV)", 100, 0., 1.);
192  weight->setAxisTitle("weight of PV Tracks (p_{T} > 1 GeV) per Event", 1);
193  weight->setAxisTitle("Number of Event", 2);
194 
195  sumpt = iBooker.book1D("sumpt", "#Sum p_{T} of PV tracks (p_{T} > 1 GeV)", 100,-0.5,249.5);
196  chi2ndf = iBooker.book1D("chi2ndf", "PV tracks (p_{T} > 1 GeV) #chi^{2}/ndof", 100, 0., 20. );
197  chi2prob = iBooker.book1D("chi2prob","PV tracks (p_{T} > 1 GeV) #chi^{2} probability",100, 0., 1. );
198  dxy = iBooker.book1D("dxy", "PV tracks (p_{T} > 1 GeV) d_{xy} (cm)", DxyBin, DxyMin, DxyMax);
199  dxy2 = iBooker.book1D("dxyzoom", "PV tracks (p_{T} > 1 GeV) d_{xy} (cm)", DxyBin, DxyMin/5., DxyMax/5.);
200  dz = iBooker.book1D("dz", "PV tracks (p_{T} > 1 GeV) d_{z} (cm)", DzBin, DzMin, DzMax );
201  dxyErr = iBooker.book1D("dxyErr", "PV tracks (p_{T} > 1 GeV) d_{xy} error (cm)", 100, 0., 0.2 );
202  dzErr = iBooker.book1D("dzErr", "PV tracks (p_{T} > 1 GeV) d_{z} error(cm)", 100, 0., 1. );
203 
204  dxyVsPhi_pt1 = iBooker.bookProfile("dxyVsPhi_pt1", "PV tracks (p_{T} > 1 GeV) d_{xy} (cm) VS track #phi",PhiBin, PhiMin, PhiMax, DxyBin, DxyMin, DxyMax,"");
205  dxyVsPhi_pt1->setAxisTitle("PV track (p_{T} > 1 GeV) #phi", 1);
206  dxyVsPhi_pt1->setAxisTitle("PV track (p_{T} > 1 GeV) d_{xy}",2);
207  dzVsPhi_pt1 = iBooker.bookProfile("dzVsPhi_pt1", "PV tracks (p_{T} > 1 GeV) d_{z} (cm) VS track #phi", PhiBin, PhiMin, PhiMax, DzBin, DzMin, DzMax, "");
208  dzVsPhi_pt1->setAxisTitle("PV track (p_{T} > 1 GeV) #phi", 1);
209  dzVsPhi_pt1->setAxisTitle("PV track (p_{T} > 1 GeV) d_{z}",2);
210 
211  dxyVsEta_pt1 = iBooker.bookProfile("dxyVsEta_pt1", "PV tracks (p_{T} > 1 GeV) d_{xy} (cm) VS track #eta",EtaBin, EtaMin, EtaMax, DxyBin, DxyMin, DxyMax,"");
212  dxyVsEta_pt1->setAxisTitle("PV track (p_{T} > 1 GeV) #eta", 1);
213  dxyVsEta_pt1->setAxisTitle("PV track (p_{T} > 1 GeV) d_{xy}",2);
214  dzVsEta_pt1 = iBooker.bookProfile("dzVsEta_pt1", "PV tracks (p_{T} > 1 GeV) d_{z} (cm) VS track #eta", EtaBin, EtaMin, EtaMax, DzBin, DzMin, DzMax, "");
215  dzVsEta_pt1->setAxisTitle("PV track (p_{T} > 1 GeV) #eta", 1);
216  dzVsEta_pt1->setAxisTitle("PV track (p_{T} > 1 GeV) d_{z}",2);
217 
218  dxyVsPhi_pt10 = iBooker.bookProfile("dxyVsPhi_pt10", "PV tracks (p_{T} > 1 GeV) d_{xy} (cm) VS track #phi",PhiBin, PhiMin, PhiMax, DxyBin, DxyMin, DxyMax,"");
219  dxyVsPhi_pt10->setAxisTitle("PV track (p_{T} > 10 GeV) #phi", 1);
220  dxyVsPhi_pt10->setAxisTitle("PV track (p_{T} > 10 GeV) d_{xy}",2);
221  dzVsPhi_pt10 = iBooker.bookProfile("dzVsPhi_pt10", "PV tracks (p_{T} > 10 GeV) d_{z} (cm) VS track #phi", PhiBin, PhiMin, PhiMax, DzBin, DzMin, DzMax, "");
222  dzVsPhi_pt10->setAxisTitle("PV track (p_{T} > 10 GeV) #phi", 1);
223  dzVsPhi_pt10->setAxisTitle("PV track (p_{T} > 10 GeV) d_{z}",2);
224 
225  dxyVsEta_pt10 = iBooker.bookProfile("dxyVsEta_pt10", "PV tracks (p_{T} > 10 GeV) d_{xy} (cm) VS track #eta",EtaBin, EtaMin, EtaMax, DxyBin, DxyMin, DxyMax,"");
226  dxyVsEta_pt10->setAxisTitle("PV track (p_{T} > 10 GeV) #eta", 1);
227  dxyVsEta_pt10->setAxisTitle("PV track (p_{T} > 10 GeV) d_{xy}",2);
228  dzVsEta_pt10 = iBooker.bookProfile("dzVsEta_pt10", "PV tracks (p_{T} > 10 GeV) d_{z} (cm) VS track #eta", EtaBin, EtaMin, EtaMax, DzBin, DzMin, DzMax, "");
229  dzVsEta_pt10->setAxisTitle("PV track (p_{T} > 10 GeV) #eta", 1);
230  dzVsEta_pt10->setAxisTitle("PV track (p_{T} > 10 GeV) d_{z}",2);
231 
232 }
233 
234 
235 PrimaryVertexMonitor::~PrimaryVertexMonitor()
236 {}
237 
238 void PrimaryVertexMonitor::analyze(const edm::Event& iEvent, const edm::EventSetup& iSetup)
239 {
240  Handle<reco::VertexCollection> recVtxs;
241  iEvent.getByToken(vertexToken_, recVtxs);
242 
243  Handle<VertexScore> scores;
244  iEvent.getByToken(scoreToken_, scores);
245 
246 
247  edm::Handle<reco::BeamSpot> beamSpotHandle;
248  iEvent.getByToken(beamspotToken_,beamSpotHandle);
249 
250  //
251  // check for absent products and simply "return" in that case
252  //
253  if (recVtxs.isValid() == false || beamSpotHandle.isValid()== false){
254  edm::LogWarning("PrimaryVertexMonitor")
255  <<" Some products not available in the event: VertexCollection "
256  <<vertexInputTag_<<" "
257  <<recVtxs.isValid() <<" BeamSpot "
258  <<beamSpotInputTag_<<" "
259  <<beamSpotHandle.isValid()<<". Skipping plots for this event";
260  return;
261  }
262 
263  BeamSpot beamSpot = *beamSpotHandle;
264 
265  nbvtx->Fill(recVtxs->size()*1.);
266  int ng=0;
267  for (auto const & vx : (*recVtxs) )
268  if (vx.isValid() && !vx.isFake() && vx.ndof()>=4.) ++ng;
269  nbgvtx->Fill(ng*1.);
270 
271  vertexPlots(recVtxs->front(), beamSpot, 1);
272 
273  if (scores.isValid() && (*scores).size()>0) {
274  auto pvScore = (*scores).get(0);
275  score[1]->Fill(std::sqrt(pvScore));
276  for (unsigned int i=1; i<(*scores).size(); ++i) score[0]->Fill(std::sqrt((*scores).get(i)));
277  }
278 
279  // fill PV tracks MEs (as now, for alignment)
280  pvTracksPlots(recVtxs->front());
281 
282  for(reco::VertexCollection::const_iterator v=recVtxs->begin()+1;
283  v!=recVtxs->end(); ++v){
284  vertexPlots(*v, beamSpot, 0);
285  }
286 
287  // Beamline plots:
288  bsX->Fill(beamSpot.x0());
289  bsY->Fill(beamSpot.y0());
290  bsZ->Fill(beamSpot.z0());
291  bsSigmaZ->Fill(beamSpot.sigmaZ());
292  bsDxdz->Fill(beamSpot.dxdz());
293  bsDydz->Fill(beamSpot.dydz());
294  bsBeamWidthX->Fill(beamSpot.BeamWidthX()*10000);
295  bsBeamWidthY->Fill(beamSpot.BeamWidthY()*10000);
296  // bsType->Fill(beamSpot.type());
297 
298 }
299 
300 void
301 PrimaryVertexMonitor::pvTracksPlots(const Vertex & v)
302 {
303 
304  const math::XYZPoint myVertex(v.position().x(),v.position().y(),v.position().z());
305 
306  if ( !v.isValid() ) return;
307  if ( v.isFake() ) return;
308 
309  if ( v.tracksSize() == 0 ) {
310  ntracks -> Fill ( 0 );
311  return;
312  }
313 
314  size_t nTracks = 0;
315  float sumPT = 0.;
316  for (reco::Vertex::trackRef_iterator t = v.tracks_begin(); t != v.tracks_end(); t++) {
317 
318  bool isHighPurity = (**t).quality(reco::TrackBase::highPurity);
319  if ( !isHighPurity ) continue;
320 
321  float pt = (**t).pt();
322  if ( pt < 1. ) continue;
323 
324  nTracks++;
325 
326  float eta = (**t).eta();
327  float phi = (**t).phi();
328 
329  float w = v.trackWeight(*t);
330  float chi2NDF = (**t).normalizedChi2();
331  float chi2Prob = TMath::Prob((**t).chi2(),(int)(**t).ndof());
332  float Dxy = (**t).dxy(myVertex);
333  float Dz = (**t).dz(myVertex);
334  float DxyErr = (**t).dxyError();
335  float DzErr = (**t).dzError();
336 
337  sumPT += pt*pt;
338 
339 
340  // fill MEs
341  weight -> Fill (w);
342  chi2ndf -> Fill (chi2NDF);
343  chi2prob -> Fill (chi2Prob);
344  dxy -> Fill (Dxy);
345  dxy2 -> Fill (Dxy);
346  dz -> Fill (Dz);
347  dxyErr -> Fill (DxyErr);
348  dzErr -> Fill (DzErr);
349 
350  dxyVsPhi_pt1 -> Fill (phi,Dxy);
351  dzVsPhi_pt1 -> Fill (phi,Dz);
352  dxyVsEta_pt1 -> Fill (eta,Dxy);
353  dzVsEta_pt1 -> Fill (eta,Dz);
354 
355  if ( pt < 10. ) continue;
356  dxyVsPhi_pt10 -> Fill (phi,Dxy);
357  dzVsPhi_pt10 -> Fill (phi,Dz);
358  dxyVsEta_pt10 -> Fill (eta,Dxy);
359  dzVsEta_pt10 -> Fill (eta,Dz);
360  }
361  ntracks -> Fill (float(nTracks));
362  sumpt -> Fill (sumPT);
363 
364 }
365 
366 void PrimaryVertexMonitor::vertexPlots(const Vertex & v, const BeamSpot& beamSpot, int i)
367 {
368 
369  if (i < 0 || i > 1) return;
370  if (!v.isValid()) type[i]->Fill(2.);
371  else if (v.isFake()) type[i]->Fill(1.);
372  else type[i]->Fill(0.);
373 
374  if (v.isValid() && !v.isFake()) {
375  float weight = 0;
376  for(reco::Vertex::trackRef_iterator t = v.tracks_begin();
377  t!=v.tracks_end(); t++) weight+= v.trackWeight(*t);
378  trksWeight[i]->Fill(weight);
379  nbtksinvtx[i]->Fill(v.tracksSize());
380  ntracksVsZ[i]->Fill(v.position().z()- beamSpot.z0(),v.tracksSize());
381 
382  vtxchi2[i]->Fill(v.chi2());
383  vtxndf[i]->Fill(v.ndof());
384  vtxprob[i]->Fill(ChiSquaredProbability(v.chi2() ,v.ndof()));
385 
386  xrec[i]->Fill(v.position().x());
387  yrec[i]->Fill(v.position().y());
388  zrec[i]->Fill(v.position().z());
389 
390  float xb = beamSpot.x0() + beamSpot.dxdz() * (v.position().z() - beamSpot.z0());
391  float yb = beamSpot.y0() + beamSpot.dydz() * (v.position().z() - beamSpot.z0());
392  xDiff[i]->Fill((v.position().x() - xb)*10000);
393  yDiff[i]->Fill((v.position().y() - yb)*10000);
394 
395  xerr[i]->Fill(v.xError()*10000);
396  yerr[i]->Fill(v.yError()*10000);
397  zerr[i]->Fill(v.zError()*10000);
398  xerrVsTrks[i]->Fill(weight, v.xError()*10000);
399  yerrVsTrks[i]->Fill(weight, v.yError()*10000);
400  zerrVsTrks[i]->Fill(weight, v.zError()*10000);
401 
402  nans[i]->Fill(1.,edm::isNotFinite(v.position().x())*1.);
403  nans[i]->Fill(2.,edm::isNotFinite(v.position().y())*1.);
404  nans[i]->Fill(3.,edm::isNotFinite(v.position().z())*1.);
405 
406  int index = 3;
407  for (int k = 0; k != 3; k++) {
408  for (int j = k; j != 3; j++) {
409  index++;
410  nans[i]->Fill(index*1., edm::isNotFinite(v.covariance(k, j))*1.);
411  // in addition, diagonal element must be positive
412  if (j == k && v.covariance(k, j) < 0) {
413  nans[i]->Fill(index*1., 1.);
414  }
415  }
416  }
417  }
418 }
419 
420 //define this as a plug-in
421 DEFINE_FWK_MODULE(PrimaryVertexMonitor);
type
type
Definition: HCALResponse.h:21
PrimaryVertexMonitor::xerrVsTrks
MonitorElement * xerrVsTrks[2]
Definition: PrimaryVertexMonitor.h:59
edm::ParameterSet::getParameter
T getParameter(std::string const &) const
reco::BeamSpot::z0
double z0() const
z coordinate
Definition: BeamSpot.h:68
PrimaryVertexMonitor::xerr
MonitorElement * xerr[2]
Definition: PrimaryVertexMonitor.h:58
i
int i
Definition: DBlmapReader.cc:9
nTracks
const unsigned int nTracks(const reco::Vertex &sv)
Definition: TemplatedVertexArbitrator.h:48
PrimaryVertexMonitor::vtxndf
MonitorElement * vtxndf[2]
Definition: PrimaryVertexMonitor.h:61
PrimaryVertexMonitor::dxyVsEta_pt1
MonitorElement * dxyVsEta_pt1
Definition: PrimaryVertexMonitor.h:68
SiPixelRawToDigiRegional_cfi.beamSpot
tuple beamSpot
Definition: SiPixelRawToDigiRegional_cfi.py:11
PrimaryVertexMonitor::dzVsPhi_pt10
MonitorElement * dzVsPhi_pt10
Definition: PrimaryVertexMonitor.h:69
PrimaryVertexMonitor::dzVsPhi_pt1
MonitorElement * dzVsPhi_pt1
Definition: PrimaryVertexMonitor.h:67
w
const double w
Definition: UKUtility.cc:23
PrimaryVertexMonitor::pvTracksPlots
void pvTracksPlots(const reco::Vertex &v)
Definition: PrimaryVertexMonitor.cc:301
reco::Vertex::tracks_end
trackRef_iterator tracks_end() const
last iterator over tracks
Definition: Vertex.cc:44
DQMStore::IBooker::bookProfile
MonitorElement * bookProfile(Args &&...args)
Definition: DQMStore.h:157
PrimaryVertexMonitor::vertexPlots
void vertexPlots(const reco::Vertex &v, const reco::BeamSpot &beamSpot, int i)
Definition: PrimaryVertexMonitor.cc:366
edm::Event::getByToken
bool getByToken(EDGetToken token, Handle< PROD > &result) const
Definition: Event.h:462
reco::Vertex::zError
double zError() const
error on z
Definition: Vertex.h:111
DEFINE_FWK_MODULE
#define DEFINE_FWK_MODULE(type)
Definition: MakerMacros.h:17
PrimaryVertexMonitor::bsSigmaZ
MonitorElement * bsSigmaZ
Definition: PrimaryVertexMonitor.h:63
PrimaryVertexMonitor::beamSpotInputTag_
edm::InputTag beamSpotInputTag_
Definition: PrimaryVertexMonitor.h:46
reco::Vertex::isValid
bool isValid() const
Tells whether the vertex is valid.
Definition: Vertex.h:60
relval_2017.k
list k
Definition: relval_2017.py:23
PrimaryVertexMonitor::scoreToken_
edm::EDGetTokenT< VertexScore > scoreToken_
Definition: PrimaryVertexMonitor.h:44
NULL
#define NULL
Definition: scimark2.h:8
PrimaryVertexMonitor::chi2ndf
MonitorElement * chi2ndf
Definition: PrimaryVertexMonitor.h:65
reco::Vertex::covariance
double covariance(int i, int j) const
(i, j)-th element of error matrix, i, j = 0, ... 2
Definition: Vertex.h:116
PrimaryVertexMonitor::dxyVsPhi_pt10
MonitorElement * dxyVsPhi_pt10
Definition: PrimaryVertexMonitor.h:69
PrimaryVertexMonitor::xDiff
MonitorElement * xDiff[2]
Definition: PrimaryVertexMonitor.h:58
dt_dqm_sourceclient_common_cff.reco
tuple reco
Definition: dt_dqm_sourceclient_common_cff.py:107
reco::Vertex::position
const Point & position() const
position
Definition: Vertex.h:99
PrimaryVertexMonitor::dzVsEta_pt1
MonitorElement * dzVsEta_pt1
Definition: PrimaryVertexMonitor.h:68
AlCaHLTBitMon_QueryRunRegistry.string
string string
Definition: AlCaHLTBitMon_QueryRunRegistry.py:255
PrimaryVertexMonitor::ntracksVsZ
MonitorElement * ntracksVsZ[2]
Definition: PrimaryVertexMonitor.h:60
PrimaryVertexMonitor::conf_
edm::ParameterSet conf_
Definition: PrimaryVertexMonitor.h:48
MonitorElement::Fill
void Fill(long long x)
Definition: MonitorElement.h:170
PrimaryVertexMonitor::zerrVsTrks
MonitorElement * zerrVsTrks[2]
Definition: PrimaryVertexMonitor.h:59
PrimaryVertexMonitor::nans
MonitorElement * nans[2]
Definition: PrimaryVertexMonitor.h:61
reco::BeamSpot::dydz
double dydz() const
dydz slope
Definition: BeamSpot.h:84
iEvent
int iEvent
Definition: GenABIO.cc:230
edm::isNotFinite
bool isNotFinite(T x)
Definition: isFinite.h:10
PrimaryVertexMonitor::yerrVsTrks
MonitorElement * yerrVsTrks[2]
Definition: PrimaryVertexMonitor.h:59
PrimaryVertexMonitor::bookHistograms
virtual void bookHistograms(DQMStore::IBooker &, edm::Run const &, edm::EventSetup const &) override
Definition: PrimaryVertexMonitor.cc:63
PrimaryVertexMonitor::dxyVsEta_pt10
MonitorElement * dxyVsEta_pt10
Definition: PrimaryVertexMonitor.h:70
HcalObjRepresent::Fill
void Fill(HcalDetId &id, double val, std::vector< TH2F > &depth)
Definition: HcalObjRepresent.h:611
PrimaryVertexMonitor::zerr
MonitorElement * zerr[2]
Definition: PrimaryVertexMonitor.h:58
mathSSE::sqrt
T sqrt(T t)
Definition: SSEVec.h:18
HLT_25ns10e33_v2_cff.EtaMax
tuple EtaMax
Definition: HLT_25ns10e33_v2_cff.py:2840
PrimaryVertexMonitor::yrec
MonitorElement * yrec[2]
Definition: PrimaryVertexMonitor.h:58
DQMStore::IBooker::book1D
MonitorElement * book1D(Args &&...args)
Definition: DQMStore.h:115
reco::Vertex::chi2
double chi2() const
chi-squares
Definition: Vertex.h:88
j
int j
Definition: DBlmapReader.cc:9
ChiSquaredProbability
float ChiSquaredProbability(double chiSquared, double nrDOF)
Definition: ChiSquaredProbability.cc:13
reco::BeamSpot::BeamWidthX
double BeamWidthX() const
beam width X
Definition: BeamSpot.h:86
PrimaryVertexMonitor::zrec
MonitorElement * zrec[2]
Definition: PrimaryVertexMonitor.h:58
Gflash::EtaMin
const double EtaMin[kNumberCalorimeter]
Definition: GflashNameSpace.h:32
reco::Vertex::trackWeight
float trackWeight(const TREF &r) const
returns the weight with which a Track has contributed to the vertex-fit.
Definition: Vertex.h:73
PrimaryVertexMonitor::vtxchi2
MonitorElement * vtxchi2[2]
Definition: PrimaryVertexMonitor.h:61
PrimaryVertexMonitor::beamspotToken_
edm::EDGetTokenT< reco::BeamSpot > beamspotToken_
Definition: PrimaryVertexMonitor.h:42
PrimaryVertexMonitor::vertexToken_
edm::EDGetTokenT< reco::VertexCollection > vertexToken_
Definition: PrimaryVertexMonitor.h:41
edm::HandleBase::isValid
bool isValid() const
Definition: HandleBase.h:75
reco::BeamSpot::dxdz
double dxdz() const
dxdz slope
Definition: BeamSpot.h:82
reco::Vertex::ndof
double ndof() const
Definition: Vertex.h:95
PrimaryVertexMonitor::dzVsEta_pt10
MonitorElement * dzVsEta_pt10
Definition: PrimaryVertexMonitor.h:70
reco::Vertex::xError
double xError() const
error on x
Definition: Vertex.h:107
DQMStore::IBooker::setCurrentFolder
void setCurrentFolder(const std::string &fullpath)
Definition: DQMStore.cc:276
reco::Vertex::isFake
bool isFake() const
Definition: Vertex.h:64
PrimaryVertexMonitor::chi2prob
MonitorElement * chi2prob
Definition: PrimaryVertexMonitor.h:65
DQMStore::IBooker::book2D
MonitorElement * book2D(Args &&...args)
Definition: DQMStore.h:133
math::XYZPoint
XYZPointD XYZPoint
point in space with cartesian internal representation
Definition: Point3D.h:12
PrimaryVertexMonitor::nbtksinvtx
MonitorElement * nbtksinvtx[2]
Definition: PrimaryVertexMonitor.h:56
reco::BeamSpot::sigmaZ
double sigmaZ() const
sigma z
Definition: BeamSpot.h:80
MonitorElement::getTH1F
TH1F * getTH1F(void) const
Definition: MonitorElement.cc:1457
reco::BeamSpot::BeamWidthY
double BeamWidthY() const
beam width Y
Definition: BeamSpot.h:88
edm::InputTag::label
std::string const & label() const
Definition: InputTag.h:36
PrimaryVertexMonitor::yerr
MonitorElement * yerr[2]
Definition: PrimaryVertexMonitor.h:58
PrimaryVertexMonitor::analyze
virtual void analyze(const edm::Event &, const edm::EventSetup &) override
Definition: PrimaryVertexMonitor.cc:238
reco::Vertex::trackRef_iterator
std::vector< TrackBaseRef >::const_iterator trackRef_iterator
The iteratator for the vector&lt;TrackRef&gt;
Definition: Vertex.h:37
reco::BeamSpot::y0
double y0() const
y coordinate
Definition: BeamSpot.h:66
PrimaryVertexMonitor::xrec
MonitorElement * xrec[2]
Definition: PrimaryVertexMonitor.h:58
PrimaryVertexMonitor::vtxprob
MonitorElement * vtxprob[2]
Definition: PrimaryVertexMonitor.h:61
PrimaryVertexMonitor::ntracks
MonitorElement * ntracks
Definition: PrimaryVertexMonitor.h:65
reco::Vertex::tracks_begin
trackRef_iterator tracks_begin() const
first iterator over tracks
Definition: Vertex.cc:39
PrimaryVertexMonitor::bsBeamWidthY
MonitorElement * bsBeamWidthY
Definition: PrimaryVertexMonitor.h:63
histoStyle.weight
int weight
Definition: histoStyle.py:50
PrimaryVertexMonitor::score
MonitorElement * score[2]
Definition: PrimaryVertexMonitor.h:56
MonitorElement::setAxisTitle
void setAxisTitle(const std::string &title, int axis=1)
set x-, y- or z-axis title (axis=1, 2, 3 respectively)
Definition: MonitorElement.cc:947
PrimaryVertexMonitor::yDiff
MonitorElement * yDiff[2]
Definition: PrimaryVertexMonitor.h:58
reco::Vertex::yError
double yError() const
error on y
Definition: Vertex.h:109
PrimaryVertexMonitor::trksWeight
MonitorElement * trksWeight[2]
Definition: PrimaryVertexMonitor.h:56
reco::Vertex::tracksSize
size_t tracksSize() const
number of tracks
Definition: Vertex.cc:34
PrimaryVertexMonitor::PrimaryVertexMonitor
PrimaryVertexMonitor(const edm::ParameterSet &pSet)
Definition: PrimaryVertexMonitor.cc:16
edm::Run
Definition: Run.h:43
PrimaryVertexMonitor::bsBeamWidthX
MonitorElement * bsBeamWidthX
Definition: PrimaryVertexMonitor.h:63
PrimaryVertexMonitor::dxyVsPhi_pt1
MonitorElement * dxyVsPhi_pt1
Definition: PrimaryVertexMonitor.h:67
reco::BeamSpot::x0
double x0() const
x coordinate
Definition: BeamSpot.h:64
Generated for CMSSW Reference Manual by doxygen 1.8.5