Skip to content

Commit f55fe11

Browse files
alibuildalibuild
committed
Please consider the following formatting changes
1 parent ed53e44 commit f55fe11

File tree

1 file changed

+75
-65
lines changed

1 file changed

+75
-65
lines changed

‎PWGCF/EbyEFluctuations/Tasks/kaonIsospinFluctuations.cxx‎

Lines changed: 75 additions & 65 deletions
Original file line numberDiff line numberDiff line change
@@ -85,17 +85,17 @@ struct KaonIsospinFluctuations {
8585
Configurable<float> cfgK0sLowPt{"cfgK0sLowPt", 0.1, "cfgK0sLowPt"};
8686
Configurable<double> cfgK0sHighPt{"cfgK0sHighPt", 1.5, "cfgK0sHighPt"};
8787
Configurable<double> cfgK0sRapitidy{"cfgK0sRapitidy", 0.5, "cfgK0sRapitidy"};
88-
Configurable<double> cfgK0sARMcut{"cfgK0sARMcut", 0.2, "cfgK0sARMcut"};
88+
Configurable<double> cfgK0sARMcut{"cfgK0sARMcut", 0.2, "cfgK0sARMcut"};
8989
} k0sSelCut;
9090

91-
// Histogram Configurables
91+
// Histogram Configurables
9292
struct : ConfigurableGroup {
9393
Configurable<int> centBins{"centBins", 1020, "No of bins in centrality axis"};
9494
Configurable<double> centBinsxLow{"centBinsxLow", -1.0, "centBinsxLow"};
9595
Configurable<double> centBinsxUp{"centBinsxUp", 101.0, "centBinsxUp"};
9696
Configurable<int> centAxisType{"centAxisType", 0, "centAxisType"};
9797
} cfgCentAxis;
98-
//Track Configurables
98+
// Track Configurables
9999
struct : ConfigurableGroup {
100100
Configurable<int> cfgTrkTpcNClsCrossedRows{"cfgTrkTpcNClsCrossedRows", 70, "cfgTrkTpcNClsCrossedRows"};
101101
Configurable<float> cfgTrkdcaXY{"cfgTrkdcaXY", 0.2, "cfgTrkdcaXY"};
@@ -177,7 +177,7 @@ struct KaonIsospinFluctuations {
177177
Configurable<bool> doFWDPtDependentCheck{"doFWDPtDependentCheck", 1, "doFWDPtDependentCheck"};
178178
Configurable<float> cfgFWDPtCut{"cfgFWDPtCut", 0.2, "cfgFWDPtCut"};
179179
Configurable<std::vector<int>> cfgFinalParticleIdList{"cfgFinalParticleIdList", {11, -11, 13, -13, 15, -15, 211, -211, 321, -321, 2212, -2212}, "cfgFinalParticleIdList"};
180-
Configurable<std::vector<int>> cfgNonFinalParticleIdList{"cfgNonFinalParticleIdList", {11, -11, 13, -13, 15, -15, 211, -211, 321, -321, 2212, -2212}, "cfgNonFinalParticleIdList"};
180+
Configurable<std::vector<int>> cfgNonFinalParticleIdList{"cfgNonFinalParticleIdList", {11, -11, 13, -13, 15, -15, 211, -211, 321, -321, 2212, -2212}, "cfgNonFinalParticleIdList"};
181181

182182
void init(InitContext const&)
183183
{
@@ -239,7 +239,7 @@ struct KaonIsospinFluctuations {
239239
const AxisSpec axisParticleCount3 = {1060, -10, 1050, "particleCount"};
240240

241241
const AxisSpec axisArmenterosAlpha = {100, -1.0, 1.0, "ArmenterosAlpha"};
242-
const AxisSpec axisArmenterosQt = {150, 0, 0.3, "ArmenterosQt"};
242+
const AxisSpec axisArmenterosQt = {150, 0, 0.3, "ArmenterosQt"};
243243

244244
HistogramConfigSpec histPDcaXY({HistType::kTH2F, {axisP, axisDcaXY}});
245245
HistogramConfigSpec histPtDcaXY({HistType::kTH2F, {axisPt, axisDcaXY}});
@@ -295,7 +295,7 @@ struct KaonIsospinFluctuations {
295295
recoV0s.add("v0Table/Full/h16_rapidity", "rapidity", kTH1F, {axisRapidity});
296296
recoV0s.add("v0Table/Full/h17_alpha", "alpha", kTH1F, {axisArmenterosAlpha});
297297
recoV0s.add("v0Table/Full/h18_qtarm", "qtarm", kTH1F, {axisArmenterosQt});
298-
recoV0s.add("v0Table/Full/h19_alpha_qtarm","alpha_qtarm",kTH2F, {axisArmenterosAlpha, axisArmenterosQt});
298+
recoV0s.add("v0Table/Full/h19_alpha_qtarm", "alpha_qtarm", kTH2F, {axisArmenterosAlpha, axisArmenterosQt});
299299
recoV0s.add("v0Table/Full/h20_pt_eta", "pt_eta", kTH2F, {axisPt, axisEta});
300300

301301
// K0s-Daughter Info
@@ -345,9 +345,15 @@ struct KaonIsospinFluctuations {
345345
recoV0s.addClone("v0Table/Full/Pi/tpcId/", "v0Table/Full/Pi/tpctofId/"); // for identification using tof+tpc
346346
recoV0s.addClone("v0Table/Full/Pi/tpcId/", "v0Table/Full/Pi/NoId/"); // for unidentified case // to observe and debug
347347

348-
if (cfgFillV0TablePostK0sCheck ){recoV0s.addClone("v0Table/Full/", "v0Table/postK0sCheck/"); }
349-
if (cfgFillV0TablePostMassCut ){recoV0s.addClone("v0Table/Full/", "v0Table/postMassCut/"); }
350-
if (cfgFillV0TablePostSelectionCut){recoV0s.addClone("v0Table/Full/", "v0Table/postSelectionCut/");}
348+
if (cfgFillV0TablePostK0sCheck) {
349+
recoV0s.addClone("v0Table/Full/", "v0Table/postK0sCheck/");
350+
}
351+
if (cfgFillV0TablePostMassCut) {
352+
recoV0s.addClone("v0Table/Full/", "v0Table/postMassCut/");
353+
}
354+
if (cfgFillV0TablePostSelectionCut) {
355+
recoV0s.addClone("v0Table/Full/", "v0Table/postSelectionCut/");
356+
}
351357

352358
recoV0s.add("v0Table/postSelectionCut/hTrueV0TagCount", "hTrueV0TagCount", {HistType::kTH1F, {{12, -2, 10}}}); // 001 = Kaon, 010 = Lambda, 100 = AnitLambda
353359
recoV0s.add("v0Table/postSelectionCut/nCommonPionOfDifferentK0s", "nCommonPionOfDifferentK0s", {HistType::kTH1D, {{44, -2, 20}}});
@@ -398,7 +404,7 @@ struct KaonIsospinFluctuations {
398404
recoK0s.add("recoK0s/PreSel/h16_rapidity", "rapidity", kTH1F, {axisRapidity});
399405
recoK0s.add("recoK0s/PreSel/h17_alpha", "alpha", kTH1F, {axisArmenterosAlpha});
400406
recoK0s.add("recoK0s/PreSel/h18_qtarm", "qtarm", kTH1F, {axisArmenterosQt});
401-
recoK0s.add("recoK0s/PreSel/h19_alpha_qtarm","alpha_qtarm",kTH2F, {axisArmenterosAlpha, axisArmenterosQt});
407+
recoK0s.add("recoK0s/PreSel/h19_alpha_qtarm", "alpha_qtarm", kTH2F, {axisArmenterosAlpha, axisArmenterosQt});
402408
recoK0s.add("recoK0s/PreSel/h20_pt_eta", "pt_eta", kTH2F, {axisPt, axisEta});
403409

404410
// K0s-Daughter Info
@@ -590,7 +596,7 @@ struct KaonIsospinFluctuations {
590596
genAnalysis.print();
591597
}
592598

593-
enum RejectionTagEnum{
599+
enum RejectionTagEnum {
594600
kPassed = 0,
595601
kFailTpcNClsCrossedRows,
596602
kFailTrkdcaXY,
@@ -1291,10 +1297,9 @@ struct KaonIsospinFluctuations {
12911297
bool selK0s(T v0)
12921298
{
12931299
if (k0sSelCut.cfgK0sMLow < v0.mK0Short() && v0.mK0Short() < k0sSelCut.cfgK0sMHigh &&
1294-
k0sSelCut.cfgK0sLowPt < v0.pt() && v0.pt() < k0sSelCut.cfgK0sHighPt &&
1295-
std::abs(v0.rapidity(MassK0Short)) < k0sSelCut.cfgK0sRapitidy &&
1296-
v0.qtarm() > (k0sSelCut.cfgK0sARMcut * std::abs(v0.alpha()))
1297-
) {
1300+
k0sSelCut.cfgK0sLowPt < v0.pt() && v0.pt() < k0sSelCut.cfgK0sHighPt &&
1301+
std::abs(v0.rapidity(MassK0Short)) < k0sSelCut.cfgK0sRapitidy &&
1302+
v0.qtarm() > (k0sSelCut.cfgK0sARMcut * std::abs(v0.alpha()))) {
12981303
return true;
12991304
} else {
13001305
return false;
@@ -1322,27 +1327,32 @@ struct KaonIsospinFluctuations {
13221327
bool checkTrackSelection(const T& track, int& rejectionTag)
13231328
{
13241329
if (track.tpcNClsCrossedRows() < cfgTrackCuts.cfgTrkTpcNClsCrossedRows) {
1325-
rejectionTag = kFailTpcNClsCrossedRows; return false;
1330+
rejectionTag = kFailTpcNClsCrossedRows;
1331+
return false;
13261332
}
13271333
if (std::fabs(track.dcaXY()) > cfgTrackCuts.cfgTrkdcaXY) {
1328-
rejectionTag = kFailTrkdcaXY; return false;
1334+
rejectionTag = kFailTrkdcaXY;
1335+
return false;
13291336
}
13301337
if (!track.isGlobalTrack()) {
1331-
rejectionTag = kFailGlobalTrack; return false;
1338+
rejectionTag = kFailGlobalTrack;
1339+
return false;
13321340
}
1333-
if (cfgTrackCuts.cfgDoVGselTrackCheck){
1334-
if(!selTrackForId(track)) {
1335-
rejectionTag = kFailVGSelCheck; return false;
1341+
if (cfgTrackCuts.cfgDoVGselTrackCheck) {
1342+
if (!selTrackForId(track)) {
1343+
rejectionTag = kFailVGSelCheck;
1344+
return false;
13361345
}
13371346
}
13381347
return true;
13391348
}
13401349

13411350
template <typename T>
13421351
bool checkTrackInList(const T& track, const std::vector<int64_t>& vecList, int& rejectionTag, const int& listTagValue)
1343-
{
1344-
if(std::binary_search(vecList.begin(), vecList.end(), track.globalIndex())){
1345-
rejectionTag = listTagValue; return true; //Binary Search is fastest search in a sorted array.
1352+
{
1353+
if (std::binary_search(vecList.begin(), vecList.end(), track.globalIndex())) {
1354+
rejectionTag = listTagValue;
1355+
return true; // Binary Search is fastest search in a sorted array.
13461356
}
13471357
return false;
13481358
}
@@ -1439,7 +1449,7 @@ struct KaonIsospinFluctuations {
14391449
recoTracks.fill(HIST(HistRegDire[Mode]) + HIST("h16_tpcInnerParam_dcaZ"), track.tpcInnerParam(), track.dcaZ());
14401450
recoTracks.fill(HIST(HistRegDire[Mode]) + HIST("h17_tofExpMom_dcaZ"), track.tofExpMom(), track.dcaZ());
14411451

1442-
recoTracks.fill(HIST(HistRegDire[Mode]) + HIST("h20_pt_eta"), track.pt(), track.eta());
1452+
recoTracks.fill(HIST(HistRegDire[Mode]) + HIST("h20_pt_eta"), track.pt(), track.eta());
14431453
// momemtum
14441454
recoTracks.fill(HIST(HistRegDire[Mode]) + HIST("h20_p_pt"), track.p(), track.pt());
14451455
recoTracks.fill(HIST(HistRegDire[Mode]) + HIST("h21_p_tpcInnerParam"), track.p(), track.tpcInnerParam());
@@ -1514,8 +1524,8 @@ struct KaonIsospinFluctuations {
15141524
histReg.fill(HIST(HistRegDire[Mode]) + HIST("h16_rapidity"), v0.rapidity(MassK0Short));
15151525
histReg.fill(HIST(HistRegDire[Mode]) + HIST("h17_alpha"), v0.alpha());
15161526
histReg.fill(HIST(HistRegDire[Mode]) + HIST("h18_qtarm"), v0.qtarm());
1517-
histReg.fill(HIST(HistRegDire[Mode]) + HIST("h19_alpha_qtarm"),v0.alpha(),v0.qtarm());
1518-
histReg.fill(HIST(HistRegDire[Mode]) + HIST("h20_pt_eta"),v0.pt(), v0.eta());
1527+
histReg.fill(HIST(HistRegDire[Mode]) + HIST("h19_alpha_qtarm"), v0.alpha(), v0.qtarm());
1528+
histReg.fill(HIST(HistRegDire[Mode]) + HIST("h20_pt_eta"), v0.pt(), v0.eta());
15191529

15201530
if (posPiIdMethod == kTPCidentified) {
15211531
fillV0DaughterQA<Mode, kPi, tpcId>(histReg, posDaughterTrack, MassProton);
@@ -1542,7 +1552,7 @@ struct KaonIsospinFluctuations {
15421552
histReg.fill(HIST(HistRegDire[Mode]) + HIST(PidDire[pidMode]) + HIST("h14_eta"), mcTrack.eta());
15431553
histReg.fill(HIST(HistRegDire[Mode]) + HIST(PidDire[pidMode]) + HIST("h15_phi"), mcTrack.phi());
15441554
histReg.fill(HIST(HistRegDire[Mode]) + HIST(PidDire[pidMode]) + HIST("h16_rapidity"), mcTrack.y());
1545-
histReg.fill(HIST(HistRegDire[Mode]) + HIST(PidDire[pidMode]) + HIST("h20_pt_eta"), mcTrack.pt(), mcTrack.eta());
1555+
histReg.fill(HIST(HistRegDire[Mode]) + HIST(PidDire[pidMode]) + HIST("h20_pt_eta"), mcTrack.pt(), mcTrack.eta());
15461556
}
15471557

15481558
template <typename T, typename U, typename H>
@@ -1595,7 +1605,7 @@ struct KaonIsospinFluctuations {
15951605
} // End of K0s block
15961606
}
15971607

1598-
template<typename T>
1608+
template <typename T>
15991609
void executeSortPairDaughters(const auto& posDauList, const auto& negDauList, auto& fullDauList, const T& hist)
16001610
{
16011611
findRepeatedEntries(posDauList, hist);
@@ -1660,12 +1670,12 @@ struct KaonIsospinFluctuations {
16601670
recoAnalysis.fill(HIST("recoAnalysis/RejectedTrack_RejectionTag"), rejectionTag);
16611671
isAcceptedTrack = false;
16621672
return; // for continue;
1663-
} else if(checkTrackInList(track, fullDauList, rejectionTag, kFailK0ShortDaughter)){
1673+
} else if (checkTrackInList(track, fullDauList, rejectionTag, kFailK0ShortDaughter)) {
16641674
recoAnalysis.fill(HIST("recoAnalysis/RejectedTrack_RejectionTag"), rejectionTag);
16651675
if (track.signed1Pt() > 0) {
1666-
nRejectedPiPlus++; //DOEFFCORR
1676+
nRejectedPiPlus++; // DOEFFCORR
16671677
} else if (track.signed1Pt() < 0) {
1668-
nRejectedPiMinus++; //DOEFFCORR
1678+
nRejectedPiMinus++; // DOEFFCORR
16691679
}
16701680
isAcceptedTrack = false;
16711681
return; // for continue;
@@ -1677,7 +1687,6 @@ struct KaonIsospinFluctuations {
16771687
nTrack++;
16781688
}
16791689

1680-
16811690
template <typename T>
16821691
void executeTrackAnalysisPart(const T& track, const int& trackIdTag,
16831692
const int& idMethodPi, const bool& trackIsPion, float& nPiMinus, float& nPiPlus,
@@ -1687,11 +1696,11 @@ struct KaonIsospinFluctuations {
16871696
const int& idMethodDe, const bool& trackIsDeuteron, float& nDePlus, float& nDeMinus)
16881697
{
16891698
if (trackIsPion) {
1690-
if (idMethodPi == kTPCidentified){
1699+
if (idMethodPi == kTPCidentified) {
16911700
fillIdentificationQA<recoAnalysisDir, kPi, tpcId, true>(recoAnalysis, track);
1692-
} else if (idMethodPi == kTPCTOFidentified){
1701+
} else if (idMethodPi == kTPCTOFidentified) {
16931702
fillIdentificationQA<recoAnalysisDir, kPi, tpctofId, true>(recoAnalysis, track);
1694-
} else if (idMethodPi == kUnidentified){
1703+
} else if (idMethodPi == kUnidentified) {
16951704
fillIdentificationQA<recoAnalysisDir, kPi, NoId, true>(recoAnalysis, track);
16961705
}
16971706
if (track.sign() > 0) {
@@ -1701,7 +1710,7 @@ struct KaonIsospinFluctuations {
17011710
}
17021711
}
17031712
if (trackIsKaon) {
1704-
if (idMethodKa == kTPCidentified){
1713+
if (idMethodKa == kTPCidentified) {
17051714
fillIdentificationQA<recoAnalysisDir, kKa, tpcId, true>(recoAnalysis, track);
17061715
} else if (idMethodKa == kTPCTOFidentified) {
17071716
fillIdentificationQA<recoAnalysisDir, kKa, tpctofId, true>(recoAnalysis, track);
@@ -1710,16 +1719,16 @@ struct KaonIsospinFluctuations {
17101719
}
17111720
if (track.sign() > 0) {
17121721
nKaPlus++;
1713-
} else if (track.sign() < 0) {
1722+
} else if (track.sign() < 0) {
17141723
nKaMinus++;
17151724
}
17161725
}
17171726
if (trackIsProton) {
1718-
if (idMethodPr == kTPCidentified){
1727+
if (idMethodPr == kTPCidentified) {
17191728
fillIdentificationQA<recoAnalysisDir, kPr, tpcId, true>(recoAnalysis, track);
17201729
} else if (idMethodPr == kTPCTOFidentified) {
17211730
fillIdentificationQA<recoAnalysisDir, kPr, tpctofId, true>(recoAnalysis, track);
1722-
} else if (idMethodPr == kUnidentified){
1731+
} else if (idMethodPr == kUnidentified) {
17231732
fillIdentificationQA<recoAnalysisDir, kPr, NoId, true>(recoAnalysis, track);
17241733
}
17251734
if (track.sign() > 0) {
@@ -1729,11 +1738,11 @@ struct KaonIsospinFluctuations {
17291738
}
17301739
}
17311740
if (trackIsElectron) {
1732-
if (idMethodEl == kTPCidentified){
1741+
if (idMethodEl == kTPCidentified) {
17331742
fillIdentificationQA<recoAnalysisDir, kEl, tpcId, true>(recoAnalysis, track);
17341743
} else if (idMethodEl == kTPCTOFidentified) {
17351744
fillIdentificationQA<recoAnalysisDir, kEl, tpctofId, true>(recoAnalysis, track);
1736-
} else if (idMethodEl == kUnidentified){
1745+
} else if (idMethodEl == kUnidentified) {
17371746
fillIdentificationQA<recoAnalysisDir, kEl, NoId, true>(recoAnalysis, track);
17381747
}
17391748
if (track.sign() > 0) {
@@ -1743,9 +1752,9 @@ struct KaonIsospinFluctuations {
17431752
}
17441753
}
17451754
if (trackIsDeuteron) {
1746-
if (idMethodDe == kTPCidentified){
1755+
if (idMethodDe == kTPCidentified) {
17471756
fillIdentificationQA<recoAnalysisDir, kDe, tpcId, true>(recoAnalysis, track);
1748-
} else if (idMethodDe == kTPCTOFidentified){
1757+
} else if (idMethodDe == kTPCTOFidentified) {
17491758
fillIdentificationQA<recoAnalysisDir, kDe, tpctofId, true>(recoAnalysis, track);
17501759
} else if (idMethodDe == kUnidentified) {
17511760
fillIdentificationQA<recoAnalysisDir, kDe, NoId, true>(recoAnalysis, track);
@@ -2398,13 +2407,14 @@ struct KaonIsospinFluctuations {
23982407
}
23992408
PROCESS_SWITCH(KaonIsospinFluctuations, processGen, "Process for Gen", false);
24002409

2401-
template<typename T>
2402-
void getV0MCount(const T& mcTrack, float& multV0M){
2403-
if((-3.7 < mcTrack.eta() && mcTrack.eta() < -1.7 ) || (2.8 < mcTrack.eta() && mcTrack.eta() < 5.1)){
2404-
if(doFWDPtDependentCheck){
2405-
if(mcTrack.pt() > cfgFWDPtCut){
2410+
template <typename T>
2411+
void getV0MCount(const T& mcTrack, float& multV0M)
2412+
{
2413+
if ((-3.7 < mcTrack.eta() && mcTrack.eta() < -1.7) || (2.8 < mcTrack.eta() && mcTrack.eta() < 5.1)) {
2414+
if (doFWDPtDependentCheck) {
2415+
if (mcTrack.pt() > cfgFWDPtCut) {
24062416
multV0M++; // V0C: at -3.7 < η < -1.7 (backward direction).
2407-
// V0A: at 2.8 < η < 5.1 (forward direction).
2417+
// V0A: at 2.8 < η < 5.1 (forward direction).
24082418
}
24092419
} else {
24102420
multV0M++;
@@ -2424,8 +2434,8 @@ struct KaonIsospinFluctuations {
24242434
for (const auto& mcColl : mcCollisions) {
24252435
centrality = -1;
24262436

2427-
if(cfgVtxZCheck){
2428-
if(std::abs(mcColl.posZ()) >= cutZvertex){
2437+
if (cfgVtxZCheck) {
2438+
if (std::abs(mcColl.posZ()) >= cutZvertex) {
24292439
continue;
24302440
}
24312441
}
@@ -2452,32 +2462,33 @@ struct KaonIsospinFluctuations {
24522462

24532463
for (const auto& mcTrack : mcTracksTablePerMcColl) {
24542464

2455-
if(cfgCountFinalParticles){
2456-
if(!mcTrack.has_daughters() && std::binary_search(finalParticleIdList.begin(), finalParticleIdList.end(), mcTrack.pdgCode())){
2465+
if (cfgCountFinalParticles) {
2466+
if (!mcTrack.has_daughters() && std::binary_search(finalParticleIdList.begin(), finalParticleIdList.end(), mcTrack.pdgCode())) {
24572467
getV0MCount(mcTrack, multV0M);
24582468
}
2459-
}
2460-
if(cfgCountNonFinalParticles){
2461-
if(mcTrack.has_daughters() && std::binary_search(nonFinalParticleIdList.begin(), nonFinalParticleIdList.end(), mcTrack.pdgCode())){
2469+
}
2470+
if (cfgCountNonFinalParticles) {
2471+
if (mcTrack.has_daughters() && std::binary_search(nonFinalParticleIdList.begin(), nonFinalParticleIdList.end(), mcTrack.pdgCode())) {
24622472
getV0MCount(mcTrack, multV0M);
24632473
}
24642474
}
24652475

2466-
if(cfgCountPhysicalPrimAndFinalParticles){
2467-
if(!mcTrack.has_daughters() && mcTrack.isPhysicalPrimary()){
2468-
if(!(std::abs(mcTrack.pdgCode()) == 12 || std::abs(mcTrack.pdgCode()) == 14 ||std::abs(mcTrack.pdgCode()) == 16)){
2469-
//Removed invisible neutrinos;
2476+
if (cfgCountPhysicalPrimAndFinalParticles) {
2477+
if (!mcTrack.has_daughters() && mcTrack.isPhysicalPrimary()) {
2478+
if (!(std::abs(mcTrack.pdgCode()) == 12 || std::abs(mcTrack.pdgCode()) == 14 || std::abs(mcTrack.pdgCode()) == 16)) {
2479+
// Removed invisible neutrinos;
24702480
getV0MCount(mcTrack, multV0M);
24712481
}
2472-
}}
2482+
}
2483+
}
24732484

24742485
if (!mcTrack.isPhysicalPrimary()) {
24752486
continue;
24762487
}
24772488

24782489
if (mcTrack.pdgCode() == kK0Short &&
2479-
k0sSelCut.cfgK0sLowPt < mcTrack.pt() && mcTrack.pt() < k0sSelCut.cfgK0sHighPt &&
2480-
std::abs(mcTrack.y()) < k0sSelCut.cfgK0sRapitidy) {
2490+
k0sSelCut.cfgK0sLowPt < mcTrack.pt() && mcTrack.pt() < k0sSelCut.cfgK0sHighPt &&
2491+
std::abs(mcTrack.y()) < k0sSelCut.cfgK0sRapitidy) {
24812492
nK0s++;
24822493
fillGenTrackQA<genAnalysisDir, kK0s>(genAnalysis, mcTrack);
24832494
}
@@ -2534,7 +2545,6 @@ struct KaonIsospinFluctuations {
25342545
} // collision loop is over
25352546
}
25362547
PROCESS_SWITCH(KaonIsospinFluctuations, processSim, "Process for Sim", false);
2537-
25382548
};
25392549

25402550
WorkflowSpec defineDataProcessing(ConfigContext const& cfgc)

0 commit comments

Comments
 (0)