HafHistogram/interface/HHistogram.h

#ifndef HHISTOGRAM_H
#define HHISTOGRAM_H

/** \class HHistogram  Ganzhur/HafHistogram/interface/HHistogram.h
 *
 * Description:
 * Analysis code of the CMS experiment;
 * Original version is a part of HAF package developed for CMS: UserCode/HAF  
 * Histogram class
 *
 * \author Marcel Kunze,        Ruhr-University Bochum, Germany 
 * \author Serguei Ganjour,     CEA-Saclay/IRFU, FR
 *
 * \version $Id: HHistogram.h,v 1.2 2009/11/19 09:40:31 ganzhur Exp $
 *
 */

#include "TH1.h"
#include "TH2.h"

class HCandList;
class HCandidate;

class TObject;
class TH1;
class HHistID;

class HHistogram : public TNamed { 
    
public:
    
    // Constructors:
    // 1-D:
    HHistogram( const char*, const char*, Int_t, Axis_t, Axis_t );
    // 2-D:
    HHistogram( const char*, const char*, Int_t, Axis_t, Axis_t, Int_t, Axis_t, Axis_t );
    // Profile:
    HHistogram( const char*, const char*, Int_t, Axis_t, Axis_t, Axis_t, Axis_t );
    
    // Destructor:
    virtual ~HHistogram();
    
    // Functions to fill a histo:
    void Accumulate( Axis_t, Stat_t weight= 1.0 );
    void Accumulate1( Axis_t, Stat_t weight= 1.0 );
    void Accumulate( Axis_t x, Axis_t y, Stat_t weight );
    void Accumulate2( Axis_t, Axis_t, Stat_t weight= 1.0 );
    
    // clear all contents
    void Reset();
    
    // accessors:
    const char* Title() const { return GetTitle(); };
    Float_t GetContents( Int_t, Int_t nbinsY= 0 ) const;
    
    
    // error on a content of a bin
    Float_t GetErrors( Int_t, Int_t nbinsY= 0 ) const;
    
    
    // contents of a bin
    Int_t GetEntries() const;
    // total number of entries
    Float_t GetWtSum() const;
    // weighted sum
    Int_t GetNbins( Int_t ) const;
    // number of bins in a dimension
    Float_t GetLow( Int_t ) const;
    // low edge of a dimension
    Float_t GetHigh( Int_t ) const;
    // high edge of a dimension
    Float_t GetAvg( Int_t ) const; // average X * wtsum
    // average of theDim
    Float_t GetCovar( Int_t, Int_t dim2= 0 ) const;
    // covariance between the two dimensions
    Int_t GetType() const;
    
    // Return the HHistID:
    HHistID GetHistID() const;
    
    // Test if input pointer matches pointer to ROOT histo:
    Bool_t PtrIsEqual( TObject* ptr ) const;
    
    HHistogram& operator<<(Axis_t x);

    // Extra functionality
    void Fill( Axis_t x ) { Accumulate1(x); }
    void Fill( Axis_t x, Axis_t y ) { Accumulate2(x,y); }
    void SetFillColor(Color_t fcolor) { histp->SetFillColor(fcolor); }
    virtual void Draw(Option_t *option="") { histp->Draw(option); }

private:
    
    // Satisfy Scotts weird function
    void setEntries( Int_t ) {}; 
    
    
    // Data membrs:
    TH1* histp;
    
public:
    ClassDef(HHistogram,1) //T histogram
};

//----------------Spezialization-----------------

class HMassHistogram : public HHistogram {
    
public:
    //Constructor
    HMassHistogram(const Text_t *name,const Text_t *title,Int_t nbins,Axis_t xlow,Axis_t xup); 
    void Accumulate( Axis_t, Stat_t weight= 1.0 );
    //Destructor
    virtual ~HMassHistogram();          
private:
public:
    ClassDef(HMassHistogram,1) //Mass histogram
}; 

class HEnergyHistogram : public HHistogram {
    
public:
    //Constructor
    HEnergyHistogram(const Text_t *name,const Text_t *title,Int_t nbins,Axis_t xlow,Axis_t xup);
    void Accumulate( Axis_t, Stat_t weight= 1.0 );
    //Destructor
    virtual ~HEnergyHistogram();          
private:
public:
    ClassDef(HEnergyHistogram,1) //Energy histogram
}; 

