Test histogram for TOF nSigma with MC

Author: BanajitBarman

PWGLF/Tasks/Nuspex/spectraTOF.cxx

@@ -591,6 +591,7 @@ struct tofSpectra {
       for (int par = 2; par <= 4; par++) {
         for (int i = 0; i < NpCharge; i++) {
           hMCpdg_nsigmaTPC[par - 2][i] = histos.add<TH3>(Form("test_mclabels/nsigmatpc/%s/%s/pdg_%i", (i < Np) ? "pos" : "neg", pN[par], PDGs[i % Np]), Form("True %s (%i) in %s selection", pTCharge[i], PDGs[i], (i < Np) ? pTCharge[par] : pTCharge[par + Np]), kTH3D, {ptAxis, nsigmaTPCAxisOccupancy, multAxis});
+          hMCpdg_nsigmaTOF[par - 2][i] = histos.add<TH3>(Form("test_mclabels/nsigmatof/%s/%s/pdg_%i", (i < Np) ? "pos" : "neg", pN[par], PDGs[i % Np]), Form("True %s (%i) in %s selection", pTCharge[i], PDGs[i], (i < Np) ? pTCharge[par] : pTCharge[par + Np]), kTH3D, {ptAxis, nsigmaTOFAxisOccupancy, multAxis});
         }
       }
     }
@@ -2817,14 +2818,17 @@ struct tofSpectra {
       int pdgCode = mcParticle.pdgCode();
       static_for<2, 4>([&](auto par) {
         const auto& nsigmaTPCpar = o2::aod::pidutils::tpcNSigma<par>(track);
+        const auto& nsigmaTOFpar = o2::aod::pidutils::tofNSigma<par>(track);
         bool isTPCpar = std::abs(nsigmaTPCpar) < trkselOptions.cfgCutNsigma;
+        bool isTOFpar = std::abs(nsigmaTOFpar) < trkselOptions.cfgCutNsigma;
         // Precompute rapidity values to avoid redundant calculations
         double rapiditypar = std::abs(track.rapidity(PID::getMass(par)));
         // TPC Selection and histogram filling
         if (isTPCpar && rapiditypar <= trkselOptions.cfgCutY) {
           static_for<0, 17>([&](auto i) {
             if (pdgCode == PDGs[i]) {
               hMCpdg_nsigmaTPC[par - 2][i]->Fill(track.pt(), nsigmaTPCpar, multiplicity);
+              hMCpdg_nsigmaTOF[par - 2][i]->Fill(track.pt(), nsigmaTOFpar, multiplicity);
             }
           });
         }