PWGLF: NucleiTask - Add cluster cut (#1619)

* PWGLF: NucleiTask - Add cluster cut

* Minor format fix

* Fix MegaLinter issue

Author: giovannimalfattore

PWGLF/DataModel/LFNucleiTables.h

@@ -84,6 +84,7 @@ DECLARE_SOA_COLUMN(Beta, beta, float);
 // TPC and ITS QA
 DECLARE_SOA_COLUMN(TPCNClsCrossedRows, tpcNClsCrossedRows, int16_t);
 DECLARE_SOA_COLUMN(TPCCrossedRowsOverFindableCls, tpcCrossedRowsOverFindableCls, float);
+DECLARE_SOA_COLUMN(TPCNClsFound, tpcNClsFound, int16_t);
 DECLARE_SOA_COLUMN(TPCChi2Ncl, tpcChi2NCl, float);
 DECLARE_SOA_COLUMN(ITSChi2NCl, itsChi2NCl, float);
 // For MC
@@ -126,6 +127,7 @@ DECLARE_SOA_TABLE(LfCandNucleusFull, "AOD", "LFNUCL",
                   full::Sign,
                   full::TPCNClsCrossedRows,
                   full::TPCCrossedRowsOverFindableCls,
+                  full::TPCNClsFound,
                   full::TPCChi2Ncl,
                   full::ITSChi2NCl);
 DECLARE_SOA_TABLE(LfCandNucleusMC, "AOD", "LFNUCLMC",

PWGLF/TableProducer/LFTreeCreatorNuclei.cxx

@@ -20,6 +20,12 @@
 
 #include "PWGLF/DataModel/LFNucleiTables.h"
 
+#include <cmath>
+
+#include <TLorentzVector.h>
+#include <TMath.h>
+#include <TObjArray.h>
+
 #include "ReconstructionDataFormats/Track.h"
 #include "Framework/runDataProcessing.h"
 #include "Framework/AnalysisTask.h"
@@ -35,11 +41,6 @@
 #include "Common/DataModel/TrackSelectionTables.h"
 #include "Common/Core/trackUtilities.h"
 
-#include <TLorentzVector.h>
-#include <TMath.h>
-#include <TObjArray.h>
-#include <cmath>
-
 using namespace o2;
 using namespace o2::framework;
 using namespace o2::framework::expressions;
@@ -133,6 +134,7 @@ struct LfTreeCreatorNuclei {
         track.sign(),
         track.tpcNClsCrossedRows(),
         track.tpcCrossedRowsOverFindableCls(),
+        track.tpcNClsFound(),
         track.tpcChi2NCl(),
         track.itsChi2NCl());
 

PWGLF/Tasks/LFNucleiBATask.cxx

@@ -46,6 +46,8 @@ struct LFNucleiBATask {
   Configurable<float> etaCut{"etaCut", 0.8f, "Value of the eta selection for spectra (default 0.8)"};
   Configurable<float> yCut{"yCut", 0.5f, "Value of the rapidity selection for spectra (default 0.5)"};
   Configurable<float> cfgCutVertex{"cfgCutVSertex", 10.0f, "Accepted z-vertex range"};
+  Configurable<float> cfgCutTPCXRows{"cfgCutTPCXRows", -1.f, "Minimum number of crossed TPC rows"};
+  Configurable<float> cfgCutTPCClusters{"cfgCutTPCClusters", -1.f, "Minimum number of found TPC clusters"};
   ConfigurableAxis binsPt{"binsPt", {VARIABLE_WIDTH, 0.0, 0.1, 0.12, 0.14, 0.16, 0.18, 0.2, 0.25, 0.3, 0.35, 0.4, 0.45, 0.5, 0.55, 0.6, 0.65, 0.7, 0.75, 0.8, 0.85, 0.9, 0.95, 1.0, 1.1, 1.2, 1.3, 1.4, 1.5, 1.6, 1.7, 1.8, 2.0, 2.2, 2.4, 2.6, 2.8, 3.0, 3.2, 3.4, 3.6, 3.8, 4.0, 4.2, 4.4, 4.6, 4.8, 5.0}, ""};
   ConfigurableAxis binsdEdx{"binsdEdx", {1000, 0.f, 1000.f}, ""};
   Configurable<bool> enableEvTimeSplitting{"enableEvTimeSplitting", true, "Flag to enable histograms splitting depending on the Event Time used"};
@@ -654,7 +656,11 @@ struct LFNucleiBATask {
           continue;
         }
       }
-      // QUI DCA
+
+      if (track.tpcNClsCrossedRows() < cfgCutTPCXRows)
+        continue;
+      if (track.tpcNClsFound() < cfgCutTPCClusters)
+        continue;
 
       // Tracks DCA histos fill
       histos.fill(HIST("tracks/hDCAxy"), track.dcaXY());
@@ -714,7 +720,7 @@ struct LFNucleiBATask {
           continue;
         }
       }
-      // QUI TUTTO IL RESTO
+
       // LOG(info)<<"\n collisionId ============>"<<track.collisionId();
 
       // QA histos fill
@@ -753,13 +759,6 @@ struct LFNucleiBATask {
         }
       }
 
-      // Calculate on fly NSigmaTPC* calibrated from BB parameters for light antinuclei
-      /*
-      To be included
-
-      return NSigmaTPC[De,Tr,He,Al]Cal
-      */
-
       // Tracks histos fill
       histos.fill(HIST("tracks/h1Eta"), track.eta());
       histos.fill(HIST("tracks/h1VarPhi"), track.phi());