Skip to content

Commit ebe05bb

Browse files
SuJeong-Jialibuild
SuJeong-Ji
and
alibuild
authored
[PWGLF] Hot fix and clean up secondary cuts in chk892flow.cxx (#9535)
Co-authored-by: ALICE Action Bot <alibuild@cern.ch>
1 parent c6bb3ed commit ebe05bb

File tree

1 file changed

+49
-41
lines changed

1 file changed

+49
-41
lines changed

‎PWGLF/Tasks/Resonances/chk892flow.cxx‎

Lines changed: 49 additions & 41 deletions
Original file line numberDiff line numberDiff line change
@@ -186,16 +186,16 @@ struct chk892flow {
186186
Configurable<bool> cSecondaryArmenterosCut{"cSecondaryArmenterosCut", true, "cut on Armenteros-Podolanski graph"};
187187

188188
Configurable<bool> cfgByPassDauPIDSelection{"cfgByPassDauPIDSelection", true, "Bypass Daughters PID selection"};
189-
Configurable<float> cSecondaryDauDCAMax{"cSecondaryDauDCAMax", 1., "Maximum DCA Secondary daughters to PV"};
190-
Configurable<float> cSecondaryDauPosDCAtoPVMin{"cSecondaryDauPosDCAtoPVMin", 0.0, "Minimum DCA Secondary positive daughters to PV"};
191-
Configurable<float> cSecondaryDauNegDCAtoPVMin{"cSecondaryDauNegDCAtoPVMin", 0.0, "Minimum DCA Secondary negative daughters to PV"};
189+
Configurable<float> cSecondaryDauDCAMax{"cSecondaryDauDCAMax", 0.2, "Maximum DCA Secondary daughters to PV"};
190+
// Configurable<float> cSecondaryDauPosDCAtoPVMin{"cSecondaryDauPosDCAtoPVMin", 0.0, "Minimum DCA Secondary positive daughters to PV"};
191+
// Configurable<float> cSecondaryDauNegDCAtoPVMin{"cSecondaryDauNegDCAtoPVMin", 0.0, "Minimum DCA Secondary negative daughters to PV"};
192192

193-
Configurable<float> cSecondaryPtMin{"cSecondaryPtMin", 0.f, "Minimum transverse momentum of Secondary"};
193+
// Configurable<float> cSecondaryPtMin{"cSecondaryPtMin", 0.f, "Minimum transverse momentum of Secondary"};
194194
Configurable<float> cSecondaryRapidityMax{"cSecondaryRapidityMax", 0.5, "Maximum rapidity of Secondary"};
195-
Configurable<float> cSecondaryRadiusMin{"cSecondaryRadiusMin", 1.2, "Minimum transverse radius of Secondary"};
196-
Configurable<float> cSecondaryCosPAMin{"cSecondaryCosPAMin", 0.995, "Mininum cosine pointing angle of Secondary"};
197-
Configurable<float> cSecondaryDCAtoPVMax{"cSecondaryDCAtoPVMax", 0.3, "Maximum DCA Secondary to PV"};
198-
Configurable<float> cSecondaryProperLifetimeMax{"cSecondaryProperLifetimeMax", 20, "Maximum Secondary Lifetime"};
195+
// Configurable<float> cSecondaryRadiusMin{"cSecondaryRadiusMin", 1.2, "Minimum transverse radius of Secondary"};
196+
Configurable<float> cSecondaryCosPAMin{"cSecondaryCosPAMin", 0.998, "Mininum cosine pointing angle of Secondary"};
197+
Configurable<float> cSecondaryDCAtoPVMax{"cSecondaryDCAtoPVMax", 0.4, "Maximum DCA Secondary to PV"};
198+
Configurable<float> cSecondaryProperLifetimeMax{"cSecondaryProperLifetimeMax", 20., "Maximum Secondary Lifetime"};
199199
Configurable<float> cSecondaryparamArmenterosCut{"cSecondaryparamArmenterosCut", 0.2, "parameter for Armenteros Cut"};
200200
Configurable<float> cSecondaryMassWindow{"cSecondaryMassWindow", 0.03, "Secondary inv mass selciton window"};
201201

@@ -339,27 +339,27 @@ struct chk892flow {
339339

340340
// K0s
341341
histos.add("QA/before/hDauDCASecondary", "DCA of daughters of secondary resonance", HistType::kTH1D, {dcaAxis});
342-
histos.add("QA/before/hDauPosDCAtoPVSecondary", "Pos DCA to PV of daughters secondary resonance", HistType::kTH1D, {dcaAxis});
343-
histos.add("QA/before/hDauNegDCAtoPVSecondary", "Neg DCA to PV of daughters secondary resonance", HistType::kTH1D, {dcaAxis});
344-
histos.add("QA/before/hpT_Secondary", "pT distribution of secondary resonance", HistType::kTH1D, {ptAxisQA});
342+
// histos.add("QA/before/hDauPosDCAtoPVSecondary", "Pos DCA to PV of daughters secondary resonance", HistType::kTH1D, {dcaAxis});
343+
// histos.add("QA/before/hDauNegDCAtoPVSecondary", "Neg DCA to PV of daughters secondary resonance", HistType::kTH1D, {dcaAxis});
344+
// histos.add("QA/before/hpT_Secondary", "pT distribution of secondary resonance", HistType::kTH1D, {ptAxisQA});
345345
histos.add("QA/before/hy_Secondary", "Rapidity distribution of secondary resonance", HistType::kTH1D, {yAxis});
346-
histos.add("QA/before/hRadiusSecondary", "Radius distribution of secondary resonance", HistType::kTH1D, {radiusAxis});
346+
// histos.add("QA/before/hRadiusSecondary", "Radius distribution of secondary resonance", HistType::kTH1D, {radiusAxis});
347347
histos.add("QA/before/hCPASecondary", "Cosine pointing angle distribution of secondary resonance", HistType::kTH1D, {cpaAxis});
348348
histos.add("QA/before/hDCAtoPVSecondary", "DCA to PV distribution of secondary resonance", HistType::kTH1D, {dcaAxis});
349-
histos.add("QA/before/hPropTauSecondary", "Proper Lifetime distribution of secondary resonance", HistType::kTH1D, {tauAxis});
350-
histos.add("QA/before/hPtAsymSecondary", "pT asymmetry distribution of secondary resonance", HistType::kTH1D, {AxisSpec{100, -1, 1, "Pair asymmetry"}});
349+
// histos.add("QA/before/hPropTauSecondary", "Proper Lifetime distribution of secondary resonance", HistType::kTH1D, {tauAxis});
350+
// histos.add("QA/before/hPtAsymSecondary", "pT asymmetry distribution of secondary resonance", HistType::kTH1D, {AxisSpec{100, -1, 1, "Pair asymmetry"}});
351351
histos.add("QA/before/hInvmassSecondary", "Invariant mass of unlike-sign secondary resonance", HistType::kTH1D, {invMassAxisK0s});
352352

353353
histos.add("QA/after/hDauDCASecondary", "DCA of daughters of secondary resonance", HistType::kTH1D, {dcaAxis});
354-
histos.add("QA/after/hDauPosDCAtoPVSecondary", "Pos DCA to PV of daughters secondary resonance", HistType::kTH1D, {dcaAxis});
355-
histos.add("QA/after/hDauNegDCAtoPVSecondary", "Neg DCA to PV of daughters secondary resonance", HistType::kTH1D, {dcaAxis});
356-
histos.add("QA/after/hpT_Secondary", "pT distribution of secondary resonance", HistType::kTH1D, {ptAxisQA});
354+
// histos.add("QA/after/hDauPosDCAtoPVSecondary", "Pos DCA to PV of daughters secondary resonance", HistType::kTH1D, {dcaAxis});
355+
// histos.add("QA/after/hDauNegDCAtoPVSecondary", "Neg DCA to PV of daughters secondary resonance", HistType::kTH1D, {dcaAxis});
356+
// histos.add("QA/after/hpT_Secondary", "pT distribution of secondary resonance", HistType::kTH1D, {ptAxisQA});
357357
histos.add("QA/after/hy_Secondary", "Rapidity distribution of secondary resonance", HistType::kTH1D, {yAxis});
358-
histos.add("QA/after/hRadiusSecondary", "Radius distribution of secondary resonance", HistType::kTH1D, {radiusAxis});
358+
// histos.add("QA/after/hRadiusSecondary", "Radius distribution of secondary resonance", HistType::kTH1D, {radiusAxis});
359359
histos.add("QA/after/hCPASecondary", "Cosine pointing angle distribution of secondary resonance", HistType::kTH1D, {cpaAxis});
360360
histos.add("QA/after/hDCAtoPVSecondary", "DCA to PV distribution of secondary resonance", HistType::kTH1D, {dcaAxis});
361-
histos.add("QA/after/hPropTauSecondary", "Proper Lifetime distribution of secondary resonance", HistType::kTH1D, {tauAxis});
362-
histos.add("QA/after/hPtAsymSecondary", "pT asymmetry distribution of secondary resonance", HistType::kTH1D, {AxisSpec{100, -1, 1, "Pair asymmetry"}});
361+
// histos.add("QA/after/hPropTauSecondary", "Proper Lifetime distribution of secondary resonance", HistType::kTH1D, {tauAxis});
362+
// histos.add("QA/after/hPtAsymSecondary", "pT asymmetry distribution of secondary resonance", HistType::kTH1D, {AxisSpec{100, -1, 1, "Pair asymmetry"}});
363363
histos.add("QA/after/hInvmassSecondary", "Invariant mass of unlike-sign secondary resonance", HistType::kTH1D, {invMassAxisK0s});
364364

365365
// Kstar
@@ -585,11 +585,11 @@ struct chk892flow {
585585
bool selectionK0s(CollisionType const& collision, K0sType const& candidate)
586586
{
587587
auto lDauDCA = candidate.dcaV0daughters();
588-
auto lDauPosDCAtoPV = candidate.dcapostopv();
589-
auto lDauNegDCAtoPV = candidate.dcanegtopv();
590-
auto lPt = candidate.pt();
588+
// auto lDauPosDCAtoPV = candidate.dcapostopv();
589+
// auto lDauNegDCAtoPV = candidate.dcanegtopv();
590+
// auto lPt = candidate.pt();
591591
auto lRapidity = candidate.yK0Short();
592-
auto lRadius = candidate.v0radius();
592+
// auto lRadius = candidate.v0radius();
593593
auto lDCAtoPV = candidate.dcav0topv();
594594
auto lCPA = candidate.v0cosPA();
595595
auto lPropTauK0s = candidate.distovertotmom(collision.posX(), collision.posY(), collision.posZ()) * MassK0Short;
@@ -604,6 +604,7 @@ struct chk892flow {
604604
histos.fill(HIST("QA/K0sCutCheck"), 1);
605605
returnFlag = false;
606606
}
607+
/*
607608
if (lDauPosDCAtoPV < cSecondaryDauPosDCAtoPVMin) {
608609
histos.fill(HIST("QA/K0sCutCheck"), 2);
609610
returnFlag = false;
@@ -616,14 +617,17 @@ struct chk892flow {
616617
histos.fill(HIST("QA/K0sCutCheck"), 4);
617618
returnFlag = false;
618619
}
619-
if (lRapidity > cSecondaryRapidityMax) {
620+
*/
621+
if (std::fabs(lRapidity) > cSecondaryRapidityMax) {
620622
histos.fill(HIST("QA/K0sCutCheck"), 5);
621623
returnFlag = false;
622624
}
625+
/*
623626
if (lRadius < cSecondaryRadiusMin) {
624627
histos.fill(HIST("QA/K0sCutCheck"), 6);
625628
returnFlag = false;
626629
}
630+
*/
627631
if (lDCAtoPV > cSecondaryDCAtoPVMax) {
628632
histos.fill(HIST("QA/K0sCutCheck"), 7);
629633
returnFlag = false;
@@ -664,6 +668,7 @@ struct chk892flow {
664668
histos.fill(HIST("QA/K0sCutCheck"), 1);
665669
return false;
666670
}
671+
/*
667672
if (lDauPosDCAtoPV < cSecondaryDauPosDCAtoPVMin) {
668673
histos.fill(HIST("QA/K0sCutCheck"), 2);
669674
return false;
@@ -676,14 +681,17 @@ struct chk892flow {
676681
histos.fill(HIST("QA/K0sCutCheck"), 4);
677682
return false;
678683
}
679-
if (lRapidity > cSecondaryRapidityMax) {
684+
*/
685+
if (std::fabs(lRapidity) > cSecondaryRapidityMax) {
680686
histos.fill(HIST("QA/K0sCutCheck"), 5);
681687
return false;
682688
}
689+
/*
683690
if (lRadius < cSecondaryRadiusMin) {
684691
histos.fill(HIST("QA/K0sCutCheck"), 6);
685692
return false;
686693
}
694+
*/
687695
if (lDCAtoPV > cSecondaryDCAtoPVMax) {
688696
histos.fill(HIST("QA/K0sCutCheck"), 7);
689697
return false;
@@ -835,14 +843,14 @@ struct chk892flow {
835843

836844
/// K0s
837845
auto trkkDauDCA = k0sCand.dcaV0daughters();
838-
auto trkkDauDCAPostoPV = k0sCand.dcapostopv();
839-
auto trkkDauDCANegtoPV = k0sCand.dcanegtopv();
840-
auto trkkpt = k0sCand.pt();
846+
// auto trkkDauDCAPostoPV = k0sCand.dcapostopv();
847+
// auto trkkDauDCANegtoPV = k0sCand.dcanegtopv();
848+
// auto trkkpt = k0sCand.pt();
841849
auto trkky = k0sCand.yK0Short();
842-
auto trkkRadius = k0sCand.v0radius();
850+
// auto trkkRadius = k0sCand.v0radius();
843851
auto trkkDCAtoPV = k0sCand.dcav0topv();
844852
auto trkkCPA = k0sCand.v0cosPA();
845-
auto trkkPropTau = k0sCand.distovertotmom(collision.posX(), collision.posY(), collision.posZ()) * MassK0Short;
853+
// auto trkkPropTau = k0sCand.distovertotmom(collision.posX(), collision.posY(), collision.posZ()) * MassK0Short;
846854
auto trkkMass = k0sCand.mK0Short();
847855

848856
if constexpr (!IsMix) {
@@ -866,15 +874,15 @@ struct chk892flow {
866874
histos.fill(HIST("QA/before/trknpionDCAz"), negDauTrack.dcaZ());
867875

868876
histos.fill(HIST("QA/before/hDauDCASecondary"), trkkDauDCA);
869-
histos.fill(HIST("QA/before/hDauPosDCAtoPVSecondary"), trkkDauDCAPostoPV);
870-
histos.fill(HIST("QA/before/hDauNegDCAtoPVSecondary"), trkkDauDCANegtoPV);
877+
// histos.fill(HIST("QA/before/hDauPosDCAtoPVSecondary"), trkkDauDCAPostoPV);
878+
// histos.fill(HIST("QA/before/hDauNegDCAtoPVSecondary"), trkkDauDCANegtoPV);
871879

872-
histos.fill(HIST("QA/before/hpT_Secondary"), trkkpt);
880+
// histos.fill(HIST("QA/before/hpT_Secondary"), trkkpt);
873881
histos.fill(HIST("QA/before/hy_Secondary"), trkky);
874-
histos.fill(HIST("QA/before/hRadiusSecondary"), trkkRadius);
882+
// histos.fill(HIST("QA/before/hRadiusSecondary"), trkkRadius);
875883
histos.fill(HIST("QA/before/hDCAtoPVSecondary"), trkkDCAtoPV);
876884
histos.fill(HIST("QA/before/hCPASecondary"), trkkCPA);
877-
histos.fill(HIST("QA/before/hPropTauSecondary"), trkkPropTau);
885+
// histos.fill(HIST("QA/before/hPropTauSecondary"), trkkPropTau);
878886
histos.fill(HIST("QA/before/hInvmassSecondary"), trkkMass);
879887
}
880888

@@ -909,15 +917,15 @@ struct chk892flow {
909917
histos.fill(HIST("QA/after/trknpionDCAz"), negDauTrack.dcaZ());
910918

911919
histos.fill(HIST("QA/after/hDauDCASecondary"), trkkDauDCA);
912-
histos.fill(HIST("QA/after/hDauPosDCAtoPVSecondary"), trkkDauDCAPostoPV);
913-
histos.fill(HIST("QA/after/hDauNegDCAtoPVSecondary"), trkkDauDCANegtoPV);
920+
// histos.fill(HIST("QA/after/hDauPosDCAtoPVSecondary"), trkkDauDCAPostoPV);
921+
// histos.fill(HIST("QA/after/hDauNegDCAtoPVSecondary"), trkkDauDCANegtoPV);
914922

915-
histos.fill(HIST("QA/after/hpT_Secondary"), trkkpt);
923+
// histos.fill(HIST("QA/after/hpT_Secondary"), trkkpt);
916924
histos.fill(HIST("QA/after/hy_Secondary"), trkky);
917-
histos.fill(HIST("QA/after/hRadiusSecondary"), trkkRadius);
925+
// histos.fill(HIST("QA/after/hRadiusSecondary"), trkkRadius);
918926
histos.fill(HIST("QA/after/hDCAtoPVSecondary"), trkkDCAtoPV);
919927
histos.fill(HIST("QA/after/hCPASecondary"), trkkCPA);
920-
histos.fill(HIST("QA/after/hPropTauSecondary"), trkkPropTau);
928+
// histos.fill(HIST("QA/after/hPropTauSecondary"), trkkPropTau);
921929
histos.fill(HIST("QA/after/hInvmassSecondary"), trkkMass);
922930
}
923931
k0sIndicies.push_back(k0sCand.index());

0 commit comments

Comments
 (0)