[u/mrichter/AliRoot.git] / TRD / Cal / AliTRDCalROC.h

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

/* $Id$ */

//////////////////////////////////////////////////
//                                              //
//  TRD calibration base class for one ROC      //
//                                              //
//////////////////////////////////////////////////

#include <TObject.h>

class TArrayI;
class TArrayF;
class TH2F;
class TH1F;
class TGraph2D;
class TH2D;

//_____________________________________________________________________________
class AliTRDCalROC : public TObject 
{

 public:

  AliTRDCalROC();
  AliTRDCalROC(Int_t p, Int_t c);
  AliTRDCalROC(const AliTRDCalROC &c);
  virtual      ~AliTRDCalROC();
  AliTRDCalROC &operator=(const AliTRDCalROC &c);
  virtual void  Copy(TObject &c) const;

  Int_t         GetNrows() const                   { return fNrows; };
  Int_t         GetNcols() const                   { return fNcols; };

  Int_t         GetChannel(Int_t c, Int_t r) const { return r+c*fNrows; };
  Int_t         GetNchannels() const               { return fNchannels; };

  Float_t       GetValue(Int_t ich) const          { return (Float_t) fData[ich] / 10000;  };
  Float_t       GetValue(Int_t col, Int_t row)     { return GetValue(GetChannel(col,row)); };

  void          SetValue(Int_t ich, Float_t value) { fData[ich] = (UShort_t) (value * 10000); };
  void          SetValue(Int_t col, Int_t row, Float_t value) 
                                                   { SetValue(GetChannel(col,row), value);    };

  // statistic
  //
  Double_t GetMean(AliTRDCalROC * const outlierROC=0) const; 
  Double_t GetRMS(AliTRDCalROC * const outlierROC=0) const;
  Double_t GetMedian(AliTRDCalROC * const outlierROC=0) const;
  Double_t GetLTM(Double_t *sigma=0, Double_t fraction=0.9, AliTRDCalROC * const outlierROC=0);

  // algebra
  Bool_t Add(Float_t c1);
  Bool_t Multiply(Float_t c1);
  Bool_t Add(const AliTRDCalROC * roc, Double_t c1 = 1);
  Bool_t Multiply(const AliTRDCalROC * roc);   
  Bool_t Divide(const AliTRDCalROC * roc);   

  // noise
  Bool_t Unfold();
  
  //Plots
  TH2F *   MakeHisto2D(Float_t min, Float_t max,Int_t type, Float_t mu = 1.0);
  TH1F *   MakeHisto1D(Float_t min, Float_t max,Int_t type, Float_t mu = 1.0);
  
 protected:

  Int_t     fPla;              //  Plane number
  Int_t     fCha;              //  Chamber number

  Int_t     fNrows;            //  Number of rows
  Int_t     fNcols;            //  Number of columns

  Int_t     fNchannels;        //  Number of channels
  UShort_t *fData;             //[fNchannels] Data

  ClassDef(AliTRDCalROC, 2)    //  TRD ROC calibration class

};

#endif
Commit	Line	Data
	1	#ifndef ALITRDCALROC_H
	2	#define ALITRDCALROC_H
	3	/* Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved. *
	4	* See cxx source for full Copyright notice */
	5
	6	/* $Id$ */
	7
	8	//////////////////////////////////////////////////
	9	// //
	10	// TRD calibration base class for one ROC //
	11	// //
	12	//////////////////////////////////////////////////
	13
	14	#include <TObject.h>
	15
	16	class TArrayI;
	17	class TArrayF;
	18	class TH2F;
	19	class TH1F;
	20	class TGraph2D;
	21	class TH2D;
	22
	23	//_____________________________________________________________________________
	24	class AliTRDCalROC : public TObject
	25	{
	26
	27	public:
	28
	29	AliTRDCalROC();
	30	AliTRDCalROC(Int_t p, Int_t c);
	31	AliTRDCalROC(const AliTRDCalROC &c);
	32	virtual ~AliTRDCalROC();
	33	AliTRDCalROC &operator=(const AliTRDCalROC &c);
	34	virtual void Copy(TObject &c) const;
	35
	36	Int_t GetNrows() const { return fNrows; };
	37	Int_t GetNcols() const { return fNcols; };
	38
	39	Int_t GetChannel(Int_t c, Int_t r) const { return r+c*fNrows; };
	40	Int_t GetNchannels() const { return fNchannels; };
	41
	42	Float_t GetValue(Int_t ich) const { return (Float_t) fData[ich] / 10000; };
	43	Float_t GetValue(Int_t col, Int_t row) { return GetValue(GetChannel(col,row)); };
	44
	45	void SetValue(Int_t ich, Float_t value) { fData[ich] = (UShort_t) (value * 10000); };
	46	void SetValue(Int_t col, Int_t row, Float_t value)
	47	{ SetValue(GetChannel(col,row), value); };
	48
	49	// statistic
	50	//
	51	Double_t GetMean(AliTRDCalROC * const outlierROC=0) const;
	52	Double_t GetRMS(AliTRDCalROC * const outlierROC=0) const;
	53	Double_t GetMedian(AliTRDCalROC * const outlierROC=0) const;
	54	Double_t GetLTM(Double_t sigma=0, Double_t fraction=0.9, AliTRDCalROC const outlierROC=0);
	55
	56	// algebra
	57	Bool_t Add(Float_t c1);
	58	Bool_t Multiply(Float_t c1);
	59	Bool_t Add(const AliTRDCalROC * roc, Double_t c1 = 1);
	60	Bool_t Multiply(const AliTRDCalROC * roc);
	61	Bool_t Divide(const AliTRDCalROC * roc);
	62
	63	// noise
	64	Bool_t Unfold();
	65
	66	//Plots
	67	TH2F * MakeHisto2D(Float_t min, Float_t max,Int_t type, Float_t mu = 1.0);
	68	TH1F * MakeHisto1D(Float_t min, Float_t max,Int_t type, Float_t mu = 1.0);
	69
	70	protected:
	71
	72	Int_t fPla; // Plane number
	73	Int_t fCha; // Chamber number
	74
	75	Int_t fNrows; // Number of rows
	76	Int_t fNcols; // Number of columns
	77
	78	Int_t fNchannels; // Number of channels
	79	UShort_t *fData; //[fNchannels] Data
	80
	81	ClassDef(AliTRDCalROC, 2) // TRD ROC calibration class
	82
	83	};
	84
	85	#endif