1515// / \author prottay.das@cern.ch
1616
1717#include " PWGLF/DataModel/LFCKSSpinalignmentTables.h"
18- #include " PWGLF/DataModel/LFSpincorrelationTables.h"
1918
2019#include " Common/Core/RecoDecay.h"
2120#include " Common/Core/trackUtilities.h"
3231#include < Math/GenVector/Boost.h>
3332#include < Math/Vector3D.h>
3433#include < Math/Vector4D.h>
35- #include < TLorentzVector.h>
3634#include < TMath.h>
3735
3836#include < fairlogger/Logger.h>
3937
4038#include < cmath> // for std::fabs
4139#include < deque>
42- #include < iostream>
40+ // #include <iostream>
4341#include < iterator>
4442#include < set> // <<< CHANGED: for dedup sets
4543#include < string>
@@ -62,7 +60,7 @@ struct cksspinalignder {
6260
6361 struct : ConfigurableGroup {
6462 Configurable<std::string> cfgURL{" cfgURL" , " http://alice-ccdb.cern.ch" , " Address of the CCDB to browse" };
65- Configurable<int64_t > nolaterthan{ " ccdb-no-later-than " , std::chrono::duration_cast<std::chrono::milliseconds>(std::chrono::system_clock::now ().time_since_epoch ()).count (), " Latest acceptable timestamp of creation for the object" };
63+ Configurable<int64_t > ccdbNoLaterThan{ " ccdbNoLaterThan " , std::chrono::duration_cast<std::chrono::milliseconds>(std::chrono::system_clock::now ().time_since_epoch ()).count (), " Latest acceptable timestamp of creation for the object" };
6664 } cfgCcdbParam;
6765
6866 // event sel/////////
@@ -85,8 +83,8 @@ struct cksspinalignder {
8583
8684 // Event Mixing
8785 Configurable<int > nEvtMixing{" nEvtMixing" , 5 , " Number of events to mix" };
88- ConfigurableAxis CfgVtxBins{ " CfgVtxBins " , {10 , -10 , 10 }, " Mixing bins - z-vertex" };
89- ConfigurableAxis CfgMultBins{ " CfgMultBins " , {8 , 0.0 , 80 }, " Mixing bins - centrality" };
86+ ConfigurableAxis cfgVtxBins{ " cfgVtxBins " , {10 , -10 , 10 }, " Mixing bins - z-vertex" };
87+ ConfigurableAxis cfgMultBins{ " cfgMultBins " , {8 , 0.0 , 80 }, " Mixing bins - centrality" };
9088
9189 // THnsparse bining
9290 ConfigurableAxis configThnAxisInvMass{" configThnAxisInvMass" , {50 , 1.09 , 1.14 }, " #it{M} (GeV/#it{c}^{2})" };
@@ -162,13 +160,14 @@ struct cksspinalignder {
162160 auto py = candidate.pionBachPy ();
163161 auto pz = candidate.pionBachPz ();
164162 auto p = std::sqrt (px * px + py * py + pz * pz);
165- if (p < 0.5 ) {
166- if (!candidate.pionBachTOFHit () && TMath::Abs (candidate.pionBachTPC ()) < nsigmaCutTPC) {
163+ float lowmom = 0.5 ;
164+ if (p < lowmom) {
165+ if (!candidate.pionBachTOFHit () && std::abs (candidate.pionBachTPC ()) < nsigmaCutTPC) {
167166 return true ;
168- } else if (candidate.pionBachTOFHit () && TMath::Sqrt (candidate.pionBachTPC () * candidate.pionBachTPC () + candidate.pionBachTOF () * candidate.pionBachTOF ()) < nsigmaCutTOF) {
167+ } else if (candidate.pionBachTOFHit () && std::sqrt (candidate.pionBachTPC () * candidate.pionBachTPC () + candidate.pionBachTOF () * candidate.pionBachTOF ()) < nsigmaCutTOF) {
169168 return true ;
170169 }
171- } else if (candidate.pionBachTOFHit () && TMath::Sqrt (candidate.pionBachTPC () * candidate.pionBachTPC () + candidate.pionBachTOF () * candidate.pionBachTOF ()) < nsigmaCutTOF) {
170+ } else if (candidate.pionBachTOFHit () && std::sqrt (candidate.pionBachTPC () * candidate.pionBachTPC () + candidate.pionBachTOF () * candidate.pionBachTOF ()) < nsigmaCutTOF) {
172171 return true ;
173172 }
174173 return false ;
@@ -234,13 +233,13 @@ struct cksspinalignder {
234233 // Processing Event Mixing
235234 SliceCache cache;
236235 using BinningType = ColumnBinningPolicy<aod::kshortpionevent::Posz, aod::kshortpionevent::Cent>;
237- BinningType colBinning{{CfgVtxBins, CfgMultBins }, true };
236+ BinningType colBinning{{cfgVtxBins, cfgMultBins }, true };
238237 Preslice<aod::KShortTracks> tracksPerCollisionV0 = aod::kshortpionpair::kshortpioneventId;
239238 Preslice<aod::PionTracks> tracksPerCollisionBach = aod::kshortpionpair::kshortpioneventId;
240239
241240 void processMixedData (EventCandidates const & collisions, aod::KShortTracks const & V0s, aod::PionTracks const & piontracks)
242241 {
243- for (auto & [collision1, collision2] : selfCombinations (colBinning, nEvtMixing, -1 , collisions, collisions)) {
242+ for (const auto & [collision1, collision2] : selfCombinations (colBinning, nEvtMixing, -1 , collisions, collisions)) {
244243 if (collision1.index () == collision2.index ()) {
245244 continue ;
246245 }
@@ -249,7 +248,7 @@ struct cksspinalignder {
249248
250249 auto groupV0 = V0s.sliceBy (tracksPerCollisionV0, collision1.index ());
251250 auto groupPion = piontracks.sliceBy (tracksPerCollisionBach, collision2.index ());
252- for (auto & [t1, t2] : soa::combinations (o2::soa::CombinationsFullIndexPolicy (groupV0, groupPion))) {
251+ for (const auto & [t1, t2] : soa::combinations (o2::soa::CombinationsFullIndexPolicy (groupV0, groupPion))) {
253252 if (!selectionV0 (t1))
254253 continue ;
255254 auto [kshortPtmix, kshortEtamix, kshortPhimix] = computePtEtaPhi (t1.kShortPx (), t1.kShortPy (), t1.kShortPz ());
0 commit comments