Skip to content

Commit ef1cd22

Browse files
alibuildalibuild
committed
Please consider the following formatting changes
1 parent 9ca261c commit ef1cd22

1 file changed

Lines changed: 42 additions & 38 deletions

File tree

‎PWGJE/Tasks/jetShape.cxx‎

Lines changed: 42 additions & 38 deletions
Original file line numberDiff line numberDiff line change
@@ -13,28 +13,27 @@
1313
/// \author Yuto Nishida <yuto.nishida@cern.ch>
1414
/// \brief Task for measuring the dependence of the jet shape function rho(r) on the distance r from the jet axis.
1515

16-
#include <string>
17-
#include <vector>
18-
#include <cmath>
19-
20-
#include "Framework/ASoA.h"
21-
#include "Framework/AnalysisDataModel.h"
22-
#include "Framework/AnalysisTask.h"
23-
#include "Framework/HistogramRegistry.h"
16+
#include "PWGJE/Core/FastJetUtilities.h"
17+
#include "PWGJE/Core/JetDerivedDataUtilities.h"
18+
#include "PWGJE/Core/JetUtilities.h"
19+
#include "PWGJE/DataModel/Jet.h"
2420

2521
#include "Common/Core/RecoDecay.h"
2622
#include "Common/Core/TrackSelection.h"
2723
#include "Common/Core/TrackSelectionDefaults.h"
2824
#include "Common/DataModel/EventSelection.h"
2925
#include "Common/DataModel/TrackSelectionTables.h"
3026

31-
#include "PWGJE/Core/FastJetUtilities.h"
32-
#include "PWGJE/Core/JetUtilities.h"
33-
#include "PWGJE/Core/JetDerivedDataUtilities.h"
34-
#include "PWGJE/DataModel/Jet.h"
35-
27+
#include "Framework/ASoA.h"
28+
#include "Framework/AnalysisDataModel.h"
29+
#include "Framework/AnalysisTask.h"
30+
#include "Framework/HistogramRegistry.h"
3631
#include "Framework/runDataProcessing.h"
3732

