Skip to content

Commit b3c0355

Browse files
fuchuncuifuchuncui
authored
[PWGCF] add cut QA and flow uncertainties
1 parent cc10477 commit b3c0355

File tree

1 file changed

+19
-19
lines changed

1 file changed

+19
-19
lines changed

‎PWGCF/Flow/Tasks/flowGFWOmegaXi.cxx‎

Lines changed: 19 additions & 19 deletions
Original file line numberDiff line numberDiff line change
@@ -395,7 +395,7 @@ struct FlowGFWOmegaXi {
395395
return;
396396
if (!corrconf.pTDif) {
397397
val = fGFW->Calculate(corrconf, 0, kFALSE).real() / dnx;
398-
if (abs(val) < 1)
398+
if (std::fabs(val) < 1)
399399
registry.fill(tarName, cent, val, dnx);
400400
return;
401401
}
@@ -411,7 +411,7 @@ struct FlowGFWOmegaXi {
411411
return;
412412
if (!corrconf.pTDif) {
413413
val = fGFW->Calculate(corrconf, 0, kFALSE).real() / dnx;
414-
if (abs(val) < 1)
414+
if (std::fabs(val) < 1)
415415
TProfile->Fill(cent, val, dnx);
416416
return;
417417
}
@@ -427,7 +427,7 @@ struct FlowGFWOmegaXi {
427427
if (dnx == 0)
428428
return;
429429
val = fGFW->Calculate(corrconf, ptbin - 1, kFALSE).real() / dnx;
430-
if (abs(val) < 1) {
430+
if (std::fabs(val) < 1) {
431431
registry.fill(tarName, fPtAxis->GetBinCenter(ptbin), cent, val, dnx);
432432
}
433433
return;
@@ -472,7 +472,7 @@ struct FlowGFWOmegaXi {
472472
if (dnx == 0)
473473
continue;
474474
val = fGFW->Calculate(corrconf, (ptbin - 1) + ((massbin - 1) * nptbins), kFALSE).real() / dnx;
475-
if (abs(val) < 1) {
475+
if (std::fabs(val) < 1) {
476476
registry.fill(tarName, fpt->GetBinCenter(ptbin), fMass->GetBinCenter(massbin), cent, val, dnx);
477477
}
478478
}
@@ -517,7 +517,7 @@ struct FlowGFWOmegaXi {
517517
if (dnx == 0)
518518
continue;
519519
val = fGFW->Calculate(corrconf, (ptbin - 1) + ((massbin - 1) * nptbins), kFALSE).real() / dnx;
520-
if (abs(val) < 1) {
520+
if (std::fabs(val) < 1) {
521521
TProfile3D->Fill(fpt->GetBinCenter(ptbin), fMass->GetBinCenter(massbin), cent, val, dnx);
522522
}
523523
}
@@ -608,7 +608,7 @@ struct FlowGFWOmegaXi {
608608
auto multNTracksPV = collision.multNTracksPV();
609609
auto occupancy = collision.trackOccupancyInTimeRange();
610610

611-
if (abs(vtxz) > cfgCutVertex)
611+
if (std::fabs(vtxz) > cfgCutVertex)
612612
return false;
613613
if (multNTracksPV < fMultPVCutLow->Eval(centrality))
614614
return false;
@@ -622,7 +622,7 @@ struct FlowGFWOmegaXi {
622622
return 0;
623623

624624
// V0A T0A 5 sigma cut
625-
if (abs(collision.multFV0A() - fT0AV0AMean->Eval(collision.multFT0A())) > 5 * fT0AV0ASigma->Eval(collision.multFT0A()))
625+
if (std::fabs(collision.multFV0A() - fT0AV0AMean->Eval(collision.multFT0A())) > 5 * fT0AV0ASigma->Eval(collision.multFT0A()))
626626
return 0;
627627

628628
return true;
@@ -683,17 +683,17 @@ struct FlowGFWOmegaXi {
683683
// fill QA
684684
registry.fill(HIST("hqaarm_podobefore"), v0.alpha(), v0.qtarm());
685685
// check daughter TPC and TOF
686-
if (v0.qtarm() / abs(v0.alpha()) > cfgv0_ArmPodocut &&
687-
(!cfgcheckDauTPC || (abs(v0posdau.tpcNSigmaPi()) < cfgtpcNSigmaCascPion && abs(v0negdau.tpcNSigmaPi()) < cfgtpcNSigmaCascPion)) &&
688-
(!cfgcheckDauTOF || (abs(v0posdau.tofNSigmaPi()) < cfgtofNSigmaCascPion && abs(v0negdau.tofNSigmaPi()) < cfgtofNSigmaCascPion))) {
686+
if (v0.qtarm() / std::fabs(v0.alpha()) > cfgv0_ArmPodocut &&
687+
(!cfgcheckDauTPC || (std::fabs(v0posdau.tpcNSigmaPi()) < cfgtpcNSigmaCascPion && std::fabs(v0negdau.tpcNSigmaPi()) < cfgtpcNSigmaCascPion)) &&
688+
(!cfgcheckDauTOF || (std::fabs(v0posdau.tofNSigmaPi()) < cfgtofNSigmaCascPion && std::fabs(v0negdau.tofNSigmaPi()) < cfgtofNSigmaCascPion))) {
689689
registry.fill(HIST("InvMassK0s_all"), v0.pt(), v0.mK0Short(), v0.eta(), cent);
690690
if (!setCurrentParticleWeights(weff, wacc, v0, vtxz, 1))
691691
continue;
692692
PDGCode = kK0Short;
693693
CandNum_all[0] = CandNum_all[0] + 1;
694694
registry.fill(HIST("hqaarm_podoafter"), v0.alpha(), v0.qtarm());
695-
} else if ((!cfgcheckDauTPC || (abs(v0posdau.tpcNSigmaPr()) < cfgtpcNSigmaCascProton && abs(v0negdau.tpcNSigmaPi()) < cfgtpcNSigmaCascPion)) &&
696-
(!cfgcheckDauTOF || (abs(v0posdau.tofNSigmaPr()) < cfgtofNSigmaCascProton && abs(v0negdau.tofNSigmaPi()) < cfgtofNSigmaCascPion))) {
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))) {
697697
registry.fill(HIST("InvMassLambda_all"), v0.pt(), v0.mLambda(), v0.eta(), cent);
698698
if (!setCurrentParticleWeights(weff, wacc, v0, vtxz, 2))
699699
continue;
@@ -729,9 +729,9 @@ struct FlowGFWOmegaXi {
729729
continue;
730730
if (v0.dcaV0daughters() > cfgv0_dcav0dau)
731731
continue;
732-
if (abs(v0.dcapostopv()) < cfgv0_dcadautopv)
732+
if (std::fabs(v0.dcapostopv()) < cfgv0_dcadautopv)
733733
continue;
734-
if (abs(v0.dcanegtopv()) < cfgv0_dcadautopv)
734+
if (std::fabs(v0.dcanegtopv()) < cfgv0_dcadautopv)
735735
continue;
736736
// fill QA after cut
737737
registry.fill(HIST("hqaV0radiusafter"), v0.v0radius());
@@ -740,7 +740,7 @@ struct FlowGFWOmegaXi {
740740
registry.fill(HIST("hqadcapostoPVafter"), v0.dcapostopv());
741741
registry.fill(HIST("hqadcanegtoPVafter"), v0.dcanegtopv());
742742
if (PDGCode == kK0Short) {
743-
if (abs(v0.mK0Short() - o2::constants::physics::MassK0Short) < cfgv0_mk0swindow) {
743+
if (std::fabs(v0.mK0Short() - o2::constants::physics::MassK0Short) < cfgv0_mk0swindow) {
744744
CandNum[0] = CandNum[0] + 1;
745745
registry.fill(HIST("InvMassK0s"), v0.pt(), v0.mK0Short(), v0.eta(), cent);
746746
registry.fill(HIST("hEtaPhiVtxzPOIK0s"), v0.phi(), v0.eta(), vtxz, wacc);
@@ -749,7 +749,7 @@ struct FlowGFWOmegaXi {
749749
fWeightsK0s->Fill(v0.phi(), v0.eta(), vtxz, v0.pt(), cent, 0);
750750
}
751751
} else if (PDGCode == kLambda0) {
752-
if (abs(v0.mLambda() - o2::constants::physics::MassLambda0) < cfgv0_mlambdawindow) {
752+
if (std::fabs(v0.mLambda() - o2::constants::physics::MassLambda0) < cfgv0_mlambdawindow) {
753753
CandNum[1] = CandNum[1] + 1;
754754
registry.fill(HIST("InvMassLambda"), v0.pt(), v0.mLambda(), v0.eta(), cent);
755755
registry.fill(HIST("hEtaPhiVtxzPOILambda"), v0.phi(), v0.eta(), vtxz, wacc);
@@ -769,13 +769,13 @@ struct FlowGFWOmegaXi {
769769
continue;
770770
}
771771
int PDGCode = 0;
772-
if (casc.sign() < 0 && abs(casc.yOmega()) < cfgCasc_rapidity && abs(bachelor.tpcNSigmaKa()) < cfgtpcNSigmaCascKaon && abs(posdau.tpcNSigmaPr()) < cfgtpcNSigmaCascProton && abs(negdau.tpcNSigmaPi()) < cfgtpcNSigmaCascPion && abs(bachelor.tofNSigmaKa()) < cfgtofNSigmaCascKaon && abs(posdau.tofNSigmaPr()) < cfgtofNSigmaCascProton && abs(negdau.tofNSigmaPi()) < cfgtofNSigmaCascPion) {
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) {
773773
registry.fill(HIST("InvMassOmegaMinus_all"), casc.pt(), casc.mOmega(), casc.eta(), cent);
774774
if (!setCurrentParticleWeights(weff, wacc, casc, vtxz, 4))
775775
continue;
776776
PDGCode = kOmegaMinus;
777777
CandNum_all[3] = CandNum_all[3] + 1;
778-
} else if (casc.sign() < 0 && abs(casc.yXi()) < cfgCasc_rapidity && abs(bachelor.tpcNSigmaPi()) < cfgtpcNSigmaCascPion && abs(posdau.tpcNSigmaPr()) < cfgtpcNSigmaCascProton && abs(negdau.tpcNSigmaPi()) < cfgtpcNSigmaCascPion && abs(bachelor.tofNSigmaPi()) < cfgtofNSigmaCascPion && abs(posdau.tofNSigmaPr()) < cfgtofNSigmaCascProton && abs(negdau.tofNSigmaPi()) < cfgtofNSigmaCascPion) {
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) {
779779
registry.fill(HIST("InvMassXiMinus_all"), casc.pt(), casc.mXi(), casc.eta(), cent);
780780
if (!setCurrentParticleWeights(weff, wacc, casc, vtxz, 3))
781781
continue;
@@ -799,7 +799,7 @@ struct FlowGFWOmegaXi {
799799
continue;
800800
if (casc.dcaV0daughters() > cfgcasc_dcav0dau)
801801
continue;
802-
if (abs(casc.mLambda() - o2::constants::physics::MassLambda0) > cfgcasc_mlambdawindow)
802+
if (std::fabs(casc.mLambda() - o2::constants::physics::MassLambda0) > cfgcasc_mlambdawindow)
803803
continue;
804804
// track quality check
805805
if (bachelor.tpcNClsFound() < cfgtpcclusters)

0 commit comments

Comments
 (0)