CaloClusterCorr::CaloSwCalibHitsShowerDepth Class Reference

Helper to calculate the shower depth as used in the calib hits SW cluster correction. More...

#include <CaloSwCalibHitsShowerDepth.h>

Collaboration diagram for CaloClusterCorr::CaloSwCalibHitsShowerDepth:

Static Public Member Functions
static double	depth (const float aeta, const float start_crack, const float end_crack, const CaloRec::Array< 2 > &sampling_depth, const float etamax, const xAOD::CaloCluster *cluster, MsgStream &log)
	Calculate the depth of the cluster. More...

Static Private Member Functions
static void	barrelCoefs (const float aeta, float R[4])
	Calculate the sampling depth coefficients for the barrel. More...
static bool	endcapCoefs (const float aeta, const CaloRec::Array< 2 > &sampling_depth, const float etamax, float R[4])
	Calculate the sampling depth coefficients for the endcap. More...

Detailed Description

Helper to calculate the shower depth as used in the calib hits SW cluster correction.

Broken out into a separate class so that it can be used by other components.

Definition at line 32 of file CaloSwCalibHitsShowerDepth.h.

Member Function Documentation

barrelCoefs()

void CaloClusterCorr::CaloSwCalibHitsShowerDepth::barrelCoefs ( const float  aeta, float  R[4]  )

staticprivate

Calculate the sampling depth coefficients for the barrel.

Parameters

aeta	abs(eta) of the cluster.
R[out]	The set of coefficients per layer.

Definition at line 94 of file CaloSwCalibHitsShowerDepth.cxx.

{
//----------------------------------------------------------------------
//     APPROXIMATE Longitudinal segmentation of the EM Barrel
//     Inspired from M.Seman's former ACCOMPA.
//     - PURPOSE Computes the shower longitudinal barycenter
//     - Inputs: shower eta and energies in the compartments
//----------------------------------------------------------------------
 
   static const float
     X0IN = 1.6,          // 1.7 X0 before strips
     X0S1 = 6.0,          //  6  X0 till end of strips
  //   X0S2 = 24.0,         // 24  X0 till end of S2
  //   X0S3 = 2.0,          //  2  X0 at least in S3
  //   RIN  = 150.0024,     // where strips start (cm)
  //   ROUT = 197.0482,     // where S3 ends (cm)
  //   RCLR = 0.7,          // clearance for resistive ink (cm)
     //BEND1 = 152.1000,    // first 3 bends radii (cm)
     //BEND2 = 155.9659,
     //BEND3 = 159.7202,
     ETASTEP = 0.8,       // where lead changes
     //DENSHI = 2.137,      // cm/X0 high density
     //DENSLO = 2.695,      // cm/X0  low density
     X0THI = 23.74,       // total depth high density part (not indep?)
     X0TLO = 19.177;      // total depth  low density part (not indep?)
 
   float X0T;
   float R12, R23;
 
   if ( aeta < ETASTEP ) {
     //DENS = DENSHI;
     X0T = X0THI;
     //RLO = BEND2;             // break strips between bends 2 and 3
     //RHI = BEND3;
    }else{
     //DENS = DENSLO;
     X0T = X0TLO;
     //RLO = BEND1;             // break strips between bends 1 and 2
     //RHI = BEND2;
   }
 
   R[0] = (X0IN) * cosh(aeta);
 
   //     R12 is radius at transition strips->Sampling2     
   //     R23 is radius at transition Sampling2->Sampling3  
    
   R12 = X0S1;
 
   //..   total depth of calo at given eta
 
   float RADLTOT = X0T*cosh(aeta);
 
   //..   radius where S2 ends. Make sure there are 2X0 left for S3
 
   if (aeta < 0.5){
     R23 = 22 * cosh(aeta) ;
   }else{
     R23 = 25;
   }
 
   R[1] = R[0]  + (R12 - R[0])/2 ;
   R[2] = R12 + (R23 - R12)/2;
   R[3] = R23 + (RADLTOT - R23)/2;
}

depth()

double CaloClusterCorr::CaloSwCalibHitsShowerDepth::depth ( const float  aeta, const float  start_crack, const float  end_crack, const CaloRec::Array< 2 > &  sampling_depth, const float  etamax, const xAOD::CaloCluster *  cluster, MsgStream &  log  )

static

Calculate the depth of the cluster.

Parameters

aeta	abs(eta) of the cluster.
start_crack	Eta of the start of the crack.
end_crack	Eta of the end of the crack.
sampling_depth	Array of sampling depths per bin/sampling in the EC.
etamax	Maximum eta value in sampling_depth.
cluster	Cluster for which to calculate the depth.
log	Stream for debug messages.

Note that the sampling depth is used only in the endcap. Parameters for the barrel are hardcoded. (FIXME!) Parameters are passed by reference to allow them to be changed.

Definition at line 36 of file CaloSwCalibHitsShowerDepth.cxx.

{
  float R[4];
  int si;
 
  if (aeta < start_crack) {
    barrelCoefs (aeta, R);
    si = 0;
  }
  else if (aeta > end_crack) {
    if (!endcapCoefs (aeta, sampling_depth, etamax, R))
      return 0;
    si = 1;
  }
  else
    return 0;
 
  static const CaloSampling::CaloSample samps[2][4] = {
    { CaloSampling::PreSamplerB,
      CaloSampling::EMB1,
      CaloSampling::EMB2,
      CaloSampling::EMB3 },
    { CaloSampling::PreSamplerE,
      CaloSampling::EME1,
      CaloSampling::EME2,
      CaloSampling::EME3 }
  };
 
  double rsum = 0;
  double esum = 0;
  for (int i=0; i < 4; i++) {
    float e = cluster->eSample (samps[si][i]);
    rsum += R[i] * e;
    esum += e;
  }
  
  if (esum == 0)
    return 0;
 
  double depth = rsum / esum;
 
  if (log.level() <= MSG::DEBUG) log << MSG::DEBUG << "Shower barycenter ---->> " << depth << endmsg;
 
  return depth;
}

endcapCoefs()

bool CaloClusterCorr::CaloSwCalibHitsShowerDepth::endcapCoefs ( const float  aeta, const CaloRec::Array< 2 > &  sampling_depth, const float  etamax, float  R[4]  )

staticprivate

Calculate the sampling depth coefficients for the endcap.

Parameters

aeta	abs(eta) of the cluster.
sampling_depth	Array of sampling depths per bin/sampling in the EC.
etamax	Maximum eta value in sampling_depth.
R[out]	The set of coefficients per layer.
aeta	abs(eta) of the cluster.
etamax	Maximum eta value in sampling_depth.
etamax	Maximum eta value in sampling_depth.
R[out]	The set of coefficients per layer.

Definition at line 168 of file CaloSwCalibHitsShowerDepth.cxx.

{
  unsigned int ibin = (static_cast<unsigned int> (aeta / etamax * 100)) ;
  if (ibin >= sampling_depth.size())
    return false;
 
  for (int i=0; i < 4; i++)
    R[i] = sampling_depth[ibin][i+1];
  return true;
}

---

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

• CaloSwCalibHitsShowerDepth.h
• CaloSwCalibHitsShowerDepth.cxx