[u/mrichter/AliRoot.git] / TPC / TPCbase / AliTPCCalibViewer.h

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

/* $Id: AliTPCCalibViewer.h,v */

/// \class AliTPCCalibViewer
///
/// \brief TPC calibration viewer/visualization class
///  use Tree for visualization

class TFile;
class TLegend;
class TGraph;
class TH1;
class TH1F;
#include <TTree.h>
#include <TMatrixDfwd.h>
#include <TVectorDfwd.h>
#include <TVectorFfwd.h>

class AliTPCCalPad;
class AliTPCCalROC;
class TFriendElement;


class AliTPCCalibViewer : public TObject {
public:
   AliTPCCalibViewer();
   AliTPCCalibViewer(const AliTPCCalibViewer &c);
   AliTPCCalibViewer(TTree *const tree);
   AliTPCCalibViewer(const char* fileName, const char* treeName = "calPads");
   AliTPCCalibViewer &operator = (const AliTPCCalibViewer & param);
   virtual ~AliTPCCalibViewer();
   virtual void Delete(Option_t* option = "");
   
   TString& GetAbbreviation()  { return fAbbreviation;  }
   TString& GetAppendString()  { return fAppendString; }
   void SetAbbreviation(const Char_t *abr) { fAbbreviation = abr; }
   void SetAppendString(const Char_t *str) { fAppendString = str; }
   
   virtual void     Draw(Option_t* opt="") { fTree->Draw(opt); }
   virtual Long64_t Draw(const char* varexp, const TCut& selection, Option_t* option = "", Long64_t nentries = 1000000000, Long64_t firstentry = 0) { return fTree->Draw(varexp, selection, option, nentries, firstentry); };
   virtual Long64_t Draw(const char* varexp, const char* selection, Option_t* option = "", Long64_t nentries = 1000000000, Long64_t firstentry = 0) { return fTree->Draw(varexp, selection, option, nentries, firstentry); };

   const char* AddAbbreviations(const Char_t *c, Bool_t printDrawCommand = kFALSE);
   Int_t EasyDraw(const char* drawCommand, const char* sector, const char* cuts = 0, const char* drawOptions = 0, Bool_t writeDrawCommand = kFALSE) const;   // easy drawing of data, use '~' for abbreviation of '.fElements'
   Int_t EasyDraw(const char* drawCommand, Int_t sector, const char* cuts = 0, const char* drawOptions = 0, Bool_t writeDrawCommand = kFALSE) const;   // easy drawing of data, use '~' for abbreviation of '.fElements'
   Int_t EasyDraw1D(const char* drawCommand, const char* sector, const char* cuts = 0, const char* drawOptions = 0, Bool_t writeDrawCommand = kFALSE) const;   // easy drawing of data, use '~' for abbreviation of '.fElements'
   Int_t EasyDraw1D(const char* drawCommand, Int_t sector, const char* cuts = 0, const char* drawOptions = 0, Bool_t writeDrawCommand = kFALSE) const;   // easy drawing of data, use '~' for abbreviation of '.fElements'
   void FormatHistoLabels(TH1 *histo) const; // formats title and axis labels of histo, removes '.fElements'
   
