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

#ifndef ALIITSONLINESDDINJECTORS_H
#define ALIITSONLINESDDINJECTORS_H


/* $Id: */

///////////////////////////////////////////////////////////////////
//                                                               //
// Class used for SDD injector analysis                           //
// Origin: F.Prino, Torino, prino@to.infn.it                     //
//                                                               //
///////////////////////////////////////////////////////////////////

#include "AliITSOnlineSDD.h"


class TH2F;
class TGraphErrors;
class AliITSOnlineSDDInjectors : public AliITSOnlineSDD {

 public:
  AliITSOnlineSDDInjectors();      
  AliITSOnlineSDDInjectors(Int_t nddl, Int_t ncarlos, Int_t sid);
  virtual ~AliITSOnlineSDDInjectors();

  void SetThresholds(Float_t tl, Float_t th){
    fLowThreshold=tl;
    fHighThreshold=th;
  }
  void SetInjLineRange(Int_t jlin, Int_t tbmin, Int_t tbmax){
    fTbMin[jlin]=tbmin;
    fTbMax[jlin]=tbmax;
  }
  void SetPolOrder(Int_t n){fPolOrder=n;}
  void SetMinDriftSpeed(Float_t vmin){fMinDriftSpeed=vmin;}
  void SetMaxDriftSpeed(Float_t vmax){fMaxDriftSpeed=vmax;}
  void SetMaxDriftSpeedErr(Float_t maxval){
    fMaxDriftSpeedErr=maxval;
  }
  void SetFitLimits(Int_t firstpad,Int_t lastpad){
    fFirstPadForFit=firstpad;
    fLastPadForFit=lastpad;
  }
  void SetPadStatusCutForFit(Int_t cutval=1){
    fPadStatusCutForFit=cutval;
  }
  void SetDefaults();
  

  TGraphErrors* GetTimeVsDistGraph(Int_t jpad) const;
  TGraphErrors* GetDriftSpeedGraph() const;
  Float_t* GetDriftSpeedFitParam()const{ return fParam;}
  Float_t GetDriftSpeeed(Int_t jpad) const{return fDriftSpeed[jpad];}
  Float_t GetDriftSpeedErr(Int_t jpad) const{return fDriftSpeedErr[jpad];}
  Float_t GetTimeBinZero() const{return fTbZero;}

  Int_t GetAnodeNumber(Int_t iInjPad) const;
  Int_t GetInjPadNumberFromAnode(Int_t nAnode) const;
  Int_t GetInjPadStatus(Int_t jpad) const;  
  Int_t GetAnodeStatus(Int_t nAnode) const{
    Int_t jpad=GetInjPadNumberFromAnode(nAnode);
    return GetInjPadStatus(jpad);
  }  
  Float_t GetCentroid(Int_t jpad, Int_t jlin) const {
    if(jpad<kInjPads && jlin<kInjLines) return fCentroid[jpad][jlin];
    else return -9999.;
  }
  Bool_t IsInjectorGood(Int_t jpad, Int_t jlin) const {
    if(jpad<kInjPads && jlin<kInjLines) return fGoodInj[jpad][jlin];
    else return 0;
  }
  void PrintInjectorStatus();
  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 CalcDriftSpeed(Int_t jpad);
  void CalcTimeBinZero();
  void FitDriftSpeedVsAnode();

 protected:
  void SetPositions();
 private:

  enum {kInjPads  = 33};
  enum {kInjLines = 3};

  AliITSOnlineSDDInjectors(const AliITSOnlineSDDInjectors& source);
  AliITSOnlineSDDInjectors& operator = (const AliITSOnlineSDDInjectors& source);
  static const Float_t fgkSaturation;        // ADC saturation value (1008)
  static const Float_t fgkDefaultLThreshold;  // Default for fLowThreshold
  static const Float_t fgkDefaultHThreshold;  // Default for fHighThreshold
  static const Float_t fgkDefaultMinSpeed;   // Default for fMinDriftSpeed
  static const Float_t fgkDefaultMaxSpeed;   // Default for fMaxDriftSpeed
  static const Float_t fgkDefaultMaxErr;     // Default for fMaxDriftSpeedErr
  static const Int_t   fgkDefaultPolOrder;   // Default for fPolOrder
  static const UShort_t   fgkDefaultTbMin[kInjLines];  // Defaults for fTbMin
  static const UShort_t   fgkDefaultTbMax[kInjLines];  // Defaults for fTbMax


