CMS 3D CMS Logo
List of all members | Public Member Functions | Private Member Functions | Private Attributes
recoBSVTagInfoValidationAnalyzer Class Reference
Inheritance diagram for recoBSVTagInfoValidationAnalyzer:
edm::EDAnalyzer edm::EDConsumerBase

Public Member Functions

 recoBSVTagInfoValidationAnalyzer (const edm::ParameterSet &)
 
- Public Member Functions inherited from edm::EDAnalyzer
void callWhenNewProductsRegistered (std::function< void(BranchDescription const &)> const &func)
 
 EDAnalyzer ()
 
SerialTaskQueue * globalLuminosityBlocksQueue ()
 
SerialTaskQueue * globalRunsQueue ()
 
ModuleDescription const & moduleDescription () const
 
std::string workerType () const
 
 ~EDAnalyzer () override
 
- Public Member Functions inherited from edm::EDConsumerBase
std::vector< ConsumesInfo > consumesInfo () const
 
void convertCurrentProcessAlias (std::string const &processName)
 Convert "@currentProcess" in InputTag process names to the actual current process name. More...
 
 EDConsumerBase ()
 
 EDConsumerBase (EDConsumerBase const &)=delete
 
 EDConsumerBase (EDConsumerBase &&)=default
 
ESProxyIndex const * esGetTokenIndices (edm::Transition iTrans) const
 
ProductResolverIndexAndSkipBit indexFrom (EDGetToken, BranchType, TypeID const &) const
 
void itemsMayGet (BranchType, std::vector< ProductResolverIndexAndSkipBit > &) const
 
void itemsToGet (BranchType, std::vector< ProductResolverIndexAndSkipBit > &) const
 
std::vector< ProductResolverIndexAndSkipBit > const & itemsToGetFrom (BranchType iType) const
 
void labelsForToken (EDGetToken iToken, Labels &oLabels) const
 
void modulesWhoseProductsAreConsumed (std::vector< ModuleDescription const * > &modules, ProductRegistry const &preg, std::map< std::string, ModuleDescription const * > const &labelsToDesc, std::string const &processName) const
 
EDConsumerBase const & operator= (EDConsumerBase const &)=delete
 
EDConsumerBase & operator= (EDConsumerBase &&)=default
 
bool registeredToConsume (ProductResolverIndex, bool, BranchType) const
 
bool registeredToConsumeMany (TypeID const &, BranchType) const
 
ProductResolverIndexAndSkipBit uncheckedIndexFrom (EDGetToken) const
 
void updateLookup (BranchType iBranchType, ProductResolverIndexHelper const &, bool iPrefetchMayGet)
 
void updateLookup (eventsetup::ESRecordsToProxyIndices const &)
 
virtual ~EDConsumerBase () noexcept(false)
 

Private Member Functions

void analyze (const edm::Event &, const edm::EventSetup &) override
 
void bookRecoToSim (std::string const &)
 
void bookSimToReco (std::string const &)
 
void endJob () override
 
void fillRecoToSim (std::string const &, reco::Vertex const &, TrackingVertexRef const &)
 
void fillSimToReco (std::string const &, reco::VertexBaseRef const &, TrackingVertexRef const &)
 

Private Attributes

VertexClassifierByProxy< reco::SecondaryVertexTagInfoCollection > classifier_
 
std::string dqmLabel
 
DQMStore * dqmStore_
 
std::map< std::string, MonitorElement * > HistIndex_
 
Int_t n_event
 
Int_t numberVertexClassifier_
 
Int_t rs_total_nall
 
Int_t rs_total_nbsv
 
Int_t rs_total_nbv
 
Int_t rs_total_ncv
 
Int_t rs_total_nlv
 
Int_t rs_total_nsv
 
Int_t sr_total_nall
 
Int_t sr_total_nbsv
 
Int_t sr_total_nbv
 
Int_t sr_total_ncv
 
Int_t sr_total_nlv
 
Int_t sr_total_nsv
 
edm::EDGetTokenT< reco::SecondaryVertexTagInfoCollection > svInfoToken
 
edm::InputTag svTagInfoProducer_
 
Int_t total_nfake
 
Int_t total_nmiss
 
edm::InputTag trackingTruth_
 
edm::EDGetTokenT< TrackingVertexCollection > tvToken
 

Additional Inherited Members

- Public Types inherited from edm::EDAnalyzer
typedef EDAnalyzer ModuleType
 
- Public Types inherited from edm::EDConsumerBase
typedef ProductLabels Labels
 
- Static Public Member Functions inherited from edm::EDAnalyzer
static const std::string & baseType ()
 
static void fillDescriptions (ConfigurationDescriptions &descriptions)
 
static void prevalidate (ConfigurationDescriptions &)
 
static bool wantsGlobalLuminosityBlocks ()
 
static bool wantsGlobalRuns ()
 
static bool wantsStreamLuminosityBlocks ()
 
static bool wantsStreamRuns ()
 
- Protected Member Functions inherited from edm::EDConsumerBase
template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType > consumes (edm::InputTag const &tag)
 
EDGetToken consumes (const TypeToGet &id, edm::InputTag const &tag)
 
template<BranchType B>
EDGetToken consumes (TypeToGet const &id, edm::InputTag const &tag)
 
ConsumesCollector consumesCollector ()
 Use a ConsumesCollector to gather consumes information from helper functions. More...
 
template<typename ProductType , BranchType B = InEvent>
void consumesMany ()
 
void consumesMany (const TypeToGet &id)
 
template<BranchType B>
void consumesMany (const TypeToGet &id)
 
template<typename ESProduct , typename ESRecord , Transition Tr = Transition::Event>
auto esConsumes ()
 
template<typename ESProduct , typename ESRecord , Transition Tr = Transition::Event>
auto esConsumes (ESInputTag const &tag)
 
template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType > mayConsume (edm::InputTag const &tag)
 
EDGetToken mayConsume (const TypeToGet &id, edm::InputTag const &tag)
 
template<BranchType B>
EDGetToken mayConsume (const TypeToGet &id, edm::InputTag const &tag)
 

Detailed Description

Definition at line 37 of file recoBSVTagInfoValidationAnalyzer.cc.

Constructor & Destructor Documentation

recoBSVTagInfoValidationAnalyzer::recoBSVTagInfoValidationAnalyzer ( const edm::ParameterSet &  config)
explicit

Definition at line 88 of file recoBSVTagInfoValidationAnalyzer.cc.

References bookRecoToSim(), bookSimToReco(), dqmLabel, dqmStore_, edm::ParameterSet::getParameter(), HistIndex_, mps_fire::i, n_event, VertexCategories::Names, numberVertexClassifier_, Utilities::operator, rs_total_nall, rs_total_nbsv, rs_total_nbv, rs_total_ncv, rs_total_nlv, rs_total_nsv, sr_total_nall, sr_total_nbsv, sr_total_nbv, sr_total_ncv, sr_total_nlv, sr_total_nsv, svInfoToken, total_nfake, total_nmiss, tvToken, and VertexCategories::Unknown.

