[u/mrichter/AliRoot.git] / HMPID / AliHMPIDCalib.h

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

// Class of HMPID to manage digits ---> pads
//.
//.
//.

//#include "TTreePlayer.h"
//#include <TTree.h>
#include <TH1.h>
#include <TH2.h>
//#include <TH1S.h>
#include <TMath.h>
#include <TFile.h>
#include <TString.h>
#include <THnSparse.h>
#include "AliHMPIDParam.h"
#include "AliHMPIDRawStream.h"

class TFile;
class AliHMPIDCalib: public TObject { 


public:
  AliHMPIDCalib();
  virtual ~AliHMPIDCalib();
          void Init();
          void FillPedestal(Int_t pad,Int_t q);                             //absolute pad number and the charge of the pad
          void FillErrors(Int_t nDDL,Int_t nErrType, Int_t nErr);           //Fill the errors from RawStream
          void FillDDLCnt(Int_t iddl,Int_t inDDL, Int_t outDDL);            //Fill the errors from RawStream
        Bool_t CalcPedestal(Int_t nDDL, Char_t* name, Char_t *name2,Int_t nEv);           //number of the DDL, name of the output file and the number of events processed
        
        Bool_t WriteErrors(Int_t nDDL, Char_t* name, Int_t nEv);            //number of the DDL, name of the output file and the number of events processed
         void SetRunParams(ULong_t runNum,Int_t timeStamp, Int_t ldcId);    //Set Run Parameters such as Run Number, TimeStamp, LDCid 
         void SetSigCut(Int_t nSigCut) { fSigCut=nSigCut;}                  //Set Sigma Cuts from Setter
         void SetSigCutFromFile(TString name);                              //Set Sigma Cuts from File
         void SetDeadChannelMapFromFile(TString name);
         Bool_t GetSelectedDDL()     const  {return fSelectDDL;}              //Set wether ADC histos of pads are written or not
         TH2F   *GetPedMeanMap(Int_t iDDL)  {return fPedMeanMap[iDDL];}       //Get the pedestal mean map for a DDL to send to AMORE
         TH2F   *GetPedSigMap(Int_t iDDL)   {return fPedSigMap[iDDL];}        //Get the pedestal sigma map for a DDL to send to AMORE
         TH1F   *GetPedMean(Int_t iChFee)   {return f1DPedMean[iChFee];}      //Get the pedestal mean map for a FEE channel to send to AMORE
         TH1F   *GetPedSigma(Int_t iChFee)  {return f1DPedSigma[iChFee];}     //Get the pedestal Sigma map for a FEE channel to send to AMORE
     THnSparse  *GetDeadMap()               {return fDeadMap;}                //Get the masked channel map from the DAQ database
         Int_t   GetNumMaskedPads()         {return fNumMaskedPads;}          //Get the number of masked channels
         Int_t   GetNumDeadPads()           {return fNumDeadPads;}            //Get the number of masked channels
protected: 

    Bool_t     *faddl;                                                         //check is ddl is filled
    Float_t ****fsq;                                                           //Sum of pad Q
    Float_t ****fsq2;                                                          //Sum of pad Q^2
    Int_t   ****fnpc;                                                          //# of the pad was called with non zero charge
    Int_t   ****fpedQ0;                                                        //Check how many times a pad gives 0 charge in pedestal runs
    Int_t     **fErr;                                                          // Store the numner of errors for a given error type and a given DDL
    TH1I      **fPadAdc;                                                       //Charge distribution for pads    
    Bool_t     *fIsPad;                                                        //Check if the ADC histo for the pad is booked or not
    TFile      *fFile;                                                         //ADC histo output file (one per LDC)      
    UInt_t      fLdcId;                                                        //Ldc ID 
    UInt_t      fTimeStamp;                                                    //Time Stamp
    Int_t       fRunNum;                                                       //Run Number
    Int_t       fSigCut;                                                       //n. of pedestal distribution sigmas used to create zero suppresion table                          
    Int_t      *fnDDLInStream;                                                 // if the DDL is in the raw data
    Int_t      *fnDDLOutStream;                                                // if the DDL is in the raw data
    Bool_t      fLargeHisto;                                                   //Default is kFALSE.if kTRUE then write large pad histograms with 4093 bins!!!! Only if you have 2GB of RAM!!!   
    Int_t       fSelectDDL;                                                    //Select the DDL to write for the in the large histograms. Only ONE at one time!
    THnSparse  *fDeadMap;                                                      //Dead Channel Map
    TH2F       **fPedMeanMap;                                                  //2D mean pedestal map to export to AMORE
    TH2F       **fPedSigMap;                                                   //2D pedestal sigma map to export to AMORE
    TH1F      **f1DPedMean;                                                    //1D mean pedestal map to export to AMORE
    TH1F      **f1DPedSigma;                                                   //1D pedestal sigma map to export to AMORE
    Int_t       fNumMaskedPads;                                                //Number of masked pads     
    Int_t       fNumDeadPads;                                                  //Number of currently dead channels   
            
  private:
                                           
  AliHMPIDCalib(const AliHMPIDCalib& c);                                       //dummy copy constructor
  AliHMPIDCalib &operator=(const AliHMPIDCalib& c);                            //dummy assignment operator
     
    ClassDef(AliHMPIDCalib,5)                                                  //HMPID calibration and pedestal class        
};
#endif
Commit	Line	Data
21f61e25	1	#ifndef AliHMPIDCalib_h
	2	#define AliHMPIDCalib_h
	3	/* Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved. *
	4	* See cxx source for full Copyright notice */
	5
	6	// Class of HMPID to manage digits ---> pads
	7	//.
	8	//.
	9	//.
