splitcalCluster Class Reference

#include <splitcalCluster.h>

Inheritance diagram for splitcalCluster:

Collaboration diagram for splitcalCluster:

Public Member Functions
	splitcalCluster ()
virtual	~splitcalCluster ()
virtual void	Print () const
void	SetEtaPhiE (double &eta, double &phi, double &e)
void	SetEta (double &eta)
void	SetPhi (double &phi)
void	SetEnergy (double &e)
void	SetIndex (int i)
void	SetStartPoint (const double &x, const double &y, const double &z)
void	SetEndPoint (const double &x, const double &y, const double &z)
void	SetHitIndices (const std::vector< int > &v)
void	AddHit (int hitIndex, double weight=1.0)
int	GetIndex ()
double	GetEta ()
double	GetPhi ()
double	GetEnergy ()
double	GetPx ()
double	GetPy ()
double	GetPz ()
double	GetEx ()
double	GetEy ()
double	GetEz ()
TVector3	GetStartPoint () const
TVector3	GetEndPoint () const
const std::vector< int > &	GetHitIndices () const
const std::vector< double > &	GetHitWeights () const
void	ComputeEtaPhiE (const std::vector< splitcalHit > &hits)
	splitcalCluster (const splitcalCluster &cluster)=default
splitcalCluster &	operator= (const splitcalCluster &cluster)=default

Private Member Functions
	ClassDef (splitcalCluster, 4)

Private Attributes
int	_index
double	_eta
double	_phi
double	_energy
std::array< double, 3 >	_start
std::array< double, 3 >	_end
std::vector< int >	_hitIndices
std::vector< double >	_hitWeights

Detailed Description

Definition at line 16 of file splitcalCluster.h.

Constructor & Destructor Documentation

splitcalCluster() [1/2]

splitcalCluster::splitcalCluster

(

)

Constructors

Definition at line 12 of file splitcalCluster.cxx.

{}

~splitcalCluster()

splitcalCluster::~splitcalCluster ( )

virtual

Destructor

Definition at line 93 of file splitcalCluster.cxx.

{}

splitcalCluster() [2/2]

splitcalCluster::splitcalCluster ( const splitcalCluster &  cluster )

default

Copy constructor

Member Function Documentation

AddHit()

void splitcalCluster::AddHit ( int  hitIndex, double  weight = 1.0  )

inline

Definition at line 47 of file splitcalCluster.h.

                                                 {
    _hitIndices.push_back(hitIndex);
    _hitWeights.push_back(weight);
  }

ClassDef()

splitcalCluster::ClassDef ( splitcalCluster  , 4    )

private

ComputeEtaPhiE()

void splitcalCluster::ComputeEtaPhiE

(

const std::vector< splitcalHit > &

hits

)

