BEmcCluster.h

Go to the documentation of this file. Or view the newest version in sPHENIX GitHub for file BEmcCluster.h

#ifndef CALORECO_EMCCLUSTER_H
#define CALORECO_EMCCLUSTER_H

// Name: EmcCluster.h
// Author: A. Bazilevsky (RIKEN-BNL)
// Major modifications by M. Volkov (RRC KI) Jan 27 2000

#include <TObject.h>

#include <cmath>
#include <cstdlib>
#include <vector>

// Forward declarations
class BEmcRec;

class EmcModule
{
 public:
  EmcModule();
  EmcModule(int ich_, float amp_, float tof_);

  virtual ~EmcModule() {}

  int ich;    // module id (linear)
  float amp;  // module signal
  float tof;  // module time-of-flight
};

// ///////////////////////////////////////////////////////////////////////////

class EmcCluster : public TObject
{
 public:
  EmcCluster()
    : fOwner(nullptr)
  {
  }

  EmcCluster(BEmcRec* sector)
    : fOwner(sector)
  {
  }


  EmcCluster(const std::vector<EmcModule>& hlist,
             BEmcRec* sector)
    : fHitList(hlist)
    , fOwner(sector)
  {
  }

  ~EmcCluster() override
  {
  }


  void ReInitialize(const std::vector<EmcModule>& hlist)
  {
    fHitList = hlist;
  }
  int GetNofHits() { return fHitList.size(); }
  std::vector<EmcModule> GetHitList() { return fHitList; };
  EmcModule GetMaxTower();
  //  EmcModule GetImpactTower();
  float GetTowerEnergy(int ich);
  float GetTowerEnergy(int ix, int iy);
  float GetTowerToF(int ich);
  float GetE4();
  float GetE9();
  float GetE9(int ich);
  float GetECore();
  float GetECoreCorrected();
  float GetTotalEnergy();
  void GetMoments(float& pxcg, float& pycg,
                  float& pxx, float& pxy, float& pyy);
  void GetCorrPos(float& xc, float& yc);
  void GetGlobalPos(float& xg, float& yg, float& zg);
  int GetSubClusters(std::vector<EmcCluster>* sClList, std::vector<EmcModule>* ppeaks);
  float GetProb(float& chi2, int& ndf);

 protected:
  std::vector<EmcModule> fHitList;

  BEmcRec* fOwner;  // what sector it belongs to

  // static members
  static int const fgMaxNofPeaks;
  static int const fgPeakIter;
  static float const fgEmin;

 public:
  // MV 2002/02/28 moved these functions here from #define's

  static int max(int a, int b)
  {
    return a > b ? a : b;
  }
  static float max(float a, float b)
  {
    return a > b ? a : b;
  }
  static double max(double a, double b)
  {
    return a > b ? a : b;
  }

  static int min(int a, int b)
  {
    return a < b ? a : b;
  }
  static float min(float a, float b)
  {
    return a < b ? a : b;
  }
  static double min(double a, double b)
  {
    return a < b ? a : b;
  }

  static int ABS(int x)
  {
    return abs(x);
  }
  static float ABS(float x)
  {
    return fabsf(x);
  }
  static double ABS(double x)
  {
    return fabs(x);
  }

  static int lowint(float x)
  {
    return x < 0. ? int(x - 1) : int(x);
  }
};

// ///////////////////////////////////////////////////////////////////////////

#endif