AliFMDGainDA.h

Go to the documentation of this file.

#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:
  typedef AliFMDBaseDA::Array Array;
  AliFMDGainDA() ;
  AliFMDGainDA(const AliFMDGainDA & gainDA) ;
  AliFMDGainDA& operator=(const AliFMDGainDA&) { return *this; }
  virtual ~AliFMDGainDA();
  Bool_t OpenFiles(Bool_t appendRun=false);
  void Init();
  void SetMaxPulse(Int_t highPulse = 256) {fHighPulse = highPulse; }
  void SetNumberOfStrips(Int_t nStrips) {fNumberOfStripsPerChip = nStrips;}

protected:
  void AddChannelContainer(Array* sectorArray, 
         UShort_t det, Char_t ring, 
         UShort_t sec, UShort_t strip);
  virtual void AddSectorSummary(Array* secArray, UShort_t det, 
        Char_t ring, UShort_t sector, 
        UShort_t nStrip);
  void FillChannels(AliFMDDigit* digit);
  void Analyse(UShort_t det, Char_t ring, UShort_t sec, UShort_t strip);
  void WriteHeaderToFile();
  void UpdatePulseAndADC(UShort_t det, 
       Char_t ring, 
       UShort_t sec, 
       UShort_t strip);
  void ResetPulseAndUpdateChannel();
  void FinishEvent();
  void Terminate(TFile* dummy);
  virtual void InitContainer(TDirectory* dir);

private:
  TH1S* GetChannelHistogram(UShort_t det, Char_t ring, UShort_t sec, 
          UShort_t va);
  TGraphErrors* GetChannel(UShort_t det, Char_t ring, 
         UShort_t sec, UShort_t strip);
  TH1F* GetSectorSummary(UShort_t det, Char_t   ring, UShort_t sec);
  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: