|
28 | 28 | /// |
29 | 29 | /// \author Alberto Calivà <alberto.caliva@cern.ch> |
30 | 30 |
|
31 | | -#include "Framework/AnalysisTask.h" |
32 | | -#include "Framework/Configurable.h" |
33 | | -#include "Framework/HistogramRegistry.h" |
34 | | -#include "Framework/InitContext.h" |
35 | | -#include "Framework/Logger.h" |
36 | | -#include "Framework/OutputObjHeader.h" |
37 | | -#include "Framework/runDataProcessing.h" |
| 31 | +#include <Framework/AnalysisTask.h> |
| 32 | +#include <Framework/Configurable.h> |
| 33 | +#include <Framework/HistogramRegistry.h> |
| 34 | +#include <Framework/InitContext.h> |
| 35 | +#include <Framework/Logger.h> |
| 36 | +#include <Framework/OutputObjHeader.h> |
| 37 | +#include <Framework/runDataProcessing.h> |
38 | 38 |
|
39 | 39 | #include <Math/Boost.h> |
40 | 40 | #include <Math/Vector3D.h> |
@@ -71,7 +71,7 @@ struct CoalescenceTreeProducer { |
71 | 71 |
|
72 | 72 | Preslice<aod::McParticles> mcParticlesPerMcCollision = aod::mcparticle::mcCollisionId; |
73 | 73 |
|
74 | | - OutputObj<TTree> treeBoundState{"treeBoundState"}; |
| 74 | + OutputObj<TTree> treeBoundState{"treeBoundState","treeBoundState"}; |
75 | 75 |
|
76 | 76 | int64_t eventID = 0; // Event ID |
77 | 77 | int64_t idB1 = 0, idB2 = 0, idB3 = 0; // MC particle IDs |
@@ -154,45 +154,46 @@ struct CoalescenceTreeProducer { |
154 | 154 | // Tree for pairs/triplets. |
155 | 155 | // For deuterons only the first two baryons are stored. |
156 | 156 | // For three-body states, the third-baryon branches are also created. |
157 | | - treeBoundState->Branch("eventID", &eventID); |
| 157 | + treeBoundState->Branch("eventID", &eventID, "eventID/L"); |
158 | 158 |
|
159 | | - treeBoundState->Branch("idB1", &idB1); |
160 | | - treeBoundState->Branch("idB2", &idB2); |
| 159 | + treeBoundState->Branch("idB1", &idB1, "idB1/L"); |
| 160 | + treeBoundState->Branch("idB2", &idB2, "idB2/L"); |
161 | 161 |
|
162 | | - treeBoundState->Branch("pdgB1", &pdgB1); |
163 | | - treeBoundState->Branch("pdgB2", &pdgB2); |
| 162 | + treeBoundState->Branch("pdgB1", &pdgB1, "pdgB1/I"); |
| 163 | + treeBoundState->Branch("pdgB2", &pdgB2, "pdgB2/I"); |
164 | 164 |
|
165 | | - treeBoundState->Branch("chargeB1", &chargeB1); |
166 | | - treeBoundState->Branch("chargeB2", &chargeB2); |
| 165 | + treeBoundState->Branch("chargeB1", &chargeB1, "chargeB1/I"); |
| 166 | + treeBoundState->Branch("chargeB2", &chargeB2, "chargeB2/I"); |
167 | 167 |
|
168 | | - treeBoundState->Branch("xB1", &xB1); |
169 | | - treeBoundState->Branch("yB1", &yB1); |
170 | | - treeBoundState->Branch("zB1", &zB1); |
171 | | - treeBoundState->Branch("tB1", &tB1); |
172 | | - treeBoundState->Branch("pxB1", &pxB1); |
173 | | - treeBoundState->Branch("pyB1", &pyB1); |
174 | | - treeBoundState->Branch("pzB1", &pzB1); |
| 168 | + treeBoundState->Branch("xB1", &xB1, "xB1/F"); |
| 169 | + treeBoundState->Branch("yB1", &yB1, "yB1/F"); |
| 170 | + treeBoundState->Branch("zB1", &zB1, "zB1/F"); |
| 171 | + treeBoundState->Branch("tB1", &tB1, "tB1/F"); |
| 172 | + treeBoundState->Branch("pxB1", &pxB1, "pxB1/F"); |
| 173 | + treeBoundState->Branch("pyB1", &pyB1, "pyB1/F"); |
| 174 | + treeBoundState->Branch("pzB1", &pzB1, "pzB1/F"); |
175 | 175 |
|
176 | | - treeBoundState->Branch("xB2", &xB2); |
177 | | - treeBoundState->Branch("yB2", &yB2); |
178 | | - treeBoundState->Branch("zB2", &zB2); |
179 | | - treeBoundState->Branch("tB2", &tB2); |
180 | | - treeBoundState->Branch("pxB2", &pxB2); |
181 | | - treeBoundState->Branch("pyB2", &pyB2); |
182 | | - treeBoundState->Branch("pzB2", &pzB2); |
| 176 | + treeBoundState->Branch("xB2", &xB2, "xB2/F"); |
| 177 | + treeBoundState->Branch("yB2", &yB2, "yB2/F"); |
| 178 | + treeBoundState->Branch("zB2", &zB2, "zB2/F"); |
| 179 | + treeBoundState->Branch("tB2", &tB2, "tB2/F"); |
| 180 | + treeBoundState->Branch("pxB2", &pxB2, "pxB2/F"); |
| 181 | + treeBoundState->Branch("pyB2", &pyB2, "pyB2/F"); |
| 182 | + treeBoundState->Branch("pzB2", &pzB2, "pzB2/F"); |
183 | 183 |
|
184 | 184 | if (static_cast<int>(boundStateSpecies) != kDeuteron) { |
185 | | - treeBoundState->Branch("idB3", &idB3); |
186 | | - treeBoundState->Branch("pdgB3", &pdgB3); |
187 | | - treeBoundState->Branch("chargeB3", &chargeB3); |
188 | | - |
189 | | - treeBoundState->Branch("xB3", &xB3); |
190 | | - treeBoundState->Branch("yB3", &yB3); |
191 | | - treeBoundState->Branch("zB3", &zB3); |
192 | | - treeBoundState->Branch("tB3", &tB3); |
193 | | - treeBoundState->Branch("pxB3", &pxB3); |
194 | | - treeBoundState->Branch("pyB3", &pyB3); |
195 | | - treeBoundState->Branch("pzB3", &pzB3); |
| 185 | + treeBoundState->Branch("idB3", &idB3, "idB3/L"); |
| 186 | + |
| 187 | + treeBoundState->Branch("pdgB3", &pdgB3, "pdgB3/I"); |
| 188 | + treeBoundState->Branch("chargeB3", &chargeB3, "chargeB3/I"); |
| 189 | + |
| 190 | + treeBoundState->Branch("xB3", &xB3, "xB3/F"); |
| 191 | + treeBoundState->Branch("yB3", &yB3, "yB3/F"); |
| 192 | + treeBoundState->Branch("zB3", &zB3, "zB3/F"); |
| 193 | + treeBoundState->Branch("tB3", &tB3, "tB3/F"); |
| 194 | + treeBoundState->Branch("pxB3", &pxB3, "pxB3/F"); |
| 195 | + treeBoundState->Branch("pyB3", &pyB3, "pyB3/F"); |
| 196 | + treeBoundState->Branch("pzB3", &pzB3, "pzB3/F"); |
196 | 197 | } |
197 | 198 | } |
198 | 199 |
|
|
0 commit comments