Skip to content

Commit bfff0fd

Browse files
arvindkhuntiaArvindiiti
arvindkhuntia
and
Arvindiiti
authored
[PWGJE] Added abs value for nSigmaTPC (#11022)
Co-authored-by: Arvind Khuntia <arvind.khuntia@cern.ch>
1 parent e212f8e commit bfff0fd

File tree

1 file changed

+42
-44
lines changed

1 file changed

+42
-44
lines changed

‎PWGJE/Tasks/nucleiInJets.cxx‎

Lines changed: 42 additions & 44 deletions
Original file line numberDiff line numberDiff line change
@@ -41,7 +41,6 @@
4141
#include "EventFiltering/ZorroSummary.h"
4242

4343
#include "ReconstructionDataFormats/Track.h"
44-
4544
#include "PWGLF/DataModel/LFParticleIdentification.h"
4645

4746
#include "PWGJE/Core/FastJetUtilities.h"
@@ -698,13 +697,13 @@ struct nucleiInJets {
698697
}
699698

700699
if (useTOFNsigmaPreSel && trk.hasTOF()) {
701-
if (std::fabs(trk.tofNSigmaPr()) < cfgnTPCPIDPrTOF)
700+
if (std::abs(trk.tofNSigmaPr()) < cfgnTPCPIDPrTOF)
702701
jetHist.fill(HIST("tracks/proton/h3PtVsProtonNSigmaTPCVsPtJet_jet"), trk.pt(), trk.tpcNSigmaPr(), jetPt);
703-
if (std::fabs(trk.tofNSigmaDe()) < cfgnTPCPIDDeTOF)
702+
if (std::abs(trk.tofNSigmaDe()) < cfgnTPCPIDDeTOF)
704703
jetHist.fill(HIST("tracks/deuteron/h3PtVsDeuteronNSigmaTPCVsPtJet_jet"), trk.pt(), trk.tpcNSigmaDe(), jetPt);
705-
if (std::fabs(trk.tofNSigmaHe()) < cfgnTPCPIDHeTOF)
704+
if (std::abs(trk.tofNSigmaHe()) < cfgnTPCPIDHeTOF)
706705
jetHist.fill(HIST("tracks/helium/h3PtVsHeliumNSigmaTPCVsPtJet_jet"), trk.pt(), trk.tpcNSigmaHe(), jetPt);
707-
if (std::fabs(trk.tofNSigmaTr()) < cfgnTPCPIDTrTOF)
706+
if (std::abs(trk.tofNSigmaTr()) < cfgnTPCPIDTrTOF)
708707
jetHist.fill(HIST("tracks/triton/h3PtVsTritonNSigmaTPCVsPtJet_jet"), trk.pt(), trk.tpcNSigmaTr(), jetPt);
709708
} else if (!useTOFNsigmaPreSel) {
710709
jetHist.fill(HIST("tracks/proton/h3PtVsProtonNSigmaTPCVsPtJet_jet"), trk.pt(), trk.tpcNSigmaPr(), jetPt);
@@ -720,7 +719,7 @@ struct nucleiInJets {
720719
jetHist.fill(HIST("tracks/deuteron/dca/after/hDCAxyVsPtDeuteron_jet"), trk.dcaXY(), trk.pt());
721720
jetHist.fill(HIST("tracks/deuteron/dca/after/hDCAzVsPtDeuteron_jet"), trk.dcaZ(), trk.pt());
722721
}
723-
if (cEnableTritonQA && std::abs(trk.tpcNSigmaHe()) < cfgnTPCPIDTr) {
722+
if (cEnableTritonQA && std::abs(trk.tpcNSigmaTr()) < cfgnTPCPIDTr) {
724723
jetHist.fill(HIST("tracks/triton/dca/after/hDCAxyVsPtTriton_jet"), trk.dcaXY(), trk.pt());
725724
jetHist.fill(HIST("tracks/triton/dca/after/hDCAzVsPtTriton_jet"), trk.dcaZ(), trk.pt());
726725
}
@@ -746,24 +745,24 @@ struct nucleiInJets {
746745
jetHist.fill(HIST("tracks/proton/h3TpcNsigmaTofNsigmaProtonVsPt_jet"), trk.tpcNSigmaPr(), trk.tofNSigmaPr(), trk.pt());
747746
jetHist.fill(HIST("tracks/deuteron/h3TpcNsigmaTofNsigmaDeuteronVsPt_jet"), trk.tpcNSigmaDe(), trk.tofNSigmaDe(), trk.pt());
748747
} else {
749-
if (trk.tpcNSigmaPr() < useTPCpreSel) {
748+
if (std::abs(trk.tpcNSigmaPr()) < cfgnTPCPIDPr) {
750749
jetHist.fill(HIST("tracks/proton/h2TOFmassProtonVsPt_jet"), massTOF, trk.pt());
751750
jetHist.fill(HIST("tracks/proton/h2TOFmass2ProtonVsPt_jet"), massTOF * massTOF - gMassProton * gMassProton, trk.pt());
752751
jetHist.fill(HIST("tracks/proton/h2TofNsigmaProtonVsPt_jet"), trk.tofNSigmaPr(), trk.pt());
753752
jetHist.fill(HIST("tracks/proton/h3TpcNsigmaTofNsigmaProtonVsPt_jet"), trk.tpcNSigmaPr(), trk.tofNSigmaPr(), trk.pt());
754753
}
755-
if (trk.tpcNSigmaDe() < useTPCpreSel) {
754+
if (std::abs(trk.tpcNSigmaDe()) < cfgnTPCPIDDe) {
756755
jetHist.fill(HIST("tracks/deuteron/h2TOFmassDeuteronVsPt_jet"), massTOF, trk.pt());
757756
jetHist.fill(HIST("tracks/deuteron/h2TOFmass2DeuteronVsPt_jet"), massTOF * massTOF - gMassDeuteron * gMassDeuteron, trk.pt());
758757
jetHist.fill(HIST("tracks/deuteron/h2TofNsigmaDeuteronVsPt_jet"), trk.tofNSigmaDe(), trk.pt());
759758
jetHist.fill(HIST("tracks/deuteron/h3TpcNsigmaTofNsigmaDeuteronVsPt_jet"), trk.tpcNSigmaDe(), trk.tofNSigmaDe(), trk.pt());
760759
}
761-
if (trk.tpcNSigmaTr() < useTPCpreSel) {
760+
if (std::abs(trk.tpcNSigmaTr()) < cfgnTPCPIDTr) {
762761
jetHist.fill(HIST("tracks/triton/h2TOFmassTritonVsPt_jet"), massTOF, trk.pt());
763762
jetHist.fill(HIST("tracks/triton/h2TOFmass2TritonVsPt_jet"), massTOF * massTOF - gMassTriton * gMassTriton, trk.pt());
764763
jetHist.fill(HIST("tracks/triton/h2TofNsigmaTritonVsPt_jet"), trk.tofNSigmaTr(), trk.pt());
765764
}
766-
if (trk.tpcNSigmaHe() < useTPCpreSel) {
765+
if (std::abs(trk.tpcNSigmaHe()) < cfgnTPCPIDHe) {
767766
jetHist.fill(HIST("tracks/helium/h2TOFmassHeliumVsPt_jet"), massTOF, trk.pt());
768767
jetHist.fill(HIST("tracks/helium/h2TOFmass2HeliumVsPt_jet"), massTOF * massTOF - gMassHelium * gMassHelium, trk.pt());
769768
jetHist.fill(HIST("tracks/helium/h2TofNsigmaHeliumVsPt_jet"), trk.tofNSigmaHe(), trk.pt());
@@ -807,7 +806,7 @@ struct nucleiInJets {
807806
jetHist.fill(HIST("tracks/antiDeuteron/dca/after/hDCAxyVsPtantiDeuteron_jet"), trk.dcaXY(), trk.pt());
808807
jetHist.fill(HIST("tracks/antiDeuteron/dca/after/hDCAzVsPtantiDeuteron_jet"), trk.dcaZ(), trk.pt());
809808
}
810-
if (cEnableHeliumQA && std::abs(trk.tpcNSigmaHe()) < cfgnTPCPIDHe) {
809+
if (cEnableHeliumQA && std::abs(trk.tpcNSigmaTr()) < cfgnTPCPIDTr) {
811810
jetHist.fill(HIST("tracks/antiTriton/dca/after/hDCAxyVsPtantiTriton_jet"), trk.dcaXY(), trk.pt());
812811
jetHist.fill(HIST("tracks/antiTriton/dca/after/hDCAzVsPtantiTriton_jet"), trk.dcaZ(), trk.pt());
813812
}
@@ -836,24 +835,24 @@ struct nucleiInJets {
836835
jetHist.fill(HIST("tracks/antiHelium/h2TOFmass2antiHeliumVsPt_jet"), massTOF * massTOF - gMassHelium * gMassHelium, trk.pt());
837836
jetHist.fill(HIST("tracks/antiHelium/h2TofNsigmaantiHeliumVsPt_jet"), trk.tofNSigmaHe(), trk.pt());
838837
} else {
839-
if (trk.tpcNSigmaPr() < useTPCpreSel) {
838+
if (std::abs(trk.tpcNSigmaPr()) < cfgnTPCPIDPr) {
840839
jetHist.fill(HIST("tracks/antiProton/h2TOFmassantiProtonVsPt_jet"), massTOF, trk.pt());
841840
jetHist.fill(HIST("tracks/antiProton/h2TOFmass2antiProtonVsPt_jet"), massTOF * massTOF - gMassProton * gMassProton, trk.pt());
842841
jetHist.fill(HIST("tracks/antiProton/h2TofNsigmaantiProtonVsPt_jet"), trk.tofNSigmaPr(), trk.pt());
843842
jetHist.fill(HIST("tracks/antiProton/h3TpcNsigmaTofNsigmaantiProtonVsPt_jet"), trk.tpcNSigmaPr(), trk.tofNSigmaPr(), trk.pt());
844843
}
845-
if (trk.tpcNSigmaDe() < useTPCpreSel) {
844+
if (std::abs(trk.tpcNSigmaDe()) < cfgnTPCPIDDe) {
846845
jetHist.fill(HIST("tracks/antiDeuteron/h2TOFmassantiDeuteronVsPt_jet"), massTOF, trk.pt());
847846
jetHist.fill(HIST("tracks/antiDeuteron/h2TOFmass2antiDeuteronVsPt_jet"), massTOF * massTOF - gMassDeuteron * gMassDeuteron, trk.pt());
848847
jetHist.fill(HIST("tracks/antiDeuteron/h2TofNsigmaantiDeuteronVsPt_jet"), trk.tofNSigmaDe(), trk.pt());
849848
jetHist.fill(HIST("tracks/antiDeuteron/h3TpcNsigmaTofNsigmaantiDeuteronVsPt_jet"), trk.tpcNSigmaDe(), trk.tofNSigmaDe(), trk.pt());
850849
}
851-
if (trk.tpcNSigmaTr() < useTPCpreSel) {
850+
if (std::abs(trk.tpcNSigmaTr()) < cfgnTPCPIDTr) {
852851
jetHist.fill(HIST("tracks/antiTriton/h2TOFmassantiTritonVsPt_jet"), massTOF, trk.pt());
853852
jetHist.fill(HIST("tracks/antiTriton/h2TOFmass2antiTritonVsPt_jet"), massTOF * massTOF - gMassTriton * gMassTriton, trk.pt());
854853
jetHist.fill(HIST("tracks/antiTriton/h2TofNsigmaantiTritonVsPt_jet"), trk.tofNSigmaTr(), trk.pt());
855854
}
856-
if (trk.tpcNSigmaHe() < useTPCpreSel) {
855+
if (std::abs(trk.tpcNSigmaHe()) < cfgnTPCPIDHe) {
857856
jetHist.fill(HIST("tracks/antiHelium/h2TOFmassantiHeliumVsPt_jet"), massTOF, trk.pt());
858857
jetHist.fill(HIST("tracks/antiHelium/h2TOFmass2antiHeliumVsPt_jet"), massTOF * massTOF - gMassHelium * gMassHelium, trk.pt());
859858
jetHist.fill(HIST("tracks/antiHelium/h2TofNsigmaantiHeliumVsPt_jet"), trk.tofNSigmaHe(), trk.pt());
@@ -864,10 +863,10 @@ struct nucleiInJets {
864863
jetHist.fill(HIST("tracks/antiPion/h2TofNsigmaantiPionVsPt_jet"), trk.tofNSigmaPi(), trk.pt());
865864
jetHist.fill(HIST("tracks/antiKaon/h2TofNsigmaantiKaonVsPt_jet"), trk.tofNSigmaKa(), trk.pt());
866865
} else {
867-
if (trk.tpcNSigmaPi() < useTPCpreSel) {
866+
if (std::abs(trk.tpcNSigmaPi()) < useTPCpreSel) {
868867
jetHist.fill(HIST("tracks/antiPion/h2TofNsigmaantiPionVsPt_jet"), trk.tofNSigmaPi(), trk.pt());
869868
}
870-
if (trk.tpcNSigmaKa() < useTPCpreSel) {
869+
if (std::abs(trk.tpcNSigmaKa()) < useTPCpreSel) {
871870
jetHist.fill(HIST("tracks/antiKaon/h2TofNsigmaantiKaonVsPt_jet"), trk.tofNSigmaKa(), trk.pt());
872871
}
873872
}
@@ -917,7 +916,7 @@ struct nucleiInJets {
917916
jetHist.fill(HIST("tracks/deuteron/dca/after/hDCAxyVsPtDeuteron"), trk.dcaXY(), trk.pt());
918917
jetHist.fill(HIST("tracks/deuteron/dca/after/hDCAzVsPtDeuteron"), trk.dcaZ(), trk.pt());
919918
}
920-
if (cEnableTritonQA && std::abs(trk.tpcNSigmaHe()) < cfgnTPCPIDTr) {
919+
if (cEnableTritonQA && std::abs(trk.tpcNSigmaTr()) < cfgnTPCPIDTr) {
921920
jetHist.fill(HIST("tracks/triton/dca/after/hDCAxyVsPtTriton"), trk.dcaXY(), trk.pt());
922921
jetHist.fill(HIST("tracks/triton/dca/after/hDCAzVsPtTriton"), trk.dcaZ(), trk.pt());
923922
}
@@ -943,28 +942,28 @@ struct nucleiInJets {
943942
jetHist.fill(HIST("tracks/helium/h2TOFmass2HeliumVsPt"), massTOF * massTOF - gMassHelium * gMassHelium, trk.pt());
944943
jetHist.fill(HIST("tracks/triton/h2TofNsigmaTritonVsPt"), trk.tofNSigmaTr(), trk.pt());
945944
} else {
946-
if (trk.tpcNSigmaPr() < useTPCpreSel) {
945+
if (std::abs(trk.tpcNSigmaPr()) < cfgnTPCPIDPr) {
947946
jetHist.fill(HIST("tracks/proton/h2TOFmassProtonVsPt"), massTOF, trk.pt());
948947
jetHist.fill(HIST("tracks/proton/h2TOFmass2ProtonVsPt"), massTOF * massTOF - gMassProton * gMassProton, trk.pt());
949948
jetHist.fill(HIST("tracks/proton/h2TofNsigmaProtonVsPt"), trk.tofNSigmaPr(), trk.pt());
950949
if (jetFlagPerpCone && isWithLeadingJet)
951950
jetHist.fill(HIST("tracks/perpCone/proton/h2TofNsigmaProtonVsPt"), trk.tofNSigmaPr(), trk.pt());
952951
}
953-
if (trk.tpcNSigmaDe() < useTPCpreSel) {
952+
if (std::abs(trk.tpcNSigmaDe()) < cfgnTPCPIDDe) {
954953
jetHist.fill(HIST("tracks/deuteron/h2TOFmassDeuteronVsPt"), massTOF, trk.pt());
955954
jetHist.fill(HIST("tracks/deuteron/h2TOFmass2DeuteronVsPt"), massTOF * massTOF - gMassDeuteron * gMassDeuteron, trk.pt());
956955
jetHist.fill(HIST("tracks/deuteron/h2TofNsigmaDeuteronVsPt"), trk.tofNSigmaDe(), trk.pt());
957956
if (jetFlagPerpCone && isWithLeadingJet)
958957
jetHist.fill(HIST("tracks/perpCone/deuteron/h2TofNsigmaDeuteronVsPt"), trk.tofNSigmaDe(), trk.pt());
959958
}
960-
if (trk.tpcNSigmaTr() < useTPCpreSel) {
959+
if (std::abs(trk.tpcNSigmaTr()) < cfgnTPCPIDTr) {
961960
jetHist.fill(HIST("tracks/triton/h2TOFmassTritonVsPt"), massTOF, trk.pt());
962961
jetHist.fill(HIST("tracks/triton/h2TOFmass2TritonVsPt"), massTOF * massTOF - gMassTriton * gMassTriton, trk.pt());
963962
jetHist.fill(HIST("tracks/triton/h2TofNsigmaTritonVsPt"), trk.tofNSigmaTr(), trk.pt());
964963
if (jetFlagPerpCone && isWithLeadingJet)
965964
jetHist.fill(HIST("tracks/perpCone/triton/h2TofNsigmaTritonVsPt"), trk.tofNSigmaTr(), trk.pt());
966965
}
967-
if (trk.tpcNSigmaHe() < useTPCpreSel) {
966+
if (std::abs(trk.tpcNSigmaHe()) < cfgnTPCPIDHe) {
968967
jetHist.fill(HIST("tracks/helium/h2TOFmassHeliumVsPt"), massTOF, trk.pt());
969968
jetHist.fill(HIST("tracks/helium/h2TOFmass2HeliumVsPt"), massTOF * massTOF - gMassHelium * gMassHelium, trk.pt());
970969
jetHist.fill(HIST("tracks/helium/h2TofNsigmaHeliumVsPt"), trk.tofNSigmaHe(), trk.pt());
@@ -978,7 +977,6 @@ struct nucleiInJets {
978977
jetHist.fill(HIST("tracks/antiDeuteron/h3PtVsantiDeuteronNSigmaTPC"), trk.pt(), trk.tpcNSigmaDe()); // De
979978
jetHist.fill(HIST("tracks/antiHelium/h3PtVsantiHeliumNSigmaTPC"), trk.pt(), trk.tpcNSigmaHe()); // He
980979
jetHist.fill(HIST("tracks/antiTriton/h3PtVsantiTritonNSigmaTPC"), trk.pt(), trk.tpcNSigmaTr()); // Tr
981-
982980
// perpCone
983981
if (jetFlagPerpCone && isWithLeadingJet) {
984982
// antiparticle info
@@ -1033,28 +1031,28 @@ struct nucleiInJets {
10331031
jetHist.fill(HIST("tracks/antiHelium/h2TOFmass2antiHeliumVsPt"), massTOF * massTOF - gMassHelium * gMassHelium, trk.pt());
10341032
jetHist.fill(HIST("tracks/antiTriton/h2TofNsigmaantiTritonVsPt"), trk.tofNSigmaTr(), trk.pt());
10351033
} else {
1036-
if (trk.tpcNSigmaPr() < useTPCpreSel) {
1034+
if (std::abs(trk.tpcNSigmaPr()) < cfgnTPCPIDPr) {
10371035
jetHist.fill(HIST("tracks/antiProton/h2TOFmassantiProtonVsPt"), massTOF, trk.pt());
10381036
jetHist.fill(HIST("tracks/antiProton/h2TOFmass2antiProtonVsPt"), massTOF * massTOF - gMassProton * gMassProton, trk.pt());
10391037
jetHist.fill(HIST("tracks/antiProton/h2TofNsigmaantiProtonVsPt"), trk.tofNSigmaPr(), trk.pt());
10401038
if (jetFlagPerpCone && isWithLeadingJet)
10411039
jetHist.fill(HIST("tracks/perpCone/antiProton/h2TofNsigmaantiProtonVsPt"), trk.tofNSigmaPr(), trk.pt());
10421040
}
1043-
if (trk.tpcNSigmaDe() < useTPCpreSel) {
1041+
if (std::abs(trk.tpcNSigmaDe()) < cfgnTPCPIDDe) {
10441042
jetHist.fill(HIST("tracks/antiDeuteron/h2TOFmassantiDeuteronVsPt"), massTOF, trk.pt());
10451043
jetHist.fill(HIST("tracks/antiDeuteron/h2TOFmass2antiDeuteronVsPt"), massTOF * massTOF - gMassDeuteron * gMassDeuteron, trk.pt());
10461044
jetHist.fill(HIST("tracks/antiDeuteron/h2TofNsigmaantiDeuteronVsPt"), trk.tofNSigmaDe(), trk.pt());
10471045
if (jetFlagPerpCone && isWithLeadingJet)
10481046
jetHist.fill(HIST("tracks/perpCone/antiDeuteron/h2TofNsigmaantiDeuteronVsPt"), trk.tofNSigmaDe(), trk.pt());
10491047
}
1050-
if (trk.tpcNSigmaTr() < useTPCpreSel) {
1048+
if (std::abs(trk.tpcNSigmaTr()) < cfgnTPCPIDTr) {
10511049
jetHist.fill(HIST("tracks/antiTriton/h2TOFmassantiTritonVsPt"), massTOF, trk.pt());
10521050
jetHist.fill(HIST("tracks/antiTriton/h2TOFmass2antiTritonVsPt"), massTOF * massTOF - gMassTriton * gMassTriton, trk.pt());
10531051
jetHist.fill(HIST("tracks/antiTriton/h2TofNsigmaantiTritonVsPt"), trk.tofNSigmaTr(), trk.pt());
10541052
if (jetFlagPerpCone && isWithLeadingJet)
10551053
jetHist.fill(HIST("tracks/perpCone/antiTriton/h2TofNsigmaantiTritonVsPt"), trk.tofNSigmaTr(), trk.pt());
10561054
}
1057-
if (trk.tpcNSigmaHe() < useTPCpreSel) {
1055+
if (std::abs(trk.tpcNSigmaHe()) < cfgnTPCPIDHe) {
10581056
jetHist.fill(HIST("tracks/antiHelium/h2TOFmassantiHeliumVsPt"), massTOF, trk.pt());
10591057
jetHist.fill(HIST("tracks/antiHelium/h2TOFmass2antiHeliumVsPt"), massTOF * massTOF - gMassHelium * gMassHelium, trk.pt());
10601058
jetHist.fill(HIST("tracks/antiHelium/h2TofNsigmaantiHeliumVsPt"), trk.tofNSigmaHe(), trk.pt());
@@ -1082,7 +1080,7 @@ struct nucleiInJets {
10821080
}
10831081
jetHist.fill(HIST("hNEvents"), 1.5);
10841082
}
1085-
if (fabs(collision.posZ()) > 10)
1083+
if (std::abs(collision.posZ()) > 10)
10861084
return;
10871085
jetHist.fill(HIST("hNEvents"), 2.5);
10881086
if (!jetderiveddatautilities::selectCollision(collision, jetderiveddatautilities::initialiseEventSelectionBits("sel8")))
@@ -1133,7 +1131,7 @@ struct nucleiInJets {
11331131
}
11341132
jetHist.fill(HIST("hNEvents"), 1.5);
11351133
}
1136-
if (fabs(collision.posZ()) > 10)
1134+
if (std::abs(collision.posZ()) > 10)
11371135
return;
11381136
jetHist.fill(HIST("hNEvents"), 2.5);
11391137
if (!jetderiveddatautilities::selectCollision(collision, jetderiveddatautilities::initialiseEventSelectionBits("sel8")))
@@ -1175,7 +1173,7 @@ struct nucleiInJets {
11751173
jetHist.fill(HIST("mcpJet/eventStat"), 0.5);
11761174
jetHist.fill(HIST("mcpJet/eventStat"), 1.5);
11771175

1178-
if (fabs(collision.posZ()) > 10) // bad vertex
1176+
if (std::abs(collision.posZ()) > 10) // bad vertex
11791177
return;
11801178

11811179
jetHist.fill(HIST("mcpJet/eventStat"), 2.5);
@@ -1184,7 +1182,7 @@ struct nucleiInJets {
11841182

11851183
bool INELgt0 = false;
11861184
for (const auto& mcParticle : mcParticles) {
1187-
if (fabs(mcParticle.eta()) < cfgtrkMaxEta) {
1185+
if (std::fabs(mcParticle.eta()) < cfgtrkMaxEta) {
11881186
INELgt0 = true;
11891187
break;
11901188
}
@@ -1205,9 +1203,9 @@ struct nucleiInJets {
12051203
for (const auto& mcParticle : mcParticles) {
12061204
if (!mcParticle.isPhysicalPrimary())
12071205
continue;
1208-
if (fabs(mcParticle.eta()) > cfgtrkMaxEta)
1206+
if (std::fabs(mcParticle.eta()) > cfgtrkMaxEta)
12091207
continue;
1210-
if (fabs(mcParticle.y()) > cfgtrkMaxRap)
1208+
if (std::fabs(mcParticle.y()) > cfgtrkMaxRap)
12111209
continue;
12121210

12131211
bool jetFlag = false;
@@ -1236,13 +1234,13 @@ struct nucleiInJets {
12361234
// bool jetFlag = kFALSE;
12371235
jetHist.fill(HIST("mcdJet/eventStat"), 1.5);
12381236

1239-
if (fabs(collisionJet.posZ()) > 10)
1237+
if (std::abs(collisionJet.posZ()) > 10)
12401238
return;
12411239

12421240
jetHist.fill(HIST("mcdJet/eventStat"), 2.5);
12431241
bool INELgt0 = false;
12441242
for (const auto& track : tracks) {
1245-
if (fabs(track.eta()) < cfgtrkMaxEta) {
1243+
if (std::fabs(track.eta()) < cfgtrkMaxEta) {
12461244
INELgt0 = true;
12471245
break;
12481246
}
@@ -1275,7 +1273,7 @@ struct nucleiInJets {
12751273
if (!track.has_mcParticle())
12761274
continue;
12771275
auto mcTrack = track.mcParticle_as<aod::JetParticles>();
1278-
if (fabs(mcTrack.eta()) > cfgtrkMaxEta)
1276+
if (std::fabs(mcTrack.eta()) > cfgtrkMaxEta)
12791277
continue;
12801278
if (!mcTrack.isPhysicalPrimary())
12811279
continue;
@@ -1319,15 +1317,15 @@ struct nucleiInJets {
13191317
soa::Join<aod::JetTracks, aod::JTrackPIs, aod::JMcTrackLbs> const& tracks,
13201318
JetMCPartTable const&, TrackCandidatesMC const&, aod::JetParticles const& particleTracks, aod::JMcCollisions const&)
13211319
{
1322-
if (fabs(collision.posZ()) > 10)
1320+
if (std::abs(collision.posZ()) > 10)
13231321
return;
13241322
if (!jetderiveddatautilities::selectCollision(collision, jetderiveddatautilities::initialiseEventSelectionBits("sel8")))
13251323
return;
13261324

13271325
jetHist.fill(HIST("recmatched/vertexZ"), collision.posZ());
13281326
bool INELgt0 = false;
13291327
for (const auto& track : tracks) {
1330-
if (fabs(track.eta()) < cfgtrkMaxEta) {
1328+
if (std::fabs(track.eta()) < cfgtrkMaxEta) {
13311329
INELgt0 = true;
13321330
break;
13331331
}
@@ -1403,7 +1401,7 @@ struct nucleiInJets {
14031401

14041402
for (const auto& track : tracks) {
14051403
auto completeTrack = track.track_as<TrackCandidatesMC>();
1406-
if (fabs(completeTrack.eta()) > cfgtrkMaxEta)
1404+
if (std::fabs(completeTrack.eta()) > cfgtrkMaxEta)
14071405
continue;
14081406
if (!isTrackSelected(completeTrack))
14091407
continue;
@@ -1412,7 +1410,7 @@ struct nucleiInJets {
14121410
auto mcTrack = track.mcParticle_as<aod::JetParticles>();
14131411
if (!mcTrack.isPhysicalPrimary())
14141412
continue;
1415-
if (fabs(mcTrack.y()) > cfgtrkMaxRap)
1413+
if (std::fabs(mcTrack.y()) > cfgtrkMaxRap)
14161414
continue;
14171415

14181416
bool isTpcPassed(true);
@@ -1487,9 +1485,9 @@ struct nucleiInJets {
14871485
for (const auto& mcParticle : mcParticles_per_coll) {
14881486
if (!mcParticle.isPhysicalPrimary())
14891487
continue;
1490-
if (fabs(mcParticle.eta()) > cfgtrkMaxEta)
1488+
if (std::fabs(mcParticle.eta()) > cfgtrkMaxEta)
14911489
continue;
1492-
if (fabs(mcParticle.y()) > cfgtrkMaxRap)
1490+
if (std::fabs(mcParticle.y()) > cfgtrkMaxRap)
14931491
continue;
14941492
bool jetFlagMC = false;
14951493
bool jetFlagPerpConeMC = false;
@@ -1546,7 +1544,7 @@ struct nucleiInJets {
15461544
if (!jetderiveddatautilities::selectCollision(recocoll, jetderiveddatautilities::initialiseEventSelectionBits("sel8")))
15471545
return;
15481546
}
1549-
if (fabs(collision.posZ()) > 10)
1547+
if (std::abs(collision.posZ()) > 10)
15501548
return;
15511549
jetHist.fill(HIST("genmatched/vertexZ"), collision.posZ());
15521550
std::vector<double> mcdJetPt{};
@@ -1612,7 +1610,7 @@ struct nucleiInJets {
16121610
} // leading jet only
16131611

16141612
for (const auto& mcParticle : mcParticles) {
1615-
if (fabs(mcParticle.eta()) > cfgtrkMaxEta)
1613+
if (std::fabs(mcParticle.eta()) > cfgtrkMaxEta)
16161614
continue;
16171615
// add pid later
16181616

0 commit comments

Comments
 (0)