Skip to content

Commit e2c8d48

Browse files
author
Aswathy Menon
committed
Changes after PR type 12 sept 2026
1 parent f19a6ee commit e2c8d48

6 files changed

Lines changed: 85 additions & 103 deletions

File tree

‎PWGDQ/Core/HistogramsLibrary.cxx‎

Lines changed: 29 additions & 0 deletions
Original file line numberDiff line numberDiff line change
@@ -121,6 +121,26 @@ void o2::aod::dqhistograms::DefineHistograms(HistogramManager* hm, const char* h
121121
if (subGroupStr.Contains("randomplane")) {
122122
hm->AddHistogram(histClass, "random plane angle", "random Psi", false, 100, -TMath::Pi() / 2, TMath::Pi() / 2, VarManager::kRandomPsi2);
123123
}
124+
if (subGroupStr.Contains("spherocity")) {
125+
hm->AddHistogram(histClass, "Spherocity", "Transverse Spherocity", false, 100, 0.0, 1.0, VarManager::kSpherocity);
126+
hm->AddHistogram(histClass, "SpherocityPtWeighted", "pT-weighted Spherocity", false, 100, 0.0, 1.0, VarManager::kSpherocityPtWeighted); // Spherocity vs event properties
127+
hm->AddHistogram(histClass, "Spherocity_VtxZ", "Spherocity vs Vtx Z", false, 100, 0.0, 1.0, VarManager::kSpherocity, 60, -15.0, 15.0, VarManager::kVtxZ);
128+
hm->AddHistogram(histClass, "MultITSTPC_VtxZ", "MultITSTPC vs Vtx Z", false, 200, 0.0, 200.0, VarManager::kMultNTracksITSTPC, 60, -15.0, 15.0, VarManager::kVtxZ);
129+
hm->AddHistogram(histClass, "Spherocity_VtxNContrib", "Spherocity vs Vtx n contributors", false, 100, 0.0, 1.0, VarManager::kSpherocity, 100, 0.0, 100.0, VarManager::kVtxNcontrib);
130+
hm->AddHistogram(histClass, "MultITSTPC_VtxNContrib", "MultITSTPC vs Vtx n contributors", false, 200, 0.0, 200.0, VarManager::kMultNTracksITSTPC, 100, 0.0, 100.0, VarManager::kVtxNcontrib);
131+
hm->AddHistogram(histClass, "Spherocity_MultITS", "Spherocity vs ITS only Multiplicity", false, 100, 0.0, 1.0, VarManager::kSpherocity, 250, 0.0, 250.0, VarManager::kMultNTracksITSOnly);
132+
hm->AddHistogram(histClass, "Spherocity_MultTPC", "Spherocity vs TPC only Multiplicity", false, 100, 0.0, 1.0, VarManager::kSpherocity, 250, 0.0, 5000.0, VarManager::kMultTPC);
133+
hm->AddHistogram(histClass, "Spherocity_MultITSTPC", "Spherocity vs ITS-TPC Multiplicity", false, 100, 0.0, 1.0, VarManager::kSpherocity, 250, 0.0, 250.0, VarManager::kMultNTracksITSTPC);
134+
hm->AddHistogram(histClass, "Spherocity_CentFT0C", "Spherocity vs Centrality FT0C", false, 100, 0.0, 1.0, VarManager::kSpherocity, 20, 0.0, 100.0, VarManager::kCentFT0C);
135+
hm->AddHistogram(histClass, "SpherocityPtWeighted_CentFT0C", "pT-weighted Spherocity vs Centrality FT0C", false, 100, 0.0, 1.0, VarManager::kSpherocityPtWeighted, 20, 0.0, 100.0, VarManager::kCentFT0C);
136+
hm->AddHistogram(histClass, "Spherocity_CentFT0C_3D", "Spherocity vs Centrality vs Mult", false, 50, 0.0, 1.0, VarManager::kSpherocity, 20, 0.0, 100.0, VarManager::kCentFT0C, 100, 0.0, 50000.0, VarManager::kMultTPC);
137+
138+
if (subGroupStr.Contains("mySpheroCent")) {
139+
hm->AddHistogram(histClass, "Spherocity_CentFT0C", "Spherocity vs Centrality FT0C", false, 100, 0.0, 1.0, VarManager::kSpherocity, 20, 0.0, 100.0, VarManager::kCentFT0C);
140+
hm->AddHistogram(histClass, "SpherocityPtWeighted_CentFT0C", "pT-weighted Spherocity vs Centrality FT0C", false, 100, 0.0, 1.0, VarManager::kSpherocityPtWeighted, 20, 0.0, 100.0, VarManager::kCentFT0C);
141+
hm->AddHistogram(histClass, "Spherocity_CentFT0C_3D", "Spherocity vs Centrality vs Mult", false, 50, 0.0, 1.0, VarManager::kSpherocity, 20, 0.0, 100.0, VarManager::kCentFT0C, 100, 0.0, 50000.0, VarManager::kMultTPC);
142+
}
143+
}
124144
if (subGroupStr.Contains("mult")) {
125145
if (subGroupStr.Contains("pp")) {
126146
hm->AddHistogram(histClass, "MultTPC", "MultTPC", false, 250, 0.0, 500.0, VarManager::kMultTPC);
@@ -1632,6 +1652,15 @@ void o2::aod::dqhistograms::DefineHistograms(HistogramManager* hm, const char* h
16321652
hm->AddHistogram(histClass, "Mass_Pt", "", false, 750, 0.0, 15.0, VarManager::kMass, 120, 0.0, 30.0, VarManager::kPt);
16331653
hm->AddHistogram(histClass, "Mass_Rapidity", "", false, 750, 0.0, 15.0, VarManager::kMass, 150, 2.5, 4.0, VarManager::kRap);
16341654
hm->AddHistogram(histClass, "Mass_Phi", "", false, 750, 0.0, 15.0, VarManager::kMass, 180, constants::math::PI, 2 * constants::math::PI, VarManager::kPhi);
1655+
hm->AddHistogram(histClass, "Mass", "", false, 750, 0.0, 5.0, VarManager::kMass);
1656+
if (subGroupStr.Contains("dimuon-spherocity")) {
1657+
hm->AddHistogram(histClass, "Spherocity", "Transverse Spherocity", false, 100, 0.0, 1.0, VarManager::kSpherocity);
1658+
hm->AddHistogram(histClass, "SpherocityPtWeighted", "pT-weighted Spherocity", false, 100, 0.0, 1.0, VarManager::kSpherocityPtWeighted);
1659+
hm->AddHistogram(histClass, "Mass_Spherocity", "", false, 750, 0.0, 5.0, VarManager::kMass, 100, 0.0, 1.0, VarManager::kSpherocity);
1660+
hm->AddHistogram(histClass, "Mass_SpherocityPtWeighted", "", false, 750, 0.0, 5.0, VarManager::kMass, 100, 0.0, 1.0, VarManager::kSpherocityPtWeighted);
1661+
hm->AddHistogram(histClass, "Mass_Pt_Spherocity", "", false, 750, 0.0, 5.0, VarManager::kMass, 120, 0.0, 30.0, VarManager::kPt, 100, 0.0, 1.0, VarManager::kSpherocity);
1662+
hm->AddHistogram(histClass, "Mass_Pt_SpherocityPtWeighted", "", false, 200, 0.0, 5.0, VarManager::kMass, 30, 0.0, 30.0, VarManager::kPt, 50, 0.0, 1.0, VarManager::kSpherocityPtWeighted);
1663+
}
16351664
if (subGroupStr.Contains("dimuon-tag-and-probe")) {
16361665
hm->AddHistogram(histClass, "Mass_PtProbe", "mass vs probe pT", false, 750, 0.0, 15.0, VarManager::kMass, 120, 0.0, 30.0, VarManager::kPt2); // Warning: in the tag-and-probe task t2 is always the probe
16371666
hm->AddHistogram(histClass, "Mass_EtaProbe", "mass vs probe eta", false, 750, 0.0, 15.0, VarManager::kMass, 120, 0.0, 30.0, VarManager::kEta2); // Warning: in the tag-and-probe task t2 is always the probe

‎PWGDQ/Core/VarManager.cxx‎

Lines changed: 4 additions & 4 deletions
Original file line numberDiff line numberDiff line change
@@ -816,9 +816,9 @@ void VarManager::SetDefaultVarNames()
816816
fgVariableUnits[kMultDimuonsME] = "";
817817
fgVariableNames[kCentFT0C] = "Centrality FT0C";
818818
fgVariableUnits[kCentFT0C] = "%";
819-
fgVariableNames[kSpherocity] = "Spherocity"; // transverse spherocity
819+
fgVariableNames[kSpherocity] = "Spherocity";
820820
fgVariableUnits[kSpherocity] = "";
821-
fgVariableNames[kSpherocityPtWeighted] = "SpherocityPtWeighted"; // pT-weighted transverse spherocity
821+
fgVariableNames[kSpherocityPtWeighted] = "SpherocityPtWeighted";
822822
fgVariableUnits[kSpherocityPtWeighted] = "";
823823
fgVariableNames[kMCEventGeneratorId] = "MC Generator ID";
824824
fgVariableNames[kMCEventSubGeneratorId] = "MC SubGenerator ID";
@@ -2100,8 +2100,8 @@ void VarManager::SetDefaultVarNames()
21002100
fgVarNamesMap["kVtxChi2"] = kVtxChi2;
21012101
fgVarNamesMap["kCentVZERO"] = kCentVZERO;
21022102
fgVarNamesMap["kCentFT0C"] = kCentFT0C;
2103-
fgVarNamesMap["kSpherocity"] = kSpherocity; // transverse spherocity
2104-
fgVarNamesMap["kSpherocityPtWeighted"] = kSpherocityPtWeighted; // pT-weighted transverse spherocity
2103+
fgVarNamesMap["kSpherocity"] = kSpherocity;
2104+
fgVarNamesMap["kSpherocityPtWeighted"] = kSpherocityPtWeighted;
21052105
fgVarNamesMap["kCentFT0A"] = kCentFT0A;
21062106
fgVarNamesMap["kCentFT0M"] = kCentFT0M;
21072107
fgVarNamesMap["kMultTPC"] = kMultTPC;

‎PWGDQ/Core/VarManager.h‎

Lines changed: 14 additions & 24 deletions
Original file line numberDiff line numberDiff line change
@@ -118,7 +118,7 @@ class VarManager : public TObject
118118
RapidityGapFilter = BIT(21),
119119
Fit = BIT(22),
120120
ReducedFit = BIT(23),
121-
ReducedEventSpherocity = BIT(24), // Achu-Changes
121+
ReducedEventSpherocity = BIT(24),
122122
Track = BIT(0),
123123
TrackCov = BIT(1),
124124
TrackExtra = BIT(2),
@@ -284,7 +284,7 @@ class VarManager : public TObject
284284
kMultNTracksITSTPC,
285285
kMultNTracksPVeta1,
286286
kMultNTracksPVetaHalf,
287-
kSpherocity, // transverse spherocity
287+
kSpherocity, // transverse spherocity
288288
kSpherocityPtWeighted, // pT-weighted transverse spherocity
289289
kTrackOccupancyInTimeRange,
290290
kFT0COccupancyInTimeRange,
@@ -1442,10 +1442,8 @@ class VarManager : public TObject
14421442
return RecoDecay::constrainAngle(psi1 - psi2, -o2::constants::math::PI / harmonic, harmonic);
14431443
}
14441444

1445-
// Spherocity function declaration--Start
14461445
template <typename T>
1447-
static float CalculateSpherocity(T const & tracks, float ptMin = 0.15, float ptMax = 10.0, float etaMin = -0.8, float etaMax = 0.8, int minMult = 10, bool usePtWeight = false);
1448-
//------------------------------------End
1446+
static float CalculateSpherocity(T const& tracks, float ptMin = 0.15, float ptMax = 10.0, float etaMin = -0.8, float etaMax = 0.8, int minMult = 10, bool usePtWeight = false);
14491447

14501448
template <typename T, typename T1>
14511449
static o2::dataformats::VertexBase RecalculatePrimaryVertex(T const& track0, T const& track1, const T1& collision);
@@ -8011,11 +8009,8 @@ void VarManager::FillTripletVertexingALICE3(C const& collision, T const& t1, T c
80118009
}
80128010
}
80138011
template <typename T>
8014-
float VarManager::CalculateSpherocity(T const& tracks, float ptMin, float ptMax, float etaMin, float etaMax, int minMult, bool usePtWeight)
8012+
float VarManager::CalculateSpherocity(T const& tracks, float ptMin, float ptMax, float etaMin, float etaMax, int minMult, bool usePtWeight)
80158013
{
8016-
LOGF(info, " >>> CalculateSpherocity CALLED: ptMin=%.2f, ptMax=%.2f, etaMin=%.2f, etaMax=%.2f, minMult=%d, usePtWeight=%s",
8017-
ptMin, ptMax, etaMin, etaMax, minMult, usePtWeight ? "true" : "false");
8018-
80198014
std::vector<float> pxNorm;
80208015
std::vector<float> pyNorm;
80218016
std::vector<float> ptValues;
@@ -8026,10 +8021,12 @@ float VarManager::CalculateSpherocity(T const& tracks, float ptMin, float ptMax,
80268021
for (const auto& track : tracks) {
80278022
float pt = track.pt();
80288023
float eta = track.eta();
8029-
8030-
if (pt < ptMin || pt > ptMax) continue;
8031-
if (eta < etaMin || eta > etaMax) continue;
8032-
8024+
8025+
if (pt < ptMin || pt > ptMax)
8026+
continue;
8027+
if (eta < etaMin || eta > etaMax)
8028+
continue;
8029+
80338030
float px = track.px();
80348031
float py = track.py();
80358032

@@ -8048,8 +8045,8 @@ float VarManager::CalculateSpherocity(T const& tracks, float ptMin, float ptMax,
80488045
return -1.0f;
80498046
}
80508047

8051-
// Minimize over azimuthal angle
8052-
float spherocity = 5.0f; //float prevSpherocity = 5.0f; const float convergenceThreshold = 1e-4f;
8048+
// Minimize over azimuthal angle, hence nSteps = 360 for 1 degree steps
8049+
float spherocity = 5.0f;
80538050
const int nSteps = 360;
80548051

80558052
for (int iStep = 0; iStep < nSteps; iStep++) {
@@ -8061,9 +8058,8 @@ float VarManager::CalculateSpherocity(T const& tracks, float ptMin, float ptMax,
80618058
float denominator = 0.0f;
80628059

80638060
for (int iTrk = 0; iTrk < nTracks; iTrk++) {
8064-
// Cross product magnitude: |p_T x n| = |px*ny - py*nx|
8065-
float crossProduct = std::abs(pxNorm[iTrk] * ny - pyNorm[iTrk] * nx);
8066-
8061+
float crossProduct = std::abs(pxNorm[iTrk] * ny - pyNorm[iTrk] * nx); // Cross product magnitude: |p_T x n| = |px*ny - py*nx|
8062+
80678063
if (usePtWeight) {
80688064
numerator += ptValues[iTrk] * crossProduct;
80698065
denominator += ptValues[iTrk];
@@ -8077,16 +8073,10 @@ float VarManager::CalculateSpherocity(T const& tracks, float ptMin, float ptMax,
80778073
if (sph < spherocity) {
80788074
spherocity = sph;
80798075
}
8080-
// // Early exit if converged
8081-
// if (iStep > 10 && std::abs(spherocity - prevSpherocity) < convergenceThreshold) {
8082-
// break;
8083-
// }
8084-
// prevSpherocity = spherocity;
80858076
}
80868077

80878078
spherocity *= (M_PI * M_PI / 4.0f);
80888079
return spherocity;
80898080
}
80908081

8091-
80928082
#endif // PWGDQ_CORE_VARMANAGER_H_

‎PWGDQ/DataModel/ReducedInfoTables.h‎

Lines changed: 5 additions & 17 deletions
Original file line numberDiff line numberDiff line change
@@ -69,8 +69,8 @@ DECLARE_SOA_COLUMN(NTPCoccupMeanTimeShortA, nTPCoccupMeanTimeShortA, float);
6969
DECLARE_SOA_COLUMN(NTPCoccupMeanTimeShortC, nTPCoccupMeanTimeShortC, float); //! TPC pileup mean time on C side (short time range)
7070
DECLARE_SOA_COLUMN(NTPCoccupMedianTimeShortA, nTPCoccupMedianTimeShortA, float); //! TPC pileup median time on A side (short time range)
7171
DECLARE_SOA_COLUMN(NTPCoccupMedianTimeShortC, nTPCoccupMedianTimeShortC, float); //! TPC pileup median time on C side (short time range)
72-
DECLARE_SOA_COLUMN(Spherocity, spherocity, float); //! transverse spherocity
73-
DECLARE_SOA_COLUMN(SpherocityPtWeighted, spherocityPtWeighted, float); //! pT-weighted transverse spherocity// namespace reducedevent
72+
DECLARE_SOA_COLUMN(Spherocity, spherocity, float); //! transverse spherocity
73+
DECLARE_SOA_COLUMN(SpherocityPtWeighted, spherocityPtWeighted, float); //! pT-weighted transverse spherocity
7474
DECLARE_SOA_COLUMN(DCAzBimodalityCoefficient, dcazBimodalityCoefficient, float); //! Bimodality coefficient of the DCAz distribution of the tracks in the event
7575
DECLARE_SOA_COLUMN(DCAzBimodalityCoefficientBinned, dcazBimodalityCoefficientBinned, float); //! Bimodality coefficient of the DCAz distribution of the tracks in the event, binned
7676
DECLARE_SOA_COLUMN(DCAzBimodalityCoefficientBinnedTrimmed1, dcazBimodalityCoefficientBinnedTrimmed1, float); //! Bimodality coefficient of the DCAz distribution of the tracks in the event, binned and trimmed 1
@@ -179,23 +179,11 @@ DECLARE_SOA_TABLE_VERSIONED(ReducedEventsExtended_001, "AOD", "REEXTENDED", 1, /
179179
mult::MultFDDA, mult::MultFDDC, mult::MultZNA, mult::MultZNC, mult::MultTracklets, mult::MultNTracksPV,
180180
cent::CentFT0C, cent::CentFT0A, cent::CentFT0M);
181181

182-
DECLARE_SOA_TABLE_VERSIONED(ReducedEventsExtended_002, "AOD", "REEXTENDED", 2, //! Extended event information
183-
bc::GlobalBC, evsel::Alias, evsel::Selection, timestamp::Timestamp, cent::CentRun2V0M,
184-
mult::MultTPC, mult::MultFV0A, mult::MultFV0C, mult::MultFT0A, mult::MultFT0C,
185-
mult::MultFDDA, mult::MultFDDC, mult::MultZNA, mult::MultZNC, mult::MultTracklets, mult::MultNTracksPV,
186-
cent::CentFT0C, cent::CentFT0A, cent::CentFT0M,
187-
reducedevent::Spherocity,
188-
reducedevent::SpherocityPtWeighted);
189-
190182
using ReducedEventsExtended = ReducedEventsExtended_001;
191-
// using ReducedEventsExtended = ReducedEventsExtended_002;
192-
193-
// // In ReducedInfoTables.h, add this after the ReducedEventsExtended declaration
194-
// DECLARE_SOA_TABLE(ReducedEventsSpherocity, "AOD", "RESPHEROCITY",
195-
// reducedevent::Spherocity,
196-
// reducedevent::SpherocityPtWeighted);
197183

198-
// using ReducedEventSphero = ReducedEventsSpherocity::iterator;
184+
DECLARE_SOA_TABLE(ReducedEventsSpherocity, "AOD", "RESPHEROCITY",
185+
reducedevent::Spherocity,
186+
reducedevent::SpherocityPtWeighted);
199187

200188
DECLARE_SOA_TABLE(ReducedEventsMultPV_000, "AOD", "REMULTPV", //! Multiplicity information for primary vertex
201189
mult::MultNTracksHasITS, mult::MultNTracksHasTPC, mult::MultNTracksHasTOF, mult::MultNTracksHasTRD,

‎PWGDQ/TableProducer/tableMaker_withAssoc.cxx‎

Lines changed: 3 additions & 34 deletions
Original file line numberDiff line numberDiff line change
@@ -222,7 +222,7 @@ struct TableMaker {
222222
// Q-vector related tables, to be filled only if the user selects the corresponding option
223223
Produces<ReducedEventsQvectorCentr> eventQvectorCentr;
224224
Produces<ReducedEventsQvectorCentrExtra> eventQvectorCentrExtra;
225-
Produces<ReducedEventsSpherocity> eventSpherocity; // Achu-Changes
225+
Produces<ReducedEventsSpherocity> eventSpherocity;
226226
} outTables;
227227

228228
OutputObj<THashList> fOutputList{"output"}; //! the histogram manager output list
@@ -1054,12 +1054,6 @@ struct TableMaker {
10541054
// NOTE: So far, collisions are filtered based on the user specified analysis cuts AND the filterPP or Zorro event filter.
10551055
// The collision-track associations which point to an event that is not selected for writing are discarded!
10561056

1057-
LOGF(info, "=== skimCollisions CALLED ===");
1058-
LOGF(info, "TEventFillMap = %u", TEventFillMap);
1059-
LOGF(info, "TTrackFillMap = %u", TTrackFillMap);
1060-
LOGF(info, "TTrackFillMap is zero: %s", TTrackFillMap == 0 ? "TRUE" : "FALSE");
1061-
LOGF(info, "Number of collisions = %d", static_cast<int>(collisions.size()));
1062-
10631057
VarManager::FillTimeFrame(collisions);
10641058
fCollIndexMap.clear();
10651059
int multTPC = -1.0;
@@ -1245,23 +1239,9 @@ struct TableMaker {
12451239
float spherocityPtWeighted = -1.0f;
12461240

12471241
if constexpr (static_cast<bool>(TTrackFillMap)) {
1248-
// DIAGNOSTIC: Check table sizes
1249-
LOGF(info, "=== skimCollisions DIAGNOSTIC ===");
1250-
LOGF(info, "Total collisions: %d", static_cast<int>(collisions.size()));
1251-
LOGF(info, "Total tracksBarrel: %d", static_cast<int>(tracksBarrel.size()));
1252-
LOGF(info, "Total trackAssocs: %d", static_cast<int>(trackAssocs.size()));
1253-
LOGF(info, "================================");
1254-
1255-
LOGF(info, ">>> Collision %d: TTrackFillMap check PASSED", collision.globalIndex());
1256-
12571242
auto groupedTrackIndices = trackAssocs.sliceBy(trackIndicesPerCollision, collision.globalIndex());
1258-
LOGF(info, ">>> Collision %d: groupedTrackIndices.size() = %d",
1259-
collision.globalIndex(), static_cast<int>(groupedTrackIndices.size()));
12601243

12611244
if (static_cast<int>(groupedTrackIndices.size()) >= fConfigSpherocity.fConfigSpherocityMinMult) {
1262-
LOGF(info, ">>> Collision %d: PASSED multiplicity check, collecting tracks...",
1263-
collision.globalIndex());
1264-
12651245
std::vector<typename std::decay_t<decltype(tracksBarrel.begin())>> trackVec;
12661246
trackVec.reserve(groupedTrackIndices.size());
12671247

@@ -1274,9 +1254,6 @@ struct TableMaker {
12741254
}
12751255
}
12761256

1277-
LOGF(info, ">>> Collision %d: Collected %d valid tracks",
1278-
collision.globalIndex(), validTracks);
1279-
12801257
if (validTracks >= fConfigSpherocity.fConfigSpherocityMinMult) {
12811258
// Calculate unweighted spherocity
12821259
spherocityValue = VarManager::CalculateSpherocity(
@@ -1297,21 +1274,13 @@ struct TableMaker {
12971274
fConfigSpherocity.fConfigSpherocityEtaMax,
12981275
fConfigSpherocity.fConfigSpherocityMinMult,
12991276
true);
1300-
1301-
LOGF(info, "Event %d: Spherocity = %.3f, pT-weighted = %.3f, nTracks = %d",
1302-
collision.globalIndex(), spherocityValue, spherocityPtWeighted, validTracks);
13031277
}
13041278
}
13051279
}
13061280

1307-
// CRITICAL: Fill VarManager::fgValues so histograms can be filled!
13081281
VarManager::fgValues[VarManager::kSpherocity] = spherocityValue;
13091282
VarManager::fgValues[VarManager::kSpherocityPtWeighted] = spherocityPtWeighted;
1310-
LOGF(info, "Filled VarManager: Spherocity = %.3f, pT-weighted = %.3f",
1311-
VarManager::fgValues[VarManager::kSpherocity],
1312-
VarManager::fgValues[VarManager::kSpherocityPtWeighted]);
13131283

1314-
// fHistMan->FillHistClass("Event_AfterCuts", VarManager::fgValues);
13151284
fHistMan->FillHistClass("Event_AfterCuts", dqtablemaker_helpers::varValues());
13161285

13171286
// create the event tables
@@ -1335,7 +1304,7 @@ struct TableMaker {
13351304
centFT0M = collision.centFT0M();
13361305
}
13371306
outTables.eventExtended(bc.globalBC(), collision.alias_raw(), collision.selection_raw(), bc.timestamp(), VarManager::fgValues[VarManager::kCentVZERO],
1338-
multTPC, multFV0A, multFV0C, multFT0A, multFT0C, multFDDA, multFDDC, multZNA, multZNC, multTracklets, multTracksPV, centFT0C, centFT0A, centFT0M, spherocityValue, spherocityPtWeighted);
1307+
multTPC, multFV0A, multFV0C, multFT0A, multFT0C, multFDDA, multFDDC, multZNA, multZNC, multTracklets, multTracksPV, centFT0C, centFT0A, centFT0M);
13391308
outTables.eventVtxCov(collision.covXX(), collision.covXY(), collision.covXZ(), collision.covYY(), collision.covYZ(), collision.covZZ(), collision.chi2());
13401309
outTables.eventInfo(collision.globalIndex());
13411310
outTables.eventSpherocity(spherocityValue, spherocityPtWeighted);
@@ -2005,7 +1974,7 @@ struct TableMaker {
20051974
outTables.event.reserve(collisions.size());
20061975
outTables.eventExtended.reserve(collisions.size());
20071976
outTables.eventVtxCov.reserve(collisions.size());
2008-
outTables.eventSpherocity.reserve(collisions.size()); // ADD THIS LINE//Achu-Changes
1977+
outTables.eventSpherocity.reserve(collisions.size());
20091978

20101979
skimCollisions<TEventFillMap, TTrackFillMap>(collisions, bcs, zdcs, ft0s, fv0as, fdds, trackAssocs, tracksBarrel);
20111980
if (fCollIndexMap.empty()) {

0 commit comments

Comments
 (0)