|
29 | 29 | // david.dobrigkeit.chinellato@cern.ch |
30 | 30 | // |
31 | 31 |
|
32 | | -#include <Math/Vector4D.h> |
33 | | -#include <cmath> |
34 | | -#include <array> |
35 | | -#include <cstdlib> |
36 | | -#include <map> |
37 | | -#include <string> |
38 | | -#include <vector> |
39 | | - |
40 | | -#include <TFile.h> |
41 | | -#include <TH2F.h> |
42 | | -#include <TProfile.h> |
43 | | -#include <TLorentzVector.h> |
44 | | -#include <TPDGCode.h> |
45 | | -#include "Math/Vector3D.h" |
46 | | - |
47 | | -#include "DCAFitter/DCAFitterN.h" |
48 | | -#include "DataFormatsParameters/GRPMagField.h" |
49 | | -#include "DataFormatsParameters/GRPObject.h" |
50 | | -#include "DetectorsBase/GeometryManager.h" |
51 | | -#include "DetectorsBase/Propagator.h" |
52 | | -#include "Framework/runDataProcessing.h" |
53 | | -#include "Framework/AnalysisTask.h" |
54 | | -#include "Framework/AnalysisDataModel.h" |
55 | | -#include "Framework/ASoAHelpers.h" |
56 | | -#include "Framework/O2DatabasePDGPlugin.h" |
57 | | -#include "ReconstructionDataFormats/Track.h" |
58 | | -#include "CCDB/BasicCCDBManager.h" |
59 | | -#include "CommonConstants/PhysicsConstants.h" |
60 | | -#include "Common/Core/trackUtilities.h" |
61 | | -#include "PWGLF/DataModel/LFStrangenessTables.h" |
62 | 32 | #include "PWGLF/DataModel/LFStrangenessMLTables.h" |
63 | 33 | #include "PWGLF/DataModel/LFStrangenessPIDTables.h" |
64 | 34 | #include "PWGLF/DataModel/LFStrangenessTables.h" |
|
79 | 49 |
|
80 | 50 | #include "CCDB/BasicCCDBManager.h" |
81 | 51 | #include "CommonConstants/PhysicsConstants.h" |
| 52 | +#include "DCAFitter/DCAFitterN.h" |
| 53 | +#include "DataFormatsParameters/GRPMagField.h" |
| 54 | +#include "DataFormatsParameters/GRPObject.h" |
| 55 | +#include "DetectorsBase/GeometryManager.h" |
| 56 | +#include "DetectorsBase/Propagator.h" |
82 | 57 | #include "Framework/ASoAHelpers.h" |
83 | 58 | #include "Framework/AnalysisDataModel.h" |
84 | 59 | #include "Framework/AnalysisTask.h" |
85 | 60 | #include "Framework/O2DatabasePDGPlugin.h" |
86 | 61 | #include "Framework/runDataProcessing.h" |
87 | 62 | #include "ReconstructionDataFormats/Track.h" |
88 | 63 |
|
| 64 | +#include "Math/Vector3D.h" |
89 | 65 | #include <Math/Vector4D.h> |
90 | 66 | #include <TFile.h> |
91 | 67 | #include <TH2F.h> |
@@ -389,7 +365,7 @@ struct QuarkoniaToHyperons { |
389 | 365 | ConfigurableAxis axisHypPairPhi{"axisHypPairPhi", {180, 0.0f, constants::math::TwoPI}, "Hyperon pair azimuthal angle (rad)"}; |
390 | 366 | } axes; |
391 | 367 |
|
392 | | - o2::base::MatLayerCylSet* lut; // material LUT for DCA fitter |
| 368 | + o2::base::MatLayerCylSet* lut; // material LUT for DCA fitter |
393 | 369 | o2::vertexing::DCAFitterN<2> fitter; |
394 | 370 |
|
395 | 371 | // helper object |
@@ -815,8 +791,8 @@ struct QuarkoniaToHyperons { |
815 | 791 | histos.add("QA/XiXiBar/h3dMassXiXiBarVsDCAPair", "h3dMassXiXiBarVsDCAPair", kTH3F, {axes.axisDCAHypPair, axes.axisPt, axes.axisQuarkoniumMass}); |
816 | 792 | histos.add("QA/XiXiBar/h3dMassXiXiBarVsPairCosPA", "h3dMassXiXiBarVsPairCosPA", kTH3F, {axes.axisHypPairCosPA, axes.axisPt, axes.axisQuarkoniumMass}); |
817 | 793 | histos.add("QA/XiXiBar/h3dMassXiXiBarVsPairOpAngle", "h3dMassXiXiBarVsPairOpAngle", kTH3F, {axes.axisHypPairOpAngle, axes.axisPt, axes.axisQuarkoniumMass}); |
818 | | - histos.add("QA/XiXiBar/h3dMassXiXiBarVsPairEta", "h3dMassXiXiBarVsPairEta", kTH3F, {axes.axisHypPairEta, axes.axisPt, axes.axisQuarkoniumMass}); |
819 | | - histos.add("QA/XiXiBar/h3dMassXiXiBarVsPairPhi", "h3dMassXiXiBarVsPairPhi", kTH3F, {axes.axisHypPairPhi, axes.axisPt, axes.axisQuarkoniumMass}); |
| 794 | + histos.add("QA/XiXiBar/h3dMassXiXiBarVsPairEta", "h3dMassXiXiBarVsPairEta", kTH3F, {axes.axisHypPairEta, axes.axisPt, axes.axisQuarkoniumMass}); |
| 795 | + histos.add("QA/XiXiBar/h3dMassXiXiBarVsPairPhi", "h3dMassXiXiBarVsPairPhi", kTH3F, {axes.axisHypPairPhi, axes.axisPt, axes.axisQuarkoniumMass}); |
820 | 796 | histos.add("QA/LaLaBar/h3dDeltaEtaXiXiBarVsPairEta", "h3dDeltaEtaXiXiBarVsPairEta", kTH3F, {axes.axisHypPairEta, axes.axisPt, axes.axisHypPairEta}); |
821 | 797 | histos.add("QA/LaLaBar/h3dDeltaPhiXiXiBarVsPairPhi", "h3dDeltaPhiXiXiBarVsPairPhi", kTH3F, {axes.axisHypPairPhi, axes.axisPt, axes.axisHypPairPhi}); |
822 | 798 | } |
@@ -1082,11 +1058,13 @@ struct QuarkoniaToHyperons { |
1082 | 1058 | std::array<float, 3> antiHyperonMomentum = {0.0f, 0.0f, 0.0f}; |
1083 | 1059 | float DCADau = -999.f; |
1084 | 1060 | float CosPA = -1.f; |
1085 | | - float OpAngle = -999.f; |
1086 | | - float Eta() const { |
| 1061 | + float OpAngle = -999.f; |
| 1062 | + float Eta() const |
| 1063 | + { |
1087 | 1064 | return RecoDecay::eta(std::array{hyperonMomentum[0] + antiHyperonMomentum[0], hyperonMomentum[1] + antiHyperonMomentum[1], hyperonMomentum[2] + antiHyperonMomentum[2]}); |
1088 | 1065 | } |
1089 | | - float Phi() const { |
| 1066 | + float Phi() const |
| 1067 | + { |
1090 | 1068 | return RecoDecay::phi(std::array{hyperonMomentum[0] + antiHyperonMomentum[0], hyperonMomentum[1] + antiHyperonMomentum[1]}); |
1091 | 1069 | } |
1092 | 1070 | }; |
@@ -1141,7 +1119,7 @@ struct QuarkoniaToHyperons { |
1141 | 1119 | TVector3 hyp2Momentum(pairInfo.antiHyperonMomentum[0], pairInfo.antiHyperonMomentum[1], pairInfo.antiHyperonMomentum[2]); |
1142 | 1120 | pairInfo.OpAngle = hyp1Momentum.Angle(hyp2Momentum); |
1143 | 1121 |
|
1144 | | - if (d < 1e-5f) { // Parallel or nearly parallel lines |
| 1122 | + if (d < 1e-5f) { // Parallel or nearly parallel lines |
1145 | 1123 | pairInfo.X = pairInfo.Y = pairInfo.Z = -999.f; // should we use another dummy value? Perhaps 999.f? |
1146 | 1124 | return pairInfo; |
1147 | 1125 | } |
@@ -1221,7 +1199,6 @@ struct QuarkoniaToHyperons { |
1221 | 1199 | return pairInfo; |
1222 | 1200 | } |
1223 | 1201 |
|
1224 | | - |
1225 | 1202 | template <typename TCollision> |
1226 | 1203 | bool isEventAccepted(TCollision collision, bool fillHists) |
1227 | 1204 | // check whether the collision passes our collision selections |
|
0 commit comments