Definition at line 14 of file splitcalCluster.cxx.

                                                                       {
  // Compute energy weighted average for hits in the same layer
  // This is in preparation of the linear fit to get the cluster eta and phi
 
  // maps to compute the weighted average of all the hits in the same layer
  std::map<int, double> mapLayerWeigthedX;
  std::map<int, double> mapLayerWeigthedY;
  std::map<int, double> mapLayerZ1;
  std::map<int, double> mapLayerZ2;
  std::map<int, double> mapLayerSumWeigthsX;
  std::map<int, double> mapLayerSumWeigthsY;
 
  // loop over hits to compute cluster energy sum and to compute the coordinates
  // weighted average per layer
  double energy = 0.;
  for (size_t i = 0; i < _hitIndices.size(); ++i) {
    const auto& hit = hits[_hitIndices[i]];
    double hitEnergy =
        hit.GetEnergy() * _hitWeights[i];  // Use weight from cluster
    energy += hitEnergy;
    int layer = hit.GetLayerNumber();
    // hits from high precision layers give both x and y coordinates --> use
    // if-if instead of if-else
    if (hit.IsX()) {
      if (!mapLayerWeigthedX.contains(layer)) {
        mapLayerWeigthedX[layer] = 0.;
        mapLayerSumWeigthsX[layer] = 0.;
      }
      mapLayerWeigthedX[layer] += hit.GetX() * hitEnergy;
      mapLayerSumWeigthsX[layer] += hitEnergy;
      mapLayerZ1[layer] = hit.GetZ();
    }
    if (hit.IsY()) {
      if (!mapLayerWeigthedY.contains(layer)) {
        mapLayerWeigthedY[layer] = 0.;
        mapLayerSumWeigthsY[layer] = 0.;
      }
      mapLayerWeigthedY[layer] += hit.GetY() * hitEnergy;
      mapLayerSumWeigthsY[layer] += hitEnergy;
      mapLayerZ2[layer] = hit.GetZ();
    }
  }  // end loop on hit
 
  auto const& [minLayerX, weightedMinX] = *mapLayerWeigthedX.begin();
  double minX = weightedMinX / mapLayerSumWeigthsX[minLayerX];
  double minZ1 = mapLayerZ1[minLayerX];
 
  auto const& [minLayerY, weightedMinY] = *mapLayerWeigthedY.begin();
  double minY = weightedMinY / mapLayerSumWeigthsY[minLayerY];
  double minZ2 = mapLayerZ1[minLayerY];
 
  double minZ = (minZ1 + minZ2) / 2.;
 
  SetStartPoint(minX, minY, minZ);
 
  auto const& [maxLayerX, weightedMaxX] = *mapLayerWeigthedX.rbegin();
  double maxX = weightedMaxX / mapLayerSumWeigthsX[maxLayerX];
  double maxZ1 = mapLayerZ1[maxLayerX];
 
  auto const& [maxLayerY, weightedMaxY] = *mapLayerWeigthedY.rbegin();
  double maxY = weightedMaxY / mapLayerSumWeigthsY[maxLayerY];
  double maxZ2 = mapLayerZ1[maxLayerY];
 
  double maxZ = (maxZ1 + maxZ2) / 2.;
 
  SetEndPoint(maxX, maxY, maxZ);
 
  // get direction vector from end-start vector difference
  TVector3 start(_start[0], _start[1], _start[2]);
  TVector3 end(_end[0], _end[1], _end[2]);
  TVector3 direction = end - start;
  double eta = direction.Eta();
  double phi = direction.Phi();
  SetEtaPhiE(eta, phi, energy);
}

GetEndPoint()

TVector3 splitcalCluster::GetEndPoint ( ) const

inline

Definition at line 65 of file splitcalCluster.h.

{ return TVector3(_end[0], _end[1], _end[2]); }

GetEnergy()

double splitcalCluster::GetEnergy ( )

inline

Definition at line 55 of file splitcalCluster.h.

{ return _energy; }

GetEta()

double splitcalCluster::GetEta ( )

inline

Definition at line 53 of file splitcalCluster.h.

{ return _eta; }

GetEx()

double splitcalCluster::GetEx ( )

inline

Definition at line 59 of file splitcalCluster.h.

{ return GetPx(); }

GetEy()

double splitcalCluster::GetEy ( )

inline

Definition at line 60 of file splitcalCluster.h.

{ return GetPy(); }

GetEz()

double splitcalCluster::GetEz ( )

inline

Definition at line 61 of file splitcalCluster.h.

{ return GetPz(); }

GetHitIndices()

const std::vector< int > & splitcalCluster::GetHitIndices ( ) const

inline

Definition at line 66 of file splitcalCluster.h.

{ return _hitIndices; }

GetHitWeights()

const std::vector< double > & splitcalCluster::GetHitWeights ( ) const

inline

Definition at line 67 of file splitcalCluster.h.

{ return _hitWeights; }

GetIndex()

int splitcalCluster::GetIndex ( )

inline

Definition at line 52 of file splitcalCluster.h.

{ return _index; }

GetPhi()

double splitcalCluster::GetPhi ( )

inline

Definition at line 54 of file splitcalCluster.h.

{ return _phi; }

GetPx()

double splitcalCluster::GetPx ( )

inline

Definition at line 56 of file splitcalCluster.h.

{ return _energy * sin(_eta) * cos(_phi); }

