clean error + add histograms

Author: Valerio Di Bella

PWGHF/HFC/TableProducer/correlatorDplusDplusReduced.cxx

@@ -283,7 +283,7 @@ struct HfCorrelatorDplusDplusReduced {
     }
   }
 
-  void processData(aod::Collisions const& collisions, SelectedCandidates const& candidates, aod::Tracks const&)
+  void processData(aod::Collisions const& collisions, SelectedCandidates const& candidates, aod::Tracks const&, aod::BCsWithTimestamps const&)
   {
     static int lastRunNumber = -1;
     // reserve memory
@@ -316,15 +316,15 @@ struct HfCorrelatorDplusDplusReduced {
       fillEvent(collision);
       for (const auto& candidate : candidatesInThisCollision) {
         auto prong_candidate = candidate.prong1_as<aod::Tracks>();
-        auto candidate_sign = prong_candidate.sign();
+        auto candidate_sign = -prong_candidate.sign();
         fillCandidateTable<aod::Collisions>(candidate, rowCandidateFullEvents.lastIndex(), candidate_sign);
       }
     }
   }
   PROCESS_SWITCH(HfCorrelatorDplusDplusReduced, processData, "Process data per collision", false);
 
   void processMcRec(aod::Collisions const& collisions,
-                    SelectedCandidatesMc const& candidates)
+                    SelectedCandidatesMc const& candidates, aod::Tracks const&)
   {
     // reserve memory
     rowCandidateFullEvents.reserve(collisions.size());
@@ -343,7 +343,7 @@ struct HfCorrelatorDplusDplusReduced {
       fillEvent(collision);
       for (const auto& candidate : candidatesInThisCollision) {
         auto prong_candidate = candidate.prong1_as<aod::Tracks>();
-        auto candidate_sign = prong_candidate.sign();
+        auto candidate_sign = -prong_candidate.sign();
         fillCandidateTable<aod::Collisions, true>(candidate, rowCandidateFullEvents.lastIndex(), candidate_sign);
       }
     }

PWGHF/HFC/Tasks/taskCorrelationDplusDplusReduced.cxx

@@ -55,7 +55,9 @@ struct HfTaskCorrelationDplusDplusReduced {
   void init(InitContext const&)
   {
     registry.add("hMassDplus", "D+ candidates;inv. mass (#pi#pi K) (GeV/#it{c}^{2}))", {HistType::kTH1F, {{120, 1.5848, 2.1848}}});
+    registry.add("hMassDminus", "D- candidates;inv. mass (#pi#pi K) (GeV/#it{c}^{2}))", {HistType::kTH1F, {{120, 1.5848, 2.1848}}});
     registry.add("hMassDplusMatched", "D+ matched candidates;inv. mass (#pi#pi K) (GeV/#it{c}^{2}))", {HistType::kTH1F, {{120, 1.5848, 2.1848}}});
+    registry.add("hMassDminusMatched", "D- matched candidates;inv. mass (#pi#pi K) (GeV/#it{c}^{2}))", {HistType::kTH1F, {{120, 1.5848, 2.1848}}});
     registry.add("hMassDplusminusPair", "D plus-minus pair candidates;inv. mass (#pi K) (GeV/#it{c}^{2});inv. mass (#pi K) (GeV/#it{c}^{2})", {HistType::kTH2F, {{120, 1.5848, 2.1848}, {120, 1.5848, 2.1848}}});
     registry.add("hMassDplusPair", "D plus pair candidates;inv. mass (#pi K) (GeV/#it{c}^{2});inv. mass (#pi K) (GeV/#it{c}^{2})", {HistType::kTH2F, {{120, 1.5848, 2.1848}, {120, 1.5848, 2.1848}}});
     registry.add("hMassDminusPair", "D minus pair candidates;inv. mass (#pi K) (GeV/#it{c}^{2});inv. mass (#pi K) (GeV/#it{c}^{2})", {HistType::kTH2F, {{120, 1.5848, 2.1848}, {120, 1.5848, 2.1848}}});
@@ -72,7 +74,11 @@ struct HfTaskCorrelationDplusDplusReduced {
       auto sign1 = 1;
       if (cand1.pt() < 0) {
         sign1 = -1;
+        registry.fill(HIST("hMassDminus"), mass1);
+      } else {
+        registry.fill(HIST("hMassDplus"), mass1);
       }
+
       for (auto cand2 = cand1 + 1; cand2 != localCandidates.end(); ++cand2) {
         auto mass2 = cand2.m();
         auto sign2 = 1;
@@ -100,9 +106,15 @@ struct HfTaskCorrelationDplusDplusReduced {
 
     for (const auto& cand1 : localCandidates) {
       auto mass1 = cand1.m();
-      registry.fill(HIST("hMassDplus"), mass1);
-      if (std::abs(cand1.flagMcMatchRec()) == hf_decay::hf_cand_3prong::DecayChannelMain::DplusToPiKPi)
-        registry.fill(HIST("hMassDplusMatched"), mass1);
+      if (cand1.pt() < 0) {
+        registry.fill(HIST("hMassDminus"), mass1);
+        if (std::abs(cand1.flagMcMatchRec()) == hf_decay::hf_cand_3prong::DecayChannelMain::DplusToPiKPi)
+          registry.fill(HIST("hMassDminusMatched"), mass1);
+      } else {
+        registry.fill(HIST("hMassDplus"), mass1);
+        if (std::abs(cand1.flagMcMatchRec()) == hf_decay::hf_cand_3prong::DecayChannelMain::DplusToPiKPi)
+          registry.fill(HIST("hMassDplusMatched"), mass1);
+      }
     }
   }
   PROCESS_SWITCH(HfTaskCorrelationDplusDplusReduced, processLocalDataMcRec, "Process local MC data", false);