1515// \brief This code loops over MFT tracks and collisions and fills histograms
1616// useful to compute dNdeta
1717
18- #include " PWGMM/Mult/DataModel/bestCollisionTable.h "
19-
20- #include " Common/DataModel/Centrality.h "
21- #include " Common/DataModel/EventSelection.h "
22- #include " Common/DataModel/Multiplicity.h "
23- #include " Common/DataModel/TrackSelectionTables.h "
18+ #include < chrono >
19+ # include < cmath >
20+ #include < cstdint >
21+ #include < cstdlib >
22+ #include < unordered_set >
23+ #include < vector >
2424
25- #include " CommonConstants/MathConstants.h"
2625#include " Framework/ASoAHelpers.h"
2726#include " Framework/AnalysisDataModel.h"
2827#include " Framework/AnalysisTask.h"
2928#include " Framework/Configurable.h"
3029#include " Framework/O2DatabasePDGPlugin.h"
3130#include " Framework/RuntimeError.h"
3231#include " Framework/runDataProcessing.h"
32+
33+ #include " Common/DataModel/Centrality.h"
34+ #include " Common/DataModel/EventSelection.h"
35+ #include " Common/DataModel/Multiplicity.h"
36+ #include " Common/DataModel/TrackSelectionTables.h"
37+ #include " CommonConstants/MathConstants.h"
3338#include " MathUtils/Utils.h"
3439#include " ReconstructionDataFormats/GlobalTrackID.h"
3540
41+ #include " PWGMM/Mult/DataModel/bestCollisionTable.h"
3642#include " TFile.h"
3743
38- #include < chrono>
39- #include < cmath>
40- #include < cstdint>
41- #include < cstdlib>
42- #include < unordered_set>
43- #include < vector>
44-
4544using namespace o2 ;
4645using namespace o2 ::framework;
4746using namespace o2 ::framework::expressions;
@@ -446,13 +445,13 @@ struct PseudorapidityDensityMFT {
446445 {
447446
448447 std::vector<typename std::decay_t <decltype (collisions)>::iterator> cols;
449- for (auto & bc : bcs) {
448+ for (const auto & bc : bcs) {
450449 if (!useEvSel ||
451450 (useEvSel && ((bc.selection_bit (aod::evsel::kIsBBT0A ) &&
452451 bc.selection_bit (aod::evsel::kIsBBT0C )) != 0 ))) {
453452 registry.fill (HIST (" EventSelection" ), 8 ); // added 5->12
454453 cols.clear ();
455- for (auto & collision : collisions) {
454+ for (const auto & collision : collisions) {
456455 if (collision.has_foundBC ()) {
457456 if (collision.foundBCId () == bc.globalIndex ()) {
458457 cols.emplace_back (collision);
@@ -522,7 +521,7 @@ struct PseudorapidityDensityMFT {
522521 }
523522
524523 if (tracks.size () > 0 ) {
525- for (auto & track : tracks) {
524+ for (const auto & track : tracks) {
526525
527526 float phi = track.phi ();
528527 o2::math_utils::bringTo02Pi (phi);
@@ -570,7 +569,7 @@ struct PseudorapidityDensityMFT {
570569 if ((z >= cfgVzCut1) && (z <= cfgVzCut2)) {
571570 registry.fill (HIST (" Tracks/2Danalysis/EventsNtrkZvtx_all" ), Ntrk, z);
572571 registry.fill (HIST (" EventSelection" ), 2 .);
573- for (auto & retrack : retracks) {
572+ for (const auto & retrack : retracks) {
574573 auto track = retrack.mfttrack ();
575574 if ((cfgnEta1 < track.eta ()) && (track.eta () < cfgnEta2) && track.nClusters () >= cfgnCluster && retrack.ambDegree () > 0 ) {
576575 registry.fill (HIST (" Tracks/2Danalysis/EtaZvtx" ), track.eta (), z);
@@ -602,7 +601,7 @@ struct PseudorapidityDensityMFT {
602601 registry.fill (HIST (" EventSelection" ), 12 .);
603602 }
604603 int64_t i = 0.0 , j = 0.0 , k = 0.0 ;
605- for (auto & retrack : retracks) {
604+ for (const auto & retrack : retracks) {
606605 auto track = retrack.mfttrack ();
607606 if ((cfgnEta1 < track.eta ()) && (track.eta () < cfgnEta2) && track.nClusters () >= cfgnCluster && retrack.ambDegree () > 0 ) {
608607 registry.fill (HIST (" Tracks/2Danalysis/EtaZvtx_sel8" ), track.eta (), z);
@@ -616,7 +615,7 @@ struct PseudorapidityDensityMFT {
616615 if (midtracks.size () > 0 ) {
617616 registry.fill (HIST (" Tracks/2Danalysis/EventsNtrkZvtx_sel8_inelfwdgt0" ), Ntrk, z);
618617 }
619- for (auto & retrack : retracks) {
618+ for (const auto & retrack : retracks) {
620619 auto track = retrack.mfttrack ();
621620
622621 if ((cfgnEta1 < track.eta ()) && (track.eta () < cfgnEta2) && track.nClusters () >= cfgnCluster) {
@@ -762,7 +761,7 @@ struct PseudorapidityDensityMFT {
762761
763762 registry.fill (HIST (" Events/Centrality/NtrkZvtx" ), Ntrk, z, c);
764763
765- for (auto & track : tracks) {
764+ for (const auto & track : tracks) {
766765
767766 float phi = track.phi ();
768767 o2::math_utils::bringTo02Pi (phi);
@@ -810,7 +809,7 @@ struct PseudorapidityDensityMFT {
810809 auto perCollisionMCSample = mcSample->sliceByCached (
811810 aod::mcparticle::mcCollisionId, mcCollision.globalIndex (), cache);
812811 auto nCharged = 0 ;
813- for (auto & particle : perCollisionMCSample) {
812+ for (const auto & particle : perCollisionMCSample) {
814813 auto charge = 0 .;
815814 auto p = pdg->GetParticle (particle.pdgCode ());
816815 if (p != nullptr ) {
@@ -827,7 +826,7 @@ struct PseudorapidityDensityMFT {
827826 auto perCollisionMCSampleCentral = mcSampleCentral->sliceByCached (
828827 aod::mcparticle::mcCollisionId, mcCollision.globalIndex (), cache);
829828 auto nChargedCentral = 0 ;
830- for (auto & particle : perCollisionMCSampleCentral ) {
829+ for (const auto & particle : perCollisionMCSample ) {
831830 auto charge = 0 .;
832831 auto p = pdg->GetParticle (particle.pdgCode ());
833832 if (p != nullptr ) {
@@ -852,7 +851,7 @@ struct PseudorapidityDensityMFT {
852851
853852 LOGP (debug, " MC col {} has {} reco cols" , mcCollision.globalIndex (),
854853 collisions.size ());
855- for (auto & collision : collisions) {
854+ for (const auto & collision : collisions) {
856855 registry.fill (HIST (" EventEfficiency" ), 3 .);
857856 if (!disableITSROFCut && !collision.selection_bit (aod::evsel::kNoITSROFrameBorder )) {
858857 return ;
@@ -894,7 +893,7 @@ struct PseudorapidityDensityMFT {
894893 registry.fill (HIST (" EventsSplitMult" ), nCharged);
895894 }
896895 if ((mcCollision.posZ () >= cfgVzCut1) && (mcCollision.posZ () <= cfgVzCut2)) {
897- for (auto & particle : particles) {
896+ for (const auto & particle : particles) {
898897 auto p = pdg->GetParticle (particle.pdgCode ());
899898 auto charge = 0 ;
900899 if (p != nullptr ) {
@@ -949,7 +948,7 @@ struct PseudorapidityDensityMFT {
949948
950949 float c_gen = -1 ;
951950 bool atLeastOne = false ;
952- for (auto & collision : collisions) {
951+ for (const auto & collision : collisions) {
953952 float c_rec = -1 ;
954953 if constexpr (ExColsGenCent::template contains<aod::CentFT0Cs>()) {
955954 c_rec = collision.centFT0C ();
@@ -978,7 +977,7 @@ struct PseudorapidityDensityMFT {
978977 aod::mcparticle::mcCollisionId, mcCollision.globalIndex (), cache);
979978 auto nCharged = 0 ;
980979
981- for (auto & particle : perCollisionMCSample) {
980+ for (const auto & particle : perCollisionMCSample) {
982981 auto p = pdg->GetParticle (particle.pdgCode ());
983982 auto charge = 0 ;
984983 if (p != nullptr ) {
@@ -995,7 +994,7 @@ struct PseudorapidityDensityMFT {
995994 mcCollision.posZ (), c_gen);
996995 }
997996
998- for (auto & particle : particles) {
997+ for (const auto & particle : particles) {
999998 auto p = pdg->GetParticle (particle.pdgCode ());
1000999 auto charge = 0 ;
10011000 if (p != nullptr ) {
@@ -1031,7 +1030,7 @@ struct PseudorapidityDensityMFT {
10311030 MFTTracksLabeled const & tracks, aod::McParticles const &)
10321031 {
10331032 if (!useEvSel || (useEvSel && collision.sel8 ())) {
1034- for (auto & track : tracks) {
1033+ for (const auto & track : tracks) {
10351034 if (!track.has_mcParticle ()) {
10361035 continue ;
10371036 }
0 commit comments