[u/mrichter/AliRoot.git] / TPC / AliTPCclustererMI.h

#ifndef ALITPCCLUSTERERMI_H
#define ALITPCCLUSTERERMI_H

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

/* $Id$ */


//-------------------------------------------------------
//                       TPC clusterer
//
//   Origin: Marian Ivanov  
//-------------------------------------------------------
#include <Rtypes.h>
#include <TObject.h>
#define kMAXCLUSTER 2500

class TFile;
class AliTPCParam;
class AliTPCRecoParam;
class AliTPCclusterMI;
class AliTPCClustersRow;
class AliRawReader;
class AliSimDigits;
class TTree;
class TTreeSRedirector;
class  AliRawEventHeaderBase;

class AliTPCclustererMI : public TObject{
public:
  AliTPCclustererMI(const AliTPCParam* par, const AliTPCRecoParam * recoParam = 0);
  virtual ~AliTPCclustererMI();
  virtual void Digits2Clusters();
  virtual void Digits2Clusters(AliRawReader* rawReader);
  virtual void SetOldRCUFormat(Bool_t rcuFormat = kFALSE)
    { fIsOldRCUFormat = rcuFormat; };
  virtual void SetInput(TTree * tree);  // set input tree with digits    
  virtual void SetOutput(TTree * tree); //set output tree with 
private:
  Bool_t IsMaximum(Float_t k, Int_t max, const Float_t *bins) const; 
  void MakeCluster2(Int_t k,Int_t max,Float_t *bins,UInt_t m,
   AliTPCclusterMI &c);  
  void MakeCluster(Int_t k,Int_t max,Float_t *bins,UInt_t m,
   AliTPCclusterMI &c); 
  Float_t  GetSigmaY2(Int_t iz);
  Float_t  GetSigmaZ2(Int_t iz);
  Float_t  FitMax(Float_t vmatrix[5][5], Float_t y, Float_t z, Float_t sigmay, Float_t sigmaz);
  void AddCluster(AliTPCclusterMI &c);  // add the cluster to the array
  void UnfoldCluster(Float_t * matrix[7], Float_t recmatrix[5][5], 
		     Float_t & meani, Float_t & meanj, Float_t & sum, Float_t &overlap );
  void FindClusters();
  Double_t  ProcesSignal(Float_t * signal, Int_t nchannels, Int_t id[3], Double_t &rms);
  void DumpHistos();


  Float_t * fBins;       //!digits array
  Float_t * fResBins;    //!digits array with res. after 1 finder
  Int_t fLoop;         //loop - cf in 2 loops
  Int_t fMaxBin;       //current ( for current sector)  maximal bin
  Int_t fMaxTime;      //current ( for current sector)  maximal time
  Int_t fMaxPad;       //current ( for current sector)  maximal pad
  Int_t fSector;      //!current sector
  Int_t fRow;         //!current row
  Float_t fSign;      //!current sign 
  Float_t fRx;        // current radius
  Float_t fPadWidth;  // the width of the pad
  Float_t fPadLength;  // the width of the pad
  Float_t fZWidth;     //the z bin width

  Bool_t  fPedSubtraction; // perform pedestal subtraction or not
  Bool_t  fIsOldRCUFormat; // assume old RCU raw data format
  AliRawEventHeaderBase *fEventHeader; //! event header information
  UInt_t  fTimeStamp;   // Time Stamp
  UInt_t  fEventType;   // Event Type
  TTree * fInput;   //!input  tree with digits - object not owner
  TTree * fOutput;   //!output tree with digits - object not owner
  AliTPCClustersRow * fRowCl;  //! current cluster row
  AliSimDigits * fRowDig;      //! current digits row
  const AliTPCParam * fParam;        //! tpc parameters
  Int_t fNcluster;             // number of clusters - for given row
  TObjArray * fAmplitudeHisto;          //! array of histograms of amplitudes
  TTreeSRedirector *fDebugStreamer;     //!debug streamer
  const AliTPCRecoParam  * fRecoParam;        //! reconstruction parameters
  ClassDef(AliTPCclustererMI,1)  // Time Projection Chamber digits
};

inline Bool_t AliTPCclustererMI::IsMaximum(Float_t q,Int_t max,const Float_t *bins) const {
  //is this a local maximum ?
  if (bins[-max] >= q) return kFALSE;
  if (bins[-1  ] >= q) return kFALSE; 
  if (bins[+max] > q) return kFALSE; 
  if (bins[+1  ] > q) return kFALSE; 
  if (bins[-max-1] >= q) return kFALSE;
  if (bins[+max-1] >= q) return kFALSE; 
  if (bins[+max+1] > q) return kFALSE; 
  if (bins[-max+1] >= q) return kFALSE;
  return kTRUE; 
}


//-----------------------------------------------------------------

#endif
Commit	Line	Data
1c53abe2	1	#ifndef ALITPCCLUSTERERMI_H
	2	#define ALITPCCLUSTERERMI_H
	3
722aa38c	4	/* Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved. *
	5	* See cxx source for full Copyright notice */
	6
	7	/* $Id$ */
	8
	9
	10
1c53abe2	11	//-------------------------------------------------------
	12	// TPC clusterer
	13	//
	14	// Origin: Marian Ivanov
	15	//-------------------------------------------------------
	16	#include <Rtypes.h>
	17	#include <TObject.h>
	18	#define kMAXCLUSTER 2500
	19
	20	class TFile;
	21	class AliTPCParam;
