[u/mrichter/AliRoot.git] / ITS / AliITSOnlineSDDInjectors.h

#ifndef ALIITSONLINESDDINJECTORS_H
#define ALIITSONLINESDDINJECTORS_H


///////////////////////////////////////////////////////////////////
//                                                               //
// Class used for SDD injector analysis                           //
// Origin: F.Prino, Torino, prino@to.infn.it                     //
//                                                               //
///////////////////////////////////////////////////////////////////
#include "AliITSOnlineSDD.h"

/* $Id$ */
class TH2F;
class TGraphErrors;
class AliITSOnlineSDDInjectors : public AliITSOnlineSDD {

 public:
  AliITSOnlineSDDInjectors();      
  AliITSOnlineSDDInjectors(Int_t mod, Int_t sid);
  virtual ~AliITSOnlineSDDInjectors();

  void SetSide(Int_t sid){fSide=sid;}
  void SetThreshold(Float_t thr=75.){fThreshold=thr;}
  void SetRangeLine1(Int_t tbmin=40, Int_t tbmax=90){
    fTbMin[0]=tbmin; fTbMax[0]=tbmax; 
  }
  void SetRangeLine2(Int_t tbmin=90, Int_t tbmax=140){
    fTbMin[1]=tbmin; fTbMax[1]=tbmax; 
  }
  void SetRangeLine3(Int_t tbmin=170, Int_t tbmax=220){
    fTbMin[2]=tbmin; fTbMax[2]=tbmax; 
  }
  void SetPolOrder(Int_t n=3){fPolOrder=n;}
  void SetMinDriftVel(Float_t vmin=4.){fMinDriftVel=vmin;}
  void SetMaxDriftVel(Float_t vmax=9.){fMaxDriftVel=vmax;}

  TGraphErrors* GetLineGraph(Int_t jlin);
  TGraphErrors* GetDriftVelocityGraph() const;
  Float_t* GetDriftVelFitParam()const{ return fParam;}
  Float_t GetDriftVelocity(Int_t jlin) const{return fDriftVel[jlin];}
  Float_t GetSigmaDriftVelocity(Int_t jlin) const{return fSigmaDriftVel[jlin];}
  Float_t GetTimeBinZero() const{return fTbZero;}
  Float_t GetDriftCoordinate(Float_t cAnode, Float_t cTimeBin);
  Int_t GetAnodeNumber(Int_t iInjLine) const;
  Float_t GetCentroid(Int_t injnumb, Int_t injline) const {
    if(injnumb<kNInjectors && injline<3) return fCentroid[injnumb][injline];
    else return -9999.;
  }
  Bool_t IsInjectorGood(Int_t injnumb, Int_t injline) const {
    if(injnumb<kNInjectors && injline<3) return fGoodInj[injnumb][injline];
    else return 0;
  }
  void PrintInjMap();
  void PrintCentroids();
  void WriteToASCII(Int_t evNumb, UInt_t timeStamp, Int_t optAppend=0);

  void Reset();      
  void AnalyzeEvent(TH2F* his);      
  void FindGoodInjectors();
  void FindCentroids();
  void CalcDriftVelocity(Int_t jlin);
  void CalcTimeBinZero();
  void FitDriftVelocityVsAnode();

 protected:
  void SetPositions();
 private:

  enum {
    kNInjectors = 33
  };

  AliITSOnlineSDDInjectors(const AliITSOnlineSDDInjectors& source);
  AliITSOnlineSDDInjectors& operator = (const AliITSOnlineSDDInjectors& source);
  static const Float_t fgkSaturation;
  static const Float_t fgkJitterTB;

  TH2F* fHisto;                         // histogram of module channel counts
  Float_t fTbZero;                      // Time zero for injector event
  Float_t fPosition[3];                 // Coordinates of injector lines
  UShort_t fTbMin[3];                   // Minimum time bin for each line
  UShort_t fTbMax[3];                   // Maximum time bin for each line
  Bool_t fGoodInj[kNInjectors][3];      // array of good injectors
  Float_t fCentroid[kNInjectors][3];    // array of time centroids of injectors
  Float_t fRMSCentroid[kNInjectors][3]; // array of time rms of injectors
  Float_t fDriftVel[kNInjectors];       // drift velocity
  Float_t fSigmaDriftVel[kNInjectors];  // error on drift velocity
  Float_t *fParam;                      // parameters of polinomial fit
                                        //  of drift vel. vs. anode number
  Int_t fPolOrder;                      // order of polinomial fit
  Float_t fMinDriftVel;                 // Cut value for minimum drift speed
  Float_t fMaxDriftVel;                 // Cut value for maximum drift speed
  Float_t fThreshold;                   // Threshold for injector signal

  ClassDef(AliITSOnlineSDDInjectors,1)
};
#endif
Commit	Line	Data
348f80b7	1	#ifndef ALIITSONLINESDDINJECTORS_H
	2	#define ALIITSONLINESDDINJECTORS_H
	3
	4
	5	///////////////////////////////////////////////////////////////////
	6	// //
	7	// Class used for SDD injector analysis //
	8	// Origin: F.Prino, Torino, prino@to.infn.it //
	9	// //
	10	///////////////////////////////////////////////////////////////////
