Add files via upload

Author: 15071832337

PWGLF/Tasks/Strangeness/lambdaJetpolarization.cxx

@@ -146,6 +146,8 @@ struct LambdaJetpolarization {
   Configurable<int> pdgkZeroShort{"pdgkZeroShort", 310, "pdgkZeroShort"};
   Configurable<int> pdgProton{"pdgProton", 2212, "pdgProton"};
   Configurable<int> pdgPion{"pdgPion", 211, "pdgPion"};
+  Configurable<int> pdgXiMinus{"pdgXiMinus", 3312, "pdgXiMinus"};
+  Configurable<int> pdgXiZero{"pdgXiZero", 3322, "pdgXiZero"};
 
   void init(InitContext const&)
   {
@@ -2348,7 +2350,7 @@ struct LambdaJetpolarization {
 
         if (std::abs(v0mcparticle.pdgCode()) == pdgLmbda && v0mcparticle.has_mothers()) {
           for (const auto& mcparticleMother0 : v0mcparticle.mothers_as<aod::McParticles>()) {
-            if (std::abs(mcparticleMother0.pdgCode()) == 3312 || std::abs(mcparticleMother0.pdgCode()) == 3322) {
+            if (std::abs(mcparticleMother0.pdgCode()) == pdgXiMinus || std::abs(mcparticleMother0.pdgCode()) == pdgXiZero) {
               ptMotherMC = mcparticleMother0.pt();
               pdgMother = mcparticleMother0.pdgCode();
             }
@@ -2371,10 +2373,10 @@ struct LambdaJetpolarization {
             }
 
           } else if (std::abs(v0.mLambda() - o2::constants::physics::MassLambda0) < 0.01) {
-            if (pdgMother == 3312) {
+            if (pdgMother == pdgXiMinus) {
               registryData.fill(HIST("MC/hFDVsPtLambdaVsMotherPt_DoubleCharged_MC"), v0.pt(), ptMotherMC, mcCollision.centFT0M());
             }
-            if (pdgMother == 3312 || pdgMother == 3322) {
+            if (pdgMother == pdgXiMinus || pdgMother == pdgXiZero) {
               registryData.fill(HIST("MC/hFDVsPtLambdaVsMotherPt_MCRatio_MC"), v0.pt(), ptMotherMC, mcCollision.centFT0M());
             }
           }
@@ -2426,7 +2428,7 @@ struct LambdaJetpolarization {
         float pdgMother = 0.;
         if (std::abs(v0mcparticle.pdgCode()) == pdgLmbda && v0mcparticle.has_mothers()) {
           for (const auto& mcparticleMother0 : v0mcparticle.mothers_as<aod::McParticles>()) {
-            if (std::abs(mcparticleMother0.pdgCode()) == 3312 || std::abs(mcparticleMother0.pdgCode()) == 3322) {
+            if (std::abs(mcparticleMother0.pdgCode()) == pdgXiMinus || std::abs(mcparticleMother0.pdgCode()) == pdgXiZero) {
               // ptMotherMC = mcparticleMother0.pt();
               pdgMother = mcparticleMother0.pdgCode();
             }
@@ -2489,10 +2491,10 @@ struct LambdaJetpolarization {
             }
 
           } else if (std::abs(v0.mLambda() - o2::constants::physics::MassLambda0) < 0.01) {
-            if (pdgMother == 3312) {
+            if (pdgMother == pdgXiMinus) {
               continue;
             }
-            if (pdgMother == 3312 || pdgMother == 3322) {
+            if (pdgMother == pdgXiMinus || pdgMother == pdgXiZero) {
               continue;
             }
           }