Skip to content

Commit f44f588

Browse files
committed
Merge branch 'master' into HFTagging
2 parents 111dae5 + c7cc561 commit f44f588

File tree

280 files changed

+20943
-9493
lines changed

Some content is hidden

Large Commits have some content hidden by default. Use the searchbox below for content that may be hidden.

280 files changed

+20943
-9493
lines changed

ALICE3/Core/FastTracker.cxx

Lines changed: 29 additions & 1 deletion
Original file line numberDiff line numberDiff line change
@@ -144,6 +144,34 @@ void FastTracker::AddSiliconALICE3v2(std::vector<float> pixelResolution)
144144
AddLayer("B10", 80., 250, x0OT, xrhoOT, resRPhiOT, resZOT, eff, 1);
145145
}
146146

147+
void FastTracker::AddSiliconALICE3(std::vector<float> pixelResolution)
148+
{
149+
float x0IT = 0.001; // 0.1%
150+
float x0OT = 0.01; // 1.0%
151+
float xrhoIB = 2.3292e-02; // 100 mum Si
152+
float xrhoOT = 2.3292e-01; // 1000 mum Si
153+
float eff = 1.00;
154+
155+
float resRPhiIT = pixelResolution[0];
156+
float resZIT = pixelResolution[1];
157+
float resRPhiOT = pixelResolution[2];
158+
float resZOT = pixelResolution[3];
159+
160+
AddLayer("bpipe0", 0.48, 250, 0.00042, 2.772e-02, 0.0f, 0.0f, 0.0f, 0); // 150 mum Be
161+
AddLayer("B00", 0.5, 250, x0IT, xrhoIB, resRPhiIT, resZIT, eff, 1);
162+
AddLayer("B01", 1.2, 250, x0IT, xrhoIB, resRPhiIT, resZIT, eff, 1);
163+
AddLayer("B02", 2.5, 250, x0IT, xrhoIB, resRPhiIT, resZIT, eff, 1);
164+
AddLayer("bpipe1", 3.7, 250, 0.0014, 9.24e-02, 0.0f, 0.0f, 0.0f, 0); // 500 mum Be
165+
AddLayer("B03", 7., 250, x0OT, xrhoOT, resRPhiOT, resZOT, eff, 1);
166+
AddLayer("B04", 9., 250, x0OT, xrhoOT, resRPhiOT, resZOT, eff, 1);
167+
AddLayer("B05", 12., 250, x0OT, xrhoOT, resRPhiOT, resZOT, eff, 1);
168+
AddLayer("B06", 20., 250, x0OT, xrhoOT, resRPhiOT, resZOT, eff, 1);
169+
AddLayer("B07", 30., 250, x0OT, xrhoOT, resRPhiOT, resZOT, eff, 1);
170+
AddLayer("B08", 45., 250, x0OT, xrhoOT, resRPhiOT, resZOT, eff, 1);
171+
AddLayer("B09", 60., 250, x0OT, xrhoOT, resRPhiOT, resZOT, eff, 1);
172+
AddLayer("B10", 80., 250, x0OT, xrhoOT, resRPhiOT, resZOT, eff, 1);
173+
}
174+
147175
void FastTracker::AddTPC(float phiResMean, float zResMean)
148176
{
149177
LOG(info) << " Adding standard time projection chamber";
@@ -307,7 +335,7 @@ int FastTracker::FastTrack(o2::track::TrackParCov inputTrack, o2::track::TrackPa
307335
break;
308336
}
309337
}
310-
if (firstActiveLayer <= 0) {
338+
if (firstActiveLayer < 0) {
311339
LOG(fatal) << "No active layers found in FastTracker, check layer setup";
312340
return -2; // no active layers
313341
}

ALICE3/Core/FastTracker.h

Lines changed: 3 additions & 0 deletions
Original file line numberDiff line numberDiff line change
@@ -42,9 +42,11 @@ class FastTracker
4242
// Layer and layer configuration
4343
void AddLayer(TString name, float r, float z, float x0, float xrho, float resRPhi = 0.0f, float resZ = 0.0f, float eff = 0.0f, int type = 0);
4444
DetLayer GetLayer(const int layer, bool ignoreBarrelLayers = true) const;
45+
std::vector<DetLayer> GetLayers() const { return layers; }
4546
int GetLayerIndex(const std::string& name) const;
4647
size_t GetNLayers() const { return layers.size(); }
4748
bool IsLayerInert(const int layer) const { return layers[layer].isInert(); }
49+
void ClearLayers() { layers.clear(); }
4850
void SetRadiationLength(const std::string layerName, float x0) { layers[GetLayerIndex(layerName)].setRadiationLength(x0); }
4951
void SetRadius(const std::string layerName, float r) { layers[GetLayerIndex(layerName)].setRadius(r); }
5052
void SetResolutionRPhi(const std::string layerName, float resRPhi) { layers[GetLayerIndex(layerName)].setResolutionRPhi(resRPhi); }
@@ -57,6 +59,7 @@ class FastTracker
5759

5860
void AddSiliconALICE3v4(std::vector<float> pixelResolution);
5961
void AddSiliconALICE3v2(std::vector<float> pixelResolution);
62+
void AddSiliconALICE3(std::vector<float> pixelResolution);
6063
void AddTPC(float phiResMean, float zResMean);
6164

6265
void Print();

ALICE3/TableProducer/alice3-multicharmTable.cxx

Lines changed: 17 additions & 2 deletions
Original file line numberDiff line numberDiff line change
@@ -162,7 +162,8 @@ struct alice3multicharmTable {
162162
Partition<alice3tracks> tracksPiFromXiC =
163163
((aod::a3DecayMap::decayMap & trackSelectionPiFromXiC) == trackSelectionPiFromXiC) && aod::track::signed1Pt > 0.0f && 1.0f / nabs(aod::track::signed1Pt) > minPiCPt&& nabs(aod::track::dcaXY) > piFromXiC_dcaXYconstant + piFromXiC_dcaXYpTdep* nabs(aod::track::signed1Pt) && nabs(aod::track::dcaZ) > piFromXiC_dcaZconstant + piFromXiC_dcaZpTdep* nabs(aod::track::signed1Pt);
164164

165-
Partition<alice3tracks> tracksPiFromXiCC = ((aod::a3DecayMap::decayMap & trackSelectionPiFromXiCC) == trackSelectionPiFromXiCC) && aod::track::signed1Pt > 0.0f && 1.0f / nabs(aod::track::signed1Pt) > minPiCCPt&& nabs(aod::track::dcaXY) > piFromXiCC_dcaXYconstant + piFromXiCC_dcaXYpTdep* nabs(aod::track::signed1Pt);
165+
Partition<alice3tracks> tracksPiFromXiCC =
166+
((aod::a3DecayMap::decayMap & trackSelectionPiFromXiCC) == trackSelectionPiFromXiCC) && aod::track::signed1Pt > 0.0f && 1.0f / nabs(aod::track::signed1Pt) > minPiCCPt&& nabs(aod::track::dcaXY) > piFromXiCC_dcaXYconstant + piFromXiCC_dcaXYpTdep* nabs(aod::track::signed1Pt) && nabs(aod::track::dcaZ) > piFromXiCC_dcaZconstant + piFromXiCC_dcaZpTdep* nabs(aod::track::signed1Pt);
166167

167168
// Helper struct to pass candidate information
168169
struct {
@@ -446,6 +447,13 @@ struct alice3multicharmTable {
446447
histos.add("hPi2cPt", "hPi2cPt", kTH1D, {axisPt});
447448
histos.add("hPiccPt", "hPiccPt", kTH1D, {axisPt});
448449

450+
histos.add("hPi1cDCAxy", "hPi1cDCAxy", kTH1D, {axisDCA});
451+
histos.add("hPi1cDCAz", "hPi1cDCAz", kTH1D, {axisDCA});
452+
histos.add("hPi2cDCAxy", "hPi2cDCAxy", kTH1D, {axisDCA});
453+
histos.add("hPi2cDCAz", "hPi2cDCAz", kTH1D, {axisDCA});
454+
histos.add("hPiccDCAxy", "hPiccDCAxy", kTH1D, {axisDCA});
455+
histos.add("hPiccDCAz", "hPiccDCAz", kTH1D, {axisDCA});
456+
449457
histos.add("hMinXiDecayRadius", "hMinXiDecayRadius", kTH1D, {axisRadius2DXi});
450458
histos.add("hMinXiCDecayRadius", "hMinXiCDecayRadius", kTH1D, {axisRadius});
451459
histos.add("hMinXiCCDecayRadius", "hMinXiCCDecayRadius", kTH1D, {axisRadius});
@@ -723,7 +731,7 @@ struct alice3multicharmTable {
723731
histos.fill(HIST("hDCAxyXiCC"), std::fabs(xiccdcaXY * 1e+4));
724732
histos.fill(HIST("hDCAzXiCC"), std::fabs(xiccdcaZ * 1e+4));
725733

726-
if (std::fabs(thisXiCcandidate.eta) > xiccMaxEta)
734+
if (std::fabs(thisXiCCcandidate.eta) > xiccMaxEta)
727735
continue; // not in central barrel
728736

729737
histos.fill(HIST("hCharmBuilding"), 3.0f);
@@ -780,6 +788,13 @@ struct alice3multicharmTable {
780788
piFromLa.pt(), piFromLa.eta(),
781789
piFromLa.dcaXY(), piFromLa.dcaZ(),
782790
pi1c.eta(), pi2c.eta(), picc.eta());
791+
792+
histos.fill(HIST("hPi1cDCAxy"), std::abs(pi1c.dcaXY() * 1e+4));
793+
histos.fill(HIST("hPi1cDCAz"), std::abs(pi1c.dcaZ() * 1e+4));
794+
histos.fill(HIST("hPi2cDCAxy"), std::abs(pi2c.dcaXY() * 1e+4));
795+
histos.fill(HIST("hPi2cDCAz"), std::abs(pi2c.dcaZ() * 1e+4));
796+
histos.fill(HIST("hPiccDCAxy"), std::abs(picc.dcaXY() * 1e+4));
797+
histos.fill(HIST("hPiccDCAz"), std::abs(picc.dcaZ() * 1e+4));
783798
}
784799
}
785800
histos.fill(HIST("hCombinationsXiCC"), nCombinationsCC);

ALICE3/Tasks/alice3-multicharm.cxx

Lines changed: 6 additions & 6 deletions
Original file line numberDiff line numberDiff line change
@@ -255,7 +255,7 @@ struct alice3multicharm {
255255
histos.add("hBDTScoreVsXiccPt", "hBDTScoreVsXiccPt", kTH2D, {axisPt, axisBDTScore});
256256
histos.add("h3dBDTScore", "h3dBDTScore", kTH3D, {axisPt, axisXiccMass, axisBDTScore});
257257
for (const auto& score : bdt.requiredScores.value) {
258-
histPath = std::format("MLQA/RequiredBDTScore_{}/", static_cast<int>(score * 100));
258+
histPath = std::format("MLQA/RequiredBDTScore_{}/", static_cast<int>(score * 10000));
259259
histPointers.insert({histPath + "hDCAXicDaughters", histos.add((histPath + "hDCAXicDaughters").c_str(), "hDCAXicDaughters", {kTH1D, {{axisDcaDaughters}}})});
260260
histPointers.insert({histPath + "hDCAXiccDaughters", histos.add((histPath + "hDCAXiccDaughters").c_str(), "hDCAXiccDaughters", {kTH1D, {{axisDcaDaughters}}})});
261261
histPointers.insert({histPath + "hDCAxyXi", histos.add((histPath + "hDCAxyXi").c_str(), "hDCAxyXi", {kTH1D, {{axisDCA}}})});
@@ -329,7 +329,7 @@ struct alice3multicharm {
329329

330330
for (const auto& requiredScore : bdt.requiredScores.value) {
331331
if (bdtScore > requiredScore) {
332-
histPath = std::format("MLQA/RequiredBDTScore_{}/", static_cast<int>(requiredScore * 100));
332+
histPath = std::format("MLQA/RequiredBDTScore_{}/", static_cast<int>(requiredScore * 10000));
333333
getHist(TH1, histPath + "hDCAXicDaughters")->Fill(xiccCand.xicDauDCA() * 1e+4);
334334
getHist(TH1, histPath + "hDCAXiccDaughters")->Fill(xiccCand.xiccDauDCA() * 1e+4);
335335
getHist(TH1, histPath + "hDCAxyXi")->Fill(std::fabs(xiccCand.xiDCAxy() * 1e+4));
@@ -349,13 +349,13 @@ struct alice3multicharm {
349349
getHist(TH1, histPath + "hPi2cDCAz")->Fill(xiccCand.pi2cDCAz() * 1e+4);
350350
getHist(TH1, histPath + "hPiccDCAxy")->Fill(xiccCand.piccDCAxy() * 1e+4);
351351
getHist(TH1, histPath + "hPiccDCAz")->Fill(xiccCand.piccDCAz() * 1e+4);
352-
getHist(TH1, histPath + "hPi1cDCAz")->Fill(xiccCand.pi1cPt());
353-
getHist(TH1, histPath + "hPi2cDCAz")->Fill(xiccCand.pi2cPt());
354-
getHist(TH1, histPath + "hPiccDCAz")->Fill(xiccCand.piccPt());
352+
getHist(TH1, histPath + "hPi1cPt")->Fill(xiccCand.pi1cPt());
353+
getHist(TH1, histPath + "hPi2cPt")->Fill(xiccCand.pi2cPt());
354+
getHist(TH1, histPath + "hPiccPt")->Fill(xiccCand.piccPt());
355355
getHist(TH1, histPath + "hXiccMass")->Fill(xiccCand.xiccMass());
356356
getHist(TH1, histPath + "hXicMass")->Fill(xiccCand.xicMass());
357-
getHist(TH1, histPath + "hXiccPt")->Fill(xiccCand.xiccPt());
358357
getHist(TH1, histPath + "hXicPt")->Fill(xiccCand.xicPt());
358+
getHist(TH1, histPath + "hXiccPt")->Fill(xiccCand.xiccPt());
359359
getHist(TH3, histPath + "h3dXicc")->Fill(xiccCand.xiccPt(), xiccCand.xiccEta(), xiccCand.xiccMass());
360360
}
361361
}

Common/CCDB/TriggerAliases.h

Lines changed: 3 additions & 2 deletions
Original file line numberDiff line numberDiff line change
@@ -12,10 +12,11 @@
1212
#ifndef COMMON_CCDB_TRIGGERALIASES_H_
1313
#define COMMON_CCDB_TRIGGERALIASES_H_
1414

15+
#include <Rtypes.h>
16+
1517
#include <cstdint>
1618
#include <map>
1719
#include <string>
18-
#include <Rtypes.h>
1920

2021
enum triggerAliases {
2122
kINT7 = 0,
@@ -57,7 +58,7 @@ class TriggerAliases
5758
TriggerAliases() = default;
5859
~TriggerAliases() = default;
5960

60-
void AddAlias(uint32_t aliasId, std::string classNames) { mAliasToClassNames[aliasId] = classNames; }
61+
void AddAlias(uint32_t aliasId, std::string const& classNames) { mAliasToClassNames[aliasId] = classNames; }
6162
void AddClassIdToAlias(uint32_t aliasId, int classId);
6263
const std::map<uint32_t, std::string>& GetAliasToClassNamesMap() const { return mAliasToClassNames; }
6364
const std::map<uint32_t, ULong64_t>& GetAliasToTriggerMaskMap() const { return mAliasToTriggerMask; }

Common/CCDB/ctpRateFetcher.cxx

Lines changed: 8 additions & 7 deletions
Original file line numberDiff line numberDiff line change
@@ -11,18 +11,19 @@
1111

1212
#include "ctpRateFetcher.h"
1313

14+
#include <CCDB/BasicCCDBManager.h>
15+
#include <CommonConstants/LHCConstants.h>
16+
#include <DataFormatsCTP/Configuration.h>
17+
#include <DataFormatsCTP/Scalers.h>
18+
#include <DataFormatsParameters/GRPLHCIFData.h>
19+
1420
#include <map>
21+
#include <string>
1522
#include <vector>
1623

17-
#include "CommonConstants/LHCConstants.h"
18-
#include "DataFormatsCTP/Configuration.h"
19-
#include "DataFormatsCTP/Scalers.h"
20-
#include "DataFormatsParameters/GRPLHCIFData.h"
21-
#include "CCDB/BasicCCDBManager.h"
22-
2324
namespace o2
2425
{
25-
double ctpRateFetcher::fetch(o2::ccdb::BasicCCDBManager* ccdb, uint64_t timeStamp, int runNumber, std::string sourceName, bool fCrashOnNull)
26+
double ctpRateFetcher::fetch(o2::ccdb::BasicCCDBManager* ccdb, uint64_t timeStamp, int runNumber, const std::string& sourceName, bool fCrashOnNull)
2627
{
2728
setupRun(runNumber, ccdb, timeStamp);
2829
if (sourceName.find("ZNC") != std::string::npos) {

Common/CCDB/ctpRateFetcher.h

Lines changed: 3 additions & 3 deletions
Original file line numberDiff line numberDiff line change
@@ -12,9 +12,9 @@
1212
#ifndef COMMON_CCDB_CTPRATEFETCHER_H_
1313
#define COMMON_CCDB_CTPRATEFETCHER_H_
1414

15-
#include <string>
15+
#include <CCDB/BasicCCDBManager.h>
1616

17-
#include "CCDB/BasicCCDBManager.h"
17+
#include <string>
1818

1919
namespace o2
2020
{
@@ -34,7 +34,7 @@ class ctpRateFetcher
3434
{
3535
public:
3636
ctpRateFetcher() = default;
37-
double fetch(o2::ccdb::BasicCCDBManager* ccdb, uint64_t timeStamp, int runNumber, std::string sourceName, bool fCrashOnNull = true);
37+
double fetch(o2::ccdb::BasicCCDBManager* ccdb, uint64_t timeStamp, int runNumber, const std::string& sourceName, bool fCrashOnNull = true);
3838

3939
void setManualCleanup(bool manualCleanup = true) { mManualCleanup = manualCleanup; }
4040

Common/Core/CollisionAssociation.h

Lines changed: 25 additions & 25 deletions
Original file line numberDiff line numberDiff line change
@@ -20,13 +20,13 @@
2020
#ifndef COMMON_CORE_COLLISIONASSOCIATION_H_
2121
#define COMMON_CORE_COLLISIONASSOCIATION_H_
2222

23-
#include <vector>
23+
#include <CommonConstants/LHCConstants.h>
24+
#include <Framework/ASoAHelpers.h>
25+
#include <Framework/AnalysisDataModel.h>
26+
2427
#include <memory>
2528
#include <utility>
26-
27-
#include "CommonConstants/LHCConstants.h"
28-
#include "Framework/AnalysisDataModel.h"
29-
#include "Framework/ASoAHelpers.h"
29+
#include <vector>
3030

3131
namespace o2::aod
3232
{
@@ -191,58 +191,58 @@ class CollisionAssociation
191191
for (auto& iterationWindow : trackIterationWindows) {
192192
bool iteratorMoved = false;
193193
const bool isAssignedTrackWindow = (iterationWindow.first != iterationWindow.second) ? iterationWindow.first.has_collision() : false;
194-
for (auto track = iterationWindow.first; track != iterationWindow.second; ++track) {
195-
int64_t trackBC = globalBC[track.filteredIndex()];
194+
for (auto trackInWindow = iterationWindow.first; trackInWindow != iterationWindow.second; ++trackInWindow) {
195+
int64_t trackBC = globalBC[trackInWindow.filteredIndex()];
196196
if (trackBC < 0) {
197197
continue;
198198
}
199199

200200
// Optimization to avoid looping over the full track list each time. This builds on that tracks are sorted by BCs (which they should be because collisions are sorted by BCs)
201-
const int64_t bcOffset = trackBC - (int64_t)collBC;
201+
const int64_t bcOffset = trackBC - static_cast<int64_t>(collBC);
202202
if constexpr (isCentralBarrel) {
203203
// only for blocks with collision association
204204
if (isAssignedTrackWindow) {
205205
constexpr int margin = 200;
206206
if (!iteratorMoved && bcOffset > -bcOffsetMax - margin) {
207-
iterationWindow.first.setCursor(track.filteredIndex());
207+
iterationWindow.first.setCursor(trackInWindow.filteredIndex());
208208
iteratorMoved = true;
209-
LOGP(debug, "Moving iterator begin {}", track.filteredIndex());
209+
LOGP(debug, "Moving iterator begin {}", trackInWindow.filteredIndex());
210210
} else if (bcOffset > bcOffsetMax + margin) {
211-
LOGP(debug, "Stopping iterator {}", track.filteredIndex());
211+
LOGP(debug, "Stopping iterator {}", trackInWindow.filteredIndex());
212212
break;
213213
}
214214
}
215215
}
216216

217-
int64_t bcOffsetWindow = trackBCCache[track.filteredIndex()] - (int64_t)collBC;
217+
int64_t bcOffsetWindow = trackBCCache[trackInWindow.filteredIndex()] - static_cast<int64_t>(collBC);
218218
if (std::abs(bcOffsetWindow) > bcOffsetMax) {
219219
continue;
220220
}
221221

222222
float trackTime = 0;
223223
float trackTimeRes = 0;
224224
if constexpr (isCentralBarrel) {
225-
if (mUsePvAssociation && track.isPVContributor()) {
226-
trackTime = track.collision().collisionTime(); // if PV contributor, we assume the time to be the one of the collision
227-
trackTimeRes = o2::constants::lhc::LHCBunchSpacingNS; // 1 BC
225+
if (mUsePvAssociation && trackInWindow.isPVContributor()) {
226+
trackTime = trackInWindow.collision().collisionTime(); // if PV contributor, we assume the time to be the one of the collision
227+
trackTimeRes = o2::constants::lhc::LHCBunchSpacingNS; // 1 BC
228228
} else {
229-
trackTime = track.trackTime();
230-
trackTimeRes = track.trackTimeRes();
229+
trackTime = trackInWindow.trackTime();
230+
trackTimeRes = trackInWindow.trackTimeRes();
231231
}
232232
} else {
233-
trackTime = track.trackTime();
234-
trackTimeRes = track.trackTimeRes();
233+
trackTime = trackInWindow.trackTime();
234+
trackTimeRes = trackInWindow.trackTimeRes();
235235
}
236236

237237
const float deltaTime = trackTime - collTime + bcOffset * o2::constants::lhc::LHCBunchSpacingNS;
238238
float sigmaTimeRes2 = collTimeRes2 + trackTimeRes * trackTimeRes;
239-
LOGP(debug, "collision time={}, collision time res={}, track time={}, track time res={}, bc collision={}, bc track={}, delta time={}", collTime, collision.collisionTimeRes(), track.trackTime(), track.trackTimeRes(), collBC, trackBC, deltaTime);
239+
LOGP(debug, "collision time={}, collision time res={}, track time={}, track time res={}, bc collision={}, bc track={}, delta time={}", collTime, collision.collisionTimeRes(), trackInWindow.trackTime(), trackInWindow.trackTimeRes(), collBC, trackBC, deltaTime);
240240

241241
float thresholdTime = 0.;
242242
if constexpr (isCentralBarrel) {
243-
if (mUsePvAssociation && track.isPVContributor()) {
243+
if (mUsePvAssociation && trackInWindow.isPVContributor()) {
244244
thresholdTime = trackTimeRes;
245-
} else if (TESTBIT(track.flags(), o2::aod::track::TrackTimeResIsRange)) {
245+
} else if (TESTBIT(trackInWindow.flags(), o2::aod::track::TrackTimeResIsRange)) {
246246
// the track time resolution is a range, not a gaussian resolution
247247
thresholdTime = trackTimeRes + mNumSigmaForTimeCompat * std::sqrt(collTimeRes2) + mTimeMargin;
248248
} else {
@@ -262,7 +262,7 @@ class CollisionAssociation
262262

263263
if (std::abs(deltaTime) < thresholdTime) {
264264
const auto collIdx = collision.globalIndex();
265-
const auto trackIdx = track.globalIndex();
265+
const auto trackIdx = trackInWindow.globalIndex();
266266
LOGP(debug, "Filling track id {} for coll id {}", trackIdx, collIdx);
267267
association(collIdx, trackIdx);
268268
if (mFillTableOfCollIdsPerTrack) {
@@ -278,9 +278,9 @@ class CollisionAssociation
278278
// create reverse index track to collisions if enabled
279279
if (mFillTableOfCollIdsPerTrack) {
280280
std::vector<int> empty{};
281-
for (const auto& track : tracksUnfiltered) {
281+
for (const auto& trackUnfiltered : tracksUnfiltered) {
282282

283-
const auto trackId = track.globalIndex();
283+
const auto trackId = trackUnfiltered.globalIndex();
284284
if (collsPerTrack[trackId] == nullptr) {
285285
reverseIndices(empty);
286286
} else {

Common/Core/CollisionTypeHelper.cxx

Lines changed: 2 additions & 3 deletions
Original file line numberDiff line numberDiff line change
@@ -17,9 +17,8 @@
1717

1818
#include "Common/Core/CollisionTypeHelper.h"
1919

20-
#include "DataFormatsParameters/GRPLHCIFData.h"
21-
22-
#include <fairlogger/Logger.h>
20+
#include <DataFormatsParameters/GRPLHCIFData.h>
21+
#include <Framework/Logger.h>
2322

2423
#include <string>
2524

Common/Core/EventPlaneHelper.cxx

Lines changed: 11 additions & 10 deletions
Original file line numberDiff line numberDiff line change
@@ -20,13 +20,14 @@
2020

2121
#include "Common/Core/EventPlaneHelper.h"
2222

23+
#include <TMath.h>
24+
#include <TVector3.h>
25+
2326
#include <algorithm>
27+
#include <cstdio>
2428
#include <iterator>
25-
#include <vector>
2629
#include <memory>
27-
28-
#include "TMath.h"
29-
#include "TVector3.h"
30+
#include <vector>
3031

3132
double EventPlaneHelper::GetPhiFV0(int chno, o2::fv0::Geometry* fv0geom)
3233
{
@@ -37,8 +38,8 @@ double EventPlaneHelper::GetPhiFV0(int chno, o2::fv0::Geometry* fv0geom)
3738
float offsetX = 0.;
3839
float offsetY = 0.;
3940

40-
int cellsInLeft[] = {0, 1, 2, 3, 8, 9, 10, 11, 16, 17, 18, 19, 24, 25, 26, 27,
41-
32, 40, 33, 41, 34, 42, 35, 43};
41+
const int cellsInLeft[] = {0, 1, 2, 3, 8, 9, 10, 11, 16, 17, 18, 19, 24, 25, 26, 27,
42+
32, 40, 33, 41, 34, 42, 35, 43};
4243
bool isChnoInLeft = std::find(std::begin(cellsInLeft), std::end(cellsInLeft), chno) != std::end(cellsInLeft);
4344

4445
if (isChnoInLeft) {
@@ -131,10 +132,10 @@ void EventPlaneHelper::DoCorrections(float& qx, float& qy,
131132
qy = (qy - corrections[2] * qx) / (1.0 - corrections[3] * corrections[2]);
132133

133134
// Rescaling of the Qx-Qy into a circle.
134-
if (fabs(corrections[4]) > 1e-8) {
135+
if (std::fabs(corrections[4]) > 1e-8) {
135136
qx /= corrections[4];
136137
}
137-
if (fabs(corrections[5]) > 1e-8) {
138+
if (std::fabs(corrections[5]) > 1e-8) {
138139
qy /= corrections[5];
139140
}
140141
}
@@ -153,9 +154,9 @@ void EventPlaneHelper::DoTwist(float& qx, float& qy, float lp, float lm)
153154

154155
void EventPlaneHelper::DoRescale(float& qx, float& qy, float ap, float am)
155156
{
156-
if (fabs(ap) > 1e-8)
157+
if (std::fabs(ap) > 1e-8)
157158
qx /= ap;
158-
if (fabs(am) > 1e-8)
159+
if (std::fabs(am) > 1e-8)
159160
qy /= am;
160161
}
161162

0 commit comments

Comments
 (0)