   Int_t  DrawHisto1D(const char* drawCommand,       Int_t sector, const char* cuts = 0, const char *sigmas = "2;4;6", Bool_t plotMean = kTRUE, Bool_t plotMedian = kTRUE, Bool_t plotLTM = kTRUE) const; // draws 1d histograms and superimposes mean, median, ltm and several sigma cuts
   Int_t  DrawHisto1D(const char* drawCommand, const char* sector, const char* cuts = 0, const char *sigmas = "2;4;6", Bool_t plotMean = kTRUE, Bool_t plotMedian = kTRUE, Bool_t plotLTM = kTRUE) const; // draws 1d histograms and superimposes mean, median, ltm and several sigma cuts
   Int_t     SigmaCut(const char* drawCommand,       Int_t sector, const char* cuts = 0, Float_t sigmaMax = 5, Bool_t plotMean = kTRUE, Bool_t plotMedian = kTRUE, Bool_t plotLTM = kTRUE, Bool_t pm = kFALSE, const char *sigmas = "", Float_t sigmaStep = -1) const;    // draws fraction of used pads over different sigma cuts
   Int_t  SigmaCutNew(const char* drawCommand, const char* sector, const char* cuts = 0, Float_t sigmaMax = 5, Bool_t plotMean = kTRUE, Bool_t plotMedian = kTRUE, Bool_t plotLTM = kTRUE, Bool_t pm = kFALSE, const char *sigmas = "", Float_t sigmaStep = -1) const;    // draws fraction of used pads over different sigma cuts
   Int_t     SigmaCut(const char* drawCommand, const char* sector, const char* cuts = 0, Float_t sigmaMax = 5, Bool_t plotMean = kTRUE, Bool_t plotMedian = kTRUE, Bool_t plotLTM = kTRUE, Bool_t pm = kFALSE, const char *sigmas = "", Float_t sigmaStep = -1) const;    // draws fraction of used pads over different sigma cuts
   Int_t    Integrate(const char* drawCommand, const char* sector, const char* cuts = 0, Float_t sigmaMax = 5, Bool_t plotMean = kTRUE, Bool_t plotMedian = kTRUE, Bool_t plotLTM = kTRUE, const char *sigmas = "", Float_t sigmaStep = -1) const;    // draws an integrated histogram
   Int_t    Integrate(const char* drawCommand,       Int_t sector, const char* cuts = 0, Float_t sigmaMax = 5, Bool_t plotMean = kTRUE, Bool_t plotMedian = kTRUE, Bool_t plotLTM = kTRUE, const char *sigmas = "", Float_t sigmaStep = -1) const;    // draws an integrated histogram
   Int_t IntegrateOld(const char* drawCommand, const char* sector, const char* cuts = 0, Float_t sigmaMax = 5, Bool_t plotMean = kTRUE, Bool_t plotMedian = kTRUE, Bool_t plotLTM = kTRUE, const char *sigmas = "", Float_t sigmaStep = -1) const;    // draws an integrated histogram
   
   AliTPCCalPad* GetCalPadOld(const char* desiredData, const char* cuts = "", const char* calPadName = "NoName") const;     // returns an AliTPCCalPad object containing the specified data with cuts applied
   AliTPCCalPad* GetCalPad(const char* desiredData, const char* cuts = "", const char* calPadName = "NoName") const;     // returns an AliTPCCalPad object containing the specified data with cuts applied

   AliTPCCalROC* GetCalROC(const char* desiredData, UInt_t sector, const char* cuts = "") const;  // returns an AliTPCCalROC object containing the specified data for sector with cuts applied
   
   TObjArray* GetArrayOfCalPads();
   TObjArray* GetListOfVariables(Bool_t printList = kFALSE);
   TObjArray* GetListOfNormalizationVariables(Bool_t printList = kFALSE) const;
   
   static void MakeTreeWithObjects(const char * fileName, const TObjArray *const array, const char * mapFileName = 0);
   static void MakeTree(const char * fileName, TObjArray * array, const char * mapFileName = 0, AliTPCCalPad *const outlierPad = 0, Float_t ltmFraction = 0.9);
   static void MakeTree(const char *outPutFileName, const Char_t *inputFileName, AliTPCCalPad *outlierPad = 0, Float_t ltmFraction = 0.9, const char *mapFileName = "$ALICE_ROOT/TPC/Calib/MapCalibrationObjects.root");
   static void CreateObjectList(const Char_t *filename, TObjArray *calibObjects);
   