33+
#include <cmath>
34+
#include <string>
35+
#include <vector>
36+
3837
using namespace o2;
3938
using namespace o2::framework;
4039
using namespace o2::framework::expressions;
@@ -59,7 +58,7 @@ struct JetShapeTask {
5958
Configurable<float> tpcNSigmaPiMax{"tpcNSigmaPiMax", 3.5f, "Max value of tpcNsigmaPion"};
6059

6160
HistogramRegistry registry{"registry",
62-
{{"tpcTofPi", "tpcTofPi", {HistType::kTHnSparseD, {{35, 0, pMax}, {nBinsNSigma, nSigmaMin, nSigmaMax},{nBinsDistance, 0, distanceMax}}}},
61+
{{"tpcTofPi", "tpcTofPi", {HistType::kTHnSparseD, {{35, 0, pMax}, {nBinsNSigma, nSigmaMin, nSigmaMax}, {nBinsDistance, 0, distanceMax}}}},
6362
{"tpcTofPr", "tpcTofPr", {HistType::kTHnSparseD, {{35, 0, pMax}, {nBinsNSigma, nSigmaMin, nSigmaMax}, {nBinsDistance, 0, distanceMax}}}},
6463
{"tpcTofPiOutOfJet", "tpcTofPiOutOfJet", {HistType::kTH2F, {{35, 0, pMax}, {nBinsNSigma, nSigmaMin, nSigmaMax}}}},
6564
{"tpcTofPrOutOfJet", "tpcTofPrOutOfJet", {HistType::kTH2F, {{35, 0, pMax}, {nBinsNSigma, nSigmaMin, nSigmaMax}}}},
@@ -79,7 +78,7 @@ struct JetShapeTask {
7978
{"trackEta", "trackEta", {HistType::kTH1F, {{100, -1, 1}}}},
8079
{"trackTpcNClsCrossedRows", "trackTpcNClsCrossedRows", {HistType::kTH1F, {{50, 0, 200}}}},
8180
{"trackDcaXY", "trackDcaXY", {HistType::kTH1F, {{40, -10, 10}}}},
82-
{"trackItsChi2NCl", "trackItsChi2NCl", {HistType::kTH1F, { {60, 0, 30}}}},
81+
{"trackItsChi2NCl", "trackItsChi2NCl", {HistType::kTH1F, {{60, 0, 30}}}},
8382
{"trackTpcChi2NCl", "trackTpcChi2NCl", {HistType::kTH1F, {{100, 0, 50}}}},
8483
{"trackTpcNClsFound", "trackTpcNClsFound", {HistType::kTH1F, {{100, 0, 200}}}},
8584
{"trackItsNCls", "trackItsNCls", {HistType::kTH1F, {{10, 0, 10}}}},
@@ -184,7 +183,6 @@ struct JetShapeTask {
184183
void processJetShape(soa::Filtered<soa::Join<aod::JetCollisions, aod::BkgChargedRhos>>::iterator const& collision, aod::JetTracks const& tracks, soa::Join<aod::ChargedJets, aod::ChargedJetConstituents> const& jets)
185184
{
186185

187-
188186
std::vector<float> ptDensity(distanceCategory->size() - 1, 0.f);
189187
std::vector<float> ptDensityBg1(distanceCategory->size() - 1, 0.f);
190188
std::vector<float> ptDensityBg2(distanceCategory->size() - 1, 0.f);
@@ -271,7 +269,7 @@ struct JetShapeTask {
271269
registry.fill(HIST("event/vertexz"), collision.posZ());
272270

273271
std::vector<typename std::decay_t<decltype(tracks)>::iterator> goodTracks;
274-
for (auto it = tracks.begin(); it != tracks.end() ; ++it){
272+
for (auto it = tracks.begin(); it != tracks.end(); ++it) {
275273
const auto& track = *it;
276274

277275
registry.fill(HIST("trackTpcNClsCrossedRows"), track.tpcNClsCrossedRows());
@@ -280,16 +278,23 @@ struct JetShapeTask {
280278
registry.fill(HIST("trackTpcChi2NCl"), track.tpcChi2NCl());
281279
registry.fill(HIST("trackTpcNClsFound"), track.tpcNClsFound());
282280
registry.fill(HIST("trackItsNCls"), track.itsNCls());
283-
registry.fill(HIST("trackEta"),track.eta());
284-
registry.fill(HIST("trackPhi"),track.phi());
281+
registry.fill(HIST("trackEta"), track.eta());
282+
registry.fill(HIST("trackPhi"), track.phi());
285283

286-
if (std::abs(track.eta()) > etaTrUp) continue;
287-
if (track.tpcNClsCrossedRows() < nclcrossTpcMin) continue;
288-
if (std::abs(track.dcaXY()) > dcaxyMax) continue;
289-
if (track.itsChi2NCl() > chi2ItsMax) continue;
290-
if (track.tpcChi2NCl() > chi2TpcMax) continue;
291-
if (track.tpcNClsFound() < nclTpcMin) continue;
292-
if (track.itsNCls() < nclItsMin) continue;
284+
if (std::abs(track.eta()) > etaTrUp)
285+
continue;
286+
if (track.tpcNClsCrossedRows() < nclcrossTpcMin)
287+
continue;
288+
if (std::abs(track.dcaXY()) > dcaxyMax)
289+
continue;
290+
if (track.itsChi2NCl() > chi2ItsMax)
291+
continue;
292+
if (track.tpcChi2NCl() > chi2TpcMax)
293+
continue;
294+
if (track.tpcNClsFound() < nclTpcMin)
295+
continue;
296+
if (track.itsNCls() < nclItsMin)
297+
continue;
293298

294299
goodTracks.push_back(it);
295300
}
@@ -301,7 +306,7 @@ struct JetShapeTask {
301306

302307
// tracks conditions
303308
for (const auto& track_it : goodTracks) {
304-
309+
305310
const auto& track = *track_it;
306311

307312
// PID check
@@ -337,15 +342,15 @@ struct JetShapeTask {
337342
if (distanceBg1 < jetR || distanceBg2 < jetR) {
338343
registry.fill(HIST("tpcDedxOutOfJet"), track.p(), track.tpcSignal());
339344

340-
if(std::abs(track.tofNSigmaPi()) < nSigmaTofCut){
345+
if (std::abs(track.tofNSigmaPi()) < nSigmaTofCut) {
341346
registry.fill(HIST("tpcTofPiOutOfJet"), track.p(), track.tpcNSigmaPi());
342-
if(track.tpcNSigmaPi()>tpcNSigmaPiMin && track.tpcNSigmaPi()<tpcNSigmaPiMax){
343-
registry.fill(HIST("pVsPtForPiOutOfJet"), track.p(), track.pt());
344-
}
347+
if (track.tpcNSigmaPi() > tpcNSigmaPiMin && track.tpcNSigmaPi() < tpcNSigmaPiMax) {
348+
registry.fill(HIST("pVsPtForPiOutOfJet"), track.p(), track.pt());
349+
}
345350
}
346-
if(std::abs(track.tofNSigmaPr()) < nSigmaTofCut){
351+
if (std::abs(track.tofNSigmaPr()) < nSigmaTofCut) {
347352
registry.fill(HIST("tpcTofPrOutOfJet"), track.p(), track.tpcNSigmaPr());
348-
if(track.tpcNSigmaPr()>tpcNSigmaPrMin && track.tpcNSigmaPr()<tpcNSigmaPrMax){
353+
if (track.tpcNSigmaPr() > tpcNSigmaPrMin && track.tpcNSigmaPr() < tpcNSigmaPrMax) {
349354
registry.fill(HIST("pVsPtForPrOutOfJet"), track.p(), track.pt());
350355
}
351356
}
@@ -355,16 +360,16 @@ struct JetShapeTask {
355360
registry.fill(HIST("tpcDedx"), track.p(), track.tpcSignal(), distance);
356361
registry.fill(HIST("tofBeta"), track.p(), track.beta());
357362

358-
if(std::abs(track.tofNSigmaPr()) < nSigmaTofCut){
363+
if (std::abs(track.tofNSigmaPr()) < nSigmaTofCut) {
359364
registry.fill(HIST("tpcTofPr"), track.p(), track.tpcNSigmaPr(), distance);
360-
if(track.tpcNSigmaPr()>tpcNSigmaPrMin && track.tpcNSigmaPr()<tpcNSigmaPrMax){
365+
if (track.tpcNSigmaPr() > tpcNSigmaPrMin && track.tpcNSigmaPr() < tpcNSigmaPrMax) {
361366
registry.fill(HIST("pVsPtForPr"), track.p(), track.pt(), distance);
362367
}
363368
}
364369

365-
if(std::abs(track.tofNSigmaPi()) < nSigmaTofCut){
370+
if (std::abs(track.tofNSigmaPi()) < nSigmaTofCut) {
366371
registry.fill(HIST("tpcTofPi"), track.p(), track.tpcNSigmaPi(), distance);
367-
if(track.tpcNSigmaPi()>tpcNSigmaPiMin && track.tpcNSigmaPi()<tpcNSigmaPiMax){
372+
if (track.tpcNSigmaPi() > tpcNSigmaPiMin && track.tpcNSigmaPi() < tpcNSigmaPiMax) {
368373
registry.fill(HIST("pVsPtForPi"), track.p(), track.pt(), distance);
369374
}
370375
}
@@ -375,4 +380,3 @@ struct JetShapeTask {
375380
};
376381

377382
WorkflowSpec defineDataProcessing(ConfigContext const& cfgc) { return WorkflowSpec{adaptAnalysisTask<JetShapeTask>(cfgc)}; }
378-

0 commit comments

Comments
 (0)