2222#include " Common/DataModel/PIDResponseTPC.h"
2323#include " Common/DataModel/TrackSelectionTables.h"
2424
25+ #include " PID/PIDTOFParamService.h"
2526#include < CCDB/BasicCCDBManager.h>
2627#include < CommonConstants/MathConstants.h>
2728#include < CommonConstants/PhysicsConstants.h>
3233#include < Framework/AnalysisDataModel.h>
3334#include < Framework/AnalysisHelpers.h>
3435#include < Framework/AnalysisTask.h>
36+ #include < Framework/Concepts.h>
3537#include < Framework/Configurable.h>
36- #include < Framework/DataTypes.h>
3738#include < Framework/HistogramRegistry.h>
3839#include < Framework/HistogramSpec.h>
3940#include < Framework/InitContext.h>
4041#include < Framework/OutputObjHeader.h>
4142#include < Framework/runDataProcessing.h>
43+ #include < ReconstructionDataFormats/DCA.h>
4244
4345#include < Math/Vector4D.h> // IWYU pragma: keep (do not replace with Math/Vector4Dfwd.h)
4446#include < Math/Vector4Dfwd.h>
@@ -216,11 +218,11 @@ struct skimmerPrimaryElectronFromDalitzEE {
216218 fRegistry .add (" Pair/uls/hMCutMvsPt" , " m_{ee} vs. p_{T,ee};m_{ee} (GeV/c^{2});p_{T,ee} (GeV/c)" , kTH2F , {{100 , 0 , 0.1 }, {200 , 0 , 2 }}, false );
217219 fRegistry .add (" Pair/uls/hMPhiCutMvsPt" , " m_{ee} vs. p_{T,ee};m_{ee} (GeV/c^{2});p_{T,ee} (GeV/c)" , kTH2F , {{100 , 0 , 0.1 }, {200 , 0 , 2 }}, false );
218220
219- fRegistry .add (" Pair/uls/hTrackMvsPhiV" , " m_{ee} vs. #varphi_{V};#varphi_{V} (rad.);m_{ee} (GeV/c^{2})" , kTH2F , {{180 , 0 , M_PI }, {100 , 0 , 0.1 }}, false );
220- fRegistry .add (" Pair/uls/hCheckEMvsPhiV" , " m_{ee} vs. #varphi_{V};#varphi_{V} (rad.);m_{ee} (GeV/c^{2})" , kTH2F , {{180 , 0 , M_PI }, {100 , 0 , 0.1 }}, false );
221- fRegistry .add (" Pair/uls/hMvsPhiV" , " m_{ee} vs. #varphi_{V};#varphi_{V} (rad.);m_{ee} (GeV/c^{2})" , kTH2F , {{180 , 0 , M_PI }, {100 , 0 , 0.1 }}, false );
222- fRegistry .add (" Pair/uls/hMCutMvsPhiV" , " m_{ee} vs. #varphi_{V};#varphi_{V} (rad.);m_{ee} (GeV/c^{2})" , kTH2F , {{180 , 0 , M_PI }, {100 , 0 , 0.1 }}, false );
223- fRegistry .add (" Pair/uls/hMPhiCutMvsPhiV" , " m_{ee} vs. #varphi_{V};#varphi_{V} (rad.);m_{ee} (GeV/c^{2})" , kTH2F , {{180 , 0 , M_PI }, {100 , 0 , 0.1 }}, false );
221+ fRegistry .add (" Pair/uls/hTrackMvsPhiV" , " m_{ee} vs. #varphi_{V};#varphi_{V} (rad.);m_{ee} (GeV/c^{2})" , kTH2F , {{180 , 0 , o2::constants::math:: PI }, {100 , 0 , 0.1 }}, false );
222+ fRegistry .add (" Pair/uls/hCheckEMvsPhiV" , " m_{ee} vs. #varphi_{V};#varphi_{V} (rad.);m_{ee} (GeV/c^{2})" , kTH2F , {{180 , 0 , o2::constants::math:: PI }, {100 , 0 , 0.1 }}, false );
223+ fRegistry .add (" Pair/uls/hMvsPhiV" , " m_{ee} vs. #varphi_{V};#varphi_{V} (rad.);m_{ee} (GeV/c^{2})" , kTH2F , {{180 , 0 , o2::constants::math:: PI }, {100 , 0 , 0.1 }}, false );
224+ fRegistry .add (" Pair/uls/hMCutMvsPhiV" , " m_{ee} vs. #varphi_{V};#varphi_{V} (rad.);m_{ee} (GeV/c^{2})" , kTH2F , {{180 , 0 , o2::constants::math:: PI }, {100 , 0 , 0.1 }}, false );
225+ fRegistry .add (" Pair/uls/hMPhiCutMvsPhiV" , " m_{ee} vs. #varphi_{V};#varphi_{V} (rad.);m_{ee} (GeV/c^{2})" , kTH2F , {{180 , 0 , o2::constants::math:: PI }, {100 , 0 , 0.1 }}, false );
224226
225227 fRegistry .addClone (" Pair/uls/" , " Pair/lspp/" );
226228 fRegistry .addClone (" Pair/uls/" , " Pair/lsmm/" );
@@ -377,7 +379,7 @@ struct skimmerPrimaryElectronFromDalitzEE {
377379 return false ;
378380 }
379381
380- if (!includeITSsa && ( !track.hasITS () || !track. hasTPC () )) {
382+ if (!includeITSsa && !track.hasTPC ()) {
381383 return false ;
382384 }
383385
@@ -698,7 +700,6 @@ struct skimmerPrimaryElectronFromDalitzEE {
698700 std::vector<int > acceptedPosTrackIds_per_collision;
699701 std::vector<int > acceptedNegTrackIds_per_collision;
700702 std::vector<int > acceptedTrackIds_per_collision;
701- std::vector<std::pair<int , int >> stored_trackIds;
702703 // Filter trackFilter = minpt < o2::aod::track::pt && nabs(o2::aod::track::eta) < maxeta && o2::aod::track::itsChi2NCl < maxchi2its && ncheckbit(aod::track::v001::detectorMap, (uint8_t)o2::aod::track::ITS) == true && nabs(o2::aod::track::dcaXY) < dca_xy_max && nabs(o2::aod::track::dcaZ) < dca_z_max;
703704 // Filter trackFilter
704705 // using MyFilteredTracks = soa::Filtered<MyTracks>;
@@ -718,11 +719,8 @@ struct skimmerPrimaryElectronFromDalitzEE {
718719 mTOFResponse ->processSetup (bcs.iteratorAt (0 ));
719720
720721 for (const auto & track : tracks) {
721- if (mapCollisionTime.find (track.collisionId ()) == mapCollisionTime.end ()) {
722- // LOGF(info, "track.collisionId() = %d, track.tofEvTime() = %f, track.tofEvTimeErr() = %f", track.collisionId(), track.tofEvTime(), track.tofEvTimeErr());
723- mapCollisionTime[track.collisionId ()] = track.tofEvTime ();
724- mapCollisionTimeError[track.collisionId ()] = track.tofEvTimeErr ();
725- }
722+ mapCollisionTime.try_emplace (track.collisionId (), track.tofEvTime ());
723+ mapCollisionTimeError.try_emplace (track.collisionId (), track.tofEvTimeErr ());
726724 }
727725 calculateTOFNSigmaWithReassociation<true >(collisions, bcs, tracks, trackIndices);
728726
@@ -772,11 +770,8 @@ struct skimmerPrimaryElectronFromDalitzEE {
772770
773771 // mTOFResponse->processSetup(bcs.iteratorAt(0));
774772 // for (const auto& track : tracks) {
775- // if (mapCollisionTime.find(track.collisionId()) == mapCollisionTime.end()) {
776- // // LOGF(info, "track.collisionId() = %d, track.tofEvTime() = %f, track.tofEvTimeErr() = %f", track.collisionId(), track.tofEvTime(), track.tofEvTimeErr());
777- // mapCollisionTime[track.collisionId()] = track.tofEvTime();
778- // mapCollisionTimeError[track.collisionId()] = track.tofEvTimeErr();
779- // }
773+ // mapCollisionTime.try_emplace(track.collisionId(), track.tofEvTime());
774+ // mapCollisionTimeError.try_emplace(track.collisionId(), track.tofEvTimeErr());
780775 // }
781776 // calculateTOFNSigmaWithReassociation<true>(collisions, bcs, tracks, trackIndices);
782777
@@ -840,11 +835,8 @@ struct skimmerPrimaryElectronFromDalitzEE {
840835
841836 mTOFResponse ->processSetup (bcs.iteratorAt (0 ));
842837 for (const auto & track : tracks) {
843- if (mapCollisionTime.find (track.collisionId ()) == mapCollisionTime.end ()) {
844- // LOGF(info, "track.collisionId() = %d, track.tofEvTime() = %f, track.tofEvTimeErr() = %f", track.collisionId(), track.tofEvTime(), track.tofEvTimeErr());
845- mapCollisionTime[track.collisionId ()] = track.tofEvTime ();
846- mapCollisionTimeError[track.collisionId ()] = track.tofEvTimeErr ();
847- }
838+ mapCollisionTime.try_emplace (track.collisionId (), track.tofEvTime ());
839+ mapCollisionTimeError.try_emplace (track.collisionId (), track.tofEvTimeErr ());
848840 }
849841 calculateTOFNSigmaWithReassociation<true >(collisions, bcs, tracks, trackIndices);
850842
0 commit comments