89  : classifier_(config, consumesCollector()) {
90  // Initialize counters
91  n_event = 0;
92  rs_total_nall = 0;
93  rs_total_nsv = 0;
94  rs_total_nbv = 0;
95  rs_total_nbsv = 0;
96  rs_total_ncv = 0;
97  rs_total_nlv = 0;
98  total_nfake = 0;
99 
100  sr_total_nall = 0;
101  sr_total_nsv = 0;
102  sr_total_nbv = 0;
103  sr_total_nbsv = 0;
104  sr_total_ncv = 0;
105  sr_total_nlv = 0;
106  total_nmiss = 0;
107 
108  // get the store
109  dqmStore_ = edm::Service<DQMStore>().operator->();
110  dqmLabel = "SVValidation/";
111  dqmStore_->setCurrentFolder(dqmLabel);
112 
113  // Get the track collection
114  svInfoToken = consumes<reco::SecondaryVertexTagInfoCollection>(config.getParameter<InputTag>("svTagInfoProducer"));
115  // Name of the traking pariticle collection
116  tvToken = consumes<TrackingVertexCollection>(config.getParameter<InputTag>("trackingTruth"));
117  // Number of track categories
118  numberVertexClassifier_ = VertexCategories::Unknown + 1;
119 
120  // Define histogram for counting categories
121  HistIndex_["VertexClassifier"] = dqmStore_->book1D("VertexClassifier",
122  "Frequency for the different track categories",
123  numberVertexClassifier_,
124  -0.5,
125  numberVertexClassifier_ - 0.5);
126 
127  //--- RecoToSim
128  HistIndex_["rs_All_MatchQuality"] = dqmStore_->book1D("rs_All_MatchQuality", "Quality of Match", 51, -0.01, 1.01);
129  HistIndex_["rs_All_FlightDistance2d"] =
130  dqmStore_->book1D("rs_All_FlightDistance2d", "Transverse flight distance [cm]", 100, 0, 5);
131  HistIndex_["rs_SecondaryVertex_FlightDistance2d"] =
132  dqmStore_->book1D("rs_SecondaryVertex_FlightDistance2d", "Transverse flight distance [cm]", 100, 0, 5);
133  HistIndex_["rs_BSV_FlightDistance2d"] =
134  dqmStore_->book1D("rs_BSV_FlightDistance2d", "Transverse flight distance [cm]", 100, 0, 5);
135  HistIndex_["rs_BWeakDecay_FlightDistance2d"] =
136  dqmStore_->book1D("rs_BWeakDecay_FlightDistance2d", "Transverse flight distance [cm]", 100, 0, 5);
137  HistIndex_["rs_CWeakDecay_FlightDistance2d"] =
138  dqmStore_->book1D("rs_CWeakDecay_FlightDistance2d", "Transverse flight distance [cm]", 100, 0, 5);
139  HistIndex_["rs_Light_FlightDistance2d"] =
140  dqmStore_->book1D("rs_Light_FlightDistance2d", "Transverse flight distance [cm]", 100, 0, 5);
141 
142  HistIndex_["rs_All_nRecVtx"] = dqmStore_->book1D("rs_All_nRecVtx", "Number of Vertices per event", 11, -0.5, 10.5);
143  HistIndex_["rs_SecondaryVertex_nRecVtx"] =
144  dqmStore_->book1D("rs_SecondaryVertex_nRecVtx", "Number of Vertices per event", 11, -0.5, 10.5);
145  HistIndex_["rs_BSV_nRecVtx"] = dqmStore_->book1D("rs_BSV_nRecVtx", "Number of Vertices per event", 11, -0.5, 10.5);
146  HistIndex_["rs_BWeakDecay_nRecVtx"] =
147  dqmStore_->book1D("rs_BWeakDecay_nRecVtx", "Number of Vertices per event", 11, -0.5, 10.5);
148  HistIndex_["rs_CWeakDecay_nRecVtx"] =
149  dqmStore_->book1D("rs_CWeakDecay_nRecVtx", "Number of Vertices per event", 11, -0.5, 10.5);
150  HistIndex_["rs_Light_nRecVtx"] =
151  dqmStore_->book1D("rs_Light_nRecVtx", "Number of Vertices per event", 11, -0.5, 10.5);
152 
153  //--- SimToReco
154  HistIndex_["sr_All_MatchQuality"] = dqmStore_->book1D("sr_All_MatchQuality", "Quality of Match", 51, -0.01, 1.01);
155  HistIndex_["sr_All_nRecVtx"] = dqmStore_->book1D("sr_All_nRecVtx", "Number of Vertices per event", 11, -0.5, 10.5);
156  HistIndex_["sr_SecondaryVertex_nRecVtx"] =
157  dqmStore_->book1D("sr_SecondaryVertex_nRecVtx", "Number of Vertices per event", 11, -0.5, 10.5);
158  HistIndex_["sr_BSV_nRecVtx"] = dqmStore_->book1D("sr_BSV_nRecVtx", "Number of Vertices per event", 11, -0.5, 10.5);
159  HistIndex_["sr_BWeakDecay_nRecVtx"] =
160  dqmStore_->book1D("sr_BWeakDecay_nRecVtx", "Number of Vertices per event", 11, -0.5, 10.5);
161  HistIndex_["sr_CWeakDecay_nRecVtx"] =
162  dqmStore_->book1D("sr_CWeakDecay_nRecVtx", "Number of Vertices per event", 11, -0.5, 10.5);
163  HistIndex_["sr_Light_nRecVtx"] =
164  dqmStore_->book1D("sr_Light_nRecVtx", "Number of Vertices per event", 11, -0.5, 10.5);
165 
166  // Set the proper categories names
167  for (Int_t i = 0; i < numberVertexClassifier_; ++i)
168  HistIndex_["VertexClassifier"]->getTH1F()->GetXaxis()->SetBinLabel(i + 1, VertexCategories::Names[i]);
169 
170  // book histograms
171  bookRecoToSim("rs_All");
172  bookRecoToSim("rs_SecondaryVertex");
173  bookRecoToSim("rs_BSV");
174  bookRecoToSim("rs_BWeakDecay");
175  bookRecoToSim("rs_CWeakDecay");
176  bookRecoToSim("rs_Light");
177 
178  bookSimToReco("sr_All");
179  bookSimToReco("sr_SecondaryVertex");
180  bookSimToReco("sr_BSV");
181  bookSimToReco("sr_BWeakDecay");
182  bookSimToReco("sr_CWeakDecay");
183  bookSimToReco("sr_Light");
184 }
edm::ParameterSet::getParameter
T getParameter(std::string const &) const
mps_fire.i
i
Definition: mps_fire.py:338
edm::EDConsumerBase::consumesCollector
ConsumesCollector consumesCollector()
Use a ConsumesCollector to gather consumes information from helper functions.
Definition: EDConsumerBase.cc:71
recoBSVTagInfoValidationAnalyzer::classifier_
VertexClassifierByProxy< reco::SecondaryVertexTagInfoCollection > classifier_
Definition: recoBSVTagInfoValidationAnalyzer.cc:46
recoBSVTagInfoValidationAnalyzer::tvToken
edm::EDGetTokenT< TrackingVertexCollection > tvToken
Definition: recoBSVTagInfoValidationAnalyzer.cc:85
recoBSVTagInfoValidationAnalyzer::HistIndex_
std::map< std::string, MonitorElement * > HistIndex_
Definition: recoBSVTagInfoValidationAnalyzer.cc:81
VertexCategories::Names
static const char *const Names[]
Name of the different categories.
Definition: VertexCategories.h:52
recoBSVTagInfoValidationAnalyzer::svInfoToken
edm::EDGetTokenT< reco::SecondaryVertexTagInfoCollection > svInfoToken
Definition: recoBSVTagInfoValidationAnalyzer.cc:84

