Skip to content

Commit 726a76e

Browse files
fuchuncuifuchuncui
authored
Add files via upload
1 parent 129eb94 commit 726a76e

File tree

1 file changed

+57
-39
lines changed

1 file changed

+57
-39
lines changed

‎PWGCF/Flow/Tasks/flowGFWOmegaXi.cxx‎

Lines changed: 57 additions & 39 deletions
Original file line numberDiff line numberDiff line change
@@ -143,10 +143,10 @@ struct FlowGFWOmegaXi {
143143
TF1* fT0AV0AMean = nullptr;
144144
TF1* fT0AV0ASigma = nullptr;
145145

146-
using TracksPID = soa::Join<aod::pidTPCFullPi, aod::pidTOFFullPi, aod::pidTPCFullKa, aod::pidTPCFullPr>;
146+
using TracksPID = soa::Join<aod::pidTPCPi, aod::pidTOFPi, aod::pidTPCKa, aod::pidTPCPr, aod::pidTOFKa, aod::pidTOFPr>;
147147
using aodTracks = soa::Filtered<soa::Join<aod::Tracks, aod::TrackSelection, aod::TracksExtra, TracksPID>>; // tracks filter
148148
using aodCollisions = soa::Filtered<soa::Join<aod::Collisions, aod::EvSels, aod::CentFT0Cs, aod::Mults>>; // collisions filter
149-
using DaughterTracks = soa::Join<aod::Tracks, aod::TracksExtra, aod::TracksDCA, aod::pidTPCPi, aod::pidTPCPr, aod::pidTPCKa>;
149+
using DaughterTracks = soa::Join<aod::Tracks, aod::TracksExtra, aod::TracksDCA, TracksPID>;
150150

151151
// Declare the pt, mult and phi Axis;
152152
int nPtBins = 0;
@@ -203,6 +203,8 @@ struct FlowGFWOmegaXi {
203203
// cumulant of flow
204204
registry.add("c22", ";Centrality (%) ; C_{2}{2} ", {HistType::kTProfile, {axisMultiplicity}});
205205
registry.add("c24", ";Centrality (%) ; C_{2}{2} ", {HistType::kTProfile, {axisMultiplicity}});
206+
registry.add("K0sc22", ";Centrality (%) ; C_{2}{2} ", {HistType::kTProfile, {axisMultiplicity}});
207+
registry.add("Lambdac22", ";Centrality (%) ; C_{2}{2} ", {HistType::kTProfile, {axisMultiplicity}});
206208
registry.add("c22dpt", ";Centrality (%) ; C_{2}{2}", {HistType::kTProfile2D, {cfgaxisPt, axisMultiplicity}});
207209
registry.add("c24dpt", ";Centrality (%) ; C_{2}{4}", {HistType::kTProfile2D, {cfgaxisPt, axisMultiplicity}});
208210
// pt-diff cumulant of flow
@@ -258,38 +260,50 @@ struct FlowGFWOmegaXi {
258260
fGFW->AddRegion("reffulldpt", -0.8, 0.8, nPtBins, 1);
259261
fGFW->AddRegion("refoldpt", -0.8, 0.8, nPtBins, 1);
260262
int nXiptMassBins = nXiPtBins * cfgXiMassbins;
261-
fGFW->AddRegion("poiXiP", 0.4, 0.8, nXiptMassBins, 2);
262-
fGFW->AddRegion("poiXiN", -0.8, -0.4, nXiptMassBins, 2);
263-
fGFW->AddRegion("poiXifull", -0.8, 0.8, nXiptMassBins, 2);
263+
fGFW->AddRegion("poiXiPdpt", 0.4, 0.8, nXiptMassBins, 2);
264+
fGFW->AddRegion("poiXiNdpt", -0.8, -0.4, nXiptMassBins, 2);
265+
fGFW->AddRegion("poiXifulldpt", -0.8, 0.8, nXiptMassBins, 2);
266+
fGFW->AddRegion("poiXiP", 0.4, 0.8, 1, 2);
267+
fGFW->AddRegion("poiXiN", -0.8, -0.4, 1, 2);
264268
int nOmegaptMassBins = nXiPtBins * cfgOmegaMassbins;
265-
fGFW->AddRegion("poiOmegaP", 0.4, 0.8, nOmegaptMassBins, 4);
266-
fGFW->AddRegion("poiOmegaN", -0.8, -0.4, nOmegaptMassBins, 4);
267-
fGFW->AddRegion("poiOmegafull", -0.8, 0.8, nOmegaptMassBins, 4);
269+
fGFW->AddRegion("poiOmegaPdpt", 0.4, 0.8, nOmegaptMassBins, 4);
270+
fGFW->AddRegion("poiOmegaNdpt", -0.8, -0.4, nOmegaptMassBins, 4);
271+
fGFW->AddRegion("poiOmegafulldpt", -0.8, 0.8, nOmegaptMassBins, 4);
272+
fGFW->AddRegion("poiOmegaP", 0.4, 0.8, 1, 4);
273+
fGFW->AddRegion("poiOmegaN", -0.8, -0.4, 1, 4);
268274
int nK0sptMassBins = nV0PtBins * cfgK0sMassbins;
269-
fGFW->AddRegion("poiK0sP", 0.4, 0.8, nK0sptMassBins, 8);
270-
fGFW->AddRegion("poiK0sN", -0.8, -0.4, nK0sptMassBins, 8);
271-
fGFW->AddRegion("poiK0sfull", -0.8, 0.8, nK0sptMassBins, 8);
275+
fGFW->AddRegion("poiK0sPdpt", 0.4, 0.8, nK0sptMassBins, 8);
276+
fGFW->AddRegion("poiK0sNdpt", -0.8, -0.4, nK0sptMassBins, 8);
277+
fGFW->AddRegion("poiK0sfulldpt", -0.8, 0.8, nK0sptMassBins, 8);
278+
fGFW->AddRegion("poiK0sP", 0.4, 0.8, 1, 8);
279+
fGFW->AddRegion("poiK0sN", -0.8, 0.4, 1, 8);
272280
int nLambdaptMassBins = nV0PtBins * cfgLambdaMassbins;
273-
fGFW->AddRegion("poiLambdaP", 0.4, 0.8, nLambdaptMassBins, 16);
274-
fGFW->AddRegion("poiLambdaN", -0.8, -0.4, nLambdaptMassBins, 16);
275-
fGFW->AddRegion("poiLambdafull", -0.8, 0.8, nLambdaptMassBins, 16);
281+
fGFW->AddRegion("poiLambdaPdpt", 0.4, 0.8, nLambdaptMassBins, 16);
282+
fGFW->AddRegion("poiLambdaNdpt", -0.8, -0.4, nLambdaptMassBins, 16);
283+
fGFW->AddRegion("poiLambdafulldpt", -0.8, 0.8, nLambdaptMassBins, 16);
284+
fGFW->AddRegion("poiLambdaP", 0.4, 0.8, 1, 16);
285+
fGFW->AddRegion("poiLambdaN", -0.8, -0.4, 1, 16);
276286
// pushback
277287
corrconfigs.push_back(fGFW->GetCorrelatorConfig("refP10dpt {2} refN10dpt {-2}", "Ref10Gap22dpt", kTRUE));
278288
corrconfigs.push_back(fGFW->GetCorrelatorConfig("reffulldpt reffulldpt {2 2 -2 -2}", "Ref10Gap24dpt", kTRUE));
279-
corrconfigs.push_back(fGFW->GetCorrelatorConfig("poiXiP {2} refN10 {-2}", "Xi10Gap22a", kTRUE));
280-
corrconfigs.push_back(fGFW->GetCorrelatorConfig("poiXiN {2} refP10 {-2}", "Xi10Gap22b", kTRUE));
281-
corrconfigs.push_back(fGFW->GetCorrelatorConfig("poiXifull reffull {2 2 -2 -2}", "Xi10Gap24", kTRUE));
282-
corrconfigs.push_back(fGFW->GetCorrelatorConfig("poiOmegaP {2} refN10 {-2}", "Omega10Gap22a", kTRUE));
283-
corrconfigs.push_back(fGFW->GetCorrelatorConfig("poiOmegaN {2} refP10 {-2}", "Omega10Gap22b", kTRUE));
284-
corrconfigs.push_back(fGFW->GetCorrelatorConfig("poiOmegaP reffull {2 2 -2 -2}", "Xi10Gap24", kTRUE));
285-
corrconfigs.push_back(fGFW->GetCorrelatorConfig("poiK0sP {2} refN10 {-2}", "K0short10Gap22a", kTRUE));
286-
corrconfigs.push_back(fGFW->GetCorrelatorConfig("poiK0sN {2} refP10 {-2}", "K0short10Gap22b", kTRUE));
287-
corrconfigs.push_back(fGFW->GetCorrelatorConfig("poiK0sP reffull {2 2 -2 -2}", "Xi10Gap24", kTRUE));
288-
corrconfigs.push_back(fGFW->GetCorrelatorConfig("poiLambdaP {2} refN10 {-2}", "Lambda10Gap22a", kTRUE));
289-
corrconfigs.push_back(fGFW->GetCorrelatorConfig("poiLambdaN {2} refP10 {-2}", "Lambda10Gap22b", kTRUE));
290-
corrconfigs.push_back(fGFW->GetCorrelatorConfig("poiLambdaP reffull {2 2 -2 -2}", "Xi10Gap24a", kTRUE));
289+
corrconfigs.push_back(fGFW->GetCorrelatorConfig("poiXiPdpt {2} refN10 {-2}", "Xi10Gap22a", kTRUE));
290+
corrconfigs.push_back(fGFW->GetCorrelatorConfig("poiXiNdpt {2} refP10 {-2}", "Xi10Gap22b", kTRUE));
291+
corrconfigs.push_back(fGFW->GetCorrelatorConfig("poiXifulldpt reffull {2 2 -2 -2}", "Xi10Gap24", kTRUE));
292+
corrconfigs.push_back(fGFW->GetCorrelatorConfig("poiOmegaPdpt {2} refN10 {-2}", "Omega10Gap22a", kTRUE));
293+
corrconfigs.push_back(fGFW->GetCorrelatorConfig("poiOmegaNdpt {2} refP10 {-2}", "Omega10Gap22b", kTRUE));
294+
corrconfigs.push_back(fGFW->GetCorrelatorConfig("poiOmegaPdpt reffull {2 2 -2 -2}", "Xi10Gap24", kTRUE));
295+
corrconfigs.push_back(fGFW->GetCorrelatorConfig("poiK0sPdpt {2} refN10 {-2}", "K0short10Gap22a", kTRUE));
296+
corrconfigs.push_back(fGFW->GetCorrelatorConfig("poiK0sNdpt {2} refP10 {-2}", "K0short10Gap22b", kTRUE));
297+
corrconfigs.push_back(fGFW->GetCorrelatorConfig("poiK0sPdpt reffull {2 2 -2 -2}", "Xi10Gap24", kTRUE));
298+
corrconfigs.push_back(fGFW->GetCorrelatorConfig("poiLambdaPdpt {2} refN10 {-2}", "Lambda10Gap22a", kTRUE));
299+
corrconfigs.push_back(fGFW->GetCorrelatorConfig("poiLambdaNdpt {2} refP10 {-2}", "Lambda10Gap22b", kTRUE));
300+
corrconfigs.push_back(fGFW->GetCorrelatorConfig("poiLambdaPdpt reffull {2 2 -2 -2}", "Xi10Gap24a", kTRUE));
291301
corrconfigs.push_back(fGFW->GetCorrelatorConfig("refP10 {2} refN10 {-2}", "Ref10Gap22", kFALSE));
292302
corrconfigs.push_back(fGFW->GetCorrelatorConfig("reffull reffull {2 2 -2 -2}", "Ref10Gap24", kFALSE));
303+
corrconfigs.push_back(fGFW->GetCorrelatorConfig("poiK0sP {2} refN10 {-2}", "K0s10Gap22inta", kFALSE));
304+
corrconfigs.push_back(fGFW->GetCorrelatorConfig("poiK0sN {2} refP10 {-2}", "K0s10Gap22intb", kFALSE));
305+
corrconfigs.push_back(fGFW->GetCorrelatorConfig("poiLambdaP {2} refN10 {-2}", "Lambda10Gap22inta", kFALSE));
306+
corrconfigs.push_back(fGFW->GetCorrelatorConfig("poiLambdaN {2} refP10 {-2}", "Lambda10Gap22intb", kFALSE));
293307
fGFW->CreateRegions(); // finalize the initialization
294308

295309
// used for event selection
@@ -342,10 +356,10 @@ struct FlowGFWOmegaXi {
342356
{
343357
float dnx = 0;
344358
float val = 0;
345-
dnx = fGFW->Calculate(corrconf, ptbin, kTRUE).real();
359+
dnx = fGFW->Calculate(corrconf, ptbin - 1, kTRUE).real();
346360
if (dnx == 0)
347361
return;
348-
val = fGFW->Calculate(corrconf, ptbin, kFALSE).real() / dnx;
362+
val = fGFW->Calculate(corrconf, ptbin - 1, kFALSE).real() / dnx;
349363
if (TMath::Abs(val) < 1) {
350364
registry.fill(tarName, fPtAxis->GetBinCenter(ptbin), cent, val, dnx);
351365
}
@@ -528,7 +542,7 @@ struct FlowGFWOmegaXi {
528542
continue;
529543
registry.fill(HIST("hPhi"), track.phi());
530544
registry.fill(HIST("hEta"), track.eta());
531-
registry.fill(HIST("hEtaPhiVtxzREF"), track.phi(), track.eta(), vtxz);
545+
registry.fill(HIST("hEtaPhiVtxzREF"), track.phi(), track.eta(), vtxz, wacc);
532546
registry.fill(HIST("hPt"), track.pt());
533547
int ptbin = fPtAxis->FindBin(track.pt()) - 1;
534548
if ((track.pt() > cfgCutPtMin) && (track.pt() < cfgCutPtMax)) {
@@ -543,13 +557,13 @@ struct FlowGFWOmegaXi {
543557
auto v0negdau = v0.negTrack_as<DaughterTracks>();
544558
// check tpc
545559
int PDGCode = 0;
546-
if (v0.qtarm() / TMath::Abs(v0.alpha()) > cfgv0_ArmPodocut && TMath::Abs(v0posdau.tpcNSigmaPi()) < cfgNSigmaCascPion && TMath::Abs(v0negdau.tpcNSigmaPi()) < cfgNSigmaCascPion) {
560+
if (v0.qtarm() / TMath::Abs(v0.alpha()) > cfgv0_ArmPodocut && TMath::Abs(v0posdau.tpcNSigmaPi()) < cfgNSigmaCascPion && TMath::Abs(v0negdau.tpcNSigmaPi()) < cfgNSigmaCascPion && TMath::Abs(v0posdau.tofNSigmaPi()) < cfgNSigmaCascPion && TMath::Abs(v0negdau.tofNSigmaPi()) < cfgNSigmaCascPion) {
547561
registry.fill(HIST("InvMassK0s_all"), v0.pt(), v0.mK0Short(), v0.eta(), cent);
548562
if (!setCurrentParticleWeights(weff, wacc, v0, vtxz, 1))
549563
continue;
550564
PDGCode = kK0Short;
551565
CandNum_all[0] = CandNum_all[0] + 1;
552-
} else if (v0.qtarm() / TMath::Abs(v0.alpha()) < cfgv0_ArmPodocut && TMath::Abs(v0posdau.tpcNSigmaPr()) < cfgNSigmaCascProton && TMath::Abs(v0negdau.tpcNSigmaPi()) < cfgNSigmaCascPion) {
566+
} else if (TMath::Abs(v0posdau.tpcNSigmaPr()) < cfgNSigmaCascProton && TMath::Abs(v0negdau.tpcNSigmaPi()) < cfgNSigmaCascPion && TMath::Abs(v0posdau.tofNSigmaPr()) < cfgNSigmaCascProton && TMath::Abs(v0negdau.tofNSigmaPi()) < cfgNSigmaCascPion) {
553567
registry.fill(HIST("InvMassLambda_all"), v0.pt(), v0.mLambda(), v0.eta(), cent);
554568
if (!setCurrentParticleWeights(weff, wacc, v0, vtxz, 2))
555569
continue;
@@ -576,7 +590,7 @@ struct FlowGFWOmegaXi {
576590
if (TMath::Abs(v0.mK0Short() - o2::constants::physics::MassK0Short) < cfgv0_mk0swindow) {
577591
CandNum[0] = CandNum[0] + 1;
578592
registry.fill(HIST("InvMassK0s"), v0.pt(), v0.mK0Short(), v0.eta(), cent);
579-
registry.fill(HIST("hEtaPhiVtxzPOIK0s"), v0.phi(), v0.eta(), vtxz);
593+
registry.fill(HIST("hEtaPhiVtxzPOIK0s"), v0.phi(), v0.eta(), vtxz, wacc);
580594
fGFW->Fill(v0.eta(), fV0PtAxis->FindBin(v0.pt()) - 1 + ((fK0sMass->FindBin(v0.mK0Short()) - 1) * nV0PtBins), v0.phi(), wacc * weff, 8);
581595
if (cfgOutputNUAWeights)
582596
fWeightsK0s->Fill(v0.phi(), v0.eta(), vtxz, v0.pt(), cent, 0);
@@ -585,7 +599,7 @@ struct FlowGFWOmegaXi {
585599
if (TMath::Abs(v0.mLambda() - o2::constants::physics::MassLambda0) < cfgv0_mlambdawindow) {
586600
CandNum[1] = CandNum[1] + 1;
587601
registry.fill(HIST("InvMassLambda"), v0.pt(), v0.mLambda(), v0.eta(), cent);
588-
registry.fill(HIST("hEtaPhiVtxzPOILambda"), v0.phi(), v0.eta(), vtxz);
602+
registry.fill(HIST("hEtaPhiVtxzPOILambda"), v0.phi(), v0.eta(), vtxz, wacc);
589603
fGFW->Fill(v0.eta(), fV0PtAxis->FindBin(v0.pt()) - 1 + ((fLambdaMass->FindBin(v0.mLambda()) - 1) * nV0PtBins), v0.phi(), wacc * weff, 16);
590604
if (cfgOutputNUAWeights)
591605
fWeightsLambda->Fill(v0.phi(), v0.eta(), vtxz, v0.pt(), cent, 0);
@@ -602,13 +616,13 @@ struct FlowGFWOmegaXi {
602616
continue;
603617
}
604618
int PDGCode = 0;
605-
if (casc.sign() < 0 && TMath::Abs(casc.yOmega()) < cfgCasc_rapidity && TMath::Abs(bachelor.tpcNSigmaKa()) < cfgNSigmaCascKaon && TMath::Abs(posdau.tpcNSigmaPr()) < cfgNSigmaCascProton && TMath::Abs(negdau.tpcNSigmaPi()) < cfgNSigmaCascPion) {
619+
if (casc.sign() < 0 && TMath::Abs(casc.yOmega()) < cfgCasc_rapidity && TMath::Abs(bachelor.tpcNSigmaKa()) < cfgNSigmaCascKaon && TMath::Abs(posdau.tpcNSigmaPr()) < cfgNSigmaCascProton && TMath::Abs(negdau.tpcNSigmaPi()) < cfgNSigmaCascPion && TMath::Abs(bachelor.tofNSigmaKa()) < cfgNSigmaCascKaon && TMath::Abs(posdau.tofNSigmaPr()) < cfgNSigmaCascProton && TMath::Abs(negdau.tofNSigmaPi()) < cfgNSigmaCascPion) {
606620
registry.fill(HIST("InvMassOmegaMinus_all"), casc.pt(), casc.mOmega(), casc.eta(), cent);
607621
if (!setCurrentParticleWeights(weff, wacc, casc, vtxz, 4))
608622
continue;
609623
PDGCode = kOmegaMinus;
610624
CandNum_all[3] = CandNum_all[3] + 1;
611-
} else if (casc.sign() < 0 && TMath::Abs(casc.yXi()) < cfgCasc_rapidity && TMath::Abs(bachelor.tpcNSigmaPi()) < cfgNSigmaCascPion && TMath::Abs(posdau.tpcNSigmaPr()) < cfgNSigmaCascProton && TMath::Abs(negdau.tpcNSigmaPi()) < cfgNSigmaCascPion) {
625+
} else if (casc.sign() < 0 && TMath::Abs(casc.yXi()) < cfgCasc_rapidity && TMath::Abs(bachelor.tpcNSigmaPi()) < cfgNSigmaCascPion && TMath::Abs(posdau.tpcNSigmaPr()) < cfgNSigmaCascProton && TMath::Abs(negdau.tpcNSigmaPi() ) < cfgNSigmaCascPion && TMath::Abs(bachelor.tofNSigmaPi()) < cfgNSigmaCascPion && TMath::Abs(posdau.tofNSigmaPr()) < cfgNSigmaCascProton && TMath::Abs(negdau.tofNSigmaPi() ) < cfgNSigmaCascPion) {
612626
registry.fill(HIST("InvMassXiMinus_all"), casc.pt(), casc.mXi(), casc.eta(), cent);
613627
if (!setCurrentParticleWeights(weff, wacc, casc, vtxz, 3))
614628
continue;
@@ -649,7 +663,7 @@ struct FlowGFWOmegaXi {
649663
continue;
650664
if (PDGCode == kOmegaMinus) {
651665
CandNum[3] = CandNum[3] + 1;
652-
registry.fill(HIST("hEtaPhiVtxzPOIOmega"), casc.phi(), casc.eta(), vtxz);
666+
registry.fill(HIST("hEtaPhiVtxzPOIOmega"), casc.phi(), casc.eta(), vtxz, wacc);
653667
registry.fill(HIST("InvMassOmegaMinus"), casc.pt(), casc.mOmega(), casc.eta(), cent);
654668
if ((casc.pt() < cfgCutPtPOIMax) && (casc.pt() > cfgCutPtPOIMin) && (casc.mOmega() > 1.63) && (casc.mOmega() < 1.71)) {
655669
fGFW->Fill(casc.eta(), fXiPtAxis->FindBin(casc.pt()) - 1 + ((fOmegaMass->FindBin(casc.mOmega()) - 1) * nXiPtBins), casc.phi(), wacc * weff, 4);
@@ -659,7 +673,7 @@ struct FlowGFWOmegaXi {
659673
}
660674
if (PDGCode == kXiMinus) {
661675
CandNum[2] = CandNum[2] + 1;
662-
registry.fill(HIST("hEtaPhiVtxzPOIXi"), casc.phi(), casc.eta(), vtxz);
676+
registry.fill(HIST("hEtaPhiVtxzPOIXi"), casc.phi(), casc.eta(), vtxz, wacc);
663677
registry.fill(HIST("InvMassXiMinus"), casc.pt(), casc.mXi(), casc.eta(), cent);
664678
if ((casc.pt() < cfgCutPtPOIMax) && (casc.pt() > cfgCutPtPOIMin) && (casc.mXi() > 1.30) && (casc.mXi() < 1.37)) {
665679
fGFW->Fill(casc.eta(), fXiPtAxis->FindBin(casc.pt()) - 1 + ((fXiMass->FindBin(casc.mXi()) - 1) * nXiPtBins), casc.phi(), wacc * weff, 2);
@@ -669,16 +683,20 @@ struct FlowGFWOmegaXi {
669683
}
670684
}
671685
for (int i = 0; i < 4; i++) {
672-
if (CandNum_all[i] > 1) {
686+
if (CandNum_all[i] > 0) {
673687
registry.fill(HIST("hEventCount"), 2.5, i + 0.5);
674688
}
675-
if (CandNum[i] > 1) {
689+
if (CandNum[i] > 0) {
676690
registry.fill(HIST("hEventCount"), 3.5, i + 0.5);
677691
}
678692
}
679693
// Filling cumulant with ROOT TProfile and loop for all ptBins
680694
FillProfile(corrconfigs.at(14), HIST("c22"), cent);
681695
FillProfile(corrconfigs.at(15), HIST("c24"), cent);
696+
FillProfile(corrconfigs.at(16), HIST("K0sc22"), cent);
697+
FillProfile(corrconfigs.at(17), HIST("K0sc22"), cent);
698+
FillProfile(corrconfigs.at(18), HIST("Lambdac22"), cent);
699+
FillProfile(corrconfigs.at(19), HIST("Lambdac22"), cent);
682700
for (int i = 1; i <= nPtBins; i++) {
683701
FillProfilepT(corrconfigs.at(0), HIST("c22dpt"), i, cent);
684702
FillProfilepT(corrconfigs.at(1), HIST("c24dpt"), i, cent);

0 commit comments

Comments
 (0)