diff --git a/docs/XeLikelihoods_8h_source.html b/docs/XeLikelihoods_8h_source.html index 69d7837..d9395b4 100644 --- a/docs/XeLikelihoods_8h_source.html +++ b/docs/XeLikelihoods_8h_source.html @@ -241,9 +241,9 @@
LKParameter
Definition: XeStat.h:2991
dataHandler
Definition: dataHandler.h:35
pdfLikelihood
Definition: XeLikelihoods.h:20
-
histoCompare
Definition: XePdfObjects.h:161
+
histoCompare
Definition: XePdfObjects.h:164
ProfileLikelihood
Definition: XeStat.h:3463
-
pdfComponent
Definition: XePdfObjects.h:68
+
pdfComponent
Definition: XePdfObjects.h:71
diff --git a/docs/XePdfObjects_8h_source.html b/docs/XePdfObjects_8h_source.html index d061e9a..e7a86c9 100644 --- a/docs/XePdfObjects_8h_source.html +++ b/docs/XePdfObjects_8h_source.html @@ -109,208 +109,208 @@
26 
27  public:
28 
-
30  shapeSys(TString name);
-
31 
-
32 
-
34  double getNearestLow();
+
33  shapeSys(TString name);
+
34 
35 
-
37  double getNearestHigh();
-
38 
-
39  ClassDef(shapeSys,1)
-
40 };
-
41 
-
42 
-
43 class scaleSys : public LKParameter {
+
37  double getNearestLow();
+
38 
+
40  double getNearestHigh();
+
41 
+
42  ClassDef(shapeSys,1)
+
43 };
44 
-
45  public:
-
46  // relativeUncertainty = 1 sigma relative uncertainty (example 10% uncertainty would be 0.1)
-
47  scaleSys(TString name, double relativeUncertainty);
-
49 
-
50  // gives the Normalization scale factor according to t-valued uncertainty --> (1 + t_val * relativeUncertainty )
-
51  // Uses the current t-value which is automatically linked to its pdfLikelihood
-
52  double getNormModifier();
-
53 
-
54 
-
55  //this case return a sys that is centered in zero, half a gaussian, t-value strictly positive
-
56  //and for a tvalue=0 have zero events. Histo is supposed to be normalized to 1
-
57  void setNull(){ isNull = true; setMinimum(0.);};
-
58 
-
59  private:
-
60  double relUnc;
-
61  bool isNull;
-
62 
-
63 
-
64  ClassDef(scaleSys,1)
-
65 };
+
45 
+
46 class scaleSys : public LKParameter {
+
47 
+
48  public:
+
49  // relativeUncertainty = 1 sigma relative uncertainty (example 10% uncertainty would be 0.1)
+
50  scaleSys(TString name, double relativeUncertainty);
+
52 
+
53  // gives the Normalization scale factor according to t-valued uncertainty --> (1 + t_val * relativeUncertainty )
+
54  // Uses the current t-value which is automatically linked to its pdfLikelihood
+
55  double getNormModifier();
+
56 
+
57 
+
58  //this case return a sys that is centered in zero, half a gaussian, t-value strictly positive
+
59  //and for a tvalue=0 have zero events. Histo is supposed to be normalized to 1
+
60  void setNull(){ isNull = true; setMinimum(0.);};
+
61 
+
62  private:
+
63  double relUnc;
+
64  bool isNull;
+
65 
66 
-
67 
-
68 class pdfComponent :errorHandler{
+
67  ClassDef(scaleSys,1)
+
68 };
69 
-
70  public:
-
71 
-
72  pdfComponent(TString name, TString filename);
-
73  ~pdfComponent();
+
70 
+
71 class pdfComponent :errorHandler{
+
72 
+
73  public:
74 
-
75 
-
76  void addScaleSys(scaleSys *addMe) { myScaleUnc.push_back(addMe); };
+
75  pdfComponent(TString name, TString filename);
+
76  ~pdfComponent();
77 
-
78  void addShapeSys(shapeSys *addMe) { myShapeUnc.push_back(addMe); };
-
79 
-
80  //load histogram according to the current value of the parameters
-
81  void loadHistos();
+
78 
+
79  void addScaleSys(scaleSys *addMe) { myScaleUnc.push_back(addMe); };
+
80 
+
81  void addShapeSys(shapeSys *addMe) { myShapeUnc.push_back(addMe); };
82 
-
83  //load default histogram, no sys.
-
84  void loadDefaultHisto();
+
83  //load histogram according to the current value of the parameters
+
84  void loadHistos();
85 
-
86  // returns the interpolated pdf over shape sys computed in (s1,s2),
-
87  // it is normalized to number of events and scale uncertainty are also
-
88  // taken into account
-
89  double getNormalizedDensity(double s1, double s2);
-
90 
-
91  // returns the (s1,s2) bin content of the default histo, no shape nor scale sys is applied
-
92  double getDefaultDensity(double s1, double s2);
+
86  //load default histogram, no sys.
+
87  void loadDefaultHisto();
+
88 
+
89  // returns the interpolated pdf over shape sys computed in (s1,s2),
+
90  // it is normalized to number of events and scale uncertainty are also
+
91  // taken into account
+
92  double getNormalizedDensity(double s1, double s2);
93 
-
94  // Returns the total integrated number of events taking into account shape sys
-
95  // and scale sys.
-
96  double getNormalizedEvents();
-
97 
-
98  // returns total integral of the default histo, no shape nor scale sys is applied
-
99  double getDefaultEvents();
-
100 
-
101  // Returns a copy of the interpolated histogram according to the "currentValue" of the
-
102  // shapeSys associated with it. The current value of each uncertainty can
-
103  // be set simply by shapeSys::setCurrentValue(val) method.
-
104  TH2F getInterpolatedHisto();
-
105 
-
106  // Returns a copy of the default histogram, no shape nor scale sys is applied
-
107  TH2F getDefaultHisto();
+
94  // returns the (s1,s2) bin content of the default histo, no shape nor scale sys is applied
+
95  double getDefaultDensity(double s1, double s2);
+
96 
+
97  // Returns the total integrated number of events taking into account shape sys
+
98  // and scale sys.
+
99  double getNormalizedEvents();
+
100 
+
101  // returns total integral of the default histo, no shape nor scale sys is applied
+
102  double getDefaultEvents();
+
103 
+
104  // Returns a copy of the interpolated histogram according to the "currentValue" of the
+
105  // shapeSys associated with it. The current value of each uncertainty can
+
106  // be set simply by shapeSys::setCurrentValue(val) method.
+
107  TH2F getInterpolatedHisto();
108 
-
109  // returns the grid points in histogram space that needs to be loaded
-
110  // this method need to be modified for arbitrary number of shape sys
-
111  TString getNearestHistoName(vector<bool> setOfVal);
-
112 
-
113  // returns the default name of the histogram
-
114  TString getDefaultHistoName();
+
109  // Returns a copy of the default histogram, no shape nor scale sys is applied
+
110  TH2F getDefaultHisto();
+
111 
+
112  // returns the grid points in histogram space that needs to be loaded
+
113  // this method need to be modified for arbitrary number of shape sys
+
114  TString getNearestHistoName(vector<bool> setOfVal);
115 
-
116  // returns a string with current values of parameters
-
117  TString getParamValueString();
+
116  // returns the default name of the histogram
+
117  TString getDefaultHistoName();
118 
-
119  TString getParamValueWritable();
-
120 
-
121  // returns the integral of the square between points (not bins), linearly interpolates
-
122  // the area if desn't correspond to an integer number of bins.
-
123  // Uses the default histo, does not consider shape or scale uncertainty.
-
124  double getDefaultPdfIntegral(double s1_min, double s1_max, double s2_min, double s2_max);
-
125 
-
126  //This scans the parameter space given a number of steps for each parameter
-
127  //produce a projection given min and max, produce a PDF file for comparison
-
128  //of all interpolation and saves also the TH2F to file.
-
129  void plotInterpolatedSpace(bool doProjectionX, double min, double max, int Nsteps);
-
130  //scale the pdf by VAL, this happend for the methods:
-
131  //getDefaultHisto, getInterpolatedHisto, getNormalizedDensity,
-
132  //getDefaultDensity, getNormalizedEvents, getDefaultEvents
-
133  void setScaleFactor(double val) {scaleFactor = val;};
-
134 
-
135  TString getName() {return pdf_name;};
-
136 
-
137  vector<scaleSys*> myScaleUnc; //container of scale sys
-
138  vector<shapeSys*> myShapeUnc; //container of shape sys
+
119  // returns a string with current values of parameters
+
120  TString getParamValueString();
+
121 
+
122  TString getParamValueWritable();
+
123 
+
124  // returns the integral of the square between points (not bins), linearly interpolates
+
125  // the area if desn't correspond to an integer number of bins.
+
126  // Uses the default histo, does not consider shape or scale uncertainty.
+
127  double getDefaultPdfIntegral(double s1_min, double s1_max, double s2_min, double s2_max);
+
128 
+
129  //This scans the parameter space given a number of steps for each parameter
+
130  //produce a projection given min and max, produce a PDF file for comparison
+
131  //of all interpolation and saves also the TH2F to file.
+
132  void plotInterpolatedSpace(bool doProjectionX, double min, double max, int Nsteps);
+
133  //scale the pdf by VAL, this happend for the methods:
+
134  //getDefaultHisto, getInterpolatedHisto, getNormalizedDensity,
+
135  //getDefaultDensity, getNormalizedEvents, getDefaultEvents
+
136  void setScaleFactor(double val) {scaleFactor = val;};
+
137 
+
138  TString getName() {return pdf_name;};
139 
-
140  TString suffix;
-
141  bool doExtend;
+
140  vector<scaleSys*> myScaleUnc; //container of scale sys
+
141  vector<shapeSys*> myShapeUnc; //container of shape sys
142 
-
143  private:
-
144  TFile *file;
-
145  TH2F *defaultDistro;
-
146  vector<TH2F*> histos; // contains the 2^N histo for the hyperplane interpolation of shapeSys
-
147  vector<double> InterpFactors; // contains the 2^N interpolation factors for the hyperplane interpolation of shapeSys
-
148  TString pdf_name;
-
149  vector<double> old_t_val; //contains the last value interpolated, the interpolation is lazy, doesn't ricompute it if is for the same set of values.
-
150  double scaleFactor;
-
151 
-
152 
-
153  void extendHisto(TH2F &h);
+
143  TString suffix;
+
144  bool doExtend;
+
145 
+
146  private:
+
147  TFile *file;
+
148  TH2F *defaultDistro;
+
149  vector<TH2F*> histos; // contains the 2^N histo for the hyperplane interpolation of shapeSys
+
150  vector<double> InterpFactors; // contains the 2^N interpolation factors for the hyperplane interpolation of shapeSys
+
151  TString pdf_name;
+
152  vector<double> old_t_val; //contains the last value interpolated, the interpolation is lazy, doesn't ricompute it if is for the same set of values.
+
153  double scaleFactor;
154 
-
155  ClassDef(pdfComponent,1)
-
156 };
+
155 
+
156  void extendHisto(TH2F &h);
157 
-
158 
-
159 
+
158  ClassDef(pdfComponent,1)
+
159 };
160 
-
161 class histoCompare :errorHandler{
+
161 
162 
-
163  // This class is supposed to compare a "base" histo
-
164  // with a set of other histogram. One can choose for ratio plot,
-
165  // stack plot of different component, one can set the projected axis
-
166  // and the slice on which do the projection, one can also rebin the histo.
-
167  // the histo are supposed to be TH2F.
-
168 
-
169  public:
-
170 
-
171  histoCompare();
-
172  ~histoCompare();
+
163 
+
164 class histoCompare :errorHandler{
+
165 
+
166  // This class is supposed to compare a "base" histo
+
167  // with a set of other histogram. One can choose for ratio plot,
+
168  // stack plot of different component, one can set the projected axis
+
169  // and the slice on which do the projection, one can also rebin the histo.
+
170  // the histo are supposed to be TH2F.
+
171 
+
172  public:
173 
-
174  //---------OPTIONS-----------//
-
175  int rebinX;
-
176  int rebinY;
-
177  bool projectionX;
-
178  double projectionMin; // minimum value on the projected axis
-
179  double projectionMax; // maximum value on the projected axis
-
180  int binMin; // same but in bin nuber
-
181  int binMax;
-
182  bool doStack; // True if you want to sum up histo in compareList
-
183  TString titleX;
-
184  TString titleY;
-
185  //---------------------------//
-
186 
-
187  // set histo with wich one one to compare the rest
-
188  void setBaseHisto(TH2F b, TString n=""){base = b; names[0] = n;};
-
189  // set the comparison histos
-
190  void addHistoToList(TH2F h, TString n ="") {compareList.push_back(h); names.push_back(n);};
-
191  //set Name of component
-
192  void setNameofComponent(unsigned int i, TString n);
-
193  // draw just a normal comparison plot
-
194  void compare();
-
195  // add a ratio plot to compare()
-
196  void compareWithRatio();
-
197 
-
198  void drawLegend(TH1D *baseH, vector <TH1D*> list );
-
199 
-
200  //return a string of info
-
201  TString projectionInfo();
-
202 
-
203  //------ histo holders --------//
-
204  TH2F base;
-
205  vector <TString> names;
-
206  vector <TH2F> compareList;
-
207 
-
208  TH1D *projectedBase;
-
209  vector <TH1D*> projectedList;
+
174  histoCompare();
+
175  ~histoCompare();
+
176 
+
177  //---------OPTIONS-----------//
+
178  int rebinX;
+
179  int rebinY;
+
180  bool projectionX;
+
181  double projectionMin; // minimum value on the projected axis
+
182  double projectionMax; // maximum value on the projected axis
+
183  int binMin; // same but in bin nuber
+
184  int binMax;
+
185  bool doStack; // True if you want to sum up histo in compareList
+
186  TString titleX;
+
187  TString titleY;
+
188  //---------------------------//
+
189 
+
190  // set histo with wich one one to compare the rest
+
191  void setBaseHisto(TH2F b, TString n=""){base = b; names[0] = n;};
+
192  // set the comparison histos
+
193  void addHistoToList(TH2F h, TString n ="") {compareList.push_back(h); names.push_back(n);};
+
194  //set Name of component
+
195  void setNameofComponent(unsigned int i, TString n);
+
196  // draw just a normal comparison plot
+
197  void compare();
+
198  // add a ratio plot to compare()
+
199  void compareWithRatio();
+
200 
+
201  void drawLegend(TH1D *baseH, vector <TH1D*> list );
+
202 
+
203  //return a string of info
+
204  TString projectionInfo();
+
205 
+
206  //------ histo holders --------//
+
207  TH2F base;
+
208  vector <TString> names;
+
209  vector <TH2F> compareList;
210 
-
211  private:
-
212 
-
213  void setOptions(TH1D *h, bool dataLike = true, bool isBottom = false);
-
214  void project(); // push histo into the TH1D
+
211  TH1D *projectedBase;
+
212  vector <TH1D*> projectedList;
+
213 
+
214  private:
215 
-
216 
-
217  ClassDef(histoCompare,1)
-
218 
+
216  void setOptions(TH1D *h, bool dataLike = true, bool isBottom = false);
+
217  void project(); // push histo into the TH1D
+
218 
219 
-
220 
-
221 };
+
220  ClassDef(histoCompare,1)
+
221 
222 
223 
-
224 
+
224 };
225 
226 
227 
-
228 #endif
+
228 
+
229 
+
230 
+
231 #endif
LKParameter
Definition: XeStat.h:2991
shapeSys
it represents a Shape Systematic. It inherits from LKParameter. Here you can use the LKParameter func...
Definition: XePdfObjects.h:25
errorHandler
Definition: XeCore.h:2142
-
scaleSys
Definition: XePdfObjects.h:43
-
histoCompare
Definition: XePdfObjects.h:161
-
pdfComponent
Definition: XePdfObjects.h:68
+
scaleSys
Definition: XePdfObjects.h:46
+
histoCompare
Definition: XePdfObjects.h:164
+
pdfComponent
Definition: XePdfObjects.h:71
diff --git a/docs/classshapeSys.html b/docs/classshapeSys.html index 4694a32..961ed5e 100644 --- a/docs/classshapeSys.html +++ b/docs/classshapeSys.html @@ -761,9 +761,10 @@
+

Constructor for class shapeSys.

Parameters
- +
name,nameof the sys as appears in histogram name, see pdfComponent for a detailed use case.
nameName of the sys as appears in histogram name, see pdfComponent for a detailed use case.
diff --git a/src/XePdfObjects.h b/src/XePdfObjects.h index d5f4d8f..70ccb06 100644 --- a/src/XePdfObjects.h +++ b/src/XePdfObjects.h @@ -26,7 +26,10 @@ class shapeSys : public LKParameter { public: - //! @param name, name of the sys as appears in histogram name, see pdfComponent for a detailed use case. + /** + * Constructor for class shapeSys. + * @param name: Name of the sys as appears in histogram name, see pdfComponent for a detailed use case. + */ shapeSys(TString name);