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"
1716#include " PWGHF/Core/DecayChannels.h"
17+ #include " PWGHF/DataModel/CandidateReconstructionTables.h"
1818#include " PWGHF/DataModel/CandidateSelectionTables.h"
1919#include " PWGJE/Core/JetDerivedDataUtilities.h"
2020#include " PWGJE/Core/JetUtilities.h"
3737#include < Framework/runDataProcessing.h>
3838
3939#include < TH1.h>
40- #include " TVector3.h"
41-
40+ #include < TVector3.h>
4241
4342#include < cmath>
4443#include < string>
@@ -48,29 +47,26 @@ using namespace o2;
4847using namespace o2 ::framework;
4948using namespace o2 ::framework::expressions;
5049
51-
5250consteval float getValFromBin (int bin)
5351{
5452 return static_cast <float >(bin) - 0 .5f ;
5553}
5654
57- enum BinExpColCntr {AllCollisions = 1 ,
58- Sel8ZCut = 2 };
55+ enum BinExpColCntr { AllCollisions = 1 ,
56+ Sel8ZCut = 2 };
5957
60- enum BinExpJetCntr {ChargedJets = 1 };
61- enum BinMCColCntr {All = 1 ,
62- ZCut = 2 ,
63- Matched = 3 ,
64- MatchedSel8ZCut = 4
65- };
66-
67- enum BinMCJetCntr {DetectorLevelJetInMCCollision = 1 ,
68- ParticleLevelJetInMCCollision = 2 ,
69- DetectorLevelJetWithMatchedCandidate = 3 ,
70- ParticleLevelJetWithMatchedCandidate = 4
58+ enum BinExpJetCntr { ChargedJets = 1 };
59+ enum BinMCColCntr { All = 1 ,
60+ ZCut = 2 ,
61+ Matched = 3 ,
62+ MatchedSel8ZCut = 4
7163};
7264
73-
65+ enum BinMCJetCntr { DetectorLevelJetInMCCollision = 1 ,
66+ ParticleLevelJetInMCCollision = 2 ,
67+ DetectorLevelJetWithMatchedCandidate = 3 ,
68+ ParticleLevelJetWithMatchedCandidate = 4
69+ };
7470
7571struct JetDsSpecSubs {
7672
@@ -104,16 +100,17 @@ struct JetDsSpecSubs {
104100 int trackSelection = -1 ;
105101
106102 // Filters
107- // Filter jetCuts = aod::jet::pt > jetPtMin&& aod::jet::r == nround(jetR.node() * 100.0f);
108- // Filter collisionFilter = nabs(aod::jcollision::posZ) < vertexZCut;
103+ // Filter jetCuts = aod::jet::pt > jetPtMin&& aod::jet::r == nround(jetR.node() * 100.0f);
104+ // Filter collisionFilter = nabs(aod::jcollision::posZ) < vertexZCut;
109105
110106 // =============
111107 // Histograms
112108 // =============
113109
114110 HistogramRegistry registry{
115111 " registry" ,
116- { {" h_collisions" , " event status;event status;entries" , {HistType::kTH1F , {{10 , 0.0 , 10.0 }}}},
112+ {
113+ {" h_collisions" , " event status;event status;entries" , {HistType::kTH1F , {{10 , 0.0 , 10.0 }}}},
117114 {" h_collision_counter_data" , " ;event counter;entries" , {HistType::kTH1F , {{10 , 0 ., 10 .}}}},
118115 {" h_collision_counter_mcd" , " ;event counter;entries" , {HistType::kTH1F , {{10 , 0 ., 10 .}}}},
119116 {" h_collision_counter_mcp" , " ;event counter;entries" , {HistType::kTH1F , {{10 , 0 ., 10 .}}}},
@@ -203,7 +200,7 @@ struct JetDsSpecSubs {
203200 hMcp->GetXaxis ()->SetBinLabel (2 , " Ds particles" );
204201 hMcp->GetXaxis ()->SetBinLabel (3 , " Ds jets with >=1 particle" );
205202
206- // Labels
203+ // Labels
207204 auto mcCollisionCounter = registry.get <TH1 >(HIST (" McEffCol" ));
208205 mcCollisionCounter->GetXaxis ()->SetBinLabel (BinMCColCntr::All, " mccollisions" );
209206 mcCollisionCounter->GetXaxis ()->SetBinLabel (BinMCColCntr::ZCut, " z_cut" );
@@ -215,7 +212,6 @@ struct JetDsSpecSubs {
215212 jetCounter->GetXaxis ()->SetBinLabel (BinMCJetCntr::DetectorLevelJetInMCCollision, " detector level" );
216213 jetCounter->GetXaxis ()->SetBinLabel (BinMCJetCntr::DetectorLevelJetWithMatchedCandidate, " particle matched jets" );
217214 jetCounter->GetXaxis ()->SetBinLabel (BinMCJetCntr::ParticleLevelJetWithMatchedCandidate, " detector matched jets" );
218-
219215 }
220216 // ===============
221217 // Lambda compute
@@ -298,7 +294,6 @@ struct JetDsSpecSubs {
298294 }
299295 PROCESS_SWITCH (JetDsSpecSubs, processCollisions, " collision QA" , false );
300296
301-
302297 // ==============
303298 // DATA process
304299 // ==============
@@ -316,7 +311,7 @@ struct JetDsSpecSubs {
316311 registry.fill (HIST (" h_collision_counter_data" ), 3.0 );
317312
318313 for (const auto & jet : jets) {
319- registry.fill (HIST (" h_dsjet_counter_data" ), 0.5 ); // DsChargedJets entries
314+ registry.fill (HIST (" h_dsjet_counter_data" ), 0.5 ); // DsChargedJets entries
320315
321316 registry.fill (HIST (" h_jet_pt_data" ), jet.pt ());
322317 registry.fill (HIST (" h_jet_eta_data" ), jet.eta ());
@@ -329,7 +324,6 @@ struct JetDsSpecSubs {
329324
330325 const float mjet = computeJetMass (jetTracks);
331326
332-
333327 TVector3 jetVector (jet.px (), jet.py (), jet.pz ());
334328
335329 int nDsInJet = 0 ;
@@ -384,13 +378,11 @@ struct JetDsSpecSubs {
384378 if (lambda12 >= 0 .f ) {
385379 registry.fill (HIST (" h_ds_jet_lambda12_data" ), lambda12);
386380 }
387-
388381 }
389382 }
390383 }
391384 PROCESS_SWITCH (JetDsSpecSubs, processDataChargedSubstructure, " Data charged jets" , false );
392385
393-
394386 // ==============
395387 // MC function
396388 // ==============
@@ -477,9 +469,7 @@ struct JetDsSpecSubs {
477469 mcdjet.pt (),
478470 mcd_zParallel,
479471 mcd_deltaR);
480-
481472 }
482-
483473 }
484474 }
485475 }
@@ -508,8 +498,6 @@ struct JetDsSpecSubs {
508498 mcdCandidates,
509499 mcpCandidates,
510500 jettracks);
511-
512-
513501 }
514502 PROCESS_SWITCH (JetDsSpecSubs, processMonteCarloEfficiencyDs, " Non-matched and matched MC Ds and jets" , false );
515503};
0 commit comments