|
36 | 36 | #include <Framework/OutputObjHeader.h> |
37 | 37 | #include <Framework/runDataProcessing.h> |
38 | 38 |
|
39 | | -#include <TMath.h> |
40 | 39 | #include <TF1.h> |
| 40 | +#include <TMath.h> |
41 | 41 |
|
42 | 42 | #include <chrono> |
43 | 43 | #include <cstddef> |
@@ -149,18 +149,18 @@ struct EseTableProducer { |
149 | 149 | Configurable<bool> cfgIsGoodZvtxFT0vsPV{"cfgIsGoodZvtxFT0vsPV", true, "kIsGoodZvtxFT0vsPV"}; |
150 | 150 | Configurable<bool> cfgNoCollInTimeRangeStandard{"cfgNoCollInTimeRangeStandard", true, "kNoCollInTimeRangeStandard"}; |
151 | 151 | Configurable<bool> cfgIsVertexITSTPC{"cfgIsVertexITSTPC", true, "Selects collisions with at least one ITS-TPC track"}; |
152 | | - Configurable<bool> cfgIsGoodITSLayersAll{"cfgIsGoodITSLayersAll", true, "kIsGoodITSLayersAll"}; |
| 152 | + Configurable<bool> cfgIsGoodITSLayersAll{"cfgIsGoodITSLayersAll", true, "kIsGoodITSLayersAll"}; |
153 | 153 | Configurable<std::string> cfgEfficiency{"cfgEfficiency", "", "CCDB path to efficiency object"}; |
154 | 154 |
|
155 | 155 | // Cut Configurables |
156 | | - Configurable<std::vector<double>> cfgMultGlobalCutPars{"cfgMultGlobalCutPars", std::vector<double>{2272.16, -76.6932, 1.01204, -0.00631545, 1.59868e-05, 136.336, -4.97006, 0.121199, -0.0015921, 7.66197e-06}, "Global vs FT0C multiplicity cut parameter values"}; |
| 156 | + Configurable<std::vector<double>> cfgMultGlobalCutPars{"cfgMultGlobalCutPars", std::vector<double>{2272.16, -76.6932, 1.01204, -0.00631545, 1.59868e-05, 136.336, -4.97006, 0.121199, -0.0015921, 7.66197e-06}, "Global vs FT0C multiplicity cut parameter values"}; |
157 | 157 | Configurable<std::vector<double>> cfgMultPVCutPars{"cfgMultPVCutPars", std::vector<double>{3074.43, -106.192, 1.46176, -0.00968364, 2.61923e-05, 182.128, -7.43492, 0.193901, -0.00256715, 1.22594e-05}, "PV vs FT0C multiplicity cut parameter values"}; |
158 | 158 | Configurable<std::vector<double>> cfgMultGlobalPVCutPars{"cfgMultGlobalPVCutPars", std::vector<double>{-0.223013, 0.715849, 0.664242, 0.0829653, -0.000503733, 1.21185e-06}, "Global vs PV multiplicity cut parameter values"}; |
159 | 159 | Configurable<std::string> cfgMultCorrHighCutFunction{"cfgMultCorrHighCutFunction", "[0] + [1]*x + [2]*x*x + [3]*x*x*x + [4]*x*x*x*x + 3.*([5] + [6]*x + [7]*x*x + [8]*x*x*x + [9]*x*x*x*x)", "Functional for multiplicity correlation cut"}; |
160 | 160 | Configurable<std::string> cfgMultCorrLowCutFunction{"cfgMultCorrLowCutFunction", "[0] + [1]*x + [2]*x*x + [3]*x*x*x + [4]*x*x*x*x - 3.*([5] + [6]*x + [7]*x*x + [8]*x*x*x + [9]*x*x*x*x)", "Functional for multiplicity correlation cut"}; |
161 | 161 | Configurable<std::string> cfgMultGlobalPVCorrCutFunction{"cfgMultGlobalPVCorrCutFunction", "[0] + [1]*x + 3*([2] + [3]*x + [4]*x*x + [5]*x*x*x)", "Functional for global vs pv multiplicity correlation cut"}; |
162 | 162 | struct : ConfigurableGroup { |
163 | | - Configurable<std::string> cfgMultGlobalASideCorrCutFunction{"cfgMultGlobalASideCorrCutFunction","[0]+[1]*x+[2]*x*x+[3]*x*x*x+[4]*x*x*x*x + [10]*([5]+[6]*x+[7]*x*x+[8]*x*x*x+[9]*x*x*x*x)", "Functional for global vs V0A multiplicity low correlation cut"}; |
| 163 | + Configurable<std::string> cfgMultGlobalASideCorrCutFunction{"cfgMultGlobalASideCorrCutFunction", "[0]+[1]*x+[2]*x*x+[3]*x*x*x+[4]*x*x*x*x + [10]*([5]+[6]*x+[7]*x*x+[8]*x*x*x+[9]*x*x*x*x)", "Functional for global vs V0A multiplicity low correlation cut"}; |
164 | 164 | Configurable<std::vector<double>> cfgMultGlobalV0ACutPars{"cfgMultGlobalV0ACutPars", std::vector<double>{567.785, 172.715, 0.77888, -0.00693466, 1.40564e-05, 679.853, 66.8068, -0.444332, 0.00115002, -4.92064e-07}, "Global vs FV0A multiplicity cut parameter values"}; |
165 | 165 | Configurable<std::vector<double>> cfgMultGlobalT0ACutPars{"cfgMultGlobalT0ACutPars", std::vector<double>{241.618, 61.8402, 0.348049, -0.00306078, 6.20357e-06, 315.235, 29.1491, -0.188639, 0.00044528, -9.08912e-08}, "Global vs FT0A multiplicity cut parameter values"}; |
166 | 166 | Configurable<float> cfgGlobalV0ALowSigma{"cfgGlobalV0ALowSigma", -3.0f, "Number of sigma deviations below expected value in global vs V0A correlation"}; |
@@ -189,7 +189,6 @@ struct EseTableProducer { |
189 | 189 | EseTableProducer::multGlobalV0ACutPars = cfgGlobalAsideCorrCuts.cfgMultGlobalV0ACutPars; |
190 | 190 | EseTableProducer::multGlobalT0ACutPars = cfgGlobalAsideCorrCuts.cfgMultGlobalT0ACutPars; |
191 | 191 | EseTableProducer::multGlobalPVCorrCutPars = cfgMultGlobalPVCutPars; |
192 | | - |
193 | 192 |
|
194 | 193 | AxisSpec t0cAxis = {1000, 0, 10000, "N_{ch} (T0C)"}; |
195 | 194 | AxisSpec t0aAxis = {300, 0, 30000, "N_{ch} (T0A)"}; |
@@ -449,7 +448,7 @@ struct EseTableProducer { |
449 | 448 | float centrality; |
450 | 449 | int64_t multiplicity; |
451 | 450 | }; |
452 | | - |
| 451 | + |
453 | 452 | template <typename TCollision> |
454 | 453 | bool eventSelected(TCollision collision, const int& multTrk, const float& centrality) |
455 | 454 | { |
@@ -546,11 +545,11 @@ struct EseTableProducer { |
546 | 545 | return; |
547 | 546 | } |
548 | 547 | } |
549 | | - |
| 548 | + |
550 | 549 | const XAxis xaxis{collision.centFT0C(), tracks.size()}; |
551 | 550 | if (cfgUseAdditionalEventCut && !eventSelected(collision, xaxis.multiplicity, xaxis.centrality)) |
552 | 551 | return; |
553 | | - |
| 552 | + |
554 | 553 | registry.fill(HIST("hMeanPtStat"), 0.5); |
555 | 554 | const auto centrality = collision.centFT0C(); |
556 | 555 | const auto mean = calculateMeanPt(tracks, centrality); |
|
0 commit comments