Log statistics on track to vertex assignment

shahor02 · shahor02 · commit ae53b9881ded · 2021-05-19T08:04:45.000+01:00
diff --git a/Detectors/Vertexing/src/VertexTrackMatcher.cxx b/Detectors/Vertexing/src/VertexTrackMatcher.cxx
@@ -95,7 +95,7 @@ void VertexTrackMatcher::process(const o2::globaltracking::RecoContainer& recoDa
 
   extractTracks(recoData, vcont); // extract all track t-brackets, excluding those tracks which contribute to vertex (already attached)
 
-  int ivStart = 0;
+  int ivStart = 0, nAssigned = 0, nAmbiguous = 0;
   std::vector<int> vtxList; // list of vertices which match to checked track
   for (const auto& tro : mTBrackets) {
     vtxList.clear();
@@ -116,13 +116,15 @@ void VertexTrackMatcher::process(const o2::globaltracking::RecoContainer& recoDa
       // track matches to vertex, register
       vtxList.push_back(vto.origID); // flag matching vertex
     }
-    bool ambigous = vtxList.size() > 1; // did track match to multiple vertices?
-    for (auto v : vtxList) {
-      auto& ref = tmpMap[v].emplace_back(tro.origID);
-      if (ambigous) {
-        ref.setAmbiguous();
+    if (vtxList.size() > 1) { // did track match to multiple vertices?
+      nAmbiguous++;
+      for (auto v : vtxList) {
+        tmpMap[v].emplace_back(tro.origID).setAmbiguous();
       }
     }
+    if (!vtxList.empty()) {
+      nAssigned++;
+    }
   }
 
   // build final vector of global indices
@@ -151,6 +153,7 @@ void VertexTrackMatcher::process(const o2::globaltracking::RecoContainer& recoDa
     vr.setEnd(trackIndex.size());
     LOG(INFO) << "Vertxex " << iv << " Tracks " << vr;
   }
+  LOG(INFO) << "Assigned " << nAssigned << " (" << nAmbiguous << " ambigously) out of " << mTBrackets.size() << " non-contributor tracks + " << vcont.size() << " contributors";
 }
 
 //________________________________________________________
@@ -185,5 +188,5 @@ void VertexTrackMatcher::extractTracks(const o2::globaltracking::RecoContainer&
   // sort in increasing min.time
   std::sort(mTBrackets.begin(), mTBrackets.end(), [](const TrackTBracket& a, const TrackTBracket& b) { return a.tBracket.getMin() < b.tBracket.getMin(); });
 
-  LOG(INFO) << "collected " << mTBrackets.size() << " seeds";
+  LOG(INFO) << "collected " << mTBrackets.size() << " non-contributor and " << vcont.size() << " contributor seeds";
 }

Original file line number	Diff line number	Diff line change
`@@ -95,7 +95,7 @@ void VertexTrackMatcher::process(const o2::globaltracking::RecoContainer& recoDa`
`95`	`95`
`96`	`96`	`extractTracks(recoData, vcont); // extract all track t-brackets, excluding those tracks which contribute to vertex (already attached)`
`97`	`97`
`98`		`- int ivStart = 0;`
	`98`	`+ int ivStart = 0, nAssigned = 0, nAmbiguous = 0;`
`99`	`99`	`std::vector<int> vtxList; // list of vertices which match to checked track`
`100`	`100`	`for (const auto& tro : mTBrackets) {`
`101`	`101`	`vtxList.clear();`
`@@ -116,13 +116,15 @@ void VertexTrackMatcher::process(const o2::globaltracking::RecoContainer& recoDa`
`116`	`116`	`// track matches to vertex, register`
`117`	`117`	`vtxList.push_back(vto.origID); // flag matching vertex`
`118`	`118`	`}`
`119`		`- bool ambigous = vtxList.size() > 1; // did track match to multiple vertices?`
`120`		`- for (auto v : vtxList) {`
`121`		`- auto& ref = tmpMap[v].emplace_back(tro.origID);`
`122`		`- if (ambigous) {`
`123`		`- ref.setAmbiguous();`
	`119`	`+ if (vtxList.size() > 1) { // did track match to multiple vertices?`
	`120`	`+ nAmbiguous++;`
	`121`	`+ for (auto v : vtxList) {`
	`122`	`+ tmpMap[v].emplace_back(tro.origID).setAmbiguous();`
`124`	`123`	`}`
`125`	`124`	`}`
	`125`	`+ if (!vtxList.empty()) {`
	`126`	`+ nAssigned++;`
	`127`	`+ }`
`126`	`128`	`}`
`127`	`129`
`128`	`130`	`// build final vector of global indices`
`@@ -151,6 +153,7 @@ void VertexTrackMatcher::process(const o2::globaltracking::RecoContainer& recoDa`
`151`	`153`	`vr.setEnd(trackIndex.size());`
`152`	`154`	`LOG(INFO) << "Vertxex " << iv << " Tracks " << vr;`
`153`	`155`	`}`
	`156`	`+ LOG(INFO) << "Assigned " << nAssigned << " (" << nAmbiguous << " ambigously) out of " << mTBrackets.size() << " non-contributor tracks + " << vcont.size() << " contributors";`
`154`	`157`	`}`
`155`	`158`
`156`	`159`	`//________________________________________________________`
`@@ -185,5 +188,5 @@ void VertexTrackMatcher::extractTracks(const o2::globaltracking::RecoContainer&`
`185`	`188`	`// sort in increasing min.time`
`186`	`189`	`std::sort(mTBrackets.begin(), mTBrackets.end(), [](const TrackTBracket& a, const TrackTBracket& b) { return a.tBracket.getMin() < b.tBracket.getMin(); });`
`187`	`190`
`188`		`- LOG(INFO) << "collected " << mTBrackets.size() << " seeds";`
	`191`	`+ LOG(INFO) << "collected " << mTBrackets.size() << " non-contributor and " << vcont.size() << " contributor seeds";`
`189`	`192`	`}`