99// granted to it by virtue of its status as an Intergovernmental Organization
1010// or submit itself to any jurisdiction.
1111// /
12+ // / \file taskFlattenicityD0Lc.cxx
1213// / \brief Analysis of D0/Lambda_c yield as a function of flattenicity
1314// / \author Laszlo Gyulai, laszlo.gyulai@cern.ch
1415
@@ -61,11 +62,11 @@ enum CandTypeSel {
6162} // namespace
6263
6364static const int nCellsFV0 = 48 ;
64- static const int CinnerFV0 = 32 ;
65+ static const int nInnerCellsFV0 = 32 ;
6566std::array<float , nCellsFV0> rhoLatticeFV0{};
6667std::array<float , nCellsFV0> fv0AmplitudeWoCalib{};
6768std::array<float , nCellsFV0> calib = {1.01697 , 1.122 , 1.03854 , 1.108 , 1.11634 , 1.14971 , 1.19321 , 1.06866 , 0.954675 , 0.952695 , 0.969853 , 0.957557 , 0.989784 , 1.01549 , 1.02182 , 0.976005 , 1.01865 , 1.06871 , 1.06264 , 1.02969 , 1.07378 , 1.06622 , 1.15057 , 1.0433 , 0.83654 , 0.847178 , 0.890027 , 0.920814 , 0.888271 , 1.04662 , 0.8869 , 0.856348 , 0.863181 , 0.906312 , 0.902166 , 1.00122 , 1.03303 , 0.887866 , 0.892437 , 0.906278 , 0.884976 , 0.864251 , 0.917221 , 1.10618 , 1.04028 , 0.893184 , 0.915734 , 0.892676 };
68- std::map<int , int > ChannelsToRings = {{0 , 0 }, {1 , 1 }, {2 , 2 }, {3 , 3 }, {4 , 7 }, {5 , 6 }, {6 , 5 }, {7 , 4 }, {8 , 8 }, {9 , 9 }, {10 , 10 }, {11 , 11 }, {12 , 15 }, {13 , 14 }, {14 , 13 }, {15 , 12 }, {16 , 16 }, {17 , 17 }, {18 , 18 }, {19 , 19 }, {20 , 23 }, {21 , 22 }, {22 , 21 }, {23 , 20 }, {24 , 24 }, {25 , 25 }, {26 , 26 }, {27 , 27 }, {28 , 31 }, {29 , 30 }, {30 , 29 }, {31 , 28 }, {32 , 32 }, {33 , 34 }, {34 , 36 }, {35 , 38 }, {36 , 47 }, {37 , 45 }, {38 , 43 }, {39 , 41 }, {40 , 33 }, {41 , 35 }, {42 , 37 }, {43 , 39 }, {44 , 46 }, {45 , 44 }, {46 , 42 }, {47 , 40 }};
69+ std::map<int , int > channelsToRings = {{0 , 0 }, {1 , 1 }, {2 , 2 }, {3 , 3 }, {4 , 7 }, {5 , 6 }, {6 , 5 }, {7 , 4 }, {8 , 8 }, {9 , 9 }, {10 , 10 }, {11 , 11 }, {12 , 15 }, {13 , 14 }, {14 , 13 }, {15 , 12 }, {16 , 16 }, {17 , 17 }, {18 , 18 }, {19 , 19 }, {20 , 23 }, {21 , 22 }, {22 , 21 }, {23 , 20 }, {24 , 24 }, {25 , 25 }, {26 , 26 }, {27 , 27 }, {28 , 31 }, {29 , 30 }, {30 , 29 }, {31 , 28 }, {32 , 32 }, {33 , 34 }, {34 , 36 }, {35 , 38 }, {36 , 47 }, {37 , 45 }, {38 , 43 }, {39 , 41 }, {40 , 33 }, {41 , 35 }, {42 , 37 }, {43 , 39 }, {44 , 46 }, {45 , 44 }, {46 , 42 }, {47 , 40 }};
6970
7071struct HfTaskFlattenicityD0Lc {
7172 Configurable<int > selectionFlagD0{" selectionFlagD0" , 1 , " Selection Flag for D0" };
@@ -478,7 +479,7 @@ struct HfTaskFlattenicityD0Lc {
478479 }
479480
480481 const float flat = fillFlat<true >(collision, 0 );
481- const float flat_calibrated = fillFlat<true >(collision, 1 );
482+ const float flatCalibrated = fillFlat<true >(collision, 1 );
482483
483484 const auto thisCollId = collision.globalIndex ();
484485
@@ -638,7 +639,7 @@ struct HfTaskFlattenicityD0Lc {
638639 auto collision = candidate.template collision_as <CollType>();
639640
640641 const float flat = fillFlat<false >(collision, 0 );
641- const float flat_calibrated = fillFlat<false >(collision, 1 );
642+ const float flatCalibrated = fillFlat<false >(collision, 1 );
642643
643644 float massD0{0 .f }, massD0bar{0 .f };
644645 massD0 = HfHelper::invMassD0ToPiK (candidate);
@@ -851,7 +852,7 @@ struct HfTaskFlattenicityD0Lc {
851852 const auto & groupedLcCandidates = candidatesLc.sliceBy (candLcPerCollision, thisCollId);
852853
853854 const float flat = fillFlat<true >(collision, 0 );
854- const float flat_calibrated = fillFlat<true >(collision, 1 );
855+ const float flatCalibrated = fillFlat<true >(collision, 1 );
855856
856857 for (const auto & candidate : groupedLcCandidates) {
857858 if (!(candidate.hfflag () & 1 << aod::hf_cand_3prong::DecayType::LcToPKPi)) {
@@ -1061,14 +1062,14 @@ struct HfTaskFlattenicityD0Lc {
10611062 for (std::size_t ich = 0 ; ich < fv0.channel ().size (); ich++) {
10621063 float amplCh = fv0.amplitude ()[ich];
10631064 int chv0 = fv0.channel ()[ich];
1064- int chv0phi = ChannelsToRings .at (chv0);
1065+ int chv0phi = channelsToRings .at (chv0);
10651066 if (amplCh > 0.0 ) {
10661067 if (chv0phi > 0.0 ) {
10671068 fv0AmplitudeWoCalib[chv0phi] = amplCh;
10681069 if (ifCalib) {
10691070 amplCh *= calib[chv0phi];
10701071 }
1071- if (chv0 < CinnerFV0 ) {
1072+ if (chv0 < nInnerCellsFV0 ) {
10721073 rhoLatticeFV0[chv0phi] += amplCh;
10731074 } else {
10741075 rhoLatticeFV0[chv0phi] += amplCh / 2 .;
@@ -1100,11 +1101,9 @@ struct HfTaskFlattenicityD0Lc {
11001101 int entries = signals.size ();
11011102 float flat{-1 };
11021103 float mRho {0 };
1103- float mRho_debug {0 };
11041104 for (int iCell = 0 ; iCell < entries; ++iCell) {
11051105 if (signals[iCell] > 0.0 ) {
11061106 mRho += 1.0 * signals[iCell];
1107- mRho_debug += 1.0 * signals[iCell];
11081107 }
11091108 }
11101109 mRho /= (1.0 * entries);
@@ -1132,4 +1131,4 @@ struct HfTaskFlattenicityD0Lc {
11321131WorkflowSpec defineDataProcessing (ConfigContext const & cfgc)
11331132{
11341133 return WorkflowSpec{adaptAnalysisTask<HfTaskFlattenicityD0Lc>(cfgc)};
1135- }
1134+ }
0 commit comments