merge master into kfstudy

ALICE3/TableProducer/OTF/onTheFlyTOFPID.cxx

@@ -276,10 +276,10 @@ struct OnTheFlyTOFPID {
 
       if (scalarProduct1 > scalarProduct2) {
         modulus = std::hypot(point1[0] - trcCircle.xC, point1[1] - trcCircle.yC) * std::hypot(startPoint[0] - trcCircle.xC, startPoint[1] - trcCircle.yC);
-        cosAngle = (point1[0] - trcCircle.xC) * (startPoint[0] - trcCircle.xC) + (point1[1] - trcCircle.yC) * (startPoint[0] - trcCircle.yC);
+        cosAngle = (point1[0] - trcCircle.xC) * (startPoint[0] - trcCircle.xC) + (point1[1] - trcCircle.yC) * (startPoint[1] - trcCircle.yC);
       } else {
         modulus = std::hypot(point2[0] - trcCircle.xC, point2[1] - trcCircle.yC) * std::hypot(startPoint[0] - trcCircle.xC, startPoint[1] - trcCircle.yC);
-        cosAngle = (point2[0] - trcCircle.xC) * (startPoint[0] - trcCircle.xC) + (point2[1] - trcCircle.yC) * (startPoint[0] - trcCircle.yC);
+        cosAngle = (point2[0] - trcCircle.xC) * (startPoint[0] - trcCircle.xC) + (point2[1] - trcCircle.yC) * (startPoint[1] - trcCircle.yC);
       }
       cosAngle /= modulus;
       length = trcCircle.rC * TMath::ACos(cosAngle);

CODEOWNERS

@@ -29,7 +29,7 @@
 /PWGDQ @alibuild @iarsene @dsekihat @feisenhu @lucamicheletti93
 /PWGEM @alibuild @mikesas @rbailhac @feisenhu
 /PWGEM/Dilepton @alibuild @mikesas @rbailhac @dsekihat @ivorobye @feisenhu
-/PWGEM/PhotonMeson @alibuild @mikesas @rbailhac @m-c-danisch @novitzky @mhemmer-cern
+/PWGEM/PhotonMeson @alibuild @mikesas @rbailhac @m-c-danisch @novitzky @mhemmer-cern @dsekihat
 /PWGHF @alibuild @vkucera @fcolamar @fgrosa @fcatalan92 @mfaggin @mmazzilli @deepathoms @nzardosh @NicoleBastid
 /PWGLF @alibuild @lramona @alcaliva @lbariogl @chiarapinto @BongHwi @lbarnby @mbombara @iravasen @njacazio @ChiaraDeMartin95 @skundu692
 /PWGMM @alibuild @aalkin

Common/DataModel/PIDResponse.h

@@ -46,6 +46,13 @@ void packInTable(const float& valueToBin, T& table)
   }
 }
 
+// Function to unpack a binned value into a float
+template <typename binningType>
+float unPackInTable(const typename binningType::binned_t& valueToUnpack)
+{
+  return binningType::bin_width * static_cast<float>(valueToUnpack);
+}
+
 // Checkers for TOF PID hypothesis availability (runtime)
 template <class T>
 using hasTOFEl = decltype(std::declval<T&>().tofNSigmaEl());
