TPC/AliTPCclustererKr.h

   1 #ifndef ALITPCCLUSTERERKR_H\r
   2 #define ALITPCCLUSTERERKR_H\r
   3 /* Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved. *\r
   4  * See cxx source for full Copyright notice                               */\r
   5 \r
   6 /* $Id: AliTPCclustererKr.h,v 1.8 2008/02/07 16:07:15 matyja Exp $ */\r
   7 \r
   8 //-------------------------------------------------------\r
   9 //                    TPC Kr Clusterer Class\r
  10 //\r
  11 //   Origin: Adam Matyja, INP PAN, adam.matyja@ifj.edu.pl\r
  12 //-------------------------------------------------------\r
  13 \r
  14 #include "TObject.h"\r
  15 \r
  16 class AliTPCclusterKr;\r
  17 class AliPadMax;\r
  18 class AliSimDigits;\r
  19 class AliTPCv4;\r
  20 class AliTPCParam;\r
  21 class AliTPCDigitsArray;\r
  22 class AliTPCvtpr;\r
  23 class AliTPCClustersRow;\r
  24 class TTree;\r
  25 class TH1F;\r
  26 class TH2F;\r
  27 \r
  28 class AliTPCTransform;\r
  29 \r
  30 //used in raw data finder\r
  31 class AliTPCROC;\r
  32 class AliTPCCalPad;\r
  33 class AliTPCAltroMapping;\r
  34 class AliTPCcalibDB;\r
  35 class AliTPCRawStream;\r
  36 class AliTPCRecoParam;\r
  37 class AliTPCReconstructor;\r
  38 class AliRawReader;\r
  39 class AliTPCCalROC;\r
  40 class TTreeSRedirector;\r
  41 class AliRawEventHeaderBase;\r
  42 //_____________________________________________________________________________\r
  43 class AliTPCclustererKr: public TObject{\r
  44 public:\r
  45   AliTPCclustererKr();\r
  46   AliTPCclustererKr(const AliTPCclustererKr &param);//copy constructor\r
  47   AliTPCclustererKr &operator = (const AliTPCclustererKr & param); \r
  48   virtual ~AliTPCclustererKr();\r
  49 \r
  50   //finders\r
  51   virtual Int_t FinderIO();//for MC\r
  52   virtual Int_t FinderIO(AliRawReader* rawReader);//for data\r
  53   virtual Int_t FinderIOold(AliRawReader* rawReader);//for data\r
  54   virtual Int_t FindClusterKrIO();//main routine for finding clusters\r
  55   virtual void CleanSector(Int_t sector); // clean isolated digits\r
  56 \r
  57   //other\r
  58   void GetXY(Int_t sec,Int_t row,Int_t pad,Double_t& xGlob,Double_t& yGlob);//give XY coordinate of the pad\r
  59 \r
  60   virtual void SetInput(TTree * tree);  //set input tree with digits    \r
  61   virtual void SetOutput(TTree * tree); //set output tree with clusters\r
  62   virtual void SetParam(AliTPCParam *param){fParam=param;}//set TPC parameters\r
  63   virtual void SetDigArr(AliTPCDigitsArray *digarr){fDigarr=digarr;}//set current array of digits\r
  64   virtual void SetRecoParam(AliTPCRecoParam *recoParam=0);//set reconstruction parameters\r
  65 \r
  66   virtual void SetTimeStamp(UInt_t timestamp){ fTimeStamp = timestamp; }\r
  67   virtual void SetRun(UInt_t run){ fRun = run; }\r
  68 \r
  69   //setters for cluster finder parameters\r
  70   virtual void SetZeroSup(Int_t v){fZeroSup=v;}//set zero suppresion parameter\r
  71   virtual void SetFirstBin(Int_t v){fFirstBin=v;}//set first considered timebin\r
  72   virtual void SetLastBin(Int_t v){fLastBin=v;}//set last considered timebin\r
  73   virtual void SetMaxNoiseAbs(Float_t v){fMaxNoiseAbs=v;}//set maximal noise value\r
  74   virtual void SetMaxNoiseSigma(Float_t v){fMaxNoiseSigma=v;}//set maximal noise sigma\r
  75 \r
  76   virtual void SetMinAdc(Int_t v){v<=0?fMinAdc=1:fMinAdc=v;}//set fMinAdc\r
  77   virtual void SetMinTimeBins(Int_t v){fMinTimeBins=v;}//set fMinTimeBins\r
  78 //  virtual void SetMaxPadRange(Int_t v){fMaxPadRange=v;}//set fMaxPadRange\r
  79 //  virtual void SetMaxRowRange(Int_t v){fMaxRowRange=v;}//set fMaxRowRange\r
  80   virtual void SetMaxTimeRange(Int_t v){fMaxTimeRange=v;}//set fMaxTimeRange\r
  81   virtual void SetValueToSize(Float_t v){fValueToSize=v;}//set fValueToSize\r
  82 \r
  83   virtual void SetMaxPadRangeCm(Double_t v){fMaxPadRangeCm=v;}//set fMaxPadRangeCm\r
  84   virtual void SetMaxRowRangeCm(Double_t v){fMaxRowRangeCm=v;}//set fMaxRowRangeCm\r
  85 \r
  86   virtual void SetIsolCut(Short_t v){fIsolCut=v;}\r
  87 \r
  88   virtual void SetDebugLevel(Int_t debug){fDebugLevel=debug;}\r
  89   //debug = 0 to 71 -sector number to  print\r
  90   // = 72 - all sectors\r
  91   // = 73 - inners\r
  92   // = 74 - outers\r
  93 \r
  94   virtual void SetHistoRow(TH1F *histo)   {fHistoRow   =histo;}\r
  95   virtual void SetHistoPad(TH1F *histo)   {fHistoPad   =histo;}\r
  96   virtual void SetHistoTime(TH1F *histo)  {fHistoTime  =histo;}\r
  97   virtual void SetHistoRowPad(TH2F *histo){fHistoRowPad=histo;}\r
  98 \r
  99   //getters for cluster finder parameters\r
 100   Int_t GetZeroSup() const {return fZeroSup;}//get zero suppresion parameter\r
 101   Int_t GetFirstBin() const {return fFirstBin;}//get first considered timebin\r
 102   Int_t GetLastBin() const {return fLastBin;}//get last considered timebin\r
 103   Float_t GetMaxNoiseAbs() const {return fMaxNoiseAbs;}//get maximal noise value\r
 104   Float_t GetMaxNoiseSigma() const {return fMaxNoiseSigma;}//get maximal noise sigma\r
 105 \r
 106   Int_t GetMinAdc() const {return fMinAdc;}//get fMinAdc\r
 107   Int_t GetMinTimeBins() const {return fMinTimeBins;}//get fMinTimeBins\r
 108 //  Int_t GetMaxPadRange() const {return fMaxPadRange;}//get fMaxPadRange\r
 109 //  Int_t GetMaxRowRange() const {return fMaxRowRange;}//get fMaxRowRange\r
 110   Int_t GetMaxTimeRange() const {return fMaxTimeRange;}//get fMaxTimeRange\r
 111   Float_t GetValueToSize() const {return fValueToSize;}//get fValueToSize\r
 112 \r
 113   Double_t GetMaxPadRangeCm() const {return fMaxPadRangeCm;}//get fMaxPadRangeCm\r
 114   Double_t GetMaxRowRangeCm() const {return fMaxRowRangeCm;}//get fMaxRowRangeCm\r
 115 \r
 116   Short_t GetIsolCut() const {return fIsolCut;}\r
 117 \r
 118   Int_t GetDebugLevel() const {return fDebugLevel;}\r
 119   TH1F * GetHistoRow(){return fHistoRow;} \r
 120   TH1F * GetHistoPad(){return fHistoPad;}\r
 121   TH1F * GetHistoTime(){return fHistoTime;}\r
 122   TH2F * GetHistoRowPad(){return fHistoRowPad;}\r
 123 \r
 124   UInt_t GetTimeStamp() const {return fTimeStamp;}\r
 125   UInt_t GetRun() const {return fRun;}\r
 126 \r
 127 private:\r
 128   void MakeClusters(TObjArray * maximaInSector, Int_t iSec, Int_t &clusterCounter);\r
 129   Bool_t fRawData; //flag =0 for MC =1 for real data\r
 130 \r
 131   TTree * fInput;   //!input  tree with digits - object not owner\r
 132   TTreeSRedirector * fOutput;   //!output tree with clusters - object not owner\r
 133 \r
 134   AliTPCParam * fParam;//!TPC parameters\r
 135   AliTPCDigitsArray *fDigarr;//! pointer to current array if digits\r
 136 \r
 137   //only for raw data :)\r
 138   const AliTPCRecoParam  * fRecoParam;        //! reconstruction parameters\r
 139 \r
 140   //cluster finder parameters\r
 141   Int_t fZeroSup;//zero suppresion parameter = 2 def.\r
 142   Int_t fFirstBin;//first considered time bin used by cluster finder = 60 def.\r
 143   Int_t fLastBin;//last considered time bin used by cluster finder = 950 def.\r
 144   Float_t fMaxNoiseAbs;// maximal noise value on pad used in cluster finder = 2 def.\r
 145   Float_t fMaxNoiseSigma;// maximal noise sigma on pad used in cluster finder = 3 def.\r
 146 \r
 147   Int_t fMinAdc;//minimal value of acd count in each timebin = 3 def.\r
 148   Int_t fMinTimeBins;//minimal value of time bins one after each other = 2 def.\r
 149 //  Int_t fMaxPadRange;//maximal pad range from maximum = 4 def.\r
 150 //  Int_t fMaxRowRange;//maximal row range from maximum = 3 def.\r
 151   Int_t fMaxTimeRange;//maximal time bin range from maximum = 7 def.\r
 152   Float_t fValueToSize;//ratio cluster value to cluster size = 3.5 def.\r
 153 \r
 154   Double_t fMaxPadRangeCm;//maximal pad range in cm from maximum = 2.5cm def.\r
 155   Double_t fMaxRowRangeCm;//maximal row range in cm from maximum = 3.5cm def.\r
 156 \r
 157   Short_t fIsolCut;//isolation cut in 3D = 5 def.\r
 158 \r
 159   Int_t fDebugLevel;//! debug level variable\r
 160   TH1F *fHistoRow;//!debug histo for rows\r
 161   TH1F *fHistoPad;//!debug histo for pads\r
 162   TH1F *fHistoTime;//!debug histo for timebins\r
 163   TH2F *fHistoRowPad;//!debug histo for rows and pads\r
 164 \r
 165   UInt_t fTimeStamp; //!time stamp from event header\r
 166   UInt_t fRun; //!run from event header\r
 167   ClassDef(AliTPCclustererKr,8)  // Time Projection Chamber Kr clusters\r
 168 };\r
 169 \r
 170 \r
 171 #endif\r