Feature; implemented manual track cuts

Author: Rafael Manhart

PWGLF/Tasks/NucleiHistTask.cxx

@@ -37,6 +37,7 @@ using namespace o2;
 using namespace o2::framework;
 using namespace o2::framework::expressions;
 
+
 namespace o2::aod
 {
 namespace NucleiTableHist
@@ -236,9 +237,23 @@ struct NucleiHistTask {
   Configurable<float> cfgCutEta{"cfgCutEta", 0.8f, "Eta range for tracks"};
   Configurable<float> nsigmacutLow{"nsigmacutLow", -3.0, "Value of the Nsigma cut"};
   Configurable<float> nsigmacutHigh{"nsigmacutHigh", +3.0, "Value of the Nsigma cut"};
+  
+  // Replacement for globalTrack filter
+  Configurable<float> minNCrossedRowsTPC{"minNCrossedRowsTPC", 70.0f, "min number of crossed rows TPC"};
+  Configurable<float> minRatioCrossedRowsTPC{"minRatioCrossedRowsTPC", 0.8f, "min ratio of crossed rows over findable clusters TPC"};
+  Configurable<float> maxChi2TPC{"maxChi2TPC", 4.0f, "max chi2 per cluster TPC"};
+  Configurable<float> maxChi2ITS{"maxChi2ITS", 36.0f, "max chi2 per cluster ITS"};
+  Configurable<float> reqTPCrefit{"reqTPCrefit", true, "require TPC refit"};
+  Configurable<float> reqITSrefit{"reqITSrefit", true, "require ITS refit"};
+  Configurable<float> maxDCA_Z{"maxDCA_Z", 2.0f, "max DCA to vertex z"};
+  Configurable<float> maxDCA_XY{"maxDCA_XY", 0.5f, "max DCA to vertex xy"};
+  Configurable<float> minReqClusterITS{"minReqClusterITS", 1.0, "min number of clusters required in ITS"};
+  Configurable<float> pTmin{"pTmin", 0.1f, "min pT"};
+  Configurable<float> pTmax{"pTmax", 1e+10f, "max pT"};
+
 
   Filter collisionFilter = nabs(aod::collision::posZ) < cfgCutVertex;
-  Filter trackFilter = (nabs(aod::track::eta) < cfgCutEta) && (requireGlobalTrackInFilter());
+  Filter trackFilter = (nabs(aod::track::eta) < cfgCutEta);// && (requireGlobalTrackInFilter());
 
   using TrackCandidates = soa::Filtered<soa::Join<aod::Tracks, aod::TracksExtra, aod::TracksDCA, aod::pidTPCLfFullPr, aod::pidTOFFullPr, aod::pidTPCLfFullDe, aod::pidTOFFullDe, aod::pidTPCLfFullHe, aod::pidTOFFullHe, aod::TrackSelection, aod::TrackSelectionExtension, aod::TOFSignal, aod::pidTOFmass, aod::pidTOFbeta>>; // aod::ReducedTracks, aod::CentFV0As
 
@@ -257,6 +272,29 @@ struct NucleiHistTask {
     spectra.fill(HIST("histRecVtxZData"), collision.posZ());
 
     for (auto track : tracks) { // start loop over tracks
+    
+      float TPCnumberCls = track.tpcNClsFound();
+      float RatioCrossedRowsOverFindableTPC = track.tpcCrossedRowsOverFindableCls();
+      float Chi2perClusterTPC = track.tpcChi2NCl();
+      float Chi2perClusterITS = track.itsChi2NCl();
+    
+      if ( TPCnumberCls < minNCrossedRowsTPC || RatioCrossedRowsOverFindableTPC < minRatioCrossedRowsTPC || Chi2perClusterTPC > maxChi2TPC || Chi2perClusterITS > maxChi2ITS || !(track.passedTPCRefit()) || !(track.passedITSRefit()) ||  track.itsNCls() < minReqClusterITS || track.pt() < pTmin || track.pt() > pTmax) {
+        continue;
+      }
+      
+      if (track.sign() > 0) {
+        spectra.fill(HIST("histDcaVsPtData_particle"), track.pt(), track.dcaXY());
+        spectra.fill(HIST("histDcaZVsPtData_particle"), track.pt(), track.dcaZ());
+      }
+      
+      if (track.sign() < 0) {
+        spectra.fill(HIST("histDcaVsPtData_antiparticle"), track.pt(), track.dcaXY());
+        spectra.fill(HIST("histDcaZVsPtData_antiparticle"), track.pt(), track.dcaZ());
+      }
+      
+      if (TMath::Abs(track.dcaXY()) > maxDCA_XY || TMath::Abs(track.dcaZ()) > maxDCA_Z) {
+        continue;
+      }
 
       // cut on rapidity
       TLorentzVector lorentzVector_proton{};
@@ -283,13 +321,12 @@ struct NucleiHistTask {
       spectra.fill(HIST("histNClusterITS"), track.pt(), track.itsNCls());
       spectra.fill(HIST("histChi2TPC"), track.pt(), track.tpcChi2NCl());
       spectra.fill(HIST("histChi2ITS"), track.pt(), track.itsChi2NCl());
+      
 
       if (track.sign() > 0) {
         proton_erg.fill(HIST("histTpcNsigmaData"), track.pt(), nSigmaProton);
         deuteron_reg.fill(HIST("histTpcNsigmaData"), track.pt(), nSigmaDeut);
         Helium3_reg.fill(HIST("histTpcNsigmaData"), track.pt() * 2.0, nSigmaHe3);
-        spectra.fill(HIST("histDcaVsPtData_particle"), track.pt(), track.dcaXY());
-        spectra.fill(HIST("histDcaZVsPtData_particle"), track.pt(), track.dcaZ());
         /*
                 if (track.centFV0A() > 0.0 && track.centFV0A() < 5.0) {
                   proton_erg.fill(HIST("histTpcNsigmaData_cent_0-5"), track.pt(), nSigmaProton);
@@ -329,8 +366,6 @@ struct NucleiHistTask {
         aproton_erg.fill(HIST("histTpcNsigmaData"), track.pt(), nSigmaProton);
         adeuteron_reg.fill(HIST("histTpcNsigmaData"), track.pt(), nSigmaDeut);
         aHelium3_reg.fill(HIST("histTpcNsigmaData"), track.pt() * 2.0, nSigmaHe3);
-        spectra.fill(HIST("histDcaVsPtData_antiparticle"), track.pt(), track.dcaXY());
-        spectra.fill(HIST("histDcaZVsPtData_antiparticle"), track.pt(), track.dcaZ());
         /*
                 if (track.centFV0A() > 0.0 && track.centFV0A() < 5.0) {
                   aproton_erg.fill(HIST("histTpcNsigmaData_cent_0-5"), track.pt(), nSigmaProton);