MtdValid: more plots for track multiplicity

Author: kdeleo

Validation/MtdValidation/plugins/Primary4DVertexValidation.cc

@@ -334,7 +334,9 @@ class Primary4DVertexValidation : public DQMEDAnalyzer {
 
   // histogram declaration
   MonitorElement* meUnAssocTracks_;
+  MonitorElement* meUnAssocTracksFake_;
   MonitorElement* meFractionUnAssocTracks_;
+  MonitorElement* meFractionUnAssocTracksFake_;
   MonitorElement* meTrackEffPtTot_;
   MonitorElement* meTrackMatchedTPEffPtTot_;
   MonitorElement* meTrackMatchedTPEffPtMtd_;
@@ -383,6 +385,7 @@ class Primary4DVertexValidation : public DQMEDAnalyzer {
   MonitorElement* meSimPVTvsZ_;
 
   MonitorElement* meVtxTrackMult_;
+  MonitorElement* meVtxTrackMultAllVtx_;
   MonitorElement* meVtxTrackW_;
   MonitorElement* meVtxTrackWnt_;
   MonitorElement* meVtxTrackRecLVMult_;
@@ -590,7 +593,10 @@ void Primary4DVertexValidation::bookHistograms(DQMStore::IBooker& ibook,
   ibook.setCurrentFolder(folder_);
   // --- histograms booking
   meUnAssocTracks_ = ibook.book1D("UnAssocTracks", "Unassociated tracks", 160, 0.5, 4.5);
+  meUnAssocTracksFake_ = ibook.book1D("UnAssocTracksFake", "Unassociated fake tracks", 160, 0.5, 4.5);
   meFractionUnAssocTracks_ = ibook.book1D("FractionUnAssocTracks", "Fraction Unassociated tracks", 160, 0.0, 1.);
+  meFractionUnAssocTracksFake_ =
+      ibook.book1D("FractionUnAssocTracksFake", "Fraction Unassociated fake tracks", 160, 0.0, 1.);
   meTrackEffPtTot_ = ibook.book1D("EffPtTot", "Pt of tracks associated to LV; track pt [GeV] ", 110, 0., 11.);
   meTrackEffEtaTot_ = ibook.book1D("EffEtaTot", "Eta of tracks associated to LV; track eta ", 66, 0., 3.3);
   meTrackMatchedTPEffPtTot_ =
@@ -727,6 +733,7 @@ void Primary4DVertexValidation::bookHistograms(DQMStore::IBooker& ibook,
   meSimPVTvsZ_ = ibook.bookProfile("simPVTvsZ", "PV Time vs Z", 30, -15., 15., 30, -0.75, 0.75);
 
   meVtxTrackMult_ = ibook.book1D("VtxTrackMult", "Log10(Vertex track multiplicity)", 80, 0.5, 2.5);
+  meVtxTrackMultAllVtx_ = ibook.book1D("VtxTrackMultAllVtx", "Log10(Vertex track multiplicity all vtx)", 80, 0.5, 2.5);
   meVtxTrackW_ = ibook.book1D("VtxTrackW", "Vertex track weight (all)", 50, 0., 1.);
   meVtxTrackWnt_ = ibook.book1D("VtxTrackWnt", "Vertex track Wnt", 50, 0., 1.);
   meVtxTrackRecLVMult_ =
@@ -2095,6 +2102,7 @@ void Primary4DVertexValidation::analyze(const edm::Event& iEvent, const edm::Eve
   std::vector<reco::TransientTrack>&& seltks = theTrackFilter->select(t_tks);
 
   int unassociatedCount = 0;
+  int unassociatedCountFake = 0;
   for (std::vector<reco::TransientTrack>::const_iterator itk = seltks.begin(); itk != seltks.end(); itk++) {
     reco::TrackBaseRef trackref = (*itk).trackBaseRef();
     bool isAssociated = false;
@@ -2112,12 +2120,19 @@ void Primary4DVertexValidation::analyze(const edm::Event& iEvent, const edm::Eve
 
     if (!isAssociated) {
       unassociatedCount++;
+      auto found = r2s_->find(trackref);
+      if (found == r2s_->end())
+        unassociatedCountFake++;
     }
   }
   double fraction = double(unassociatedCount) / (seltks.size());
   meUnAssocTracks_->Fill(log10(unassociatedCount));
   meFractionUnAssocTracks_->Fill(fraction);
 
+  double fractionFake = double(unassociatedCountFake) / (seltks.size());
+  meUnAssocTracksFake_->Fill(log10(unassociatedCountFake));
+  meFractionUnAssocTracksFake_->Fill(fractionFake);
+
   // Loop on tracks
   for (unsigned int iv = 0; iv < recopv.size(); iv++) {
     if (recopv.at(iv).ndof > selNdof_) {
@@ -2662,6 +2677,17 @@ void Primary4DVertexValidation::analyze(const edm::Event& iEvent, const edm::Eve
   for (unsigned int ir = 0; ir < recopv.size(); ir++) {
     if (recopv.at(ir).ndof > selNdof_) {
       meRecPVZ_->Fill(recopv.at(ir).z);
+      const reco::Vertex* vertex = recopv.at(ir).recVtx;
+      unsigned int nt = 0;
+      for (auto iTrack = vertex->tracks_begin(); iTrack != vertex->tracks_end(); ++iTrack) {
+        if (trackAssoc[*iTrack] == -1) {
+          edm::LogWarning("mtdTracks") << "Extended track not associated";
+          continue;
+        }
+        nt++;
+      }
+      meVtxTrackMultAllVtx_->Fill(log10(nt));
+
       if (recopv.at(ir).recVtx->tError() > 0.) {
         meRecPVT_->Fill(recopv.at(ir).recVtx->t());
       }