diff --git a/docs/XeLikelihoods_8h_source.html b/docs/XeLikelihoods_8h_source.html
index 82e9abd..69d7837 100644
--- a/docs/XeLikelihoods_8h_source.html
+++ b/docs/XeLikelihoods_8h_source.html
@@ -241,9 +241,9 @@
<div class="ttc" id="classLKParameter_html"><div class="ttname"><a href="classLKParameter.html">LKParameter</a></div><div class="ttdef"><b>Definition:</b> XeStat.h:2991</div></div>
<div class="ttc" id="classdataHandler_html"><div class="ttname"><a href="classdataHandler.html">dataHandler</a></div><div class="ttdef"><b>Definition:</b> dataHandler.h:35</div></div>
<div class="ttc" id="classpdfLikelihood_html"><div class="ttname"><a href="classpdfLikelihood.html">pdfLikelihood</a></div><div class="ttdef"><b>Definition:</b> XeLikelihoods.h:20</div></div>
-<div class="ttc" id="classhistoCompare_html"><div class="ttname"><a href="classhistoCompare.html">histoCompare</a></div><div class="ttdef"><b>Definition:</b> XePdfObjects.h:160</div></div>
+<div class="ttc" id="classhistoCompare_html"><div class="ttname"><a href="classhistoCompare.html">histoCompare</a></div><div class="ttdef"><b>Definition:</b> XePdfObjects.h:161</div></div>
<div class="ttc" id="classProfileLikelihood_html"><div class="ttname"><a href="classProfileLikelihood.html">ProfileLikelihood</a></div><div class="ttdef"><b>Definition:</b> XeStat.h:3463</div></div>
-<div class="ttc" id="classpdfComponent_html"><div class="ttname"><a href="classpdfComponent.html">pdfComponent</a></div><div class="ttdef"><b>Definition:</b> XePdfObjects.h:67</div></div>
+<div class="ttc" id="classpdfComponent_html"><div class="ttname"><a href="classpdfComponent.html">pdfComponent</a></div><div class="ttdef"><b>Definition:</b> XePdfObjects.h:68</div></div>
</div><!-- fragment --></div><!-- contents -->
<!-- start footer part -->
<hr class="footer"/><address class="footer"><small>
diff --git a/docs/XePdfObjects_8h_source.html b/docs/XePdfObjects_8h_source.html
index c3bf5f7..d061e9a 100644
--- a/docs/XePdfObjects_8h_source.html
+++ b/docs/XePdfObjects_8h_source.html
@@ -105,212 +105,212 @@
<div class="line"><a name="l00016"></a><span class="lineno"> 16</span> <span class="keyword">using namespace </span>std;</div>
<div class="line"><a name="l00017"></a><span class="lineno"> 17</span> </div>
<div class="line"><a name="l00018"></a><span class="lineno"> 18</span> </div>
-<div class="line"><a name="l00024"></a><span class="lineno"><a class="line" href="classshapeSys.html"> 24</a></span> <span class="keyword">class </span><a class="code" href="classshapeSys.html">shapeSys</a> : <span class="keyword">public</span> <a class="code" href="classLKParameter.html">LKParameter</a> {</div>
-<div class="line"><a name="l00025"></a><span class="lineno"> 25</span> </div>
-<div class="line"><a name="l00026"></a><span class="lineno"> 26</span>  <span class="keyword">public</span>:</div>
-<div class="line"><a name="l00027"></a><span class="lineno"> 27</span> </div>
-<div class="line"><a name="l00029"></a><span class="lineno"> 29</span>  <a class="code" href="classshapeSys.html">shapeSys</a>(TString name);</div>
-<div class="line"><a name="l00030"></a><span class="lineno"> 30</span> </div>
+<div class="line"><a name="l00025"></a><span class="lineno"><a class="line" href="classshapeSys.html"> 25</a></span> <span class="keyword">class </span><a class="code" href="classshapeSys.html">shapeSys</a> : <span class="keyword">public</span> <a class="code" href="classLKParameter.html">LKParameter</a> {</div>
+<div class="line"><a name="l00026"></a><span class="lineno"> 26</span> </div>
+<div class="line"><a name="l00027"></a><span class="lineno"> 27</span>  <span class="keyword">public</span>:</div>
+<div class="line"><a name="l00028"></a><span class="lineno"> 28</span> </div>
+<div class="line"><a name="l00030"></a><span class="lineno"> 30</span>  <a class="code" href="classshapeSys.html">shapeSys</a>(TString name);</div>
<div class="line"><a name="l00031"></a><span class="lineno"> 31</span> </div>
-<div class="line"><a name="l00033"></a><span class="lineno"> 33</span>  <span class="keywordtype">double</span> getNearestLow(); </div>
-<div class="line"><a name="l00034"></a><span class="lineno"> 34</span> </div>
-<div class="line"><a name="l00036"></a><span class="lineno"> 36</span>  <span class="keywordtype">double</span> getNearestHigh();</div>
-<div class="line"><a name="l00037"></a><span class="lineno"> 37</span>  </div>
-<div class="line"><a name="l00038"></a><span class="lineno"> 38</span>  ClassDef(<a class="code" href="classshapeSys.html">shapeSys</a>,1)</div>
-<div class="line"><a name="l00039"></a><span class="lineno"> 39</span> };</div>
-<div class="line"><a name="l00040"></a><span class="lineno"> 40</span> </div>
+<div class="line"><a name="l00032"></a><span class="lineno"> 32</span> </div>
+<div class="line"><a name="l00034"></a><span class="lineno"> 34</span>  <span class="keywordtype">double</span> getNearestLow(); </div>
+<div class="line"><a name="l00035"></a><span class="lineno"> 35</span> </div>
+<div class="line"><a name="l00037"></a><span class="lineno"> 37</span>  <span class="keywordtype">double</span> getNearestHigh();</div>
+<div class="line"><a name="l00038"></a><span class="lineno"> 38</span>  </div>
+<div class="line"><a name="l00039"></a><span class="lineno"> 39</span>  ClassDef(<a class="code" href="classshapeSys.html">shapeSys</a>,1)</div>
+<div class="line"><a name="l00040"></a><span class="lineno"> 40</span> };</div>
<div class="line"><a name="l00041"></a><span class="lineno"> 41</span> </div>
-<div class="line"><a name="l00042"></a><span class="lineno"><a class="line" href="classscaleSys.html"> 42</a></span> <span class="keyword">class </span><a class="code" href="classscaleSys.html">scaleSys</a> : <span class="keyword">public</span> <a class="code" href="classLKParameter.html">LKParameter</a> {</div>
-<div class="line"><a name="l00043"></a><span class="lineno"> 43</span> </div>
-<div class="line"><a name="l00044"></a><span class="lineno"> 44</span>  <span class="keyword">public</span>:</div>
-<div class="line"><a name="l00045"></a><span class="lineno"> 45</span>  <span class="comment">// relativeUncertainty = 1 sigma relative uncertainty (example 10% uncertainty would be 0.1) </span></div>
-<div class="line"><a name="l00046"></a><span class="lineno"> 46</span>  <a class="code" href="classscaleSys.html">scaleSys</a>(TString name, <span class="keywordtype">double</span> relativeUncertainty);</div>
-<div class="line"><a name="l00048"></a><span class="lineno"> 48</span> </div>
-<div class="line"><a name="l00049"></a><span class="lineno"> 49</span>  <span class="comment">// gives the Normalization scale factor according to t-valued uncertainty --> (1 + t_val * relativeUncertainty ) </span></div>
-<div class="line"><a name="l00050"></a><span class="lineno"> 50</span>  <span class="comment">// Uses the current t-value which is automatically linked to its pdfLikelihood </span></div>
-<div class="line"><a name="l00051"></a><span class="lineno"> 51</span>  <span class="keywordtype">double</span> getNormModifier();</div>
-<div class="line"><a name="l00052"></a><span class="lineno"> 52</span> </div>
-<div class="line"><a name="l00053"></a><span class="lineno"> 53</span>  </div>
-<div class="line"><a name="l00054"></a><span class="lineno"> 54</span>  <span class="comment">//this case return a sys that is centered in zero, half a gaussian, t-value strictly positive </span></div>
-<div class="line"><a name="l00055"></a><span class="lineno"> 55</span>  <span class="comment">//and for a tvalue=0 have zero events. Histo is supposed to be normalized to 1</span></div>
-<div class="line"><a name="l00056"></a><span class="lineno"> 56</span>  <span class="keywordtype">void</span> setNull(){ isNull = <span class="keyword">true</span>; setMinimum(0.);};</div>
-<div class="line"><a name="l00057"></a><span class="lineno"> 57</span> </div>
-<div class="line"><a name="l00058"></a><span class="lineno"> 58</span>  <span class="keyword">private</span>:</div>
-<div class="line"><a name="l00059"></a><span class="lineno"> 59</span>  <span class="keywordtype">double</span> relUnc;</div>
-<div class="line"><a name="l00060"></a><span class="lineno"> 60</span>  <span class="keywordtype">bool</span> isNull;</div>
-<div class="line"><a name="l00061"></a><span class="lineno"> 61</span>  </div>
-<div class="line"><a name="l00062"></a><span class="lineno"> 62</span> </div>
-<div class="line"><a name="l00063"></a><span class="lineno"> 63</span>  ClassDef(<a class="code" href="classscaleSys.html">scaleSys</a>,1)</div>
-<div class="line"><a name="l00064"></a><span class="lineno"> 64</span> };</div>
-<div class="line"><a name="l00065"></a><span class="lineno"> 65</span> </div>
+<div class="line"><a name="l00042"></a><span class="lineno"> 42</span> </div>
+<div class="line"><a name="l00043"></a><span class="lineno"><a class="line" href="classscaleSys.html"> 43</a></span> <span class="keyword">class </span><a class="code" href="classscaleSys.html">scaleSys</a> : <span class="keyword">public</span> <a class="code" href="classLKParameter.html">LKParameter</a> {</div>
+<div class="line"><a name="l00044"></a><span class="lineno"> 44</span> </div>
+<div class="line"><a name="l00045"></a><span class="lineno"> 45</span>  <span class="keyword">public</span>:</div>
+<div class="line"><a name="l00046"></a><span class="lineno"> 46</span>  <span class="comment">// relativeUncertainty = 1 sigma relative uncertainty (example 10% uncertainty would be 0.1) </span></div>
+<div class="line"><a name="l00047"></a><span class="lineno"> 47</span>  <a class="code" href="classscaleSys.html">scaleSys</a>(TString name, <span class="keywordtype">double</span> relativeUncertainty);</div>
+<div class="line"><a name="l00049"></a><span class="lineno"> 49</span> </div>
+<div class="line"><a name="l00050"></a><span class="lineno"> 50</span>  <span class="comment">// gives the Normalization scale factor according to t-valued uncertainty --> (1 + t_val * relativeUncertainty ) </span></div>
+<div class="line"><a name="l00051"></a><span class="lineno"> 51</span>  <span class="comment">// Uses the current t-value which is automatically linked to its pdfLikelihood </span></div>
+<div class="line"><a name="l00052"></a><span class="lineno"> 52</span>  <span class="keywordtype">double</span> getNormModifier();</div>
+<div class="line"><a name="l00053"></a><span class="lineno"> 53</span> </div>
+<div class="line"><a name="l00054"></a><span class="lineno"> 54</span>  </div>
+<div class="line"><a name="l00055"></a><span class="lineno"> 55</span>  <span class="comment">//this case return a sys that is centered in zero, half a gaussian, t-value strictly positive </span></div>
+<div class="line"><a name="l00056"></a><span class="lineno"> 56</span>  <span class="comment">//and for a tvalue=0 have zero events. Histo is supposed to be normalized to 1</span></div>
+<div class="line"><a name="l00057"></a><span class="lineno"> 57</span>  <span class="keywordtype">void</span> setNull(){ isNull = <span class="keyword">true</span>; setMinimum(0.);};</div>
+<div class="line"><a name="l00058"></a><span class="lineno"> 58</span> </div>
+<div class="line"><a name="l00059"></a><span class="lineno"> 59</span>  <span class="keyword">private</span>:</div>
+<div class="line"><a name="l00060"></a><span class="lineno"> 60</span>  <span class="keywordtype">double</span> relUnc;</div>
+<div class="line"><a name="l00061"></a><span class="lineno"> 61</span>  <span class="keywordtype">bool</span> isNull;</div>
+<div class="line"><a name="l00062"></a><span class="lineno"> 62</span>  </div>
+<div class="line"><a name="l00063"></a><span class="lineno"> 63</span> </div>
+<div class="line"><a name="l00064"></a><span class="lineno"> 64</span>  ClassDef(<a class="code" href="classscaleSys.html">scaleSys</a>,1)</div>
+<div class="line"><a name="l00065"></a><span class="lineno"> 65</span> };</div>
<div class="line"><a name="l00066"></a><span class="lineno"> 66</span> </div>
-<div class="line"><a name="l00067"></a><span class="lineno"><a class="line" href="classpdfComponent.html"> 67</a></span> <span class="keyword">class </span><a class="code" href="classpdfComponent.html">pdfComponent</a> :<a class="code" href="classerrorHandler.html">errorHandler</a>{</div>
-<div class="line"><a name="l00068"></a><span class="lineno"> 68</span> </div>
-<div class="line"><a name="l00069"></a><span class="lineno"> 69</span>  <span class="keyword">public</span>:</div>
-<div class="line"><a name="l00070"></a><span class="lineno"> 70</span> </div>
-<div class="line"><a name="l00071"></a><span class="lineno"> 71</span>  <a class="code" href="classpdfComponent.html">pdfComponent</a>(TString name, TString filename);</div>
-<div class="line"><a name="l00072"></a><span class="lineno"> 72</span>  ~<a class="code" href="classpdfComponent.html">pdfComponent</a>();</div>
-<div class="line"><a name="l00073"></a><span class="lineno"> 73</span> </div>
+<div class="line"><a name="l00067"></a><span class="lineno"> 67</span> </div>
+<div class="line"><a name="l00068"></a><span class="lineno"><a class="line" href="classpdfComponent.html"> 68</a></span> <span class="keyword">class </span><a class="code" href="classpdfComponent.html">pdfComponent</a> :<a class="code" href="classerrorHandler.html">errorHandler</a>{</div>
+<div class="line"><a name="l00069"></a><span class="lineno"> 69</span> </div>
+<div class="line"><a name="l00070"></a><span class="lineno"> 70</span>  <span class="keyword">public</span>:</div>
+<div class="line"><a name="l00071"></a><span class="lineno"> 71</span> </div>
+<div class="line"><a name="l00072"></a><span class="lineno"> 72</span>  <a class="code" href="classpdfComponent.html">pdfComponent</a>(TString name, TString filename);</div>
+<div class="line"><a name="l00073"></a><span class="lineno"> 73</span>  ~<a class="code" href="classpdfComponent.html">pdfComponent</a>();</div>
<div class="line"><a name="l00074"></a><span class="lineno"> 74</span> </div>
-<div class="line"><a name="l00075"></a><span class="lineno"> 75</span>  <span class="keywordtype">void</span> addScaleSys(<a class="code" href="classscaleSys.html">scaleSys</a> *addMe) { myScaleUnc.push_back(addMe); };</div>
-<div class="line"><a name="l00076"></a><span class="lineno"> 76</span> </div>
-<div class="line"><a name="l00077"></a><span class="lineno"> 77</span>  <span class="keywordtype">void</span> addShapeSys(<a class="code" href="classshapeSys.html">shapeSys</a> *addMe) { myShapeUnc.push_back(addMe); };</div>
-<div class="line"><a name="l00078"></a><span class="lineno"> 78</span> </div>
-<div class="line"><a name="l00079"></a><span class="lineno"> 79</span>  <span class="comment">//load histogram according to the current value of the parameters</span></div>
-<div class="line"><a name="l00080"></a><span class="lineno"> 80</span>  <span class="keywordtype">void</span> loadHistos();</div>
-<div class="line"><a name="l00081"></a><span class="lineno"> 81</span> </div>
-<div class="line"><a name="l00082"></a><span class="lineno"> 82</span>  <span class="comment">//load default histogram, no sys.</span></div>
-<div class="line"><a name="l00083"></a><span class="lineno"> 83</span>  <span class="keywordtype">void</span> loadDefaultHisto();</div>
-<div class="line"><a name="l00084"></a><span class="lineno"> 84</span> </div>
-<div class="line"><a name="l00085"></a><span class="lineno"> 85</span>  <span class="comment">// returns the interpolated pdf over shape sys computed in (s1,s2),</span></div>
-<div class="line"><a name="l00086"></a><span class="lineno"> 86</span>  <span class="comment">// it is normalized to number of events and scale uncertainty are also </span></div>
-<div class="line"><a name="l00087"></a><span class="lineno"> 87</span>  <span class="comment">// taken into account</span></div>
-<div class="line"><a name="l00088"></a><span class="lineno"> 88</span>  <span class="keywordtype">double</span> getNormalizedDensity(<span class="keywordtype">double</span> s1, <span class="keywordtype">double</span> s2);</div>
-<div class="line"><a name="l00089"></a><span class="lineno"> 89</span> </div>
-<div class="line"><a name="l00090"></a><span class="lineno"> 90</span>  <span class="comment">// returns the (s1,s2) bin content of the default histo, no shape nor scale sys is applied</span></div>
-<div class="line"><a name="l00091"></a><span class="lineno"> 91</span>  <span class="keywordtype">double</span> getDefaultDensity(<span class="keywordtype">double</span> s1, <span class="keywordtype">double</span> s2);</div>
-<div class="line"><a name="l00092"></a><span class="lineno"> 92</span> </div>
-<div class="line"><a name="l00093"></a><span class="lineno"> 93</span>  <span class="comment">// Returns the total integrated number of events taking into account shape sys</span></div>
-<div class="line"><a name="l00094"></a><span class="lineno"> 94</span>  <span class="comment">// and scale sys.</span></div>
-<div class="line"><a name="l00095"></a><span class="lineno"> 95</span>  <span class="keywordtype">double</span> getNormalizedEvents();</div>
-<div class="line"><a name="l00096"></a><span class="lineno"> 96</span>  </div>
-<div class="line"><a name="l00097"></a><span class="lineno"> 97</span>  <span class="comment">// returns total integral of the default histo, no shape nor scale sys is applied</span></div>
-<div class="line"><a name="l00098"></a><span class="lineno"> 98</span>  <span class="keywordtype">double</span> getDefaultEvents();</div>
-<div class="line"><a name="l00099"></a><span class="lineno"> 99</span> </div>
-<div class="line"><a name="l00100"></a><span class="lineno"> 100</span>  <span class="comment">// Returns a copy of the interpolated histogram according to the "currentValue" of the</span></div>
-<div class="line"><a name="l00101"></a><span class="lineno"> 101</span>  <span class="comment">// shapeSys associated with it. The current value of each uncertainty can</span></div>
-<div class="line"><a name="l00102"></a><span class="lineno"> 102</span>  <span class="comment">// be set simply by shapeSys::setCurrentValue(val) method.</span></div>
-<div class="line"><a name="l00103"></a><span class="lineno"> 103</span>  TH2F getInterpolatedHisto();</div>
-<div class="line"><a name="l00104"></a><span class="lineno"> 104</span> </div>
-<div class="line"><a name="l00105"></a><span class="lineno"> 105</span>  <span class="comment">// Returns a copy of the default histogram, no shape nor scale sys is applied</span></div>
-<div class="line"><a name="l00106"></a><span class="lineno"> 106</span>  TH2F getDefaultHisto();</div>
-<div class="line"><a name="l00107"></a><span class="lineno"> 107</span> </div>
-<div class="line"><a name="l00108"></a><span class="lineno"> 108</span>  <span class="comment">// returns the grid points in histogram space that needs to be loaded</span></div>
-<div class="line"><a name="l00109"></a><span class="lineno"> 109</span>  <span class="comment">// this method need to be modified for arbitrary number of shape sys </span></div>
-<div class="line"><a name="l00110"></a><span class="lineno"> 110</span>  TString getNearestHistoName(vector<bool> setOfVal);</div>
-<div class="line"><a name="l00111"></a><span class="lineno"> 111</span> </div>
-<div class="line"><a name="l00112"></a><span class="lineno"> 112</span>  <span class="comment">// returns the default name of the histogram</span></div>
-<div class="line"><a name="l00113"></a><span class="lineno"> 113</span>  TString getDefaultHistoName();</div>
-<div class="line"><a name="l00114"></a><span class="lineno"> 114</span> </div>
-<div class="line"><a name="l00115"></a><span class="lineno"> 115</span>  <span class="comment">// returns a string with current values of parameters</span></div>
-<div class="line"><a name="l00116"></a><span class="lineno"> 116</span>  TString getParamValueString();</div>
-<div class="line"><a name="l00117"></a><span class="lineno"> 117</span> </div>
-<div class="line"><a name="l00118"></a><span class="lineno"> 118</span>  TString getParamValueWritable();</div>
-<div class="line"><a name="l00119"></a><span class="lineno"> 119</span> </div>
-<div class="line"><a name="l00120"></a><span class="lineno"> 120</span>  <span class="comment">// returns the integral of the square between points (not bins), linearly interpolates</span></div>
-<div class="line"><a name="l00121"></a><span class="lineno"> 121</span>  <span class="comment">// the area if desn't correspond to an integer number of bins.</span></div>
-<div class="line"><a name="l00122"></a><span class="lineno"> 122</span>  <span class="comment">// Uses the default histo, does not consider shape or scale uncertainty.</span></div>
-<div class="line"><a name="l00123"></a><span class="lineno"> 123</span>  <span class="keywordtype">double</span> getDefaultPdfIntegral(<span class="keywordtype">double</span> s1_min, <span class="keywordtype">double</span> s1_max, <span class="keywordtype">double</span> s2_min, <span class="keywordtype">double</span> s2_max);</div>
-<div class="line"><a name="l00124"></a><span class="lineno"> 124</span> </div>
-<div class="line"><a name="l00125"></a><span class="lineno"> 125</span>  <span class="comment">//This scans the parameter space given a number of steps for each parameter</span></div>
-<div class="line"><a name="l00126"></a><span class="lineno"> 126</span>  <span class="comment">//produce a projection given min and max, produce a PDF file for comparison</span></div>
-<div class="line"><a name="l00127"></a><span class="lineno"> 127</span>  <span class="comment">//of all interpolation and saves also the TH2F to file.</span></div>
-<div class="line"><a name="l00128"></a><span class="lineno"> 128</span>  <span class="keywordtype">void</span> plotInterpolatedSpace(<span class="keywordtype">bool</span> doProjectionX, <span class="keywordtype">double</span> min, <span class="keywordtype">double</span> max, <span class="keywordtype">int</span> Nsteps);</div>
-<div class="line"><a name="l00129"></a><span class="lineno"> 129</span>  <span class="comment">//scale the pdf by VAL, this happend for the methods:</span></div>
-<div class="line"><a name="l00130"></a><span class="lineno"> 130</span>  <span class="comment">//getDefaultHisto, getInterpolatedHisto, getNormalizedDensity, </span></div>
-<div class="line"><a name="l00131"></a><span class="lineno"> 131</span>  <span class="comment">//getDefaultDensity, getNormalizedEvents, getDefaultEvents</span></div>
-<div class="line"><a name="l00132"></a><span class="lineno"> 132</span>  <span class="keywordtype">void</span> setScaleFactor(<span class="keywordtype">double</span> val) {scaleFactor = val;}; </div>
-<div class="line"><a name="l00133"></a><span class="lineno"> 133</span> </div>
-<div class="line"><a name="l00134"></a><span class="lineno"> 134</span>  TString getName() {<span class="keywordflow">return</span> pdf_name;};</div>
-<div class="line"><a name="l00135"></a><span class="lineno"> 135</span> </div>
-<div class="line"><a name="l00136"></a><span class="lineno"> 136</span>  vector<scaleSys*> myScaleUnc; <span class="comment">//container of scale sys</span></div>
-<div class="line"><a name="l00137"></a><span class="lineno"> 137</span>  vector<shapeSys*> myShapeUnc; <span class="comment">//container of shape sys</span></div>
-<div class="line"><a name="l00138"></a><span class="lineno"> 138</span> </div>
-<div class="line"><a name="l00139"></a><span class="lineno"> 139</span>  TString suffix;</div>
-<div class="line"><a name="l00140"></a><span class="lineno"> 140</span>  <span class="keywordtype">bool</span> doExtend;</div>
-<div class="line"><a name="l00141"></a><span class="lineno"> 141</span> </div>
-<div class="line"><a name="l00142"></a><span class="lineno"> 142</span>  <span class="keyword">private</span>:</div>
-<div class="line"><a name="l00143"></a><span class="lineno"> 143</span>  TFile *file;</div>
-<div class="line"><a name="l00144"></a><span class="lineno"> 144</span>  TH2F *defaultDistro;</div>
-<div class="line"><a name="l00145"></a><span class="lineno"> 145</span>  vector<TH2F*> histos; <span class="comment">// contains the 2^N histo for the hyperplane interpolation of shapeSys</span></div>
-<div class="line"><a name="l00146"></a><span class="lineno"> 146</span>  vector<double> InterpFactors; <span class="comment">// contains the 2^N interpolation factors for the hyperplane interpolation of shapeSys </span></div>
-<div class="line"><a name="l00147"></a><span class="lineno"> 147</span>  TString pdf_name;</div>
-<div class="line"><a name="l00148"></a><span class="lineno"> 148</span>  vector<double> old_t_val; <span class="comment">//contains the last value interpolated, the interpolation is lazy, doesn't ricompute it if is for the same set of values.</span></div>
-<div class="line"><a name="l00149"></a><span class="lineno"> 149</span>  <span class="keywordtype">double</span> scaleFactor;</div>
-<div class="line"><a name="l00150"></a><span class="lineno"> 150</span> </div>
+<div class="line"><a name="l00075"></a><span class="lineno"> 75</span> </div>
+<div class="line"><a name="l00076"></a><span class="lineno"> 76</span>  <span class="keywordtype">void</span> addScaleSys(<a class="code" href="classscaleSys.html">scaleSys</a> *addMe) { myScaleUnc.push_back(addMe); };</div>
+<div class="line"><a name="l00077"></a><span class="lineno"> 77</span> </div>
+<div class="line"><a name="l00078"></a><span class="lineno"> 78</span>  <span class="keywordtype">void</span> addShapeSys(<a class="code" href="classshapeSys.html">shapeSys</a> *addMe) { myShapeUnc.push_back(addMe); };</div>
+<div class="line"><a name="l00079"></a><span class="lineno"> 79</span> </div>
+<div class="line"><a name="l00080"></a><span class="lineno"> 80</span>  <span class="comment">//load histogram according to the current value of the parameters</span></div>
+<div class="line"><a name="l00081"></a><span class="lineno"> 81</span>  <span class="keywordtype">void</span> loadHistos();</div>
+<div class="line"><a name="l00082"></a><span class="lineno"> 82</span> </div>
+<div class="line"><a name="l00083"></a><span class="lineno"> 83</span>  <span class="comment">//load default histogram, no sys.</span></div>
+<div class="line"><a name="l00084"></a><span class="lineno"> 84</span>  <span class="keywordtype">void</span> loadDefaultHisto();</div>
+<div class="line"><a name="l00085"></a><span class="lineno"> 85</span> </div>
+<div class="line"><a name="l00086"></a><span class="lineno"> 86</span>  <span class="comment">// returns the interpolated pdf over shape sys computed in (s1,s2),</span></div>
+<div class="line"><a name="l00087"></a><span class="lineno"> 87</span>  <span class="comment">// it is normalized to number of events and scale uncertainty are also </span></div>
+<div class="line"><a name="l00088"></a><span class="lineno"> 88</span>  <span class="comment">// taken into account</span></div>
+<div class="line"><a name="l00089"></a><span class="lineno"> 89</span>  <span class="keywordtype">double</span> getNormalizedDensity(<span class="keywordtype">double</span> s1, <span class="keywordtype">double</span> s2);</div>
+<div class="line"><a name="l00090"></a><span class="lineno"> 90</span> </div>
+<div class="line"><a name="l00091"></a><span class="lineno"> 91</span>  <span class="comment">// returns the (s1,s2) bin content of the default histo, no shape nor scale sys is applied</span></div>
+<div class="line"><a name="l00092"></a><span class="lineno"> 92</span>  <span class="keywordtype">double</span> getDefaultDensity(<span class="keywordtype">double</span> s1, <span class="keywordtype">double</span> s2);</div>
+<div class="line"><a name="l00093"></a><span class="lineno"> 93</span> </div>
+<div class="line"><a name="l00094"></a><span class="lineno"> 94</span>  <span class="comment">// Returns the total integrated number of events taking into account shape sys</span></div>
+<div class="line"><a name="l00095"></a><span class="lineno"> 95</span>  <span class="comment">// and scale sys.</span></div>
+<div class="line"><a name="l00096"></a><span class="lineno"> 96</span>  <span class="keywordtype">double</span> getNormalizedEvents();</div>
+<div class="line"><a name="l00097"></a><span class="lineno"> 97</span>  </div>
+<div class="line"><a name="l00098"></a><span class="lineno"> 98</span>  <span class="comment">// returns total integral of the default histo, no shape nor scale sys is applied</span></div>
+<div class="line"><a name="l00099"></a><span class="lineno"> 99</span>  <span class="keywordtype">double</span> getDefaultEvents();</div>
+<div class="line"><a name="l00100"></a><span class="lineno"> 100</span> </div>
+<div class="line"><a name="l00101"></a><span class="lineno"> 101</span>  <span class="comment">// Returns a copy of the interpolated histogram according to the "currentValue" of the</span></div>
+<div class="line"><a name="l00102"></a><span class="lineno"> 102</span>  <span class="comment">// shapeSys associated with it. The current value of each uncertainty can</span></div>
+<div class="line"><a name="l00103"></a><span class="lineno"> 103</span>  <span class="comment">// be set simply by shapeSys::setCurrentValue(val) method.</span></div>
+<div class="line"><a name="l00104"></a><span class="lineno"> 104</span>  TH2F getInterpolatedHisto();</div>
+<div class="line"><a name="l00105"></a><span class="lineno"> 105</span> </div>
+<div class="line"><a name="l00106"></a><span class="lineno"> 106</span>  <span class="comment">// Returns a copy of the default histogram, no shape nor scale sys is applied</span></div>
+<div class="line"><a name="l00107"></a><span class="lineno"> 107</span>  TH2F getDefaultHisto();</div>
+<div class="line"><a name="l00108"></a><span class="lineno"> 108</span> </div>
+<div class="line"><a name="l00109"></a><span class="lineno"> 109</span>  <span class="comment">// returns the grid points in histogram space that needs to be loaded</span></div>
+<div class="line"><a name="l00110"></a><span class="lineno"> 110</span>  <span class="comment">// this method need to be modified for arbitrary number of shape sys </span></div>
+<div class="line"><a name="l00111"></a><span class="lineno"> 111</span>  TString getNearestHistoName(vector<bool> setOfVal);</div>
+<div class="line"><a name="l00112"></a><span class="lineno"> 112</span> </div>
+<div class="line"><a name="l00113"></a><span class="lineno"> 113</span>  <span class="comment">// returns the default name of the histogram</span></div>
+<div class="line"><a name="l00114"></a><span class="lineno"> 114</span>  TString getDefaultHistoName();</div>
+<div class="line"><a name="l00115"></a><span class="lineno"> 115</span> </div>
+<div class="line"><a name="l00116"></a><span class="lineno"> 116</span>  <span class="comment">// returns a string with current values of parameters</span></div>
+<div class="line"><a name="l00117"></a><span class="lineno"> 117</span>  TString getParamValueString();</div>
+<div class="line"><a name="l00118"></a><span class="lineno"> 118</span> </div>
+<div class="line"><a name="l00119"></a><span class="lineno"> 119</span>  TString getParamValueWritable();</div>
+<div class="line"><a name="l00120"></a><span class="lineno"> 120</span> </div>
+<div class="line"><a name="l00121"></a><span class="lineno"> 121</span>  <span class="comment">// returns the integral of the square between points (not bins), linearly interpolates</span></div>
+<div class="line"><a name="l00122"></a><span class="lineno"> 122</span>  <span class="comment">// the area if desn't correspond to an integer number of bins.</span></div>
+<div class="line"><a name="l00123"></a><span class="lineno"> 123</span>  <span class="comment">// Uses the default histo, does not consider shape or scale uncertainty.</span></div>
+<div class="line"><a name="l00124"></a><span class="lineno"> 124</span>  <span class="keywordtype">double</span> getDefaultPdfIntegral(<span class="keywordtype">double</span> s1_min, <span class="keywordtype">double</span> s1_max, <span class="keywordtype">double</span> s2_min, <span class="keywordtype">double</span> s2_max);</div>
+<div class="line"><a name="l00125"></a><span class="lineno"> 125</span> </div>
+<div class="line"><a name="l00126"></a><span class="lineno"> 126</span>  <span class="comment">//This scans the parameter space given a number of steps for each parameter</span></div>
+<div class="line"><a name="l00127"></a><span class="lineno"> 127</span>  <span class="comment">//produce a projection given min and max, produce a PDF file for comparison</span></div>
+<div class="line"><a name="l00128"></a><span class="lineno"> 128</span>  <span class="comment">//of all interpolation and saves also the TH2F to file.</span></div>
+<div class="line"><a name="l00129"></a><span class="lineno"> 129</span>  <span class="keywordtype">void</span> plotInterpolatedSpace(<span class="keywordtype">bool</span> doProjectionX, <span class="keywordtype">double</span> min, <span class="keywordtype">double</span> max, <span class="keywordtype">int</span> Nsteps);</div>
+<div class="line"><a name="l00130"></a><span class="lineno"> 130</span>  <span class="comment">//scale the pdf by VAL, this happend for the methods:</span></div>
+<div class="line"><a name="l00131"></a><span class="lineno"> 131</span>  <span class="comment">//getDefaultHisto, getInterpolatedHisto, getNormalizedDensity, </span></div>
+<div class="line"><a name="l00132"></a><span class="lineno"> 132</span>  <span class="comment">//getDefaultDensity, getNormalizedEvents, getDefaultEvents</span></div>
+<div class="line"><a name="l00133"></a><span class="lineno"> 133</span>  <span class="keywordtype">void</span> setScaleFactor(<span class="keywordtype">double</span> val) {scaleFactor = val;}; </div>
+<div class="line"><a name="l00134"></a><span class="lineno"> 134</span> </div>
+<div class="line"><a name="l00135"></a><span class="lineno"> 135</span>  TString getName() {<span class="keywordflow">return</span> pdf_name;};</div>
+<div class="line"><a name="l00136"></a><span class="lineno"> 136</span> </div>
+<div class="line"><a name="l00137"></a><span class="lineno"> 137</span>  vector<scaleSys*> myScaleUnc; <span class="comment">//container of scale sys</span></div>
+<div class="line"><a name="l00138"></a><span class="lineno"> 138</span>  vector<shapeSys*> myShapeUnc; <span class="comment">//container of shape sys</span></div>
+<div class="line"><a name="l00139"></a><span class="lineno"> 139</span> </div>
+<div class="line"><a name="l00140"></a><span class="lineno"> 140</span>  TString suffix;</div>
+<div class="line"><a name="l00141"></a><span class="lineno"> 141</span>  <span class="keywordtype">bool</span> doExtend;</div>
+<div class="line"><a name="l00142"></a><span class="lineno"> 142</span> </div>
+<div class="line"><a name="l00143"></a><span class="lineno"> 143</span>  <span class="keyword">private</span>:</div>
+<div class="line"><a name="l00144"></a><span class="lineno"> 144</span>  TFile *file;</div>
+<div class="line"><a name="l00145"></a><span class="lineno"> 145</span>  TH2F *defaultDistro;</div>
+<div class="line"><a name="l00146"></a><span class="lineno"> 146</span>  vector<TH2F*> histos; <span class="comment">// contains the 2^N histo for the hyperplane interpolation of shapeSys</span></div>
+<div class="line"><a name="l00147"></a><span class="lineno"> 147</span>  vector<double> InterpFactors; <span class="comment">// contains the 2^N interpolation factors for the hyperplane interpolation of shapeSys </span></div>
+<div class="line"><a name="l00148"></a><span class="lineno"> 148</span>  TString pdf_name;</div>
+<div class="line"><a name="l00149"></a><span class="lineno"> 149</span>  vector<double> old_t_val; <span class="comment">//contains the last value interpolated, the interpolation is lazy, doesn't ricompute it if is for the same set of values.</span></div>
+<div class="line"><a name="l00150"></a><span class="lineno"> 150</span>  <span class="keywordtype">double</span> scaleFactor;</div>
<div class="line"><a name="l00151"></a><span class="lineno"> 151</span> </div>
-<div class="line"><a name="l00152"></a><span class="lineno"> 152</span>  <span class="keywordtype">void</span> extendHisto(TH2F &h);</div>
-<div class="line"><a name="l00153"></a><span class="lineno"> 153</span> </div>
-<div class="line"><a name="l00154"></a><span class="lineno"> 154</span>  ClassDef(<a class="code" href="classpdfComponent.html">pdfComponent</a>,1)</div>
-<div class="line"><a name="l00155"></a><span class="lineno"> 155</span> };</div>
-<div class="line"><a name="l00156"></a><span class="lineno"> 156</span> </div>
+<div class="line"><a name="l00152"></a><span class="lineno"> 152</span> </div>
+<div class="line"><a name="l00153"></a><span class="lineno"> 153</span>  <span class="keywordtype">void</span> extendHisto(TH2F &h);</div>
+<div class="line"><a name="l00154"></a><span class="lineno"> 154</span> </div>
+<div class="line"><a name="l00155"></a><span class="lineno"> 155</span>  ClassDef(<a class="code" href="classpdfComponent.html">pdfComponent</a>,1)</div>
+<div class="line"><a name="l00156"></a><span class="lineno"> 156</span> };</div>
<div class="line"><a name="l00157"></a><span class="lineno"> 157</span> </div>
<div class="line"><a name="l00158"></a><span class="lineno"> 158</span> </div>
<div class="line"><a name="l00159"></a><span class="lineno"> 159</span> </div>
-<div class="line"><a name="l00160"></a><span class="lineno"><a class="line" href="classhistoCompare.html"> 160</a></span> <span class="keyword">class </span><a class="code" href="classhistoCompare.html">histoCompare</a> :<a class="code" href="classerrorHandler.html">errorHandler</a>{</div>
-<div class="line"><a name="l00161"></a><span class="lineno"> 161</span> </div>
-<div class="line"><a name="l00162"></a><span class="lineno"> 162</span>  <span class="comment">// This class is supposed to compare a "base" histo </span></div>
-<div class="line"><a name="l00163"></a><span class="lineno"> 163</span>  <span class="comment">// with a set of other histogram. One can choose for ratio plot,</span></div>
-<div class="line"><a name="l00164"></a><span class="lineno"> 164</span>  <span class="comment">// stack plot of different component, one can set the projected axis</span></div>
-<div class="line"><a name="l00165"></a><span class="lineno"> 165</span>  <span class="comment">// and the slice on which do the projection, one can also rebin the histo.</span></div>
-<div class="line"><a name="l00166"></a><span class="lineno"> 166</span>  <span class="comment">// the histo are supposed to be TH2F.</span></div>
-<div class="line"><a name="l00167"></a><span class="lineno"> 167</span> </div>
-<div class="line"><a name="l00168"></a><span class="lineno"> 168</span>  <span class="keyword">public</span>:</div>
-<div class="line"><a name="l00169"></a><span class="lineno"> 169</span> </div>
-<div class="line"><a name="l00170"></a><span class="lineno"> 170</span>  <a class="code" href="classhistoCompare.html">histoCompare</a>();</div>
-<div class="line"><a name="l00171"></a><span class="lineno"> 171</span>  ~<a class="code" href="classhistoCompare.html">histoCompare</a>();</div>
-<div class="line"><a name="l00172"></a><span class="lineno"> 172</span> </div>
-<div class="line"><a name="l00173"></a><span class="lineno"> 173</span>  <span class="comment">//---------OPTIONS-----------// </span></div>
-<div class="line"><a name="l00174"></a><span class="lineno"> 174</span>  <span class="keywordtype">int</span> rebinX;</div>
-<div class="line"><a name="l00175"></a><span class="lineno"> 175</span>  <span class="keywordtype">int</span> rebinY;</div>
-<div class="line"><a name="l00176"></a><span class="lineno"> 176</span>  <span class="keywordtype">bool</span> projectionX;</div>
-<div class="line"><a name="l00177"></a><span class="lineno"> 177</span>  <span class="keywordtype">double</span> projectionMin; <span class="comment">// minimum value on the projected axis</span></div>
-<div class="line"><a name="l00178"></a><span class="lineno"> 178</span>  <span class="keywordtype">double</span> projectionMax; <span class="comment">// maximum value on the projected axis</span></div>
-<div class="line"><a name="l00179"></a><span class="lineno"> 179</span>  <span class="keywordtype">int</span> binMin; <span class="comment">// same but in bin nuber</span></div>
-<div class="line"><a name="l00180"></a><span class="lineno"> 180</span>  <span class="keywordtype">int</span> binMax; </div>
-<div class="line"><a name="l00181"></a><span class="lineno"> 181</span>  <span class="keywordtype">bool</span> doStack; <span class="comment">// True if you want to sum up histo in compareList</span></div>
-<div class="line"><a name="l00182"></a><span class="lineno"> 182</span>  TString titleX;</div>
-<div class="line"><a name="l00183"></a><span class="lineno"> 183</span>  TString titleY;</div>
-<div class="line"><a name="l00184"></a><span class="lineno"> 184</span>  <span class="comment">//---------------------------//</span></div>
-<div class="line"><a name="l00185"></a><span class="lineno"> 185</span> </div>
-<div class="line"><a name="l00186"></a><span class="lineno"> 186</span>  <span class="comment">// set histo with wich one one to compare the rest</span></div>
-<div class="line"><a name="l00187"></a><span class="lineno"> 187</span>  <span class="keywordtype">void</span> setBaseHisto(TH2F b, TString n=<span class="stringliteral">""</span>){base = b; names[0] = n;};</div>
-<div class="line"><a name="l00188"></a><span class="lineno"> 188</span>  <span class="comment">// set the comparison histos</span></div>
-<div class="line"><a name="l00189"></a><span class="lineno"> 189</span>  <span class="keywordtype">void</span> addHistoToList(TH2F h, TString n =<span class="stringliteral">""</span>) {compareList.push_back(h); names.push_back(n);};</div>
-<div class="line"><a name="l00190"></a><span class="lineno"> 190</span>  <span class="comment">//set Name of component</span></div>
-<div class="line"><a name="l00191"></a><span class="lineno"> 191</span>  <span class="keywordtype">void</span> setNameofComponent(<span class="keywordtype">unsigned</span> <span class="keywordtype">int</span> i, TString n);</div>
-<div class="line"><a name="l00192"></a><span class="lineno"> 192</span>  <span class="comment">// draw just a normal comparison plot</span></div>
-<div class="line"><a name="l00193"></a><span class="lineno"> 193</span>  <span class="keywordtype">void</span> compare();</div>
-<div class="line"><a name="l00194"></a><span class="lineno"> 194</span>  <span class="comment">// add a ratio plot to compare()</span></div>
-<div class="line"><a name="l00195"></a><span class="lineno"> 195</span>  <span class="keywordtype">void</span> compareWithRatio();</div>
-<div class="line"><a name="l00196"></a><span class="lineno"> 196</span> </div>
-<div class="line"><a name="l00197"></a><span class="lineno"> 197</span>  <span class="keywordtype">void</span> drawLegend(TH1D *baseH, vector <TH1D*> list );</div>
-<div class="line"><a name="l00198"></a><span class="lineno"> 198</span>  </div>
-<div class="line"><a name="l00199"></a><span class="lineno"> 199</span>  <span class="comment">//return a string of info</span></div>
-<div class="line"><a name="l00200"></a><span class="lineno"> 200</span>  TString projectionInfo();</div>
-<div class="line"><a name="l00201"></a><span class="lineno"> 201</span> </div>
-<div class="line"><a name="l00202"></a><span class="lineno"> 202</span>  <span class="comment">//------ histo holders --------//</span></div>
-<div class="line"><a name="l00203"></a><span class="lineno"> 203</span>  TH2F base;</div>
-<div class="line"><a name="l00204"></a><span class="lineno"> 204</span>  vector <TString> names; </div>
-<div class="line"><a name="l00205"></a><span class="lineno"> 205</span>  vector <TH2F> compareList;</div>
-<div class="line"><a name="l00206"></a><span class="lineno"> 206</span> </div>
-<div class="line"><a name="l00207"></a><span class="lineno"> 207</span>  TH1D *projectedBase;</div>
-<div class="line"><a name="l00208"></a><span class="lineno"> 208</span>  vector <TH1D*> projectedList;</div>
-<div class="line"><a name="l00209"></a><span class="lineno"> 209</span> </div>
-<div class="line"><a name="l00210"></a><span class="lineno"> 210</span>  <span class="keyword">private</span>:</div>
-<div class="line"><a name="l00211"></a><span class="lineno"> 211</span> </div>
-<div class="line"><a name="l00212"></a><span class="lineno"> 212</span>  <span class="keywordtype">void</span> setOptions(TH1D *h, <span class="keywordtype">bool</span> dataLike = <span class="keyword">true</span>, <span class="keywordtype">bool</span> isBottom = <span class="keyword">false</span>);</div>
-<div class="line"><a name="l00213"></a><span class="lineno"> 213</span>  <span class="keywordtype">void</span> project(); <span class="comment">// push histo into the TH1D</span></div>
-<div class="line"><a name="l00214"></a><span class="lineno"> 214</span> </div>
+<div class="line"><a name="l00160"></a><span class="lineno"> 160</span> </div>
+<div class="line"><a name="l00161"></a><span class="lineno"><a class="line" href="classhistoCompare.html"> 161</a></span> <span class="keyword">class </span><a class="code" href="classhistoCompare.html">histoCompare</a> :<a class="code" href="classerrorHandler.html">errorHandler</a>{</div>
+<div class="line"><a name="l00162"></a><span class="lineno"> 162</span> </div>
+<div class="line"><a name="l00163"></a><span class="lineno"> 163</span>  <span class="comment">// This class is supposed to compare a "base" histo </span></div>
+<div class="line"><a name="l00164"></a><span class="lineno"> 164</span>  <span class="comment">// with a set of other histogram. One can choose for ratio plot,</span></div>
+<div class="line"><a name="l00165"></a><span class="lineno"> 165</span>  <span class="comment">// stack plot of different component, one can set the projected axis</span></div>
+<div class="line"><a name="l00166"></a><span class="lineno"> 166</span>  <span class="comment">// and the slice on which do the projection, one can also rebin the histo.</span></div>
+<div class="line"><a name="l00167"></a><span class="lineno"> 167</span>  <span class="comment">// the histo are supposed to be TH2F.</span></div>
+<div class="line"><a name="l00168"></a><span class="lineno"> 168</span> </div>
+<div class="line"><a name="l00169"></a><span class="lineno"> 169</span>  <span class="keyword">public</span>:</div>
+<div class="line"><a name="l00170"></a><span class="lineno"> 170</span> </div>
+<div class="line"><a name="l00171"></a><span class="lineno"> 171</span>  <a class="code" href="classhistoCompare.html">histoCompare</a>();</div>
+<div class="line"><a name="l00172"></a><span class="lineno"> 172</span>  ~<a class="code" href="classhistoCompare.html">histoCompare</a>();</div>
+<div class="line"><a name="l00173"></a><span class="lineno"> 173</span> </div>
+<div class="line"><a name="l00174"></a><span class="lineno"> 174</span>  <span class="comment">//---------OPTIONS-----------// </span></div>
+<div class="line"><a name="l00175"></a><span class="lineno"> 175</span>  <span class="keywordtype">int</span> rebinX;</div>
+<div class="line"><a name="l00176"></a><span class="lineno"> 176</span>  <span class="keywordtype">int</span> rebinY;</div>
+<div class="line"><a name="l00177"></a><span class="lineno"> 177</span>  <span class="keywordtype">bool</span> projectionX;</div>
+<div class="line"><a name="l00178"></a><span class="lineno"> 178</span>  <span class="keywordtype">double</span> projectionMin; <span class="comment">// minimum value on the projected axis</span></div>
+<div class="line"><a name="l00179"></a><span class="lineno"> 179</span>  <span class="keywordtype">double</span> projectionMax; <span class="comment">// maximum value on the projected axis</span></div>
+<div class="line"><a name="l00180"></a><span class="lineno"> 180</span>  <span class="keywordtype">int</span> binMin; <span class="comment">// same but in bin nuber</span></div>
+<div class="line"><a name="l00181"></a><span class="lineno"> 181</span>  <span class="keywordtype">int</span> binMax; </div>
+<div class="line"><a name="l00182"></a><span class="lineno"> 182</span>  <span class="keywordtype">bool</span> doStack; <span class="comment">// True if you want to sum up histo in compareList</span></div>
+<div class="line"><a name="l00183"></a><span class="lineno"> 183</span>  TString titleX;</div>
+<div class="line"><a name="l00184"></a><span class="lineno"> 184</span>  TString titleY;</div>
+<div class="line"><a name="l00185"></a><span class="lineno"> 185</span>  <span class="comment">//---------------------------//</span></div>
+<div class="line"><a name="l00186"></a><span class="lineno"> 186</span> </div>
+<div class="line"><a name="l00187"></a><span class="lineno"> 187</span>  <span class="comment">// set histo with wich one one to compare the rest</span></div>
+<div class="line"><a name="l00188"></a><span class="lineno"> 188</span>  <span class="keywordtype">void</span> setBaseHisto(TH2F b, TString n=<span class="stringliteral">""</span>){base = b; names[0] = n;};</div>
+<div class="line"><a name="l00189"></a><span class="lineno"> 189</span>  <span class="comment">// set the comparison histos</span></div>
+<div class="line"><a name="l00190"></a><span class="lineno"> 190</span>  <span class="keywordtype">void</span> addHistoToList(TH2F h, TString n =<span class="stringliteral">""</span>) {compareList.push_back(h); names.push_back(n);};</div>
+<div class="line"><a name="l00191"></a><span class="lineno"> 191</span>  <span class="comment">//set Name of component</span></div>
+<div class="line"><a name="l00192"></a><span class="lineno"> 192</span>  <span class="keywordtype">void</span> setNameofComponent(<span class="keywordtype">unsigned</span> <span class="keywordtype">int</span> i, TString n);</div>
+<div class="line"><a name="l00193"></a><span class="lineno"> 193</span>  <span class="comment">// draw just a normal comparison plot</span></div>
+<div class="line"><a name="l00194"></a><span class="lineno"> 194</span>  <span class="keywordtype">void</span> compare();</div>
+<div class="line"><a name="l00195"></a><span class="lineno"> 195</span>  <span class="comment">// add a ratio plot to compare()</span></div>
+<div class="line"><a name="l00196"></a><span class="lineno"> 196</span>  <span class="keywordtype">void</span> compareWithRatio();</div>
+<div class="line"><a name="l00197"></a><span class="lineno"> 197</span> </div>
+<div class="line"><a name="l00198"></a><span class="lineno"> 198</span>  <span class="keywordtype">void</span> drawLegend(TH1D *baseH, vector <TH1D*> list );</div>
+<div class="line"><a name="l00199"></a><span class="lineno"> 199</span>  </div>
+<div class="line"><a name="l00200"></a><span class="lineno"> 200</span>  <span class="comment">//return a string of info</span></div>
+<div class="line"><a name="l00201"></a><span class="lineno"> 201</span>  TString projectionInfo();</div>
+<div class="line"><a name="l00202"></a><span class="lineno"> 202</span> </div>
+<div class="line"><a name="l00203"></a><span class="lineno"> 203</span>  <span class="comment">//------ histo holders --------//</span></div>
+<div class="line"><a name="l00204"></a><span class="lineno"> 204</span>  TH2F base;</div>
+<div class="line"><a name="l00205"></a><span class="lineno"> 205</span>  vector <TString> names; </div>
+<div class="line"><a name="l00206"></a><span class="lineno"> 206</span>  vector <TH2F> compareList;</div>
+<div class="line"><a name="l00207"></a><span class="lineno"> 207</span> </div>
+<div class="line"><a name="l00208"></a><span class="lineno"> 208</span>  TH1D *projectedBase;</div>
+<div class="line"><a name="l00209"></a><span class="lineno"> 209</span>  vector <TH1D*> projectedList;</div>
+<div class="line"><a name="l00210"></a><span class="lineno"> 210</span> </div>
+<div class="line"><a name="l00211"></a><span class="lineno"> 211</span>  <span class="keyword">private</span>:</div>
+<div class="line"><a name="l00212"></a><span class="lineno"> 212</span> </div>
+<div class="line"><a name="l00213"></a><span class="lineno"> 213</span>  <span class="keywordtype">void</span> setOptions(TH1D *h, <span class="keywordtype">bool</span> dataLike = <span class="keyword">true</span>, <span class="keywordtype">bool</span> isBottom = <span class="keyword">false</span>);</div>
+<div class="line"><a name="l00214"></a><span class="lineno"> 214</span>  <span class="keywordtype">void</span> project(); <span class="comment">// push histo into the TH1D</span></div>
<div class="line"><a name="l00215"></a><span class="lineno"> 215</span> </div>
-<div class="line"><a name="l00216"></a><span class="lineno"> 216</span>  ClassDef(<a class="code" href="classhistoCompare.html">histoCompare</a>,1)</div>
-<div class="line"><a name="l00217"></a><span class="lineno"> 217</span>  </div>
-<div class="line"><a name="l00218"></a><span class="lineno"> 218</span> </div>
+<div class="line"><a name="l00216"></a><span class="lineno"> 216</span> </div>
+<div class="line"><a name="l00217"></a><span class="lineno"> 217</span>  ClassDef(<a class="code" href="classhistoCompare.html">histoCompare</a>,1)</div>
+<div class="line"><a name="l00218"></a><span class="lineno"> 218</span>  </div>
<div class="line"><a name="l00219"></a><span class="lineno"> 219</span> </div>
-<div class="line"><a name="l00220"></a><span class="lineno"> 220</span> };</div>
-<div class="line"><a name="l00221"></a><span class="lineno"> 221</span> </div>
+<div class="line"><a name="l00220"></a><span class="lineno"> 220</span> </div>
+<div class="line"><a name="l00221"></a><span class="lineno"> 221</span> };</div>
<div class="line"><a name="l00222"></a><span class="lineno"> 222</span> </div>
<div class="line"><a name="l00223"></a><span class="lineno"> 223</span> </div>
<div class="line"><a name="l00224"></a><span class="lineno"> 224</span> </div>
<div class="line"><a name="l00225"></a><span class="lineno"> 225</span> </div>
<div class="line"><a name="l00226"></a><span class="lineno"> 226</span> </div>
-<div class="line"><a name="l00227"></a><span class="lineno"> 227</span> <span class="preprocessor">#endif</span></div>
+<div class="line"><a name="l00227"></a><span class="lineno"> 227</span> </div>
+<div class="line"><a name="l00228"></a><span class="lineno"> 228</span> <span class="preprocessor">#endif</span></div>
<div class="ttc" id="classLKParameter_html"><div class="ttname"><a href="classLKParameter.html">LKParameter</a></div><div class="ttdef"><b>Definition:</b> XeStat.h:2991</div></div>
-<div class="ttc" id="classshapeSys_html"><div class="ttname"><a href="classshapeSys.html">shapeSys</a></div><div class="ttdef"><b>Definition:</b> XePdfObjects.h:24</div></div>
+<div class="ttc" id="classshapeSys_html"><div class="ttname"><a href="classshapeSys.html">shapeSys</a></div><div class="ttdoc">it represents a Shape Systematic. It inherits from LKParameter. Here you can use the LKParameter func...</div><div class="ttdef"><b>Definition:</b> XePdfObjects.h:25</div></div>
<div class="ttc" id="classerrorHandler_html"><div class="ttname"><a href="classerrorHandler.html">errorHandler</a></div><div class="ttdef"><b>Definition:</b> XeCore.h:2142</div></div>
-<div class="ttc" id="classscaleSys_html"><div class="ttname"><a href="classscaleSys.html">scaleSys</a></div><div class="ttdef"><b>Definition:</b> XePdfObjects.h:42</div></div>
-<div class="ttc" id="classhistoCompare_html"><div class="ttname"><a href="classhistoCompare.html">histoCompare</a></div><div class="ttdef"><b>Definition:</b> XePdfObjects.h:160</div></div>
-<div class="ttc" id="classpdfComponent_html"><div class="ttname"><a href="classpdfComponent.html">pdfComponent</a></div><div class="ttdef"><b>Definition:</b> XePdfObjects.h:67</div></div>
+<div class="ttc" id="classscaleSys_html"><div class="ttname"><a href="classscaleSys.html">scaleSys</a></div><div class="ttdef"><b>Definition:</b> XePdfObjects.h:43</div></div>
+<div class="ttc" id="classhistoCompare_html"><div class="ttname"><a href="classhistoCompare.html">histoCompare</a></div><div class="ttdef"><b>Definition:</b> XePdfObjects.h:161</div></div>
+<div class="ttc" id="classpdfComponent_html"><div class="ttname"><a href="classpdfComponent.html">pdfComponent</a></div><div class="ttdef"><b>Definition:</b> XePdfObjects.h:68</div></div>
</div><!-- fragment --></div><!-- contents -->
<!-- start footer part -->
<hr class="footer"/><address class="footer"><small>
diff --git a/docs/annotated.html b/docs/annotated.html
index 1146b96..13fb13e 100644
--- a/docs/annotated.html
+++ b/docs/annotated.html
@@ -166,7 +166,7 @@
<tr id="row_75_"><td class="entry"><img src="ftv2node.png" alt="o" width="16" height="22" /><img src="ftv2cl.png" alt="C" width="24" height="22" /><a class="el" href="classSelectionCutS1.html" target="_self">SelectionCutS1</a></td><td class="desc"></td></tr>
<tr id="row_76_" class="even"><td class="entry"><img src="ftv2node.png" alt="o" width="16" height="22" /><img src="ftv2cl.png" alt="C" width="24" height="22" /><a class="el" href="classSelectionCutUnsmearedS1.html" target="_self">SelectionCutUnsmearedS1</a></td><td class="desc"></td></tr>
<tr id="row_77_"><td class="entry"><img src="ftv2node.png" alt="o" width="16" height="22" /><img src="ftv2cl.png" alt="C" width="24" height="22" /><a class="el" href="classSensitivityBands.html" target="_self">SensitivityBands</a></td><td class="desc"></td></tr>
-<tr id="row_78_" class="even"><td class="entry"><img src="ftv2node.png" alt="o" width="16" height="22" /><img src="ftv2cl.png" alt="C" width="24" height="22" /><a class="el" href="classshapeSys.html" target="_self">shapeSys</a></td><td class="desc"></td></tr>
+<tr id="row_78_" class="even"><td class="entry"><img src="ftv2node.png" alt="o" width="16" height="22" /><img src="ftv2cl.png" alt="C" width="24" height="22" /><a class="el" href="classshapeSys.html" target="_self">shapeSys</a></td><td class="desc">It represents a Shape Systematic. It inherits from <a class="el" href="classLKParameter.html">LKParameter</a>. Here you can use the <a class="el" href="classLKParameter.html">LKParameter</a> functions to set initial val, min, max and setp </td></tr>
<tr id="row_79_"><td class="entry"><img src="ftv2node.png" alt="o" width="16" height="22" /><img src="ftv2cl.png" alt="C" width="24" height="22" /><a class="el" href="classSigmaParameter.html" target="_self">SigmaParameter</a></td><td class="desc"></td></tr>
<tr id="row_80_" class="even"><td class="entry"><img src="ftv2node.png" alt="o" width="16" height="22" /><img src="ftv2cl.png" alt="C" width="24" height="22" /><a class="el" href="classSigmaRange.html" target="_self">SigmaRange</a></td><td class="desc"></td></tr>
<tr id="row_81_"><td class="entry"><img src="ftv2node.png" alt="o" width="16" height="22" /><img src="ftv2cl.png" alt="C" width="24" height="22" /><a class="el" href="classSignalAndBackground.html" target="_self">SignalAndBackground</a></td><td class="desc"></td></tr>
diff --git a/docs/classLKParameter.html b/docs/classLKParameter.html
index be94f1a..8bba89d 100644
--- a/docs/classLKParameter.html
+++ b/docs/classLKParameter.html
@@ -105,7 +105,7 @@
<area href="classXeCore.html" alt="XeCore" shape="rect" coords="0,0,134,24"/>
<area href="classXeCore.html" alt="XeCore" shape="rect" coords="144,0,278,24"/>
<area href="classscaleSys.html" alt="scaleSys" shape="rect" coords="252,224,386,248"/>
-<area href="classshapeSys.html" alt="shapeSys" shape="rect" coords="252,280,386,304"/>
+<area href="classshapeSys.html" title="it represents a Shape Systematic. It inherits from LKParameter. Here you can use the LKParameter func..." alt="shapeSys" shape="rect" coords="252,280,386,304"/>
<area href="classSigmaParameter.html" alt="SigmaParameter" shape="rect" coords="252,336,386,360"/>
<area href="classTEfficiencyParameter.html" alt="TEfficiencyParameter" shape="rect" coords="252,392,386,416"/>
<area href="classTGaussParameter.html" alt="TGaussParameter" shape="rect" coords="252,448,386,472"/>
diff --git a/docs/classshapeSys.html b/docs/classshapeSys.html
index 5dc83d4..4694a32 100644
--- a/docs/classshapeSys.html
+++ b/docs/classshapeSys.html
@@ -90,6 +90,9 @@
</div><!--header-->
<div class="contents">
+<p>it represents a Shape Systematic. It inherits from <a class="el" href="classLKParameter.html">LKParameter</a>. Here you can use the <a class="el" href="classLKParameter.html">LKParameter</a> functions to set initial val, min, max and setp
+ <a href="classshapeSys.html#details">More...</a></p>
+
<p><code>#include <<a class="el" href="XePdfObjects_8h_source.html">XePdfObjects.h</a>></code></p>
<div class="dynheader">
Inheritance diagram for shapeSys:</div>
@@ -743,7 +746,7 @@
<tr class="separator:a77dd420f3efa9017becd6970480adfd7 inherit pro_static_attribs_classXeObject"><td class="memSeparator" colspan="2"> </td></tr>
</table>
<a name="details" id="details"></a><h2 class="groupheader">Detailed Description</h2>
-<div class="textblock"><p>Class representing a Shape Systematic, it inherits from <a class="el" href="classLKParameter.html">LKParameter</a>. here you can use the <a class="el" href="classLKParameter.html">LKParameter</a> functions to set initial val, min, max and setp </p>
+<div class="textblock"><p>it represents a Shape Systematic. It inherits from <a class="el" href="classLKParameter.html">LKParameter</a>. Here you can use the <a class="el" href="classLKParameter.html">LKParameter</a> functions to set initial val, min, max and setp </p>
</div><h2 class="groupheader">Constructor & Destructor Documentation</h2>
<a class="anchor" id="a1816a4031694074a2734f73a57eb93c1"></a>
<div class="memitem">
diff --git a/docs/hierarchy.html b/docs/hierarchy.html
index 633c57d..9ae2c01 100644
--- a/docs/hierarchy.html
+++ b/docs/hierarchy.html
@@ -197,7 +197,7 @@
<tr id="row_3_0_9_3_1_4_" style="display:none;"><td class="entry"><img src="ftv2vertline.png" alt="|" width="16" height="22" /><img src="ftv2blank.png" alt=" " width="16" height="22" /><img src="ftv2vertline.png" alt="|" width="16" height="22" /><img src="ftv2vertline.png" alt="|" width="16" height="22" /><img src="ftv2blank.png" alt=" " width="16" height="22" /><img src="ftv2lastnode.png" alt="\" width="16" height="22" /><img src="ftv2cl.png" alt="C" width="24" height="22" /><a class="el" href="classS1S2PL.html" target="_self">S1S2PL</a></td><td class="desc"></td></tr>
<tr id="row_3_0_9_4_" style="display:none;"><td class="entry"><img src="ftv2vertline.png" alt="|" width="16" height="22" /><img src="ftv2blank.png" alt=" " width="16" height="22" /><img src="ftv2vertline.png" alt="|" width="16" height="22" /><img id="arr_3_0_9_4_" src="ftv2pnode.png" alt="o" width="16" height="22" onclick="toggleFolder('3_0_9_4_')"/><img src="ftv2cl.png" alt="C" width="24" height="22" /><a class="el" href="classLKParameter.html" target="_self">LKParameter</a></td><td class="desc"></td></tr>
<tr id="row_3_0_9_4_0_" style="display:none;"><td class="entry"><img src="ftv2vertline.png" alt="|" width="16" height="22" /><img src="ftv2blank.png" alt=" " width="16" height="22" /><img src="ftv2vertline.png" alt="|" width="16" height="22" /><img src="ftv2vertline.png" alt="|" width="16" height="22" /><img src="ftv2node.png" alt="o" width="16" height="22" /><img src="ftv2cl.png" alt="C" width="24" height="22" /><a class="el" href="classscaleSys.html" target="_self">scaleSys</a></td><td class="desc"></td></tr>
-<tr id="row_3_0_9_4_1_" style="display:none;"><td class="entry"><img src="ftv2vertline.png" alt="|" width="16" height="22" /><img src="ftv2blank.png" alt=" " width="16" height="22" /><img src="ftv2vertline.png" alt="|" width="16" height="22" /><img src="ftv2vertline.png" alt="|" width="16" height="22" /><img src="ftv2node.png" alt="o" width="16" height="22" /><img src="ftv2cl.png" alt="C" width="24" height="22" /><a class="el" href="classshapeSys.html" target="_self">shapeSys</a></td><td class="desc"></td></tr>
+<tr id="row_3_0_9_4_1_" style="display:none;"><td class="entry"><img src="ftv2vertline.png" alt="|" width="16" height="22" /><img src="ftv2blank.png" alt=" " width="16" height="22" /><img src="ftv2vertline.png" alt="|" width="16" height="22" /><img src="ftv2vertline.png" alt="|" width="16" height="22" /><img src="ftv2node.png" alt="o" width="16" height="22" /><img src="ftv2cl.png" alt="C" width="24" height="22" /><a class="el" href="classshapeSys.html" target="_self">shapeSys</a></td><td class="desc">It represents a Shape Systematic. It inherits from <a class="el" href="classLKParameter.html">LKParameter</a>. Here you can use the <a class="el" href="classLKParameter.html">LKParameter</a> functions to set initial val, min, max and setp </td></tr>
<tr id="row_3_0_9_4_2_" style="display:none;"><td class="entry"><img src="ftv2vertline.png" alt="|" width="16" height="22" /><img src="ftv2blank.png" alt=" " width="16" height="22" /><img src="ftv2vertline.png" alt="|" width="16" height="22" /><img src="ftv2vertline.png" alt="|" width="16" height="22" /><img src="ftv2node.png" alt="o" width="16" height="22" /><img src="ftv2cl.png" alt="C" width="24" height="22" /><a class="el" href="classSigmaParameter.html" target="_self">SigmaParameter</a></td><td class="desc"></td></tr>
<tr id="row_3_0_9_4_3_" style="display:none;"><td class="entry"><img src="ftv2vertline.png" alt="|" width="16" height="22" /><img src="ftv2blank.png" alt=" " width="16" height="22" /><img src="ftv2vertline.png" alt="|" width="16" height="22" /><img src="ftv2vertline.png" alt="|" width="16" height="22" /><img src="ftv2node.png" alt="o" width="16" height="22" /><img src="ftv2cl.png" alt="C" width="24" height="22" /><a class="el" href="classTEfficiencyParameter.html" target="_self">TEfficiencyParameter</a></td><td class="desc"></td></tr>
<tr id="row_3_0_9_4_4_" style="display:none;"><td class="entry"><img src="ftv2vertline.png" alt="|" width="16" height="22" /><img src="ftv2blank.png" alt=" " width="16" height="22" /><img src="ftv2vertline.png" alt="|" width="16" height="22" /><img src="ftv2vertline.png" alt="|" width="16" height="22" /><img src="ftv2node.png" alt="o" width="16" height="22" /><img src="ftv2cl.png" alt="C" width="24" height="22" /><a class="el" href="classTGaussParameter.html" target="_self">TGaussParameter</a></td><td class="desc"></td></tr>
diff --git a/src/XePdfObjects.h b/src/XePdfObjects.h
index f33959e..d5f4d8f 100644
--- a/src/XePdfObjects.h
+++ b/src/XePdfObjects.h
@@ -17,8 +17,9 @@ using namespace std;
/**
- * Class representing a Shape Systematic, it inherits from LKParameter.
- * here you can use the LKParameter functions to set initial val,
+ * \class shapeSys
+ * \brief it represents a Shape Systematic.
+ * It inherits from LKParameter. Here you can use the LKParameter functions to set initial val,
* min, max and setp
*/
class shapeSys : public LKParameter {