|
14 | 14 | // This code loops over v0 photons for studying material budget. |
15 | 15 | // Please write to: daiki.sekihata@cern.ch |
16 | 16 |
|
17 | | -#include <cstring> |
18 | | -#include <iterator> |
| 17 | +#include "PWGEM/Dilepton/Utils/MCUtilities.h" |
| 18 | +#include "PWGEM/PhotonMeson/Core/CutsLibrary.h" |
| 19 | +#include "PWGEM/PhotonMeson/Core/HistogramsLibrary.h" |
| 20 | +#include "PWGEM/PhotonMeson/Core/PairCut.h" |
| 21 | +#include "PWGEM/PhotonMeson/Core/V0PhotonCut.h" |
| 22 | +#include "PWGEM/PhotonMeson/DataModel/gammaTables.h" |
| 23 | +#include "PWGEM/PhotonMeson/Utils/MCUtilities.h" |
| 24 | +#include "PWGEM/PhotonMeson/Utils/PairUtilities.h" |
| 25 | + |
| 26 | +#include "Common/Core/RecoDecay.h" |
| 27 | +#include "Common/DataModel/Centrality.h" |
| 28 | +#include "Common/DataModel/EventSelection.h" |
| 29 | +#include "Common/DataModel/TrackSelectionTables.h" |
19 | 30 |
|
20 | | -#include "TString.h" |
21 | | -#include "Math/Vector4D.h" |
22 | | -#include "Framework/runDataProcessing.h" |
23 | | -#include "Framework/AnalysisTask.h" |
24 | | -#include "Framework/AnalysisDataModel.h" |
25 | 31 | #include "Framework/ASoAHelpers.h" |
| 32 | +#include "Framework/AnalysisDataModel.h" |
| 33 | +#include "Framework/AnalysisTask.h" |
| 34 | +#include "Framework/runDataProcessing.h" |
26 | 35 | #include "ReconstructionDataFormats/Track.h" |
27 | | -#include "Common/DataModel/TrackSelectionTables.h" |
28 | | -#include "Common/DataModel/EventSelection.h" |
29 | | -#include "Common/DataModel/Centrality.h" |
30 | | -#include "Common/DataModel/PIDResponse.h" |
31 | | -#include "Common/Core/RecoDecay.h" |
32 | | -#include "PWGEM/PhotonMeson/DataModel/gammaTables.h" |
33 | | -#include "PWGEM/PhotonMeson/Utils/PairUtilities.h" |
34 | | -#include "PWGEM/PhotonMeson/Utils/MCUtilities.h" |
35 | | -#include "PWGEM/PhotonMeson/Core/V0PhotonCut.h" |
36 | | -#include "PWGEM/PhotonMeson/Core/PairCut.h" |
37 | | -#include "PWGEM/PhotonMeson/Core/CutsLibrary.h" |
38 | | -#include "PWGEM/PhotonMeson/Core/HistogramsLibrary.h" |
39 | | -#include "PWGEM/Dilepton/Utils/MCUtilities.h" |
| 36 | + |
| 37 | +#include "Math/Vector4D.h" |
| 38 | +#include "TString.h" |
| 39 | + |
| 40 | +#include <cstring> |
| 41 | +#include <iterator> |
40 | 42 |
|
41 | 43 | using namespace o2; |
42 | 44 | using namespace o2::aod; |
@@ -127,8 +129,8 @@ struct MaterialBudgetMC { |
127 | 129 | THashList* list_pair_subsys_paircut = reinterpret_cast<THashList*>(list_pair_subsys_photoncut->FindObject(pair_cut_name.data())); |
128 | 130 | o2::aod::pwgem::photon::histogram::DefineHistograms(list_pair_subsys_paircut, "material_budget_study", "Pair"); |
129 | 131 | } // end of cut3 loop pair cut |
130 | | - } // end of probecut loop |
131 | | - } // end of tagcut loop |
| 132 | + } // end of probecut loop |
| 133 | + } // end of tagcut loop |
132 | 134 | } |
133 | 135 |
|
134 | 136 | static constexpr std::string_view pairnames[9] = {"PCMPCM", "PHOSPHOS", "EMCEMC", "PCMPHOS", "PCMEMC", "PCMDalitzEE", "PCMDalitzMuMu", "PHOSEMC", "DalitzEEDalitzEE"}; |
@@ -280,7 +282,7 @@ struct MaterialBudgetMC { |
280 | 282 | reinterpret_cast<THnSparseF*>(list_v0->FindObject(cut.GetName())->FindObject("hs_conv_point"))->Fill(value); |
281 | 283 |
|
282 | 284 | } // end of photon loop |
283 | | - } // end of cut loop |
| 285 | + } // end of cut loop |
284 | 286 |
|
285 | 287 | } // end of collision loop |
286 | 288 | } |
@@ -372,11 +374,11 @@ struct MaterialBudgetMC { |
372 | 374 | value[5] = eta_cp2; |
373 | 375 | reinterpret_cast<THnSparseF*>(list_pair_ss->FindObject(Form("%s_%s", tagcut.GetName(), probecut.GetName()))->FindObject(paircut.GetName())->FindObject("hs_conv_point_same"))->Fill(value); |
374 | 376 | } // end of pair cut loop |
375 | | - } // end of g2 loop |
376 | | - } // end of g1 loop |
377 | | - } // end of probecut loop |
378 | | - } // end of tagcut loop |
379 | | - } // end of collision loop |
| 377 | + } // end of g2 loop |
| 378 | + } // end of g1 loop |
| 379 | + } // end of probecut loop |
| 380 | + } // end of tagcut loop |
| 381 | + } // end of collision loop |
380 | 382 | } |
381 | 383 |
|
382 | 384 | Partition<MyCollisions> grouped_collisions = cfgCentMin < o2::aod::cent::centFT0M && o2::aod::cent::centFT0M < cfgCentMax; // this goes to same event. |
|
0 commit comments