[PWGLF] fix missing abs() and config-axes bugs

PWGLF/Tasks/Strangeness/cascadecorrelations.cxx

@@ -15,38 +15,40 @@
 // Author: Rik Spijkers (rik.spijkers@cern.ch)
 //
 
-#include <Math/Vector4D.h>
-#include <cmath>
-#include <array>
-#include <cstdlib>
-#include <string>
-#include <utility>
+#include "PWGLF/DataModel/LFStrangenessTables.h"
+#include "PWGLF/Utils/inelGt.h"
 
-#include "Framework/runDataProcessing.h"
-#include "Framework/AnalysisTask.h"
-#include "Framework/AnalysisDataModel.h"
-#include "Framework/ASoAHelpers.h"
-#include "Framework/O2DatabasePDGPlugin.h"
-#include "ReconstructionDataFormats/Track.h"
 #include "Common/Core/RecoDecay.h"
-#include "Common/Core/trackUtilities.h"
-#include "PWGLF/DataModel/LFStrangenessTables.h"
 #include "Common/Core/TrackSelection.h"
-#include "Common/DataModel/TrackSelectionTables.h"
-#include "Common/DataModel/EventSelection.h"
+#include "Common/Core/trackUtilities.h"
 #include "Common/DataModel/Centrality.h"
+#include "Common/DataModel/EventSelection.h"
 #include "Common/DataModel/Multiplicity.h"
 #include "Common/DataModel/PIDResponse.h"
-#include "CCDB/BasicCCDBManager.h"
+#include "Common/DataModel/TrackSelectionTables.h"
 #include "EventFiltering/Zorro.h"
-#include "PWGLF/Utils/inelGt.h"
 
+#include "CCDB/BasicCCDBManager.h"
+#include "Framework/ASoAHelpers.h"
+#include "Framework/AnalysisDataModel.h"
+#include "Framework/AnalysisTask.h"
+#include "Framework/O2DatabasePDGPlugin.h"
+#include "Framework/runDataProcessing.h"
+#include "ReconstructionDataFormats/Track.h"
+
+#include <Math/Vector4D.h>
 #include <TFile.h>
-#include <TList.h>
 #include <TH2F.h>
-#include <TProfile.h>
+#include <TList.h>
 #include <TLorentzVector.h>
 #include <TPDGCode.h>
+#include <TProfile.h>
+
+#include <array>
+#include <cmath>
+#include <cstdlib>
+#include <string>
+#include <utility>
 // #include <TDatabasePDG.h>
 
 using namespace o2;
@@ -275,7 +277,7 @@
     if (!gen.isPhysicalPrimary())
       return;
     int genpdg = gen.pdgCode();
     if ((flag < 3 && TMath::Abs(genpdg) == 3312) || (flag > 1 && TMath::Abs(genpdg) == 3334)) {
       // if casc is consistent with Xi and has matched gen Xi OR cand is consistent with Omega and has matched gen omega
       // have to do this in case we reco true Xi with only Omega hypothesis (or vice versa) (very unlikely)
       registry.fill(HIST("truerec/hV0Radius"), rec.v0radius());
@@ -366,18 +368,18 @@
         casc.v0cosPA(pvx, pvy, pvz) < v0setting_cospa ||
         casc.casccosPA(pvx, pvy, pvz) < cascadesetting_cospa ||
         casc.dcav0topv(pvx, pvy, pvz) < cascadesetting_mindcav0topv ||
         TMath::Abs(casc.mLambda() - 1.115683) > cascadesetting_v0masswindow)
       return 0; // It failed at least one topo selection
 
     registry.fill(HIST("hSelectionStatus"), 3); // passes topo
     // registry.fill(HIST("hMassXi3"), casc.mXi(), casc.pt());
 
     if (TMath::Abs(posTrack.eta()) > etaTracks || TMath::Abs(negTrack.eta()) > etaTracks || TMath::Abs(bachTrack.eta()) > etaTracks)
       return 0;
 
     registry.fill(HIST("hSelectionStatus"), 4); // passes track eta
 
     if (TMath::Abs(casc.eta()) > etaCascades)
       return 0;
 
     registry.fill(HIST("hSelectionStatus"), 5); // passes candidate eta
@@ -388,25 +390,25 @@
     // Lambda check
     if (casc.sign() < 0) {
       // Proton check:
       if (TMath::Abs(posTrack.tpcNSigmaPr()) > tpcNsigmaProton)
         return 0;
       // Pion check:
       if (TMath::Abs(negTrack.tpcNSigmaPi()) > tpcNsigmaPion)
         return 0;
     } else {
       // Proton check:
       if (TMath::Abs(negTrack.tpcNSigmaPr()) > tpcNsigmaProton)
         return 0;
       // Pion check:
       if (TMath::Abs(posTrack.tpcNSigmaPi()) > tpcNsigmaPion)
         return 0;
     }
     registry.fill(HIST("hSelectionStatus"), 6); // passes V0 daughters PID
     // registry.fill(HIST("hMassXi4"), casc.mXi(), casc.pt());
 
     // Bachelor check
     if (TMath::Abs(bachTrack.tpcNSigmaPi()) < tpcNsigmaBachelor) {
       if (TMath::Abs(bachTrack.tpcNSigmaKa()) < tpcNsigmaBachelor) {
         // consistent with both!
         registry.fill(HIST("hSelectionStatus"), 7); // passes bach PID
         // registry.fill(HIST("hMassXi5"), casc.mXi(), casc.pt());
@@ -536,8 +538,8 @@
   ConfigurableAxis dcaAxis = {"dcaAxis", {100, 0.0f, 2.0f}, "cm"};
   ConfigurableAxis multiplicityAxis{"multiplicityAxis", {100, 0, 100}, "Multiplicity (MultFT0M?)"};
   ConfigurableAxis invLambdaMassAxis{"invLambdaMassAxis", {100, 1.07f, 1.17f}, "Inv. Mass (GeV/c^{2})"};
-  AxisSpec deltaYAxis{40, -2 * maxRapidity, 2 * maxRapidity, "#Delta y"};
-  AxisSpec rapidityAxis{100, -maxRapidity, maxRapidity, "y"};
+  ConfigurableAxis deltaYAxis{"deltaYAxis", {40, -2.f, 2.f}, "#Delta y"};
+  ConfigurableAxis rapidityAxis{"rapidityAxis", {100, -1.f, 1.f}, "y"};
   AxisSpec selectionFlagAxis{4, -0.5f, 3.5f, "Selection flag of casc candidate"};
   AxisSpec itsClustersAxis{8, -0.5, 7.5, "number of ITS clusters"};
   AxisSpec tpcRowsAxis{160, -0.5, 159.5, "TPC crossed rows"};
@@ -1124,7 +1126,7 @@
 
       if (!trigger.isPhysicalPrimary() || !assoc.isPhysicalPrimary())
         continue; // require the cascades to be primaries
-      if (trigger.eta() > etaGenCascades)
+      if (TMath::Abs(trigger.eta()) > etaGenCascades)
         continue; // only apply eta cut to trigger - trigger normalization still valid without introducing 2-particle-acceptance effects
 
       double dphi = RecoDecay::constrainAngle(trigger.phi() - assoc.phi(), -PIHalf);