@@ -78,7 +78,7 @@ std::array<std::shared_ptr<TH2>, NpCharge> hDecayLengthMCNotHF; // Decay Length
7878
7979std::array<std::shared_ptr<TH2>, NpCharge> hPtNumTOFMatchWithPIDSignalPrm; // Pt distribution of particles with a hit in the TOF and a compatible signal
8080
81- std::array<std::array<std::shared_ptr<TH3>, NpCharge>, 3 > hMCpdg_nsigmaTPC; // 2D array of nsigmaTPC histograms [Selection: pi,K,p][True PDG: 18 species]
81+ std::array<std::array<std::shared_ptr<TH3>, NpCharge>, 3 > hMCpdg_nsigmaTPC; // 2D array of nsigmaTPC histograms [Selection: pi,K,p][True PDG: 18 species]
8282
8383// Spectra task
8484struct tofSpectra {
@@ -587,12 +587,12 @@ struct tofSpectra {
587587 histos.add (" MC/MultiplicityMCINELgt0" , " MC multiplicity" , kTH1D , {multAxis});
588588 histos.add (" MC/MultiplicityMCINELgt1" , " MC multiplicity" , kTH1D , {multAxis});
589589 }
590-
591- if (doprocessTrackMCLabels){
592- for (int par = 2 ; par <= 4 ; par++){
593- for (int i = 0 ; i < NpCharge; i++){
594- hMCpdg_nsigmaTPC[par- 2 ][i] = histos.add <TH3>(Form (" test_mclabels/nsigmatpc/%s/%s/pdg_%i" , (i < Np) ? " pos" : " neg" , pN[par], PDGs[i% Np]), Form (" True %s (%i) in %s selection" , pTCharge[i], PDGs[i], (i < Np) ? pTCharge[par] : pTCharge[par+ Np]), kTH3D , {ptAxis, nsigmaTPCAxisOccupancy, multAxis});
595- }
590+
591+ if (doprocessTrackMCLabels) {
592+ for (int par = 2 ; par <= 4 ; par++) {
593+ for (int i = 0 ; i < NpCharge; i++) {
594+ hMCpdg_nsigmaTPC[par - 2 ][i] = histos.add <TH3>(Form (" test_mclabels/nsigmatpc/%s/%s/pdg_%i" , (i < Np) ? " pos" : " neg" , pN[par], PDGs[i % Np]), Form (" True %s (%i) in %s selection" , pTCharge[i], PDGs[i], (i < Np) ? pTCharge[par] : pTCharge[par + Np]), kTH3D , {ptAxis, nsigmaTPCAxisOccupancy, multAxis});
595+ }
596596 }
597597 }
598598
@@ -2795,15 +2795,15 @@ struct tofSpectra {
27952795 }
27962796 }
27972797 PROCESS_SWITCH (tofSpectra, processMCgen_RecoEvs, " process generated MC (reconstructed events)" , false );
2798-
2798+
27992799 void processTrackMCLabels (CollisionCandidates::iterator const & collisions,
2800- soa::Join<TrackCandidates,
2801- aod::pidTPCFullPi, aod::pidTPCFullKa, aod::pidTPCFullPr,
2802- aod::pidTOFFullPi, aod::pidTOFFullKa, aod::pidTOFFullPr> const & tracks,
2803- aod::McTrackLabels const & mcTrackLabels, aod::McParticles const & mcParticles)
2800+ soa::Join<TrackCandidates,
2801+ aod::pidTPCFullPi, aod::pidTPCFullKa, aod::pidTPCFullPr,
2802+ aod::pidTOFFullPi, aod::pidTOFFullKa, aod::pidTOFFullPr> const & tracks,
2803+ aod::McTrackLabels const & mcTrackLabels, aod::McParticles const & mcParticles)
28042804 {
28052805 const float multiplicity = getMultiplicity (collisions);
2806-
2806+
28072807 for (const auto & track : tracks) {
28082808 if (!track.has_collision ()) {
28092809 continue ;
@@ -2818,17 +2818,17 @@ struct tofSpectra {
28182818 const auto & mcLabel = mcTrackLabels.iteratorAt (track.globalIndex ());
28192819 const auto & mcParticle = mcParticles.iteratorAt (mcLabel.mcParticleId ());
28202820 int pdgCode = mcParticle.pdgCode ();
2821-
2822- static_for<2 ,4 >([&](auto par){
2821+
2822+ static_for<2 , 4 >([&](auto par) {
28232823 const auto & nsigmaTPCpar = o2::aod::pidutils::tpcNSigma<par>(track);
28242824 bool isTPCpar = std::abs (nsigmaTPCpar) < trkselOptions.cfgCutNsigma ;
28252825 // Precompute rapidity values to avoid redundant calculations
28262826 double rapiditypar = std::abs (track.rapidity (PID::getMass (par)));
2827- // TPC Selection and histogram filling
2828- if (isTPCpar && rapiditypar <= trkselOptions.cfgCutY ){
2829- static_for<0 , 17 >([&](auto i){
2830- if (pdgCode == PDGs[i]){
2831- hMCpdg_nsigmaTPC[par- 2 ][i]->Fill (track.pt (), nsigmaTPCpar, multiplicity);
2827+ // TPC Selection and histogram filling
2828+ if (isTPCpar && rapiditypar <= trkselOptions.cfgCutY ) {
2829+ static_for<0 , 17 >([&](auto i) {
2830+ if (pdgCode == PDGs[i]) {
2831+ hMCpdg_nsigmaTPC[par - 2 ][i]->Fill (track.pt (), nsigmaTPCpar, multiplicity);
28322832 }
28332833 });
28342834 }
0 commit comments