MiniTreeDataSet_8hpp_source.html 20 KB

1234567891011121314151617181920212223242526272829303132333435363738394041424344454647484950515253545556575859606162636465666768697071727374757677787980818283848586
  1. <!DOCTYPE html PUBLIC "-//W3C//DTD XHTML 1.0 Transitional//EN" "http://www.w3.org/TR/xhtml1/DTD/xhtml1-transitional.dtd">
  2. <html xmlns="http://www.w3.org/1999/xhtml">
  3. <head>
  4. <meta http-equiv="Content-Type" content="text/xhtml;charset=UTF-8"/>
  5. <meta http-equiv="X-UA-Compatible" content="IE=9"/>
  6. <meta name="generator" content="Doxygen 1.8.13"/>
  7. <meta name="viewport" content="width=device-width, initial-scale=1"/>
  8. <title>TTTT Analysis: /home/caleb/Sources/TTTT/analysis/MiniTreeDataSet.hpp Source File</title>
  9. <link href="tabs.css" rel="stylesheet" type="text/css"/>
  10. <script type="text/javascript" src="jquery.js"></script>
  11. <script type="text/javascript" src="dynsections.js"></script>
  12. <link href="search/search.css" rel="stylesheet" type="text/css"/>
  13. <script type="text/javascript" src="search/searchdata.js"></script>
  14. <script type="text/javascript" src="search/search.js"></script>
  15. <link href="doxygen.css" rel="stylesheet" type="text/css" />
  16. </head>
  17. <body>
  18. <div id="top"><!-- do not remove this div, it is closed by doxygen! -->
  19. <div id="titlearea">
  20. <table cellspacing="0" cellpadding="0">
  21. <tbody>
  22. <tr style="height: 56px;">
  23. <td id="projectalign" style="padding-left: 0.5em;">
  24. <div id="projectname">TTTT Analysis
  25. &#160;<span id="projectnumber">0.1</span>
  26. </div>
  27. </td>
  28. </tr>
  29. </tbody>
  30. </table>
  31. </div>
  32. <!-- end header part -->
  33. <!-- Generated by Doxygen 1.8.13 -->
  34. <script type="text/javascript">
  35. var searchBox = new SearchBox("searchBox", "search",false,'Search');
  36. </script>
  37. <script type="text/javascript" src="menudata.js"></script>
  38. <script type="text/javascript" src="menu.js"></script>
  39. <script type="text/javascript">
  40. $(function() {
  41. initMenu('',true,false,'search.php','Search');
  42. $(document).ready(function() { init_search(); });
  43. });
  44. </script>
  45. <div id="main-nav"></div>
  46. <!-- window showing the filter options -->
  47. <div id="MSearchSelectWindow"
  48. onmouseover="return searchBox.OnSearchSelectShow()"
  49. onmouseout="return searchBox.OnSearchSelectHide()"
  50. onkeydown="return searchBox.OnSearchSelectKey(event)">
  51. </div>
  52. <!-- iframe showing the search results (closed by default) -->
  53. <div id="MSearchResultsWindow">
  54. <iframe src="javascript:void(0)" frameborder="0"
  55. name="MSearchResults" id="MSearchResults">
  56. </iframe>
  57. </div>
  58. <div id="nav-path" class="navpath">
  59. <ul>
  60. <li class="navelem"><a class="el" href="dir_b8678fa8510b7ff9a55ffd4d18fd5e47.html">analysis</a></li> </ul>
  61. </div>
  62. </div><!-- top -->
  63. <div class="header">
  64. <div class="headertitle">
  65. <div class="title">MiniTreeDataSet.hpp</div> </div>
  66. </div><!--header-->
  67. <div class="contents">
  68. <a href="MiniTreeDataSet_8hpp.html">Go to the documentation of this file.</a><div class="fragment"><div class="line"><a name="l00001"></a><span class="lineno"> 1</span>&#160;</div><div class="line"><a name="l00031"></a><span class="lineno"> 31</span>&#160;<span class="preprocessor">#ifndef minitreedataset_h</span></div><div class="line"><a name="l00032"></a><span class="lineno"> 32</span>&#160;<span class="preprocessor">#define minitreedataset_h</span></div><div class="line"><a name="l00033"></a><span class="lineno"> 33</span>&#160;</div><div class="line"><a name="l00034"></a><span class="lineno"> 34</span>&#160;<span class="preprocessor">#include &lt;string&gt;</span></div><div class="line"><a name="l00035"></a><span class="lineno"> 35</span>&#160;<span class="preprocessor">#include &lt;tuple&gt;</span></div><div class="line"><a name="l00036"></a><span class="lineno"> 36</span>&#160;</div><div class="line"><a name="l00037"></a><span class="lineno"> 37</span>&#160;<span class="preprocessor">#include &quot;filval/filval.hpp&quot;</span></div><div class="line"><a name="l00038"></a><span class="lineno"> 38</span>&#160;<span class="preprocessor">#include &quot;filval_root/filval_root.hpp&quot;</span></div><div class="line"><a name="l00039"></a><span class="lineno"> 39</span>&#160;<span class="preprocessor">#include &quot;<a class="code" href="MiniTree_8hpp.html">MiniTree.hpp</a>&quot;</span></div><div class="line"><a name="l00040"></a><span class="lineno"> 40</span>&#160;</div><div class="line"><a name="l00041"></a><span class="lineno"> 41</span>&#160;<span class="keyword">using namespace </span><a class="code" href="namespacestd.html">std</a>;</div><div class="line"><a name="l00042"></a><span class="lineno"> 42</span>&#160;<span class="keyword">using namespace </span><a class="code" href="namespacefv.html">fv</a>;</div><div class="line"><a name="l00043"></a><span class="lineno"> 43</span>&#160;<span class="keyword">using namespace </span><a class="code" href="namespacefv_1_1root.html">fv::root</a>;</div><div class="line"><a name="l00044"></a><span class="lineno"> 44</span>&#160;</div><div class="line"><a name="l00045"></a><span class="lineno"> 45</span>&#160;<span class="keyword">class </span>MiniTreeDataSet : <span class="keyword">public</span> DataSet,</div><div class="line"><a name="l00046"></a><span class="lineno"> 46</span>&#160; <span class="keyword">public</span> MiniTree{</div><div class="line"><a name="l00047"></a><span class="lineno"> 47</span>&#160; <span class="keyword">private</span>:</div><div class="line"><a name="l00048"></a><span class="lineno"> 48</span>&#160; std::string input_filename;</div><div class="line"><a name="l00049"></a><span class="lineno"> 49</span>&#160; std::string output_filename;</div><div class="line"><a name="l00050"></a><span class="lineno"> 50</span>&#160; TFile* input_file;</div><div class="line"><a name="l00051"></a><span class="lineno"> 51</span>&#160; TFile* output_file;</div><div class="line"><a name="l00052"></a><span class="lineno"> 52</span>&#160; <span class="keywordtype">long</span> next_entry;</div><div class="line"><a name="l00053"></a><span class="lineno"> 53</span>&#160; <span class="keywordtype">long</span> nentries;</div><div class="line"><a name="l00054"></a><span class="lineno"> 54</span>&#160; <span class="keywordtype">bool</span> load_next(){</div><div class="line"><a name="l00055"></a><span class="lineno"> 55</span>&#160; <span class="keywordflow">if</span> (next_entry &gt;= nentries) <span class="keywordflow">return</span> <span class="keyword">false</span>;</div><div class="line"><a name="l00056"></a><span class="lineno"> 56</span>&#160; fChain-&gt;GetEntry(next_entry);</div><div class="line"><a name="l00057"></a><span class="lineno"> 57</span>&#160; ++next_entry;</div><div class="line"><a name="l00058"></a><span class="lineno"> 58</span>&#160; <span class="keywordflow">return</span> <span class="keyword">true</span>;</div><div class="line"><a name="l00059"></a><span class="lineno"> 59</span>&#160; }</div><div class="line"><a name="l00060"></a><span class="lineno"> 60</span>&#160; <span class="keywordtype">int</span> get_events(){</div><div class="line"><a name="l00061"></a><span class="lineno"> 61</span>&#160; <span class="keywordflow">return</span> nentries;</div><div class="line"><a name="l00062"></a><span class="lineno"> 62</span>&#160; }</div><div class="line"><a name="l00063"></a><span class="lineno"> 63</span>&#160; <span class="keywordtype">int</span> get_current_event(){</div><div class="line"><a name="l00064"></a><span class="lineno"> 64</span>&#160; <span class="keywordflow">return</span> next_entry-1;</div><div class="line"><a name="l00065"></a><span class="lineno"> 65</span>&#160; }</div><div class="line"><a name="l00066"></a><span class="lineno"> 66</span>&#160;</div><div class="line"><a name="l00067"></a><span class="lineno"> 67</span>&#160; <span class="keyword">public</span>:</div><div class="line"><a name="l00068"></a><span class="lineno"> 68</span>&#160; MiniTreeDataSet(<span class="keyword">const</span> std::string&amp; input_filename, <span class="keyword">const</span> std::string output_filename)</div><div class="line"><a name="l00069"></a><span class="lineno"> 69</span>&#160; :input_filename(input_filename),</div><div class="line"><a name="l00070"></a><span class="lineno"> 70</span>&#160; output_filename(output_filename),</div><div class="line"><a name="l00071"></a><span class="lineno"> 71</span>&#160; next_entry(0) {</div><div class="line"><a name="l00072"></a><span class="lineno"> 72</span>&#160; input_file = TFile::Open(input_filename.c_str());</div><div class="line"><a name="l00073"></a><span class="lineno"> 73</span>&#160; Init((TTree*) input_file-&gt;Get(<span class="stringliteral">&quot;tree&quot;</span>));</div><div class="line"><a name="l00074"></a><span class="lineno"> 74</span>&#160; nentries = fChain-&gt;GetEntriesFast();</div><div class="line"><a name="l00075"></a><span class="lineno"> 75</span>&#160; output_file = TFile::Open(output_filename.c_str(), <span class="stringliteral">&quot;RECREATE&quot;</span>);</div><div class="line"><a name="l00076"></a><span class="lineno"> 76</span>&#160; this-&gt;fChain-&gt;SetBranchStatus(<span class="stringliteral">&quot;*&quot;</span>, <span class="keyword">false</span>);</div><div class="line"><a name="l00077"></a><span class="lineno"> 77</span>&#160; }</div><div class="line"><a name="l00078"></a><span class="lineno"> 78</span>&#160;</div><div class="line"><a name="l00079"></a><span class="lineno"> 79</span>&#160; ~MiniTreeDataSet(){</div><div class="line"><a name="l00080"></a><span class="lineno"> 80</span>&#160; input_file-&gt;Close();</div><div class="line"><a name="l00081"></a><span class="lineno"> 81</span>&#160; output_file-&gt;Close();</div><div class="line"><a name="l00082"></a><span class="lineno"> 82</span>&#160; }</div><div class="line"><a name="l00083"></a><span class="lineno"> 83</span>&#160;</div><div class="line"><a name="l00084"></a><span class="lineno"> 84</span>&#160; <span class="keyword">template</span> &lt;<span class="keyword">typename</span> T&gt;</div><div class="line"><a name="l00085"></a><span class="lineno"> 85</span>&#160; <a class="code" href="classfv_1_1Value.html">Value&lt;T&gt;</a>* track_branch(<span class="keyword">const</span> std::string&amp; bname){</div><div class="line"><a name="l00086"></a><span class="lineno"> 86</span>&#160; TBranch* branch = fChain-&gt;GetBranch(bname.c_str());</div><div class="line"><a name="l00087"></a><span class="lineno"> 87</span>&#160; <span class="keywordflow">if</span> (branch == <span class="keyword">nullptr</span>){</div><div class="line"><a name="l00088"></a><span class="lineno"> 88</span>&#160; CRITICAL(<span class="stringliteral">&quot;Branch: &quot;</span> &lt;&lt; bname &lt;&lt; <span class="stringliteral">&quot; does not exist in input tree.&quot;</span>, -1);</div><div class="line"><a name="l00089"></a><span class="lineno"> 89</span>&#160; }</div><div class="line"><a name="l00090"></a><span class="lineno"> 90</span>&#160; T* bref = (T*) branch-&gt;GetAddress();</div><div class="line"><a name="l00091"></a><span class="lineno"> 91</span>&#160; fChain-&gt;SetBranchStatus(bname.c_str(), <span class="keyword">true</span>);</div><div class="line"><a name="l00092"></a><span class="lineno"> 92</span>&#160; INFO(<span class="stringliteral">&quot;Registering branch \&quot;&quot;</span> &lt;&lt; bname</div><div class="line"><a name="l00093"></a><span class="lineno"> 93</span>&#160; &lt;&lt; <span class="stringliteral">&quot;\&quot; with address &quot;</span> &lt;&lt; bref</div><div class="line"><a name="l00094"></a><span class="lineno"> 94</span>&#160; &lt;&lt; <span class="stringliteral">&quot; and type &quot;</span> &lt;&lt; <span class="keyword">typeid</span>(bref).name());</div><div class="line"><a name="l00095"></a><span class="lineno"> 95</span>&#160; <span class="keywordflow">return</span> <span class="keyword">new</span> <a class="code" href="classfv_1_1ObservedValue.html">ObservedValue&lt;T&gt;</a>(bname, bref);</div><div class="line"><a name="l00096"></a><span class="lineno"> 96</span>&#160; }</div><div class="line"><a name="l00097"></a><span class="lineno"> 97</span>&#160;</div><div class="line"><a name="l00098"></a><span class="lineno"> 98</span>&#160; <span class="keyword">template</span> &lt;<span class="keyword">typename</span> T&gt;</div><div class="line"><a name="l00099"></a><span class="lineno"> 99</span>&#160; <a class="code" href="classfv_1_1Value.html">Value&lt;T*&gt;</a>* track_branch_ptr(<span class="keyword">const</span> std::string&amp; bname){</div><div class="line"><a name="l00100"></a><span class="lineno"> 100</span>&#160; TBranch* branch = fChain-&gt;GetBranch(bname.c_str());</div><div class="line"><a name="l00101"></a><span class="lineno"> 101</span>&#160; <span class="keywordflow">if</span> (branch == <span class="keyword">nullptr</span>){</div><div class="line"><a name="l00102"></a><span class="lineno"> 102</span>&#160; CRITICAL(<span class="stringliteral">&quot;Branch: &quot;</span> &lt;&lt; bname &lt;&lt; <span class="stringliteral">&quot; does not exist in input tree.&quot;</span>, -1);</div><div class="line"><a name="l00103"></a><span class="lineno"> 103</span>&#160; }</div><div class="line"><a name="l00104"></a><span class="lineno"> 104</span>&#160; T* bref = (T*) branch-&gt;GetAddress();</div><div class="line"><a name="l00105"></a><span class="lineno"> 105</span>&#160; fChain-&gt;SetBranchStatus(bname.c_str(), <span class="keyword">true</span>);</div><div class="line"><a name="l00106"></a><span class="lineno"> 106</span>&#160; INFO(<span class="stringliteral">&quot;Registering pointer branch \&quot;&quot;</span> &lt;&lt; bname</div><div class="line"><a name="l00107"></a><span class="lineno"> 107</span>&#160; &lt;&lt; <span class="stringliteral">&quot;\&quot; with address &quot;</span> &lt;&lt; bref</div><div class="line"><a name="l00108"></a><span class="lineno"> 108</span>&#160; &lt;&lt; <span class="stringliteral">&quot; and type &quot;</span> &lt;&lt; <span class="keyword">typeid</span>(bref).name());</div><div class="line"><a name="l00109"></a><span class="lineno"> 109</span>&#160; <span class="keywordflow">return</span> <span class="keyword">new</span> <a class="code" href="classfv_1_1PointerValue.html">PointerValue&lt;T&gt;</a>(bname, bref);</div><div class="line"><a name="l00110"></a><span class="lineno"> 110</span>&#160; }</div><div class="line"><a name="l00111"></a><span class="lineno"> 111</span>&#160;</div><div class="line"><a name="l00112"></a><span class="lineno"> 112</span>&#160; <span class="keyword">template</span> &lt;<span class="keyword">typename</span> T&gt;</div><div class="line"><a name="l00113"></a><span class="lineno"> 113</span>&#160; decltype(<span class="keyword">auto</span>) track_branch_vec(const <a class="code" href="namespacestd.html">std</a>::<span class="keywordtype">string</span>&amp; size_bname, const <a class="code" href="namespacestd.html">std</a>::<span class="keywordtype">string</span>&amp; value_bname){</div><div class="line"><a name="l00114"></a><span class="lineno"> 114</span>&#160; track_branch_ptr&lt;T&gt;(value_bname);</div><div class="line"><a name="l00115"></a><span class="lineno"> 115</span>&#160; <span class="keywordflow">return</span> wrapper_vector&lt;T&gt;(lookup&lt;int&gt;(size_bname), lookup&lt;T*&gt;(value_bname), value_bname);</div><div class="line"><a name="l00116"></a><span class="lineno"> 116</span>&#160; }</div><div class="line"><a name="l00117"></a><span class="lineno"> 117</span>&#160;</div><div class="line"><a name="l00118"></a><span class="lineno"> 118</span>&#160; <span class="keywordtype">void</span> save_all(){</div><div class="line"><a name="l00119"></a><span class="lineno"> 119</span>&#160; output_file-&gt;cd();</div><div class="line"><a name="l00120"></a><span class="lineno"> 120</span>&#160; <span class="keywordflow">for</span>(<span class="keyword">auto</span> container : containers){</div><div class="line"><a name="l00121"></a><span class="lineno"> 121</span>&#160; container.second-&gt;save_as(<span class="stringliteral">&quot;outfile&quot;</span>, SaveOption::ROOT);</div><div class="line"><a name="l00122"></a><span class="lineno"> 122</span>&#160; }</div><div class="line"><a name="l00123"></a><span class="lineno"> 123</span>&#160; <span class="comment">// Save the value names for each container to enable looking up</span></div><div class="line"><a name="l00124"></a><span class="lineno"> 124</span>&#160; <span class="comment">// what values are plotted</span></div><div class="line"><a name="l00125"></a><span class="lineno"> 125</span>&#160; std::map&lt;string,string&gt; value_lookup = this-&gt;get_container_name_value_map();</div><div class="line"><a name="l00126"></a><span class="lineno"> 126</span>&#160; gDirectory-&gt;WriteObjectAny(&amp;value_lookup, <span class="stringliteral">&quot;std::map&lt;std::string,std::string&gt;&quot;</span>, <span class="stringliteral">&quot;_value_lookup&quot;</span>);</div><div class="line"><a name="l00127"></a><span class="lineno"> 127</span>&#160;</div><div class="line"><a name="l00128"></a><span class="lineno"> 128</span>&#160; std::map&lt;string,string&gt; fn_impl_lookup = this-&gt;get_function_name_impl_map();</div><div class="line"><a name="l00129"></a><span class="lineno"> 129</span>&#160; gDirectory-&gt;WriteObjectAny(&amp;fn_impl_lookup, <span class="stringliteral">&quot;std::map&lt;std::string,std::string&gt;&quot;</span>, <span class="stringliteral">&quot;_function_impl_lookup&quot;</span>);</div><div class="line"><a name="l00130"></a><span class="lineno"> 130</span>&#160; }</div><div class="line"><a name="l00131"></a><span class="lineno"> 131</span>&#160;};</div><div class="line"><a name="l00132"></a><span class="lineno"> 132</span>&#160;<span class="preprocessor">#endif // minitreedataset_h</span></div><div class="ttc" id="namespacestd_html"><div class="ttname"><a href="namespacestd.html">std</a></div></div>
  69. <div class="ttc" id="classfv_1_1ObservedValue_html"><div class="ttname"><a href="classfv_1_1ObservedValue.html">fv::ObservedValue</a></div><div class="ttdoc">A value supplied by the dataset, not derived. </div><div class="ttdef"><b>Definition:</b> value.hpp:448</div></div>
  70. <div class="ttc" id="namespacefv_html"><div class="ttname"><a href="namespacefv.html">fv</a></div><div class="ttdoc">The namespace containing all filval classes and functions. </div><div class="ttdef"><b>Definition:</b> api.hpp:38</div></div>
  71. <div class="ttc" id="classfv_1_1Value_html"><div class="ttname"><a href="classfv_1_1Value.html">fv::Value</a></div><div class="ttdoc">A templated value. </div><div class="ttdef"><b>Definition:</b> value.hpp:415</div></div>
  72. <div class="ttc" id="namespacefv_1_1root_html"><div class="ttname"><a href="namespacefv_1_1root.html">fv::root</a></div><div class="ttdef"><b>Definition:</b> api.hpp:8</div></div>
  73. <div class="ttc" id="MiniTree_8hpp_html"><div class="ttname"><a href="MiniTree_8hpp.html">MiniTree.hpp</a></div></div>
  74. <div class="ttc" id="classfv_1_1PointerValue_html"><div class="ttname"><a href="classfv_1_1PointerValue.html">fv::PointerValue</a></div><div class="ttdoc">A Value of a pointer. </div><div class="ttdef"><b>Definition:</b> value.hpp:1208</div></div>
  75. </div><!-- fragment --></div><!-- contents -->
  76. <!-- start footer part -->
  77. <hr class="footer"/><address class="footer"><small>
  78. Generated by &#160;<a href="http://www.doxygen.org/index.html">
  79. <img class="footer" src="doxygen.png" alt="doxygen"/>
  80. </a> 1.8.13
  81. </small></address>
  82. </body>
  83. </html>