PWGHF: Bug fix to the creator that was pushed (#225)

* Bug fix to the creator that was pushed

* Format changes

Author: pchrista

PWGHF/TableProducer/HFCandidateCreatorLb.cxx

@@ -105,6 +105,9 @@ struct HFCandidateCreatorLb {
       if (lcCand.isSelLcpiKp() >= d_selectionFlagLc) {
         hMassLcToPKPi->Fill(InvMassLcpiKp(lcCand), lcCand.pt());
       }
+      hPtLc->Fill(lcCand.pt());
+      hCPALc->Fill(lcCand.cpa());
+
       auto track0 = lcCand.index0_as<aod::BigTracks>();
       auto track1 = lcCand.index1_as<aod::BigTracks>();
       auto track2 = lcCand.index2_as<aod::BigTracks>();
@@ -138,12 +141,13 @@ struct HFCandidateCreatorLb {
         if (trackPion.pt() < ptPionMin) {
           continue;
         }
-        if (trackPion.sign() * charge < 0) {
+        if (trackPion.sign() > 0) {
           continue;
         }
         if (trackPion.globalIndex() == index0Lc || trackPion.globalIndex() == index1Lc || trackPion.globalIndex() == index2Lc) {
           continue;
         }
+        hPtPion->Fill(trackPion.pt());
         array<float, 3> pvecPion;
         auto trackParVarPi = getTrackParCov(trackPion);
 
@@ -168,6 +172,9 @@ struct HFCandidateCreatorLb {
         trackLc.propagateToDCA(primaryVertex, magneticField, &impactParameter0);
         trackParVarPi.propagateToDCA(primaryVertex, magneticField, &impactParameter1);
 
+        hCovSVXX->Fill(covMatrixPCA[0]);
+        hCovPVXX->Fill(covMatrixPV[0]);
+
         // get uncertainty of the decay length
         double phi, theta;
         getPointDirection(array{collision.posX(), collision.posY(), collision.posZ()}, secondaryVertexLb, phi, theta);

PWGHF/TableProducer/HFLbToLcPiCandidateSelector.cxx

@@ -80,8 +80,8 @@ struct HfLbToLcPiCandidateSelector {
       return false;
     }
 
-    //d0(Lc)xd0(pi)
-    if (hfCandLb.impactParameterProduct() > cuts->get(pTBin, "Imp. Par. Product")) {
+    //Lc+ pt
+    if (hfCandLc.pt() < cuts->get(pTBin, "pT Lc+")) {
       return false;
     }
 
@@ -103,16 +103,20 @@ struct HfLbToLcPiCandidateSelector {
     }
 
     //Lb CPA cut
-    //if (hfCandLb.cpa() < cuts->get(pTBin, "CPA")) {
-    //return false;
-    //}
+    if (hfCandLb.cpa() < cuts->get(pTBin, "CPA")) {
+      return false;
+    }
 
-    //d0 of Lc and pi
-    if ((std::abs(hfCandLb.impactParameter0()) > cuts->get(pTBin, "d0 Lc+")) ||
-        (std::abs(hfCandLb.impactParameter1()) > cuts->get(pTBin, "d0 Pi"))) {
+    //d0 of pi
+    if (std::abs(hfCandLb.impactParameter1()) < cuts->get(pTBin, "d0 Pi")) {
       return false;
     }
 
+    //d0 of Lc+
+    //if(std::abs(hfCandLb.impactParameter0()) > cuts->get(pTBin, "d0 Lc+")){
+    //return false;
+    //}
+
     return true;
   }
 

PWGHF/Tasks/taskLb.cxx

@@ -72,9 +72,8 @@ struct HfTaskLb {
   Filter filterSelectCandidates = (aod::hf_selcandidate_lb::isSelLbToLcPi >= selectionFlagLb);
 
   void process(soa::Join<aod::Collisions, aod::CentV0Ms>::iterator const& collision, soa::Filtered<soa::Join<aod::HfCandLb, aod::HFSelLbToLcPiCandidate>> const& candidates, soa::Join<aod::HfCandProng3, aod::HFSelLcCandidate>, aod::BigTracks)
-  //void process(soa::Filtered<soa::Join<aod::HfCandLb, aod::HFSelLbToLcPiCandidate>> const& candidates, soa::Join<aod::HfCandProng3, aod::HFSelLcCandidate>, aod::BigTracks)
   {
-    float centrality = collision.centV0M(); //FIXME: Add centrality properly;removed for the time being (taking forever to run)
+    float centrality = collision.centV0M();
     registry.fill(HIST("hCentrality"), centrality);
 
     for (auto& candidate : candidates) {