Adding interaction time, time and amplitude for each PMT (Alla)

[u/mrichter/AliRoot.git] / RICH / AliRICHCluster.h

#ifndef AliRICHCluster_h
#define AliRICHCluster_h
/* Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved. *
 * See cxx source for full Copyright notice                               */

#include "AliRICHParam.h"    //DigAdd(), Dig()

class AliRICHCluster :public TObject
{
public:
  enum EClusterStatus {kFormed,kCoG,kUnfolded,kEmpty=-1}; 
                       AliRICHCluster(                                     ):TObject(),fQdc(-1),fCham(-1),fX(-1),fY(-1),fStatus(kEmpty   ),fDigs(0) {}
                       AliRICHCluster(Int_t c,Double_t x,Double_t y,Int_t q):TObject(),fQdc(q ),fCham(c) ,fX(x ),fY(y ),fStatus(kUnfolded),fDigs(0) {}
  virtual             ~AliRICHCluster(                                     )                                                                        {DigDel();}
//framework part                   
         void          Print  (Option_t *opt=""                                  )const;                  //overloaded TObject::Print() to print cluster info
  static void          FitFunc(Int_t &, Double_t *, Double_t &, Double_t *, Int_t);                       //fit function to be used by MINUIT
//private part  
         void          CoG      (                                         );                                                      //calculates center of gravity
         Int_t         C        (                                         )const{return fCham;                                  } //chamber number
  inline void          DigAdd   (AliRICHDigit *pDig                       );                                                      //add new digit ot the cluster
         void          DigDel   (                                         )     {if(fDigs) {delete fDigs;fDigs=0;}              } //deletes the list of digits (not digits!) 
         void          DistXY   (const TVector2 &p,Double_t &x,Double_t &y)const{ x=p.X()-fX; y=p.Y()-fY;                       } //distance in x to given point 
         AliRICHDigit* Dig      (Int_t i                                  )const{return (AliRICHDigit*)fDigs->At(i);            } //pointer to i-th digit
         TObjArray*    Digits   (                                         )const{return fDigs;                                  } //list of digits  
         TVector3      Lors2Mars()                            const{return AliRICHParam::Instance()->Lors2Mars(fCham,fX,fY);     } //cluster position in MARS
         void          Reset    (                                         )     {DigDel();fQdc=fCham=-1;fX=fY=-1;fStatus=kEmpty;} //cleans the cluster
         Int_t         Solve    (TClonesArray *pCluLst,Bool_t isUnfold    );                                                      //solve cluster: MINUIT fit or CoG
         Int_t         Size     (                                         )const{return (fDigs)?fDigs->GetEntriesFast():0;      } //number of pads in cluster
         Int_t         Q        (                                         )const{return fQdc;                                   } //cluster charge in QDC channels 
         Double_t      X        (                                         )const{return fX;                                     } //cluster x position in LRS
         Double_t      Y        (                                         )const{return fY;                                     } //cluster y position in LRS 
//test part         
  static void          Test     (Double_t x,Double_t y,Double_t e=0,Bool_t isUnfold=kTRUE);                                        //test by hit (x [cm] , y [cm] , e [GeV])
  static void          Test     (                                                        );                                        //test by predifined patterns
protected:
  Int_t         fQdc;         //QDC value
  Int_t         fCham;        //10*chamber number+sector number 
  Double_t      fX;           //local x postion, cm 
  Double_t      fY;           //local y postion, cm  
  Int_t         fStatus;      //flag to mark the quality of the cluster   
  TObjArray    *fDigs;        //! list of digits forming this cluster
  ClassDef(AliRICHCluster,3)  //RICH cluster class       
};//class AliRICHCluster
//++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
void AliRICHCluster::DigAdd(AliRICHDigit *pDig)
{
// Adds a given digit to the list of digits belonging to this cluster    
  if(!fDigs) {fQdc=0;fDigs = new TObjArray;}
  fDigs->Add(pDig);
  fQdc+=(Int_t)pDig->Qdc(); 
  fCham=pDig->C();
  fStatus=kFormed;
}
//++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
#endif