194b0609	22	class AliTPCRecoParam;
1c53abe2	23	class AliTPCclusterMI;
1c53abe2	24	class AliTPCClustersRow;
f8aae377	25	class AliRawReader;
1c53abe2	26	class AliSimDigits;
1c53abe2	27	class TTree;
65c045f0	28	class TTreeSRedirector;
ef4ad662	29	class AliRawEventHeaderBase;
1c53abe2	30
	31	class AliTPCclustererMI : public TObject{
	32	public:
194b0609	33	AliTPCclustererMI(const AliTPCParam* par, const AliTPCRecoParam * recoParam = 0);
65c045f0	34	virtual ~AliTPCclustererMI();
f8aae377	35	virtual void Digits2Clusters();
f8aae377	36	virtual void Digits2Clusters(AliRawReader* rawReader);
c8c88e5d	37	virtual void SetOldRCUFormat(Bool_t rcuFormat = kFALSE)
c8c88e5d	38	{ fIsOldRCUFormat = rcuFormat; };
1c53abe2	39	virtual void SetInput(TTree * tree); // set input tree with digits
	40	virtual void SetOutput(TTree * tree); //set output tree with
	41	private:
97f127bb	42	Bool_t IsMaximum(Float_t k, Int_t max, const Float_t *bins) const;
13116aec	43	void MakeCluster2(Int_t k,Int_t max,Float_t *bins,UInt_t m,
1c53abe2	44	AliTPCclusterMI &c);
13116aec	45	void MakeCluster(Int_t k,Int_t max,Float_t *bins,UInt_t m,
1c53abe2	46	AliTPCclusterMI &c);
	47	Float_t GetSigmaY2(Int_t iz);
	48	Float_t GetSigmaZ2(Int_t iz);
	49	Float_t FitMax(Float_t vmatrix[5][5], Float_t y, Float_t z, Float_t sigmay, Float_t sigmaz);
	50	void AddCluster(AliTPCclusterMI &c); // add the cluster to the array
13116aec	51	void UnfoldCluster(Float_t * matrix[7], Float_t recmatrix[5][5],
1c53abe2	52	Float_t & meani, Float_t & meanj, Float_t & sum, Float_t &overlap );
f8aae377	53	void FindClusters();
ef4ad662	54	Double_t ProcesSignal(Float_t * signal, Int_t nchannels, Int_t id[3], Double_t &rms);
65c045f0	55	void DumpHistos();
1c53abe2	56
1c53abe2	57
13116aec	58	Float_t * fBins; //!digits array
13116aec	59	Float_t * fResBins; //!digits array with res. after 1 finder
1c53abe2	60	Int_t fLoop; //loop - cf in 2 loops
97f127bb	61	Int_t fMaxBin; //current ( for current sector) maximal bin
	62	Int_t fMaxTime; //current ( for current sector) maximal time
	63	Int_t fMaxPad; //current ( for current sector) maximal pad
1c53abe2	64	Int_t fSector; //!current sector
508541c7	65	Int_t fRow; //!current row
1c53abe2	66	Float_t fSign; //!current sign
	67	Float_t fRx; // current radius
	68	Float_t fPadWidth; // the width of the pad
	69	Float_t fPadLength; // the width of the pad
	70	Float_t fZWidth; //the z bin width
	71
b6f060dc	72	Bool_t fPedSubtraction; // perform pedestal subtraction or not
c8c88e5d	73	Bool_t fIsOldRCUFormat; // assume old RCU raw data format
ef4ad662	74	AliRawEventHeaderBase *fEventHeader; //! event header information
	75	UInt_t fTimeStamp; // Time Stamp
	76	UInt_t fEventType; // Event Type
1c53abe2	77	TTree * fInput; //!input tree with digits - object not owner
	78	TTree * fOutput; //!output tree with digits - object not owner
	79	AliTPCClustersRow * fRowCl; //! current cluster row
	80	AliSimDigits * fRowDig; //! current digits row
	81	const AliTPCParam * fParam; //! tpc parameters
	82	Int_t fNcluster; // number of clusters - for given row
65c045f0	83	TObjArray * fAmplitudeHisto; //! array of histograms of amplitudes
65c045f0	84	TTreeSRedirector *fDebugStreamer; //!debug streamer
194b0609	85	const AliTPCRecoParam * fRecoParam; //! reconstruction parameters
1c53abe2	86	ClassDef(AliTPCclustererMI,1) // Time Projection Chamber digits
	87	};
	88
97f127bb	89	inline Bool_t AliTPCclustererMI::IsMaximum(Float_t q,Int_t max,const Float_t *bins) const {
1c53abe2	90	//is this a local maximum ?
	91	if (bins[-max] >= q) return kFALSE;
	92	if (bins[-1 ] >= q) return kFALSE;
	93	if (bins[+max] > q) return kFALSE;
	94	if (bins[+1 ] > q) return kFALSE;
	95	if (bins[-max-1] >= q) return kFALSE;
	96	if (bins[+max-1] >= q) return kFALSE;
	97	if (bins[+max+1] > q) return kFALSE;
	98	if (bins[-max+1] >= q) return kFALSE;
	99	return kTRUE;
	100	}
	101
	102
	103
	104	//-----------------------------------------------------------------
	105
	106	#endif
	107
	108