@@ -481,7 +488,7 @@ DECLARE_SOA_COLUMN(TOFNSigmaAl, tofNSigmaAl, float); //! Nsigma separation with
 // Macro to convert the stored Nsigmas to floats
 #define DEFINE_UNWRAP_NSIGMA_COLUMN(COLUMN, COLUMN_NAME) \
   DECLARE_SOA_DYNAMIC_COLUMN(COLUMN, COLUMN_NAME,        \
-                             [](binning::binned_t nsigma_binned) -> float { return binning::bin_width * static_cast<float>(nsigma_binned); });
+                             [](binning::binned_t nsigma_binned) -> float { return o2::aod::pidutils::unPackInTable<binning>(nsigma_binned); });
 
 namespace pidtof_tiny
 {

Common/TableProducer/trackPropagation.cxx

@@ -151,7 +151,7 @@ struct TrackPropagation {
 
     for (auto& track : tracks) {
       if constexpr (fillCovMat) {
-        if (fillTracksDCACov) {
+        if (fillTracksDCA || fillTracksDCACov) {
           mDcaInfoCov.set(999, 999, 999, 999, 999);
         }
         setTrackParCov(track, mTrackParCov);

DPG/Tasks/AOTTrack/perfK0sResolution.cxx

@@ -46,6 +46,16 @@ struct perfK0sResolution {
     const AxisSpec etaAxis{etaBins, "#eta"};
     const AxisSpec phiAxis{phiBins, "#phi"};
 
+    int nProc = 0;
+    if (doprocessData) {
+      LOG(info) << "processData enabled";
+      nProc++;
+    }
+    if (doprocessMC) {
+      LOG(info) << "processMC enabled";
+      nProc++;
+    }
+
     registry.add("h2_masspT", "h2_masspT", {HistType::kTH2F, {mAxis, pTAxis}});
     registry.add("h2_masseta", "h2_masseta", {HistType::kTH2F, {mAxis, etaAxis}});
     registry.add("h2_massphi", "h2_massphi", {HistType::kTH2F, {mAxis, phiAxis}});
@@ -62,26 +72,70 @@ struct perfK0sResolution {
   Configurable<float> nSigTPC{"nSigTPC", 10., "nSigTPC"};
   Configurable<int> trdSelectionPos{"trdSelectionPos", 0, "Flag for the TRD selection on positive daughters: -1 no TRD, 0 no selection, 1 TRD"};
   Configurable<int> trdSelectionNeg{"trdSelectionNeg", 0, "Flag for the TRD selection on negative daughters: -1 no TRD, 0 no selection, 1 TRD"};
+  Configurable<int> tofSelectionPos{"tofSelectionPos", 0, "Flag for the TOF selection on positive daughters: -1 no TOF, 0 no selection, 1 TOF"};
+  Configurable<int> tofSelectionNeg{"tofSelectionNeg", 0, "Flag for the TOF selection on negative daughters: -1 no TOF, 0 no selection, 1 TOF"};
   Configurable<bool> eventSelection{"eventSelection", true, "event selection"};
 
   template <typename T1, typename T2, typename C>
   bool acceptV0(const T1& v0, const T2& ntrack, const T2& ptrack, const C& collision)
   {
     // Apply selections on V0
-    if (TMath::Abs(v0.yK0Short()) > rapidity)
+    if (TMath::Abs(v0.yK0Short()) > rapidity) {
       return false;
-    if (v0.v0cosPA() < v0setting_cospa)
+    }
+    if (v0.v0cosPA() < v0setting_cospa) {
       return false;
-    if (v0.v0radius() < v0setting_radius)
+    }
+    if (v0.v0radius() < v0setting_radius) {
       return false;
-    if (v0.distovertotmom(collision.posX(), collision.posY(), collision.posZ()) * pid_constants::sMasses[PID::K0] > 2.684 * v0lifetime)
+    }
+    if (v0.distovertotmom(collision.posX(), collision.posY(), collision.posZ()) * pid_constants::sMasses[PID::K0] > 2.684 * v0lifetime) {
       return false;
+    }
 
     // Apply selections on V0 daughters
-    if (!ntrack.hasTPC() || !ptrack.hasTPC())
+    if (!ntrack.hasTPC() || !ptrack.hasTPC()) {
       return false;
-    if (ntrack.tpcNSigmaPi() > nSigTPC || ptrack.tpcNSigmaPi() > nSigTPC)
+    }
+    if (ntrack.tpcNSigmaPi() > nSigTPC || ptrack.tpcNSigmaPi() > nSigTPC) {
       return false;
+    }
+    // TOF selection
+    switch (tofSelectionPos) {
+      case -1:
+        if (ptrack.hasTOF()) {
+          return false;
+        }
+        break;
+      case 0:
+        break;
+      case 1:
+        if (!ptrack.hasTOF()) {
+          return false;
+        }
+        break;
+      default:
+        LOG(fatal) << "Invalid TRD selection for positive daughter";
+        break;
+    }
+    switch (tofSelectionNeg) {
+      case -1:
+        if (ntrack.hasTOF()) {
+          return false;
+        }
+        break;
+      case 0:
+        break;
+      case 1:
+        if (!ntrack.hasTOF()) {
+          return false;
+        }
+        break;
+      default:
+        LOG(fatal) << "Invalid TRD selection for negative daughter";
+        break;
+    }
+    // TRD selection
     switch (trdSelectionPos) {
       case -1:
         if (ptrack.hasTRD()) {
@@ -121,7 +175,7 @@ struct perfK0sResolution {
 
   Filter v0Filter = nabs(aod::v0data::dcapostopv) > v0setting_dcapostopv&& nabs(aod::v0data::dcanegtopv) > v0setting_dcanegtopv&& aod::v0data::dcaV0daughters < v0setting_dcav0dau;
 
-  void process(SelectedCollisions::iterator const& collision, soa::Filtered<aod::V0Datas> const& fullV0s, PIDTracks const& tracks)
+  void processData(SelectedCollisions::iterator const& collision, soa::Filtered<aod::V0Datas> const& fullV0s, PIDTracks const& tracks)
   {
     if (eventSelection && !collision.sel8())
       return;
@@ -138,10 +192,34 @@ struct perfK0sResolution {
       registry.fill(HIST("h2_massphi"), v0.mK0Short(), v0.phi());
     }
   }
+  PROCESS_SWITCH(perfK0sResolution, processData, "Process data", true);
+
+  void processMC(SelectedCollisions::iterator const& collision, soa::Filtered<aod::V0Datas> const& fullV0s, soa::Join<PIDTracks, aod::McTrackLabels> const& tracks, aod::McParticles const&)
+  {
+    if (eventSelection && !collision.sel8())
+      return;
+
+    for (auto& v0 : fullV0s) {
+
+      const auto& posTrack = v0.posTrack_as<soa::Join<PIDTracks, aod::McTrackLabels>>();
+      const auto& negTrack = v0.negTrack_as<soa::Join<PIDTracks, aod::McTrackLabels>>();
+      if (!acceptV0(v0, negTrack, posTrack, collision))
+        continue;
+      if (!posTrack.has_mcParticle()) {
+        continue;
+      }
+      if (!negTrack.has_mcParticle()) {
+        continue;
+      }
+      if (posTrack.mcParticle().pdgCode() != PID::Pion || negTrack.mcParticle().pdgCode() != PID::Pion) {
+        continue;
+      }
+      registry.fill(HIST("h2_masspT"), v0.mK0Short(), v0.pt());
+      registry.fill(HIST("h2_masseta"), v0.mK0Short(), v0.eta());
+      registry.fill(HIST("h2_massphi"), v0.mK0Short(), v0.phi());
+    }
+  }
+  PROCESS_SWITCH(perfK0sResolution, processMC, "Process MC", false);
 };
 
-WorkflowSpec defineDataProcessing(ConfigContext const& cfgc)
-{
-  return WorkflowSpec{
-    adaptAnalysisTask<perfK0sResolution>(cfgc)};
-}
+WorkflowSpec defineDataProcessing(ConfigContext const& cfgc) { return WorkflowSpec{adaptAnalysisTask<perfK0sResolution>(cfgc)}; }

DPG/Tasks/AOTTrack/qaEfficiencyV0s.cxx

@@ -36,45 +36,53 @@
 #include "THashList.h"
 
 using namespace o2::framework;
+static constexpr int nSpecies = 2; // One per PDG
+// static constexpr const char* particleTitle[nSpecies] = {"K0s", "-K0s"};
+static constexpr int PDGs[nSpecies] = {kK0Short, -kK0Short};
+int pdgToIndex(int pdg)
+{
+  for (int i = 0; i < nSpecies; i++) {
+    if (pdg == PDGs[i]) {
+      return i;
+    }
+  }
+  return -1;
+}
+static constexpr int kHistoPtNum = 0;
+static constexpr int kHistoPtDen = 1;
+static constexpr int kHistoTot = 2;
+
+std::shared_ptr<TH1> histograms[nSpecies][kHistoTot] = {{nullptr}};
+std::shared_ptr<TH1> histogramsPrm[nSpecies][kHistoTot] = {{nullptr}};
+std::shared_ptr<TH1> histogramsPrmRap[nSpecies][kHistoTot] = {{nullptr}};
 
 struct QaEfficiencyV0s {
   // Particle information
-  static constexpr int nSpecies = 1; // One per PDG
-  static constexpr const char* particleTitle[nSpecies] = {"K0s"};
-  static constexpr int PDGs[nSpecies] = {kK0Short};
-  // Particle only selection
-  Configurable<bool> doK0s{"do-k0s", false, "Flag to run with the PDG code of k0s"};
   Configurable<float> rapidityCut{"rapidityCut", 0.5, "Rapidity cut"};
   ConfigurableAxis ptBins{"ptBins", {200, 0.f, 5.f}, "Pt binning"};
   OutputObj<THashList> listEfficiencyMC{"EfficiencyMC"};
-  // Histograms
-
   HistogramRegistry registry{"Histos", {}, OutputObjHandlingPolicy::AnalysisObject};
-
   void init(InitContext&)
   {
     const AxisSpec axisPt{ptBins, "#it{p}_{T} GeV/#it{c}"};
-    auto h = registry.add<TH1>("Pos/PtNum_310", "Pos/Num_310", kTH1F, {axisPt});
-    registry.add("Pos/PtDen_310", "Pos/Den_310", kTH1F, {axisPt});
-    registry.add("Neg/PtNum_310", "Neg/Num_310", kTH1F, {axisPt});
-    registry.add("Neg/PtDen_310", "Neg/Den_310", kTH1F, {axisPt});
-    registry.add("Prm/Pos/PtNum_310", "Prm/Pos/Num_310", kTH1F, {axisPt});
-    registry.add("Prm/Pos/PtDen_310", "Prm/Pos/Den_310", kTH1F, {axisPt});
-    registry.add("Prm/Neg/PtNum_310", "Prm/Neg/Num_310", kTH1F, {axisPt});
-    registry.add("Prm/Neg/PtDen_310", "Prm/Neg/Den_310", kTH1F, {axisPt});
-    registry.add("PrmRap/Pos/PtNum_310", "PrmRap/Pos/Num_310", kTH1F, {axisPt});
-    registry.add("PrmRap/Pos/PtDen_310", "PrmRap/Pos/Den_310", kTH1F, {axisPt});
-    registry.add("PrmRap/Neg/PtNum_310", "PrmRap/Neg/Num_310", kTH1F, {axisPt});
-    registry.add("PrmRap/Neg/PtDen_310", "PrmRap/Neg/Den_310", kTH1F, {axisPt});
-
-    TAxis* axis = h->GetXaxis();
     listEfficiencyMC.setObject(new THashList);
-    listEfficiencyMC->Add(new TEfficiency(Form("efficiencyPt_pdg%d", PDGs[0]), Form("efficiencyPt_pdg%d", PDGs[0]), axis->GetNbins(), axis->GetXmin(), axis->GetXmax()));
-    listEfficiencyMC->Add(new TEfficiency(Form("efficiencyPt_pdg-%d", PDGs[0]), Form("efficiencyPt_pdg-%d", PDGs[0]), axis->GetNbins(), axis->GetXmin(), axis->GetXmax()));
-    listEfficiencyMC->Add(new TEfficiency(Form("efficiencyPtPrm_pdg%d", PDGs[0]), Form("efficiencyPtPrm_pdg%d", PDGs[0]), axis->GetNbins(), axis->GetXmin(), axis->GetXmax()));
-    listEfficiencyMC->Add(new TEfficiency(Form("efficiencyPtPrm_pdg-%d", PDGs[0]), Form("efficiencyPtPrm_pdg-%d", PDGs[0]), axis->GetNbins(), axis->GetXmin(), axis->GetXmax()));
-    listEfficiencyMC->Add(new TEfficiency(Form("efficiencyPtPrmRap_pdg%d", PDGs[0]), Form("efficiencyPtPrmRap_pdg%d", PDGs[0]), axis->GetNbins(), axis->GetXmin(), axis->GetXmax()));
-    listEfficiencyMC->Add(new TEfficiency(Form("efficiencyPtPrmRap_pdg-%d", PDGs[0]), Form("efficiencyPtPrmRap_pdg-%d", PDGs[0]), axis->GetNbins(), axis->GetXmin(), axis->GetXmax()));
+
+    for (int i = 0; i < nSpecies; i++) {
+      // MC efficiency (PDG code
+      histograms[i][kHistoPtNum] = registry.add<TH1>(Form("Pt/Num_%i", PDGs[i]), Form("Num %i", PDGs[i]), kTH1F, {axisPt});
+      histograms[i][kHistoPtDen] = registry.add<TH1>(Form("Pt/Den_%i", PDGs[i]), Form("Den %i", PDGs[i]), kTH1F, {axisPt});
+
+      histogramsPrm[i][kHistoPtNum] = registry.add<TH1>(Form("Pt/Prm/Num_%i", PDGs[i]), Form("Pt Prm Num %i", PDGs[i]), kTH1F, {axisPt});
+      histogramsPrm[i][kHistoPtDen] = registry.add<TH1>(Form("Pt/Prm/Den_%i", PDGs[i]), Form("Pt Prm Den %i", PDGs[i]), kTH1F, {axisPt});
+
+      histogramsPrmRap[i][kHistoPtNum] = registry.add<TH1>(Form("Pt/Prm/Rap/PtNum_%i", PDGs[i]), Form("Pt Prm Rap Num %i", PDGs[i]), kTH1F, {axisPt});
+      histogramsPrmRap[i][kHistoPtDen] = registry.add<TH1>(Form("Pt/Prm/Rap/PtDen_%i", PDGs[i]), Form("Pt Prm Rap Den %i", PDGs[i]), kTH1F, {axisPt});
+
+      TAxis* axis = histograms[i][0]->GetXaxis();
+      listEfficiencyMC->Add(new TEfficiency(Form("efficiencyPt_pdg%d", PDGs[i]), Form("efficiencyPt_pdg%d", PDGs[i]), axis->GetNbins(), axis->GetXmin(), axis->GetXmax()));
+      listEfficiencyMC->Add(new TEfficiency(Form("efficiencyPtPrm_pdg%d", PDGs[i]), Form("efficiencyPtPrm_pdg%d", PDGs[i]), axis->GetNbins(), axis->GetXmin(), axis->GetXmax()));
+      listEfficiencyMC->Add(new TEfficiency(Form("efficiencyPtPrmRap_pdg%d", PDGs[i]), Form("efficiencyPtPrmRap_pdg%d", PDGs[i]), axis->GetNbins(), axis->GetXmin(), axis->GetXmax()));
+    }
   }
 
   // MC process
@@ -83,62 +91,44 @@ struct QaEfficiencyV0s {
   void process(o2::aod::McV0Labels const& V0s,
                o2::aod::McParticles const& mcParticles)
   {
-    for (auto const& v0 : V0s) {
+    for (auto const& v0 : V0s) { // Numerator
       if (v0.has_mcParticle()) {
         auto mcparticle = v0.mcParticle();
-        if (mcparticle.pdgCode() == PDGs[0]) {
-          registry.fill(HIST("Pos/PtNum_310"), mcparticle.pt());
-          if (mcparticle.isPhysicalPrimary()) {
-            registry.fill(HIST("Prm/Pos/PtNum_310"), mcparticle.pt());
-            if (TMath::Abs(mcparticle.y()) < rapidityCut) {
-              registry.fill(HIST("PrmRap/Pos/PtNum_310"), mcparticle.pt());
-            }
-          }
-        } else if (mcparticle.pdgCode() == -PDGs[0]) {
-          registry.fill(HIST("Neg/PtNum_310"), mcparticle.pt());
-          if (mcparticle.isPhysicalPrimary()) {
-            registry.fill(HIST("Prm/Neg/PtNum_310"), mcparticle.pt());
-            if (TMath::Abs(mcparticle.y()) < rapidityCut) {
-              registry.fill(HIST("PrmRap/Neg/PtNum_310"), mcparticle.pt());
-            }
-          }
+        const auto index = pdgToIndex(mcparticle.pdgCode());
+        if (index < 0) {
+          continue;
         }
-      }
-    }
-    for (auto const& mcparticle : mcParticles) {
-      if (mcparticle.pdgCode() == PDGs[0]) {
-        registry.fill(HIST("Pos/PtDen_310"), mcparticle.pt());
+        histograms[index][kHistoPtNum]->Fill(mcparticle.pt());
         if (mcparticle.isPhysicalPrimary()) {
-          registry.fill(HIST("Prm/Pos/PtDen_310"), mcparticle.pt());
+          histogramsPrm[index][kHistoPtNum]->Fill(mcparticle.pt());
           if (TMath::Abs(mcparticle.y()) < rapidityCut) {
-            registry.fill(HIST("PrmRap/Pos/PtDen_310"), mcparticle.pt());
+            histogramsPrmRap[index][kHistoPtNum]->Fill(mcparticle.pt());
           }
         }
-      } else if (mcparticle.pdgCode() == -PDGs[0]) {
-        registry.fill(HIST("Neg/PtDen_310"), mcparticle.pt());
-        if (mcparticle.isPhysicalPrimary()) {
-          registry.fill(HIST("Prm/Neg/PtDen_310"), mcparticle.pt());
-          if (TMath::Abs(mcparticle.y()) < rapidityCut) {
-            registry.fill(HIST("PrmRap/Neg/PtDen_310"), mcparticle.pt());
-          }
+      }
+    }
+    for (auto const& mcparticle : mcParticles) { // Denominator
+      const auto index = pdgToIndex(mcparticle.pdgCode());
+      if (index < 0) {
+        continue;
+      }
+      histograms[index][kHistoPtDen]->Fill(mcparticle.pt());
+      if (mcparticle.isPhysicalPrimary()) {
+        histogramsPrm[index][kHistoPtDen]->Fill(mcparticle.pt());
+        if (TMath::Abs(mcparticle.y()) < rapidityCut) {
+          histogramsPrmRap[index][kHistoPtDen]->Fill(mcparticle.pt());
         }
       }
     }
 
-    static_cast<TEfficiency*>(listEfficiencyMC->At(0))->SetPassedHistogram(*registry.get<TH1>(HIST("Pos/PtNum_310")), "f");
-    static_cast<TEfficiency*>(listEfficiencyMC->At(0))->SetTotalHistogram(*registry.get<TH1>(HIST("Pos/PtDen_310")), "f");
-    static_cast<TEfficiency*>(listEfficiencyMC->At(1))->SetPassedHistogram(*registry.get<TH1>(HIST("Neg/PtNum_310")), "f");
-    static_cast<TEfficiency*>(listEfficiencyMC->At(1))->SetTotalHistogram(*registry.get<TH1>(HIST("Neg/PtDen_310")), "f");
-
-    static_cast<TEfficiency*>(listEfficiencyMC->At(2))->SetPassedHistogram(*registry.get<TH1>(HIST("Prm/Pos/PtNum_310")), "f");
-    static_cast<TEfficiency*>(listEfficiencyMC->At(2))->SetTotalHistogram(*registry.get<TH1>(HIST("Prm/Pos/PtDen_310")), "f");
-    static_cast<TEfficiency*>(listEfficiencyMC->At(3))->SetPassedHistogram(*registry.get<TH1>(HIST("Prm/Neg/PtNum_310")), "f");
-    static_cast<TEfficiency*>(listEfficiencyMC->At(3))->SetTotalHistogram(*registry.get<TH1>(HIST("Prm/Neg/PtDen_310")), "f");
-
-    static_cast<TEfficiency*>(listEfficiencyMC->At(4))->SetPassedHistogram(*registry.get<TH1>(HIST("PrmRap/Pos/PtNum_310")), "f");
-    static_cast<TEfficiency*>(listEfficiencyMC->At(4))->SetTotalHistogram(*registry.get<TH1>(HIST("PrmRap/Pos/PtDen_310")), "f");
-    static_cast<TEfficiency*>(listEfficiencyMC->At(5))->SetPassedHistogram(*registry.get<TH1>(HIST("PrmRap/Neg/PtNum_310")), "f");
-    static_cast<TEfficiency*>(listEfficiencyMC->At(5))->SetTotalHistogram(*registry.get<TH1>(HIST("PrmRap/Neg/PtDen_310")), "f");
+    for (int i = 0; i < nSpecies; i++) {
+      static_cast<TEfficiency*>(listEfficiencyMC->At(i * nSpecies))->SetPassedHistogram(*histograms[i][kHistoPtNum].get(), "f");
+      static_cast<TEfficiency*>(listEfficiencyMC->At(i * nSpecies))->SetTotalHistogram(*histograms[i][kHistoPtDen].get(), "f");
+      static_cast<TEfficiency*>(listEfficiencyMC->At(i * nSpecies + 1))->SetPassedHistogram(*histogramsPrm[i][kHistoPtNum].get(), "f");
+      static_cast<TEfficiency*>(listEfficiencyMC->At(i * nSpecies + 1))->SetTotalHistogram(*histogramsPrm[i][kHistoPtDen].get(), "f");
+      static_cast<TEfficiency*>(listEfficiencyMC->At(i * nSpecies + 2))->SetPassedHistogram(*histogramsPrmRap[i][kHistoPtNum].get(), "f");
+      static_cast<TEfficiency*>(listEfficiencyMC->At(i * nSpecies + 2))->SetTotalHistogram(*histogramsPrmRap[i][kHistoPtDen].get(), "f");
+    }
   }
 };