GetPy()

double splitcalCluster::GetPy ( )

inline

Definition at line 57 of file splitcalCluster.h.

{ return _energy * sin(_eta) * sin(_phi); }

GetPz()

double splitcalCluster::GetPz ( )

inline

Definition at line 58 of file splitcalCluster.h.

{ return _energy * cos(_eta); }

GetStartPoint()

TVector3 splitcalCluster::GetStartPoint ( ) const

inline

Definition at line 62 of file splitcalCluster.h.

                                 {
    return TVector3(_start[0], _start[1], _start[2]);
  }

operator=()

splitcalCluster & splitcalCluster::operator= ( const splitcalCluster &  cluster )

default

Print()

void splitcalCluster::Print ( ) const

virtual

Methods

Definition at line 97 of file splitcalCluster.cxx.

                                  {
  std::cout << "-I- splitcalCluster: " << std::endl;
  std::cout << "    (eta,phi,energy) = " << _eta << " ,  " << _phi << " ,  "
            << _energy << std::endl;
  std::cout << "    start(x,y,z) = " << _start[0] << " ,  " << _start[1]
            << " ,  " << _start[2] << std::endl;
  std::cout << "    end(x,y,z) = " << _end[0] << " ,  " << _end[1] << " ,  "
            << _end[2] << std::endl;
  std::cout << "------- " << std::endl;
}

SetEndPoint()

void splitcalCluster::SetEndPoint ( const double &  x, const double &  y, const double &  z  )

inline

Definition at line 41 of file splitcalCluster.h.

                                                                      {
    _end[0] = x;
    _end[1] = y;
    _end[2] = z;
  }

SetEnergy()

void splitcalCluster::SetEnergy ( double &  e )

inline

Definition at line 34 of file splitcalCluster.h.

{ _energy = e; }

SetEta()

void splitcalCluster::SetEta ( double &  eta )

inline

Definition at line 32 of file splitcalCluster.h.

{ _eta = eta; }

SetEtaPhiE()

void splitcalCluster::SetEtaPhiE ( double &  eta, double &  phi, double &  e  )

inline

Definition at line 27 of file splitcalCluster.h.

                                                       {
    _eta = eta;
    _phi = phi;
    _energy = e;
  }

SetHitIndices()

void splitcalCluster::SetHitIndices ( const std::vector< int > &  v )

inline

Definition at line 46 of file splitcalCluster.h.

{ _hitIndices = v; }

SetIndex()

void splitcalCluster::SetIndex ( int  i )

inline

Definition at line 35 of file splitcalCluster.h.

{ _index = i; }

SetPhi()

void splitcalCluster::SetPhi ( double &  phi )

inline

Definition at line 33 of file splitcalCluster.h.

{ _phi = phi; }

SetStartPoint()

void splitcalCluster::SetStartPoint ( const double &  x, const double &  y, const double &  z  )

inline

Definition at line 36 of file splitcalCluster.h.

                                                                        {
    _start[0] = x;
    _start[1] = y;
    _start[2] = z;
  }

Member Data Documentation

_end

std::array<double, 3> splitcalCluster::_end

private

Definition at line 79 of file splitcalCluster.h.

_energy

double splitcalCluster::_energy

private

Definition at line 77 of file splitcalCluster.h.

_eta

double splitcalCluster::_eta

private

Definition at line 77 of file splitcalCluster.h.

_hitIndices

std::vector<int> splitcalCluster::_hitIndices

private

Definition at line 80 of file splitcalCluster.h.

_hitWeights

std::vector<double> splitcalCluster::_hitWeights

private

Definition at line 81 of file splitcalCluster.h.

_index

int splitcalCluster::_index

private

Definition at line 76 of file splitcalCluster.h.

_phi

double splitcalCluster::_phi

private

Definition at line 77 of file splitcalCluster.h.

_start

std::array<double, 3> splitcalCluster::_start

private

Definition at line 78 of file splitcalCluster.h.

---

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

• splitcal/splitcalCluster.h
• splitcal/splitcalCluster.cxx