AliTPCcalibDButil.h

Go to the documentation of this file.

#ifndef ALITPCCALIBDBUTIL_H
#define ALITPCCALIBDBUTIL_H

/* Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved. *
 * See cxx source for full Copyright notice                               */


#include <TObject.h>
#include <TArrayI.h>
#include <TVectorDfwd.h>
#include <TVectorFfwd.h>
#include <TMatrixDfwd.h>

class TGraph;
class TMap;
class AliDCSSensorArray;
class AliTPCcalibDB;
class AliTPCCalPad;
class AliTPCCalROC;
class AliTPCmapper;
class AliTPCCalibRaw;
class AliCDBEntry;
class AliDCSSensor;
class AliDCSSensorArray;
class AliTPCSensorTempArray;
class AliTPCdataQA;
class TGraphErrors;
class TTreeSRedirector;
class AliTPCCalROC;
class AliSplineFit;

class AliTPCcalibDButil : public TObject
{
public:
  AliTPCcalibDButil();
  virtual ~AliTPCcalibDButil();

  void UpdateFromCalibDB();
  //data processing functions
  void ProcessCEdata(const char* fitFormula, TVectorD &fitResultsA, TVectorD &fitResultsC,
                     Int_t &noutliersCE, Double_t & chi2A, Double_t &chi2C, AliTPCCalPad * const outCE=0);
  void ProcessCEgraphs(TVectorD &vecTEntries, TVectorD &vecTMean, TVectorD &vecTRMS, TVectorD &vecTMedian,
                       TVectorD &vecQEntries, TVectorD &vecQMean, TVectorD &vecQRMS, TVectorD &vecQMedian,
                       Float_t &driftTimeA, Float_t &driftTimeC );
  void ProcessNoiseData(TVectorD &vNoiseMean, TVectorD &vNoiseMeanSenRegions,
                        TVectorD &vNoiseRMS, TVectorD &vNoiseRMSSenRegions,
                        Int_t &nonMaskedZero, Int_t &nNaN);
  void ProcessQAData(TVectorD &vQaOcc, TVectorD &vQaQtot, TVectorD &vQaQmax);
  void ProcessPulser(TVectorD &vMeanTime);
  void ProcessALTROConfig(Int_t &nMasked);
  void ProcessGoofie(TVectorD & vecEntries, TVectorD & vecMedian, TVectorD &vecMean, TVectorD &vecRMS);
  
  //processing functions using reference data
  void ProcessPedestalVariations(TVectorF &pedestalDeviations);
  void ProcessNoiseVariations(TVectorF &noiseDeviations);
  void ProcessPulserVariations(TVectorF &pulserQdeviations, Float_t &varQMean, Int_t &npadsOutOneTB, Int_t &npadsOffAdd);
  
  //getter preprocess information
  Int_t GetNPulserOutliers() const {return fNpulserOutliers;}
  Float_t GetMeanAltro(const AliTPCCalROC *roc, const Int_t row, const Int_t pad, AliTPCCalROC * const rocOut=0x0);
  AliTPCCalPad *GetPulserOutlierMap() const {return fPulserOutlier;}

  //getters ref data
  TMap *GetReferenceMap() const {return fRefMap;}
  Int_t GetReferenceRun(const char* type) const;
  const char* GetRefValidity() const {return fRefValidity.Data();}
  