Member Function Documentation

void recoBSVTagInfoValidationAnalyzer::analyze ( const edm::Event &  event,
const edm::EventSetup &  setup 
)
overrideprivate

Definition at line 186 of file recoBSVTagInfoValidationAnalyzer.cc.

References VertexCategories::BWeakDecay, classifier_, gather_cfg::cout, VertexCategories::CWeakDecay, VertexClassifierByProxy< Collection >::evaluate(), VertexCategories::Fake, fillRecoToSim(), fillSimToReco(), HistIndex_, VertexClassifier::history(), mps_fire::i, VertexCategories::is(), n_event, VertexClassifierByProxy< Collection >::newEvent(), numberVertexClassifier_, VertexHistory::quality(), VertexCategories::Reconstructed, VertexHistory::recoVertex(), rs_total_nall, rs_total_nbsv, rs_total_nbv, rs_total_ncv, rs_total_nlv, rs_total_nsv, VertexCategories::SecondaryVertex, HistoryBase::simVertex(), sr_total_nall, sr_total_nbsv, sr_total_nbv, sr_total_ncv, sr_total_nlv, sr_total_nsv, svInfoToken, total_nfake, total_nmiss, and tvToken.

186  {
187  ++n_event;
188 
189  std::cout << "*** Analyzing " << event.id() << " n_event = " << n_event << std::endl << std::endl;
190 
191  // Set the classifier for a new event
192  classifier_.newEvent(event, setup);
193 
194  // Vertex collection
195  edm::Handle<reco::SecondaryVertexTagInfoCollection> svTagInfoCollection;
196  event.getByToken(svInfoToken, svTagInfoCollection);
197  // Get a constant reference to the track history associated to the classifier
198  VertexHistory const &tracer = classifier_.history();
199 
200  cout << "* Event " << n_event << " ; svTagInfoCollection->size() = " << svTagInfoCollection->size() << endl;
201 
202  int rs_nall = 0;
203  int rs_nsv = 0;
204  int rs_nbv = 0;
205  int rs_nbsv = 0;
206  int rs_ncv = 0;
207  int rs_nlv = 0;
208  int nfake = 0;
209 
210  int sr_nall = 0;
211  int sr_nsv = 0;
212  int sr_nbv = 0;
213  int sr_nbsv = 0;
214  int sr_ncv = 0;
215  int sr_nlv = 0;
216  int nmiss = 0;
217 
218  // Loop over the svTagInfo collection.
219  for (std::size_t index = 0; index < svTagInfoCollection->size(); ++index) {
220  reco::SecondaryVertexTagInfoRef svTagInfo(svTagInfoCollection, index);
221 
222  // Loop over the vertexes in svTagInfo
223  for (std::size_t vindex = 0; vindex < svTagInfo->nVertices(); ++vindex) {
224  // Classify the vertices
225  classifier_.evaluate(svTagInfo, vindex);
226 
227  // quality of the match
228  double rs_quality = tracer.quality();
229 
230  // Fill the histogram with the categories
231  for (Int_t i = 0; i != numberVertexClassifier_; ++i) {
232  if (classifier_.is((VertexCategories::Category)i)) {
233  HistIndex_["VertexClassifier"]->Fill(i);
234  }
235  }
236  if (!classifier_.is(VertexCategories::Fake)) {
237  HistIndex_["rs_All_MatchQuality"]->Fill(rs_quality);
238  fillRecoToSim("rs_All", svTagInfo->secondaryVertex(vindex), tracer.simVertex());
239  HistIndex_["rs_All_FlightDistance2d"]->Fill(svTagInfo->flightDistance(vindex, true).value());
240  rs_nall++;
241 
242  if (classifier_.is(VertexCategories::SecondaryVertex)) {
243  fillRecoToSim("rs_SecondaryVertex", svTagInfo->secondaryVertex(vindex), tracer.simVertex());
244  HistIndex_["rs_SecondaryVertex_FlightDistance2d"]->Fill(svTagInfo->flightDistance(vindex, true).value());
245  rs_nsv++;
246  }
247 
248  if (classifier_.is(VertexCategories::BWeakDecay)) {
249  fillRecoToSim("rs_BWeakDecay", svTagInfo->secondaryVertex(vindex), tracer.simVertex());
250  HistIndex_["rs_BWeakDecay_FlightDistance2d"]->Fill(svTagInfo->flightDistance(vindex, true).value());
251  rs_nbv++;
252 
253  if (classifier_.is(VertexCategories::SecondaryVertex)) {
254  fillRecoToSim("rs_BSV", svTagInfo->secondaryVertex(vindex), tracer.simVertex());
255  HistIndex_["rs_BSV_FlightDistance2d"]->Fill(svTagInfo->flightDistance(vindex, true).value());
256  rs_nbsv++;
257  }
258  } // BWeakDecay
259 
260  else if (classifier_.is(VertexCategories::CWeakDecay)) {
261  fillRecoToSim("rs_CWeakDecay", svTagInfo->secondaryVertex(vindex), tracer.simVertex());
262  HistIndex_["rs_CWeakDecay_FlightDistance2d"]->Fill(svTagInfo->flightDistance(vindex, true).value());
263  rs_ncv++;
264 
265  } else {
266  fillRecoToSim("rs_Light", svTagInfo->secondaryVertex(vindex), tracer.simVertex());
267  HistIndex_["rs_Light_FlightDistance2d"]->Fill(svTagInfo->flightDistance(vindex, true).value());
268  rs_nlv++;
269  }
270  } // end if classifier
271 
272  else {
273  cout << " VertexCategories::Fake!!" << endl;
274  nfake++;
275  }
276 
277  } // end loop over vertices in svTagInfo
278 
279  } // loop over svTagInfo
280 
281  HistIndex_["rs_All_nRecVtx"]->Fill(rs_nall);
282  HistIndex_["rs_SecondaryVertex_nRecVtx"]->Fill(rs_nsv);
283  HistIndex_["rs_BWeakDecay_nRecVtx"]->Fill(rs_nbv);
284  HistIndex_["rs_BSV_nRecVtx"]->Fill(rs_nbsv);
285  HistIndex_["rs_CWeakDecay_nRecVtx"]->Fill(rs_ncv);
286  HistIndex_["rs_Light_nRecVtx"]->Fill(rs_nlv);
287  cout << endl;
288 
289  //----------------------------------------------------------------
290  // SIM TO RECO!
291 
292  // Vertex collection
293  edm::Handle<TrackingVertexCollection> TVCollection;
294  event.getByToken(tvToken, TVCollection);
295  // Loop over the TV collection.
296  for (std::size_t index = 0; index < TVCollection->size(); ++index) {
297  TrackingVertexRef trackingVertex(TVCollection, index);
298 
299  classifier_.evaluate(trackingVertex);
300 
301  double sr_quality = tracer.quality();
302 
303  if (classifier_.is(VertexCategories::Reconstructed)) {
304  HistIndex_["sr_All_MatchQuality"]->Fill(sr_quality);
305  fillSimToReco("sr_All", tracer.recoVertex(), trackingVertex);
306  sr_nall++;
307 
308  if (classifier_.is(VertexCategories::SecondaryVertex)) {
309  fillSimToReco("sr_SecondaryVertex", tracer.recoVertex(), trackingVertex);
310  sr_nsv++;
311  }
312 
313  if (classifier_.is(VertexCategories::BWeakDecay)) {
314  fillSimToReco("sr_BWeakDecay", tracer.recoVertex(), trackingVertex);
315  sr_nbv++;
316 
317  if (classifier_.is(VertexCategories::SecondaryVertex)) {
318  fillSimToReco("sr_BSV", tracer.recoVertex(), trackingVertex);
319  sr_nbsv++;
320  }
321 
322  } // BWeakDecay
323 
324  else if (classifier_.is(VertexCategories::CWeakDecay)) {
325  fillSimToReco("sr_CWeakDecay", tracer.recoVertex(), trackingVertex);
326  sr_ncv++;
327  }
328 
329  else {
330  fillSimToReco("sr_Light", tracer.recoVertex(), trackingVertex);
331  sr_nlv++;
332  }
333 
334  } // Reconstructed
335  else {
336  // cout << "##### Not reconstructed!" << endl;
337  nmiss++;
338  }
339 
340  } // TVCollection.size()
341 
342  HistIndex_["sr_All_nRecVtx"]->Fill(sr_nall);
343  HistIndex_["sr_SecondaryVertex_nRecVtx"]->Fill(sr_nsv);
344  HistIndex_["sr_BWeakDecay_nRecVtx"]->Fill(sr_nbv);
345  HistIndex_["sr_BSV_nRecVtx"]->Fill(sr_nbsv);
346  HistIndex_["sr_CWeakDecay_nRecVtx"]->Fill(sr_ncv);
347  HistIndex_["sr_Light_nRecVtx"]->Fill(rs_nlv);
348 
349  rs_total_nall += rs_nall;
350  rs_total_nsv += rs_nsv;
351  rs_total_nbv += rs_nbv;
352  rs_total_nbsv += rs_nbsv;
353  rs_total_ncv += rs_ncv;
354  rs_total_nlv += rs_nlv;
355  total_nfake += nfake;
356 
357  sr_total_nall += sr_nall;
358  sr_total_nsv += sr_nsv;
359  sr_total_nbv += sr_nbv;
360  sr_total_nbsv += sr_nbsv;
361  sr_total_ncv += sr_ncv;
362  sr_total_nlv += sr_nlv;
363  total_nmiss += nmiss;
364 }
VertexHistory
This class traces the simulated and generated history of a given track.
Definition: VertexHistory.h:18
mps_fire.i
i
Definition: mps_fire.py:338
VertexClassifier::history
VertexHistory const & history() const
Returns a reference to the vertex history used in the classification.
Definition: VertexClassifier.h:41
VertexHistory::recoVertex
const reco::VertexBaseRef & recoVertex() const
Return a reference to the reconstructed track.
Definition: VertexHistory.h:56
VertexClassifierByProxy::newEvent
void newEvent(edm::Event const &event, edm::EventSetup const &config) override
Pre-process event information (for accessing reconstraction information).
Definition: VertexClassifierByProxy.h:24
recoBSVTagInfoValidationAnalyzer::fillSimToReco
void fillSimToReco(std::string const &, reco::VertexBaseRef const &, TrackingVertexRef const &)
Definition: recoBSVTagInfoValidationAnalyzer.cc:589
VertexCategories::Category
Category
Categories available to vertexes.
Definition: VertexCategories.h:11
recoBSVTagInfoValidationAnalyzer::classifier_
VertexClassifierByProxy< reco::SecondaryVertexTagInfoCollection > classifier_
Definition: recoBSVTagInfoValidationAnalyzer.cc:46
recoBSVTagInfoValidationAnalyzer::tvToken
edm::EDGetTokenT< TrackingVertexCollection > tvToken
Definition: recoBSVTagInfoValidationAnalyzer.cc:85
recoBSVTagInfoValidationAnalyzer::HistIndex_
std::map< std::string, MonitorElement * > HistIndex_
Definition: recoBSVTagInfoValidationAnalyzer.cc:81
VertexCategories::is
bool is(Category category) const
Returns track flag for a given category.
Definition: VertexCategories.h:61
HistoryBase::simVertex
const TrackingVertexRef & simVertex() const
Return the initial tracking vertex from the history.
Definition: HistoryBase.h:70
VertexClassifierByProxy::evaluate
VertexClassifierByProxy< Collection > const & evaluate(TrackingVertexRef const &vertex)
Classify the TrackingVertex in categories.
Definition: VertexClassifierByProxy.h:32
gather_cfg.cout
cout
Definition: gather_cfg.py:144
VertexHistory::quality
double quality() const
Return the quality of the match.
Definition: VertexHistory.h:59
recoBSVTagInfoValidationAnalyzer::fillRecoToSim
void fillRecoToSim(std::string const &, reco::Vertex const &, TrackingVertexRef const &)
Definition: recoBSVTagInfoValidationAnalyzer.cc:496
recoBSVTagInfoValidationAnalyzer::svInfoToken
edm::EDGetTokenT< reco::SecondaryVertexTagInfoCollection > svInfoToken
Definition: recoBSVTagInfoValidationAnalyzer.cc:84
void recoBSVTagInfoValidationAnalyzer::bookRecoToSim ( std::string const &  prefix)
private