class HMomentumHistogram : public HHistogram {
    
public:
    //Constructor
    HMomentumHistogram(const Text_t *name,const Text_t *title,Int_t nbins,Axis_t xlow,Axis_t xup);
    void Accumulate( Axis_t, Stat_t weight= 1.0 );
    //Destructor
    virtual ~HMomentumHistogram();        
private:
public:
    ClassDef(HMomentumHistogram,1) //Momentum histogram
};

class HEoverPHistogram : public HHistogram {
    
public:
    //Constructor
    HEoverPHistogram(const Text_t *name,const Text_t *title,Int_t nbinsx,Axis_t xlow,Axis_t xup
        ,Int_t nbinsy,Axis_t ylow,Axis_t yup);
    void Accumulate( Axis_t x, Axis_t y, Stat_t weight=1.0 );
    //Destructor
    virtual ~HEoverPHistogram();        
private:
public:
    ClassDef(HEoverPHistogram,1) //Energy vs. momentum histogram
};

class HMoverPHistogram : public HHistogram {
    
public:
    //Constructor
    HMoverPHistogram(const Text_t *name,const Text_t *title,Int_t nbinsx,Axis_t xlow,Axis_t xup
        ,Int_t nbinsy,Axis_t ylow,Axis_t yup);
    void Accumulate( Axis_t x, Axis_t y, Stat_t weight=1.0 );
    //Destructor
    virtual ~HMoverPHistogram();        
private:
public:
    ClassDef(HMoverPHistogram,1) //Mass vs. momentum histogram
};

class HDalitzPlot : public HHistogram {
    
public:
    //Constructor
    HDalitzPlot(const Text_t *name,const Text_t *title,Int_t nbinsx,Axis_t xlow,Axis_t xup
        ,Int_t nbinsy,Axis_t ylow,Axis_t yup);
    //Destructor
    virtual ~HDalitzPlot();        
private:
public:
    ClassDef(HDalitzPlot,1) //Generate a Dalitz plot
};

#endif

