Skip to content

Commit 0510b42

Browse files
fuchuncuifuchuncui
authored
[PWGCF] Updated event selection and particle reconstruction in flowGFWOmegaXi.cxx
1 parent 66bce9e commit 0510b42

File tree

1 file changed

+30
-30
lines changed

1 file changed

+30
-30
lines changed

‎PWGCF/Flow/Tasks/flowGFWOmegaXi.cxx‎

Lines changed: 30 additions & 30 deletions
Original file line numberDiff line numberDiff line change
@@ -167,10 +167,10 @@ struct FlowGFWOmegaXi {
167167
TF1* fT0AV0AMean = nullptr;
168168
TF1* fT0AV0ASigma = nullptr;
169169

170-
using tracksPID = soa::Join<aod::pidTPCPi, aod::pidTPCKa, aod::pidTPCPr, aod::pidTOFPi, aod::pidTOFKa, aod::pidTOFPr>;
171-
using aodTracks = soa::Filtered<soa::Join<aod::Tracks, aod::TrackSelection, aod::TracksExtra, tracksPID>>; // tracks filter
172-
using aodCollisions = soa::Filtered<soa::Join<aod::Collisions, aod::EvSels, aod::CentFT0Cs, aod::Mults>>; // collisions filter
173-
using daughterTracks = soa::Join<aod::Tracks, aod::TracksExtra, aod::TracksDCA, tracksPID>;
170+
using TracksPID = soa::Join<aod::pidTPCPi, aod::pidTPCKa, aod::pidTPCPr, aod::pidTOFPi, aod::pidTOFKa, aod::pidTOFPr>;
171+
using AodTracks = soa::Filtered<soa::Join<aod::Tracks, aod::TrackSelection, aod::TracksExtra, TracksPID>>; // tracks filter
172+
using AodCollisions = soa::Filtered<soa::Join<aod::Collisions, aod::EvSels, aod::CentFT0Cs, aod::Mults>>; // collisions filter
173+
using DaughterTracks = soa::Join<aod::Tracks, aod::TracksExtra, aod::TracksDCA, TracksPID>;
174174

175175
// Declare the pt, mult and phi Axis;
176176
int nPtBins = 0;
@@ -636,15 +636,15 @@ struct FlowGFWOmegaXi {
636636
return true;
637637
}
638638

639-
void process(aodCollisions::iterator const& collision, aod::BCsWithTimestamps const&, aodTracks const& tracks, aod::CascDataExt const& Cascades, aod::V0Datas const& V0s, daughterTracks const&)
639+
void process(AodCollisions::iterator const& collision, aod::BCsWithTimestamps const&, AodTracks const& tracks, aod::CascDataExt const& Cascades, aod::V0Datas const& V0s, DaughterTracks const&)
640640
{
641-
int Ntot = tracks.size();
642-
int CandNum_all[4] = {0, 0, 0, 0};
643-
int CandNum[4] = {0, 0, 0, 0};
641+
int nTot = tracks.size();
642+
int candNum_all[4] = {0, 0, 0, 0};
643+
int candNum[4] = {0, 0, 0, 0};
644644
for (int i = 0; i < 4; i++) {
645645
registry.fill(HIST("hEventCount"), 0.5, i + 0.5);
646646
}
647-
if (Ntot < 1)
647+
if (nTot < 1)
648648
return;
649649
fGFW->Clear();
650650
const auto cent = collision.centFT0C();
@@ -656,7 +656,7 @@ struct FlowGFWOmegaXi {
656656
loadCorrections(bc.timestamp());
657657
float vtxz = collision.posZ();
658658
registry.fill(HIST("hVtxZ"), vtxz);
659-
registry.fill(HIST("hMult"), Ntot);
659+
registry.fill(HIST("hMult"), nTot);
660660
registry.fill(HIST("hCent"), collision.centFT0C());
661661
for (int i = 0; i < 4; i++) {
662662
registry.fill(HIST("hEventCount"), 1.5, i + 0.5);
@@ -685,8 +685,8 @@ struct FlowGFWOmegaXi {
685685
}
686686
// fill GFW of V0 flow
687687
for (const auto& v0 : V0s) {
688-
auto v0posdau = v0.posTrack_as<daughterTracks>();
689-
auto v0negdau = v0.negTrack_as<daughterTracks>();
688+
auto v0posdau = v0.posTrack_as<DaughterTracks>();
689+
auto v0negdau = v0.negTrack_as<DaughterTracks>();
690690
// check tpc
691691
bool isK0s = false;
692692
bool isLambda = false;
@@ -701,7 +701,7 @@ struct FlowGFWOmegaXi {
701701
if (!setCurrentParticleWeights(weff, wacc, v0, vtxz, 1))
702702
continue;
703703
isK0s = true;
704-
CandNum_all[0] = CandNum_all[0] + 1;
704+
candNum_all[0] = candNum_all[0] + 1;
705705
registry.fill(HIST("hqaarm_podoafter"), v0.alpha(), v0.qtarm());
706706
}
707707
// Lambda and antiLambda
@@ -712,15 +712,15 @@ struct FlowGFWOmegaXi {
712712
if (!setCurrentParticleWeights(weff, wacc, v0, vtxz, 2))
713713
continue;
714714
isLambda = true;
715-
CandNum_all[1] = CandNum_all[1] + 1;
715+
candNum_all[1] = candNum_all[1] + 1;
716716
} else if (std::fabs(v0.y()) < 0.5 && std::fabs(v0.mLambda() - o2::constants::physics::MassLambda) < cfgv0_mlambdawindow &&
717717
(!cfgcheckDauTPC || (std::fabs(v0negdau.tpcNSigmaPr()) < cfgtpcNSigmaCascProton && std::fabs(v0posdau.tpcNSigmaPi()) < cfgtpcNSigmaCascPion)) &&
718718
(!cfgcheckDauTOF || ((std::fabs(v0negdau.tofNSigmaPr()) < cfgtofNSigmaCascProton || v0posdau.pt() < 0.4) && (std::fabs(v0posdau.tofNSigmaPi()) < cfgtofNSigmaCascPion || v0negdau.pt() < 0.4)))) {
719719
registry.fill(HIST("InvMassLambda_all"), v0.pt(), v0.mLambda(), v0.eta(), cent);
720720
if (!setCurrentParticleWeights(weff, wacc, v0, vtxz, 2))
721721
continue;
722722
isLambda = true;
723-
CandNum_all[1] = CandNum_all[1] + 1;
723+
candNum_all[1] = candNum_all[1] + 1;
724724
}
725725
// fill QA before cut
726726
registry.fill(HIST("hqaV0radiusbefore"), v0.v0radius());
@@ -762,15 +762,15 @@ struct FlowGFWOmegaXi {
762762
registry.fill(HIST("hqadcapostoPVafter"), v0.dcapostopv());
763763
registry.fill(HIST("hqadcanegtoPVafter"), v0.dcanegtopv());
764764
if (isK0s) {
765-
CandNum[0] = CandNum[0] + 1;
765+
candNum[0] = candNum[0] + 1;
766766
registry.fill(HIST("InvMassK0s"), v0.pt(), v0.mK0Short(), v0.eta(), cent);
767767
registry.fill(HIST("hEtaPhiVtxzPOIK0s"), v0.phi(), v0.eta(), vtxz, wacc);
768768
fGFW->Fill(v0.eta(), fV0PtAxis->FindBin(v0.pt()) - 1 + ((fK0sMass->FindBin(v0.mK0Short()) - 1) * nV0PtBins), v0.phi(), wacc * weff, 8);
769769
if (cfgOutputNUAWeights)
770770
fWeightsK0s->Fill(v0.phi(), v0.eta(), vtxz, v0.pt(), cent, 0);
771771
}
772772
if (isLambda) {
773-
CandNum[1] = CandNum[1] + 1;
773+
candNum[1] = candNum[1] + 1;
774774
registry.fill(HIST("InvMassLambda"), v0.pt(), v0.mLambda(), v0.eta(), cent);
775775
registry.fill(HIST("hEtaPhiVtxzPOILambda"), v0.phi(), v0.eta(), vtxz, wacc);
776776
fGFW->Fill(v0.eta(), fV0PtAxis->FindBin(v0.pt()) - 1 + ((fLambdaMass->FindBin(v0.mLambda()) - 1) * nV0PtBins), v0.phi(), wacc * weff, 16);
@@ -780,9 +780,9 @@ struct FlowGFWOmegaXi {
780780
}
781781
// fill GFW of casc flow
782782
for (const auto& casc : Cascades) {
783-
auto bachelor = casc.bachelor_as<daughterTracks>();
784-
auto posdau = casc.posTrack_as<daughterTracks>();
785-
auto negdau = casc.negTrack_as<daughterTracks>();
783+
auto bachelor = casc.bachelor_as<DaughterTracks>();
784+
auto posdau = casc.posTrack_as<DaughterTracks>();
785+
auto negdau = casc.negTrack_as<DaughterTracks>();
786786
// check TPC
787787
if (cfgcheckDauTPC && (!posdau.hasTPC() || !negdau.hasTPC() || !bachelor.hasTPC())) {
788788
continue;
@@ -797,15 +797,15 @@ struct FlowGFWOmegaXi {
797797
if (!setCurrentParticleWeights(weff, wacc, casc, vtxz, 4))
798798
continue;
799799
isOmega = true;
800-
CandNum_all[3] = CandNum_all[3] + 1;
800+
candNum_all[3] = candNum_all[3] + 1;
801801
} else if (casc.sign() < 0 && (casc.mOmega() > 1.63) && (casc.mOmega() < 1.71) && std::fabs(casc.yOmega()) < cfgCasc_rapidity &&
802802
(!cfgcheckDauTPC || (std::fabs(bachelor.tpcNSigmaKa()) < cfgtpcNSigmaCascKaon && std::fabs(negdau.tpcNSigmaPr()) < cfgtpcNSigmaCascProton && std::fabs(posdau.tpcNSigmaPi()) < cfgtpcNSigmaCascPion)) &&
803803
(!cfgcheckDauTOF || ((std::fabs(bachelor.tofNSigmaKa()) < cfgtofNSigmaCascKaon || bachelor.pt() < 0.4) && (std::fabs(negdau.tofNSigmaPr()) < cfgtofNSigmaCascProton || negdau.pt() < 0.4) && (std::fabs(posdau.tofNSigmaPi()) < cfgtofNSigmaCascPion || posdau.pt() < 0.4)))) {
804804
registry.fill(HIST("InvMassOmega_all"), casc.pt(), casc.mOmega(), casc.eta(), cent);
805805
if (!setCurrentParticleWeights(weff, wacc, casc, vtxz, 4))
806806
continue;
807807
isOmega = true;
808-
CandNum_all[3] = CandNum_all[3] + 1;
808+
candNum_all[3] = candNum_all[3] + 1;
809809
}
810810
// Xi and antiXi
811811
if (casc.sign() < 0 && (casc.mXi() > 1.30) && (casc.mXi() < 1.37) && std::fabs(casc.yXi()) < cfgCasc_rapidity &&
@@ -815,15 +815,15 @@ struct FlowGFWOmegaXi {
815815
if (!setCurrentParticleWeights(weff, wacc, casc, vtxz, 3))
816816
continue;
817817
isXi = true;
818-
CandNum_all[2] = CandNum_all[2] + 1;
818+
candNum_all[2] = candNum_all[2] + 1;
819819
} else if (casc.sign() < 0 && (casc.mXi() > 1.30) && (casc.mXi() < 1.37) && std::fabs(casc.yXi()) < cfgCasc_rapidity &&
820820
(!cfgcheckDauTPC || (std::fabs(bachelor.tpcNSigmaPi()) < cfgtpcNSigmaCascPion && std::fabs(posdau.tpcNSigmaPr()) < cfgtpcNSigmaCascProton && std::fabs(negdau.tpcNSigmaPi()) < cfgtpcNSigmaCascPion)) &&
821821
(!cfgcheckDauTOF || ((std::fabs(bachelor.tofNSigmaPi()) < cfgtofNSigmaCascPion || bachelor.pt() < 0.4) && (std::fabs(negdau.tofNSigmaPr()) < cfgtofNSigmaCascProton || negdau.pt() < 0.4) && (std::fabs(posdau.tofNSigmaPi()) < cfgtofNSigmaCascPion || posdau.pt() < 0.4)))) {
822822
registry.fill(HIST("InvMassXi_all"), casc.pt(), casc.mXi(), casc.eta(), cent);
823823
if (!setCurrentParticleWeights(weff, wacc, casc, vtxz, 3))
824824
continue;
825825
isXi = true;
826-
CandNum_all[2] = CandNum_all[2] + 1;
826+
candNum_all[2] = candNum_all[2] + 1;
827827
}
828828
// topological cut
829829
if (casc.cascradius() < cfgcasc_radius)
@@ -856,15 +856,15 @@ struct FlowGFWOmegaXi {
856856
if (negdau.itsNCls() < cfgitsclusters)
857857
continue;
858858
if (isOmega) {
859-
CandNum[3] = CandNum[3] + 1;
859+
candNum[3] = candNum[3] + 1;
860860
registry.fill(HIST("hEtaPhiVtxzPOIOmega"), casc.phi(), casc.eta(), vtxz, wacc);
861861
registry.fill(HIST("InvMassOmega"), casc.pt(), casc.mOmega(), casc.eta(), cent);
862862
fGFW->Fill(casc.eta(), fXiPtAxis->FindBin(casc.pt()) - 1 + ((fOmegaMass->FindBin(casc.mOmega()) - 1) * nXiPtBins), casc.phi(), wacc * weff, 4);
863863
if (cfgOutputNUAWeights)
864864
fWeightsOmega->Fill(casc.phi(), casc.eta(), vtxz, casc.pt(), cent, 0);
865865
}
866866
if (isXi) {
867-
CandNum[2] = CandNum[2] + 1;
867+
candNum[2] = candNum[2] + 1;
868868
registry.fill(HIST("hEtaPhiVtxzPOIXi"), casc.phi(), casc.eta(), vtxz, wacc);
869869
registry.fill(HIST("InvMassXi"), casc.pt(), casc.mXi(), casc.eta(), cent);
870870
fGFW->Fill(casc.eta(), fXiPtAxis->FindBin(casc.pt()) - 1 + ((fXiMass->FindBin(casc.mXi()) - 1) * nXiPtBins), casc.phi(), wacc * weff, 2);
@@ -873,10 +873,10 @@ struct FlowGFWOmegaXi {
873873
}
874874
}
875875
for (int i = 0; i < 4; i++) {
876-
if (CandNum_all[i] > 0) {
876+
if (candNum_all[i] > 0) {
877877
registry.fill(HIST("hEventCount"), 2.5, i + 0.5);
878878
}
879-
if (CandNum[i] > 0) {
879+
if (candNum[i] > 0) {
880880
registry.fill(HIST("hEventCount"), 3.5, i + 0.5);
881881
}
882882
}
@@ -910,9 +910,9 @@ struct FlowGFWOmegaXi {
910910
}
911911
// Fill subevents flow
912912
TRandom3* fRdm = new TRandom3(0);
913-
double Eventrdm = 10 * fRdm->Rndm();
913+
double eventrdm = 10 * fRdm->Rndm();
914914
for (int j = 1; j <= 10; j++) {
915-
if (Eventrdm > (j - 1) && Eventrdm < j)
915+
if (eventrdm > (j - 1) && eventrdm < j)
916916
continue;
917917
fillProfile(corrconfigs.at(15), refc22[j - 1], cent);
918918
fillProfile(corrconfigs.at(16), refc24[j - 1], cent);

0 commit comments

Comments
 (0)