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
    ~BadChannelAna();                // dtor
    BadChannelAna(TString period, TString train, TString trigger, Int_t runNumber,Int_t trial, TString workDir, TString listName, Bool_t runByRun=0);

    void Run(Bool_t mergeOnly=0);

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

protected:

    void Init();
    TString MergeRuns();
    void LoadExternalBadMap();
    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* BuildHitAndEnergyMeanScaled(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 FlagAsBad_Time(Int_t crit, TH1F* inhisto, Double_t nSig = 3);

    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();
    void PrintCellInfo(Int_t number);

    //Settings for analysed period
    Int_t   fCurrentRunNumber;            
      TString fPeriod;                      
    TString fTrainNo;                     
    TString fTrigger;                     
    Int_t   fNoOfCells;                   
    Int_t   fCellStartDCal;               
    Int_t   fStartCell;                   
    Int_t   fStartCellSM[21];             
    Double_t fEndLowerBound;              

    //Genergal paths
    TString fAnalysisOutput;              
    TString fAnalysisOutputRbR;           
    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;            
    TString fExternalBadMapName;          
    Bool_t  fTestRoutine;                 
      Bool_t  fPrint;                       
    Bool_t  fTrackCellRecord;             
    Bool_t  fRunBRunMap;                  

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

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

    //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, 2);
};
#endif