   TFriendElement* AddReferenceTree(const char* filename, const char* treename = "calPads", const char* refname = "R");
   TFriendElement* AddFriend(const char* treename, const char* filename) {return fTree->AddFriend(treename, filename);};
   TFriendElement* AddFriend(TTree* tree, const char* alias, Bool_t warn=kFALSE) {return fTree->AddFriend(tree, alias, warn);};
   TFriendElement* AddFriend(const char* treename, TFile* file) {return fTree->AddFriend(treename, file);};
   TTree * GetTree() const { return fTree;}
   static void MakeCalPadAliases(TTree * tree);
   TString* Fit(const char* drawCommand, const char* formula, const char* cuts, Double_t & chi2, TVectorD &fitParam, TMatrixD &covMatrix);

  // 
  // Array tools
  // 
  static Double_t GetLTM(Int_t n, const Double_t *const array, Double_t *const sigma = 0, Double_t fraction = 0.9);
  static Int_t GetBin(Float_t value, Int_t nbins, Double_t binLow, Double_t binUp);
  static TH1F* SigmaCut(Int_t n, const Float_t *array, Float_t mean, Float_t sigma, Int_t nbins, Float_t binLow, Float_t binUp, Float_t sigmaMax, Float_t sigmaStep = -1, Bool_t pm = kFALSE);
  static TH1F* SigmaCut(TH1F *const histogram, Float_t mean, Float_t sigma, Float_t sigmaMax, Float_t sigmaStep = -1, Bool_t pm = kFALSE);
  static TH1F* Integrate(TH1F *const histogram, Float_t mean = 0, Float_t sigma = 0, Float_t sigmaMax = 0, Float_t sigmaStep = -1);
  static TH1F* Integrate(Int_t n, const Float_t *const array, Int_t nbins, Float_t binLow, Float_t binUp, Float_t mean = 0, Float_t sigma = 0, Float_t sigmaMax = 0, Float_t sigmaStep = -1);
  
  static TH1F* SigmaCut(Int_t n, const Double_t *array, Double_t mean, Double_t sigma, Int_t nbins, const Double_t *xbins, Double_t sigmaMax);
   
   
protected:
   TTree* fTree;     ///< tree containing visualization data (e.g. written by AliTPCCalPad::MakeTree(...)
   TFile* fFile;     ///< file that contains a calPads tree (e.g. written by AliTPCCalPad::MakeTree(...)
   TObjArray* fListOfObjectsToBeDeleted;  ///< Objects, that will be deleted when the destructor ist called
   Bool_t fTreeMustBeDeleted;  ///< decides weather the tree must be deleted in destructor or not
   TString fAbbreviation;       ///< the abreviation for '.fElements'
   TString fAppendString;      ///< '.fElements', stored in a TStrig
   
   void DrawLines(TH1F *cutHistoMean, TVectorF nsigma, TLegend *legend, Int_t color, Bool_t pm) const;
   void DrawLines(TGraph *graph, TVectorF nsigma, TLegend *legend, Int_t color, Bool_t pm) const;
  
   
   /// \cond CLASSIMP
   ClassDef(AliTPCCalibViewer,1)    //  TPC calibration viewer class
   /// \endcond
};

#endif
Commit	Line	Data
39bcd65d	1	#ifndef ALITPCCALIBVIEWER_H
	2	#define ALITPCCALIBVIEWER_H
	3	/* Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved. *
	4	* See cxx source for full Copyright notice */
	5
	6	/* $Id: AliTPCCalibViewer.h,v */
	7
7d855b04	8	/// \class AliTPCCalibViewer
	9	///
	10	/// \brief TPC calibration viewer/visualization class
	11	/// use Tree for visualization
39bcd65d	12
2811495d	13	class TFile;
72d0ab7e	14	class TLegend;
4f3934a1	15	class TGraph;
661f340b	16	class TH1;
661f340b	17	class TH1F;
2811495d	18	#include <TTree.h>
2811495d	19	#include <TMatrixDfwd.h>
661f340b	20	#include <TVectorDfwd.h>
661f340b	21	#include <TVectorFfwd.h>
2811495d	22
	23	class AliTPCCalPad;
	24	class AliTPCCalROC;
	25	class TFriendElement;
