Changes to the producer

Author: zchochul

PWGCF/FemtoWorld/DataModel/FemtoWorldDerived.h

@@ -74,6 +74,7 @@ DECLARE_SOA_INDEX_COLUMN(FemtoWorldCollision, femtoWorldCollision);
 DECLARE_SOA_COLUMN(Pt, pt, float);                    //! p_T (GeV/c)
 DECLARE_SOA_COLUMN(Eta, eta, float);                  //! Eta
 DECLARE_SOA_COLUMN(Phi, phi, float);                  //! Phi
+DECLARE_SOA_COLUMN(Mass, mass, int8_t);                                                  //! Mass of the particle
 DECLARE_SOA_COLUMN(PartType, partType, uint8_t);      //! Type of the particle, according to femtoworldparticle::ParticleType
 DECLARE_SOA_COLUMN(Cut, cut, cutContainerType);       //! Bit-wise container for the different selection criteria
 DECLARE_SOA_COLUMN(PIDCut, pidcut, cutContainerType); //! Bit-wise container for the different PID selection criteria \todo since bit-masking cannot be done yet with filters we use a second field for the PID
@@ -134,6 +135,7 @@ DECLARE_SOA_TABLE(FemtoWorldParticles, "AOD", "FEMTOWORLDPARTS",
                   femtoworldparticle::Eta,
                   femtoworldparticle::Phi,
                   track::P,
+                  femtoworldparticle::Mass,
                   femtoworldparticle::PartType,
                   femtoworldparticle::Cut,
                   femtoworldparticle::PIDCut,

PWGCF/FemtoWorld/TableProducer/femtoWorldProducerTask.cxx

@@ -797,15 +797,7 @@ struct femtoWorldProducerTask {
             continue;
           }
         }
-        /*else if (p1.pt() == NAN){
-         continue;
-       }*/
-        // LOGF(info, "p1.pt() = %f \n p1.p() = %f \n", p1.pt(), p1.p());
-        // LOGF(info, "p2.pt() = %f \n p2.p() = %f \n", p2.pt(), p2.p());
-        //  float phiPx = p1.px() + p2.px();
-        //  float phiPy = p1.py() + p2.py();
-        //  float phiPz = p1.pz() + p2.pz();
-
+       
         TLorentzVector part1Vec;
         TLorentzVector part2Vec;
         float mMassOne = TDatabasePDG::Instance()->GetParticle(ConfPDGCodePartOne)->Mass();
@@ -821,6 +813,7 @@ struct femtoWorldProducerTask {
         float phiPhi = sumVec.Phi();
         float phiPt = sumVec.Pt();
         float phiP = sumVec.P();
+        float phiM = sumVec.M();
 
         // LOGF(info, "FIRST DAUGHTER: \n pT = %f \n eta = %f \n phi = %f", p1.pt(), p1.eta(), p1.phi());
 
@@ -839,6 +832,7 @@ struct femtoWorldProducerTask {
                       p1.eta(),
                       p1.phi(),
                       p1.p(),
+                      mMassOne,
                       aod::femtoworldparticle::ParticleType::kPhiChild,
                       cutContainerV0.at(femtoWorldV0Selection::V0ContainerPosition::kPosCuts),
                       cutContainerV0.at(femtoWorldV0Selection::V0ContainerPosition::kPosPID),
@@ -890,6 +884,7 @@ struct femtoWorldProducerTask {
                       p2.eta(),
                       p2.phi(),
                       p2.p(),
+                      mMassTwo,
                       aod::femtoworldparticle::ParticleType::kPhiChild,
                       cutContainerV0.at(femtoWorldV0Selection::V0ContainerPosition::kNegCuts),
                       cutContainerV0.at(femtoWorldV0Selection::V0ContainerPosition::kNegPID),
@@ -937,6 +932,7 @@ struct femtoWorldProducerTask {
                       phiEta,
                       phiPhi,
                       phiP,
+                      phiM,
                       aod::femtoworldparticle::ParticleType::kPhi,
                       cutContainerV0.at(femtoWorldV0Selection::V0ContainerPosition::kV0),
                       0,

PWGCF/FemtoWorld/Tasks/femtoWorldPairTaskTrackPhi.cxx

@@ -134,17 +134,18 @@ struct femtoWorldPairTaskTrackPhi {
     eventHisto.fillQA(col);
     /// Histogramming same event
     for (auto& part : groupPartsOne) {
+      /*
       if (part.p() > cfgCutTable->get("PartOne", "MaxP") || part.pt() > cfgCutTable->get("PartOne", "MaxPt")) {
         continue;
       }
       if (!isFullPIDSelected(part.pidcut(), part.p(), cfgCutTable->get("PartOne", "PIDthr"), vPIDPartOne, cfgNspecies, kNsigma, cfgCutTable->get("PartOne", "nSigmaTPC"), cfgCutTable->get("PartOne", "nSigmaTPCTOF"))) {
         continue;
-      }
-      LOGF(info, "TEST ================================================================= \n");
+      }*/
+      LOGF(info, "TEST ==============================TRACKONEFILL=================================== \n");
       trackHistoPartOne.fillQA(part);
     }
     for (auto& part : groupPartsTwo) {
-
+      //LOGF(info, "TEST ==============================TRACKPHIFILL=================================== \n");
       trackHistoPartTwo.fillQA(part);
     }
     /// Now build the combinations