@@ -210,7 +210,7 @@ struct ThreePartCorr {
210210 }
211211
212212 if (TrackFilters (track)) {
213- A_PID = TrackPID (track);
213+ A_PID = TrackPID (track);
214214 QARegistry.fill (HIST (" hTrackPt" ), track.pt ());
215215 QARegistry.fill (HIST (" hTrackEta" ), track.eta ());
216216 QARegistry.fill (HIST (" hTrackPhi" ), track.phi ());
@@ -364,64 +364,64 @@ struct ThreePartCorr {
364364 }
365365
366366 if (TrackFilters (track)) {
367- auto particle = track.mcParticle ();
368- if (particle.isPhysicalPrimary ()) {
369-
370- // Efficiency - Reconstructed
371- MCRegistry.fill (HIST (" hReconstructed" ), track.pt ());
372- if (particle.pdgCode () == kPiPlus ) { // Pos pions
373- MCRegistry.fill (HIST (" hRecPionP" ), track.pt ());
374- } else if (particle.pdgCode () == kPiMinus ) { // Neg pions
375- MCRegistry.fill (HIST (" hRecPionN" ), track.pt ());
376- } else if (particle.pdgCode () == kKPlus ) { // Pos kaons
377- MCRegistry.fill (HIST (" hRecKaonP" ), track.pt ());
378- } else if (particle.pdgCode () == kKMinus ) { // Neg kaons
379- MCRegistry.fill (HIST (" hRecKaonN" ), track.pt ());
380- } else if (particle.pdgCode () == kProton ) { // Pos protons
381- MCRegistry.fill (HIST (" hRecProtonP" ), track.pt ());
382- } else if (particle.pdgCode () == kProtonBar ) { // Neg protons
383- MCRegistry.fill (HIST (" hRecProtonN" ), track.pt ());
384- }
385-
386- // Purity (PID)
387- A_PID = TrackPID (track);
388-
389- if (track.sign () > 0 ) { // Positive tracks
390- if (A_PID[0 ] == 0.0 ) { // Pions
391- MCRegistry.fill (HIST (" hSelectPionP" ), track.pt ());
392- if (particle.pdgCode () == kPiPlus ) {
393- MCRegistry.fill (HIST (" hTrueSelectPionP" ), track.pt ());
394- }
395- } else if (A_PID[0 ] == 1.0 ) { // Kaons
396- MCRegistry.fill (HIST (" hSelectKaonP" ), track.pt ());
397- if (particle.pdgCode () == kKPlus ) {
398- MCRegistry.fill (HIST (" hTrueSelectKaonP" ), track.pt ());
399- }
400- } else if (A_PID[0 ] == 2.0 ) { // Protons
401- MCRegistry.fill (HIST (" hSelectProtonP" ), track.pt ());
402- if (particle.pdgCode () == kProton ) {
403- MCRegistry.fill (HIST (" hTrueSelectProtonP" ), track.pt ());
404- }
405- }
406- } else if (track.sign () < 0 ) { // Negative tracks
407- if (A_PID[0 ] == 0.0 ) { // Pions
408- MCRegistry.fill (HIST (" hSelectPionN" ), track.pt ());
409- if (particle.pdgCode () == kPiMinus ) {
410- MCRegistry.fill (HIST (" hTrueSelectPionN" ), track.pt ());
411- }
412- } else if (A_PID[0 ] == 1.0 ) { // Kaons
413- MCRegistry.fill (HIST (" hSelectKaonN" ), track.pt ());
414- if (particle.pdgCode () == kKMinus ) {
415- MCRegistry.fill (HIST (" hTrueSelectKaonN" ), track.pt ());
416- }
417- } else if (A_PID[0 ] == 2.0 ) { // Protons
418- MCRegistry.fill (HIST (" hSelectProtonN" ), track.pt ());
419- if (particle.pdgCode () == kProtonBar ) {
420- MCRegistry.fill (HIST (" hTrueSelectProtonN" ), track.pt ());
421- }
422- }
423- }
424- }
367+ auto particle = track.mcParticle ();
368+ if (particle.isPhysicalPrimary ()) {
369+
370+ // Efficiency - Reconstructed
371+ MCRegistry.fill (HIST (" hReconstructed" ), track.pt ());
372+ if (particle.pdgCode () == kPiPlus ) { // Pos pions
373+ MCRegistry.fill (HIST (" hRecPionP" ), track.pt ());
374+ } else if (particle.pdgCode () == kPiMinus ) { // Neg pions
375+ MCRegistry.fill (HIST (" hRecPionN" ), track.pt ());
376+ } else if (particle.pdgCode () == kKPlus ) { // Pos kaons
377+ MCRegistry.fill (HIST (" hRecKaonP" ), track.pt ());
378+ } else if (particle.pdgCode () == kKMinus ) { // Neg kaons
379+ MCRegistry.fill (HIST (" hRecKaonN" ), track.pt ());
380+ } else if (particle.pdgCode () == kProton ) { // Pos protons
381+ MCRegistry.fill (HIST (" hRecProtonP" ), track.pt ());
382+ } else if (particle.pdgCode () == kProtonBar ) { // Neg protons
383+ MCRegistry.fill (HIST (" hRecProtonN" ), track.pt ());
384+ }
385+
386+ // Purity (PID)
387+ A_PID = TrackPID (track);
388+
389+ if (track.sign () > 0 ) { // Positive tracks
390+ if (A_PID[0 ] == 0.0 ) { // Pions
391+ MCRegistry.fill (HIST (" hSelectPionP" ), track.pt ());
392+ if (particle.pdgCode () == kPiPlus ) {
393+ MCRegistry.fill (HIST (" hTrueSelectPionP" ), track.pt ());
394+ }
395+ } else if (A_PID[0 ] == 1.0 ) { // Kaons
396+ MCRegistry.fill (HIST (" hSelectKaonP" ), track.pt ());
397+ if (particle.pdgCode () == kKPlus ) {
398+ MCRegistry.fill (HIST (" hTrueSelectKaonP" ), track.pt ());
399+ }
400+ } else if (A_PID[0 ] == 2.0 ) { // Protons
401+ MCRegistry.fill (HIST (" hSelectProtonP" ), track.pt ());
402+ if (particle.pdgCode () == kProton ) {
403+ MCRegistry.fill (HIST (" hTrueSelectProtonP" ), track.pt ());
404+ }
405+ }
406+ } else if (track.sign () < 0 ) { // Negative tracks
407+ if (A_PID[0 ] == 0.0 ) { // Pions
408+ MCRegistry.fill (HIST (" hSelectPionN" ), track.pt ());
409+ if (particle.pdgCode () == kPiMinus ) {
410+ MCRegistry.fill (HIST (" hTrueSelectPionN" ), track.pt ());
411+ }
412+ } else if (A_PID[0 ] == 1.0 ) { // Kaons
413+ MCRegistry.fill (HIST (" hSelectKaonN" ), track.pt ());
414+ if (particle.pdgCode () == kKMinus ) {
415+ MCRegistry.fill (HIST (" hTrueSelectKaonN" ), track.pt ());
416+ }
417+ } else if (A_PID[0 ] == 2.0 ) { // Protons
418+ MCRegistry.fill (HIST (" hSelectProtonN" ), track.pt ());
419+ if (particle.pdgCode () == kProtonBar ) {
420+ MCRegistry.fill (HIST (" hTrueSelectProtonN" ), track.pt ());
421+ }
422+ }
423+ }
424+ }
425425 }
426426 }
427427 // End of the Monte-Carlo reconstructed QA
@@ -503,7 +503,7 @@ struct ThreePartCorr {
503503
504504 return ID;
505505 }
506-
506+
507507 // ================================================================================================================================================================================================================
508508
509509 template <class V0Cand >
@@ -534,69 +534,69 @@ struct ThreePartCorr {
534534 }
535535
536536 if (TrackPID (Track)[0 ] == 0.0 ) { // Pions
537- if (TMath::Abs (Track.tpcNSigmaPi ()) > 4.0 ) {
538- return kFALSE ;
537+ if (TMath::Abs (Track.tpcNSigmaPi ()) > 4.0 ) {
538+ return kFALSE ;
539539 }
540540 if (Track.pt () < 0.3 ) {
541- return kFALSE ;
541+ return kFALSE ;
542542 } else if (Track.pt () > 0.3 && Track.pt () < 1.5 ) {
543- if (TMath::Abs (TrackPID (Track)[1 ]) > 4.0 ) {
544- return kFALSE ;
545- }
543+ if (TMath::Abs (TrackPID (Track)[1 ]) > 4.0 ) {
544+ return kFALSE ;
545+ }
546546 } else if (Track.pt () > 1.5 && Track.pt () < 2.3 ) {
547- if (TrackPID (Track)[1 ] < -4.0 || TrackPID (Track)[1 ] > 0.0 ) {
548- return kFALSE ;
549- }
547+ if (TrackPID (Track)[1 ] < -4.0 || TrackPID (Track)[1 ] > 0.0 ) {
548+ return kFALSE ;
549+ }
550550 } else if (Track.pt () > 2.3 ) {
551- return kFALSE ;
551+ return kFALSE ;
552552 }
553-
553+
554554 } else if (TrackPID (Track)[0 ] == 1.0 ) { // Kaons
555- if (TMath::Abs (Track.tpcNSigmaKa ()) > 4.0 ) {
556- return kFALSE ;
555+ if (TMath::Abs (Track.tpcNSigmaKa ()) > 4.0 ) {
556+ return kFALSE ;
557557 }
558558 if (Track.pt () < 0.5 ) {
559- return kFALSE ;
559+ return kFALSE ;
560560 } else if (Track.pt () > 0.5 && Track.pt () < 1.5 ) {
561- if (TMath::Abs (TrackPID (Track)[1 ]) > 4.0 ) {
562- return kFALSE ;
563- }
561+ if (TMath::Abs (TrackPID (Track)[1 ]) > 4.0 ) {
562+ return kFALSE ;
563+ }
564564 } else if (Track.pt () > 1.5 && Track.pt () < 2.0 ) {
565- if (TrackPID (Track)[1 ] < -2.0 || TrackPID (Track)[1 ] > 4.0 ) {
566- return kFALSE ;
567- }
565+ if (TrackPID (Track)[1 ] < -2.0 || TrackPID (Track)[1 ] > 4.0 ) {
566+ return kFALSE ;
567+ }
568568 } else if (Track.pt () > 2.0 && Track.pt () < 2.5 ) {
569- if (TrackPID (Track)[1 ] < 0.0 || TrackPID (Track)[1 ] > 4.0 ) {
570- return kFALSE ;
571- }
569+ if (TrackPID (Track)[1 ] < 0.0 || TrackPID (Track)[1 ] > 4.0 ) {
570+ return kFALSE ;
571+ }
572572 } else if (Track.pt () > 2.5 ) {
573- return kFALSE ;
573+ return kFALSE ;
574574 }
575575
576576 } else if (TrackPID (Track)[0 ] == 2.0 ) { // Protons
577- if (TMath::Abs (Track.tpcNSigmaPr ()) > 4.0 ) {
578- return kFALSE ;
577+ if (TMath::Abs (Track.tpcNSigmaPr ()) > 4.0 ) {
578+ return kFALSE ;
579579 }
580580 if (Track.pt () < 0.5 ) {
581- return kFALSE ;
581+ return kFALSE ;
582582 } else if (Track.pt () > 0.5 && Track.pt () < 0.7 ) {
583- if (TrackPID (Track)[1 ] < -2.0 || TrackPID (Track)[1 ] > 4.0 ) {
584- return kFALSE ;
585- }
583+ if (TrackPID (Track)[1 ] < -2.0 || TrackPID (Track)[1 ] > 4.0 ) {
584+ return kFALSE ;
585+ }
586586 } else if (Track.pt () > 0.7 && Track.pt () < 2.5 ) {
587- if (TMath::Abs (TrackPID (Track)[1 ]) > 4.0 ) {
588- return kFALSE ;
589- }
587+ if (TMath::Abs (TrackPID (Track)[1 ]) > 4.0 ) {
588+ return kFALSE ;
589+ }
590590 } else if (Track.pt () > 2.5 ) {
591- if (TrackPID (Track)[1 ] < -2.0 || TrackPID (Track)[1 ] > 4.0 ) {
592- return kFALSE ;
593- }
591+ if (TrackPID (Track)[1 ] < -2.0 || TrackPID (Track)[1 ] > 4.0 ) {
592+ return kFALSE ;
593+ }
594594 }
595595 }
596596
597597 return kTRUE ;
598598 }
599-
599+
600600 template <class V0Cand , class TrackCand >
601601 Bool_t CorrelationFilters (const V0Cand& V0, const TrackCand& Track)
602602 {
0 commit comments