1313// / \brief Ds-tagged jet analysis with substructure histogram outputs
1414// / \author Monalisa Melo <monalisa.melo@cern.ch>, Universidade de São Paulo
1515
16+ #include " PWGHF/DataModel/CandidateReconstructionTables.h"
17+ #include " PWGHF/DataModel/CandidateSelectionTables.h"
18+ #include " PWGJE/Core/JetDerivedDataUtilities.h"
19+ #include " PWGJE/Core/JetUtilities.h"
20+ #include " PWGJE/DataModel/Jet.h"
21+ #include " PWGJE/DataModel/JetReducedData.h"
22+ #include " PWGJE/DataModel/JetSubstructure.h"
23+
24+ #include " Common/Core/RecoDecay.h"
25+
1626#include < Framework/ASoA.h>
1727#include < Framework/AnalysisDataModel.h>
1828#include < Framework/AnalysisTask.h>
2333#include < Framework/InitContext.h>
2434#include < Framework/runDataProcessing.h>
2535
26- #include " Common/Core/RecoDecay.h"
27- #include " PWGJE/Core/JetUtilities.h"
28- #include " PWGJE/Core/JetDerivedDataUtilities.h"
29-
30- #include " PWGJE/DataModel/Jet.h"
31- #include " PWGJE/DataModel/JetReducedData.h"
32- #include " PWGJE/DataModel/JetSubstructure.h"
33-
34- #include " PWGHF/DataModel/CandidateReconstructionTables.h"
35- #include " PWGHF/DataModel/CandidateSelectionTables.h"
36+ #include < Math/Vector4D.h>
37+ #include < TH1.h>
38+ #include < TVector3.h>
3639
3740#include < cmath>
38- #include < vector>
3941#include < string>
40-
41- #include < TH1.h>
42- #include < TVector3.h>
43- #include < Math/Vector4D.h>
42+ #include < vector>
4443
4544using namespace o2 ;
4645using namespace o2 ::framework;
@@ -73,7 +72,7 @@ struct JetDsSpecSubs {
7372 int trackSelection = -1 ;
7473
7574 // Filters
76- Filter jetCuts = aod::jet::pt > jetPtMin && aod::jet::r == nround(jetR.node() * 100 .0f );
75+ Filter jetCuts = aod::jet::pt > jetPtMin&& aod::jet::r == nround(jetR.node() * 100 .0f );
7776 Filter collisionFilter = nabs(aod::jcollision::posZ) < vertexZCut;
7877
7978 // =============
@@ -108,31 +107,11 @@ struct JetDsSpecSubs {
108107 {" h_ds_jet_lambda12" , " ;#lambda_{2}^{1};entries" , {HistType::kTH1F , {{200 , 0 ., 1 .}}}},
109108 {" h_ds_jet_mass" , " ;m_{jet};entries" , {HistType::kTH1F , {{200 , 0 ., 50 .}}}},
110109
111- {" hSparse_ds" , " ;m_{D_{S}};p_{T,D_{S}};p_{T,jet};z_{||};#DeltaR" ,
112- {HistType::kTHnSparseF ,
113- {
114- {200 , 1.7 , 2.2 },
115- {200 , 0 ., 100 .},
116- {200 , 0 ., 100 .},
117- {200 , 0 ., 2 .},
118- {200 , 0 ., 1 .}
119- }}},
120-
121- {" h2_response_jet_pt" , " ;p_{T}^{det};p_{T}^{part}" ,
122- {HistType::kTH2F ,
123- {
124- {200 , 0 ., 100 .},
125- {200 , 0 ., 100 .}
126- }}},
127-
128- {" h2_response_lambda11" , " ;#lambda_{1}^{1,det};#lambda_{1}^{1,part}" ,
129- {HistType::kTH2F ,
130- {
131- {200 , 0 ., 1 .},
132- {200 , 0 ., 1 .}
133- }}}
134- }
135- };
110+ {" hSparse_ds" , " ;m_{D_{S}};p_{T,D_{S}};p_{T,jet};z_{||};#DeltaR" , {HistType::kTHnSparseF , {{200 , 1.7 , 2.2 }, {200 , 0 ., 100 .}, {200 , 0 ., 100 .}, {200 , 0 ., 2 .}, {200 , 0 ., 1 .}}}},
111+
112+ {" h2_response_jet_pt" , " ;p_{T}^{det};p_{T}^{part}" , {HistType::kTH2F , {{200 , 0 ., 100 .}, {200 , 0 ., 100 .}}}},
113+
114+ {" h2_response_lambda11" , " ;#lambda_{1}^{1,det};#lambda_{1}^{1,part}" , {HistType::kTH2F , {{200 , 0 ., 1 .}, {200 , 0 ., 1 .}}}}}};
136115
137116 // ========
138117 // INIT
@@ -153,7 +132,7 @@ struct JetDsSpecSubs {
153132 // Lambda compute
154133 // ===============
155134
156- template <typename JET , typename TRACKS >
135+ template <typename JET , typename TRACKS >
157136 float computeLambda (JET const & jet,
158137 TRACKS const & tracks,
159138 float alpha,
@@ -164,7 +143,7 @@ struct JetDsSpecSubs {
164143 }
165144
166145 float sum = 0 .f ;
167-
146+
168147 for (auto const & trk : tracks) {
169148 const float dr = jetutilities::deltaR (jet, trk);
170149 sum += std::pow (trk.pt (), kappa) * std::pow (dr, alpha);
@@ -182,7 +161,7 @@ struct JetDsSpecSubs {
182161 // Jet Mass compute
183162 // =================
184163
185- template <typename TRACKS >
164+ template <typename TRACKS >
186165 float computeJetMass (TRACKS const & tracks)
187166 {
188167 double sumPx = 0.0 , sumPy = 0.0 , sumPz = 0.0 , sumE = 0.0 ;
@@ -240,7 +219,7 @@ struct JetDsSpecSubs {
240219 // DATA analysis
241220 // ===============
242221
243- template <typename JETS , typename CANDS >
222+ template <typename JETS , typename CANDS >
244223 void analyseData (JETS const & jets,
245224 CANDS const &)
246225 {
@@ -306,7 +285,7 @@ struct JetDsSpecSubs {
306285 // MCD analysis
307286 // ==============
308287
309- template <typename JETS >
288+ template <typename JETS >
310289 void analyseMCD (JETS const & jets)
311290 {
312291 for (const auto & jet : jets) {
@@ -320,7 +299,7 @@ struct JetDsSpecSubs {
320299
321300 const float mjet = computeJetMass (jetTracks);
322301
323- TVector3 jetVector (jet.px (), jet.py (),jet.pz ());
302+ TVector3 jetVector (jet.px (), jet.py (), jet.pz ());
324303
325304 for (const auto & ds : jet.template candidates_as <DsCandidatesMCD>()) {
326305
@@ -348,7 +327,7 @@ struct JetDsSpecSubs {
348327 // MC particle level
349328 // ==================
350329
351- template <typename JETS >
330+ template <typename JETS >
352331 void analyseMCP (JETS const & jets)
353332 {
354333 for (const auto & jet : jets) {
@@ -443,4 +422,4 @@ struct JetDsSpecSubs {
443422WorkflowSpec defineDataProcessing (ConfigContext const & cfgc)
444423{
445424 return WorkflowSpec{adaptAnalysisTask<JetDsSpecSubs>(cfgc)};
446- }
425+ }
0 commit comments