Definition at line 366 of file recoBSVTagInfoValidationAnalyzer.cc.

References dqmStore_, HistIndex_, dataset::name, and AlCaHLTBitMon_QueryRunRegistry::string.

Referenced by recoBSVTagInfoValidationAnalyzer().

366  {
367  // Book pull histograms
368 
369  std::string name = prefix + "_Pullx";
370  HistIndex_[name] = dqmStore_->book1D(name.c_str(), name.c_str(), 100, -10., 10.);
371  name = prefix + "_Pully";
372  HistIndex_[name] = dqmStore_->book1D(name.c_str(), name.c_str(), 100, -10., 10.);
373  name = prefix + "_Pullz";
374  HistIndex_[name] = dqmStore_->book1D(name.c_str(), name.c_str(), 100, -10., 10.);
375 
376  name = prefix + "_Resx";
377  HistIndex_[name] = dqmStore_->book1D(name.c_str(), name.c_str(), 100, -0.05, 0.05);
378  name = prefix + "_Resy";
379  HistIndex_[name] = dqmStore_->book1D(name.c_str(), name.c_str(), 100, -0.05, 0.05);
380  name = prefix + "_Resz";
381  HistIndex_[name] = dqmStore_->book1D(name.c_str(), name.c_str(), 100, -0.05, 0.05);
382 
383  name = prefix + "_Chi2Norm";
384  HistIndex_[name] = dqmStore_->book1D(name.c_str(), name.c_str(), 100, 0, 10.);
385  name = prefix + "_Chi2Prob";
386  HistIndex_[name] = dqmStore_->book1D(name.c_str(), name.c_str(), 100, 0., 1.);
387 
388  name = prefix + "_nRecTrks";
389  HistIndex_[name] = dqmStore_->book1D(name.c_str(), name.c_str(), 501, -0.5, 500.5);
390 
391  name = prefix + "_AverageTrackWeight";
392  HistIndex_[name] = dqmStore_->book1D(name.c_str(), name.c_str(), 100, -0.1, 1.1);
393 
394  name = prefix + "_Mass";
395  HistIndex_[name] = dqmStore_->book1D(name.c_str(), name.c_str(), 65, 0., 6.5);
396 
397  name = prefix + "_RecPt";
398  HistIndex_[name] = dqmStore_->book1D(name.c_str(), name.c_str(), 2000, 0., 1000.);
399 
400  name = prefix + "_RecEta";
401  HistIndex_[name] = dqmStore_->book1D(name.c_str(), name.c_str(), 200, -3., 3.);
402 
403  name = prefix + "_RecCharge";
404  HistIndex_[name] = dqmStore_->book1D(name.c_str(), name.c_str(), 21, -0.5, 20.5);
405 
406  name = prefix + "_RecTrackPt";
407  HistIndex_[name] = dqmStore_->book1D(name.c_str(), name.c_str(), 2000, 0., 1000.);
408 
409  name = prefix + "_RecTrackEta";
410  HistIndex_[name] = dqmStore_->book1D(name.c_str(), name.c_str(), 200, -3., 3.);
411 
412  name = prefix + "_nSimTrks";
413  HistIndex_[name] = dqmStore_->book1D(name.c_str(), name.c_str(), 501, -0.5, 500.5);
414 
415  name = prefix + "_SimPt";
416  HistIndex_[name] = dqmStore_->book1D(name.c_str(), name.c_str(), 2000, 0., 1000.);
417 
418  name = prefix + "_SimEta";
419  HistIndex_[name] = dqmStore_->book1D(name.c_str(), name.c_str(), 200, -3., 3.);
420 
421  name = prefix + "_SimCharge";
422  HistIndex_[name] = dqmStore_->book1D(name.c_str(), name.c_str(), 21, -0.5, 20.5);
423 
424  name = prefix + "_SimTrackPt";
425  HistIndex_[name] = dqmStore_->book1D(name.c_str(), name.c_str(), 500, 0., 500.);
426 
427  name = prefix + "_SimTrackEta";
428  HistIndex_[name] = dqmStore_->book1D(name.c_str(), name.c_str(), 200, -3., 3.);
429 }
AlCaHLTBitMon_QueryRunRegistry.string
string
Definition: AlCaHLTBitMon_QueryRunRegistry.py:256
recoBSVTagInfoValidationAnalyzer::HistIndex_
std::map< std::string, MonitorElement * > HistIndex_
Definition: recoBSVTagInfoValidationAnalyzer.cc:81
void recoBSVTagInfoValidationAnalyzer::bookSimToReco ( std::string const &  prefix)
private

