[u/mrichter/AliRoot.git] / FMD / AliFMDGainDA.h

#ifndef ALIFMDGAINDA_H
#define ALIFMDGAINDA_H
/* Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights
 * reserved. 
 *
 * See cxx source for full Copyright notice                               
 */
// This class implements the Gain detector algorithm (DA) for the FMD. 
// The gain is the response of the VA chips to a known pulse and has to be 
// calculated strip-by-strip, pulse-by-pulse. 


#include "AliFMDBaseDA.h"
#include "TH1F.h"
#include "TObjArray.h"
// #include "TGraphErrors.h"
class TGraphErrors;
class TH1S;

class AliFMDGainDA: public AliFMDBaseDA 
{
public:
  /** 
   * Constructor 
   * 
   */
  AliFMDGainDA() ;
  /** 
   * Copy constructor 
   * 
   * @param gainDA  Object to copy from
   */  
  AliFMDGainDA(const AliFMDGainDA & gainDA) ;
  /** 
   * Assignment operator 
   * 
   * @param gainDA Object to assign from
   */  
  AliFMDGainDA& operator=(const AliFMDGainDA&) { return *this; }
  /** 
   * Destructor 
   * 
   */
  virtual ~AliFMDGainDA();
  /** 
   * Initialize 
   * 
   */
  void Init();
  /** 
   * Set the maximum pulse size 
   * 
   * @param highPulse Maximum pulse size
   */
  void SetMaxPulse(Int_t highPulse = 256) {fHighPulse = highPulse; }
  /** 
   * Set the number of strips per input channel 
   * 
   * @param nStrips Number of strips per channel
   */
  void SetNumberOfStrips(Int_t nStrips) {fNumberOfStripsPerChip = nStrips;}

protected:
  /** 
   * Make a container for a channel 
   * 
   * @param sectorArray Sectors 
   * @param det         Detector # 
   * @param ring        Ring identifier 
   * @param sec         Sector number
   * @param strip       Strip number
   */ 
  void AddChannelContainer(TObjArray* sectorArray, 
			   UShort_t det, Char_t ring, 
			   UShort_t sec, UShort_t strip);
  /** 
   * Fill channel histogram 
   * 
   * @param digit Digit to fill from
   */
  void FillChannels(AliFMDDigit* digit);
  /** 
   * Analyse the result
   * 
   * @param det   Detector # 
   * @param ring  Ring identifier 
   * @param sec   Sector number
   * @param strip Strip number
   * @param h     Summary histogram with bins for sector and strip
   */
  void Analyse(UShort_t det, Char_t ring, UShort_t sec, UShort_t strip);
  /**
   * Write out the header 
   */
  void WriteHeaderToFile();
  /** 
   * Go to the next sample
   * 
   * @param det         Detector # 
   * @param ring        Ring identifier 
   * @param sec         Sector number
   * @param strip       Strip number
   */
  void UpdatePulseAndADC(UShort_t det, 
			 Char_t ring, 
			 UShort_t sec, 
			 UShort_t strip);
  /** 
   * Reset all 
   * 
   */
  void ResetPulseAndUpdateChannel();
  /** 
   * End of event 
   * 
   */
  void FinishEvent();
  /** 
   * End of job
   * 
   * @param dummy Not used
   */
  void Terminate(TFile* dummy);

private:
  /** 
   * Get the channel histogram
   * 
   * @param det         Detector # 
   * @param ring        Ring identifier 
   * @param sec         Sector number
   * @param strip       Strip number
   * 
   * @return Histogram
   */  
  TH1S* GetChannelHistogram(UShort_t det, Char_t ring, UShort_t sec, 
			    UShort_t strip);
  /** 
   * Get strip graph
   * 
   * @param det         Detector # 
   * @param ring        Ring identifier 
   * @param sec         Sector number
   * @param strip       Strip number
   * 
   * @return Graph
   */
  TGraphErrors* GetChannel(UShort_t det, Char_t ring, 
			   UShort_t sec, UShort_t strip);
  TObjArray fGainArray;          // Array of gains 
  Int_t     fHighPulse;          // Highest pulse
  TArrayS   fEventsPerChannel;   // # of events per pulse step
  TArrayS   fCurrentPulse;       // The current pulse size 
  TArrayS   fCurrentChannel;     // The current strip number
  Int_t     fNumberOfStripsPerChip; // Number of strips
  
  TH1F      fSummaryGains;         // Summary histogram 
  Int_t     fCurrentSummaryStrip;  // Current strip for summary

  void  MakeSummary(UShort_t det, Char_t ring);

  TH2* fGainFMD1i; // AMORE DQM histogram
  TH2* fGainFMD2i; // AMORE DQM histogram
  TH2* fGainFMD2o; // AMORE DQM histogram
  TH2* fGainFMD3i; // AMORE DQM histogram
  TH2* fGainFMD3o; // AMORE DQM histogram
  TH2* fChi2FMD1i; // AMORE DQM histogram
  TH2* fChi2FMD2i; // AMORE DQM histogram
  TH2* fChi2FMD2o; // AMORE DQM histogram
  TH2* fChi2FMD3i; // AMORE DQM histogram
  TH2* fChi2FMD3o; // AMORE DQM histogram
  