0e5e647a	11	#include "AliITSOnlineSDD.h"
348f80b7	12
c668f182	13	/* $Id$ */
348f80b7	14	class TH2F;
	15	class TGraphErrors;
	16	class AliITSOnlineSDDInjectors : public AliITSOnlineSDD {
	17
	18	public:
	19	AliITSOnlineSDDInjectors();
	20	AliITSOnlineSDDInjectors(Int_t mod, Int_t sid);
	21	virtual ~AliITSOnlineSDDInjectors();
	22
	23	void SetSide(Int_t sid){fSide=sid;}
	24	void SetThreshold(Float_t thr=75.){fThreshold=thr;}
	25	void SetRangeLine1(Int_t tbmin=40, Int_t tbmax=90){
	26	fTbMin[0]=tbmin; fTbMax[0]=tbmax;
	27	}
	28	void SetRangeLine2(Int_t tbmin=90, Int_t tbmax=140){
	29	fTbMin[1]=tbmin; fTbMax[1]=tbmax;
	30	}
	31	void SetRangeLine3(Int_t tbmin=170, Int_t tbmax=220){
	32	fTbMin[2]=tbmin; fTbMax[2]=tbmax;
	33	}
	34	void SetPolOrder(Int_t n=3){fPolOrder=n;}
	35	void SetMinDriftVel(Float_t vmin=4.){fMinDriftVel=vmin;}
	36	void SetMaxDriftVel(Float_t vmax=9.){fMaxDriftVel=vmax;}
	37
	38	TGraphErrors* GetLineGraph(Int_t jlin);
beb262b4	39	TGraphErrors* GetDriftVelocityGraph() const;
348f80b7	40	Float_t* GetDriftVelFitParam()const{ return fParam;}
	41	Float_t GetDriftVelocity(Int_t jlin) const{return fDriftVel[jlin];}
	42	Float_t GetSigmaDriftVelocity(Int_t jlin) const{return fSigmaDriftVel[jlin];}
	43	Float_t GetTimeBinZero() const{return fTbZero;}
	44	Float_t GetDriftCoordinate(Float_t cAnode, Float_t cTimeBin);
beb262b4	45	Int_t GetAnodeNumber(Int_t iInjLine) const;
c668f182	46	Float_t GetCentroid(Int_t injnumb, Int_t injline) const {
	47	if(injnumb<kNInjectors && injline<3) return fCentroid[injnumb][injline];
	48	else return -9999.;
	49	}
	50	Bool_t IsInjectorGood(Int_t injnumb, Int_t injline) const {
	51	if(injnumb<kNInjectors && injline<3) return fGoodInj[injnumb][injline];
	52	else return 0;
	53	}
348f80b7	54	void PrintInjMap();
348f80b7	55	void PrintCentroids();
4c82df4c	56	void WriteToASCII(Int_t evNumb, UInt_t timeStamp, Int_t optAppend=0);
348f80b7	57
	58	void Reset();
	59	void AnalyzeEvent(TH2F* his);
	60	void FindGoodInjectors();
	61	void FindCentroids();
	62	void CalcDriftVelocity(Int_t jlin);
	63	void CalcTimeBinZero();
	64	void FitDriftVelocityVsAnode();
	65
	66	protected:
	67	void SetPositions();
	68	private:
	69
	70	enum {
	71	kNInjectors = 33
	72	};
	73
	74	AliITSOnlineSDDInjectors(const AliITSOnlineSDDInjectors& source);
	75	AliITSOnlineSDDInjectors& operator = (const AliITSOnlineSDDInjectors& source);
	76	static const Float_t fgkSaturation;
	77	static const Float_t fgkJitterTB;
	78
beb262b4	79	TH2F* fHisto; // histogram of module channel counts
	80	Float_t fTbZero; // Time zero for injector event
	81	Float_t fPosition[3]; // Coordinates of injector lines
	82	UShort_t fTbMin[3]; // Minimum time bin for each line
	83	UShort_t fTbMax[3]; // Maximum time bin for each line
	84	Bool_t fGoodInj[kNInjectors][3]; // array of good injectors
	85	Float_t fCentroid[kNInjectors][3]; // array of time centroids of injectors
	86	Float_t fRMSCentroid[kNInjectors][3]; // array of time rms of injectors
	87	Float_t fDriftVel[kNInjectors]; // drift velocity
	88	Float_t fSigmaDriftVel[kNInjectors]; // error on drift velocity
	89	Float_t *fParam; // parameters of polinomial fit
	90	// of drift vel. vs. anode number
	91	Int_t fPolOrder; // order of polinomial fit
	92	Float_t fMinDriftVel; // Cut value for minimum drift speed
	93	Float_t fMaxDriftVel; // Cut value for maximum drift speed
	94	Float_t fThreshold; // Threshold for injector signal
348f80b7	95
	96	ClassDef(AliITSOnlineSDDInjectors,1)
	97	};
	98	#endif