Definition at line 431 of file recoBSVTagInfoValidationAnalyzer.cc.

References dqmStore_, HistIndex_, dataset::name, and AlCaHLTBitMon_QueryRunRegistry::string.

Referenced by recoBSVTagInfoValidationAnalyzer().

431  {
432  // Book pull histograms
433 
434  std::string name = prefix + "_Pullx";
435  HistIndex_[name] = dqmStore_->book1D(name.c_str(), name.c_str(), 100, -10., 10.);
436  name = prefix + "_Pully";
437  HistIndex_[name] = dqmStore_->book1D(name.c_str(), name.c_str(), 100, -10., 10.);
438  name = prefix + "_Pullz";
439  HistIndex_[name] = dqmStore_->book1D(name.c_str(), name.c_str(), 100, -10., 10.);
440 
441  name = prefix + "_Resx";
442  HistIndex_[name] = dqmStore_->book1D(name.c_str(), name.c_str(), 100, -0.05, 0.05);
443  name = prefix + "_Resy";
444  HistIndex_[name] = dqmStore_->book1D(name.c_str(), name.c_str(), 100, -0.05, 0.05);
445  name = prefix + "_Resz";
446  HistIndex_[name] = dqmStore_->book1D(name.c_str(), name.c_str(), 100, -0.05, 0.05);
447 
448  name = prefix + "_Chi2Norm";
449  HistIndex_[name] = dqmStore_->book1D(name.c_str(), name.c_str(), 100, 0, 10.);
450  name = prefix + "_Chi2Prob";
451  HistIndex_[name] = dqmStore_->book1D(name.c_str(), name.c_str(), 100, 0., 1.);
452 
453  name = prefix + "_nRecTrks";
454  HistIndex_[name] = dqmStore_->book1D(name.c_str(), name.c_str(), 501, -0.5, 500.5);
455 
456  name = prefix + "_AverageTrackWeight";
457  HistIndex_[name] = dqmStore_->book1D(name.c_str(), name.c_str(), 100, -0.1, 1.1);
458 
459  name = prefix + "_Mass";
460  HistIndex_[name] = dqmStore_->book1D(name.c_str(), name.c_str(), 65, 0., 6.5);
461 
462  name = prefix + "_RecPt";
463  HistIndex_[name] = dqmStore_->book1D(name.c_str(), name.c_str(), 2000, 0., 1000.);
464 
465  name = prefix + "_RecEta";
466  HistIndex_[name] = dqmStore_->book1D(name.c_str(), name.c_str(), 200, -3., 3.);
467 
468  name = prefix + "_RecCharge";
469  HistIndex_[name] = dqmStore_->book1D(name.c_str(), name.c_str(), 21, -0.5, 20.5);
470 
471  name = prefix + "_RecTrackPt";
472  HistIndex_[name] = dqmStore_->book1D(name.c_str(), name.c_str(), 2000, 0., 1000.);
473 
474  name = prefix + "_RecTrackEta";
475  HistIndex_[name] = dqmStore_->book1D(name.c_str(), name.c_str(), 200, -3., 3.);
476 
477  name = prefix + "_nSimTrks";
478  HistIndex_[name] = dqmStore_->book1D(name.c_str(), name.c_str(), 501, -0.5, 500.5);
479 
480  name = prefix + "_SimPt";
481  HistIndex_[name] = dqmStore_->book1D(name.c_str(), name.c_str(), 2000, 0., 1000.);
482 
483  name = prefix + "_SimEta";
484  HistIndex_[name] = dqmStore_->book1D(name.c_str(), name.c_str(), 200, -3., 3.);
485 
486  name = prefix + "_SimCharge";
487  HistIndex_[name] = dqmStore_->book1D(name.c_str(), name.c_str(), 21, -0.5, 20.5);
488 
489  name = prefix + "_SimTrackPt";
490  HistIndex_[name] = dqmStore_->book1D(name.c_str(), name.c_str(), 500, 0., 500.);
491 
492  name = prefix + "_SimTrackEta";
493  HistIndex_[name] = dqmStore_->book1D(name.c_str(), name.c_str(), 200, -3., 3.);
494 }
AlCaHLTBitMon_QueryRunRegistry.string
string
Definition: AlCaHLTBitMon_QueryRunRegistry.py:256
recoBSVTagInfoValidationAnalyzer::HistIndex_
std::map< std::string, MonitorElement * > HistIndex_
Definition: recoBSVTagInfoValidationAnalyzer.cc:81
void recoBSVTagInfoValidationAnalyzer::endJob ( void  )
overrideprivatevirtual

Reimplemented from edm::EDAnalyzer.

Definition at line 682 of file recoBSVTagInfoValidationAnalyzer.cc.

References gather_cfg::cout, DEFINE_FWK_MODULE, n_event, rs_total_nall, rs_total_nbsv, rs_total_nbv, rs_total_ncv, rs_total_nlv, rs_total_nsv, sr_total_nall, sr_total_nbsv, sr_total_nbv, sr_total_ncv, sr_total_nlv, sr_total_nsv, and total_nfake.

