//#include <TTree.h>
#include <TH1.h>
#include <TH2.h>
-#include <TH1S.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 {
AliHMPIDCalib();
virtual ~AliHMPIDCalib();
void Init();
-// void FillPedestal(Int_t nDDL,Int_t row, Int_t dil,Int_t adr,Int_t q);
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);
- Bool_t CalcPedestal(ULong_t runNum, Int_t nDDL, Char_t* name, Int_t nEv); //number of the DDL, name of the output file and the number of events processed
- Bool_t WriteErrors(ULong_t runNum, Int_t nDDL, Char_t* name, Int_t nEv); //number of the DDL, name of the output file and the number of events processed
-
- enum {
- kNRows = 24, // Number of rows (starting from 1 !)//was25
- kNDILOGICAdd = 10, // Number of DILOGIC addresses in a row (starting from 1 !) //was11
- kNPadAdd = 48, // Number of pad row
- kNDDL = 14
- };
+ 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[11]; //check is ddl is filled
- Float_t fsq[kNDDL][kNRows][kNDILOGICAdd][kNPadAdd]; //Sum of pad Q
- Float_t fsq2[kNDDL][kNRows][kNDILOGICAdd][kNPadAdd]; //Sum of pad Q^2
- Int_t fNumOfErr[kNDDL][AliHMPIDRawStream::kSumErr]; // Store the numner of errors for a given error type and a given DDL
- ClassDef(AliHMPIDCalib,2) //HMPID calibration and pedestal class
+
+ 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
+