diff --git a/PWGLF/Tasks/Nuspex/spectraTOF.cxx b/PWGLF/Tasks/Nuspex/spectraTOF.cxx index aed8e70c050..4084e53d73d 100644 --- a/PWGLF/Tasks/Nuspex/spectraTOF.cxx +++ b/PWGLF/Tasks/Nuspex/spectraTOF.cxx @@ -273,7 +273,7 @@ struct tofSpectra { histos.add("track/TPC/tpcCrossedRowsOverFindableCls", "crossed TPC rows over findable clusters;crossed rows / findable clusters TPC", kTH2D, {{60, 0.7, 1.3}, chargeAxis}); histos.add("track/TPC/tpcChi2NCl", "chi2 per cluster in TPC;chi2 / cluster TPC", kTH2D, {{100, 0, 10}, chargeAxis}); histos.add("Vertex/histGenVtxMC", "MC generated vertex z position", HistType::kTH1F, {{400, -40., +40., "z position (cm)"}}); - histos.add("Centrality/ImpParm", "Centrality", HistType::kTH1F, {{binsPercentile, "Centrality FT0C"}}); + histos.add("Centrality/ImpParm", "Centrality", HistType::kTH1F, {{binsPercentile, "Centrality FT0C"}}); histos.addClone("track/ITS/itsNCls", "track/selected/ITS/itsNCls"); histos.addClone("track/ITS/itsChi2NCl", "track/selected/ITS/itsChi2NCl"); @@ -338,7 +338,7 @@ struct tofSpectra { histos.add("MC/fake/neg", "Fake negative tracks", kTH1D, {ptAxis}); histos.add("MC/no_collision/pos", "No collision pos track", kTH1D, {ptAxis}); histos.add("MC/no_collision/neg", "No collision neg track", kTH1D, {ptAxis}); - histos.add("MC/test/pi/pos/prm/pt/den", "generated MC #pi^{+}", kTHnSparseD, {ptAxis, impParamAxis}); + histos.add("MC/test/pi/pos/prm/pt/den", "generated MC #pi^{+}", kTHnSparseD, {ptAxis, impParamAxis}); histos.add("MC/test/pi/neg/prm/pt/den", "generated MC #pi^{-}", kTHnSparseD, {ptAxis, impParamAxis}); histos.add("MC/test/ka/pos/prm/pt/den", "generated MC K^{+}", kTHnSparseD, {ptAxis, impParamAxis}); histos.add("MC/test/ka/neg/prm/pt/den", "generated MC K^{-}", kTHnSparseD, {ptAxis, impParamAxis}); @@ -565,8 +565,8 @@ struct tofSpectra { histos.add(hdcazprmgoodevs[i].data(), pTCharge[i], kTH2D, {ptAxis, dcaZAxis}); } } - //histos.add("Data/pos/pt/its_tpc_trd_tof", "pos ITS-TPC-TRD-TOF", kTH1D, {ptAxis}); - // Mismatched info + // histos.add("Data/pos/pt/its_tpc_trd_tof", "pos ITS-TPC-TRD-TOF", kTH1D, {ptAxis}); + // Mismatched info histos.add(hpt_mism_its_prm[i].data(), pTCharge[i], kTH1D, {ptAxis}); histos.add(hpt_mism_tpc_prm[i].data(), pTCharge[i], kTH1D, {ptAxis}); histos.add(hpt_mism_trd_prm[i].data(), pTCharge[i], kTH1D, {ptAxis}); @@ -1357,11 +1357,11 @@ struct tofSpectra { //************************************RD************************************************** auto mcCollision = collision.mcCollision(); - auto multiplicity = 0. ; + auto multiplicity = 0.; if (includeCentralityMC) { - multiplicity = mcCollision.impactParameter(); + multiplicity = mcCollision.impactParameter(); } else { - multiplicity = collision.centFT0C(); + multiplicity = collision.centFT0C(); } //************************************RD************************************************** @@ -1488,43 +1488,43 @@ struct tofSpectra { template void process_MC_gen(const McCollisionType& mcCollision, const McParticlesType& mcParticles) { - histos.fill(HIST("Vertex/histGenVtxMC"), mcCollision.posZ()); + histos.fill(HIST("Vertex/histGenVtxMC"), mcCollision.posZ()); histos.fill(HIST("Centrality/ImpParm"), mcCollision.impactParameter()); -const float multiplicity = mcCollision.impactParameter(); + const float multiplicity = mcCollision.impactParameter(); for (const auto& mcParticleGen : mcParticles) { if (!mcParticleGen.isPhysicalPrimary()) continue; int pdgCode = mcParticleGen.pdgCode(); - float pt = mcParticleGen.pt(); - float absY = std::abs(mcParticleGen.y()); - // Apply rapidity cut - if (absY > trkselOptions.cfgCutY) { - continue; - } + float pt = mcParticleGen.pt(); + float absY = std::abs(mcParticleGen.y()); + // Apply rapidity cut + if (absY > trkselOptions.cfgCutY) { + continue; + } // Fill histograms based on particle type - switch (pdgCode) { - case 2212: - histos.fill(HIST("MC/test/pr/pos/prm/pt/den"), pt, multiplicity); - break; - case -2212: - histos.fill(HIST("MC/test/pr/neg/prm/pt/den"), pt, multiplicity); - break; - case 211: - histos.fill(HIST("MC/test/pi/pos/prm/pt/den"), pt, multiplicity); - break; - case -211: - histos.fill(HIST("MC/test/pi/neg/prm/pt/den"), pt, multiplicity); - break; - case 321: - histos.fill(HIST("MC/test/ka/pos/prm/pt/den"), pt, multiplicity); - break; - case -321: - histos.fill(HIST("MC/test/ka/neg/prm/pt/den"), pt, multiplicity); - break; - default: - break; - } + switch (pdgCode) { + case 2212: + histos.fill(HIST("MC/test/pr/pos/prm/pt/den"), pt, multiplicity); + break; + case -2212: + histos.fill(HIST("MC/test/pr/neg/prm/pt/den"), pt, multiplicity); + break; + case 211: + histos.fill(HIST("MC/test/pi/pos/prm/pt/den"), pt, multiplicity); + break; + case -211: + histos.fill(HIST("MC/test/pi/neg/prm/pt/den"), pt, multiplicity); + break; + case 321: + histos.fill(HIST("MC/test/ka/pos/prm/pt/den"), pt, multiplicity); + break; + case -321: + histos.fill(HIST("MC/test/ka/neg/prm/pt/den"), pt, multiplicity); + break; + default: + break; + } } } template @@ -1963,7 +1963,7 @@ const float multiplicity = mcCollision.impactParameter(); } } PROCESS_SWITCH(tofSpectra, processMC, "Process MC", false); -void processMCgen(aod::McCollision const& mcCollision, aod::McParticles const& mcParticles) + void processMCgen(aod::McCollision const& mcCollision, aod::McParticles const& mcParticles) { process_MC_gen(mcCollision, mcParticles); }