682  {
683  std::cout << std::endl;
684  std::cout << " ====== Total Number of analyzed events: " << n_event << " ====== " << std::endl;
685  std::cout << " ====== Total Number of R2S All: " << rs_total_nall << " ====== " << std::endl;
686  std::cout << " ====== Total Number of R2S SecondaryVertex: " << rs_total_nsv << " ====== " << std::endl;
687  std::cout << " ====== Total Number of R2S BWeakDecay: " << rs_total_nbv << " ====== " << std::endl;
688  std::cout << " ====== Total Number of R2S BWeakDecay::SecondaryVertex: " << rs_total_nbsv << " ====== " << std::endl;
689  std::cout << " ====== Total Number of R2S CWeakDecay: " << rs_total_ncv << " ====== " << std::endl;
690  std::cout << " ====== Total Number of R2S Light: " << rs_total_nlv << " ====== " << std::endl;
691  std::cout << std::endl;
692  std::cout << " ====== Total Number of S2R All: " << sr_total_nall << " ====== " << std::endl;
693  std::cout << " ====== Total Number of S2R SecondaryVertex: " << sr_total_nsv << " ====== " << std::endl;
694  std::cout << " ====== Total Number of S2R BWeakDecay: " << sr_total_nbv << " ====== " << std::endl;
695  std::cout << " ====== Total Number of S2R BWeakDecay::SecondaryVertex: " << sr_total_nbsv << " ====== " << std::endl;
696  std::cout << " ====== Total Number of S2R CWeakDecay: " << sr_total_ncv << " ====== " << std::endl;
697  std::cout << " ====== Total Number of S2R Light: " << sr_total_nlv << " ====== " << std::endl;
698  std::cout << std::endl;
699  std::cout << " ====== Total Number of Fake Vertices: " << total_nfake << " ====== " << std::endl;
700 }
gather_cfg.cout
cout
Definition: gather_cfg.py:144
void recoBSVTagInfoValidationAnalyzer::fillRecoToSim ( std::string const &  prefix,
reco::Vertex const &  vertex,
TrackingVertexRef const &  simVertex 
)
private

Definition at line 496 of file recoBSVTagInfoValidationAnalyzer.cc.

References ALCARECOTkAlJpsiMuMu_cff::charge, reco::Vertex::chi2(), ChiSquaredProbability(), PVValHelper::eta, HistIndex_, ResonanceBuilder::mass, reco::Vertex::ndof(), reco::Vertex::normalizedChi2(), AlCaHLTBitMon_ParallelJobs::p, EnergyCorrector::pt, PVValHelper::resz, mathSSE::sqrt(), reco::Vertex::tracks_begin(), reco::Vertex::tracks_end(), reco::Vertex::tracksSize(), reco::Vertex::trackWeight(), mps_merge::weight, reco::Vertex::x(), reco::Vertex::xError(), reco::Vertex::y(), reco::Vertex::yError(), reco::Vertex::z(), and reco::Vertex::zError().

Referenced by analyze().

498  {
499  double pullx = (vertex.x() - simVertex->position().x()) / vertex.xError();
500  double pully = (vertex.y() - simVertex->position().y()) / vertex.yError();
501  double pullz = (vertex.z() - simVertex->position().z()) / vertex.zError();
502 
503  double resx = vertex.x() - simVertex->position().x();
504  double resy = vertex.y() - simVertex->position().y();
505  double resz = vertex.z() - simVertex->position().z();
506 
507  double chi2norm = vertex.normalizedChi2();
508  double chi2prob = ChiSquaredProbability(vertex.chi2(), vertex.ndof());
509 
510  double sum_weight = 0.;
511  double weight = 0.;
512  double tracksize = vertex.tracksSize();
513  math::XYZVector momentum;
514  math::XYZTLorentzVector sum;
515  int charge = 0;
516  double thePiMass = 0.13957;
517  for (reco::Vertex::trackRef_iterator recDaughter = vertex.tracks_begin(); recDaughter != vertex.tracks_end();
518  ++recDaughter) {
519  ROOT::Math::LorentzVector<ROOT::Math::PxPyPzM4D<double>> vec;
520 
521  vec.SetPx((**recDaughter).px());
522  vec.SetPy((**recDaughter).py());
523  vec.SetPz((**recDaughter).pz());
524  vec.SetM(thePiMass);
525 
526  sum += vec;
527 
528  weight = vertex.trackWeight(*recDaughter);
529  sum_weight += weight;
530 
531  math::XYZVector p;
532  p = (**recDaughter).momentum();
533  momentum += p;
534 
535  charge += (*recDaughter)->charge();
536 
537  HistIndex_[prefix + "_RecTrackPt"]->Fill((*recDaughter)->pt());
538  HistIndex_[prefix + "_RecTrackEta"]->Fill((*recDaughter)->eta());
539  } // end loop to recDaughters
540  // cout << " average sum of weights = " <<
541  // sum_weight/tracksize << endl;
542 
543  double mass = sum.M();
544  double pt = sqrt(momentum.Perp2());
545  double eta = momentum.Eta();
546 
547  math::XYZVector simmomentum;
548  int simcharge = 0;
549  for (TrackingVertex::tp_iterator simDaughter = simVertex->daughterTracks_begin();
550  simDaughter != simVertex->daughterTracks_end();
551  ++simDaughter) {
552  math::XYZVector p;
553  p = (**simDaughter).momentum();
554  simmomentum += p;
555 
556  simcharge += (*simDaughter)->charge();
557 
558  HistIndex_[prefix + "_SimTrackPt"]->Fill((*simDaughter)->pt());
559  HistIndex_[prefix + "_SimTrackEta"]->Fill((*simDaughter)->eta());
560  }
561 
562  double simpt = sqrt(simmomentum.Perp2());
563  double simeta = simmomentum.Eta();
564 
565  // cout << "[fillRecoToSim] vertex.tracksSize() = " << vertex.tracksSize() <<
566  // " ; simVertex->nDaughterTracks() = " << simVertex->nDaughterTracks() <<
567  // endl;
568 
569  HistIndex_[prefix + "_nRecTrks"]->Fill(vertex.tracksSize());
570  HistIndex_[prefix + "_nSimTrks"]->Fill(simVertex->nDaughterTracks());
571  HistIndex_[prefix + "_Pullx"]->Fill(pullx);
572  HistIndex_[prefix + "_Pully"]->Fill(pully);
573  HistIndex_[prefix + "_Pullz"]->Fill(pullz);
574  HistIndex_[prefix + "_Resx"]->Fill(resx);
575  HistIndex_[prefix + "_Resy"]->Fill(resy);
576  HistIndex_[prefix + "_Resz"]->Fill(resz);
577  HistIndex_[prefix + "_AverageTrackWeight"]->Fill(sum_weight / tracksize);
578  HistIndex_[prefix + "_Chi2Norm"]->Fill(chi2norm);
579  HistIndex_[prefix + "_Chi2Prob"]->Fill(chi2prob);
580  HistIndex_[prefix + "_RecPt"]->Fill(pt);
581  HistIndex_[prefix + "_RecEta"]->Fill(eta);
582  HistIndex_[prefix + "_RecCharge"]->Fill(charge);
583  HistIndex_[prefix + "_Mass"]->Fill(mass);
584  HistIndex_[prefix + "_SimPt"]->Fill(simpt);
585  HistIndex_[prefix + "_SimEta"]->Fill(simeta);
586  HistIndex_[prefix + "_SimCharge"]->Fill(simcharge);
587 }
AlCaHLTBitMon_ParallelJobs.p
p
Definition: AlCaHLTBitMon_ParallelJobs.py:153
ecalDetailedTimeRecHit_cfi.simVertex
simVertex
Definition: ecalDetailedTimeRecHit_cfi.py:17
weight
Definition: weight.py:1
math::XYZTLorentzVector
XYZTLorentzVectorD XYZTLorentzVector
Lorentz vector with cylindrical internal representation using pseudorapidity.
Definition: LorentzVector.h:29
mathSSE::sqrt
T sqrt(T t)
Definition: SSEVec.h:18
recoBSVTagInfoValidationAnalyzer::HistIndex_
std::map< std::string, MonitorElement * > HistIndex_
Definition: recoBSVTagInfoValidationAnalyzer.cc:81
ChiSquaredProbability
float ChiSquaredProbability(double chiSquared, double nrDOF)
Definition: ChiSquaredProbability.cc:13
ALCARECOTkAlJpsiMuMu_cff.charge
charge
Definition: ALCARECOTkAlJpsiMuMu_cff.py:47
math::XYZVector
XYZVectorD XYZVector
spatial vector with cartesian internal representation
Definition: Vector3D.h:30
reco::Vertex::trackRef_iterator
std::vector< TrackBaseRef >::const_iterator trackRef_iterator
The iteratator for the vector<TrackRef>
Definition: Vertex.h:37
mps_merge.weight
weight
Definition: mps_merge.py:88
void recoBSVTagInfoValidationAnalyzer::fillSimToReco ( std::string const &  prefix,
reco::VertexBaseRef const &  vertex,
TrackingVertexRef const &  simVertex 
)
private

