@@ -343,12 +343,12 @@ struct HfTaskElectronWeakBoson {
343343 // histogram for UE
344344 registry.add (" hRho" , " rho UE density" , kTH1F , {axisE});
345345 registry.add (" hSumERC" , " RC sumE" , kTH1F , {axisE});
346- registry.add (" hUE " , " UE vs. centrality" , kTH2F , {{axisCentrality}, {axisE}});
346+ registry.add (" hEnergyUE " , " UE vs. centrality" , kTH2F , {{axisCentrality}, {axisE}});
347347 }
348348
349349 double getIsolatedCluster (const o2::aod::EMCALCluster& cluster,
350350 const SelectedClusters& clusters,
351- float UE )
351+ float EnergyUE )
352352 {
353353 double energySum = 0.0 ;
354354 double energySum_excl = 0.0 ;
@@ -374,10 +374,10 @@ struct HfTaskElectronWeakBoson {
374374 }
375375 energySum_excl = energySum - cluster.energy ();
376376 if (energySum > 0 ) {
377- isoEnergy = (energySum_excl - UE ) / cluster.energy ();
377+ isoEnergy = (energySum_excl - EnergyUE ) / cluster.energy ();
378378 }
379379
380- // LOG(info) <<"clustE = " << cluster.energy() << " ; energySum = " << energySum << " ; nclust in Cone = " << nclustSum - 1 << " ; UE = " << UE << " ; isoEnergy = " << isoEnergy;
380+ // LOG(info) <<"clustE = " << cluster.energy() << " ; energySum = " << energySum << " ; nclust in Cone = " << nclustSum - 1 << " ; UE = " << EnergyUE << " ; isoEnergy = " << isoEnergy;
381381 registry.fill (HIST (" hIsolationEnergy" ), cluster.energy (), isoEnergy);
382382
383383 return (isoEnergy);
@@ -416,8 +416,8 @@ struct HfTaskElectronWeakBoson {
416416 }
417417 float estimateRhoRC (const SelectedClusters& clusters)
418418 {
419- const float R = rIsolation;
420- const float A = o2::constants::math::PI * R * R ;
419+ const float randomConeR = rIsolation;
420+ const float randomConeArea = o2::constants::math::PI * randomConeR * randomConeR ;
421421
422422 std::vector<float > sumErc;
423423 sumErc.reserve (nRandomCones);
@@ -427,7 +427,7 @@ struct HfTaskElectronWeakBoson {
427427 float etarc = rnd->Uniform (-etaEmcMax, etaEmcMax); // in EMCal acceptance
428428 float phirc = rnd->Uniform (phiEmcMin, phiEmcMax); // in EMCal acceptance
429429
430- float s = 0 ;
430+ float energySumRC = 0 ;
431431
432432 for (auto & c : clusters) {
433433 if (c.energy () > rcHardE)
@@ -436,13 +436,13 @@ struct HfTaskElectronWeakBoson {
436436 double dPhirc = phirc - c.phi ();
437437 dPhirc = RecoDecay::constrainAngle (dPhirc, -o2::constants::math::PI);
438438 double const deltaRrc = std::sqrt (dEtarc * dEtarc + dPhirc * dPhirc);
439- if (deltaRrc < R ) {
440- s += c.energy ();
439+ if (deltaRrc < randomConeR ) {
440+ energySumRC += c.energy ();
441441 }
442442 }
443443
444- registry.fill (HIST (" hSumERC" ), s );
445- sumErc.push_back (s );
444+ registry.fill (HIST (" hSumERC" ), energySumRC );
445+ sumErc.push_back (energySumRC );
446446 }
447447
448448 if (sumErc.empty ())
@@ -454,9 +454,9 @@ struct HfTaskElectronWeakBoson {
454454
455455 float median = sumErc[sumErc.size () / 2 ];
456456 // LOG(info) << "median = " << median;
457- registry.fill (HIST (" hRho" ), median / A );
457+ registry.fill (HIST (" hRho" ), median / randomConeArea );
458458
459- return median / A ;
459+ return median / randomConeArea ;
460460 }
461461
462462 void recoMassZee (const KFParticle& kfpIsoEle,
@@ -621,13 +621,13 @@ struct HfTaskElectronWeakBoson {
621621
622622 // UE estimate
623623 float rho = 0 .f ;
624- float UE = 0 .f ;
624+ float EnergyUE = 0 .f ;
625625
626626 if (useUEsub) {
627627 rho = estimateRhoRC (emcClusters);
628- UE = rho * static_cast <float >(o2::constants::math::PI * rIsolation * rIsolation);
629- registry.fill (HIST (" hUE " ), centrality, UE );
630- // LOG(info) << "UE = " << UE ;
628+ EnergyUE = rho * static_cast <float >(o2::constants::math::PI * rIsolation * rIsolation);
629+ registry.fill (HIST (" hEnergyUE " ), centrality, EnergyUE );
630+ // LOG(info) << "UE = " << EnergyUE ;
631631 }
632632
633633 // track loop
@@ -754,7 +754,7 @@ struct HfTaskElectronWeakBoson {
754754 eop = energyEmc / match.track_as <TrackEle>().p ();
755755 // LOG(info) << "eop = " << eop;
756756
757- isoEnergy = getIsolatedCluster (cluster, emcClusters, UE );
757+ isoEnergy = getIsolatedCluster (cluster, emcClusters, EnergyUE );
758758
759759 if (match.track_as <TrackEle>().pt () > ptTHnThresh && isTHnElectron) {
760760 registry.fill (HIST (" hTHnElectrons" ), match.track_as <TrackEle>().pt (), match.track_as <TrackEle>().tpcNSigmaEl (), m02Emc, eop, isoEnergy, isoMomentum, trackCount, track.eta (), track.tpcSignal ());
0 commit comments