1313// / \brief QA for MC productions which contain hypertriton 3body decay process, including special checks for TOF PID
1414// / \author Yuanzhe Wang <yuanzhe.wang@cern.ch>
1515
16- #include < cmath>
17- #include < array>
18- #include < cstdlib>
19- #include < string>
20- #include < vector>
21- #include < algorithm>
22- #include < unordered_set>
23-
2416#include " PWGLF/DataModel/LFPIDTOFGenericTables.h"
17+ #include " PWGLF/DataModel/LFStrangenessTables.h"
2518#include " PWGLF/Utils/pidTOFGeneric.h"
2619
27- #include " CommonDataFormat/InteractionRecord.h"
28- #include " CommonDataFormat/IRFrame.h"
29- #include " Framework/runDataProcessing.h"
30- #include " Framework/AnalysisTask.h"
31- #include " Framework/AnalysisDataModel.h"
32- #include " Framework/ASoAHelpers.h"
33- #include " ReconstructionDataFormats/Track.h"
3420#include " Common/Core/RecoDecay.h"
35- #include " Common/Core/trackUtilities.h"
36- #include " PWGLF/DataModel/LFStrangenessTables.h"
3721#include " Common/Core/TrackSelection.h"
38- #include " Common/DataModel/TrackSelectionTables.h"
39- #include " Common/TableProducer/PID/pidTOFBase.h"
40- #include " Common/DataModel/EventSelection.h"
22+ #include " Common/Core/trackUtilities.h"
4123#include " Common/DataModel/Centrality.h"
24+ #include " Common/DataModel/EventSelection.h"
4225#include " Common/DataModel/PIDResponse.h"
43- #include " CommonConstants/PhysicsConstants.h"
26+ #include " Common/DataModel/TrackSelectionTables.h"
27+ #include " Common/TableProducer/PID/pidTOFBase.h"
28+
4429#include " CCDB/BasicCCDBManager.h"
30+ #include " CommonConstants/PhysicsConstants.h"
31+ #include " CommonDataFormat/IRFrame.h"
32+ #include " CommonDataFormat/InteractionRecord.h"
33+ #include " Framework/ASoAHelpers.h"
34+ #include " Framework/AnalysisDataModel.h"
35+ #include " Framework/AnalysisTask.h"
36+ #include " Framework/runDataProcessing.h"
37+ #include " ReconstructionDataFormats/Track.h"
38+
39+ #include < algorithm>
40+ #include < array>
41+ #include < cmath>
42+ #include < cstdlib>
43+ #include < string>
44+ #include < unordered_set>
45+ #include < vector>
4546
4647using namespace o2 ;
4748using namespace o2 ::framework;
@@ -59,7 +60,6 @@ namespace
5960constexpr float kCSPEED = TMath::C() * 1 .0e2f * 1 .0e-12f ; // c in cm/ps
6061} // namespace
6162
62-
6363template <class TMCTrackTo , typename TMCParticle>
6464bool is3bodyDecayedH3L (TMCParticle const & particle)
6565{
@@ -328,7 +328,7 @@ struct Hypertriton3bodyQa {
328328 }
329329 registry.fill (HIST (" hEventCounter" ), 2.5 );
330330
331- std::vector<int > indicesProton, indicesPion, indicesDeuteron; // indices for daughter tracks
331+ std::vector<int > indicesProton, indicesPion, indicesDeuteron; // indices for daughter tracks
332332 std::unordered_set<int64_t > globalIDProton, globalIDPion, globalIDDeuteron; // check duplicated daughters
333333 int itrack = -1 ;
334334
@@ -815,16 +815,15 @@ struct Hypertriton3bodyMcParticleCheck {
815815 if (isMatter) {
816816 registry.fill (HIST (" hMcHypertritonCounter" ), 3.5 );
817817 registry.fill (HIST (" hMcHypertritonCounter" ), 4.5 );
818- }
819- else {
818+ } else {
820819 registry.fill (HIST (" hMcHypertritonCounter" ), 3.5 );
821820 registry.fill (HIST (" hMcHypertritonCounter" ), 5.5 );
822821 }
823822 double hypertritonMCMass = RecoDecay::m (std::array{std::array{dauProtonMom[0 ], dauProtonMom[1 ], dauProtonMom[2 ]}, std::array{dauPionMom[0 ], dauPionMom[1 ], dauPionMom[2 ]}, std::array{dauDeuteronMom[0 ], dauDeuteronMom[1 ], dauDeuteronMom[2 ]}}, std::array{o2::constants::physics::MassProton, o2::constants::physics::MassPionCharged, o2::constants::physics::MassDeuteron});
824823 registry.fill (HIST (" hMcRecoInvMass" ), hypertritonMCMass);
825824
826- mcLifetime = RecoDecay::sqrtSumOfSquares (dauDeuteronPos[0 ] - mcparticle.vx (), dauDeuteronPos[1 ] - mcparticle.vy (), dauDeuteronPos[2 ] - mcparticle.vz ()) * o2::constants::physics::MassHyperTriton / mcparticle.p ();
827- registry.fill (HIST (" h3dMCDecayedHypertriton" ), mcparticle.y (), mcparticle.pt (), mcLifetime);
825+ mcLifetime = RecoDecay::sqrtSumOfSquares (dauDeuteronPos[0 ] - mcparticle.vx (), dauDeuteronPos[1 ] - mcparticle.vy (), dauDeuteronPos[2 ] - mcparticle.vz ()) * o2::constants::physics::MassHyperTriton / mcparticle.p ();
826+ registry.fill (HIST (" h3dMCDecayedHypertriton" ), mcparticle.y (), mcparticle.pt (), mcLifetime);
828827
829828 double diffTrackR = dauDeuteronTrackR - std::min (dauPionTrackR, dauProtonTrackR);
830829 registry.fill (HIST (" hDiffDaughterR" ), diffTrackR);
0 commit comments