Definition at line 589 of file recoBSVTagInfoValidationAnalyzer.cc.

References ALCARECOTkAlJpsiMuMu_cff::charge, reco::Vertex::chi2(), ChiSquaredProbability(), PVValHelper::eta, HistIndex_, ResonanceBuilder::mass, reco::Vertex::ndof(), reco::Vertex::normalizedChi2(), AlCaHLTBitMon_ParallelJobs::p, EnergyCorrector::pt, PVValHelper::resz, mathSSE::sqrt(), reco::Vertex::tracks_begin(), reco::Vertex::tracks_end(), reco::Vertex::tracksSize(), reco::Vertex::trackWeight(), mps_merge::weight, reco::Vertex::x(), reco::Vertex::xError(), reco::Vertex::y(), reco::Vertex::yError(), reco::Vertex::z(), and reco::Vertex::zError().

Referenced by analyze().

591  {
592  double pullx = (vertex->x() - simVertex->position().x()) / vertex->xError();
593  double pully = (vertex->y() - simVertex->position().y()) / vertex->yError();
594  double pullz = (vertex->z() - simVertex->position().z()) / vertex->zError();
595 
596  double resx = vertex->x() - simVertex->position().x();
597  double resy = vertex->y() - simVertex->position().y();
598  double resz = vertex->z() - simVertex->position().z();
599 
600  double chi2norm = vertex->normalizedChi2();
601  double chi2prob = ChiSquaredProbability(vertex->chi2(), vertex->ndof());
602 
603  double sum_weight = 0.;
604  double weight = 0.;
605  double tracksize = vertex->tracksSize();
606  math::XYZVector momentum;
607  math::XYZTLorentzVector sum;
608  int charge = 0;
609  double thePiMass = 0.13957;
610  for (reco::Vertex::trackRef_iterator recDaughter = vertex->tracks_begin(); recDaughter != vertex->tracks_end();
611  ++recDaughter) {
612  ROOT::Math::LorentzVector<ROOT::Math::PxPyPzM4D<double>> vec;
613 
614  vec.SetPx((**recDaughter).px());
615  vec.SetPy((**recDaughter).py());
616  vec.SetPz((**recDaughter).pz());
617  vec.SetM(thePiMass);
618 
619  sum += vec;
620 
621  weight = vertex->trackWeight(*recDaughter);
622  sum_weight += weight;
623 
624  math::XYZVector p;
625  p = (**recDaughter).momentum();
626  momentum += p;
627 
628  charge += (*recDaughter)->charge();
629 
630  HistIndex_[prefix + "_RecTrackPt"]->Fill((*recDaughter)->pt());
631  HistIndex_[prefix + "_RecTrackEta"]->Fill((*recDaughter)->eta());
632  }
633  // cout << " average sum of weights = " <<
634  // sum_weight/tracksize << endl;
635 
636  double mass = sum.M();
637  double pt = sqrt(momentum.Perp2());
638  double eta = momentum.Eta();
639 
640  math::XYZVector simmomentum;
641  int simcharge = 0;
642  for (TrackingVertex::tp_iterator simDaughter = simVertex->daughterTracks_begin();
643  simDaughter != simVertex->daughterTracks_end();
644  ++simDaughter) {
645  math::XYZVector p;
646  p = (**simDaughter).momentum();
647  simmomentum += p;
648 
649  simcharge += (*simDaughter)->charge();
650 
651  HistIndex_[prefix + "_SimTrackPt"]->Fill((*simDaughter)->pt());
652  HistIndex_[prefix + "_SimTrackEta"]->Fill((*simDaughter)->eta());
653  }
654 
655  double simpt = sqrt(simmomentum.Perp2());
656  double simeta = simmomentum.Eta();
657 
658  // cout << "[fillSimToReco] vertex->tracksSize() = " << vertex->tracksSize()
659  // << " ; simVertex->nDaughterTracks() = " << simVertex->nDaughterTracks() <<
660  // endl;
661 
662  HistIndex_[prefix + "_nRecTrks"]->Fill(vertex->tracksSize());
663  HistIndex_[prefix + "_nSimTrks"]->Fill(simVertex->nDaughterTracks());
664  HistIndex_[prefix + "_Pullx"]->Fill(pullx);
665  HistIndex_[prefix + "_Pully"]->Fill(pully);
666  HistIndex_[prefix + "_Pullz"]->Fill(pullz);
667  HistIndex_[prefix + "_Resx"]->Fill(resx);
668  HistIndex_[prefix + "_Resy"]->Fill(resy);
669  HistIndex_[prefix + "_Resz"]->Fill(resz);
670  HistIndex_[prefix + "_AverageTrackWeight"]->Fill(sum_weight / tracksize);
671  HistIndex_[prefix + "_Chi2Norm"]->Fill(chi2norm);
672  HistIndex_[prefix + "_Chi2Prob"]->Fill(chi2prob);
673  HistIndex_[prefix + "_RecPt"]->Fill(pt);
674  HistIndex_[prefix + "_RecEta"]->Fill(eta);
675  HistIndex_[prefix + "_RecCharge"]->Fill(charge);
676  HistIndex_[prefix + "_Mass"]->Fill(mass);
677  HistIndex_[prefix + "_SimPt"]->Fill(simpt);
678  HistIndex_[prefix + "_SimEta"]->Fill(simeta);
679  HistIndex_[prefix + "_SimCharge"]->Fill(simcharge);
680 }
AlCaHLTBitMon_ParallelJobs.p
p
Definition: AlCaHLTBitMon_ParallelJobs.py:153
ecalDetailedTimeRecHit_cfi.simVertex
simVertex
Definition: ecalDetailedTimeRecHit_cfi.py:17
weight
Definition: weight.py:1
math::XYZTLorentzVector
XYZTLorentzVectorD XYZTLorentzVector
Lorentz vector with cylindrical internal representation using pseudorapidity.
Definition: LorentzVector.h:29
mathSSE::sqrt
T sqrt(T t)
Definition: SSEVec.h:18
recoBSVTagInfoValidationAnalyzer::HistIndex_
std::map< std::string, MonitorElement * > HistIndex_
Definition: recoBSVTagInfoValidationAnalyzer.cc:81
ChiSquaredProbability
float ChiSquaredProbability(double chiSquared, double nrDOF)
Definition: ChiSquaredProbability.cc:13
ALCARECOTkAlJpsiMuMu_cff.charge
charge
Definition: ALCARECOTkAlJpsiMuMu_cff.py:47
math::XYZVector
XYZVectorD XYZVector
spatial vector with cartesian internal representation
Definition: Vector3D.h:30
reco::Vertex::trackRef_iterator
std::vector< TrackBaseRef >::const_iterator trackRef_iterator
The iteratator for the vector<TrackRef>
Definition: Vertex.h:37
mps_merge.weight
weight
Definition: mps_merge.py:88

