Removed histogram filling for reassigned MFT tracks rather fill them in tables

Removed histogram filling for reassigned and assigned MFT track ambiguity degrees. Added checks for Cellular Automaton and LTF track-finding algorithm requirements.

Author: abmodak

PWGCF/TwoParticleCorrelations/TableProducer/longrangeMaker.cxx

@@ -142,6 +142,7 @@ struct LongrangeMaker {
   struct : ConfigurableGroup {
     Configurable<bool> cfgUseChi2Cut{"cfgUseChi2Cut", false, "Use condition on MFT track: chi2/Nclusters"};
     Configurable<bool> cfgRequireCA{"cfgRequireCA", false, "Use Cellular Automaton track-finding algorithm"};
+    Configurable<bool> cfgRequireLTF{"cfgRequireLTF", false, "Use LTF track-finding algorithm"};
     Configurable<bool> useMftPtCut{"useMftPtCut", true, "Choose to apply MFT track pT cut"};
     Configurable<int> cfgMftCluster{"cfgMftCluster", 5, "cut on MFT Cluster"};
     Configurable<float> cfgMftEtaMax{"cfgMftEtaMax", -2.4f, "Maximum MFT eta cut"};
@@ -277,9 +278,6 @@ struct LongrangeMaker {
     xMftBestTrk->SetBinLabel(4, "DCAxy selection");
     xMftBestTrk->SetBinLabel(5, "DCAz selection");
 
-    histos.add("ReassignedMFTtrackAmbDegree", "ReassignedMFTtrackAmbDegree", kTH1D, {cfgAxis.axisMFTAmbDegree});
-    histos.add("AssignedMFTtrackAmbDegree", "AssignedMFTtrackAmbDegree", kTH1D, {cfgAxis.axisMFTAmbDegree});
-
     histos.add("FT0A_Amp", "FT0A_Amp", kTH1D, {cfgAxis.axisAmplitude});
     histos.add("FT0A_Amp_gaincorrected", "FT0A_Amp_gaincorrected", kTH1D, {cfgAxis.axisAmplitude});
     histos.add("FT0A_Channel_vs_Amp", "FT0A_Channel_vs_Amp", kTH2D, {cfgAxis.axisChannel, cfgAxis.axisAmplitude});
@@ -483,12 +481,6 @@ struct LongrangeMaker {
         continue;
       }
       auto phi = itrack.phi();
-      if (itrack.collisionId() != reassoMftTrack.bestCollisionId()) {
-        histos.fill(HIST("ReassignedMFTtrackAmbDegree"), reassoMftTrack.ambDegree());
-      }
-      if (itrack.collisionId() == reassoMftTrack.bestCollisionId()) {
-        histos.fill(HIST("AssignedMFTtrackAmbDegree"), reassoMftTrack.ambDegree());
-      }
       o2::math_utils::bringTo02Pi(phi);
       lrmfttracks(lrcollision.lastIndex(),
                   reassoMftTrack.ambDegree(),
@@ -497,7 +489,9 @@ struct LongrangeMaker {
                   phi,
                   itrack.nClusters(),
                   reassoMftTrack.bestDCAXY(),
-                  reassoMftTrack.bestDCAZ());
+                  reassoMftTrack.bestDCAZ(),
+                  itrack.isCA(),
+                  itrack.collisionId() != reassoMftTrack.bestCollisionId());
     }
 
     // v0 loop
@@ -664,20 +658,16 @@ struct LongrangeMaker {
         }
         auto phi = itrack.phi();
         o2::math_utils::bringTo02Pi(phi);
-        if (itrack.collisionId() != reassoMftTrack.bestCollisionId()) {
-          histos.fill(HIST("ReassignedMFTtrackAmbDegree"), reassoMftTrack.ambDegree());
-        }
-        if (itrack.collisionId() == reassoMftTrack.bestCollisionId()) {
-          histos.fill(HIST("AssignedMFTtrackAmbDegree"), reassoMftTrack.ambDegree());
-        }
         upclrmfttracks(upclrcollision.lastIndex(),
                        reassoMftTrack.ambDegree(),
                        itrack.pt(),
                        itrack.eta(),
                        phi,
                        itrack.nClusters(),
                        reassoMftTrack.bestDCAXY(),
-                       reassoMftTrack.bestDCAZ());
+                       reassoMftTrack.bestDCAZ(),
+                       itrack.isCA(),
+                       itrack.collisionId() != reassoMftTrack.bestCollisionId());
       }
 
       // v0 loop
@@ -869,20 +859,16 @@ struct LongrangeMaker {
         }
         auto phi = itrack.phi();
         o2::math_utils::bringTo02Pi(phi);
-        if (itrack.collisionId() != reassoMftTrack.bestCollisionId()) {
-          histos.fill(HIST("ReassignedMFTtrackAmbDegree"), reassoMftTrack.ambDegree());
-        }
-        if (itrack.collisionId() == reassoMftTrack.bestCollisionId()) {
-          histos.fill(HIST("AssignedMFTtrackAmbDegree"), reassoMftTrack.ambDegree());
-        }
         lrmfttracks(lrcollision.lastIndex(),
                     reassoMftTrack.ambDegree(),
                     itrack.pt(),
                     itrack.eta(),
                     phi,
                     itrack.nClusters(),
                     reassoMftTrack.bestDCAXY(),
-                    reassoMftTrack.bestDCAZ());
+                    reassoMftTrack.bestDCAZ(),
+                    itrack.isCA(),
+                    itrack.collisionId() != reassoMftTrack.bestCollisionId());
       }
 
       for (const auto& particle : mcparticles) {
@@ -1034,6 +1020,12 @@ struct LongrangeMaker {
         if (!isMftTrackSelected(track)) {
           continue;
         }
+        if (cfgmfttrksel.cfgRequireCA && !track.isCA()) {
+          continue;
+        }
+        if (cfgmfttrksel.cfgRequireLTF && track.isCA()) {
+          continue;
+        }
         if (!track.has_mcParticle())
           continue;
         auto particle = track.mcParticle();
@@ -1247,9 +1239,6 @@ struct LongrangeMaker {
     if constexpr (fillHis) {
       histos.fill(HIST("hMftTrkSel"), 5);
     }
-    if (cfgmfttrksel.cfgRequireCA && !track.isCA()) {
-      return false;
-    }
     if constexpr (fillHis) {
       histos.fill(HIST("hMftTrkSel"), 6);
     }