1616// / \date 2025-04-10
1717
1818#include " PWGCF/Femto/FemtoNuclei/DataModel/PionNucleiTables.h"
19+ #include " PWGCF/FemtoWorld/Core/FemtoWorldMath.h"
1920#include " PWGLF/DataModel/EPCalibrationTables.h"
2021#include " PWGLF/DataModel/LFHypernucleiTables.h"
2122#include " PWGLF/Utils/svPoolCreator.h"
22- #include " PWGCF/FemtoWorld/Core/FemtoWorldMath.h"
2323
2424#include " Common/Core/PID/PIDTOF.h"
2525#include " Common/Core/PID/TPCPIDResponse.h"
@@ -149,7 +149,7 @@ struct PiNucandidate {
149149
150150 // collision information
151151 int32_t collisionID = 0 ;
152- float cent =1 .f;
152+ float cent = 1 .f;
153153};
154154
155155struct PiNucleiFemto {
@@ -283,7 +283,7 @@ struct PiNucleiFemto {
283283 {" hkStaVsmTVsCent_US_M" , " ;kStar (GeV/c);mT (GeV/#it{c}^{2});Centrality" , {HistType::kTH3F , {{300 , 0 .0f , 3 .0f }, {100 , 0 .2f , 3 .2f }, {100 , 0 .0f , 100 .0f }}}},
284284 {" hkStaVsmTVsCent_US_A" , " ;kStar (GeV/c);mT (GeV/#it{c}^{2});Centrality" , {HistType::kTH3F , {{300 , 0 .0f , 3 .0f }, {100 , 0 .2f , 3 .2f }, {100 , 0 .0f , 100 .0f }}}},
285285 {" hCollIDVsCentEachPion" , " ;CollisionID;Centrality" , {HistType::kTH2F , {{4000 , 0 .0f , 4000 .0f }, {100 , 0 .0f , 100 .0f }}}},
286- {" hNHypsPerPrevColl" , " Number of V0Hypers in previous collision used for mixing;N_{V0Hypers};Entries" ,{HistType::kTH2F , {{4000 , 0 .0f , 4000 .0f }, {50 , -0.5 , 49.5 }}}},
286+ {" hNHypsPerPrevColl" , " Number of V0Hypers in previous collision used for mixing;N_{V0Hypers};Entries" , {HistType::kTH2F , {{4000 , 0 .0f , 4000 .0f }, {50 , -0.5 , 49.5 }}}},
287287 {" hkStar_LS_M" , " ;kStar (GeV/c)" , {HistType::kTH1F , {{300 , 0 .0f , 3 .0f }}}},
288288 {" hkStar_LS_A" , " ;kStar (GeV/c)" , {HistType::kTH1F , {{300 , 0 .0f , 3 .0f }}}},
289289 {" hkStar_US_M" , " ;kStar (GeV/c)" , {HistType::kTH1F , {{300 , 0 .0f , 3 .0f }}}},
@@ -295,7 +295,6 @@ struct PiNucleiFemto {
295295 false ,
296296 true };
297297
298-
299298 int numOfCentBins = 40 ;
300299 int numOfVertexZBins = 30 ;
301300 float Vz_low = -10 .0f ;
@@ -485,14 +484,14 @@ struct PiNucleiFemto {
485484 mQaRegistry .fill (HIST (" h2NsigmaPiTPC_preselection" ), candidate.tpcInnerParam (), tpcNSigmaPi);
486485 if (std::abs (candidate.pt ()) < settingCutPiptMin || std::abs (candidate.pt ()) > settingCutPiptMax)
487486 return false ;
488- // reject protons
487+ // reject protons
489488 if (std::abs (candidate.tpcNSigmaPr ()) < settingCutNsigTPCPrMin)
490489 return false ;
491490 mQaRegistry .fill (HIST (" h2NsigmaPiPrTPC" ), candidate.tpcNSigmaPr ());
492491 if (candidate.hasTOF () && std::abs (candidate.tofNSigmaPr ()) < settingCutNsigTOFPrMin)
493492 return false ;
494493 mQaRegistry .fill (HIST (" h2NsigmaPiPrTOF" ), candidate.tofNSigmaPr ());
495-
494+
496495 if (candidate.hasTOF () && candidate.tpcInnerParam () >= settingCutPinMinTOFPi) {
497496 auto tofNSigmaPi = candidate.tofNSigmaPi ();
498497 auto combNsigma = std::sqrt (tofNSigmaPi * tofNSigmaPi + tpcNSigmaPi * tpcNSigmaPi);
@@ -528,7 +527,7 @@ struct PiNucleiFemto {
528527 {
529528 float tpcInnerParam = candidate.tpcInnerParam ();
530529 mQaRegistry .fill (HIST (" h2dEdx" ), candidate.sign () * tpcInnerParam, candidate.tpcSignal ());
531-
530+
532531 float DeDCAxyMin = 0.015 + 0.0305 / TMath::Power (candidate.pt (), 1.1 );
533532 if (abs (candidate.dcaXY ()) > DeDCAxyMin || abs (candidate.dcaXY ()) > settingCutDeDCAzMin)
534533 return false ;
@@ -804,7 +803,7 @@ struct PiNucleiFemto {
804803 bool filledAllOnce = false ;
805804 // LOG(info) << "Number of tracks: " << tracks.size();
806805 for (const auto & track0 : tracks) {
807-
806+
808807 mQaRegistry .fill (HIST (" hTrackSel" ), Selections::kNoCuts );
809808
810809 if (!selectTrack (track0)) {
@@ -835,7 +834,7 @@ struct PiNucleiFemto {
835834 if (!selectTrack (track1) || !selectionPIDPion (track1)) {
836835 continue ;
837836 }
838-
837+
839838 if (!filledAllOnce) {
840839 mQaRegistry .fill (HIST (" hCollIDVsCentEachPion" ), track1.collisionId (), cent);
841840 mQaRegistry .fill (HIST (" hSinglePiPt" ), track1.pt () * track1.sign ());
@@ -849,7 +848,7 @@ struct PiNucleiFemto {
849848 trackPair.collBracket = collBracket;
850849 mTrackPairs .push_back (trackPair);
851850 }
852- filledAllOnce = true ;
851+ filledAllOnce = true ;
853852 }
854853 }
855854
@@ -1237,23 +1236,23 @@ PROCESS_SWITCH(PiNucleiFemto, processMixedEventHyper, "Process Mixed event", fal
12371236 mQaRegistry .fill (HIST (" hSkipReasons" ), 0 );
12381237 continue ;
12391238 }
1240-
1239+
12411240 auto hypdTablepreviousCollision = V0Hypers.sliceBy (hypPerCol, c2.globalIndex ());
12421241 hypdTablepreviousCollision.bindExternalIndices (&V0Hypers);
1243-
1242+
12441243 if (hypdTablepreviousCollision.size () == 0 ) {
12451244 mQaRegistry .fill (HIST (" hSkipReasons" ), 1 );
12461245 continue ;
12471246 }
1248-
1247+
12491248 auto firstHyp = hypdTablepreviousCollision.iteratorAt (0 );
12501249 int poolIndexHyp = where_pool (firstHyp.zPrimVtx (), firstHyp.centralityFT0C ());
1251-
1250+
12521251 if (poolIndexHyp != poolIndexPi) {
1253- mQaRegistry .fill (HIST (" hSkipReasons" ), 2 );
1254- continue ;
1252+ mQaRegistry .fill (HIST (" hSkipReasons" ), 2 );
1253+ continue ;
12551254 }
1256-
1255+
12571256 auto trackTableThisCollision = pitracks.sliceBy (mPerCol , collision.globalIndex ());
12581257 trackTableThisCollision.bindExternalIndices (&pitracks);
12591258
@@ -1265,7 +1264,7 @@ PROCESS_SWITCH(PiNucleiFemto, processMixedEventHyper, "Process Mixed event", fal
12651264 pool.events .pop_front ();
12661265 }
12671266 pool.events .push_back ({collision.globalIndex ()});
1268- }
1267+ }
12691268 fillPairsHyper (collisions, pitracks, V0Hypers, /* isMixedEvent*/ true );
12701269 }
12711270 PROCESS_SWITCH (PiNucleiFemto, processMixedEventHyperPool, " Process Mixed event" , false );
0 commit comments