@@ -53,12 +53,16 @@ struct QAHistTask {
5353 void init (o2::framework::InitContext&)
5454 {
5555
56- if ((process_proton == true && (process_deuteron == true || process_triton == true || process_He3 == true || process_He4 == true )) || (process_deuteron == true && (process_triton == true || process_He3 == true || process_He4 == true )) || (process_triton == true && (process_He3 == true || process_He4 == true )) || (process_He3 == true && process_He4 == true )) {
56+ if ((process_pion == true && (process_kaon == true || process_proton == true || process_deuteron == true || process_triton == true || process_He3 == true || process_He4 == true )) || (process_kaon == true && (process_proton == true || process_deuteron == true || process_triton == true || process_He3 == true || process_He4 == true )) || ( process_proton == true && (process_deuteron == true || process_triton == true || process_He3 == true || process_He4 == true )) || (process_deuteron == true && (process_triton == true || process_He3 == true || process_He4 == true )) || (process_triton == true && (process_He3 == true || process_He4 == true )) || (process_He3 == true && process_He4 == true )) {
5757 LOG (fatal) << " ++++++++ Can't enable more than one species at a time, use subwagons for that purpose. ++++++++" ;
5858 }
5959
6060 std::string species;
6161
62+ if (process_pion)
63+ species = " pi" ;
64+ if (process_kaon)
65+ species = " ka" ;
6266 if (process_proton)
6367 species = " p" ;
6468 if (process_deuteron)
@@ -198,6 +202,8 @@ struct QAHistTask {
198202 }
199203
200204 // Configurables
205+ Configurable<bool > process_pion{" process_pion" , false , " 0: disabled, 1: enabled" };
206+ Configurable<bool > process_kaon{" process_kaon" , false , " 0: disabled, 1: enabled" };
201207 Configurable<bool > process_proton{" process_proton" , false , " 0: disabled, 1: enabled" };
202208 Configurable<bool > process_deuteron{" process_deuteron" , false , " 0: disabled, 1: enabled" };
203209 Configurable<bool > process_triton{" process_triton" , false , " 0: disabled, 1: enabled" };
@@ -255,6 +261,10 @@ struct QAHistTask {
255261
256262 float nSigmaSpecies = 999.0 ;
257263
264+ if (process_pion)
265+ nSigmaSpecies = track.tpcNSigmaPi ();
266+ if (process_kaon)
267+ nSigmaSpecies = track.tpcNSigmaKa ();
258268 if (process_proton)
259269 nSigmaSpecies = track.tpcNSigmaPr ();
260270 if (process_deuteron)
@@ -360,7 +370,7 @@ struct QAHistTask {
360370 QA_reg.fill (HIST (" histTofSignalData" ), track.tpcInnerParam () * track.sign (), track.beta ());
361371 }
362372
363- // fill QA histograms (proton)
373+ // fill QA histograms
364374 if (TMath::Abs (nSigmaSpecies) < nsigmacut) {
365375 if (track.sign () > 0 ) {
366376 QA_species_pos.fill (HIST (" histDcaVsPtData" ), track.pt (), track.dcaXY ());
@@ -483,11 +493,11 @@ struct QAHistTask {
483493
484494 using EventCandidatesDataCent = soa::Filtered<soa::Join<aod::Collisions, aod::EvSels, aod::CentFT0Cs>>;
485495
486- using TrackCandidatesData = soa::Filtered<soa::Join<aod::Tracks, aod::TracksExtra, aod::TracksDCA, aod::pidTPCLfFullPr, aod::pidTOFFullPr, aod::pidTPCLfFullDe, aod::pidTOFFullDe, aod::pidTPCLfFullTr, aod::pidTOFFullTr, aod::pidTPCLfFullHe, aod::pidTOFFullHe, aod::pidTPCLfFullAl, aod::pidTOFFullAl, aod::TrackSelection, aod::TrackSelectionExtension, aod::TOFSignal, aod::pidTOFmass, aod::pidTOFbeta>>;
496+ using TrackCandidatesData = soa::Filtered<soa::Join<aod::Tracks, aod::TracksExtra, aod::TracksDCA, aod::pidTPCLfFullPi, aod::pidTOFFullPi, aod::pidTPCLfFullKa, aod::pidTOFFullKa, aod:: pidTPCLfFullPr, aod::pidTOFFullPr, aod::pidTPCLfFullDe, aod::pidTOFFullDe, aod::pidTPCLfFullTr, aod::pidTOFFullTr, aod::pidTPCLfFullHe, aod::pidTOFFullHe, aod::pidTPCLfFullAl, aod::pidTOFFullAl, aod::TrackSelection, aod::TrackSelectionExtension, aod::TOFSignal, aod::pidTOFmass, aod::pidTOFbeta>>;
487497
488498 using EventCandidatesMC = soa::Filtered<soa::Join<aod::Collisions, aod::EvSels, aod::CentFT0Cs>>;
489499
490- using TrackCandidatesMC = soa::Filtered<soa::Join<aod::Tracks, aod::TracksExtra, aod::TracksDCA, aod::pidTPCLfFullPr, aod::pidTOFFullPr, aod::pidTPCLfFullDe, aod::pidTOFFullDe, aod::pidTPCLfFullTr, aod::pidTOFFullTr, aod::pidTPCLfFullHe, aod::pidTOFFullHe, aod::pidTPCLfFullAl, aod::pidTOFFullAl, aod::TrackSelection, aod::TrackSelectionExtension, aod::TOFSignal, aod::pidTOFmass, aod::pidTOFbeta>>;
500+ using TrackCandidatesMC = soa::Filtered<soa::Join<aod::Tracks, aod::TracksExtra, aod::TracksDCA, aod::pidTPCLfFullPi, aod::pidTOFFullPi, aod::pidTPCLfFullKa, aod::pidTOFFullKa, aod:: pidTPCLfFullPr, aod::pidTOFFullPr, aod::pidTPCLfFullDe, aod::pidTOFFullDe, aod::pidTPCLfFullTr, aod::pidTOFFullTr, aod::pidTPCLfFullHe, aod::pidTOFFullHe, aod::pidTPCLfFullAl, aod::pidTOFFullAl, aod::TrackSelection, aod::TrackSelectionExtension, aod::TOFSignal, aod::pidTOFmass, aod::pidTOFbeta>>;
491501
492502 void processData (EventCandidatesData::iterator const & event, TrackCandidatesData const & tracks)
493503 {
0 commit comments