WIP: fix MCLinking in TDR branch for LHCb::Tracks

@sponce , @raaij , @sstahl The solution I adopted to restore MCLinking for the LHCb::Tracks -> vector conversion, is to append in the MCSequence an algorithm reading all the tracks in the Track container and make a copy of them in /Copy/Track. This allow to run the PrChecker without any issue. Maybe this can be extended to other algorithm which do not allow to merge TDR to master. @gligorov , you may need this hack to setup a timing & efficiency test for CPU using the exact same software stack. To - do -list

• Make the TrackConverter (move to functional) a 1-to-1 conversion One path In , One path Out (KeyedTrack) chained by the PrTrackAssociator using single containers. The final MCChecking sequence will be Per track type produced, i.e. MCCHeckingVelo, MCCheckingUpstream, MCCheckingForward, where optionally, whether no KeyedContainer is produced, a converter is placed in the sequence.
• Let TrackSys and UpgradeChecking sequence to know where converted Keyed tracks are located and which tracks have to be converted (this depends on the signature of the HLT algorithms, i.e. if they produce LHCb::Tracks or vectorLHCb::Track

Now, you can pass TrackSys().TracksToConvert = ["Velo","Upstream","ForwardFast"] to convert from vector<Track> to LHCb::Tracks when doing MC association. Also TrackSys().VeloUpgradeOnly = True will execute VeloOnly ( and Decoding detectors) in the sequence. The last can be useful concerning @cattanem suggestions in Brunel!344 (closed)

Pr/PrVeloUT/src/PrVeloUT.cpp

@@ -32,6 +32,7 @@ StatusCode PrVeloUT::initialize() {
     m_veloUTTime = m_timerTool->addTimer( "Internal VeloUT Tracking" );
     m_timerTool->decreaseIndent();
   }
+  counter("#tracks");
   return StatusCode::SUCCESS;
 }
 
@@ -60,17 +61,20 @@ std::vector<LHCb::Track> PrVeloUT::operator()(const std::vector<LHCb::Track>& in
   boost::optional<LHCb::Track> outVeloUT;
   counter("#seeds") += inputTracks.size();
 
-  // list tracks to be ignored as they are not respecting the Impact Parameter cut
-  std::map<const LHCb::Track*, bool> cutTracks;
-  for (auto& vertex : vertices) {
-    for (const LHCb::Track* track : vertex.tracks()) {
-      cutTracks[track] = acceptTrackForIPCut(*track, vertex, m_minIP);
+  auto eventState = m_veloUTTool->createState();
+  bool killtrack = false;
+  for(const auto& veloTr: inputTracks) {
+    //you  want to process only those velo input track having IP > m_minIP w.r.t. to all PV,
+    //i.e., if it has for at least 1 PV IP < m_minIP, the track is not propagated further
+
+    killtrack = false;
+    if( m_doIPCut){
+      for( const auto& vertex : vertices){
+	if( !acceptTrackForIPCut( veloTr, vertex, m_minIP)  ){ killtrack = true; break;}
+      }
     }
-  }
+    if( m_doIPCut && killtrack) continue;
 
-  auto eventState = m_veloUTTool->createState();
-  for(const auto* veloTr: inputTracks) {
-    if (!cutTracks[veloTr]) continue;
     m_veloUTTool->tracksFromTrack(veloTr, outVeloUT, eventState).ignore();
     if(outVeloUT){
       outputTracks.push_back(std::move(outVeloUT.value()));