BadChannelAna.h

Go to the documentation of this file.

#ifndef BADCHANNELANA_H
#define BADCHANNELANA_H



/* Copyright(c) 1998-2016, ALICE Experiment at CERN, All rights reserved. *
 * See cxx source for full Copyright notice
*/
 
 
#include <Riostream.h>
#include <TString.h>
#include <TArrayD.h>

class TH1;
class TH1F;
class TH1D;
class TH2F;
class TFile;
class TList;

class AliCalorimeterUtils;


class BadChannelAna : public TObject {
  
public:
      BadChannelAna() ;                // default ctor
    virtual ~BadChannelAna()  { ; }  // virtual dtor
    BadChannelAna(TString period, TString train, TString trigger, Int_t runNumber,Int_t trial, TString workDir, TString listName);

    void Run(Bool_t mergeOnly=0);

      //Setters
    void SetExternalMergedFile(TString inputName)        {fExternalFileName = inputName;}
      void SetQAChecks(Bool_t inputBool)                   {fTestRoutine      = inputBool;}
      void SetPrintOutput(Bool_t inputBool)                {fPrint            = inputBool;}
      void SetStartEndCell(Int_t start, Int_t end)         {fStartCell = start; fNoOfCells = end;}
      void AddManualMasking(std::vector<Int_t> cellVector) {fManualMask.swap(cellVector) ;}
    void AddPeriodAnalysis(Int_t criteria, Double_t nsigma, Double_t emin, Double_t emax);


protected:

    void Init();
    TString MergeRuns();
    void BCAnalysis();
    void PeriodAnalysis(Int_t criterum=7, Double_t nsigma = 4.0, Double_t emin=0.1, Double_t emax=2.0);

    TH1F* BuildHitAndEnergyMean(Int_t crit, Double_t emin = 0.1, Double_t emax=2.);
    TH1F* BuildTimeMean(Int_t crit, Double_t tmin, Double_t tmax);

    void FlagAsDead();
    void FlagAsBad(Int_t crit, TH1F* inhisto, Double_t nsigma = 4., Double_t dnbins = 200);

    void SummarizeResultsByFlag();
    void SummarizeResults();
    TH1D *BuildMeanFromGood(Int_t warmIn=0);
    Bool_t CheckDistribution(TH1* ratio, TH1* reference);
    Bool_t IsCoveredByTRD(Int_t row, Int_t collumn);
    void SaveBadCellsToPDF(Int_t version, TString pdfName);
    void PlotFlaggedCells2D(Int_t flagBegin,Int_t flagEnd=-1);
    void SaveHistoToFile();

    //Test Test
    //Settings for analysed period
    Int_t   fCurrentRunNumber;            
      TString fPeriod;                      
    TString fTrainNo;                     
    TString fTrigger;                     
    Int_t   fNoOfCells;                   
    Int_t   fCellStartDCal;               
    Int_t   fStartCell;                   

    //Genergal paths
    TString fAnalysisOutput;              
    TString fAnalysisInput;               
    TString fRunList;                     
    TString fRunListFileName;             
    TString fWorkdir;                     

    //
    TString fQADirect;                    
    TString fMergedFileName;              
    std::vector<TArrayD> fAnalysisVector; 

    //Things to be individualized by setters
    Int_t   fTrial;                       
    TString fExternalFileName;            
      Bool_t  fTestRoutine;                 
      Bool_t  fPrint;                       

    //histogram settings
    Int_t fNMaxCols;                      
    Int_t fNMaxRows;                      
    Int_t fNMaxColsAbs;                   
    Int_t fNMaxRowsAbs;                   

    //arrays to store information
    Double_t fnEventsInRange;
    Int_t *fFlag;                         
    Int_t fCriterionCounter;              
    Bool_t *fWarmCell;                    
    std::vector<Int_t> fManualMask;       

    //Calorimeter information for the investigated runs
    AliCalorimeterUtils* fCaloUtils;      

    TFile* fRootFile;                     
    TH2F* fCellAmplitude;                 
    TH2F* fCellTime;                      
    TH1F* fProcessedEvents;               
    TH1F* fhCellFlag;                     
    TH1F* fhCellWarm;                     

    TList* fOutputListBad;                
    TList* fOutputListBadRatio;           
    TList* fOutputListGood;               
    TList* fOutputListGoodRatio;          
    
private:
    BadChannelAna           (const BadChannelAna&); // not implemented
    BadChannelAna &operator=(const BadChannelAna&); // not implemented

    TH1F* fAvgNHitPerEvVsCellId;          
    TH1F* fAvgEngPerHitVsCellId;          
    
  
  ClassDef(BadChannelAna, 1);
};
#endif