@@ -1041,9 +1041,18 @@ struct lambdak0sflattenicity {
10411041 continue ;
10421042 }
10431043
1044+ float etap = mcParticle.eta ();
1045+ float phip = mcParticle.phi ();
1046+ isegment = 0 ;
10441047 for (int ieta = 0 ; ieta < 5 ; ieta++) {
10451048 etamax = maxEtaFV0 - ieta * detaFV0;
1049+ if (ieta == 0 ) {
1050+ etamax = maxEtaFV0;
1051+ }
10461052 etamin = maxEtaFV0 - (ieta + 1 ) * detaFV0;
1053+ if (ieta == 4 ) {
1054+ etamin = minEtaFV0;
1055+ }
10471056 nsectors = 8 ;
10481057 if (ieta == 4 ) {
10491058 nsectors = 16 ;
@@ -1052,7 +1061,7 @@ struct lambdak0sflattenicity {
10521061 minphi = iphi * 2.0 * TMath::Pi () / nsectors;
10531062 maxphi = (iphi + 1 ) * 2.0 * TMath::Pi () / nsectors;
10541063 dphi = TMath::Abs (maxphi - minphi);
1055- if (mcParticle. eta () >= etamin && mcParticle. eta () < etamax && mcParticle. phi () >= minphi && mcParticle. phi () < maxphi) {
1064+ if (etap >= etamin && etap < etamax && phip >= minphi && phip < maxphi) {
10561065 RhoLatticeFV0AMC[isegment] += 1.0 / TMath::Abs (dphi * detaFV0);
10571066 }
10581067 isegment++;
@@ -1061,7 +1070,7 @@ struct lambdak0sflattenicity {
10611070 }
10621071
10631072 flattenicity =
1064- GetFlatenicity ({RhoLatticeFV0AMC.data (), RhoLatticeFV0AMC.size ()});
1073+ 1.0 - GetFlatenicity ({RhoLatticeFV0AMC.data (), RhoLatticeFV0AMC.size ()});
10651074 return flattenicity;
10661075 }
10671076 // ====================== Flattenicity estimation ends =====================
@@ -1530,7 +1539,7 @@ struct lambdak0sflattenicity {
15301539 PROCESS_SWITCH (lambdak0sflattenicity, processDataRun3, " Process Run 3 Data" ,
15311540 false );
15321541 PROCESS_SWITCH (lambdak0sflattenicity, processRecMC,
1533- " Process Run 3 mc, reconstructed" , false );
1542+ " Process Run 3 mc, reconstructed" , true );
15341543 PROCESS_SWITCH (lambdak0sflattenicity, processGenMC,
15351544 " Process Run 3 mc, generated" , true );
15361545};
0 commit comments