cluster.h

Go to the documentation of this file.

#ifndef CLUSTER_H
#define CLUSTER_H

#include <list.h>

struct hit { int ich; float amp; float tof; };

class emshower;
class peakarea;

class cluster {
protected:
     list<hit> HitList;
public:
     cluster() {};
     cluster( list<hit> hlist ) { HitList = hlist; }
     void ReInitialize( list<hit> hlist ) { HitList = hlist; }
     int GetNofHits() { return HitList.size(); }
     void GetHits(hit* phit, int n);
     void GetHitList(list<hit>* plist) { *plist = HitList; };
     hit GetMaxTower();
     hit GetImpactTower();
     float GetTowerEnergy( int ich );
     float GetTowerToF( int ich );
     float GetE9( int ich );
     float GetTotalEnergy();
     void GetMoments( float* pxcg, float* pycg, 
                      float* pxx, float* pxy, float* pyy );
     void GetCorrPos( float* pxc, float* pyc );
     void GetGlobalPos( float* pxg, float* pyg, float* pzg );
     void GetErrors( float* pde, float* pdx, float* pdy, float* pdz);
     void GetChar( float* pe, 
                   float* pxcg, float* pysg, 
                   float* pxc, float* pyc, 
                   float* pxg, float* pyg, float* pzg,
                   float* pxx, float* pxy, float* pyy,
                   float* pde, float* pdx, float* pdy, float* pdz );
     int GetPeaks( peakarea*, hit* );
};

class peakarea: public cluster{
public:
     peakarea(){}
     peakarea( list<hit> hlist ):cluster(hlist) {}
     float GetChi2();
     void GetCGmin( float* pxcgmin, float* pycgmin );
     void GetChar( float* pe, 
                   float* pxcg, float* pysg, 
                   float* pxcgmin, float* pysgmin, 
                   float* pxc, float* pyc, 
                   float* pxg, float* pyg, float* pzg,
                   float* pxx, float* pxy, float* pyy, 
                   float* pchi,
                   float* pde, float* pdx, float* pdy, float* pdz );
     int GetGammas( emshower* );
};

class emshower { 
     float energy;
     float xcg;
     float ycg;
     float chisq;
public:
     emshower() { energy=0; xcg=0; ycg=0; }
     emshower( float e, float x, float y, float chi ) { energy=e; xcg=x; ycg=y; chisq=chi; }
     void ReInitialize( float e, float x, float y, float chi ) { energy=e; xcg=x; ycg=y; chisq=chi; }
     float GetTotalEnergy() { return energy; }
     void GetCG( float* px, float* py ) { *px=xcg; *py=ycg; }
     void GetCorrPos( float*, float* );
     void GetGlobalPos( float*, float*, float* );
     float GetChi2() { return chisq; }
     void GetErrors( float* pde, float* pdx, float* pdy, float* pdz);
     void GetChar( float* pe, 
                   float* pxcg, float* pysg, 
                   float* pxc, float* pyc, 
                   float* pxg, float* pyg, float* pzg,
                   float* pchi,
                   float* pde, float* pdx, float* pdy, float* pdz );
};

struct SecGeom {
    short nz;
    short ny;
    float Tower_zSize; // Tower size in Z dir
    float Tower_ySize; // Tower size in Y dir
    float xyz0[3];    // Coordinate of the 0th tower
    float nxyz[3];    // Normal unit vector inside the tower
    short zFlag;      // if zFlag=1, tower counting in z-dir is reversed
};

extern float EMCcell(float, float, float);
extern int Find_Clusters( list<hit>, list<cluster>* );
extern void SetGeometry( SecGeom, float* );
extern void SetThreshold( float );
extern void SetChi2Limit( int );
extern float Chi2Correct( float Chi2, int ND );

#endif