]> git.uio.no Git - u/mrichter/AliRoot.git/blobdiff - HMPID/AliHMPIDCalib.h
Adding more bins in QA (Alis)
[u/mrichter/AliRoot.git] / HMPID / AliHMPIDCalib.h
index 1354d4084dc9887079e6a3eefcac7fe68b7b40ba..937fe1c4fa69be3f76741bea3d6948a7b37f19ad 100644 (file)
@@ -7,10 +7,20 @@
 //.
 //.
 //.
-#include "TTreePlayer.h"
-#include <TTree.h>
+
+//#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 { 
 
 
@@ -18,19 +28,57 @@ public:
   AliHMPIDCalib();
   virtual ~AliHMPIDCalib();
           void Init();
-          void FillPedestal(Int_t ddl,Int_t row, Int_t dil,Int_t adr,Int_t q);
-        Bool_t CalcPedestal(Int_t ddl, Char_t* name);
-
+          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: 
-    TTree   *fPedTree;                                                            //Pedestal Tree
-    Int_t   fa;                                                                   //DILOGIC address
-    Int_t   fd;                                                                   //DILOGIC number
-    Int_t   fr;                                                                   //DILOGIC row
-    Bool_t  faddl[11];                                                               //check is ddl is filled
-    Int_t   fq;                                                                   //Qdc value
-    TH1F   *fPedHisto;                                                            //                                                      
-    ClassDef(AliHMPIDCalib,1)                                                     //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
+