|
13 | 13 | /// \brief Analysis of D0/Lambda_c yield as a function of flattenicity |
14 | 14 | /// \author Laszlo Gyulai, laszlo.gyulai@cern.ch |
15 | 15 |
|
| 16 | +#include "RecoDecay.h" |
16 | 17 | #include "PWGHF/Core/CentralityEstimation.h" |
17 | 18 | #include "PWGHF/Core/DecayChannels.h" |
18 | 19 | #include "PWGHF/Core/HfHelper.h" |
|
23 | 24 | #include "PWGHF/DataModel/TrackIndexSkimmingTables.h" |
24 | 25 | #include "PWGHF/Utils/utilsEvSelHf.h" |
25 | 26 |
|
26 | | -#include "Common/DataModel/Centrality.h" |
27 | 27 | #include "Common/DataModel/EventSelection.h" |
28 | 28 | #include "Common/DataModel/Multiplicity.h" |
29 | | -#include "Common/DataModel/PIDResponseTOF.h" |
30 | | -#include "Common/DataModel/PIDResponseTPC.h" |
31 | 29 | #include "Common/DataModel/TrackSelectionTables.h" |
32 | 30 |
|
33 | 31 | #include <CCDB/BasicCCDBManager.h> |
| 32 | +#include <CommonConstants/MathConstants.h> |
| 33 | +#include <CommonConstants/PhysicsConstants.h> |
| 34 | +#include <DataFormatsFIT/Triggers.h> |
34 | 35 | #include <Framework/ASoA.h> |
35 | 36 | #include <Framework/AnalysisDataModel.h> |
| 37 | +#include <Framework/AnalysisHelpers.h> |
36 | 38 | #include <Framework/AnalysisTask.h> |
| 39 | +#include <Framework/Configurable.h> |
37 | 40 | #include <Framework/HistogramRegistry.h> |
38 | 41 | #include <Framework/HistogramSpec.h> |
| 42 | +#include <Framework/InitContext.h> |
39 | 43 | #include <Framework/runDataProcessing.h> |
40 | | - |
41 | | -#include <algorithm> |
| 44 | +#include <Rtypes.h> |
| 45 | +#include <THnSparse.h> |
| 46 | +#include <TPDGCode.h> |
| 47 | + |
| 48 | +#include <cstdint> |
| 49 | +#include <bitset> |
| 50 | +#include <cstddef> |
| 51 | +#include <cmath> |
42 | 52 | #include <map> |
43 | 53 | #include <string> |
44 | 54 | #include <vector> |
@@ -1039,6 +1049,7 @@ struct HfTaskFlattenicityD0Lc { |
1039 | 1049 | { |
1040 | 1050 | rhoLatticeFV0.fill(0); |
1041 | 1051 | rhoLatticeFV0Calibrated.fill(0); |
| 1052 | + float flat = 9999; |
1042 | 1053 | if (collision.has_foundFV0()) { |
1043 | 1054 | auto fv0 = collision.foundFV0(); |
1044 | 1055 | std::bitset<8> fV0Triggers = fv0.triggerMask(); |
@@ -1068,13 +1079,10 @@ struct HfTaskFlattenicityD0Lc { |
1068 | 1079 | registry.fill(HIST("Flattenicity_calibrated"), 1 - flattenicityFV0Calibrated); |
1069 | 1080 | registry.fill(HIST("Flattenicity"), 1 - flattenicityFV0); |
1070 | 1081 | } |
1071 | | - return 1. - flattenicityFV0; |
1072 | | - } else { |
1073 | | - return 9999; |
1074 | | - } |
1075 | | - } else { |
1076 | | - return 9999; |
| 1082 | + flat = 1. - flattenicityFV0; |
| 1083 | + } |
1077 | 1084 | } |
| 1085 | + return flat; |
1078 | 1086 | } |
1079 | 1087 |
|
1080 | 1088 | template <typename T, std::size_t S> |
|
0 commit comments