  ClassDef(AliFMDGainDA,0) // Detector algorithm for gain runs

};
#endif
// Local Variables: 
//  mode: C++ 
// End:
Commit	Line	Data
	1	#ifndef ALIFMDGAINDA_H
	2	#define ALIFMDGAINDA_H
	3	/* Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights
	4	* reserved.
	5	*
	6	* See cxx source for full Copyright notice
	7	*/
	8	// This class implements the Gain detector algorithm (DA) for the FMD.
	9	// The gain is the response of the VA chips to a known pulse and has to be
	10	// calculated strip-by-strip, pulse-by-pulse.
	11
	12
	13	#include "AliFMDBaseDA.h"
	14	#include "TH1F.h"
	15	#include "TObjArray.h"
	16	// #include "TGraphErrors.h"
	17	class TGraphErrors;
	18	class TH1S;
	19
	20	class AliFMDGainDA: public AliFMDBaseDA
	21	{
	22	public:
	23	/**
	24	* Constructor
	25	*
	26	*/
	27	AliFMDGainDA() ;
	28	/**
	29	* Copy constructor
	30	*
	31	* @param gainDA Object to copy from
	32	*/
	33	AliFMDGainDA(const AliFMDGainDA & gainDA) ;
	34	/**
	35	* Assignment operator
	36	*
	37	* @param gainDA Object to assign from
	38	*/
	39	AliFMDGainDA& operator=(const AliFMDGainDA&) { return *this; }
	40	/**
	41	* Destructor
	42	*
	43	*/
	44	virtual ~AliFMDGainDA();
	45	/**
	46	* Initialize
	47	*
	48	*/
	49	void Init();
	50	/**
	51	* Set the maximum pulse size
	52	*
	53	* @param highPulse Maximum pulse size
	54	*/
	55	void SetMaxPulse(Int_t highPulse = 256) {fHighPulse = highPulse; }
	56	/**
	57	* Set the number of strips per input channel
	58	*
	59	* @param nStrips Number of strips per channel
	60	*/
	61	void SetNumberOfStrips(Int_t nStrips) {fNumberOfStripsPerChip = nStrips;}
	62
	63	protected:
	64	/**
	65	* Make a container for a channel
	66	*
	67	* @param sectorArray Sectors
	68	* @param det Detector #
	69	* @param ring Ring identifier
	70	* @param sec Sector number
	71	* @param strip Strip number
	72	*/
	73	void AddChannelContainer(TObjArray* sectorArray,
	74	UShort_t det, Char_t ring,
	75	UShort_t sec, UShort_t strip);
	76	/**
	77	* Fill channel histogram
	78	*
	79	* @param digit Digit to fill from
	80	*/
	81	void FillChannels(AliFMDDigit* digit);
	82	/**
	83	* Analyse the result
	84	*
	85	* @param det Detector #
	86	* @param ring Ring identifier
	87	* @param sec Sector number
	88	* @param strip Strip number
	89	* @param h Summary histogram with bins for sector and strip
	90	*/
	91	void Analyse(UShort_t det, Char_t ring, UShort_t sec, UShort_t strip);
	92	/**
	93	* Write out the header
	94	*/
	95	void WriteHeaderToFile();
	96	/**
	97	* Go to the next sample
	98	*
	99	* @param det Detector #
	100	* @param ring Ring identifier
	101	* @param sec Sector number
	102	* @param strip Strip number
	103	*/
	104	void UpdatePulseAndADC(UShort_t det,
	105	Char_t ring,
	106	UShort_t sec,
	107	UShort_t strip);
	108	/**
	109	* Reset all
	110	*
	111	*/
	112	void ResetPulseAndUpdateChannel();
	113	/**
	114	* End of event
	115	*
	116	*/
	117	void FinishEvent();
	118	/**
	119	* End of job
	120	*
	121	* @param dummy Not used
	122	*/
	123	void Terminate(TFile* dummy);
	124
	125	private:
	126	/**
	127	* Get the channel histogram
	128	*
	129	* @param det Detector #
	130	* @param ring Ring identifier
	131	* @param sec Sector number
	132	* @param strip Strip number
	133	*
	134	* @return Histogram
	135	*/
	136	TH1S* GetChannelHistogram(UShort_t det, Char_t ring, UShort_t sec,
	137	UShort_t strip);
	138	/**
	139	* Get strip graph
	140	*
	141	* @param det Detector #
	142	* @param ring Ring identifier
	143	* @param sec Sector number
	144	* @param strip Strip number
	145	*
	146	* @return Graph
	147	*/
	148	TGraphErrors* GetChannel(UShort_t det, Char_t ring,
	149	UShort_t sec, UShort_t strip);
	150	TObjArray fGainArray; // Array of gains
	151	Int_t fHighPulse; // Highest pulse
	152	TArrayS fEventsPerChannel; // # of events per pulse step
	153	TArrayS fCurrentPulse; // The current pulse size
	154	TArrayS fCurrentChannel; // The current strip number
	155	Int_t fNumberOfStripsPerChip; // Number of strips
	156
	157	TH1F fSummaryGains; // Summary histogram
	158	Int_t fCurrentSummaryStrip; // Current strip for summary
	159
	160	void MakeSummary(UShort_t det, Char_t ring);
	161
	162	TH2* fGainFMD1i; // AMORE DQM histogram
	163	TH2* fGainFMD2i; // AMORE DQM histogram
	164	TH2* fGainFMD2o; // AMORE DQM histogram
	165	TH2* fGainFMD3i; // AMORE DQM histogram
	166	TH2* fGainFMD3o; // AMORE DQM histogram
	167	TH2* fChi2FMD1i; // AMORE DQM histogram
	168	TH2* fChi2FMD2i; // AMORE DQM histogram
	169	TH2* fChi2FMD2o; // AMORE DQM histogram
	170	TH2* fChi2FMD3i; // AMORE DQM histogram
	171	TH2* fChi2FMD3o; // AMORE DQM histogram
	172
	173	ClassDef(AliFMDGainDA,0) // Detector algorithm for gain runs
	174
	175	};
	176	#endif
	177	// Local Variables:
	178	// mode: C++
	179	// End: