diff --git a/PWGLF/Tasks/Nuspex/hyperhelium4sigmaAnalysis.cxx b/PWGLF/Tasks/Nuspex/hyperhelium4sigmaAnalysis.cxx index 9a964cf0d62..f9b62c10a6d 100644 --- a/PWGLF/Tasks/Nuspex/hyperhelium4sigmaAnalysis.cxx +++ b/PWGLF/Tasks/Nuspex/hyperhelium4sigmaAnalysis.cxx @@ -34,16 +34,17 @@ using CollisionsFull = soa::Join; using FullTracksExtIU = soa::Join; using MCLabeledTracksIU = soa::Join; -namespace{ - constexpr int kITSLayers = 7; - constexpr int kITSInnerBarrelLayers = 3; - constexpr int kITSOuterBarrelLayers = 4; - std::shared_ptr hMotherCounter; - std::shared_ptr hDauAlphaCounter; - std::shared_ptr hDauTritonCounter; - std::shared_ptr hDauProtonCounter; - std::shared_ptr hDauPionCounter; -} +namespace +{ +constexpr int kITSLayers = 7; +constexpr int kITSInnerBarrelLayers = 3; +constexpr int kITSOuterBarrelLayers = 4; +std::shared_ptr hMotherCounter; +std::shared_ptr hDauAlphaCounter; +std::shared_ptr hDauTritonCounter; +std::shared_ptr hDauProtonCounter; +std::shared_ptr hDauPionCounter; +} // namespace //-------------------------------Check the decay channel of H4S------------------------------- enum Channel { @@ -59,7 +60,7 @@ Channel getDecayChannelH4S(TMCParticle const& particle, std::vector& list) if (std::abs(particle.pdgCode()) != o2::constants::physics::Pdg::kHyperHelium4Sigma) { return kNDecayChannel; } - + // list: charged, charged or empty, neutral list.clear(); list.resize(3, -1); @@ -127,7 +128,7 @@ Channel getDecayChannelH4S(TMCParticle const& particle, std::vector& list) } //-------------------------------------------------------------- -//check if the mcparticle is daughter of hyperhelium4sigma +// check if the mcparticle is daughter of hyperhelium4sigma template bool isDaughterTrack(TMCParticle const& mcparticle, int pdgcode = o2::constants::physics::Pdg::kHyperHelium4Sigma) { @@ -334,7 +335,7 @@ struct Hyperhelium4sigmaQa { hMotherCounter->GetXaxis()->SetBinLabel(4, "has collision"); hMotherCounter->GetXaxis()->SetBinLabel(5, "ITSonly"); hMotherCounter->GetXaxis()->SetBinLabel(6, "ITS hits"); - hMotherCounter->GetXaxis()->SetBinLabel(7, "ITS IR"); + hMotherCounter->GetXaxis()->SetBinLabel(7, "ITS IR"); hMotherCounter->GetXaxis()->SetBinLabel(8, "ITS chi2"); hMotherCounter->GetXaxis()->SetBinLabel(9, "pt"); recoQAHist.add("hTrueMotherRVsDiffPt", ";#Delta p_{T} (GeV/#it{c});R (cm);", HistType::kTH2F, {diffPtAxis, radiusAxis}); @@ -396,7 +397,7 @@ struct Hyperhelium4sigmaQa { // qa for mother track selection template bool motherTrackCheck(const TTrack& track, const std::shared_ptr hist) - { + { hist->Fill(1); if (std::abs(track.eta()) > etaMax) { @@ -404,12 +405,12 @@ struct Hyperhelium4sigmaQa { } hist->Fill(2); - if (!track.has_collision()){ + if (!track.has_collision()) { return false; } hist->Fill(3); - if (!track.hasITS() || track.hasTPC() || track.hasTOF()){ + if (!track.hasITS() || track.hasTPC() || track.hasTOF()) { return false; } hist->Fill(4); @@ -465,7 +466,7 @@ struct Hyperhelium4sigmaQa { if (track.hasTOF()) { return; - } + } hist->Fill(6); }