  TH2F* fHisto;                              // histogram of channel counts
  Float_t fTbZero;                           // Time zero for injector event
  Float_t fPosition[kInjLines];              // Coordinates of injector lines
  UShort_t fTbMin[kInjLines];                // Minimum time bin for each line
  UShort_t fTbMax[kInjLines];                // Maximum time bin for each line
  Bool_t fGoodInj[kInjPads][kInjLines];      // array of good injectors
  Float_t fCentroid[kInjPads][kInjLines];    // array of time bin centroids
  Float_t fRMSCentroid[kInjPads][kInjLines]; // array of time rms of injectors
  Float_t fDriftSpeed[kInjPads];             // drift speed
  Float_t fDriftSpeedErr[kInjPads];          // error on drift speed
  Float_t *fParam;                           // parameters of polinomial fit to
                                             // drift speed vs. anode number
  Int_t fPolOrder;                   // order of polinomial fit
  Float_t fMinDriftSpeed;            // Minimum value for drift speed
  Float_t fMaxDriftSpeed;            // Maximum value for drift speed
  Float_t fMaxDriftSpeedErr;         // Maximum value for error on drift speed
  Float_t fLowThreshold;             // Low threshold for injector signal
  Float_t fHighThreshold;            // High threshold for injector signal

  Int_t fFirstPadForFit;             // first injector pad used in fit
  Int_t fLastPadForFit;              // last injector pad used in fit
  Int_t fPadStatusCutForFit;         // minimum value of pad status for fit


  ClassDef(AliITSOnlineSDDInjectors,3)
};
#endif
Commit	Line	Data
348f80b7	1	#ifndef ALIITSONLINESDDINJECTORS_H
	2	#define ALIITSONLINESDDINJECTORS_H
	3
	4
9f026db8	5	/* $Id: */
9f026db8	6
348f80b7	7	///////////////////////////////////////////////////////////////////
	8	// //
	9	// Class used for SDD injector analysis //
	10	// Origin: F.Prino, Torino, prino@to.infn.it //
	11	// //
	12	///////////////////////////////////////////////////////////////////
