Skip to content

Commit fb43714

Browse files
author
Prottay Das
committed
detailed QA plots separately for lambda and antilambda
1 parent 23c1fcf commit fb43714

File tree

1 file changed

+43
-15
lines changed

1 file changed

+43
-15
lines changed

‎PWGLF/Tasks/Strangeness/lambdapolsp.cxx‎

Lines changed: 43 additions & 15 deletions
Original file line numberDiff line numberDiff line change
@@ -206,7 +206,11 @@ struct lambdapolsp {
206206
ConfigurableAxis axisDca{"axisDca", {100, -5, 5}, "dca axis"};
207207
ConfigurableAxis axisLT{"axisLT", {50, 0, 50}, "lifetime axis"};
208208
ConfigurableAxis axisCR{"axisCR", {40, 0, 200}, "CR axis"};
209+
ConfigurableAxis axisnsig{"axisnsig", {100, -10, 10}, "nsigma axis"};
210+
ConfigurableAxis axispt{"axispt", {100, 0, 10}, "pt axis"};
209211
Configurable<bool> filldist{"filldist", true, "fill topo distr"};
212+
Configurable<float> lowmasscut{"lowmasscut", 1.11, "low mass cut"};
213+
Configurable<float> highmasscut{"highmasscut", 1.12, "high mass cut"};
210214
} distGrp;
211215

212216
RCTFlagsChecker rctChecker;
@@ -399,15 +403,29 @@ struct lambdapolsp {
399403
}
400404

401405
if (distGrp.filldist) {
402-
histos.add("hcosine", "hcosine", HistType::kTH1D, {distGrp.axisCosine}, true);
403-
histos.add("hdcabwv0daugh", "hdcabwv0daugh", HistType::kTH1D, {distGrp.axisDca}, true);
404-
histos.add("hlifetime", "hlifetime", HistType::kTH1D, {distGrp.axisLT}, true);
405-
histos.add("hradius", "hradius", HistType::kTH1D, {distGrp.axisRadius}, true);
406+
histos.add("hcosinelambda", "hcosinelambda", HistType::kTH1D, {distGrp.axisCosine}, true);
407+
histos.add("hdcabwv0daughlambda", "hdcabwv0daughlambda", HistType::kTH1D, {distGrp.axisDca}, true);
408+
histos.add("hlifetimelambda", "hlifetimelambda", HistType::kTH1D, {distGrp.axisLT}, true);
409+
histos.add("hradiuslambda", "hradiuslambda", HistType::kTH1D, {distGrp.axisRadius}, true);
406410
histos.add("hdcaposlambda", "hdcaposlambda", HistType::kTH1D, {distGrp.axisDca}, true);
407411
histos.add("hdcaneglambda", "hdcaneglambda", HistType::kTH1D, {distGrp.axisDca}, true);
412+
histos.add("htpcCRlambda", "htpcCRlambda", HistType::kTH1D, {distGrp.axisCR}, true);
413+
histos.add("htpcposlambda", "htpcposlambda", HistType::kTH1D, {distGrp.axisnsig}, true);
414+
histos.add("htpcneglambda", "htpcneglambda", HistType::kTH1D, {distGrp.axisnsig}, true);
415+
histos.add("hptposlambda", "hptposlambda", HistType::kTH1D, {distGrp.axispt}, true);
416+
histos.add("hptneglambda", "hptneglambda", HistType::kTH1D, {distGrp.axispt}, true);
417+
418+
histos.add("hcosineantilambda", "hcosineantilambda", HistType::kTH1D, {distGrp.axisCosine}, true);
419+
histos.add("hdcabwv0daughantilambda", "hdcabwv0daughantilambda", HistType::kTH1D, {distGrp.axisDca}, true);
420+
histos.add("hlifetimeantilambda", "hlifetimeantilambda", HistType::kTH1D, {distGrp.axisLT}, true);
421+
histos.add("hradiusantilambda", "hradiusantilambda", HistType::kTH1D, {distGrp.axisRadius}, true);
408422
histos.add("hdcaposantilambda", "hdcaposantilambda", HistType::kTH1D, {distGrp.axisDca}, true);
409423
histos.add("hdcanegantilambda", "hdcanegantilambda", HistType::kTH1D, {distGrp.axisDca}, true);
410-
histos.add("htpcCR", "htpcCR", HistType::kTH1D, {distGrp.axisCR}, true);
424+
histos.add("htpcCRantilambda", "htpcCRantilambda", HistType::kTH1D, {distGrp.axisCR}, true);
425+
histos.add("htpcposantilambda", "htpcposantilambda", HistType::kTH1D, {distGrp.axisnsig}, true);
426+
histos.add("htpcnegantilambda", "htpcnegantilambda", HistType::kTH1D, {distGrp.axisnsig}, true);
427+
histos.add("hptposantilambda", "hptposantilambda", HistType::kTH1D, {distGrp.axispt}, true);
428+
histos.add("hptnegantilambda", "hptnegantilambda", HistType::kTH1D, {distGrp.axispt}, true);
411429
}
412430

