Skip to content

Commit 6e19be8

Browse files
alibuildalibuild
committed
Please consider the following formatting changes
1 parent a20670a commit 6e19be8

File tree

1 file changed

+41
-36
lines changed

1 file changed

+41
-36
lines changed

‎PWGCF/Flow/Tasks/flowGfwOmegaXi.cxx‎

Lines changed: 41 additions & 36 deletions
Original file line numberDiff line numberDiff line change
@@ -14,40 +14,45 @@
1414
/// \since Sep/13/2024
1515
/// \brief This task is to caculate V0s and cascades flow by GenericFramework
1616

17-
#include <CCDB/BasicCCDBManager.h>
18-
#include <vector>
19-
#include <string>
20-
#include <cmath>
21-
#include <memory>
22-
#include "Framework/runDataProcessing.h"
23-
#include "Framework/AnalysisTask.h"
24-
#include "Framework/ASoAHelpers.h"
25-
#include "Framework/RunningWorkflowInfo.h"
26-
#include "Framework/HistogramRegistry.h"
27-
#include "Common/DataModel/EventSelection.h"
17+
#include "GFW.h"
18+
#include "GFWCumulant.h"
19+
#include "GFWPowerArray.h"
20+
#include "GFWWeights.h"
21+
22+
#include "PWGLF/DataModel/LFStrangenessTables.h"
23+
#include "PWGMM/Mult/DataModel/Index.h"
24+
25+
#include "Common/Core/EventPlaneHelper.h"
2826
#include "Common/Core/TrackSelection.h"
29-
#include "Common/DataModel/TrackSelectionTables.h"
27+
#include "Common/Core/trackUtilities.h"
3028
#include "Common/DataModel/Centrality.h"
29+
#include "Common/DataModel/EventSelection.h"
30+
#include "Common/DataModel/Multiplicity.h"
3131
#include "Common/DataModel/PIDResponse.h"
3232
#include "Common/DataModel/PIDResponseITS.h"
33-
#include "Common/DataModel/Multiplicity.h"
34-
#include "GFWPowerArray.h"
35-
#include "GFW.h"
36-
#include "GFWCumulant.h"
37-
#include "GFWWeights.h"
3833
#include "Common/DataModel/Qvectors.h"
39-
#include "Common/Core/EventPlaneHelper.h"
40-
#include "ReconstructionDataFormats/Track.h"
34+
#include "Common/DataModel/TrackSelectionTables.h"
35+
4136
#include "CommonConstants/PhysicsConstants.h"
42-
#include "Common/Core/trackUtilities.h"
43-
#include "PWGLF/DataModel/LFStrangenessTables.h"
44-
#include "PWGMM/Mult/DataModel/Index.h"
37+
#include "Framework/ASoAHelpers.h"
38+
#include "Framework/AnalysisTask.h"
39+
#include "Framework/HistogramRegistry.h"
40+
#include "Framework/RunningWorkflowInfo.h"
41+
#include "Framework/runDataProcessing.h"
42+
#include "ReconstructionDataFormats/Track.h"
43+
#include <CCDB/BasicCCDBManager.h>
44+
4545
#include "TList.h"
46-
#include <TProfile.h>
47-
#include <TRandom3.h>
4846
#include <TF1.h>
4947
#include <TF2.h>
5048
#include <TPDGCode.h>
49+
#include <TProfile.h>
50+
#include <TRandom3.h>
51+
52+
#include <cmath>
53+
#include <memory>
54+
#include <string>
55+
#include <vector>
5156