3e60bb92	10
	11	//#include "TTreePlayer.h"
	12	//#include <TTree.h>
21f61e25	13	#include <TH1.h>
0f1281b2	14	#include <TH2.h>
606697a8	15	//#include <TH1S.h>
3e60bb92	16	#include <TMath.h>
3c8e86a0	17	#include <TFile.h>
cdd8604a	18	#include <TString.h>
2ac899f2	19	#include <THnSparse.h>
0f1281b2	20	#include "AliHMPIDParam.h"
0f1281b2	21	#include "AliHMPIDRawStream.h"
21f61e25	22
3c8e86a0	23	class TFile;
21f61e25	24	class AliHMPIDCalib: public TObject {
	25
	26
	27	public:
	28	AliHMPIDCalib();
	29	virtual ~AliHMPIDCalib();
	30	void Init();
0f1281b2	31	void FillPedestal(Int_t pad,Int_t q); //absolute pad number and the charge of the pad
3c8e86a0	32	void FillErrors(Int_t nDDL,Int_t nErrType, Int_t nErr); //Fill the errors from RawStream
ee812b5b	33	void FillDDLCnt(Int_t iddl,Int_t inDDL, Int_t outDDL); //Fill the errors from RawStream
4969a7fa	34	Bool_t CalcPedestal(Int_t nDDL, Char_t* name, Char_t *name2,Int_t nEv); //number of the DDL, name of the output file and the number of events processed
cdd8604a	35
3c8e86a0	36	Bool_t WriteErrors(Int_t nDDL, Char_t* name, Int_t nEv); //number of the DDL, name of the output file and the number of events processed
3c8e86a0	37	void SetRunParams(ULong_t runNum,Int_t timeStamp, Int_t ldcId); //Set Run Parameters such as Run Number, TimeStamp, LDCid
b38ac33a	38	void SetSigCut(Int_t nSigCut) { fSigCut=nSigCut;} //Set Sigma Cuts from Setter
4969a7fa	39	void SetSigCutFromFile(TString name); //Set Sigma Cuts from File
2ac899f2	40	void SetDeadChannelMapFromFile(TString name);
2ac899f2	41	Bool_t GetSelectedDDL() const {return fSelectDDL;} //Set wether ADC histos of pads are written or not
2a7784f8	42	TH2F *GetPedMeanMap(Int_t iDDL) {return fPedMeanMap[iDDL];} //Get the pedestal mean map for a DDL to send to AMORE
	43	TH2F *GetPedSigMap(Int_t iDDL) {return fPedSigMap[iDDL];} //Get the pedestal sigma map for a DDL to send to AMORE
	44	TH1F *GetPedMean(Int_t iChFee) {return f1DPedMean[iChFee];} //Get the pedestal mean map for a FEE channel to send to AMORE
	45	TH1F *GetPedSigma(Int_t iChFee) {return f1DPedSigma[iChFee];} //Get the pedestal Sigma map for a FEE channel to send to AMORE
	46	THnSparse *GetDeadMap() {return fDeadMap;} //Get the masked channel map from the DAQ database
	47	Int_t GetNumMaskedPads() {return fNumMaskedPads;} //Get the number of masked channels
	48	Int_t GetNumDeadPads() {return fNumDeadPads;} //Get the number of masked channels
21f61e25	49	protected:
3c8e86a0	50
ee812b5b	51	Bool_t *faddl; //check is ddl is filled
	52	Float_t ****fsq; //Sum of pad Q
	53	Float_t ****fsq2; //Sum of pad Q^2
	54	Int_t ****fnpc; //# of the pad was called with non zero charge
	55	Int_t ****fpedQ0; //Check how many times a pad gives 0 charge in pedestal runs
	56	Int_t **fErr; // Store the numner of errors for a given error type and a given DDL
	57	TH1I **fPadAdc; //Charge distribution for pads
	58	Bool_t *fIsPad; //Check if the ADC histo for the pad is booked or not
	59	TFile *fFile; //ADC histo output file (one per LDC)
	60	UInt_t fLdcId; //Ldc ID
	61	UInt_t fTimeStamp; //Time Stamp
	62	Int_t fRunNum; //Run Number
	63	Int_t fSigCut; //n. of pedestal distribution sigmas used to create zero suppresion table
ee812b5b	64	Int_t *fnDDLInStream; // if the DDL is in the raw data
	65	Int_t *fnDDLOutStream; // if the DDL is in the raw data
	66	Bool_t fLargeHisto; //Default is kFALSE.if kTRUE then write large pad histograms with 4093 bins!!!! Only if you have 2GB of RAM!!!
	67	Int_t fSelectDDL; //Select the DDL to write for the in the large histograms. Only ONE at one time!
2ac899f2	68	THnSparse *fDeadMap; //Dead Channel Map
	69	TH2F **fPedMeanMap; //2D mean pedestal map to export to AMORE
	70	TH2F **fPedSigMap; //2D pedestal sigma map to export to AMORE
2a7784f8	71	TH1F **f1DPedMean; //1D mean pedestal map to export to AMORE
2ac899f2	72	TH1F **f1DPedSigma; //1D pedestal sigma map to export to AMORE
2a7784f8	73	Int_t fNumMaskedPads; //Number of masked pads
	74	Int_t fNumDeadPads; //Number of currently dead channels
	75
2ac899f2	76	private:
	77
	78	AliHMPIDCalib(const AliHMPIDCalib& c); //dummy copy constructor
	79	AliHMPIDCalib &operator=(const AliHMPIDCalib& c); //dummy assignment operator
	80
2a7784f8	81	ClassDef(AliHMPIDCalib,5) //HMPID calibration and pedestal class
21f61e25	82	};
21f61e25	83	#endif
4969a7fa	84