Member Data Documentation

VertexClassifierByProxy<reco::SecondaryVertexTagInfoCollection> recoBSVTagInfoValidationAnalyzer::classifier_
private

Definition at line 46 of file recoBSVTagInfoValidationAnalyzer.cc.

Referenced by analyze().

std::string recoBSVTagInfoValidationAnalyzer::dqmLabel
private

Definition at line 53 of file recoBSVTagInfoValidationAnalyzer.cc.

Referenced by recoBSVTagInfoValidationAnalyzer().

DQMStore* recoBSVTagInfoValidationAnalyzer::dqmStore_
private

Definition at line 52 of file recoBSVTagInfoValidationAnalyzer.cc.

Referenced by bookRecoToSim(), bookSimToReco(), and recoBSVTagInfoValidationAnalyzer().

std::map<std::string, MonitorElement *> recoBSVTagInfoValidationAnalyzer::HistIndex_
private

Definition at line 81 of file recoBSVTagInfoValidationAnalyzer.cc.

Referenced by analyze(), bookRecoToSim(), bookSimToReco(), fillRecoToSim(), fillSimToReco(), and recoBSVTagInfoValidationAnalyzer().

Int_t recoBSVTagInfoValidationAnalyzer::n_event
private

Definition at line 55 of file recoBSVTagInfoValidationAnalyzer.cc.

Referenced by analyze(), endJob(), and recoBSVTagInfoValidationAnalyzer().

Int_t recoBSVTagInfoValidationAnalyzer::numberVertexClassifier_
private

Definition at line 48 of file recoBSVTagInfoValidationAnalyzer.cc.

Referenced by analyze(), and recoBSVTagInfoValidationAnalyzer().

Int_t recoBSVTagInfoValidationAnalyzer::rs_total_nall
private

Definition at line 56 of file recoBSVTagInfoValidationAnalyzer.cc.

Referenced by analyze(), endJob(), and recoBSVTagInfoValidationAnalyzer().

Int_t recoBSVTagInfoValidationAnalyzer::rs_total_nbsv
private

Definition at line 59 of file recoBSVTagInfoValidationAnalyzer.cc.

Referenced by analyze(), endJob(), and recoBSVTagInfoValidationAnalyzer().

Int_t recoBSVTagInfoValidationAnalyzer::rs_total_nbv
private

Definition at line 58 of file recoBSVTagInfoValidationAnalyzer.cc.

Referenced by analyze(), endJob(), and recoBSVTagInfoValidationAnalyzer().

Int_t recoBSVTagInfoValidationAnalyzer::rs_total_ncv
private

Definition at line 60 of file recoBSVTagInfoValidationAnalyzer.cc.

Referenced by analyze(), endJob(), and recoBSVTagInfoValidationAnalyzer().

Int_t recoBSVTagInfoValidationAnalyzer::rs_total_nlv
private

Definition at line 61 of file recoBSVTagInfoValidationAnalyzer.cc.

Referenced by analyze(), endJob(), and recoBSVTagInfoValidationAnalyzer().

Int_t recoBSVTagInfoValidationAnalyzer::rs_total_nsv
private

Definition at line 57 of file recoBSVTagInfoValidationAnalyzer.cc.

Referenced by analyze(), endJob(), and recoBSVTagInfoValidationAnalyzer().

Int_t recoBSVTagInfoValidationAnalyzer::sr_total_nall
private

Definition at line 64 of file recoBSVTagInfoValidationAnalyzer.cc.

Referenced by analyze(), endJob(), and recoBSVTagInfoValidationAnalyzer().

Int_t recoBSVTagInfoValidationAnalyzer::sr_total_nbsv
private

Definition at line 67 of file recoBSVTagInfoValidationAnalyzer.cc.

Referenced by analyze(), endJob(), and recoBSVTagInfoValidationAnalyzer().

Int_t recoBSVTagInfoValidationAnalyzer::sr_total_nbv
private

Definition at line 66 of file recoBSVTagInfoValidationAnalyzer.cc.

Referenced by analyze(), endJob(), and recoBSVTagInfoValidationAnalyzer().

Int_t recoBSVTagInfoValidationAnalyzer::sr_total_ncv
private

Definition at line 68 of file recoBSVTagInfoValidationAnalyzer.cc.

Referenced by analyze(), endJob(), and recoBSVTagInfoValidationAnalyzer().

Int_t recoBSVTagInfoValidationAnalyzer::sr_total_nlv
private

Definition at line 69 of file recoBSVTagInfoValidationAnalyzer.cc.

Referenced by analyze(), endJob(), and recoBSVTagInfoValidationAnalyzer().

Int_t recoBSVTagInfoValidationAnalyzer::sr_total_nsv
private

Definition at line 65 of file recoBSVTagInfoValidationAnalyzer.cc.

Referenced by analyze(), endJob(), and recoBSVTagInfoValidationAnalyzer().

edm::EDGetTokenT<reco::SecondaryVertexTagInfoCollection> recoBSVTagInfoValidationAnalyzer::svInfoToken
private

Definition at line 84 of file recoBSVTagInfoValidationAnalyzer.cc.

Referenced by analyze(), and recoBSVTagInfoValidationAnalyzer().

edm::InputTag recoBSVTagInfoValidationAnalyzer::svTagInfoProducer_
private

Definition at line 50 of file recoBSVTagInfoValidationAnalyzer.cc.

Int_t recoBSVTagInfoValidationAnalyzer::total_nfake
private

Definition at line 62 of file recoBSVTagInfoValidationAnalyzer.cc.

Referenced by analyze(), endJob(), and recoBSVTagInfoValidationAnalyzer().

Int_t recoBSVTagInfoValidationAnalyzer::total_nmiss
private

Definition at line 70 of file recoBSVTagInfoValidationAnalyzer.cc.

Referenced by analyze(), and recoBSVTagInfoValidationAnalyzer().

edm::InputTag recoBSVTagInfoValidationAnalyzer::trackingTruth_
private

Definition at line 49 of file recoBSVTagInfoValidationAnalyzer.cc.

edm::EDGetTokenT<TrackingVertexCollection> recoBSVTagInfoValidationAnalyzer::tvToken
private

Definition at line 85 of file recoBSVTagInfoValidationAnalyzer.cc.

Referenced by analyze(), and recoBSVTagInfoValidationAnalyzer().

Generated for CMSSW Reference Manual by doxygen 1.8.11