Skip to content

Commit 3247dfd

Browse files
ChiaraDeMartin95ChiaraDeMartin95
authored
Merge pull request #28 from alibuild/alibot-cleanup-12205
Please consider the following formatting changes to #12205
2 parents 14ccf35 + 955541b commit 3247dfd

File tree

1 file changed

+44
-41
lines changed

1 file changed

+44
-41
lines changed

‎PWGLF/TableProducer/Strangeness/cascadeflow.cxx‎

Lines changed: 44 additions & 41 deletions
Original file line numberDiff line numberDiff line change
@@ -1244,7 +1244,7 @@ struct cascadeFlow {
12441244
int counter = 0;
12451245
IsCascAccepted(casc, negExtra, posExtra, bachExtra, counter);
12461246
histos.fill(HIST("hCascade"), counter);
1247-
1247+
12481248
// ML selections
12491249
bool isSelectedCasc[nParticles]{false, false};
12501250

@@ -1506,11 +1506,13 @@ struct cascadeFlow {
15061506
int counterALambda = 0;
15071507
bool isLambdaCandidate = 0;
15081508
bool isALambdaCandidate = 0;
1509-
if (isLambdaAccepted(negExtra, posExtra, counterLambda)) isLambdaCandidate = 1;
1510-
if (isAntiLambdaAccepted(negExtra, posExtra, counterALambda)) isALambdaCandidate = 1;
1509+
if (isLambdaAccepted(negExtra, posExtra, counterLambda))
1510+
isLambdaCandidate = 1;
1511+
if (isAntiLambdaAccepted(negExtra, posExtra, counterALambda))
1512+
isALambdaCandidate = 1;
15111513
histos.fill(HIST("hLambdaDauSel"), counterLambda);
15121514
histos.fill(HIST("hALambdaDauSel"), counterALambda);
1513-
1515+
15141516
// pt cut
15151517
if (v0.pt() < V0Configs.MinPtV0 || v0.pt() > V0Configs.MaxPtV0) {
15161518
continue;
@@ -1527,26 +1529,28 @@ struct cascadeFlow {
15271529
isSelectedV0[1] = true;
15281530

15291531
int chargeIndex = -1;
1530-
if (isSelectedV0[0] && !isSelectedV0[1]) { //Lambdas
1531-
histos.fill(HIST("hLambdaCandidate"), 0);
1532-
chargeIndex = 0;
1532+
if (isSelectedV0[0] && !isSelectedV0[1]) { // Lambdas
1533+
histos.fill(HIST("hLambdaCandidate"), 0);
1534+
chargeIndex = 0;
15331535
}
1534-
if (isSelectedV0[1] && !isSelectedV0[0]) { //AntiLambdas
1535-
histos.fill(HIST("hLambdaCandidate"), 1);
1536-
chargeIndex = 1;
1536+
if (isSelectedV0[1] && !isSelectedV0[0]) { // AntiLambdas
1537+
histos.fill(HIST("hLambdaCandidate"), 1);
1538+
chargeIndex = 1;
15371539
}
15381540
if (isSelectedV0[0] && isSelectedV0[1]) {
1539-
histos.fill(HIST("hLambdaCandidate"), 2);
1540-
if (v0.mLambda() > V0Configs.MinMassLambda && v0.mLambda() < V0Configs.MaxMassLambda && v0.mAntiLambda() > V0Configs.MinMassLambda && v0.mAntiLambda() < V0Configs.MaxMassLambda) {
1541-
histos.fill(HIST("hLambdaCandidate"), 3);
1542-
continue; //in case of ambiguity between Lambda and AntiLambda, I skip the particle
1543-
}
1544-
if (v0.mLambda() > V0Configs.MinMassLambda && v0.mLambda() < V0Configs.MaxMassLambda) chargeIndex = 0;
1545-
else if (v0.mAntiLambda() > V0Configs.MinMassLambda && v0.mAntiLambda() < V0Configs.MaxMassLambda) chargeIndex = 1;
1546-
else {
1547-
histos.fill(HIST("hLambdaCandidate"), 4);
1548-
continue; //in case of ambiguity between Lambda and AntiLambda, I skip the particle
1549-
}
1541+
histos.fill(HIST("hLambdaCandidate"), 2);
1542+
if (v0.mLambda() > V0Configs.MinMassLambda && v0.mLambda() < V0Configs.MaxMassLambda && v0.mAntiLambda() > V0Configs.MinMassLambda && v0.mAntiLambda() < V0Configs.MaxMassLambda) {
1543+
histos.fill(HIST("hLambdaCandidate"), 3);
1544+
continue; // in case of ambiguity between Lambda and AntiLambda, I skip the particle
1545+
}
1546+
if (v0.mLambda() > V0Configs.MinMassLambda && v0.mLambda() < V0Configs.MaxMassLambda)
1547+
chargeIndex = 0;
1548+
else if (v0.mAntiLambda() > V0Configs.MinMassLambda && v0.mAntiLambda() < V0Configs.MaxMassLambda)
1549+
chargeIndex = 1;
1550+
else {
1551+
histos.fill(HIST("hLambdaCandidate"), 4);
1552+
continue; // in case of ambiguity between Lambda and AntiLambda, I skip the particle
1553+
}
15501554
}
15511555
if (!isSelectedV0[0] && !isSelectedV0[1])
15521556
continue;
@@ -1583,14 +1587,13 @@ struct cascadeFlow {
15831587
double cos2ThetaLambda = 0;
15841588
double cosThetaLambda = 0;
15851589
if (chargeIndex == 0) {
1586-
pzs2Lambda = cosThetaStarProton[0] * std::sin(2 * (v0.phi() - psiT0C)) / lambdav2::AlphaLambda[0] / meanCos2ThetaProtonFromLambda;
1587-
cos2ThetaLambda = cosThetaStarProton[0] * cosThetaStarProton[0];
1588-
cosThetaLambda = cosThetaStarProton[0] / cascadev2::AlphaLambda[0] / meanCos2ThetaProtonFromLambda;
1589-
}
1590-
else {
1591-
pzs2Lambda = cosThetaStarProton[1] * std::sin(2 * (v0.phi() - psiT0C)) / lambdav2::AlphaLambda[1] / meanCos2ThetaProtonFromLambda;
1592-
cos2ThetaLambda = cosThetaStarProton[1] * cosThetaStarProton[1];
1593-
cosThetaLambda = cosThetaStarProton[1] / cascadev2::AlphaLambda[1] / meanCos2ThetaProtonFromLambda;
1590+
pzs2Lambda = cosThetaStarProton[0] * std::sin(2 * (v0.phi() - psiT0C)) / lambdav2::AlphaLambda[0] / meanCos2ThetaProtonFromLambda;
1591+
cos2ThetaLambda = cosThetaStarProton[0] * cosThetaStarProton[0];
1592+
cosThetaLambda = cosThetaStarProton[0] / cascadev2::AlphaLambda[0] / meanCos2ThetaProtonFromLambda;
1593+
} else {
1594+
pzs2Lambda = cosThetaStarProton[1] * std::sin(2 * (v0.phi() - psiT0C)) / lambdav2::AlphaLambda[1] / meanCos2ThetaProtonFromLambda;
1595+
cos2ThetaLambda = cosThetaStarProton[1] * cosThetaStarProton[1];
1596+
cosThetaLambda = cosThetaStarProton[1] / cascadev2::AlphaLambda[1] / meanCos2ThetaProtonFromLambda;
15941597
}
15951598

15961599
histos.fill(HIST("hv2CEPvsFT0C"), coll.centFT0C(), v2CEP);
@@ -1599,19 +1602,19 @@ struct cascadeFlow {
15991602
histos.fill(HIST("hlambdaminuspsiT0C"), lambdaminuspsiT0C);
16001603

16011604
if (fillingConfigs.isFillTHNLambda) {
1602-
if (fillingConfigs.isFillTHN_V2)
1603-
histos.get<THn>(HIST("hLambdaV2"))->Fill(coll.centFT0C(), chargeIndex, v0.pt(), v0.mLambda(), v2CEP);
1604-
if (fillingConfigs.isFillTHN_Pz){
1605-
histos.get<THn>(HIST("hLambdaPzs2"))->Fill(coll.centFT0C(), chargeIndex, v0.pt(), v0.mLambda(), pzs2Lambda);
1606-
}
1607-
if (fillingConfigs.isFillTHN_Acc)
1608-
histos.get<THn>(HIST("hLambdaCos2Theta"))->Fill(coll.centFT0C(), chargeIndex, v0.eta(), v0.pt(), v0.mLambda(), cos2ThetaLambda);
1605+
if (fillingConfigs.isFillTHN_V2)
1606+
histos.get<THn>(HIST("hLambdaV2"))->Fill(coll.centFT0C(), chargeIndex, v0.pt(), v0.mLambda(), v2CEP);
1607+
if (fillingConfigs.isFillTHN_Pz) {
1608+
histos.get<THn>(HIST("hLambdaPzs2"))->Fill(coll.centFT0C(), chargeIndex, v0.pt(), v0.mLambda(), pzs2Lambda);
1609+
}
1610+
if (fillingConfigs.isFillTHN_Acc)
1611+
histos.get<THn>(HIST("hLambdaCos2Theta"))->Fill(coll.centFT0C(), chargeIndex, v0.eta(), v0.pt(), v0.mLambda(), cos2ThetaLambda);
16091612
}
16101613
if (fillingConfigs.isFillTHNLambda_PzVsPsi) {
1611-
if (fillingConfigs.isFillTHN_Pz)
1612-
histos.get<THn>(HIST("hLambdaPzVsPsi"))->Fill(coll.centFT0C(), chargeIndex, v0.pt(), v0.mLambda(), cosThetaLambda, 2 * lambdaminuspsiT0C);
1613-
if (fillingConfigs.isFillTHN_Acc)
1614-
histos.get<THn>(HIST("hLambdaCos2ThetaVsPsi"))->Fill(coll.centFT0C(), chargeIndex, v0.eta(), v0.pt(), v0.mLambda(), cos2ThetaLambda, 2 * lambdaminuspsiT0C);
1614+
if (fillingConfigs.isFillTHN_Pz)
1615+
histos.get<THn>(HIST("hLambdaPzVsPsi"))->Fill(coll.centFT0C(), chargeIndex, v0.pt(), v0.mLambda(), cosThetaLambda, 2 * lambdaminuspsiT0C);
1616+
if (fillingConfigs.isFillTHN_Acc)
1617+
histos.get<THn>(HIST("hLambdaCos2ThetaVsPsi"))->Fill(coll.centFT0C(), chargeIndex, v0.eta(), v0.pt(), v0.mLambda(), cos2ThetaLambda, 2 * lambdaminuspsiT0C);
16151618
}
16161619
}
16171620
}
@@ -1681,7 +1684,7 @@ struct cascadeFlow {
16811684
int counter = 0;
16821685
IsCascAccepted(casc, negExtra, posExtra, bachExtra, counter);
16831686
histos.fill(HIST("hCascade"), counter);
1684-
1687+
16851688
// ML selections
16861689
bool isSelectedCasc[nParticles]{false, false};
16871690

0 commit comments

Comments
 (0)