ViewVC Help
View File | Revision Log | Show Annotations | Root Listing
root/cvsroot/UserCode/MitCommon/DataFormats/src/TH2DAsymErr.cc
Revision: 1.4
Committed: Sat Mar 27 13:20:00 2010 UTC (15 years, 1 month ago) by sixie
Content type: text/plain
Branch: MAIN
CVS Tags: Mit_032, Mit_031, Mit_025c_branch2, Mit_025c_branch1, Mit_030, Mit_029c, Mit_030_pre1, Mit_029a, Mit_029, Mit_029_pre1, Mit_028a, Mit_025c_branch0, Mit_028, Mit_027a, Mit_027, Mit_026, Mit_025e, Mit_025d, Mit_025c, Mit_025b, Mit_025a, Mit_025, Mit_025pre2, Mit_024b, Mit_025pre1, Mit_024a, Mit_024, Mit_023, Mit_022a, Mit_022, Mit_020d, TMit_020d, Mit_020c, Mit_021, Mit_021pre2, Mit_021pre1, Mit_020b, Mit_020a, Mit_020, Mit_020pre1, Mit_018, Mit_017, Mit_017pre3, Mit_017pre2, Mit_017pre1, V07-05-00, Mit_016, Mit_015b, Mit_015a, Mit_015, Mit_014e, Mit_014d, Mit_014c, Mit_014b, ConvRejection-10-06-09, Mit_014a, Mit_014, Mit_014pre3, Mit_014pre2, Mit_014pre1, Mit_013d, Mit_013c, Mit_013b, Mit_013a, HEAD
Branch point for: Mit_025c_branch
Changes since 1.3: +27 -1 lines
Log Message: 
add some accessors to get bin content by bin number

File Contents

# User Rev Content
1 sixie 1.4 // $Id: TH2DAsymErr.cc,v 1.3 2009/10/21 15:05:22 sixie Exp $
2 phedex 1.1
3     #include "MitCommon/DataFormats/interface/TH2DAsymErr.h"
4     #include <TList.h>
5 sixie 1.3 #include <iostream>
6 phedex 1.1
7     ClassImp(mithep::TH2DAsymErr)
8    
9     using namespace mithep;
10 sixie 1.3 using namespace std;
11 phedex 1.1
12     //--------------------------------------------------------------------------------------------------
13 loizides 1.2 TH2DAsymErr::TH2DAsymErr(const char *name, const char *title,
14     Int_t nbinsx, Double_t xlow, Double_t xup
15     ,Int_t nbinsy, Double_t ylow, Double_t yup) :
16     TH2D(name, title, nbinsx, xlow, xup , nbinsy, ylow, yup),
17     fStatErrorLow(fNcells),
18     fStatErrorHigh(fNcells),
19     fSysErrorLow(fNcells),
20     fSysErrorHigh(fNcells)
21 phedex 1.1 {
22     // Constructor.
23     }
24    
25 loizides 1.2 TH2DAsymErr::TH2DAsymErr(const TH2D &h2d) :
26     TH2D(h2d),
27     fStatErrorLow(fNcells),
28     fStatErrorHigh(fNcells),
29     fSysErrorLow(fNcells),
30     fSysErrorHigh(fNcells)
31 phedex 1.1 {
32     // Constructor.
33     }
34    
35     //--------------------------------------------------------------------------------------------------
36 loizides 1.2 Double_t TH2DAsymErr::GetError(Double_t x, Double_t y, EErrType t)
37 phedex 1.1 {
38 loizides 1.2 // Get error corresponding to x and y for given error type.
39 phedex 1.1
40 sixie 1.3 Int_t nx = fXaxis.GetNbins()+2;
41     Int_t binx = fXaxis.FindFixBin(x);
42     Int_t biny = fYaxis.FindFixBin(y);
43     Int_t bin = binx + nx*biny;
44    
45 loizides 1.2 switch (t) {
46     case kStatErrLow:
47     return fStatErrorLow.fArray[bin];
48     break;
49     case kStatErrHigh:
50     return fStatErrorHigh.fArray[bin];
51     break;
52     case kSysErrLow:
53     return fSysErrorLow.fArray[bin];
54     break;
55     case kSysErrHigh:
56     return fSysErrorHigh.fArray[bin];
57     break;
58     }
59     return 0;
60 phedex 1.1 }
61    
62     //--------------------------------------------------------------------------------------------------
63 sixie 1.4 Double_t TH2DAsymErr::GetBinError(Int_t b, Int_t c, EErrType t)
64     {
65     // Get error corresponding to bin (b,c) for given error type.
66    
67     Int_t nx = fXaxis.GetNbins()+2;
68     Int_t bin = b + nx*c;
69    
70     switch (t) {
71     case kStatErrLow:
72     return fStatErrorLow.fArray[bin];
73     break;
74     case kStatErrHigh:
75     return fStatErrorHigh.fArray[bin];
76     break;
77     case kSysErrLow:
78     return fSysErrorLow.fArray[bin];
79     break;
80     case kSysErrHigh:
81     return fSysErrorHigh.fArray[bin];
82     break;
83     }
84     return 0;
85     }
86    
87    
88     //--------------------------------------------------------------------------------------------------
89 loizides 1.2 void TH2DAsymErr::SetBinContent(Int_t binx, Int_t biny, Double_t value,
90     Double_t statErrorLow, Double_t statErrorHigh,
91     Double_t sysErrorLow, Double_t sysErrorHigh)
92 phedex 1.1 {
93 loizides 1.2 // Set bin content for given bin and value.
94 phedex 1.1
95 loizides 1.2 Int_t bin = GetBin(binx,biny);
96     SetBinContent(bin,value);
97 phedex 1.1
98 loizides 1.2 fStatErrorLow.fArray[bin] = statErrorLow;
99     fStatErrorHigh.fArray[bin] = statErrorHigh;
100     fSysErrorLow.fArray[bin] = sysErrorLow;
101     fSysErrorHigh.fArray[bin] = sysErrorHigh;
102 sixie 1.3 return;
103 phedex 1.1 }
104    
105     //--------------------------------------------------------------------------------------------------
106 loizides 1.2 void TH2DAsymErr::SetBinError(Int_t binx, Int_t biny,
107     Double_t statErrorLow, Double_t statErrorHigh,
108     Double_t sysErrorLow, Double_t sysErrorHigh)
109 phedex 1.1 {
110 loizides 1.2 // Set bin error for asymetric errors.
111 phedex 1.1
112 loizides 1.2 Int_t bin = GetBin(binx,biny);
113 phedex 1.1
114 loizides 1.2 fStatErrorLow.fArray[bin] = statErrorLow;
115     fStatErrorHigh.fArray[bin] = statErrorHigh;
116     fSysErrorLow.fArray[bin] = sysErrorLow;
117     fSysErrorHigh.fArray[bin] = sysErrorHigh;
118 sixie 1.3 return;
119 phedex 1.1 }