-#ifndef AliTPCCALIBDBUTIL_H
-#define AliTPCCALIBDBUTIL_H
+#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 <TVectorD.h>
-#include <TVectorF.h>
+#include <TVectorDfwd.h>
+#include <TVectorFfwd.h>
+#include <TMatrixDfwd.h>
class TGraph;
class TMap;
class TGraphErrors;
class TTreeSRedirector;
class AliTPCCalROC;
+class AliSplineFit;
+
class AliTPCcalibDButil : public TObject
{
public:
void UpdateFromCalibDB();
//data processing functions
void ProcessCEdata(const char* fitFormula, TVectorD &fitResultsA, TVectorD &fitResultsC,
- Int_t &noutliersCE, Double_t & chi2A, Double_t &chi2C, AliTPCCalPad *outCE=0);
+ 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);
//getter preprocess information
Int_t GetNPulserOutliers() const {return fNpulserOutliers;}
- Float_t GetMeanAltro(const AliTPCCalROC *roc, const Int_t row, const Int_t pad, AliTPCCalROC *rocOut=0x0);
+ 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;}
- const Int_t GetReferenceRun(const char* type) const;
+ Int_t GetReferenceRun(const char* type) const;
const char* GetRefValidity() const {return fRefValidity.Data();}
AliTPCCalPad* GetRefPadNoise() const {return fRefPadNoise;}
//setters for pad by pad information
- void SetPulserData(AliTPCCalPad *tmean, AliTPCCalPad *trms=0x0, AliTPCCalPad *qmean=0x0)
+ void SetPulserData(AliTPCCalPad * const tmean, AliTPCCalPad * const trms=0x0, AliTPCCalPad * const qmean=0x0)
{fPulserTmean=tmean; fPulserTrms=trms; fPulserQmean=qmean;}
- void SetCEData(AliTPCCalPad *tmean, AliTPCCalPad *trms=0x0, AliTPCCalPad *qmean=0x0)
+ void SetCEData(AliTPCCalPad *const tmean, AliTPCCalPad *const trms=0x0, AliTPCCalPad *const qmean=0x0)
{fCETmean=tmean; fCETrms=trms; fCEQmean=qmean;}
- void SetNoisePedestal(AliTPCCalPad *noise, AliTPCCalPad *pedestal=0x0)
+ void SetNoisePedestal(AliTPCCalPad *const noise, AliTPCCalPad *const pedestal=0x0)
{fPadNoise=noise; fPedestals=pedestal;}
- void SetALTROData(AliTPCCalPad *masked)
+ void SetALTROData(AliTPCCalPad *const masked)
{fALTROMasked=masked;}
- void SetGoofieArray(AliDCSSensorArray *arr) {fGoofieArray=arr;}
+ 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 *tmean, AliTPCCalPad *trms=0x0, AliTPCCalPad *qmean=0x0)
+ void SetRefPulserData(AliTPCCalPad *const tmean, AliTPCCalPad *const trms=0x0, AliTPCCalPad *const qmean=0x0)
{fRefPulserTmean=tmean; fRefPulserTrms=trms; fRefPulserQmean=qmean;}
- void SetRefCEData(AliTPCCalPad *tmean, AliTPCCalPad *trms=0x0, AliTPCCalPad *qmean=0x0)
+ void SetRefCEData(AliTPCCalPad *const tmean, AliTPCCalPad *const trms=0x0, AliTPCCalPad *const qmean=0x0)
{fRefCETmean=tmean; fRefCETrms=trms; fRefCEQmean=qmean;}
- void SetRefNoisePedestal(AliTPCCalPad *noise, AliTPCCalPad *pedestal=0x0)
+ void SetRefNoisePedestal(AliTPCCalPad *const noise, AliTPCCalPad *const pedestal=0x0)
{fRefPadNoise=noise; fRefPedestals=pedestal;}
- void SetRefALTROData(AliTPCCalPad *masked)
+ void SetRefALTROData(AliTPCCalPad *const masked)
{fRefALTROMasked=masked;}
//creation of derived pad by pad calibration data
//
// create outlyer maps
//
- AliTPCCalPad *CreateCEOutlyerMap(Int_t &noutliersCE, AliTPCCalPad *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 *pulserOut=0, Float_t cutTime=3, Float_t cutnRMSQ=5, Float_t cutnRMSrms=5);
+ 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);
//
// graph tools
//
static Double_t GetLaserTime0(Int_t run, Int_t timeStamp, Int_t deltaT, Int_t side);
- static TGraph* FilterGraphMedian(TGraph * graph, Float_t sigmaCut, Double_t &medianY);
+ 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 *graph, Double_t xref);
+ static Double_t EvalGraphConst(TGraph * const graph, Double_t xref);
+ static Double_t EvalGraphConst(AliSplineFit *graph, Double_t xref);
+
//
// Filter sensors
//
//
// Filter AliRelAlignmentKalman - Alignment/Drift velocity
//
- static TMatrixD* MakeStatRelKalman(TObjArray *array, Float_t minFraction, Int_t minStat, Float_t maxvd);
- static TObjArray *SmoothRelKalman(TObjArray *array,TMatrixD & stat, Bool_t direction, Float_t sigmaCut);
- static TObjArray *SmoothRelKalman(TObjArray *arrayP, TObjArray *arrayM);
- static void FilterCE(Double_t deltaT=100, Double_t cutAbs=10, Double_t cutSigma=4., TTreeSRedirector *pcstream=0);
- static void FilterTracks(Int_t run, Double_t cutSigma=20., TTreeSRedirector *pcstream=0);
+ static TMatrixD* MakeStatRelKalman(TObjArray * const array, Float_t minFraction, Int_t minStat, Float_t maxvd);
+ static TObjArray *SmoothRelKalman(TObjArray * const array,const TMatrixD & stat, 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 *pcstream=0);
+ 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);
+
private:
AliTPCcalibDB *fCalibDB; //pointer to calibDB object
AliTPCCalPad *fPadNoise; //noise information
//
// filters
-public:
TArrayI fRuns; // run list with OCDB info
TArrayI fRunsStart; // start time for given run
TArrayI fRunsStop; // stop time for given run
-private:
+
AliTPCcalibDButil (const AliTPCcalibDButil& );
AliTPCcalibDButil& operator= (const AliTPCcalibDButil& );
- 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);
- const Int_t GetCurrentReferenceRun(const char* type);
- AliCDBEntry* GetRefEntry(const char* cdbPath);
+
+ static Double_t GetVDriftTPCLaserTracksCommon(Double_t &dist, Int_t timeStamp, Double_t deltaT, Int_t side, TObjArray * const array);
+
ClassDef(AliTPCcalibDButil,0)
};
git://git.uio.no / u / mrichter / AliRoot.git / blobdiff