//#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 TH1I;
class TFile;
class AliHMPIDCalib: public TObject {
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
- Bool_t CalcPedestal(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 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 InitHisto(Int_t q,Int_t histocnt,Char_t* name); //Init the pad histograms
- void FillHisto(Int_t histocnt,Int_t q); //Fill the ADC histograms
- void InitFile(Int_t nDDL); //Init the ADC histo output file (one per LDC)
- void CloseFile(Int_t nDDL); //Close the file
void SetRunParams(ULong_t runNum,Int_t timeStamp, Int_t ldcId); //Set Run Parameters such as Run Number, TimeStamp, LDCid
- inline void SetSigCut(Int_t nSigCut) { fSigCut=nSigCut;} //Set Sigma Cuts from Setter
- void SetSigCutFromFile(Char_t* name); //Set Sigma Cuts from File
- inline void SetWriteHistoPads(Bool_t isOn) {fWritePads=isOn;} //Set wether ADC histos of pads are written or not
- inline Bool_t GetWritePads() {return fWritePads;} //Set wether ADC histos of pads are written or not
+ 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[AliHMPIDRawStream::kNDDL+1][AliHMPIDRawStream::kNRows+1][AliHMPIDRawStream::kNDILOGICAdd+1][AliHMPIDRawStream::kNPadAdd+1]; //Sum of pad Q
- Float_t fsq2[AliHMPIDRawStream::kNDDL+1][AliHMPIDRawStream::kNRows+1][AliHMPIDRawStream::kNDILOGICAdd+1][AliHMPIDRawStream::kNPadAdd+1]; //Sum of pad Q^2
- Int_t fErr[AliHMPIDRawStream::kNDDL+1][AliHMPIDRawStream::kSumErr+1]; // 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
- Bool_t fWritePads; //Select wether to write ADC pad histograms or not
- 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
+
git://git.uio.no / u / mrichter / AliRoot.git / blobdiff