Skip to content

Commit 267de0c

Browse files
fuchuncuifuchuncui
authored
PWGCF: Add dead ITS staves cut
1 parent 3ec537f commit 267de0c

File tree

1 file changed

+85
-50
lines changed

1 file changed

+85
-50
lines changed

‎PWGCF/Flow/Tasks/flowGFWOmegaXi.cxx‎

Lines changed: 85 additions & 50 deletions
Original file line numberDiff line numberDiff line change
@@ -145,8 +145,8 @@ struct FlowGFWOmegaXi {
145145
OutputObj<GFWWeights> fWeightsREF{GFWWeights("weightsREF")};
146146
OutputObj<GFWWeights> fWeightsK0s{GFWWeights("weightsK0s")};
147147
OutputObj<GFWWeights> fWeightsLambda{GFWWeights("weightsLambda")};
148-
OutputObj<GFWWeights> fWeightsXi{GFWWeights("weightsXiMinus")};
149-
OutputObj<GFWWeights> fWeightsOmega{GFWWeights("weightsOmegaMinus")};
148+
OutputObj<GFWWeights> fWeightsXi{GFWWeights("weightsXi")};
149+
OutputObj<GFWWeights> fWeightsOmega{GFWWeights("weightsOmega")};
150150

151151
// define global variables
152152
GFW* fGFW = new GFW(); // GFW class used from main src
@@ -202,6 +202,7 @@ struct FlowGFWOmegaXi {
202202

203203
// Add some output objects to the histogram registry
204204
registry.add("hPhi", "", {HistType::kTH1D, {cfgaxisPhi}});
205+
registry.add("hPhicorr", "", {HistType::kTH1D, {cfgaxisPhi}});
205206
registry.add("hEta", "", {HistType::kTH1D, {cfgaxisEta}});
206207
registry.add("hVtxZ", "", {HistType::kTH1D, {cfgaxisVertex}});
207208
registry.add("hMult", "", {HistType::kTH1D, {{3000, 0.5, 3000.5}}});
@@ -266,10 +267,10 @@ struct FlowGFWOmegaXi {
266267
Lambdac24[i - 1] = registry.add<TProfile3D>(Form("Jackknife/Lambda/Lambdac24dpt_%d", i), ";pt ; C_{2}{4} ", {HistType::kTProfile3D, {cfgaxisPtV0, cfgaxisLambdaMassforflow, axisMultiplicity}});
267268
}
268269
// InvMass(GeV) of casc and v0
269-
registry.add("InvMassXiMinus_all", "", {HistType::kTHnSparseF, {cfgaxisPtXi, axisXiminusMass, cfgaxisEta, axisMultiplicity}});
270-
registry.add("InvMassOmegaMinus_all", "", {HistType::kTHnSparseF, {cfgaxisPtXi, axisOmegaminusMass, cfgaxisEta, axisMultiplicity}});
271-
registry.add("InvMassOmegaMinus", "", {HistType::kTHnSparseF, {cfgaxisPtXi, axisOmegaminusMass, cfgaxisEta, axisMultiplicity}});
272-
registry.add("InvMassXiMinus", "", {HistType::kTHnSparseF, {cfgaxisPtXi, axisXiminusMass, cfgaxisEta, axisMultiplicity}});
270+
registry.add("InvMassXi_all", "", {HistType::kTHnSparseF, {cfgaxisPtXi, axisXiminusMass, cfgaxisEta, axisMultiplicity}});
271+
registry.add("InvMassOmega_all", "", {HistType::kTHnSparseF, {cfgaxisPtXi, axisOmegaminusMass, cfgaxisEta, axisMultiplicity}});
272+
registry.add("InvMassOmega", "", {HistType::kTHnSparseF, {cfgaxisPtXi, axisOmegaminusMass, cfgaxisEta, axisMultiplicity}});
273+
registry.add("InvMassXi", "", {HistType::kTHnSparseF, {cfgaxisPtXi, axisXiminusMass, cfgaxisEta, axisMultiplicity}});
273274
registry.add("InvMassK0s_all", "", {HistType::kTHnSparseF, {cfgaxisPtV0, axisK0sMass, cfgaxisEta, axisMultiplicity}});
274275
registry.add("InvMassLambda_all", "", {HistType::kTHnSparseF, {cfgaxisPtV0, axisLambdaMass, cfgaxisEta, axisMultiplicity}});
275276
registry.add("InvMassK0s", "", {HistType::kTHnSparseF, {cfgaxisPtV0, axisK0sMass, cfgaxisEta, axisMultiplicity}});
@@ -308,6 +309,7 @@ struct FlowGFWOmegaXi {
308309
fGFW->AddRegion("poiP10dpt", 0.4, 0.8, nPtBins, 32);
309310
fGFW->AddRegion("poifulldpt", -0.8, 0.8, nPtBins, 32);
310311
fGFW->AddRegion("poioldpt", -0.8, 0.8, nPtBins, 1);
312+
311313
int nXiptMassBins = nXiPtBins * cfgXiMassbins;
312314
fGFW->AddRegion("poiXiPdpt", 0.4, 0.8, nXiptMassBins, 2);
313315
fGFW->AddRegion("poiXiNdpt", -0.8, -0.4, nXiptMassBins, 2);
@@ -598,6 +600,10 @@ struct FlowGFWOmegaXi {
598600
// no collisions in specified time range
599601
return 0;
600602
}
603+
if (!collision.selection_bit(o2::aod::evsel::kIsGoodITSLayersAll)) {
604+
// cut time intervals with dead ITS staves
605+
return 0;
606+
}
601607
float vtxz = -999;
602608
if (collision.numContrib() > 1) {
603609
vtxz = collision.posZ();
@@ -661,6 +667,7 @@ struct FlowGFWOmegaXi {
661667
if (!setCurrentParticleWeights(weff, wacc, track, vtxz, 0))
662668
continue;
663669
registry.fill(HIST("hPhi"), track.phi());
670+
registry.fill(HIST("hPhicorr"), track.phi(), wacc);
664671
registry.fill(HIST("hEta"), track.eta());
665672
registry.fill(HIST("hEtaPhiVtxzREF"), track.phi(), track.eta(), vtxz, wacc);
666673
registry.fill(HIST("hPt"), track.pt());
@@ -679,25 +686,38 @@ struct FlowGFWOmegaXi {
679686
auto v0posdau = v0.posTrack_as<DaughterTracks>();
680687
auto v0negdau = v0.negTrack_as<DaughterTracks>();
681688
// check tpc
682-
int PDGCode = 0;
689+
bool isK0s = false;
690+
bool isLambda = false;
683691
// fill QA
684692
registry.fill(HIST("hqaarm_podobefore"), v0.alpha(), v0.qtarm());
685693
// check daughter TPC and TOF
686-
if (v0.qtarm() / std::fabs(v0.alpha()) > cfgv0_ArmPodocut &&
694+
// K0short
695+
if (v0.qtarm() / std::fabs(v0.alpha()) > cfgv0_ArmPodocut && std::fabs(v0.y()) < 0.5 && std::fabs(v0.mK0Short() - o2::constants::physics::MassK0Short) < cfgv0_mk0swindow &&
687696
(!cfgcheckDauTPC || (std::fabs(v0posdau.tpcNSigmaPi()) < cfgtpcNSigmaCascPion && std::fabs(v0negdau.tpcNSigmaPi()) < cfgtpcNSigmaCascPion)) &&
688-
(!cfgcheckDauTOF || (std::fabs(v0posdau.tofNSigmaPi()) < cfgtofNSigmaCascPion && std::fabs(v0negdau.tofNSigmaPi()) < cfgtofNSigmaCascPion))) {
697+
(!cfgcheckDauTOF || ((std::fabs(v0posdau.tofNSigmaPi()) < cfgtofNSigmaCascPion || v0posdau.pt() < 0.4) && (std::fabs(v0negdau.tofNSigmaPi()) < cfgtofNSigmaCascPion || v0negdau.pt() < 0.4)))) {
689698
registry.fill(HIST("InvMassK0s_all"), v0.pt(), v0.mK0Short(), v0.eta(), cent);
690699
if (!setCurrentParticleWeights(weff, wacc, v0, vtxz, 1))
691700
continue;
692-
PDGCode = kK0Short;
701+
isK0s = true;
693702
CandNum_all[0] = CandNum_all[0] + 1;
694703
registry.fill(HIST("hqaarm_podoafter"), v0.alpha(), v0.qtarm());
695-
} else if ((!cfgcheckDauTPC || (std::fabs(v0posdau.tpcNSigmaPr()) < cfgtpcNSigmaCascProton && std::fabs(v0negdau.tpcNSigmaPi()) < cfgtpcNSigmaCascPion)) &&
696-
(!cfgcheckDauTOF || (std::fabs(v0posdau.tofNSigmaPr()) < cfgtofNSigmaCascProton && std::fabs(v0negdau.tofNSigmaPi()) < cfgtofNSigmaCascPion))) {
704+
}
705+
// Lambda and antiLambda
706+
if (std::fabs(v0.y()) < 0.5 && std::fabs(v0.mLambda() - o2::constants::physics::MassLambda) < cfgv0_mlambdawindow &&
707+
(!cfgcheckDauTPC || (std::fabs(v0posdau.tpcNSigmaPr()) < cfgtpcNSigmaCascProton && std::fabs(v0negdau.tpcNSigmaPi()) < cfgtpcNSigmaCascPion)) &&
708+
(!cfgcheckDauTOF || ((std::fabs(v0posdau.tofNSigmaPr()) < cfgtofNSigmaCascProton || v0posdau.pt() < 0.4) && (std::fabs(v0negdau.tofNSigmaPi()) < cfgtofNSigmaCascPion || v0negdau.pt() < 0.4)))) {
709+
registry.fill(HIST("InvMassLambda_all"), v0.pt(), v0.mLambda(), v0.eta(), cent);
710+
if (!setCurrentParticleWeights(weff, wacc, v0, vtxz, 2))
711+
continue;
712+
isLambda = true;
713+
CandNum_all[1] = CandNum_all[1] + 1;
714+
} else if (std::fabs(v0.y()) < 0.5 && std::fabs(v0.mLambda() - o2::constants::physics::MassLambda) < cfgv0_mlambdawindow &&
715+
(!cfgcheckDauTPC || (std::fabs(v0negdau.tpcNSigmaPr()) < cfgtpcNSigmaCascProton && std::fabs(v0posdau.tpcNSigmaPi()) < cfgtpcNSigmaCascPion)) &&
716+
(!cfgcheckDauTOF || ((std::fabs(v0negdau.tofNSigmaPr()) < cfgtofNSigmaCascProton || v0posdau.pt() < 0.4) && (std::fabs(v0posdau.tofNSigmaPi()) < cfgtofNSigmaCascPion || v0negdau.pt() < 0.4)))) {
697717
registry.fill(HIST("InvMassLambda_all"), v0.pt(), v0.mLambda(), v0.eta(), cent);
698718
if (!setCurrentParticleWeights(weff, wacc, v0, vtxz, 2))
699719
continue;
700-
PDGCode = kLambda0;
720+
isLambda = true;
701721
CandNum_all[1] = CandNum_all[1] + 1;
702722
}
703723
// fill QA before cut
@@ -739,24 +759,21 @@ struct FlowGFWOmegaXi {
739759
registry.fill(HIST("hqadcaV0dauafter"), v0.dcaV0daughters());
740760
registry.fill(HIST("hqadcapostoPVafter"), v0.dcapostopv());
741761
registry.fill(HIST("hqadcanegtoPVafter"), v0.dcanegtopv());
742-
if (PDGCode == kK0Short) {
743-
if (std::fabs(v0.mK0Short() - o2::constants::physics::MassK0Short) < cfgv0_mk0swindow) {
744-
CandNum[0] = CandNum[0] + 1;
745-
registry.fill(HIST("InvMassK0s"), v0.pt(), v0.mK0Short(), v0.eta(), cent);
746-
registry.fill(HIST("hEtaPhiVtxzPOIK0s"), v0.phi(), v0.eta(), vtxz, wacc);
747-
fGFW->Fill(v0.eta(), fV0PtAxis->FindBin(v0.pt()) - 1 + ((fK0sMass->FindBin(v0.mK0Short()) - 1) * nV0PtBins), v0.phi(), wacc * weff, 8);
748-
if (cfgOutputNUAWeights)
749-
fWeightsK0s->Fill(v0.phi(), v0.eta(), vtxz, v0.pt(), cent, 0);
750-
}
751-
} else if (PDGCode == kLambda0) {
752-
if (std::fabs(v0.mLambda() - o2::constants::physics::MassLambda0) < cfgv0_mlambdawindow) {
753-
CandNum[1] = CandNum[1] + 1;
754-
registry.fill(HIST("InvMassLambda"), v0.pt(), v0.mLambda(), v0.eta(), cent);
755-
registry.fill(HIST("hEtaPhiVtxzPOILambda"), v0.phi(), v0.eta(), vtxz, wacc);
756-
fGFW->Fill(v0.eta(), fV0PtAxis->FindBin(v0.pt()) - 1 + ((fLambdaMass->FindBin(v0.mLambda()) - 1) * nV0PtBins), v0.phi(), wacc * weff, 16);
757-
if (cfgOutputNUAWeights)
758-
fWeightsLambda->Fill(v0.phi(), v0.eta(), vtxz, v0.pt(), cent, 0);
759-
}
762+
if (isK0s) {
763+
CandNum[0] = CandNum[0] + 1;
764+
registry.fill(HIST("InvMassK0s"), v0.pt(), v0.mK0Short(), v0.eta(), cent);
765+
registry.fill(HIST("hEtaPhiVtxzPOIK0s"), v0.phi(), v0.eta(), vtxz, wacc);
766+
fGFW->Fill(v0.eta(), fV0PtAxis->FindBin(v0.pt()) - 1 + ((fK0sMass->FindBin(v0.mK0Short()) - 1) * nV0PtBins), v0.phi(), wacc * weff, 8);
767+
if (cfgOutputNUAWeights)
768+
fWeightsK0s->Fill(v0.phi(), v0.eta(), vtxz, v0.pt(), cent, 0);
769+
}
770+
if (isLambda) {
771+
CandNum[1] = CandNum[1] + 1;
772+
registry.fill(HIST("InvMassLambda"), v0.pt(), v0.mLambda(), v0.eta(), cent);
773+
registry.fill(HIST("hEtaPhiVtxzPOILambda"), v0.phi(), v0.eta(), vtxz, wacc);
774+
fGFW->Fill(v0.eta(), fV0PtAxis->FindBin(v0.pt()) - 1 + ((fLambdaMass->FindBin(v0.mLambda()) - 1) * nV0PtBins), v0.phi(), wacc * weff, 16);
775+
if (cfgOutputNUAWeights)
776+
fWeightsLambda->Fill(v0.phi(), v0.eta(), vtxz, v0.pt(), cent, 0);
760777
}
761778
}
762779
// fill GFW of casc flow
@@ -768,21 +785,43 @@ struct FlowGFWOmegaXi {
768785
if (cfgcheckDauTPC && (!posdau.hasTPC() || !negdau.hasTPC() || !bachelor.hasTPC())) {
769786
continue;
770787
}
771-
int PDGCode = 0;
772-
if (casc.sign() < 0 && std::fabs(casc.yOmega()) < cfgCasc_rapidity && std::fabs(bachelor.tpcNSigmaKa()) < cfgtpcNSigmaCascKaon && std::fabs(posdau.tpcNSigmaPr()) < cfgtpcNSigmaCascProton && std::fabs(negdau.tpcNSigmaPi()) < cfgtpcNSigmaCascPion && std::fabs(bachelor.tofNSigmaKa()) < cfgtofNSigmaCascKaon && std::fabs(posdau.tofNSigmaPr()) < cfgtofNSigmaCascProton && std::fabs(negdau.tofNSigmaPi()) < cfgtofNSigmaCascPion) {
773-
registry.fill(HIST("InvMassOmegaMinus_all"), casc.pt(), casc.mOmega(), casc.eta(), cent);
788+
bool isOmega = false;
789+
bool isXi = false;
790+
// Omega and antiOmega
791+
if (casc.sign() < 0 && (casc.mOmega() > 1.63) && (casc.mOmega() < 1.71) && std::fabs(casc.yOmega()) < cfgCasc_rapidity &&
792+
(!cfgcheckDauTPC || (std::fabs(bachelor.tpcNSigmaKa()) < cfgtpcNSigmaCascKaon && std::fabs(posdau.tpcNSigmaPr()) < cfgtpcNSigmaCascProton && std::fabs(negdau.tpcNSigmaPi()) < cfgtpcNSigmaCascPion)) &&
793+
(!cfgcheckDauTOF || ((std::fabs(bachelor.tofNSigmaKa()) < cfgtofNSigmaCascKaon || bachelor.pt() < 0.4) && (std::fabs(posdau.tofNSigmaPr()) < cfgtofNSigmaCascProton || posdau.pt() < 0.4) && (std::fabs(negdau.tofNSigmaPi()) < cfgtofNSigmaCascPion || negdau.pt() < 0.4)))) {
794+
registry.fill(HIST("InvMassOmega_all"), casc.pt(), casc.mOmega(), casc.eta(), cent);
774795
if (!setCurrentParticleWeights(weff, wacc, casc, vtxz, 4))
775796
continue;
776-
PDGCode = kOmegaMinus;
797+
isOmega = true;
777798
CandNum_all[3] = CandNum_all[3] + 1;
778-
} else if (casc.sign() < 0 && std::fabs(casc.yXi()) < cfgCasc_rapidity && std::fabs(bachelor.tpcNSigmaPi()) < cfgtpcNSigmaCascPion && std::fabs(posdau.tpcNSigmaPr()) < cfgtpcNSigmaCascProton && std::fabs(negdau.tpcNSigmaPi()) < cfgtpcNSigmaCascPion && std::fabs(bachelor.tofNSigmaPi()) < cfgtofNSigmaCascPion && std::fabs(posdau.tofNSigmaPr()) < cfgtofNSigmaCascProton && std::fabs(negdau.tofNSigmaPi()) < cfgtofNSigmaCascPion) {
779-
registry.fill(HIST("InvMassXiMinus_all"), casc.pt(), casc.mXi(), casc.eta(), cent);
799+
} else if (casc.sign() < 0 && (casc.mOmega() > 1.63) && (casc.mOmega() < 1.71) && std::fabs(casc.yOmega()) < cfgCasc_rapidity &&
800+
(!cfgcheckDauTPC || (std::fabs(bachelor.tpcNSigmaKa()) < cfgtpcNSigmaCascKaon && std::fabs(negdau.tpcNSigmaPr()) < cfgtpcNSigmaCascProton && std::fabs(posdau.tpcNSigmaPi()) < cfgtpcNSigmaCascPion)) &&
801+
(!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)))) {
802+
registry.fill(HIST("InvMassOmega_all"), casc.pt(), casc.mOmega(), casc.eta(), cent);
803+
if (!setCurrentParticleWeights(weff, wacc, casc, vtxz, 4))
804+
continue;
805+
isOmega = true;
806+
CandNum_all[3] = CandNum_all[3] + 1;
807+
}
808+
// Xi and antiXi
809+
if (casc.sign() < 0 && (casc.mXi() > 1.30) && (casc.mXi() < 1.37) && std::fabs(casc.yXi()) < cfgCasc_rapidity &&
810+
(!cfgcheckDauTPC || (std::fabs(bachelor.tpcNSigmaPi()) < cfgtpcNSigmaCascPion && std::fabs(posdau.tpcNSigmaPr()) < cfgtpcNSigmaCascProton && std::fabs(negdau.tpcNSigmaPi()) < cfgtpcNSigmaCascPion)) &&
811+
(!cfgcheckDauTOF || ((std::fabs(bachelor.tofNSigmaPi()) < cfgtofNSigmaCascPion || bachelor.pt() < 0.4) && (std::fabs(posdau.tofNSigmaPr()) < cfgtofNSigmaCascProton || posdau.pt() < 0.4) && (std::fabs(negdau.tofNSigmaPi()) < cfgtofNSigmaCascPion || negdau.pt() < 0.4)))) {
812+
registry.fill(HIST("InvMassXi_all"), casc.pt(), casc.mXi(), casc.eta(), cent);
780813
if (!setCurrentParticleWeights(weff, wacc, casc, vtxz, 3))
781814
continue;
782-
PDGCode = kXiMinus;
815+
isXi = true;
816+
CandNum_all[2] = CandNum_all[2] + 1;
817+
} else if (casc.sign() < 0 && (casc.mXi() > 1.30) && (casc.mXi() < 1.37) && std::fabs(casc.yXi()) < cfgCasc_rapidity &&
818+
(!cfgcheckDauTPC || (std::fabs(bachelor.tpcNSigmaPi()) < cfgtpcNSigmaCascPion && std::fabs(posdau.tpcNSigmaPr()) < cfgtpcNSigmaCascProton && std::fabs(negdau.tpcNSigmaPi()) < cfgtpcNSigmaCascPion)) &&
819+
(!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)))) {
820+
registry.fill(HIST("InvMassXi_all"), casc.pt(), casc.mXi(), casc.eta(), cent);
821+
if (!setCurrentParticleWeights(weff, wacc, casc, vtxz, 3))
822+
continue;
823+
isXi = true;
783824
CandNum_all[2] = CandNum_all[2] + 1;
784-
} else {
785-
continue;
786825
}
787826
// topological cut
788827
if (casc.cascradius() < cfgcasc_radius)
@@ -814,23 +853,19 @@ struct FlowGFWOmegaXi {
814853
continue;
815854
if (negdau.itsNCls() < cfgitsclusters)
816855
continue;
817-
if (PDGCode == kOmegaMinus) {
856+
if (isOmega) {
818857
CandNum[3] = CandNum[3] + 1;
819858
registry.fill(HIST("hEtaPhiVtxzPOIOmega"), casc.phi(), casc.eta(), vtxz, wacc);
820-
registry.fill(HIST("InvMassOmegaMinus"), casc.pt(), casc.mOmega(), casc.eta(), cent);
821-
if ((casc.pt() < cfgCutPtPOIMax) && (casc.pt() > cfgCutPtPOIMin) && (casc.mOmega() > 1.63) && (casc.mOmega() < 1.71)) {
822-
fGFW->Fill(casc.eta(), fXiPtAxis->FindBin(casc.pt()) - 1 + ((fOmegaMass->FindBin(casc.mOmega()) - 1) * nXiPtBins), casc.phi(), wacc * weff, 4);
823-
}
859+
registry.fill(HIST("InvMassOmega"), casc.pt(), casc.mOmega(), casc.eta(), cent);
860+
fGFW->Fill(casc.eta(), fXiPtAxis->FindBin(casc.pt()) - 1 + ((fOmegaMass->FindBin(casc.mOmega()) - 1) * nXiPtBins), casc.phi(), wacc * weff, 4);
824861
if (cfgOutputNUAWeights)
825862
fWeightsOmega->Fill(casc.phi(), casc.eta(), vtxz, casc.pt(), cent, 0);
826863
}
827-
if (PDGCode == kXiMinus) {
864+
if (isXi) {
828865
CandNum[2] = CandNum[2] + 1;
829866
registry.fill(HIST("hEtaPhiVtxzPOIXi"), casc.phi(), casc.eta(), vtxz, wacc);
830-
registry.fill(HIST("InvMassXiMinus"), casc.pt(), casc.mXi(), casc.eta(), cent);
831-
if ((casc.pt() < cfgCutPtPOIMax) && (casc.pt() > cfgCutPtPOIMin) && (casc.mXi() > 1.30) && (casc.mXi() < 1.37)) {
832-
fGFW->Fill(casc.eta(), fXiPtAxis->FindBin(casc.pt()) - 1 + ((fXiMass->FindBin(casc.mXi()) - 1) * nXiPtBins), casc.phi(), wacc * weff, 2);
833-
}
867+
registry.fill(HIST("InvMassXi"), casc.pt(), casc.mXi(), casc.eta(), cent);
868+
fGFW->Fill(casc.eta(), fXiPtAxis->FindBin(casc.pt()) - 1 + ((fXiMass->FindBin(casc.mXi()) - 1) * nXiPtBins), casc.phi(), wacc * weff, 2);
834869
if (cfgOutputNUAWeights)
835870
fWeightsXi->Fill(casc.phi(), casc.eta(), vtxz, casc.pt(), cent, 0);
836871
}

0 commit comments

Comments
 (0)