fixedErrors

Author: baasingh

PWGCF/EbyEFluctuations/Tasks/FactorialMomentsTask.cxx

@@ -135,7 +135,7 @@ struct FactorialMomentsTask {
     OutputObjHandlingPolicy::AnalysisObject,
     true};
   static const Int_t nBins = 52;
-  constexpr double kMinCharge = 1e-6;
+  double kMinCharge = 1e-6;
   static const Int_t nfqOrder = 6;
   Int_t countSamples = 0;
   Int_t testc1 = 0, testc2 = 0, testc3 = 0;
@@ -195,11 +195,11 @@ struct FactorialMomentsTask {
       for (Int_t iM = 0; iM < nBins; ++iM) {
         auto mHistsR = std::get<std::shared_ptr<TH2>>(histos.add(Form("bin%i/Reset/mEtaPhi%i", iPt + 1, iM), Form("#eta#phi_%i for bin %.2f-%.2f;#eta;#phi", iM, ptCuts.value[2 * iPt], ptCuts.value[2 * iPt + 1]), HistType::kTH2F, {{binningM[iM], -0.8, 0.8}, {binningM[iM], 0, 2 * TMath::Pi()}}));
         mHistArrReset.push_back(mHistsR);
-        for (Int_t iq = 0; iq nfqOrder; ++iq) {
+        for (Int_t iq = 0; iq < nfqOrder; ++iq) {
           tmpFqErr[iq][iPt][iM] = new TH1D(Form("tmpFqErr%i%i%i", iq, iPt, iM), Form("tmpFqErr%i%i%i", iq, iPt, iM), 100, 0, 10);
         }
       }
-      for (Int_t i = 0; i nfqOrder; ++i) {
+      for (Int_t i = 0; i < nfqOrder; ++i) {
         auto mHistFq = std::get<std::shared_ptr<TH1>>(histos.add(Form("mFinalFq%i_bin%i", i + 2, iPt + 1), Form("Final F_%i for bin %.2f-%.2f;M", i + 2, ptCuts.value[2 * iPt], ptCuts.value[2 * iPt + 1]), HistType::kTH1F, {{nBins, -0.5, nBins - 0.5}}));
         mFqBinFinal.push_back(mHistFq);
         auto mHistAv = std::get<std::shared_ptr<TH1>>(histos.add(Form("mFinalAvBin%i_bin%i", i + 2, iPt + 1), Form("Final AvBin_%i for bin %.2f-%.2f;M", i + 2, ptCuts.value[2 * iPt], ptCuts.value[2 * iPt + 1]), HistType::kTH1F, {{nBins, -0.5, nBins - 0.5}}));
@@ -255,7 +255,7 @@ struct FactorialMomentsTask {
             double binconVal = 0;
             binconVal = hist[iPt * nBins + iM]->GetBinContent(iEta, iPhi);
             binContent += binconVal;
-            for (Int_t iq = 0; iq nfqOrder; ++iq) {
+            for (Int_t iq = 0; iq < nfqOrder; ++iq) {
               Double_t fqBin = 0;
               if (binconVal >= iq + 2) {
                 fqBin = TMath::Factorial(binconVal) / (TMath::Factorial(binconVal - (iq + 2)));
@@ -268,7 +268,7 @@ struct FactorialMomentsTask {
           }
         }
         binConEvent[iPt][iM] = binContent / (TMath::Power(binningM[iM], 2));
-        for (Int_t iq = 0; iq nfqOrder; ++iq) {
+        for (Int_t iq = 0; iq < nfqOrder; ++iq) {
           if (sumfqBin[iq] > 0) {
             fqEvent[iq][iPt][iM] = sumfqBin[iq] / (TMath::Power(binningM[iM], 2));
             fqEventSampled[iq][iPt][iM] += fqEvent[iq][iPt][iM];
@@ -321,13 +321,13 @@ struct FactorialMomentsTask {
     histos.fill(HIST("mCentFV0A"), coll.centFV0A());
     histos.fill(HIST("mCentFT0A"), coll.centFT0A());
     histos.fill(HIST("mCentFT0C"), coll.centFT0C());
-    for (Int_t const& h : mHistArrReset) {
+    for (auto const& h : mHistArrReset) {
       h->Reset();
     }
     countTracks = {0, 0, 0, 0, 0};
     fqEvent = {{{{{0, 0, 0, 0, 0, 0}}}}};
     binConEvent = {{{0, 0, 0, 0, 0}}};
-    for (Int_t const& track : tracks) {
+    for (auto const& track : tracks) {
       if (track.hasTPC())
       // if (track.hasITS())
       // if (track.isGlobalTrack())
@@ -386,15 +386,15 @@ struct FactorialMomentsTask {
     histos.fill(HIST("mVertexY"), coll.posY());
     histos.fill(HIST("mVertexZ"), coll.posZ());
     histos.fill(HIST("mCentFT0M"), coll.centRun2V0M());
-    for (Int_t const& h : mHistArrReset) {
+    for (auto const& h : mHistArrReset) {
       h->Reset();
     }
     countTracks = {0, 0, 0, 0, 0};
     fqEvent = {{{{{0, 0, 0, 0, 0, 0}}}}};
     binConEvent = {{{0, 0, 0, 0, 0}}};
-    for (Int_t const& track : tracks) {
+    for (auto const& track : tracks) {
       double recoCharge = (track.sign() != 0) ? track.sign() : 0.;
-      if (std::abs(track.eta()) centralEta && track.isGlobalTrack() && std::abs(recoCharge) >= kMinCharge) {
+      if (std::abs(track.eta()) < centralEta && track.isGlobalTrack() && std::abs(recoCharge) >= kMinCharge) {
         histos.fill(HIST("mCollID"), track.collisionId());
         histos.fill(HIST("mNFindableClsTPC"), track.tpcNClsFindable());
         histos.fill(HIST("mNClsTPC"), track.tpcNClsFound());
@@ -416,7 +416,7 @@ struct FactorialMomentsTask {
 
     auto mccollision = coll.mcCollision_as<aod::McCollisions>();
     auto mcParts = mcParticles.sliceBy(perMcCollision, coll.mcCollision().globalIndex());
-    for (Int_t const& mc : mcParts) {
+    for (auto const& mc : mcParts) {
       int pdgCode = mc.pdgCode();
       auto pdgInfo = pdg->GetParticle(pdgCode);
       if (!pdgInfo) {
@@ -425,7 +425,7 @@ struct FactorialMomentsTask {
       double charge = pdgInfo->Charge();
       double physCharge = charge / 3.0;
       histos.fill(HIST("mChargeBefore"), physCharge);
-      if (mc.isPhysicalPrimary() && std::abs(mc.eta()) centralEta && std::abs(physCharge) >= kMinCharge) {
+      if (mc.isPhysicalPrimary() && std::abs(mc.eta()) < centralEta && std::abs(physCharge) >= kMinCharge) {
         histos.fill(HIST("mChargeAfter"), physCharge);
         histos.fill(HIST("mEta"), mc.eta());
         histos.fill(HIST("mPt"), mc.pt());
@@ -461,13 +461,13 @@ struct FactorialMomentsTask {
     histos.fill(HIST("mVertexZ"), coll.posZ());
     histos.fill(HIST("mCentFT0M"), coll.centRun2V0M());
 
-    for (Int_t const& h : mHistArrReset) {
+    for (auto const& h : mHistArrReset) {
       h->Reset();
     }
     countTracks = {0, 0, 0, 0, 0};
     fqEvent = {{{{{0, 0, 0, 0, 0, 0}}}}};
     binConEvent = {{{0, 0, 0, 0, 0}}};
-    for (Int_t const& track : tracks) {
+    for (auto const& track : tracks) {
       if ((track.pt() < ptMin) || (!track.isGlobalTrack()) || (track.tpcNClsFindable() < minTPCCls)) {
         continue;
       }