  AliTPCCalPad* GetRefPadNoise() const {return fRefPadNoise;}
  AliTPCCalPad* GetRefPedestals() const {return fRefPedestals;}
  AliTPCCalPad* GetRefPedestalMasked() const {return fRefPedestalMasked;}
  AliTPCCalPad* GetRefPulserTmean() const {return fRefPulserTmean;}
  AliTPCCalPad* GetRefPulserTrms() const {return fRefPulserTrms;}
  AliTPCCalPad* GetRefPulserQmean() const {return fRefPulserQmean;}
  AliTPCCalPad* GetRefPulserOutlier() const {return fRefPulserOutlier;}
  AliTPCCalPad* GetRefPulserMasked() const {return fRefPulserMasked;}
  AliTPCCalPad* GetRefCETmean() const {return fRefCETmean;}
  AliTPCCalPad* GetRefCETrms() const {return fRefCETrms;}
  AliTPCCalPad* GetRefCEQmean() const {return fRefCEQmean;}
  AliTPCCalPad* GetRefCEMasked() const {return fRefCEMasked;}
  AliTPCCalPad* GetRefALTROFPED() const {return fRefALTROFPED;}
  AliTPCCalPad* GetRefALTROZsThr() const {return fRefALTROZsThr;}
  AliTPCCalPad* GetRefALTROAcqStart() const {return fRefALTROAcqStart;}
  AliTPCCalPad* GetRefALTROAcqStop() const {return fRefALTROAcqStop;}
  AliTPCCalPad* GetRefALTROMasked() const {return fRefALTROMasked;}
  
  
  //setters for pad by pad information
  void SetPulserData(AliTPCCalPad * const tmean, AliTPCCalPad * const trms=0x0, AliTPCCalPad * const qmean=0x0)
                {fPulserTmean=tmean; fPulserTrms=trms; fPulserQmean=qmean;}
  void SetCEData(AliTPCCalPad *const tmean, AliTPCCalPad *const trms=0x0, AliTPCCalPad *const qmean=0x0)
                {fCETmean=tmean; fCETrms=trms; fCEQmean=qmean;}
  void SetNoisePedestal(AliTPCCalPad *const noise, AliTPCCalPad *const pedestal=0x0)
                {fPadNoise=noise; fPedestals=pedestal;}
  void SetALTROData(AliTPCCalPad *const masked)
                {fALTROMasked=masked;}
  void SetGoofieArray(AliDCSSensorArray *const arr) {fGoofieArray=arr;}
  
  //setters for pad by pad information
  void SetRefFile(const char* filename);
  void SetReferenceRun(Int_t run=-1);
  void UpdateRefDataFromOCDB();
  void SetRefPulserData(AliTPCCalPad *const tmean, AliTPCCalPad *const trms=0x0, AliTPCCalPad *const qmean=0x0)
                {fRefPulserTmean=tmean; fRefPulserTrms=trms; fRefPulserQmean=qmean;}
  void SetRefCEData(AliTPCCalPad *const tmean, AliTPCCalPad *const trms=0x0, AliTPCCalPad *const qmean=0x0)
                {fRefCETmean=tmean; fRefCETrms=trms; fRefCEQmean=qmean;}
  void SetRefNoisePedestal(AliTPCCalPad *const noise, AliTPCCalPad *const pedestal=0x0)
                {fRefPadNoise=noise; fRefPedestals=pedestal;}
  void SetRefALTROData(AliTPCCalPad *const masked)
                {fRefALTROMasked=masked;}
  
  //creation of derived pad by pad calibration data
  AliTPCCalPad *CreatePadTime0(Int_t model, Double_t &gyA, Double_t &gyC, Double_t &chi2A, Double_t &chi2C);
  //
  // create outlyer maps
  //
  AliTPCCalPad *CreateCEOutlyerMap(Int_t &noutliersCE, AliTPCCalPad * const ceOut=0, Float_t minSignal=10, Float_t cutTrmsMin=0.9, Float_t cutTrmsMax=1.2, Float_t cutMaxDistT=0.7);
  AliTPCCalPad *CreatePulserOutlyerMap(Int_t &noutliersPulser, AliTPCCalPad * const pulserOut=0, Float_t cutTime=3, Float_t cutnRMSQ=5, Float_t cutnRMSrms=5);
  //
  AliTPCCalPad *CreatePadTime0CE(TVectorD &fitResultsA, TVectorD&fitResultsC, Int_t &nOut, Double_t &chi2A, Double_t &chi2C, const char *dumpfile=0);
  //