413431
ccdb->setURL(cfgCcdbParam.cfgURL);
@@ -1168,23 +1186,24 @@ struct lambdapolsp {
11681186
}
11691187
} else {
11701188

1171-
if (distGrp.filldist) {
1172-
histos.fill(HIST("hcosine"), v0.v0cosPA());
1173-
histos.fill(HIST("hdcabwv0daugh"), v0.dcaV0daughters());
1174-
histos.fill(HIST("hlifetime"), TMath::Abs(v0.distovertotmom(collision.posX(), collision.posY(), collision.posZ()) * massLambda));
1175-
histos.fill(HIST("hradius"), v0.v0radius());
1176-
histos.fill(HIST("htpcCR"), postrack.tpcNClsCrossedRows());
1177-
}
1178-
11791189
if (LambdaTag) {
11801190
Lambda = Proton + AntiPion;
11811191
tagb = 0;
11821192
int binx = accprofileL->GetXaxis()->FindBin(v0.eta());
11831193
int biny = accprofileL->GetYaxis()->FindBin(v0.pt());
11841194
double acvalue = accprofileL->GetBinContent(binx, biny);
1185-
if (distGrp.filldist) {
1195+
if (distGrp.filldist && aLambdaTag == 0 && Lambda.M() > distGrp.lowmasscut && Lambda.M() < distGrp.highmasscut) {
1196+
histos.fill(HIST("hcosinelambda"), v0.v0cosPA());
1197+
histos.fill(HIST("hdcabwv0daughlambda"), v0.dcaV0daughters());
1198+
histos.fill(HIST("hlifetimelambda"), TMath::Abs(v0.distovertotmom(collision.posX(), collision.posY(), collision.posZ()) * massLambda));
1199+
histos.fill(HIST("hradiuslambda"), v0.v0radius());
1200+
histos.fill(HIST("htpcCRlambda"), postrack.tpcNClsCrossedRows());
11861201
histos.fill(HIST("hdcaposlambda"), v0.dcapostopv());
11871202
histos.fill(HIST("hdcaneglambda"), v0.dcanegtopv());
1203+
histos.fill(HIST("htpcposlambda"), postrack.tpcNSigmaPr());
1204+
histos.fill(HIST("htpcneglambda"), negtrack.tpcNSigmaPi());
1205+
histos.fill(HIST("hptposlambda"), Proton.Pt());
1206+
histos.fill(HIST("hptneglambda"), AntiPion.Pt());
11881207
}
11891208
fillHistograms(taga, tagb, Lambda, Proton, psiZDCC, psiZDCA, psiZDC, centrality, v0.mLambda(), v0.pt(), v0.eta(), acvalue, 1.0);
11901209
}
@@ -1196,9 +1215,18 @@ struct lambdapolsp {
11961215
int binx = accprofileAL->GetXaxis()->FindBin(v0.eta());
11971216
int biny = accprofileAL->GetYaxis()->FindBin(v0.pt());
11981217
double acvalue = accprofileAL->GetBinContent(binx, biny);
1199-
if (distGrp.filldist) {
1218+
if (distGrp.filldist && LambdaTag == 0 && AntiLambda.M() > distGrp.lowmasscut && AntiLambda.M() < distGrp.highmasscut) {
1219+
histos.fill(HIST("hcosineantilambda"), v0.v0cosPA());
1220+
histos.fill(HIST("hdcabwv0daughantilambda"), v0.dcaV0daughters());
1221+
histos.fill(HIST("hlifetimeantilambda"), TMath::Abs(v0.distovertotmom(collision.posX(), collision.posY(), collision.posZ()) * massLambda));
1222+
histos.fill(HIST("hradiusantilambda"), v0.v0radius());
1223+
histos.fill(HIST("htpcCRantilambda"), postrack.tpcNClsCrossedRows());
12001224
histos.fill(HIST("hdcaposantilambda"), v0.dcapostopv());
12011225
histos.fill(HIST("hdcanegantilambda"), v0.dcanegtopv());
1226+
histos.fill(HIST("htpcposantilambda"), postrack.tpcNSigmaPi());
1227+
histos.fill(HIST("htpcnegantilambda"), negtrack.tpcNSigmaPr());
1228+
histos.fill(HIST("hptposantilambda"), Pion.Pt());
1229+
histos.fill(HIST("hptnegantilambda"), AntiProton.Pt());
12021230
}
12031231
fillHistograms(taga, tagb, AntiLambda, AntiProton, psiZDCC, psiZDCA, psiZDC, centrality, v0.mAntiLambda(), v0.pt(), v0.eta(), acvalue, wgtvalue);
12041232
}

0 commit comments

Comments
 (0)