Skip to content

Commit 7b55421

Browse files
author
Prottay Das
committed
updated
1 parent aaadf6a commit 7b55421

File tree

1 file changed

+48
-20
lines changed

1 file changed

+48
-20
lines changed

‎PWGLF/Tasks/Strangeness/lambdapolsp.cxx‎

Lines changed: 48 additions & 20 deletions
Original file line numberDiff line numberDiff line change
@@ -875,6 +875,7 @@ struct lambdapolsp {
875875
if (rctCut.requireRCTFlagChecker && !rctChecker(collision)) {
876876
return;
877877
}
878+
878879
// currentRunNumber = collision.foundBC_as<BCsRun3>().runNumber();
879880
auto bc = collision.foundBC_as<BCsRun3>();
880881

@@ -902,7 +903,7 @@ struct lambdapolsp {
902903
modqyZDCC = qyZDCC;
903904
}
904905

905-
auto psiZDC = TMath::ATan2((modqyZDCC - modqyZDCA), (modqxZDCC - modqxZDCA)); // full event plane
906+
auto psiZDC = TMath::ATan2((modqyZDCC - modqyZDCA), (modqxZDCC - modqxZDCA)); // full event plane*/
906907
/*if (useonlypsis) {
907908
psiZDC = psiZDCC - psiZDCA;
908909
}*/
@@ -1123,6 +1124,7 @@ struct lambdapolsp {
11231124
accprofileL = ccdb->getForTimeStamp<TProfile2D>(ConfAccPathL.value, bc.timestamp());
11241125
accprofileAL = ccdb->getForTimeStamp<TProfile2D>(ConfAccPathAL.value, bc.timestamp());
11251126
}
1127+
double acvalue = 1.0;
11261128
int binxwgt;
11271129
double wgtvalue;
11281130
if (useyldwgt) {
@@ -1186,30 +1188,42 @@ struct lambdapolsp {
11861188
if (LambdaTag) {
11871189
Lambda = Proton + AntiPion;
11881190
tagb = 0;
1189-
int binx = accprofileL->GetXaxis()->FindBin(v0.eta());
1190-
int biny = accprofileL->GetYaxis()->FindBin(v0.pt());
1191-
double acvalue = accprofileL->GetBinContent(binx, biny);
1191+
if (useAccCorr) {
1192+
int binx = accprofileL->GetXaxis()->FindBin(v0.eta());
1193+
int biny = accprofileL->GetYaxis()->FindBin(v0.pt());
1194+
acvalue = accprofileL->GetBinContent(binx, biny);
1195+
} else {
1196+
acvalue = 1.0;
1197+
}
1198+
11921199
fillHistograms(taga, tagb, Lambda, Proton, psiZDCC, psiZDCA, psiZDC, centrality, v0.mLambda(), v0.pt(), desbinvalue, acvalue, 1.0);
11931200
}
11941201

11951202
tagb = aLambdaTag;
11961203
if (aLambdaTag) {
11971204
AntiLambda = AntiProton + Pion;
11981205
taga = 0;
1199-
int binx = accprofileAL->GetXaxis()->FindBin(v0.eta());
1200-
int biny = accprofileAL->GetYaxis()->FindBin(v0.pt());
1201-
double acvalue = accprofileAL->GetBinContent(binx, biny);
1206+
if (useAccCorr) {
1207+
int binx = accprofileL->GetXaxis()->FindBin(v0.eta());
1208+
int biny = accprofileL->GetYaxis()->FindBin(v0.pt());
1209+
acvalue = accprofileL->GetBinContent(binx, biny);
1210+
} else {
1211+
acvalue = 1.0;
1212+
}
12021213
fillHistograms(taga, tagb, AntiLambda, AntiProton, psiZDCC, psiZDCA, psiZDC, centrality, v0.mAntiLambda(), v0.pt(), desbinvalue, acvalue, 1.0);
12031214
}
12041215
}
12051216
} else {
1206-
12071217
if (LambdaTag) {
12081218
Lambda = Proton + AntiPion;
12091219
tagb = 0;
1210-
int binx = accprofileL->GetXaxis()->FindBin(v0.eta());
1211-
int biny = accprofileL->GetYaxis()->FindBin(v0.pt());
1212-
double acvalue = accprofileL->GetBinContent(binx, biny);
1220+
if (useAccCorr) {
1221+
int binx = accprofileL->GetXaxis()->FindBin(v0.eta());
1222+
int biny = accprofileL->GetYaxis()->FindBin(v0.pt());
1223+
acvalue = accprofileL->GetBinContent(binx, biny);
1224+
} else {
1225+
acvalue = 1.0;
1226+
}
12131227
if (distGrp.filldist && aLambdaTag == 0 && Lambda.M() > distGrp.lowmasscut && Lambda.M() < distGrp.highmasscut) {
12141228
histos.fill(HIST("hcosinelambda"), v0.v0cosPA());
12151229
histos.fill(HIST("hdcabwv0daughlambda"), v0.dcaV0daughters());
@@ -1230,9 +1244,13 @@ struct lambdapolsp {
12301244
if (aLambdaTag) {
12311245
AntiLambda = AntiProton + Pion;
12321246
taga = 0;
1233-
int binx = accprofileAL->GetXaxis()->FindBin(v0.eta());
1234-
int biny = accprofileAL->GetYaxis()->FindBin(v0.pt());
1235-
double acvalue = accprofileAL->GetBinContent(binx, biny);
1247+
if (useAccCorr) {
1248+
int binx = accprofileL->GetXaxis()->FindBin(v0.eta());
1249+
int biny = accprofileL->GetYaxis()->FindBin(v0.pt());
1250+
acvalue = accprofileL->GetBinContent(binx, biny);
1251+
} else {
1252+
acvalue = 1.0;
1253+
}
12361254
if (distGrp.filldist && LambdaTag == 0 && AntiLambda.M() > distGrp.lowmasscut && AntiLambda.M() < distGrp.highmasscut) {
12371255
histos.fill(HIST("hcosineantilambda"), v0.v0cosPA());
12381256
histos.fill(HIST("hdcabwv0daughantilambda"), v0.dcaV0daughters());
@@ -1500,9 +1518,14 @@ struct lambdapolsp {
15001518
if (analyzeLambda && LambdaTag) {
15011519
Lambda = Proton + AntiPion;
15021520
tagb = 0;
1503-
int binx = accprofileL->GetXaxis()->FindBin(v0.eta());
1504-
int biny = accprofileL->GetYaxis()->FindBin(v0.pt());
1505-
double acvalue = accprofileL->GetBinContent(binx, biny);
1521+
double acvalue = 1.0;
1522+
if (useAccCorr) {
1523+
int binx = accprofileL->GetXaxis()->FindBin(v0.eta());
1524+
int biny = accprofileL->GetYaxis()->FindBin(v0.pt());
1525+
acvalue = accprofileL->GetBinContent(binx, biny);
1526+
} else {
1527+
acvalue = 1.0;
1528+
}
15061529
// double acvalue = 1.0;
15071530
fillHistograms(taga, tagb, Lambda, Proton, psiZDCC, psiZDCA, psiZDC, centrality, v0.mLambda(), v0.pt(), v0.eta(), acvalue, 1.0);
15081531
}
@@ -1511,9 +1534,14 @@ struct lambdapolsp {
15111534
if (analyzeLambda && aLambdaTag) {
15121535
AntiLambda = AntiProton + Pion;
15131536
taga = 0;
1514-
int binx = accprofileAL->GetXaxis()->FindBin(v0.eta());
1515-
int biny = accprofileAL->GetYaxis()->FindBin(v0.pt());
1516-
double acvalue = accprofileAL->GetBinContent(binx, biny);
1537+
double acvalue = 1.0;
1538+
if (useAccCorr) {
1539+
int binx = accprofileAL->GetXaxis()->FindBin(v0.eta());
1540+
int biny = accprofileAL->GetYaxis()->FindBin(v0.pt());
1541+
acvalue = accprofileAL->GetBinContent(binx, biny);
1542+
} else {
1543+
acvalue = 1.0;
1544+
}
15171545
// double acvalue = 1.0;
15181546
fillHistograms(taga, tagb, AntiLambda, AntiProton, psiZDCC, psiZDCA, psiZDC, centrality, v0.mAntiLambda(), v0.pt(), v0.eta(), acvalue, wgtvalue);
15191547
}

0 commit comments

Comments
 (0)