5257
using namespace o2;
5358
using namespace o2::framework;
@@ -302,7 +307,7 @@ struct FlowGfwOmegaXi {
302307

303308
// QA
304309
if (cfgOutputQA) {
305-
// V0 QA
310+
// V0 QA
306311
registry.add("QAhisto/V0/hqaV0radiusbefore", "", {HistType::kTH1D, {{200, 0, 200}}});
307312
registry.add("QAhisto/V0/hqaV0radiusafter", "", {HistType::kTH1D, {{200, 0, 200}}});
308313
registry.add("QAhisto/V0/hqaV0cosPAbefore", "", {HistType::kTH1D, {{1000, 0.95, 1}}});
@@ -503,7 +508,7 @@ struct FlowGfwOmegaXi {
503508
fMultPVCutHigh = new TF1("fMultPVCutHigh", "[0]+[1]*x+[2]*x*x+[3]*x*x*x+[4]*x*x*x*x + 3.5*([5]+[6]*x+[7]*x*x+[8]*x*x*x+[9]*x*x*x*x)", 0, 100);
504509
fMultPVCutHigh->SetParameters(3074.43, -106.192, 1.46176, -0.00968364, 2.61923e-05, 182.128, -7.43492, 0.193901, -0.00256715, 1.22594e-05);
505510
}
506-
511+
507512
fT0AV0AMean = new TF1("fT0AV0AMean", "[0]+[1]*x", 0, 200000);
508513
fT0AV0AMean->SetParameters(-1601.0581, 9.417652e-01);
509514
fT0AV0ASigma = new TF1("fT0AV0ASigma", "[0]+[1]*x+[2]*x*x+[3]*x*x*x+[4]*x*x*x*x", 0, 200000);
@@ -905,7 +910,7 @@ struct FlowGfwOmegaXi {
905910
// check daughter TPC and TOF
906911
// K0short
907912
if (v0.pt() > trkQualityOpts.cfgCutPtK0sMin.value && v0.pt() < trkQualityOpts.cfgCutPtK0sMax.value) {
908-
if (v0.qtarm() / std::fabs(v0.alpha()) > v0BuilderOpts.cfgv0_ArmPodocut.value &&
913+
if (v0.qtarm() / std::fabs(v0.alpha()) > v0BuilderOpts.cfgv0_ArmPodocut.value &&
909914
std::fabs(v0.mK0Short() - o2::constants::physics::MassK0Short) < v0BuilderOpts.cfgv0_mk0swindow.value &&
910915
(std::fabs(v0posdau.tpcNSigmaPi()) < cfgNSigma[0] && std::fabs(v0negdau.tpcNSigmaPi()) < cfgNSigma[0]) &&
911916
((std::fabs(v0posdau.tofNSigmaPi()) < cfgNSigma[3] || v0posdau.pt() < lowpt) && (std::fabs(v0negdau.tofNSigmaPi()) < cfgNSigma[3] || v0negdau.pt() < lowpt)) &&
@@ -928,9 +933,9 @@ struct FlowGfwOmegaXi {
928933
isLambda = true;
929934
candNumAll[1] = candNumAll[1] + 1;
930935
} else if (std::fabs(v0.mLambda() - o2::constants::physics::MassLambda) < v0BuilderOpts.cfgv0_mlambdawindow.value &&
931-
(std::fabs(v0negdau.tpcNSigmaPr()) < cfgNSigma[1] && std::fabs(v0posdau.tpcNSigmaPi()) < cfgNSigma[0]) &&
932-
((std::fabs(v0negdau.tofNSigmaPr()) < cfgNSigma[4] || v0negdau.pt() < lowpt) && (std::fabs(v0posdau.tofNSigmaPi()) < cfgNSigma[3] || v0posdau.pt() < lowpt)) &&
933-
((std::fabs(itsResponse.nSigmaITS<o2::track::PID::Proton>(v0posdau)) < cfgNSigma[7]) || v0posdau.pt() < lowpt) && ((std::fabs(itsResponse.nSigmaITS<o2::track::PID::Pion>(v0negdau)) < cfgNSigma[6]) || v0negdau.pt() < lowpt)) {
936+
(std::fabs(v0negdau.tpcNSigmaPr()) < cfgNSigma[1] && std::fabs(v0posdau.tpcNSigmaPi()) < cfgNSigma[0]) &&
937+
((std::fabs(v0negdau.tofNSigmaPr()) < cfgNSigma[4] || v0negdau.pt() < lowpt) && (std::fabs(v0posdau.tofNSigmaPi()) < cfgNSigma[3] || v0posdau.pt() < lowpt)) &&
938+
((std::fabs(itsResponse.nSigmaITS<o2::track::PID::Proton>(v0posdau)) < cfgNSigma[7]) || v0posdau.pt() < lowpt) && ((std::fabs(itsResponse.nSigmaITS<o2::track::PID::Pion>(v0negdau)) < cfgNSigma[6]) || v0negdau.pt() < lowpt)) {
934939
registry.fill(HIST("InvMassLambda_all"), v0.pt(), v0.mLambda(), v0.eta(), cent);
935940
isLambda = true;
936941
candNumAll[1] = candNumAll[1] + 1;
@@ -1049,9 +1054,9 @@ struct FlowGfwOmegaXi {
10491054
isOmega = true;
10501055
candNumAll[3] = candNumAll[3] + 1;
10511056
} else if (casc.sign() > 0 && std::fabs(casc.yOmega()) < cfgCasc_rapidity &&
1052-
(std::fabs(bachelor.tpcNSigmaKa()) < cfgNSigma[2] && std::fabs(negdau.tpcNSigmaPr()) < cfgNSigma[1] && std::fabs(posdau.tpcNSigmaPi()) < cfgNSigma[0]) &&
1053-
((std::fabs(bachelor.tofNSigmaKa()) < cfgNSigma[5] || bachelor.pt() < lowpt) && (std::fabs(negdau.tofNSigmaPr()) < cfgNSigma[4] || negdau.pt() < lowpt) && (std::fabs(posdau.tofNSigmaPi()) < cfgNSigma[3] || posdau.pt() < lowpt)) &&
1054-
((std::fabs(itsResponse.nSigmaITS<o2::track::PID::Kaon>(bachelor)) < cfgNSigma[8]) || bachelor.pt() < lowpt) && ((std::fabs(itsResponse.nSigmaITS<o2::track::PID::Proton>(posdau)) < cfgNSigma[7]) || posdau.pt() < lowpt) && ((std::fabs(itsResponse.nSigmaITS<o2::track::PID::Pion>(negdau)) < cfgNSigma[6]) || negdau.pt() < lowpt)) {
1057+
(std::fabs(bachelor.tpcNSigmaKa()) < cfgNSigma[2] && std::fabs(negdau.tpcNSigmaPr()) < cfgNSigma[1] && std::fabs(posdau.tpcNSigmaPi()) < cfgNSigma[0]) &&
1058+
((std::fabs(bachelor.tofNSigmaKa()) < cfgNSigma[5] || bachelor.pt() < lowpt) && (std::fabs(negdau.tofNSigmaPr()) < cfgNSigma[4] || negdau.pt() < lowpt) && (std::fabs(posdau.tofNSigmaPi()) < cfgNSigma[3] || posdau.pt() < lowpt)) &&
1059+
((std::fabs(itsResponse.nSigmaITS<o2::track::PID::Kaon>(bachelor)) < cfgNSigma[8]) || bachelor.pt() < lowpt) && ((std::fabs(itsResponse.nSigmaITS<o2::track::PID::Proton>(posdau)) < cfgNSigma[7]) || posdau.pt() < lowpt) && ((std::fabs(itsResponse.nSigmaITS<o2::track::PID::Pion>(negdau)) < cfgNSigma[6]) || negdau.pt() < lowpt)) {
10551060
registry.fill(HIST("InvMassOmega_all"), casc.pt(), casc.mOmega(), casc.eta(), cent);
10561061
isOmega = true;
10571062
candNumAll[3] = candNumAll[3] + 1;
@@ -1067,9 +1072,9 @@ struct FlowGfwOmegaXi {
10671072
isXi = true;
10681073
candNumAll[2] = candNumAll[2] + 1;
10691074
} else if (casc.sign() > 0 && std::fabs(casc.yXi()) < cfgCasc_rapidity &&
1070-
(std::fabs(bachelor.tpcNSigmaPi()) < cfgNSigma[0] && std::fabs(negdau.tpcNSigmaPr()) < cfgNSigma[1] && std::fabs(posdau.tpcNSigmaPi()) < cfgNSigma[0]) &&
1071-
((std::fabs(bachelor.tofNSigmaPi()) < cfgNSigma[3] || bachelor.pt() < lowpt) && (std::fabs(negdau.tofNSigmaPr()) < cfgNSigma[4] || negdau.pt() < lowpt) && (std::fabs(posdau.tofNSigmaPi()) < cfgNSigma[3] || posdau.pt() < lowpt)) &&
1072-
((std::fabs(itsResponse.nSigmaITS<o2::track::PID::Pion>(bachelor)) < cfgNSigma[6]) || bachelor.pt() < lowpt) && ((std::fabs(itsResponse.nSigmaITS<o2::track::PID::Proton>(posdau)) < cfgNSigma[7]) || posdau.pt() < lowpt) && ((std::fabs(itsResponse.nSigmaITS<o2::track::PID::Pion>(negdau)) < cfgNSigma[6]) || negdau.pt() < lowpt)) {
1075+
(std::fabs(bachelor.tpcNSigmaPi()) < cfgNSigma[0] && std::fabs(negdau.tpcNSigmaPr()) < cfgNSigma[1] && std::fabs(posdau.tpcNSigmaPi()) < cfgNSigma[0]) &&
1076+
((std::fabs(bachelor.tofNSigmaPi()) < cfgNSigma[3] || bachelor.pt() < lowpt) && (std::fabs(negdau.tofNSigmaPr()) < cfgNSigma[4] || negdau.pt() < lowpt) && (std::fabs(posdau.tofNSigmaPi()) < cfgNSigma[3] || posdau.pt() < lowpt)) &&
1077+
((std::fabs(itsResponse.nSigmaITS<o2::track::PID::Pion>(bachelor)) < cfgNSigma[6]) || bachelor.pt() < lowpt) && ((std::fabs(itsResponse.nSigmaITS<o2::track::PID::Proton>(posdau)) < cfgNSigma[7]) || posdau.pt() < lowpt) && ((std::fabs(itsResponse.nSigmaITS<o2::track::PID::Pion>(negdau)) < cfgNSigma[6]) || negdau.pt() < lowpt)) {
10731078
registry.fill(HIST("InvMassXi_all"), casc.pt(), casc.mXi(), casc.eta(), cent);
10741079
isXi = true;
10751080
candNumAll[2] = candNumAll[2] + 1;

0 commit comments

Comments
 (0)