BSEmcClusterRadiusAlgo.h

Go to the documentation of this file.

#ifndef BSEMCCLUSTERRADIUSALGO_H
#define BSEMCCLUSTERRADIUSALGO_H

#include "RtypesCore.h"

class BSEmcClusterRadiusAlgo {
 public:
  BSEmcClusterRadiusAlgo(){};
  ~BSEmcClusterRadiusAlgo(){};

  Double_t GetModifiedRadius(Double_t t_clusradius, Int_t t_n = 1) const
  {
    if (t_n == 0) {
      t_n = 1;
    }
    return t_clusradius + t_n * fSingleCrystalRadius;
  }
  // Add single crystal radius to precluster radius, because distance to crystal centre is taken and not its size.
  // If this correction is not implemented, preclusters consisting of a single crystal will get a radius of 0 and will never be merged to another cluster.
  // Added multiplier n to allow the use of the diameter, or manually increase the radius of a cluster
  Double_t GetModifiedXRadius(Double_t xrad, Int_t n = 1) const { return xrad + n * fSingleCrystalRadius; }
  Double_t GetModifiedYRadius(Double_t yrad, Int_t n = 1) const { return yrad + n * fSingleCrystalRadius; }

  void SetSingleCrystalRadius(Double_t t_radius) { fSingleCrystalRadius = t_radius; }
  Double_t GetSingleCrystalRadius() const { return fSingleCrystalRadius; }

 private:
  Double_t fSingleCrystalRadius{1.84}; // cm = single crystal radius = sqrt(2)*1.3 = 1.8385, 1.3 = 0.5*crystal size
};

#endif /*BSEMCCLUSTERRADIUSALGO_H*/