@@ -66,8 +66,8 @@ struct ThreePartCorr {
6666 using MyFilteredMCParticles = soa::Filtered<aod::McParticles>;
6767 using MyFilteredMCRecCollision = soa::Filtered<soa::Join<aod::Collisions, aod::EvSels, aod::McCollisionLabels>>::iterator;
6868 using MyFilteredMCTracks = soa::Filtered<soa::Join<aod::Tracks, aod::TracksExtra, aod::TrackSelection, aod::McTrackLabels,
69- aod::pidTPCPi, aod::pidTPCKa, aod::pidTPCPr,
70- aod::pidTOFPi, aod::pidTOFKa, aod::pidTOFPr, aod::pidTOFbeta>>;
69+ aod::pidTPCPi, aod::pidTPCKa, aod::pidTPCPr,
70+ aod::pidTOFPi, aod::pidTOFKa, aod::pidTOFPr, aod::pidTOFbeta>>;
7171
7272 // Mixed-events binning policy
7373 SliceCache cache;
@@ -320,8 +320,8 @@ struct ThreePartCorr {
320320 for (const auto & particle : particles) {
321321 if (particle.isPhysicalPrimary ()) {
322322
323- // Efficiency - Generated
324- MCRegistry.fill (HIST (" hGenerated" ), particle.pt ());
323+ // Efficiency - Generated
324+ MCRegistry.fill (HIST (" hGenerated" ), particle.pt ());
325325 if (particle.pdgCode () == kPiPlus ) { // Pos pions
326326 MCRegistry.fill (HIST (" hGenPionP" ), particle.pt ());
327327 } else if (particle.pdgCode () == kPiMinus ) { // Neg pions
@@ -356,8 +356,8 @@ struct ThreePartCorr {
356356 auto particle = track.mcParticle ();
357357 if (particle.isPhysicalPrimary ()) {
358358
359- // Efficiency - Reconstructed
360- MCRegistry.fill (HIST (" hReconstructed" ), track.pt ());
359+ // Efficiency - Reconstructed
360+ MCRegistry.fill (HIST (" hReconstructed" ), track.pt ());
361361 if (particle.pdgCode () == kPiPlus ) { // Pos pions
362362 MCRegistry.fill (HIST (" hRecPionP" ), track.pt ());
363363 } else if (particle.pdgCode () == kPiMinus ) { // Neg pions
@@ -372,46 +372,46 @@ struct ThreePartCorr {
372372 MCRegistry.fill (HIST (" hRecProtonN" ), track.pt ());
373373 }
374374
375- // Purity
376- A_PID = TrackPID (track);
377- if (A_PID[1 ] < 4.0 ) {
378-
379- if (track.sign () > 0 ) { // Positive tracks
380- if (A_PID[0 ] == 0.0 ) { // Pions
381- MCRegistry.fill (HIST (" hSelectPionP" ), track.pt ());
382- if (particle.pdgCode () == kPiPlus ) {
383- MCRegistry.fill (HIST (" hTrueSelectPionP" ), track.pt ());
384- }
385- } else if (A_PID[0 ] == 1.0 ) { // Kaons
386- MCRegistry.fill (HIST (" hSelectKaonP" ), track.pt ());
387- if (particle.pdgCode () == kKPlus ) {
388- MCRegistry.fill (HIST (" hTrueSelectKaonP" ), track.pt ());
389- }
390- } else if (A_PID[0 ] == 2.0 ) { // Protons
391- MCRegistry.fill (HIST (" hSelectProtonP" ), track.pt ());
392- if (particle.pdgCode () == kProton ) {
393- MCRegistry.fill (HIST (" hTrueSelectProtonP" ), track.pt ());
394- }
395- }
396- } else if (track.sign () < 0 ) { // Negative tracks
397- if (A_PID[0 ] == 0.0 ) { // Pions
398- MCRegistry.fill (HIST (" hSelectPionN" ), track.pt ());
399- if (particle.pdgCode () == kPiMinus ) {
400- MCRegistry.fill (HIST (" hTrueSelectPionN" ), track.pt ());
401- }
402- } else if (A_PID[0 ] == 1.0 ) { // Kaons
403- MCRegistry.fill (HIST (" hSelectKaonN" ), track.pt ());
404- if (particle.pdgCode () == kKMinus ) {
405- MCRegistry.fill (HIST (" hTrueSelectKaonN" ), track.pt ());
406- }
407- } else if (A_PID[0 ] == 2.0 ) { // Protons
408- MCRegistry.fill (HIST (" hSelectProtonN" ), track.pt ());
409- if (particle.pdgCode () == kProtonBar ) {
410- MCRegistry.fill (HIST (" hTrueSelectProtonN" ), track.pt ());
411- }
412- }
413- }
414- }
375+ // Purity
376+ A_PID = TrackPID (track);
377+ if (A_PID[1 ] < 4.0 ) {
378+
379+ if (track.sign () > 0 ) { // Positive tracks
380+ if (A_PID[0 ] == 0.0 ) { // Pions
381+ MCRegistry.fill (HIST (" hSelectPionP" ), track.pt ());
382+ if (particle.pdgCode () == kPiPlus ) {
383+ MCRegistry.fill (HIST (" hTrueSelectPionP" ), track.pt ());
384+ }
385+ } else if (A_PID[0 ] == 1.0 ) { // Kaons
386+ MCRegistry.fill (HIST (" hSelectKaonP" ), track.pt ());
387+ if (particle.pdgCode () == kKPlus ) {
388+ MCRegistry.fill (HIST (" hTrueSelectKaonP" ), track.pt ());
389+ }
390+ } else if (A_PID[0 ] == 2.0 ) { // Protons
391+ MCRegistry.fill (HIST (" hSelectProtonP" ), track.pt ());
392+ if (particle.pdgCode () == kProton ) {
393+ MCRegistry.fill (HIST (" hTrueSelectProtonP" ), track.pt ());
394+ }
395+ }
396+ } else if (track.sign () < 0 ) { // Negative tracks
397+ if (A_PID[0 ] == 0.0 ) { // Pions
398+ MCRegistry.fill (HIST (" hSelectPionN" ), track.pt ());
399+ if (particle.pdgCode () == kPiMinus ) {
400+ MCRegistry.fill (HIST (" hTrueSelectPionN" ), track.pt ());
401+ }
402+ } else if (A_PID[0 ] == 1.0 ) { // Kaons
403+ MCRegistry.fill (HIST (" hSelectKaonN" ), track.pt ());
404+ if (particle.pdgCode () == kKMinus ) {
405+ MCRegistry.fill (HIST (" hTrueSelectKaonN" ), track.pt ());
406+ }
407+ } else if (A_PID[0 ] == 2.0 ) { // Protons
408+ MCRegistry.fill (HIST (" hSelectProtonN" ), track.pt ());
409+ if (particle.pdgCode () == kProtonBar ) {
410+ MCRegistry.fill (HIST (" hTrueSelectProtonN" ), track.pt ());
411+ }
412+ }
413+ }
414+ }
415415 }
416416 }
417417 // End of the Monte-Carlo reconstructed QA
0 commit comments