CaloClusterKineHelper Class Reference

#include <CaloClusterKineHelper.h>

Collaboration diagram for CaloClusterKineHelper:

Static Public Member Functions
static void	calculateKine (xAOD::CaloCluster *clu, const bool useweight=true, const bool updateLayers=true, const bool useGPUCriteria=false)
	Helper class to calculate cluster kinematics based on cells.

Detailed Description

Definition at line 13 of file CaloClusterKineHelper.h.

Member Function Documentation

calculateKine()

void CaloClusterKineHelper::calculateKine ( xAOD::CaloCluster * clu, const bool useweight = true, const bool updateLayers = true, const bool useGPUCriteria = false )

static

Helper class to calculate cluster kinematics based on cells.

Author

Walter Lampl Walte.nosp@m.r.La.nosp@m.mpl@c.nosp@m.ern..nosp@m.ch Calculate cluster kinematics from contained cells

Definition at line 223 of file CaloClusterKineHelper.cxx.

                                                                                                                                        {
  //static CaloPhiRange range;
 
  // update global kinematics 
  //
  // for the update of the position the normalization is not a trival
  // thing. The previous implementation used the sum of weighted
  // energies as normalization. This leads to unphysical eta and phi
  // values in case negative energies are added in. The new algorithm
  // takes therefore |E| instead of E which gives the same eta and phi
  // as before for the 2 cases where all cells are negative or all
  // cells are positive. In the mixed case it will give the direction
  // of activity in the calorimeter.
 
  const CaloClusterCellLink* cccl=clu->getCellLinks();
  if (!cccl || cccl->size()==0) return; //Do nothing if no cells
 
  double posNorm = 0;
  int nBarrel = 0;
  int nEndcap = 0;
  double timeNorm = 0.;
 
  std::array<double,CaloSampling::Unknown> posSamNorm{};
 
  CellAccum accum (posNorm, posSamNorm, nBarrel, nEndcap, timeNorm, useGPUCriteria);
  //accum.theNewPhi = clu->phi(); //????
  if (useweight)
    accumCells (cccl, accum, AccumWeight(), useGPUCriteria);
  else
    accumCells (cccl, accum, AccumNoWeight(), useGPUCriteria);
  
  if ( posNorm != 0. ) {
    double inorm = 1 / posNorm;
    clu->setEta(accum.theNewEta * inorm);
    clu->setPhi(CaloPhiRange::fix(accum.theNewPhi * inorm));
  }
  else {
    clu->setEta(0);
    clu->setPhi(0);
  }
 
  clu->setE(accum.theNewEnergy);
  clu->setM(0.0);
 
  if ( timeNorm != 0. ) {
    clu->setTime(accum.theNewTime/timeNorm);
    // note that standard deviation of cell time distribution is >= 0. 
    // (no particular accommodation for <x^2> < <x><x>!) 
    //    clu->setSecondTime(std::sqrt(std::max(accum.theNewSecondTime/timeNorm-clu->time()*clu->time(),0.)));
    // consistency with other second moments: store variance!
    clu->setSecondTime(accum.theNewSecondTime/timeNorm-clu->time()*clu->time());
  } else { 
    clu->setTime(0.);      
    clu->setSecondTime(0.);
  }
 
 
  //Set Sampling pattern:
  uint32_t samplingPattern=0;
  for(int i=0;i<(int)CaloSampling::Unknown;i++) {
    if (accum.PresenceInSample[i]) 
      samplingPattern |= (0x1U<<i);
  }
  if (samplingPattern!=clu->samplingPattern()) { 
    //Clear sampling data (if there is any): Invalidated if the sampling pattern has changed 
    clu->clearSamplingData();
    clu->setSamplingPattern(samplingPattern);
  }

  // Check Sampling Variable Updates //
 
  if ( !updateLayers ) return;
 
  //First loop over samplings
  for(int i=0;i<(int)CaloSampling::Unknown;i++) {
    if ( !accum.PresenceInSample[i] ) continue; // check sampling bit
    if ( posSamNorm[i] != 0 ) {
      const double inorm = 1 / posSamNorm[i];
      accum.EtaInSample[i] *= inorm;
      accum.PhiInSample[i] = CaloPhiRange::fix (accum.PhiInSample[i] * inorm);
    }
  }//End loop over sampling
 
  //Second loop over samplings 
  //Set(Engery/Eta/Phi) methods only allowed once m_samplingPattern is final!
  for(int i=0;i<(int)CaloSampling::Unknown;i++) {
    if ( !accum.PresenceInSample[i] ) continue;
    // check sampling bit
    CaloSampling::CaloSample sam = (CaloSampling::CaloSample)i;
    clu->setEnergy(sam,accum.EnergyInSample[i]);
    // clu->setNCellsPerSampling(sam, accum.NCellsInSample[i]);
    clu->setEta(sam,accum.EtaInSample[i]);
    clu->setPhi(sam,accum.PhiInSample[i]);
 
    clu->setEmax(sam,accum.MaxEnergyInSample[i]);
    clu->setEtamax(sam,accum.EtaMaxEnergyInSample[i]);
    clu->setPhimax(sam,accum.PhiMaxEnergyInSample[i]);
  }
 
  }

---

The documentation for this class was generated from the following files:

• CaloClusterKineHelper.h
• CaloClusterKineHelper.cxx