Skip to content

Commit

Permalink
Merge pull request AliceO2Group#68 from alibuild/alibot-cleanup-6380
Browse files Browse the repository at this point in the history 
Please consider the following formatting changes to AliceO2Group#6380
  • Loading branch information
@RD0407
RD0407 committed Jun 4, 2024
2 parents aeedd1d + ee047a5 commit 453a21f
Showing 1 changed file with 38 additions and 38 deletions.
76 changes: 38 additions & 38 deletions PWGLF/Tasks/Nuspex/spectraTOF.cxx
View file
Original file line number Diff line number Diff line change
Expand Up @@ -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");
Expand Down Expand Up @@ -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});
Expand Down Expand Up @@ -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});
Expand Down Expand Up @@ -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**************************************************

Expand Down Expand Up @@ -1488,43 +1488,43 @@ struct tofSpectra {
template <typename McCollisionType, typename McParticlesType>
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 <std::size_t i, typename ParticleType>
Expand Down Expand Up @@ -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);
}
Expand Down

0 comments on commit 453a21f

Please sign in to comment.