99// granted to it by virtue of its status as an Intergovernmental Organization
1010// or submit itself to any jurisdiction.
1111//
12- // Task to select electrons from dalitz decay
13- // It can produce track and pair histograms for selected tracks
14- // It creates a bitmap with selections to be stored during skimming
12+ //
13+ // \file DalitzSelection.cxx
14+ // \author Gauthier Legras, glegras@uni-muenster.de, gauthier.legras@cern.ch
15+ // \brief Task to select electrons from dalitz decay
16+ // It can produce track and pair histograms for selected tracks
17+ // It creates a bitmap with selections to be stored during skimming
1518//
1619#include " PWGDQ/Core/AnalysisCompositeCut.h"
1720#include " PWGDQ/Core/AnalysisCut.h"
@@ -65,10 +68,10 @@ using MyBarrelTracks = soa::Join<aod::Tracks, aod::TracksExtra, aod::TrackSelect
6568 aod::pidTOFFullEl, aod::pidTOFFullPi, aod::pidTOFFullMu,
6669 aod::pidTOFFullKa, aod::pidTOFFullPr, aod::pidTOFbeta>;
6770
68- constexpr static uint32_t gkEventFillMap = VarManager::ObjTypes::Collision;
69- constexpr static uint32_t gkTrackFillMap = VarManager::ObjTypes::Track | VarManager::ObjTypes::TrackExtra | VarManager::ObjTypes::TrackDCA | VarManager::ObjTypes::TrackSelection | VarManager::ObjTypes::TrackPID;
71+ constexpr static uint32_t EventFillMap = VarManager::ObjTypes::Collision;
72+ constexpr static uint32_t TrackFillMap = VarManager::ObjTypes::Track | VarManager::ObjTypes::TrackExtra | VarManager::ObjTypes::TrackDCA | VarManager::ObjTypes::TrackSelection | VarManager::ObjTypes::TrackPID;
7073
71- struct dalitzPairing {
74+ struct DalitzSelection {
7275 Produces<o2::aod::DalitzBits> dalitzbits;
7376 Preslice<MyBarrelTracks> perCollision = aod::track::collisionId;
7477
@@ -154,7 +157,7 @@ struct dalitzPairing {
154157 TString addTrackCutsStr = fConfigCuts .fConfigTrackCutsJSON .value ;
155158 if (addTrackCutsStr != " " ) {
156159 std::vector<AnalysisCut*> addTrackCuts = dqcuts::GetCutsFromJSON (addTrackCutsStr.Data ());
157- for (auto & t : addTrackCuts) {
160+ for (const auto & t : addTrackCuts) {
158161 fTrackCuts .push_back (reinterpret_cast <AnalysisCompositeCut*>(t));
159162 }
160163 }
@@ -172,7 +175,7 @@ struct dalitzPairing {
172175 // extra cuts via JSON
173176 if (addTrackCutsProbeStr != " " && (cutNamesProbeStr.CompareTo (cutNamesStr) != 0 || addTrackCutsProbeStr.CompareTo (addTrackCutsStr) != 0 )) {
174177 std::vector<AnalysisCut*> addTrackCuts = dqcuts::GetCutsFromJSON (addTrackCutsProbeStr.Data ());
175- for (auto & t : addTrackCuts) {
178+ for (const auto & t : addTrackCuts) {
176179 fTrackCutsProbe .push_back (reinterpret_cast <AnalysisCompositeCut*>(t));
177180 }
178181 fIsTagAndProbe = true ;
@@ -190,7 +193,7 @@ struct dalitzPairing {
190193 TString addPairCutsStr = fConfigCuts .fConfigPairCutsJSON .value ;
191194 if (addPairCutsStr != " " ) {
192195 std::vector<AnalysisCut*> addPairCuts = dqcuts::GetCutsFromJSON (addPairCutsStr.Data ());
193- for (auto & t : addPairCuts) {
196+ for (const auto & t : addPairCuts) {
194197 fPairCuts .push_back (reinterpret_cast <AnalysisCompositeCut*>(t));
195198 }
196199 }
@@ -294,7 +297,7 @@ struct dalitzPairing {
294297 template <uint32_t TTrackFillMap, typename TTracks>
295298 void runTrackSelection (TTracks const & tracksBarrel)
296299 {
297- for (auto & track : tracksBarrel) {
300+ for (const auto & track : tracksBarrel) {
298301 uint8_t filterMap = uint8_t (0 );
299302 uint8_t filterMapProbe = uint8_t (0 );
300303 VarManager::FillTrack<TTrackFillMap>(track);
@@ -324,7 +327,7 @@ struct dalitzPairing {
324327 template <int TPairType, uint32_t TTrackFillMap, typename TTracks>
325328 void runDalitzPairing (TTracks const & tracks1, TTracks const & tracks2)
326329 {
327- for (auto & [track1, track2] : o2::soa::combinations (CombinationsFullIndexPolicy (tracks1, tracks2))) {
330+ for (const auto & [track1, track2] : o2::soa::combinations (CombinationsFullIndexPolicy (tracks1, tracks2))) {
328331 if (track1.globalIndex () == track2.globalIndex ()) {
329332 continue ;
330333 }
@@ -375,7 +378,7 @@ struct dalitzPairing {
375378
376379 // Fill Hists
377380 if (fConfigOptions .fQA ) {
378- for (auto & track : tracks1) {
381+ for (const auto & track : tracks1) {
379382 uint8_t filterMap = fDalitzmap [track.globalIndex ()];
380383 uint8_t filterMapProbe = fDalitzmapProbe [track.globalIndex ()];
381384 if (!filterMap && !filterMapProbe) {
@@ -405,11 +408,11 @@ struct dalitzPairing {
405408 fDalitzmap .clear ();
406409 fDalitzmapProbe .clear ();
407410
408- for (auto & collision : collisions) {
411+ for (const auto & collision : collisions) {
409412 fTrackmap .clear ();
410413 fTrackmapProbe .clear ();
411414 VarManager::ResetValues (VarManager::kNRunWiseVariables , VarManager::kNBarrelTrackVariables );
412- VarManager::FillEvent<gkEventFillMap >(collision);
415+ VarManager::FillEvent<EventFillMap >(collision);
413416 bool isEventSelected = fEventCut ->IsSelected (VarManager::fgValues);
414417
415418 if (isEventSelected) {
@@ -452,26 +455,26 @@ struct dalitzPairing {
452455 }
453456
454457 auto groupedFilteredTracks = filteredTracks.sliceBy (perCollision, collision.globalIndex ());
455- runTrackSelection<gkTrackFillMap >(groupedFilteredTracks);
456- runDalitzPairing<pairType, gkTrackFillMap >(groupedFilteredTracks, groupedFilteredTracks);
458+ runTrackSelection<TrackFillMap >(groupedFilteredTracks);
459+ runDalitzPairing<pairType, TrackFillMap >(groupedFilteredTracks, groupedFilteredTracks);
457460 }
458461 }
459462
460- for (auto & track : tracks) { // Fill dalitz bits
463+ for (const auto & track : tracks) { // Fill dalitz bits
461464 dalitzbits (fDalitzmap [track.globalIndex ()]);
462465 }
463466 }
464467
465- void processDummy (MyEvents&)
468+ void processDummy (MyEvents const &)
466469 {
467470 }
468471
469- PROCESS_SWITCH (dalitzPairing , processFullTracks, " Run Dalitz selection on AO2D tables" , false );
470- PROCESS_SWITCH (dalitzPairing , processDummy, " Do nothing" , false );
472+ PROCESS_SWITCH (DalitzSelection , processFullTracks, " Run Dalitz selection on AO2D tables" , false );
473+ PROCESS_SWITCH (DalitzSelection , processDummy, " Do nothing" , false );
471474};
472475
473476WorkflowSpec defineDataProcessing (ConfigContext const & cfgc)
474477{
475478 return WorkflowSpec{
476- adaptAnalysisTask<dalitzPairing >(cfgc)};
479+ adaptAnalysisTask<DalitzSelection >(cfgc)};
477480}
0 commit comments