Revision:	1.1
Committed:	Wed Sep 29 15:44:21 2010 UTC (14 years, 7 months ago) by yjlee
Content type:	text/plain
Branch:	MAIN
CVS Tags:	HEAD
Error occurred while calculating annotation data.
Log Message:	HafHistogram tool
#	Content
1	#ifndef HHISTOGRAM_H
2	#define HHISTOGRAM_H
3
4	/** \class HHistogram Ganzhur/HafHistogram/interface/HHistogram.h
5	*
6	* Description:
7	* Analysis code of the CMS experiment;
8	* Original version is a part of HAF package developed for CMS: UserCode/HAF
9	* Histogram class
10	*
11	* \author Marcel Kunze, Ruhr-University Bochum, Germany
12	* \author Serguei Ganjour, CEA-Saclay/IRFU, FR
13	*
14	* \version $Id: HHistogram.h,v 1.2 2009/11/19 09:40:31 ganzhur Exp $
15	*
16	*/
17
18	#include "TH1.h"
19	#include "TH2.h"
20
21	class HCandList;
22	class HCandidate;
23
24	class TObject;
25	class TH1;
26	class HHistID;
27
28	class HHistogram : public TNamed {
29
30	public:
31
32	// Constructors:
33	// 1-D:
34	HHistogram( const char, const char, Int_t, Axis_t, Axis_t );
35	// 2-D:
36	HHistogram( const char, const char, Int_t, Axis_t, Axis_t, Int_t, Axis_t, Axis_t );
37	// Profile:
38	HHistogram( const char, const char, Int_t, Axis_t, Axis_t, Axis_t, Axis_t );
39
40	// Destructor:
41	virtual ~HHistogram();
42
43	// Functions to fill a histo:
44	void Accumulate( Axis_t, Stat_t weight= 1.0 );
45	void Accumulate1( Axis_t, Stat_t weight= 1.0 );
46	void Accumulate( Axis_t x, Axis_t y, Stat_t weight );
47	void Accumulate2( Axis_t, Axis_t, Stat_t weight= 1.0 );
48
49	// clear all contents
50	void Reset();
51
52	// accessors:
53	const char* Title() const { return GetTitle(); };
54	Float_t GetContents( Int_t, Int_t nbinsY= 0 ) const;
55
56
57	// error on a content of a bin
58	Float_t GetErrors( Int_t, Int_t nbinsY= 0 ) const;
59
60
61	// contents of a bin
62	Int_t GetEntries() const;
63	// total number of entries
64	Float_t GetWtSum() const;
65	// weighted sum
66	Int_t GetNbins( Int_t ) const;
67	// number of bins in a dimension
68	Float_t GetLow( Int_t ) const;
69	// low edge of a dimension
70	Float_t GetHigh( Int_t ) const;
71	// high edge of a dimension
72	Float_t GetAvg( Int_t ) const; // average X * wtsum
73	// average of theDim
74	Float_t GetCovar( Int_t, Int_t dim2= 0 ) const;
75	// covariance between the two dimensions
76	Int_t GetType() const;
77
78	// Return the HHistID:
79	HHistID GetHistID() const;
80
81	// Test if input pointer matches pointer to ROOT histo:
82	Bool_t PtrIsEqual( TObject* ptr ) const;
83
84	HHistogram& operator<<(Axis_t x);
85
86	// Extra functionality
87	void Fill( Axis_t x ) { Accumulate1(x); }
88	void Fill( Axis_t x, Axis_t y ) { Accumulate2(x,y); }
89	void SetFillColor(Color_t fcolor) { histp->SetFillColor(fcolor); }
90	virtual void Draw(Option_t *option="") { histp->Draw(option); }
91
92	private:
93
94	// Satisfy Scotts weird function
95	void setEntries( Int_t ) {};
96
97
98	// Data membrs:
99	TH1* histp;
100
101	public:
102	ClassDef(HHistogram,1) //T histogram
103	};
104
105	//----------------Spezialization-----------------
106
107	class HMassHistogram : public HHistogram {
108
109	public:
110	//Constructor
111	HMassHistogram(const Text_t name,const Text_t title,Int_t nbins,Axis_t xlow,Axis_t xup);
112	void Accumulate( Axis_t, Stat_t weight= 1.0 );
113	//Destructor
114	virtual ~HMassHistogram();
115	private:
116	public:
117	ClassDef(HMassHistogram,1) //Mass histogram
118	};
119
120	class HEnergyHistogram : public HHistogram {
121
122	public:
123	//Constructor
124	HEnergyHistogram(const Text_t name,const Text_t title,Int_t nbins,Axis_t xlow,Axis_t xup);
125	void Accumulate( Axis_t, Stat_t weight= 1.0 );
126	//Destructor
127	virtual ~HEnergyHistogram();
128	private:
129	public:
130	ClassDef(HEnergyHistogram,1) //Energy histogram
131	};
132
133	class HMomentumHistogram : public HHistogram {
134
135	public:
136	//Constructor
137	HMomentumHistogram(const Text_t name,const Text_t title,Int_t nbins,Axis_t xlow,Axis_t xup);
138	void Accumulate( Axis_t, Stat_t weight= 1.0 );
139	//Destructor
140	virtual ~HMomentumHistogram();
141	private:
142	public:
143	ClassDef(HMomentumHistogram,1) //Momentum histogram
144	};
145
146	class HEoverPHistogram : public HHistogram {
147
148	public:
149	//Constructor
150	HEoverPHistogram(const Text_t name,const Text_t title,Int_t nbinsx,Axis_t xlow,Axis_t xup
151	,Int_t nbinsy,Axis_t ylow,Axis_t yup);
152	void Accumulate( Axis_t x, Axis_t y, Stat_t weight=1.0 );
153	//Destructor
154	virtual ~HEoverPHistogram();
155	private:
156	public:
157	ClassDef(HEoverPHistogram,1) //Energy vs. momentum histogram
158	};
159
160	class HMoverPHistogram : public HHistogram {
161
162	public:
163	//Constructor
164	HMoverPHistogram(const Text_t name,const Text_t title,Int_t nbinsx,Axis_t xlow,Axis_t xup
165	,Int_t nbinsy,Axis_t ylow,Axis_t yup);
166	void Accumulate( Axis_t x, Axis_t y, Stat_t weight=1.0 );
167	//Destructor
168	virtual ~HMoverPHistogram();
169	private:
170	public:
171	ClassDef(HMoverPHistogram,1) //Mass vs. momentum histogram
172	};
173
174	class HDalitzPlot : public HHistogram {
175
176	public:
177	//Constructor
178	HDalitzPlot(const Text_t name,const Text_t title,Int_t nbinsx,Axis_t xlow,Axis_t xup
179	,Int_t nbinsy,Axis_t ylow,Axis_t yup);
180	//Destructor
181	virtual ~HDalitzPlot();
182	private:
183	public:
184	ClassDef(HDalitzPlot,1) //Generate a Dalitz plot
185	};
186
187	#endif
188