  void UpdatePulserOutlierMap();
  void UpdateRefPulserOutlierMap();
  void PulserOutlierMap(AliTPCCalPad *pulOut, const AliTPCCalPad *pulT, const AliTPCCalPad *pulQ);

  const char* GetGUIRefTreeDefaultName();
  
  Bool_t CreateGUIRefTree(const char* filename="");
  //
  // graph tools
  //
  static Double_t GetLaserTime0(Int_t run, Int_t timeStamp, Int_t deltaT, Int_t side);
  static TGraph* FilterGraphMedian(TGraph * const graph, Float_t sigmaCut, Double_t &medianY);
  static TGraph* FilterGraphMedianAbs(TGraph * graph, Float_t cut, Double_t &medianY);
  static TGraphErrors* FilterGraphMedianErr(TGraphErrors * graph, Float_t sigmaCut,Double_t &medianY);
  //
  static void Sort(TGraph *graph);
  static void SmoothGraph(TGraph *graph, Double_t delta);
  static Int_t     GetNearest(TGraph *graph, Double_t xref, Double_t &dx, Double_t &y);
  static Double_t EvalGraphConst(TGraph * const graph, Double_t xref);
  static Double_t EvalGraphConst(AliSplineFit *graph, Double_t xref);

  //
  // Filter sensors
  //
  static Float_t FilterSensor(AliDCSSensor * sensor, Double_t ymin, Double_t ymax, Double_t maxdy, Double_t sigmaCut); 
  //
  // Filter AliRelAlignmentKalman - Alignment/Drift velocity
  //
  static TMatrixD* MakeStatRelKalman(TObjArray * const array, Float_t minFraction, Int_t minStat, Float_t maxvd);
  static TObjArray *SmoothRelKalman(TObjArray * const array,const TMatrixD * statP, Bool_t direction, Float_t sigmaCut);
  static TObjArray *SmoothRelKalman(TObjArray * const arrayP, TObjArray * const arrayM);
  static void FilterCE(Double_t deltaT=100, Double_t cutAbs=10, Double_t cutSigma=4., TTreeSRedirector * const pcstream=0);
  static void FilterTracks(Int_t run, Double_t cutSigma=20., TTreeSRedirector * const pcstream=0);
  static Float_t FilterTemperature(AliTPCSensorTempArray *tempArray, Double_t ymin=15, Double_t ymax=22, Double_t sigmaCut=5); 

  static   void FilterGoofie(AliDCSSensorArray * goofieArray, Double_t deltaT=2, Double_t cutSigma=4.,  Double_t minVdn=8.5, Double_t maxVdn=9.05, TTreeSRedirector * const pcstream=0);
  static Double_t  GetTriggerOffsetTPC(Int_t run, Int_t timeStamp, Double_t deltaT=86400, Double_t deltaTLaser=3600, Int_t valType=0);
  static Double_t  GetVDriftTPC(Double_t &dist, Int_t run, Int_t timeStamp, Double_t deltaT=86400, Double_t deltaTLaser=3600, Int_t valType=0);
  static Double_t  GetVDriftTPCLaserTracks(Double_t &dist,Int_t run, Int_t timeStamp, Double_t deltaT=43200, Int_t side=2);
  static Double_t  GetVDriftTPCLaserTracksOnline(Double_t &dist, Int_t /*run*/, Int_t timeStamp, Double_t deltaT=43200, Int_t side=2);
  static Double_t  GetVDriftTPCCE(Double_t &dist, Int_t run, Int_t timeStamp, Double_t deltaT=43200, Int_t side=2);
  static Double_t  GetVDriftTPCITS(Double_t &dist, Int_t run, Int_t timeStamp);
  static Double_t  GetTime0TPCITS(Double_t &dist, Int_t run, Int_t timeStamp);
  Int_t MakeRunList(Int_t startRun, Int_t stopRun); // find the list of usable runs
  Int_t FindRunTPC(Int_t    itime, Bool_t debug=kFALSE);