72d0ab7e	26
39bcd65d	27
	28	class AliTPCCalibViewer : public TObject {
	29	public:
	30	AliTPCCalibViewer();
	31	AliTPCCalibViewer(const AliTPCCalibViewer &c);
2811495d	32	AliTPCCalibViewer(TTree *const tree);
a6e0ebfe	33	AliTPCCalibViewer(const char* fileName, const char* treeName = "calPads");
39bcd65d	34	AliTPCCalibViewer &operator = (const AliTPCCalibViewer & param);
39bcd65d	35	virtual ~AliTPCCalibViewer();
a6d2bd0c	36	virtual void Delete(Option_t* option = "");
39bcd65d	37
8192627b	38	TString& GetAbbreviation() { return fAbbreviation; }
8192627b	39	TString& GetAppendString() { return fAppendString; }
2811495d	40	void SetAbbreviation(const Char_t *abr) { fAbbreviation = abr; }
2811495d	41	void SetAppendString(const Char_t *str) { fAppendString = str; }
832bb12c	42
39bcd65d	43	virtual void Draw(Option_t* opt="") { fTree->Draw(opt); }
	44	virtual Long64_t Draw(const char* varexp, const TCut& selection, Option_t* option = "", Long64_t nentries = 1000000000, Long64_t firstentry = 0) { return fTree->Draw(varexp, selection, option, nentries, firstentry); };
	45	virtual Long64_t Draw(const char* varexp, const char* selection, Option_t* option = "", Long64_t nentries = 1000000000, Long64_t firstentry = 0) { return fTree->Draw(varexp, selection, option, nentries, firstentry); };
	46
2811495d	47	const char* AddAbbreviations(const Char_t *c, Bool_t printDrawCommand = kFALSE);
39bcd65d	48	Int_t EasyDraw(const char* drawCommand, const char* sector, const char* cuts = 0, const char* drawOptions = 0, Bool_t writeDrawCommand = kFALSE) const; // easy drawing of data, use '~' for abbreviation of '.fElements'
	49	Int_t EasyDraw(const char* drawCommand, Int_t sector, const char* cuts = 0, const char* drawOptions = 0, Bool_t writeDrawCommand = kFALSE) const; // easy drawing of data, use '~' for abbreviation of '.fElements'
	50	Int_t EasyDraw1D(const char* drawCommand, const char* sector, const char* cuts = 0, const char* drawOptions = 0, Bool_t writeDrawCommand = kFALSE) const; // easy drawing of data, use '~' for abbreviation of '.fElements'
	51	Int_t EasyDraw1D(const char* drawCommand, Int_t sector, const char* cuts = 0, const char* drawOptions = 0, Bool_t writeDrawCommand = kFALSE) const; // easy drawing of data, use '~' for abbreviation of '.fElements'
832bb12c	52	void FormatHistoLabels(TH1 *histo) const; // formats title and axis labels of histo, removes '.fElements'
72d0ab7e	53
	54	Int_t DrawHisto1D(const char* drawCommand, Int_t sector, const char* cuts = 0, const char *sigmas = "2;4;6", Bool_t plotMean = kTRUE, Bool_t plotMedian = kTRUE, Bool_t plotLTM = kTRUE) const; // draws 1d histograms and superimposes mean, median, ltm and several sigma cuts
	55	Int_t DrawHisto1D(const char* drawCommand, const char* sector, const char* cuts = 0, const char *sigmas = "2;4;6", Bool_t plotMean = kTRUE, Bool_t plotMedian = kTRUE, Bool_t plotLTM = kTRUE) const; // draws 1d histograms and superimposes mean, median, ltm and several sigma cuts
	56	Int_t SigmaCut(const char* drawCommand, Int_t sector, const char* cuts = 0, Float_t sigmaMax = 5, Bool_t plotMean = kTRUE, Bool_t plotMedian = kTRUE, Bool_t plotLTM = kTRUE, Bool_t pm = kFALSE, const char *sigmas = "", Float_t sigmaStep = -1) const; // draws fraction of used pads over different sigma cuts
4f3934a1	57	Int_t SigmaCutNew(const char* drawCommand, const char* sector, const char* cuts = 0, Float_t sigmaMax = 5, Bool_t plotMean = kTRUE, Bool_t plotMedian = kTRUE, Bool_t plotLTM = kTRUE, Bool_t pm = kFALSE, const char *sigmas = "", Float_t sigmaStep = -1) const; // draws fraction of used pads over different sigma cuts
72d0ab7e	58	Int_t SigmaCut(const char* drawCommand, const char* sector, const char* cuts = 0, Float_t sigmaMax = 5, Bool_t plotMean = kTRUE, Bool_t plotMedian = kTRUE, Bool_t plotLTM = kTRUE, Bool_t pm = kFALSE, const char *sigmas = "", Float_t sigmaStep = -1) const; // draws fraction of used pads over different sigma cuts
	59	Int_t Integrate(const char* drawCommand, const char* sector, const char* cuts = 0, Float_t sigmaMax = 5, Bool_t plotMean = kTRUE, Bool_t plotMedian = kTRUE, Bool_t plotLTM = kTRUE, const char *sigmas = "", Float_t sigmaStep = -1) const; // draws an integrated histogram
	60	Int_t Integrate(const char* drawCommand, Int_t sector, const char* cuts = 0, Float_t sigmaMax = 5, Bool_t plotMean = kTRUE, Bool_t plotMedian = kTRUE, Bool_t plotLTM = kTRUE, const char *sigmas = "", Float_t sigmaStep = -1) const; // draws an integrated histogram
4f3934a1	61	Int_t IntegrateOld(const char* drawCommand, const char* sector, const char* cuts = 0, Float_t sigmaMax = 5, Bool_t plotMean = kTRUE, Bool_t plotMedian = kTRUE, Bool_t plotLTM = kTRUE, const char *sigmas = "", Float_t sigmaStep = -1) const; // draws an integrated histogram
39bcd65d	62
a6e0ebfe	63	AliTPCCalPad* GetCalPadOld(const char* desiredData, const char* cuts = "", const char* calPadName = "NoName") const; // returns an AliTPCCalPad object containing the specified data with cuts applied
a6e0ebfe	64	AliTPCCalPad* GetCalPad(const char* desiredData, const char* cuts = "", const char* calPadName = "NoName") const; // returns an AliTPCCalPad object containing the specified data with cuts applied
ac7becc3	65
f78da5ae	66	AliTPCCalROC* GetCalROC(const char* desiredData, UInt_t sector, const char* cuts = "") const; // returns an AliTPCCalROC object containing the specified data for sector with cuts applied
39bcd65d	67
	68	TObjArray* GetArrayOfCalPads();
	69	TObjArray* GetListOfVariables(Bool_t printList = kFALSE);
72d0ab7e	70	TObjArray* GetListOfNormalizationVariables(Bool_t printList = kFALSE) const;
39bcd65d	71
2811495d	72	static void MakeTreeWithObjects(const char * fileName, const TObjArray const array, const char mapFileName = 0);
2811495d	73	static void MakeTree(const char * fileName, TObjArray * array, const char * mapFileName = 0, AliTPCCalPad *const outlierPad = 0, Float_t ltmFraction = 0.9);
4f3934a1	74	static void MakeTree(const char outPutFileName, const Char_t inputFileName, AliTPCCalPad outlierPad = 0, Float_t ltmFraction = 0.9, const char mapFileName = "$ALICE_ROOT/TPC/Calib/MapCalibrationObjects.root");
4f3934a1	75	static void CreateObjectList(const Char_t filename, TObjArray calibObjects);
39bcd65d	76
	77	TFriendElement* AddReferenceTree(const char* filename, const char* treename = "calPads", const char* refname = "R");
	78	TFriendElement* AddFriend(const char* treename, const char* filename) {return fTree->AddFriend(treename, filename);};
	79	TFriendElement* AddFriend(TTree* tree, const char* alias, Bool_t warn=kFALSE) {return fTree->AddFriend(tree, alias, warn);};
	80	TFriendElement* AddFriend(const char* treename, TFile* file) {return fTree->AddFriend(treename, file);};
2811495d	81	TTree * GetTree() const { return fTree;}
6d49cfb3	82	static void MakeCalPadAliases(TTree * tree);
a6d2bd0c	83	TString* Fit(const char* drawCommand, const char* formula, const char* cuts, Double_t & chi2, TVectorD &fitParam, TMatrixD &covMatrix);
72d0ab7e	84
	85	//
	86	// Array tools
	87	//
2811495d	88	static Double_t GetLTM(Int_t n, const Double_t const array, Double_t const sigma = 0, Double_t fraction = 0.9);
72d0ab7e	89	static Int_t GetBin(Float_t value, Int_t nbins, Double_t binLow, Double_t binUp);
2811495d	90	static TH1F* SigmaCut(Int_t n, const Float_t *array, Float_t mean, Float_t sigma, Int_t nbins, Float_t binLow, Float_t binUp, Float_t sigmaMax, Float_t sigmaStep = -1, Bool_t pm = kFALSE);
	91	static TH1F* SigmaCut(TH1F *const histogram, Float_t mean, Float_t sigma, Float_t sigmaMax, Float_t sigmaStep = -1, Bool_t pm = kFALSE);
	92	static TH1F* Integrate(TH1F *const histogram, Float_t mean = 0, Float_t sigma = 0, Float_t sigmaMax = 0, Float_t sigmaStep = -1);
	93	static TH1F* Integrate(Int_t n, const Float_t *const array, Int_t nbins, Float_t binLow, Float_t binUp, Float_t mean = 0, Float_t sigma = 0, Float_t sigmaMax = 0, Float_t sigmaStep = -1);
72d0ab7e	94
2811495d	95	static TH1F* SigmaCut(Int_t n, const Double_t array, Double_t mean, Double_t sigma, Int_t nbins, const Double_t xbins, Double_t sigmaMax);
72d0ab7e	96
	97
	98
39bcd65d	99	protected:
7d855b04	100	TTree* fTree; ///< tree containing visualization data (e.g. written by AliTPCCalPad::MakeTree(...)
	101	TFile* fFile; ///< file that contains a calPads tree (e.g. written by AliTPCCalPad::MakeTree(...)
	102	TObjArray* fListOfObjectsToBeDeleted; ///< Objects, that will be deleted when the destructor ist called
	103	Bool_t fTreeMustBeDeleted; ///< decides weather the tree must be deleted in destructor or not
	104	TString fAbbreviation; ///< the abreviation for '.fElements'
	105	TString fAppendString; ///< '.fElements', stored in a TStrig
39bcd65d	106
72d0ab7e	107	void DrawLines(TH1F cutHistoMean, TVectorF nsigma, TLegend legend, Int_t color, Bool_t pm) const;
4f3934a1	108	void DrawLines(TGraph graph, TVectorF nsigma, TLegend legend, Int_t color, Bool_t pm) const;
72d0ab7e	109
72d0ab7e	110
7d855b04	111	/// \cond CLASSIMP
39bcd65d	112	ClassDef(AliTPCCalibViewer,1) // TPC calibration viewer class
7d855b04	113	/// \endcond
39bcd65d	114	};
	115
	116	#endif
4f3934a1	117
	118
	119