Add files via upload

Author: 15071832337

PWGLF/Tasks/Strangeness/lambdaJetpolarization.cxx

@@ -127,7 +127,7 @@ struct LambdaJetpolarization {
   Configurable<float> v0TracketaMax{"v0TracketaMax", +0.8f, "eta max track"};
   Configurable<bool> requireTPC{"requireTPC", true, "require TPC hit"};
   Configurable<float> yMin{"yMin", -0.5f, "minimum y"};
-  Configurable<float> yMax{"yMax", +0.5f, "maximum y"};
+  Configurable<float> yMC{"yMax", +0.5f, "maximum y"};
   Configurable<float> v0rejLambda{"v0rejLambda", 0.00, "V0 rej Lambda"};
   Configurable<float> v0accLambda{"v0accLambda", 0.075, "V0 acc Lambda"};
   Configurable<bool> ifinitpasslambda{"ifinitpasslambda", 0, "ifinitpasslambda"};
@@ -2207,7 +2207,7 @@ struct LambdaJetpolarization {
       auto v0sThisCollision = V0s.sliceBy(perCol, collision.globalIndex());
       auto tracksThisCollision = tracks.sliceBy(perColtrack, collision.globalIndex());
       std::vector<fastjet::PseudoJet> fjParticles;
-      for (auto& track : tracksThisCollision) {
+      for (const auto& track : tracksThisCollision) {
         registryData.fill(HIST("MC/h_MCtrack_pt"), track.pt());
         registryData.fill(HIST("MC/h_MCtrack_eta"), track.eta());
         registryData.fill(HIST("MC/h_MCtrack_phi"), track.phi());
@@ -2299,7 +2299,7 @@ struct LambdaJetpolarization {
           continue;
         }
         auto v0mcparticle = v0.mcParticle();
-        if (std::abs(v0mcparticle.y()) > 0.5f) {
+        if (std::abs(v0mcparticle.y()) > yMC) {
           continue;
         }
         if (v0.v0Type() != v0TypeSelection) {
@@ -2529,7 +2529,7 @@ struct LambdaJetpolarizationMCcount {
     registry.get<TH1>(HIST("hEventSelection"))->GetXaxis()->SetBinLabel(3, "posZ cut");
   }
 
-  Configurable<float> rapidityMCcut{"rapidityMCcut", 0.5, "rapidity cut MC count"};
+  Configurable<float> rapidityMCcut{"rapidityMCcut", 0.5f, "rapidity cut MC count"};
   Configurable<bool> eventSel8selection{"eventSel8selection", true, "event selection MC count post sel8 cut"};
   Configurable<bool> eventPosZselection{"eventPosZselection", true, "event selection MC count post posZ cut"};
 
@@ -2557,7 +2557,7 @@ struct LambdaJetpolarizationMCcount {
     }
     registry.fill(HIST("hEventSelection"), 2.5);
 
-    for (auto& mcparticle : mcParticles) {
+    for (const auto& mcparticle : mcParticles) {
       if (TMath::Abs(mcparticle.y()) < rapidityMCcut) {
         if (mcparticle.isPhysicalPrimary()) {
           if (!mcparticle.has_daughters()) {
@@ -2566,8 +2566,8 @@ struct LambdaJetpolarizationMCcount {
           if (mcparticle.pdgCode() == 310) {
             registry.fill(HIST("hK0ShortCount"), 0.5);
             registry.fill(HIST("hK0ShortCount_PtDiff"), mcparticle.pt());
-            for (auto& mcparticleDaughter0 : mcparticle.daughters_as<aod::McParticles>()) {
-              for (auto& mcparticleDaughter1 : mcparticle.daughters_as<aod::McParticles>()) {
+            for (const auto& mcparticleDaughter0 : mcparticle.daughters_as<aod::McParticles>()) {
+              for (const auto& mcparticleDaughter1 : mcparticle.daughters_as<aod::McParticles>()) {
                 if (mcparticleDaughter0.pdgCode() == 211 && mcparticleDaughter1.pdgCode() == -211) {
                   registry.fill(HIST("hK0ShortCount"), 1.5);
                 }
@@ -2577,8 +2577,8 @@ struct LambdaJetpolarizationMCcount {
           if (mcparticle.pdgCode() == 3122) {
             registry.fill(HIST("hLambdaCount"), 0.5);
             registry.fill(HIST("hLambdaCount_PtDiff"), mcparticle.pt());
-            for (auto& mcparticleDaughter0 : mcparticle.daughters_as<aod::McParticles>()) {
-              for (auto& mcparticleDaughter1 : mcparticle.daughters_as<aod::McParticles>()) {
+            for (const auto& mcparticleDaughter0 : mcparticle.daughters_as<aod::McParticles>()) {
+              for (const auto& mcparticleDaughter1 : mcparticle.daughters_as<aod::McParticles>()) {
                 if (mcparticleDaughter0.pdgCode() == -211 && mcparticleDaughter1.pdgCode() == 2212) {
                   registry.fill(HIST("hLambdaCount"), 1.5);
                 }
@@ -2588,8 +2588,8 @@ struct LambdaJetpolarizationMCcount {
           if (mcparticle.pdgCode() == -3122) {
             registry.fill(HIST("hAntiLambdaCount"), 0.5);
             registry.fill(HIST("hAntiLambdaCount_PtDiff"), mcparticle.pt());
-            for (auto& mcparticleDaughter0 : mcparticle.daughters_as<aod::McParticles>()) {
-              for (auto& mcparticleDaughter1 : mcparticle.daughters_as<aod::McParticles>()) {
+            for (const auto& mcparticleDaughter0 : mcparticle.daughters_as<aod::McParticles>()) {
+              for (const auto& mcparticleDaughter1 : mcparticle.daughters_as<aod::McParticles>()) {
                 if (mcparticleDaughter0.pdgCode() == 211 && mcparticleDaughter1.pdgCode() == -2212) {
                   registry.fill(HIST("hAntiLambdaCount"), 1.5);
                 }