|
12 | 12 | /// \brief this is a starting point for the Resonances tutorial |
13 | 13 | /// \author junlee kim |
14 | 14 |
|
| 15 | +#include "PWGLF/DataModel/ReducedHeptaQuarkTables.h" |
| 16 | + |
| 17 | +#include "Common/Core/trackUtilities.h" |
| 18 | + |
| 19 | +#include "CommonConstants/PhysicsConstants.h" |
| 20 | +#include "Framework/ASoAHelpers.h" |
| 21 | +#include "Framework/AnalysisDataModel.h" |
| 22 | +#include "Framework/AnalysisTask.h" |
| 23 | +#include "Framework/StepTHn.h" |
| 24 | +#include "Framework/runDataProcessing.h" |
15 | 25 | #include <Framework/Configurable.h> |
16 | | -#include <TLorentzVector.h> |
| 26 | + |
17 | 27 | #include <Math/GenVector/Boost.h> |
18 | | -#include <Math/Vector4D.h> |
19 | 28 | #include <Math/Vector3D.h> |
| 29 | +#include <Math/Vector4D.h> |
| 30 | +#include <TLorentzVector.h> |
20 | 31 | #include <TMath.h> |
21 | 32 | #include <TRandom3.h> |
| 33 | + |
22 | 34 | #include <fairlogger/Logger.h> |
| 35 | + |
| 36 | +#include <algorithm> |
23 | 37 | #include <iostream> |
24 | 38 | #include <iterator> |
25 | 39 | #include <string> |
26 | 40 | #include <vector> |
27 | | -#include <algorithm> |
28 | | - |
29 | | -#include "Framework/AnalysisTask.h" |
30 | | -#include "Framework/ASoAHelpers.h" |
31 | | -#include "Framework/runDataProcessing.h" |
32 | | -#include "Framework/AnalysisDataModel.h" |
33 | | -#include "Framework/StepTHn.h" |
34 | | -#include "Common/Core/trackUtilities.h" |
35 | | -#include "PWGLF/DataModel/ReducedHeptaQuarkTables.h" |
36 | | -#include "CommonConstants/PhysicsConstants.h" |
37 | 41 |
|
38 | 42 | using namespace o2; |
39 | 43 | using namespace o2::framework; |
@@ -203,15 +207,15 @@ struct heptaquark { |
203 | 207 | double cosHel = hqlambda_boost.Vect().Dot(hqphipair_boost.Vect()) / (hqlambda_boost.Vect().Mag() * hqphipair_boost.Vect().Mag()); |
204 | 208 | if (std::abs(cosHel) > cfgCollinear) |
205 | 209 | selection += 4; |
206 | | -/* |
207 | | - ROOT::Math::XYZVector rPV(col.posX(), col.posY(), col.posZ()); |
208 | | - ROOT::Math::XYZVector rSV(hq3.hqx(), hq3.hqy(), hq3.hqz()); |
209 | | - ROOT::Math::XYZVector L = rSV - rPV; |
210 | | - ROOT::Math::XYZVector exMom(ex.Px(), ex.Py(), ex.Pz()); |
211 | | - double cosPoint = L.Dot(exMom) / (L.R() * pEx.R() + 1e-9); |
212 | | - if (cosPoint < cfgCosPoint) |
213 | | - return 8; |
214 | | -*/ |
| 210 | + /* |
| 211 | + ROOT::Math::XYZVector rPV(col.posX(), col.posY(), col.posZ()); |
| 212 | + ROOT::Math::XYZVector rSV(hq3.hqx(), hq3.hqy(), hq3.hqz()); |
| 213 | + ROOT::Math::XYZVector L = rSV - rPV; |
| 214 | + ROOT::Math::XYZVector exMom(ex.Px(), ex.Py(), ex.Pz()); |
| 215 | + double cosPoint = L.Dot(exMom) / (L.R() * pEx.R() + 1e-9); |
| 216 | + if (cosPoint < cfgCosPoint) |
| 217 | + return 8; |
| 218 | + */ |
215 | 219 | return selection; |
216 | 220 | } |
217 | 221 |
|
|
0 commit comments