 AliTPCCalPad* GetRefCalPad(AliCDBEntry *entry, const char* objName);
  AliTPCCalPad* GetRefCalPad(AliCDBEntry *entry);
  AliTPCCalPad* GetAltroMasked(const char* cdbPath, const char* name);
  Bool_t HasRefChanged(const char *cdbPath);
  Int_t GetCurrentReferenceRun(const char* type) const;
  AliCDBEntry* GetRefEntry(const char* cdbPath);

  static TTree* ConnectGainTrees(TString baseDir);
  static TTree* ConnectPulserTrees(TString baseDir, TTree * tree=0);
  //static TTree* ConnectNoiseTrees(TString baseDir,  TTree * tree=0); 
  static TTree* ConnectDistortionTrees(TString baseDir, TString  selection, TTree *tMain);
  static  TTree* ConnectCalPadTrees(TString baseDir, TString pattern,  TTree *tMain, Bool_t checkAliases);

private:
  AliTPCcalibDB *fCalibDB;            
  AliTPCCalPad  *fPadNoise;           
  AliTPCCalPad  *fPedestals;          
  AliTPCCalPad  *fPulserTmean;        
  AliTPCCalPad  *fPulserTrms;         
  AliTPCCalPad  *fPulserQmean;        
  AliTPCCalPad  *fPulserOutlier;      
  AliTPCCalPad  *fCETmean;            
  AliTPCCalPad  *fCETrms;             
  AliTPCCalPad  *fCEQmean;            
  AliTPCCalPad  *fALTROMasked;        
  //
  AliTPCCalibRaw *fCalibRaw;          
  //
  AliTPCdataQA   *fDataQA;            
  //reference data
  TMap *fRefMap;                        
  TMap *fCurrentRefMap;                 
  TString fRefValidity;                 
  //  
  AliTPCCalPad  *fRefPadNoise;           
  AliTPCCalPad  *fRefPedestals;          
  AliTPCCalPad  *fRefPedestalMasked;     
  AliTPCCalPad  *fRefPulserTmean;        
  AliTPCCalPad  *fRefPulserTrms;         
  AliTPCCalPad  *fRefPulserQmean;        
  AliTPCCalPad  *fRefPulserOutlier;      
  AliTPCCalPad  *fRefPulserMasked;       
  AliTPCCalPad  *fRefCETmean;            
  AliTPCCalPad  *fRefCETrms;             
  AliTPCCalPad  *fRefCEQmean;            
  AliTPCCalPad  *fRefCEMasked;           
  AliTPCCalPad  *fRefALTROFPED;          
  AliTPCCalPad  *fRefALTROZsThr;         
  AliTPCCalPad  *fRefALTROAcqStart;      
  AliTPCCalPad  *fRefALTROAcqStop;       
  AliTPCCalPad  *fRefALTROMasked;        
  //
  AliTPCCalibRaw *fRefCalibRaw;          
  //
  AliTPCdataQA   *fRefDataQA;            
  //
  AliDCSSensorArray* fGoofieArray;    
  //
  AliTPCmapper  *fMapper;             
  Int_t fNpulserOutliers;             
  
  Float_t fIrocTimeOffset;               
  Float_t fCETmaxLimitAbs;               
  Float_t fPulTmaxLimitAbs;              
  Float_t fPulQmaxLimitAbs;              
  Float_t fPulQminLimit;                 

  //
  // helpers to get the run number for given time stamps
  //
  // filters  
  
  TArrayI fRuns;                         
  TArrayI fRunsStart;                    
  TArrayI fRunsStop;                     
  
  AliTPCcalibDButil (const AliTPCcalibDButil& );
  AliTPCcalibDButil& operator= (const AliTPCcalibDButil& );

  
  static Double_t  GetVDriftTPCLaserTracksCommon(Double_t &dist, Int_t timeStamp, Double_t deltaT, Int_t side, TObjArray * const array);
  
  
  ClassDef(AliTPCcalibDButil,0)
};


#endif