9f026db8	13
0e5e647a	14	#include "AliITSOnlineSDD.h"
348f80b7	15
9f026db8	16
348f80b7	17	class TH2F;
	18	class TGraphErrors;
	19	class AliITSOnlineSDDInjectors : public AliITSOnlineSDD {
	20
	21	public:
	22	AliITSOnlineSDDInjectors();
979b5a5f	23	AliITSOnlineSDDInjectors(Int_t nddl, Int_t ncarlos, Int_t sid);
348f80b7	24	virtual ~AliITSOnlineSDDInjectors();
348f80b7	25
9f026db8	26	void SetThresholds(Float_t tl, Float_t th){
	27	fLowThreshold=tl;
	28	fHighThreshold=th;
	29	}
	30	void SetInjLineRange(Int_t jlin, Int_t tbmin, Int_t tbmax){
	31	fTbMin[jlin]=tbmin;
	32	fTbMax[jlin]=tbmax;
	33	}
	34	void SetPolOrder(Int_t n){fPolOrder=n;}
	35	void SetMinDriftSpeed(Float_t vmin){fMinDriftSpeed=vmin;}
	36	void SetMaxDriftSpeed(Float_t vmax){fMaxDriftSpeed=vmax;}
	37	void SetMaxDriftSpeedErr(Float_t maxval){
	38	fMaxDriftSpeedErr=maxval;
348f80b7	39	}
9f026db8	40	void SetFitLimits(Int_t firstpad,Int_t lastpad){
	41	fFirstPadForFit=firstpad;
	42	fLastPadForFit=lastpad;
348f80b7	43	}
9f026db8	44	void SetPadStatusCutForFit(Int_t cutval=1){
9f026db8	45	fPadStatusCutForFit=cutval;
348f80b7	46	}
9f026db8	47	void SetDefaults();
	48
	49
	50	TGraphErrors* GetTimeVsDistGraph(Int_t jpad) const;
	51	TGraphErrors* GetDriftSpeedGraph() const;
	52	Float_t* GetDriftSpeedFitParam()const{ return fParam;}
	53	Float_t GetDriftSpeeed(Int_t jpad) const{return fDriftSpeed[jpad];}
	54	Float_t GetDriftSpeedErr(Int_t jpad) const{return fDriftSpeedErr[jpad];}
348f80b7	55	Float_t GetTimeBinZero() const{return fTbZero;}
9f026db8	56
	57	Int_t GetAnodeNumber(Int_t iInjPad) const;
	58	Int_t GetInjPadNumberFromAnode(Int_t nAnode) const;
	59	Int_t GetInjPadStatus(Int_t jpad) const;
	60	Int_t GetAnodeStatus(Int_t nAnode) const{
	61	Int_t jpad=GetInjPadNumberFromAnode(nAnode);
	62	return GetInjPadStatus(jpad);
	63	}
	64	Float_t GetCentroid(Int_t jpad, Int_t jlin) const {
	65	if(jpad<kInjPads && jlin<kInjLines) return fCentroid[jpad][jlin];
c668f182	66	else return -9999.;
c668f182	67	}
9f026db8	68	Bool_t IsInjectorGood(Int_t jpad, Int_t jlin) const {
9f026db8	69	if(jpad<kInjPads && jlin<kInjLines) return fGoodInj[jpad][jlin];
c668f182	70	else return 0;
c668f182	71	}
9f026db8	72	void PrintInjectorStatus();
348f80b7	73	void PrintCentroids();
4c82df4c	74	void WriteToASCII(Int_t evNumb, UInt_t timeStamp, Int_t optAppend=0);
348f80b7	75
	76	void Reset();
	77	void AnalyzeEvent(TH2F* his);
	78	void FindGoodInjectors();
	79	void FindCentroids();
9f026db8	80	void CalcDriftSpeed(Int_t jpad);
348f80b7	81	void CalcTimeBinZero();
9f026db8	82	void FitDriftSpeedVsAnode();
348f80b7	83
	84	protected:
	85	void SetPositions();
	86	private:
	87
9f026db8	88	enum {kInjPads = 33};
9f026db8	89	enum {kInjLines = 3};
348f80b7	90
	91	AliITSOnlineSDDInjectors(const AliITSOnlineSDDInjectors& source);
	92	AliITSOnlineSDDInjectors& operator = (const AliITSOnlineSDDInjectors& source);
9f026db8	93	static const Float_t fgkSaturation; // ADC saturation value (1008)
	94	static const Float_t fgkDefaultLThreshold; // Default for fLowThreshold
	95	static const Float_t fgkDefaultHThreshold; // Default for fHighThreshold
	96	static const Float_t fgkDefaultMinSpeed; // Default for fMinDriftSpeed
	97	static const Float_t fgkDefaultMaxSpeed; // Default for fMaxDriftSpeed
	98	static const Float_t fgkDefaultMaxErr; // Default for fMaxDriftSpeedErr
	99	static const Int_t fgkDefaultPolOrder; // Default for fPolOrder
	100	static const UShort_t fgkDefaultTbMin[kInjLines]; // Defaults for fTbMin
	101	static const UShort_t fgkDefaultTbMax[kInjLines]; // Defaults for fTbMax
	102
	103
	104	TH2F* fHisto; // histogram of channel counts
	105	Float_t fTbZero; // Time zero for injector event
	106	Float_t fPosition[kInjLines]; // Coordinates of injector lines
	107	UShort_t fTbMin[kInjLines]; // Minimum time bin for each line
	108	UShort_t fTbMax[kInjLines]; // Maximum time bin for each line
	109	Bool_t fGoodInj[kInjPads][kInjLines]; // array of good injectors
	110	Float_t fCentroid[kInjPads][kInjLines]; // array of time bin centroids
	111	Float_t fRMSCentroid[kInjPads][kInjLines]; // array of time rms of injectors
	112	Float_t fDriftSpeed[kInjPads]; // drift speed
	113	Float_t fDriftSpeedErr[kInjPads]; // error on drift speed
	114	Float_t *fParam; // parameters of polinomial fit to
	115	// drift speed vs. anode number
	116	Int_t fPolOrder; // order of polinomial fit
	117	Float_t fMinDriftSpeed; // Minimum value for drift speed
	118	Float_t fMaxDriftSpeed; // Maximum value for drift speed
	119	Float_t fMaxDriftSpeedErr; // Maximum value for error on drift speed
	120	Float_t fLowThreshold; // Low threshold for injector signal
	121	Float_t fHighThreshold; // High threshold for injector signal
	122
	123	Int_t fFirstPadForFit; // first injector pad used in fit
	124	Int_t fLastPadForFit; // last injector pad used in fit
	125	Int_t fPadStatusCutForFit; // minimum value of pad status for fit
	126
	127
	128	ClassDef(AliITSOnlineSDDInjectors,3)
348f80b7	129	};
348f80b7	130	#endif