/[PAMELA software]/PamCut/doc/Doxygen/html/CaloAxis2_8h_source.html
ViewVC logotype

Annotation of /PamCut/doc/Doxygen/html/CaloAxis2_8h_source.html

Parent Directory Parent Directory | Revision Log Revision Log

Revision 1.1.1.1 - (hide annotations) (download) (as text) (vendor branch)
Wed May 27 13:30:07 2009 UTC (15 years, 10 months ago) by pam-fi
Branch: DEV
CVS Tags: v0r00
Changes since 1.1: +0 -0 lines
File MIME type: text/html 
First import
1 pam-fi 1.1 <!DOCTYPE HTML PUBLIC "-//W3C//DTD HTML 4.01 Transitional//EN">
2     <html><head><meta http-equiv="Content-Type" content="text/html;charset=UTF-8">
3     <title>PamCut: CaloAxis2.h Source File</title>
4     <link href="tabs.css" rel="stylesheet" type="text/css">
5     <link href="doxygen.css" rel="stylesheet" type="text/css">
6     </head><body>
7     <!-- Generated by Doxygen 1.5.9 -->
8     <div class="navigation" id="top">
9     <div class="tabs">
10     <ul>
11     <li><a href="index.html"><span>Main&nbsp;Page</span></a></li>
12     <li><a href="annotated.html"><span>Classes</span></a></li>
13     <li class="current"><a href="files.html"><span>Files</span></a></li>
14     </ul>
15     </div>
16     <div class="tabs">
17     <ul>
18     <li><a href="files.html"><span>File&nbsp;List</span></a></li>
19     <li><a href="globals.html"><span>File&nbsp;Members</span></a></li>
20     </ul>
21     </div>
22     <h1>CaloAxis2.h</h1><a href="CaloAxis2_8h.html">Go to the documentation of this file.</a><div class="fragment"><pre class="fragment"><a name="l00001"></a>00001 <span class="preprocessor">#ifndef CALOAXIS2_H_</span>
23     <a name="l00002"></a>00002 <span class="preprocessor"></span><span class="preprocessor">#define CALOAXIS2_H_</span>
24     <a name="l00003"></a>00003 <span class="preprocessor"></span>
25     <a name="l00004"></a>00004 <span class="preprocessor">#if !defined(__CINT__) || defined(__MAKECINT__)</span>
26     <a name="l00005"></a>00005 <span class="preprocessor"></span>
27     <a name="l00006"></a>00006
28     <a name="l00007"></a>00007 <span class="preprocessor">#include &lt;TString.h&gt;</span>
29     <a name="l00008"></a>00008 <span class="preprocessor">#include &lt;TH1F.h&gt;</span>
30     <a name="l00009"></a>00009 <span class="preprocessor">#include &lt;TH2F.h&gt;</span>
31     <a name="l00010"></a>00010 <span class="preprocessor">#include &lt;TMath.h&gt;</span>
32     <a name="l00011"></a>00011 <span class="preprocessor">#include &lt;TLine.h&gt;</span>
33     <a name="l00012"></a>00012 <span class="preprocessor">#include &lt;TPolyMarker.h&gt;</span>
34     <a name="l00013"></a>00013 <span class="preprocessor">#include &lt;TSelector.h&gt;</span>
35     <a name="l00014"></a>00014 <span class="preprocessor">#include &lt;TFile.h&gt;</span>
36     <a name="l00015"></a>00015
37     <a name="l00016"></a>00016 <span class="preprocessor">#include &lt;stdlib.h&gt;</span>
38     <a name="l00017"></a>00017 <span class="preprocessor">#include &lt;iostream&gt;</span>
39     <a name="l00018"></a>00018 <span class="keyword">using namespace </span>std;
40     <a name="l00019"></a>00019
41     <a name="l00020"></a>00020 <span class="preprocessor">#include &lt;PamLevel2.h&gt;</span>
42     <a name="l00021"></a>00021 <span class="comment">//#include &lt;CaloEvent.h&gt;</span>
43     <a name="l00022"></a>00022
44     <a name="l00023"></a>00023 <span class="preprocessor">#endif // !defined(__CINT__) || defined(__MAKECINT__)</span>
45     <a name="l00024"></a>00024 <span class="preprocessor"></span><span class="comment">//===============================================================================</span>
46     <a name="l00025"></a>00025 <span class="comment">//</span>
47     <a name="l00026"></a>00026 <span class="comment">//</span>
48     <a name="l00027"></a>00027 <span class="comment">//</span>
49     <a name="l00028"></a>00028 <span class="comment">//</span>
50     <a name="l00029"></a>00029 <span class="comment">//===============================================================================</span>
51     <a name="l00030"></a><a class="code" href="classCaloEvent.html">00030</a> <span class="keyword">class </span><a class="code" href="classCaloEvent.html">CaloEvent</a>{
52     <a name="l00031"></a>00031
53     <a name="l00032"></a>00032 <span class="keyword">public</span>:
54     <a name="l00033"></a>00033 <span class="comment">// calo event</span>
55     <a name="l00034"></a>00034 <span class="comment">// ------------------------------</span>
56     <a name="l00035"></a>00035 <span class="comment">// hit strips</span>
57     <a name="l00036"></a>00036 <span class="comment">// ------------------------------</span>
58     <a name="l00037"></a><a class="code" href="classCaloEvent.html#593e99c95482f5e63e43eaa8c843bba2">00037</a> vector&lt;int&gt; <a class="code" href="classCaloEvent.html#593e99c95482f5e63e43eaa8c843bba2">ids</a>; <span class="comment">//strip 0-95</span>
59     <a name="l00038"></a><a class="code" href="classCaloEvent.html#3ad199ffd33f6300bff64e854c6ad652">00038</a> vector&lt;int&gt; <a class="code" href="classCaloEvent.html#3ad199ffd33f6300bff64e854c6ad652">idp</a>; <span class="comment">//plane 0-21</span>
60     <a name="l00039"></a><a class="code" href="classCaloEvent.html#688d5d2803d4688d4f1ef72cbaf42dcd">00039</a> vector&lt;float&gt; <a class="code" href="classCaloEvent.html#688d5d2803d4688d4f1ef72cbaf42dcd">q</a>; <span class="comment">// signal</span>
61     <a name="l00040"></a><a class="code" href="classCaloEvent.html#340ff2ae2167b37d0b6acc56b1475daf">00040</a> vector&lt;float&gt; <a class="code" href="classCaloEvent.html#340ff2ae2167b37d0b6acc56b1475daf">zcoord</a>;
62     <a name="l00041"></a><a class="code" href="classCaloEvent.html#b48403ffdb4114d57e4ea4058438e814">00041</a> vector&lt;float&gt; <a class="code" href="classCaloEvent.html#b48403ffdb4114d57e4ea4058438e814">xycoord</a>;
63     <a name="l00042"></a>00042
64     <a name="l00043"></a>00043
65     <a name="l00044"></a>00044 <span class="comment">// ------------------------------</span>
66     <a name="l00045"></a>00045 <span class="comment">// methods</span>
67     <a name="l00046"></a>00046 <span class="comment">// ------------------------------</span>
68     <a name="l00047"></a><a class="code" href="classCaloEvent.html#fd2bc3739bbbf688c9323d8a96e17c5f">00047</a> <span class="keywordtype">void</span> <a class="code" href="classCaloEvent.html#fd2bc3739bbbf688c9323d8a96e17c5f">Add</a>(<span class="keywordtype">int</span> iis, <span class="keywordtype">int</span> iip, <span class="keywordtype">float</span> fq, <span class="keywordtype">float</span> fz, <span class="keywordtype">float</span> fxy){
69     <a name="l00048"></a>00048 <a class="code" href="classCaloEvent.html#593e99c95482f5e63e43eaa8c843bba2">ids</a>.push_back(iis);
70     <a name="l00049"></a>00049 <a class="code" href="classCaloEvent.html#3ad199ffd33f6300bff64e854c6ad652">idp</a>.push_back(iip);
71     <a name="l00050"></a>00050 <a class="code" href="classCaloEvent.html#688d5d2803d4688d4f1ef72cbaf42dcd">q</a>.push_back(fq);
72     <a name="l00051"></a>00051 <a class="code" href="classCaloEvent.html#340ff2ae2167b37d0b6acc56b1475daf">zcoord</a>.push_back(fz);
73     <a name="l00052"></a>00052 <a class="code" href="classCaloEvent.html#b48403ffdb4114d57e4ea4058438e814">xycoord</a>.push_back(fxy);
74     <a name="l00053"></a>00053 };
75     <a name="l00054"></a>00054
76     <a name="l00055"></a>00055
77     <a name="l00056"></a>00056
78     <a name="l00057"></a><a class="code" href="classCaloEvent.html#2b2597fe0263fe128101a6ab02520130">00057</a> <span class="keywordtype">int</span> <a class="code" href="classCaloEvent.html#2b2597fe0263fe128101a6ab02520130">GetN</a>(){<span class="keywordflow">return</span> (<span class="keywordtype">int</span>)(<a class="code" href="classCaloEvent.html#688d5d2803d4688d4f1ef72cbaf42dcd">q</a>.size());};
79     <a name="l00058"></a>00058
80     <a name="l00059"></a><a class="code" href="classCaloEvent.html#6cb8e9a295bf452dbca12f3089045706">00059</a> <span class="keywordtype">void</span> <a class="code" href="classCaloEvent.html#6cb8e9a295bf452dbca12f3089045706">Reset</a>(){
81     <a name="l00060"></a>00060 <a class="code" href="classCaloEvent.html#593e99c95482f5e63e43eaa8c843bba2">ids</a>.clear();
82     <a name="l00061"></a>00061 <a class="code" href="classCaloEvent.html#3ad199ffd33f6300bff64e854c6ad652">idp</a>.clear();
83     <a name="l00062"></a>00062 <a class="code" href="classCaloEvent.html#688d5d2803d4688d4f1ef72cbaf42dcd">q</a>.clear();
84     <a name="l00063"></a>00063 <a class="code" href="classCaloEvent.html#340ff2ae2167b37d0b6acc56b1475daf">zcoord</a>.clear();
85     <a name="l00064"></a>00064 <a class="code" href="classCaloEvent.html#b48403ffdb4114d57e4ea4058438e814">xycoord</a>.clear();
86     <a name="l00065"></a>00065 };
87     <a name="l00066"></a>00066
88     <a name="l00067"></a><a class="code" href="classCaloEvent.html#aa540318b81bbed181ab8cefa5df87ac">00067</a> <span class="keywordtype">void</span> <a class="code" href="classCaloEvent.html#aa540318b81bbed181ab8cefa5df87ac">Delete</a>(){ <a class="code" href="classCaloEvent.html#6cb8e9a295bf452dbca12f3089045706">Reset</a>(); };
89     <a name="l00068"></a>00068
90     <a name="l00069"></a><a class="code" href="classCaloEvent.html#d288c6c79268b65eae838c461e11b8eb">00069</a> <a class="code" href="classCaloEvent.html#d288c6c79268b65eae838c461e11b8eb">~CaloEvent</a>(){ <a class="code" href="classCaloEvent.html#aa540318b81bbed181ab8cefa5df87ac">Delete</a>(); };
91     <a name="l00070"></a>00070
92     <a name="l00071"></a><a class="code" href="classCaloEvent.html#67cd3b4a14bcdb5d15b40a177ca37fe3">00071</a> <a class="code" href="classCaloEvent.html#67cd3b4a14bcdb5d15b40a177ca37fe3">CaloEvent</a>(){ };
93     <a name="l00072"></a>00072
94     <a name="l00073"></a><a class="code" href="classCaloEvent.html#3eadf14b2b405c65a695fac497ec3c22">00073</a> <a class="code" href="classCaloEvent.html#3eadf14b2b405c65a695fac497ec3c22">CaloEvent</a>( CaloLevel1* calo, <span class="keywordtype">int</span> view, Bool_t usemechanicalalignement){ <a class="code" href="classCaloEvent.html#538af1924be400124d003802302e71dd">Set</a>(calo,view,usemechanicalalignement); };
95     <a name="l00074"></a>00074
96     <a name="l00075"></a><a class="code" href="classCaloEvent.html#d7beea6828722ab1ecff5e25fa33210a">00075</a> <a class="code" href="classCaloEvent.html#d7beea6828722ab1ecff5e25fa33210a">CaloEvent</a>( CaloLevel1* calo, <span class="keywordtype">int</span> view){ <a class="code" href="classCaloEvent.html#538af1924be400124d003802302e71dd">Set</a>(calo,view); };
97     <a name="l00076"></a>00076
98     <a name="l00077"></a><a class="code" href="classCaloEvent.html#538af1924be400124d003802302e71dd">00077</a> <span class="keywordtype">void</span> <a class="code" href="classCaloEvent.html#538af1924be400124d003802302e71dd">Set</a>(CaloLevel1* calo, <span class="keywordtype">int</span> view, Bool_t usemechanicalalignement){
99     <a name="l00078"></a>00078
100     <a name="l00079"></a>00079 <a class="code" href="classCaloEvent.html#6cb8e9a295bf452dbca12f3089045706">Reset</a>();
101     <a name="l00080"></a>00080 <span class="comment">// cout &lt;&lt; " CaloEvent::Set()"&lt;&lt;endl;</span>
102     <a name="l00081"></a>00081 CaloStrip cstrip;
103     <a name="l00082"></a>00082 cstrip = CaloStrip(calo,usemechanicalalignement);
104     <a name="l00083"></a>00083
105     <a name="l00084"></a>00084 <span class="keywordflow">for</span>(<span class="keywordtype">int</span> ih=0;ih&lt;calo-&gt;istrip;ih++){
106     <a name="l00085"></a>00085 <span class="keywordtype">int</span> iv=-1;
107     <a name="l00086"></a>00086 <span class="keywordtype">int</span> ip=-1;
108     <a name="l00087"></a>00087 <span class="keywordtype">int</span> is=-1;
109     <a name="l00088"></a>00088 calo-&gt;DecodeEstrip(ih,iv,ip,is);
110     <a name="l00089"></a>00089 <span class="keywordflow">if</span>(iv==view){
111     <a name="l00090"></a>00090 cstrip.Set(view,ip,is);
112     <a name="l00091"></a>00091 <span class="keywordtype">float</span> fxy = 0.;
113     <a name="l00092"></a>00092 <span class="keywordflow">if</span>( !view ) fxy = cstrip.GetX();
114     <a name="l00093"></a>00093 <span class="keywordflow">else</span> fxy = cstrip.GetY();
115     <a name="l00094"></a>00094 <span class="keywordtype">float</span> fz = cstrip.GetZ();
116     <a name="l00095"></a>00095 <span class="keywordtype">float</span> fq = cstrip.GetE();
117     <a name="l00096"></a>00096 <span class="keywordflow">if</span>(fq&gt;0)<a class="code" href="classCaloEvent.html#fd2bc3739bbbf688c9323d8a96e17c5f">Add</a>(is,ip,fq,fz,fxy);
118     <a name="l00097"></a>00097 }
119     <a name="l00098"></a>00098 }
120     <a name="l00099"></a>00099
121     <a name="l00100"></a>00100 };
122     <a name="l00101"></a><a class="code" href="classCaloEvent.html#eb4e3acbb34abd7b51794c6aee4852e5">00101</a> <span class="keywordtype">void</span> <a class="code" href="classCaloEvent.html#eb4e3acbb34abd7b51794c6aee4852e5">Set</a>(CaloLevel1* calo, <span class="keywordtype">int</span> view){ <a class="code" href="classCaloEvent.html#eb4e3acbb34abd7b51794c6aee4852e5">Set</a>(calo,view,0); };
123     <a name="l00102"></a>00102
124     <a name="l00103"></a>00103
125     <a name="l00104"></a>00104 };
126     <a name="l00105"></a>00105 <span class="comment">//===============================================================================</span>
127     <a name="l00106"></a>00106 <span class="comment">//</span>
128     <a name="l00107"></a>00107 <span class="comment">//</span>
129     <a name="l00108"></a>00108 <span class="comment">//</span>
130     <a name="l00109"></a>00109 <span class="comment">//</span>
131     <a name="l00110"></a>00110 <span class="comment">//===============================================================================</span>
132     <a name="l00111"></a><a class="code" href="classCaloAxis.html">00111</a> <span class="keyword">class </span><a class="code" href="classCaloAxis.html">CaloAxis</a> {
133     <a name="l00112"></a>00112
134     <a name="l00113"></a>00113 <span class="keyword">public</span>:
135     <a name="l00114"></a>00114
136     <a name="l00115"></a><a class="code" href="classCaloAxis.html#fb48c08633030431b7e509ff1b34558d">00115</a> <a class="code" href="classCaloEvent.html">CaloEvent</a> *<a class="code" href="classCaloAxis.html#fb48c08633030431b7e509ff1b34558d">cevent</a>;
137     <a name="l00116"></a>00116
138     <a name="l00117"></a>00117 <span class="comment">// ------------------------------</span>
139     <a name="l00118"></a>00118 <span class="comment">// selected points along the axis</span>
140     <a name="l00119"></a>00119 <span class="comment">// ------------------------------</span>
141     <a name="l00120"></a>00120 <span class="comment">// fitted points</span>
142     <a name="l00121"></a><a class="code" href="classCaloAxis.html#bea775e33c9f4ea5c308eace6ac4a513">00121</a> vector&lt;float&gt; <a class="code" href="classCaloAxis.html#bea775e33c9f4ea5c308eace6ac4a513">x</a>;<span class="comment">// independent coordinates (z)</span>
143     <a name="l00122"></a><a class="code" href="classCaloAxis.html#dc1ce8c49916d35fb3bfcc29e03f8905">00122</a> vector&lt;float&gt; <a class="code" href="classCaloAxis.html#dc1ce8c49916d35fb3bfcc29e03f8905">y</a>;<span class="comment">// dependente coordinate (x/y)</span>
144     <a name="l00123"></a><a class="code" href="classCaloAxis.html#a14b06b7a9cd75a27ce585e250efdc63">00123</a> vector&lt;float&gt; <a class="code" href="classCaloAxis.html#a14b06b7a9cd75a27ce585e250efdc63">w</a>;<span class="comment">// fit weight</span>
145     <a name="l00124"></a><a class="code" href="classCaloAxis.html#06da2809e33c9fa67e151d210fdab8fe">00124</a> vector&lt;float&gt; <a class="code" href="classCaloAxis.html#06da2809e33c9fa67e151d210fdab8fe">q</a>;<span class="comment">// signal</span>
146     <a name="l00125"></a>00125
147     <a name="l00126"></a>00126 <span class="comment">// ------------------------------</span>
148     <a name="l00127"></a>00127 <span class="comment">// axis parameters</span>
149     <a name="l00128"></a>00128 <span class="comment">// ------------------------------</span>
150     <a name="l00129"></a><a class="code" href="classCaloAxis.html#e4b2ce253db4df8576ccc1f122a7e432">00129</a> Float_t <a class="code" href="classCaloAxis.html#e4b2ce253db4df8576ccc1f122a7e432">par</a>[2];
151     <a name="l00130"></a><a class="code" href="classCaloAxis.html#8ffb46af208ad9a97e6d8feb8d9a088a">00130</a> Float_t <a class="code" href="classCaloAxis.html#8ffb46af208ad9a97e6d8feb8d9a088a">covpar</a>[2][2];
152     <a name="l00131"></a>00131
153     <a name="l00132"></a>00132 <span class="comment">// ------------------------------</span>
154     <a name="l00133"></a>00133 <span class="comment">// general variables</span>
155     <a name="l00134"></a>00134 <span class="comment">// ------------------------------</span>
156     <a name="l00135"></a><a class="code" href="classCaloAxis.html#575a6d4f2ab14d0039902b0818395222">00135</a> <span class="keywordtype">float</span> <a class="code" href="classCaloAxis.html#575a6d4f2ab14d0039902b0818395222">cog</a>[22];<span class="comment">// baricenter</span>
157     <a name="l00136"></a><a class="code" href="classCaloAxis.html#2ff6400f0393e7ab05627a438c6eeca5">00136</a> <span class="keywordtype">float</span> <a class="code" href="classCaloAxis.html#2ff6400f0393e7ab05627a438c6eeca5">qpl</a>[22]; <span class="comment">// charge in each plane</span>
158     <a name="l00137"></a>00137
159     <a name="l00138"></a>00138 <span class="comment">// ------------------------------</span>
160     <a name="l00139"></a>00139 <span class="comment">// objects to draw the axis</span>
161     <a name="l00140"></a>00140 <span class="comment">// ------------------------------</span>
162     <a name="l00141"></a><a class="code" href="classCaloAxis.html#9e4372d0b9d9a3ec17fc2e1570039273">00141</a> TPolyMarker *<a class="code" href="classCaloAxis.html#9e4372d0b9d9a3ec17fc2e1570039273">sel</a>;
163     <a name="l00142"></a><a class="code" href="classCaloAxis.html#7ca3199fb8575c566821219253d02ccf">00142</a> TLine *<a class="code" href="classCaloAxis.html#7ca3199fb8575c566821219253d02ccf">axis</a>;
164     <a name="l00143"></a>00143
165     <a name="l00144"></a>00144
166     <a name="l00145"></a>00145 <span class="comment">// ------------------------------</span>
167     <a name="l00146"></a>00146 <span class="comment">// methods</span>
168     <a name="l00147"></a>00147 <span class="comment">// ------------------------------</span>
169     <a name="l00148"></a><a class="code" href="classCaloAxis.html#cf42fbb09c6992fb4292da8ec22ef551">00148</a> <span class="keywordtype">void</span> <a class="code" href="classCaloAxis.html#cf42fbb09c6992fb4292da8ec22ef551">Init</a>(){
170     <a name="l00149"></a>00149 <span class="keywordflow">for</span>(Int_t i=0; i&lt;2;i++){
171     <a name="l00150"></a>00150 <a class="code" href="classCaloAxis.html#e4b2ce253db4df8576ccc1f122a7e432">par</a>[i]=0;
172     <a name="l00151"></a>00151 <span class="keywordflow">for</span>(Int_t j=0; j&lt;2;j++)<a class="code" href="classCaloAxis.html#8ffb46af208ad9a97e6d8feb8d9a088a">covpar</a>[i][j]=0;
173     <a name="l00152"></a>00152 };
174     <a name="l00153"></a>00153 <a class="code" href="classCaloAxis.html#9e4372d0b9d9a3ec17fc2e1570039273">sel</a> = 0;
175     <a name="l00154"></a>00154 <a class="code" href="classCaloAxis.html#7ca3199fb8575c566821219253d02ccf">axis</a> = 0;
176     <a name="l00155"></a>00155 <span class="keywordflow">for</span>(Int_t i=0; i&lt;22; i++){
177     <a name="l00156"></a>00156 <a class="code" href="classCaloAxis.html#575a6d4f2ab14d0039902b0818395222">cog</a>[i]=0;
178     <a name="l00157"></a>00157 <a class="code" href="classCaloAxis.html#2ff6400f0393e7ab05627a438c6eeca5">qpl</a>[i]=0;
179     <a name="l00158"></a>00158 }
180     <a name="l00159"></a>00159 <span class="comment">// if(cevent)cevent-&gt;Reset();</span>
181     <a name="l00160"></a>00160 }
182     <a name="l00161"></a>00161
183     <a name="l00162"></a><a class="code" href="classCaloAxis.html#00e32ad98be21c697b8a0671b90e8dd1">00162</a> <span class="keywordtype">void</span> <a class="code" href="classCaloAxis.html#00e32ad98be21c697b8a0671b90e8dd1">Set</a>(CaloLevel1* calo, Int_t view){ <a class="code" href="classCaloAxis.html#fb48c08633030431b7e509ff1b34558d">cevent</a>-&gt;<a class="code" href="classCaloEvent.html#538af1924be400124d003802302e71dd">Set</a>(calo,view,0); };
184     <a name="l00163"></a><a class="code" href="classCaloAxis.html#ba39f6561558754993ca0cb221f94053">00163</a> <span class="keywordtype">void</span> <a class="code" href="classCaloAxis.html#ba39f6561558754993ca0cb221f94053">Set</a>(CaloLevel1* calo, Int_t view, Bool_t usemechanicalalignement){ <a class="code" href="classCaloAxis.html#fb48c08633030431b7e509ff1b34558d">cevent</a>-&gt;<a class="code" href="classCaloEvent.html#538af1924be400124d003802302e71dd">Set</a>(calo,view,usemechanicalalignement); };
185     <a name="l00164"></a>00164
186     <a name="l00165"></a><a class="code" href="classCaloAxis.html#5fbc947804114f961214c5c26b8fe927">00165</a> <a class="code" href="classCaloAxis.html#5fbc947804114f961214c5c26b8fe927">CaloAxis</a>(){
187     <a name="l00166"></a>00166 <a class="code" href="classCaloAxis.html#fb48c08633030431b7e509ff1b34558d">cevent</a> = <span class="keyword">new</span> <a class="code" href="classCaloEvent.html">CaloEvent</a>();
188     <a name="l00167"></a>00167 <a class="code" href="classCaloAxis.html#cf42fbb09c6992fb4292da8ec22ef551">Init</a>();
189     <a name="l00168"></a>00168 };
190     <a name="l00169"></a>00169
191     <a name="l00170"></a><a class="code" href="classCaloAxis.html#a0601f823b9ba6d1e359d4f2c61d6be5">00170</a> <a class="code" href="classCaloAxis.html#a0601f823b9ba6d1e359d4f2c61d6be5">CaloAxis</a>(CaloLevel1* calo, Int_t view){
192     <a name="l00171"></a>00171 <a class="code" href="classCaloAxis.html#fb48c08633030431b7e509ff1b34558d">cevent</a> = <span class="keyword">new</span> <a class="code" href="classCaloEvent.html">CaloEvent</a>();
193     <a name="l00172"></a>00172 <a class="code" href="classCaloAxis.html#cf42fbb09c6992fb4292da8ec22ef551">Init</a>();
194     <a name="l00173"></a>00173 <a class="code" href="classCaloAxis.html#00e32ad98be21c697b8a0671b90e8dd1">Set</a>(calo,view);
195     <a name="l00174"></a>00174 };
196     <a name="l00175"></a><a class="code" href="classCaloAxis.html#c86a443d7a47406a5bdf5b05afa7d122">00175</a> <a class="code" href="classCaloAxis.html#c86a443d7a47406a5bdf5b05afa7d122">CaloAxis</a>(CaloLevel1* calo, Int_t view, Bool_t usemechanicalalignement){
197     <a name="l00176"></a>00176 <a class="code" href="classCaloAxis.html#fb48c08633030431b7e509ff1b34558d">cevent</a> = <span class="keyword">new</span> <a class="code" href="classCaloEvent.html">CaloEvent</a>();
198     <a name="l00177"></a>00177 <a class="code" href="classCaloAxis.html#cf42fbb09c6992fb4292da8ec22ef551">Init</a>();
199     <a name="l00178"></a>00178 <a class="code" href="classCaloAxis.html#00e32ad98be21c697b8a0671b90e8dd1">Set</a>(calo,view,usemechanicalalignement);
200     <a name="l00179"></a>00179 };
201     <a name="l00180"></a>00180
202     <a name="l00181"></a>00181
203     <a name="l00182"></a><a class="code" href="classCaloAxis.html#a1b225fccc498fdf49c415497d405eac">00182</a> <span class="keywordtype">void</span> <a class="code" href="classCaloAxis.html#a1b225fccc498fdf49c415497d405eac">Reset</a>(){
204     <a name="l00183"></a>00183 <span class="comment">// cout &lt;&lt; " CaloAxis::Reset()"&lt;&lt;endl;</span>
205     <a name="l00184"></a>00184 <a class="code" href="classCaloAxis.html#bea775e33c9f4ea5c308eace6ac4a513">x</a>.clear();
206     <a name="l00185"></a>00185 <a class="code" href="classCaloAxis.html#dc1ce8c49916d35fb3bfcc29e03f8905">y</a>.clear();
207     <a name="l00186"></a>00186 <a class="code" href="classCaloAxis.html#a14b06b7a9cd75a27ce585e250efdc63">w</a>.clear();
208     <a name="l00187"></a>00187 <a class="code" href="classCaloAxis.html#06da2809e33c9fa67e151d210fdab8fe">q</a>.clear();
209     <a name="l00188"></a>00188 <span class="comment">// cout &lt;&lt; sel &lt;&lt; endl;</span>
210     <a name="l00189"></a>00189 <span class="keywordflow">if</span>(<a class="code" href="classCaloAxis.html#9e4372d0b9d9a3ec17fc2e1570039273">sel</a>)<a class="code" href="classCaloAxis.html#9e4372d0b9d9a3ec17fc2e1570039273">sel</a>-&gt;Delete();
211     <a name="l00190"></a>00190 <span class="comment">// cout &lt;&lt; axis &lt;&lt; endl;</span>
212     <a name="l00191"></a>00191 <span class="keywordflow">if</span>(<a class="code" href="classCaloAxis.html#7ca3199fb8575c566821219253d02ccf">axis</a>)<a class="code" href="classCaloAxis.html#7ca3199fb8575c566821219253d02ccf">axis</a>-&gt;Delete();
213     <a name="l00192"></a>00192 <span class="comment">// if(cevent)cevent-&gt;Reset();</span>
214     <a name="l00193"></a>00193 <a class="code" href="classCaloAxis.html#cf42fbb09c6992fb4292da8ec22ef551">Init</a>();
215     <a name="l00194"></a>00194 }
216     <a name="l00195"></a>00195
217     <a name="l00196"></a><a class="code" href="classCaloAxis.html#bed1c4a34e0db2caf824f29bbb014789">00196</a> <span class="keywordtype">void</span> <a class="code" href="classCaloAxis.html#bed1c4a34e0db2caf824f29bbb014789">Delete</a>(){ <a class="code" href="classCaloAxis.html#a1b225fccc498fdf49c415497d405eac">Reset</a>(); <span class="keyword">delete</span> <a class="code" href="classCaloAxis.html#fb48c08633030431b7e509ff1b34558d">cevent</a>; };
218     <a name="l00197"></a>00197
219     <a name="l00198"></a><a class="code" href="classCaloAxis.html#12750c5b587e571b7ab632826e9e5e98">00198</a> <a class="code" href="classCaloAxis.html#12750c5b587e571b7ab632826e9e5e98">~CaloAxis</a>(){ <a class="code" href="classCaloAxis.html#bed1c4a34e0db2caf824f29bbb014789">Delete</a>(); };
220     <a name="l00199"></a>00199
221     <a name="l00200"></a>00200
222     <a name="l00201"></a><a class="code" href="classCaloAxis.html#4953d4ace9da1118672a451d66afd01b">00201</a> <span class="keywordtype">void</span> <a class="code" href="classCaloAxis.html#4953d4ace9da1118672a451d66afd01b">Add</a>(<span class="keywordtype">float</span> xin, <span class="keywordtype">float</span> yin){
223     <a name="l00202"></a>00202 <a class="code" href="classCaloAxis.html#bea775e33c9f4ea5c308eace6ac4a513">x</a>.push_back(xin);
224     <a name="l00203"></a>00203 <a class="code" href="classCaloAxis.html#dc1ce8c49916d35fb3bfcc29e03f8905">y</a>.push_back(yin);
225     <a name="l00204"></a>00204 <a class="code" href="classCaloAxis.html#a14b06b7a9cd75a27ce585e250efdc63">w</a>.push_back(1.);
226     <a name="l00205"></a>00205 <a class="code" href="classCaloAxis.html#06da2809e33c9fa67e151d210fdab8fe">q</a>.push_back(0.);
227     <a name="l00206"></a>00206 <span class="comment">/* Int_t nin = x.size(); */</span>
228     <a name="l00207"></a>00207 <span class="comment">/* if(nin&gt;2)Fit(); */</span>
229     <a name="l00208"></a>00208 };
230     <a name="l00209"></a><a class="code" href="classCaloAxis.html#3b38f46913dc5e5ec4b522541822cfae">00209</a> <span class="keywordtype">void</span> <a class="code" href="classCaloAxis.html#3b38f46913dc5e5ec4b522541822cfae">Add</a>(<span class="keywordtype">float</span> xin, <span class="keywordtype">float</span> yin, <span class="keywordtype">float</span> qin){
231     <a name="l00210"></a>00210 <a class="code" href="classCaloAxis.html#bea775e33c9f4ea5c308eace6ac4a513">x</a>.push_back(xin);
232     <a name="l00211"></a>00211 <a class="code" href="classCaloAxis.html#dc1ce8c49916d35fb3bfcc29e03f8905">y</a>.push_back(yin);
233     <a name="l00212"></a>00212 <a class="code" href="classCaloAxis.html#06da2809e33c9fa67e151d210fdab8fe">q</a>.push_back(qin);
234     <a name="l00213"></a>00213 <a class="code" href="classCaloAxis.html#a14b06b7a9cd75a27ce585e250efdc63">w</a>.push_back(1.);
235     <a name="l00214"></a>00214 <span class="comment">/* Int_t nin = x.size(); */</span>
236     <a name="l00215"></a>00215 <span class="comment">/* if(nin&gt;2)Fit(); */</span>
237     <a name="l00216"></a>00216 };
238     <a name="l00217"></a><a class="code" href="classCaloAxis.html#93078ed5c1350094a15e601351f8500e">00217</a> <span class="keywordtype">void</span> <a class="code" href="classCaloAxis.html#93078ed5c1350094a15e601351f8500e">Add</a>(<span class="keywordtype">float</span> xin, <span class="keywordtype">float</span> yin, <span class="keywordtype">float</span> qin, <span class="keywordtype">float</span> win){
239     <a name="l00218"></a>00218 <a class="code" href="classCaloAxis.html#bea775e33c9f4ea5c308eace6ac4a513">x</a>.push_back(xin);
240     <a name="l00219"></a>00219 <a class="code" href="classCaloAxis.html#dc1ce8c49916d35fb3bfcc29e03f8905">y</a>.push_back(yin);
241     <a name="l00220"></a>00220 <a class="code" href="classCaloAxis.html#a14b06b7a9cd75a27ce585e250efdc63">w</a>.push_back(win);
242     <a name="l00221"></a>00221 <a class="code" href="classCaloAxis.html#06da2809e33c9fa67e151d210fdab8fe">q</a>.push_back(qin);
243     <a name="l00222"></a>00222 <span class="comment">/* Int_t nin = x.size(); */</span>
244     <a name="l00223"></a>00223 <span class="comment">/* if(nin&gt;2)Fit(); */</span>
245     <a name="l00224"></a>00224 };
246     <a name="l00225"></a>00225
247     <a name="l00226"></a><a class="code" href="classCaloAxis.html#58d8a61ebffbef5318b9539cd7364633">00226</a> <span class="keywordtype">int</span> <a class="code" href="classCaloAxis.html#58d8a61ebffbef5318b9539cd7364633">Fit</a>(){
248     <a name="l00227"></a>00227
249     <a name="l00228"></a>00228 <span class="comment">// cout &lt;&lt; " CaloAxis::Fit()"&lt;&lt;endl;</span>
250     <a name="l00229"></a>00229 Float_t SSS = 0;
251     <a name="l00230"></a>00230 Float_t SSX = 0;
252     <a name="l00231"></a>00231 Float_t SXX = 0;
253     <a name="l00232"></a>00232
254     <a name="l00233"></a>00233 Float_t SSY = 0;
255     <a name="l00234"></a>00234 Float_t SXY = 0;
256     <a name="l00235"></a>00235
257     <a name="l00236"></a>00236 Int_t np=0;
258     <a name="l00237"></a>00237 <span class="keywordflow">for</span>(Int_t i=0; i&lt;(int)(<a class="code" href="classCaloAxis.html#bea775e33c9f4ea5c308eace6ac4a513">x</a>.size()); i++){
259     <a name="l00238"></a>00238 SSS += <a class="code" href="classCaloAxis.html#a14b06b7a9cd75a27ce585e250efdc63">w</a>[i]*<a class="code" href="classCaloAxis.html#a14b06b7a9cd75a27ce585e250efdc63">w</a>[i];
260     <a name="l00239"></a>00239 SSX += <a class="code" href="classCaloAxis.html#bea775e33c9f4ea5c308eace6ac4a513">x</a>[i]*w[i]*w[i];
261     <a name="l00240"></a>00240 SXX += <a class="code" href="classCaloAxis.html#bea775e33c9f4ea5c308eace6ac4a513">x</a>[i]*<a class="code" href="classCaloAxis.html#bea775e33c9f4ea5c308eace6ac4a513">x</a>[i]*w[i]*w[i];
262     <a name="l00241"></a>00241 SSY += <a class="code" href="classCaloAxis.html#dc1ce8c49916d35fb3bfcc29e03f8905">y</a>[i]*w[i]*w[i];
263     <a name="l00242"></a>00242 SXY += <a class="code" href="classCaloAxis.html#bea775e33c9f4ea5c308eace6ac4a513">x</a>[i]*<a class="code" href="classCaloAxis.html#dc1ce8c49916d35fb3bfcc29e03f8905">y</a>[i]*w[i]*w[i];
264     <a name="l00243"></a>00243 <span class="keywordflow">if</span>(w[i])np++;
265     <a name="l00244"></a>00244 }
266     <a name="l00245"></a>00245 Float_t det = SSS*SXX-SSX*SSX;
267     <a name="l00246"></a>00246 <span class="keywordflow">if</span>(det==0)<span class="keywordflow">return</span> 0;
268     <a name="l00247"></a>00247 <span class="keywordflow">if</span>(np&lt;3)<span class="keywordflow">return</span> 0;
269     <a name="l00248"></a>00248 <span class="comment">// cout &lt;&lt; np &lt;&lt; " points fitted -- ok "&lt;&lt;endl;</span>
270     <a name="l00249"></a>00249 <a class="code" href="classCaloAxis.html#e4b2ce253db4df8576ccc1f122a7e432">par</a>[0] = (SSY*SXX-SXY*SSX)/det;
271     <a name="l00250"></a>00250 <a class="code" href="classCaloAxis.html#e4b2ce253db4df8576ccc1f122a7e432">par</a>[1] = (SXY*SSS-SSY*SSX)/det;
272     <a name="l00251"></a>00251 <span class="keywordflow">return</span> 1;
273     <a name="l00252"></a>00252 };
274     <a name="l00253"></a>00253
275     <a name="l00254"></a><a class="code" href="classCaloAxis.html#7c3a74ca8926244b321e7ef2d12cedd2">00254</a> <span class="keywordtype">void</span> <a class="code" href="classCaloAxis.html#7c3a74ca8926244b321e7ef2d12cedd2">Print</a>(){
276     <a name="l00255"></a>00255 <span class="comment">// useful for debug for(Int_t i=0; i&lt;(int)(x.size()); i++)if(w[i])cout &lt;&lt; x[i] &lt;&lt; " - "&lt;&lt;y[i]&lt;&lt;endl;</span>
277     <a name="l00256"></a>00256 }
278     <a name="l00257"></a>00257
279     <a name="l00258"></a><a class="code" href="classCaloAxis.html#24eb8972ce108383bf4d9ee90f0bff89">00258</a> <span class="keywordtype">float</span> <a class="code" href="classCaloAxis.html#24eb8972ce108383bf4d9ee90f0bff89">GetPar</a>(Int_t i){<span class="keywordflow">return</span> <a class="code" href="classCaloAxis.html#e4b2ce253db4df8576ccc1f122a7e432">par</a>[i];};
280     <a name="l00259"></a>00259
281     <a name="l00260"></a><a class="code" href="classCaloAxis.html#be448c5906fb90329316b50796d2c721">00260</a> <span class="keywordtype">int</span> <a class="code" href="classCaloAxis.html#be448c5906fb90329316b50796d2c721">GetN</a>(){<span class="keywordflow">return</span> (<span class="keywordtype">int</span>)(<a class="code" href="classCaloAxis.html#bea775e33c9f4ea5c308eace6ac4a513">x</a>.size());};
282     <a name="l00261"></a>00261
283     <a name="l00262"></a><a class="code" href="classCaloAxis.html#2973ca2a8afc807f56805c9420b47ed0">00262</a> <span class="keywordtype">float</span>* <a class="code" href="classCaloAxis.html#2973ca2a8afc807f56805c9420b47ed0">GetX</a>(){
284     <a name="l00263"></a>00263 <span class="keywordtype">float</span> *xout = <span class="keyword">new</span> <span class="keywordtype">float</span>[(int)(<a class="code" href="classCaloAxis.html#bea775e33c9f4ea5c308eace6ac4a513">x</a>.size())];
285     <a name="l00264"></a>00264 <span class="keywordflow">for</span>(Int_t i=0; i&lt;(int)(<a class="code" href="classCaloAxis.html#bea775e33c9f4ea5c308eace6ac4a513">x</a>.size()); i++)xout[i]=<a class="code" href="classCaloAxis.html#bea775e33c9f4ea5c308eace6ac4a513">x</a>[i];
286     <a name="l00265"></a>00265 <span class="keywordflow">return</span> xout;
287     <a name="l00266"></a>00266 }
288     <a name="l00267"></a>00267
289     <a name="l00268"></a><a class="code" href="classCaloAxis.html#c3172d5d58795cc3de1522559a231066">00268</a> <span class="keywordtype">float</span>* <a class="code" href="classCaloAxis.html#c3172d5d58795cc3de1522559a231066">GetY</a>(){
290     <a name="l00269"></a>00269 <span class="keywordtype">float</span> *yout = <span class="keyword">new</span> <span class="keywordtype">float</span>[(int)(<a class="code" href="classCaloAxis.html#dc1ce8c49916d35fb3bfcc29e03f8905">y</a>.size())];
291     <a name="l00270"></a>00270 <span class="keywordflow">for</span>(Int_t i=0; i&lt;(int)(<a class="code" href="classCaloAxis.html#dc1ce8c49916d35fb3bfcc29e03f8905">y</a>.size()); i++) yout[i]=<a class="code" href="classCaloAxis.html#dc1ce8c49916d35fb3bfcc29e03f8905">y</a>[i];
292     <a name="l00271"></a>00271 <span class="keywordflow">return</span> yout;
293     <a name="l00272"></a>00272 }
294     <a name="l00273"></a>00273
295     <a name="l00274"></a><a class="code" href="classCaloAxis.html#fa2726b36147a39e2c64d4808aa435e3">00274</a> <span class="keywordtype">float</span>* <a class="code" href="classCaloAxis.html#fa2726b36147a39e2c64d4808aa435e3">GetQ</a>(){
296     <a name="l00275"></a>00275 <span class="keywordtype">float</span> *qout = <span class="keyword">new</span> <span class="keywordtype">float</span>[(int)(<a class="code" href="classCaloAxis.html#06da2809e33c9fa67e151d210fdab8fe">q</a>.size())];
297     <a name="l00276"></a>00276 <span class="keywordflow">for</span>(<span class="keywordtype">int</span> i=0; i&lt;(int)(<a class="code" href="classCaloAxis.html#06da2809e33c9fa67e151d210fdab8fe">q</a>.size()); i++) qout[i]=<a class="code" href="classCaloAxis.html#06da2809e33c9fa67e151d210fdab8fe">q</a>[i];
298     <a name="l00277"></a>00277 <span class="keywordflow">return</span> qout;
299     <a name="l00278"></a>00278 }
300     <a name="l00279"></a>00279
301     <a name="l00280"></a><a class="code" href="classCaloAxis.html#cd5c771ea0ead8eb2428118551358eb0">00280</a> <span class="keywordtype">float</span> <a class="code" href="classCaloAxis.html#cd5c771ea0ead8eb2428118551358eb0">GetChi2</a>(){
302     <a name="l00281"></a>00281 <span class="keywordtype">float</span> chi2=0;
303     <a name="l00282"></a>00282 <span class="keywordtype">int</span> nchi=0;
304     <a name="l00283"></a>00283 <span class="keywordflow">for</span>(Int_t i=0; i&lt;(int)(<a class="code" href="classCaloAxis.html#dc1ce8c49916d35fb3bfcc29e03f8905">y</a>.size()); i++){
305     <a name="l00284"></a>00284 chi2 += <a class="code" href="classCaloAxis.html#a14b06b7a9cd75a27ce585e250efdc63">w</a>[i]*<a class="code" href="classCaloAxis.html#a14b06b7a9cd75a27ce585e250efdc63">w</a>[i]*(<a class="code" href="classCaloAxis.html#dc1ce8c49916d35fb3bfcc29e03f8905">y</a>[i]-<a class="code" href="classCaloAxis.html#e4b2ce253db4df8576ccc1f122a7e432">par</a>[0]-<a class="code" href="classCaloAxis.html#e4b2ce253db4df8576ccc1f122a7e432">par</a>[1]*<a class="code" href="classCaloAxis.html#bea775e33c9f4ea5c308eace6ac4a513">x</a>[i])*(<a class="code" href="classCaloAxis.html#dc1ce8c49916d35fb3bfcc29e03f8905">y</a>[i]-<a class="code" href="classCaloAxis.html#e4b2ce253db4df8576ccc1f122a7e432">par</a>[0]-<a class="code" href="classCaloAxis.html#e4b2ce253db4df8576ccc1f122a7e432">par</a>[1]*<a class="code" href="classCaloAxis.html#bea775e33c9f4ea5c308eace6ac4a513">x</a>[i]);
306     <a name="l00285"></a>00285 nchi += (int)w[i];
307     <a name="l00286"></a>00286 };
308     <a name="l00287"></a>00287
309     <a name="l00288"></a>00288 <span class="keywordflow">if</span>(nchi&lt;3) <span class="keywordflow">return</span> -999; <span class="comment">// original</span>
310     <a name="l00289"></a>00289 <span class="comment">// if(nchi&lt;3 || chi2==0)return -999; // modified by Sergio</span>
311     <a name="l00290"></a>00290
312     <a name="l00291"></a>00291 chi2=chi2/(nchi-2);
313     <a name="l00292"></a>00292 <span class="keywordflow">return</span> chi2;
314     <a name="l00293"></a>00293 }
315     <a name="l00294"></a>00294
316     <a name="l00295"></a><a class="code" href="classCaloAxis.html#0d817c015af5b772b4a0944aa3dbf696">00295</a> <span class="keywordtype">float</span> <a class="code" href="classCaloAxis.html#0d817c015af5b772b4a0944aa3dbf696">GetQaxis</a>(){
317     <a name="l00296"></a>00296 <span class="keywordtype">float</span> qaxis=0;
318     <a name="l00297"></a>00297 <span class="keywordflow">for</span>(Int_t i=0; i&lt;(int)(<a class="code" href="classCaloAxis.html#06da2809e33c9fa67e151d210fdab8fe">q</a>.size()); i++) qaxis +=<a class="code" href="classCaloAxis.html#06da2809e33c9fa67e151d210fdab8fe">q</a>[i];
319     <a name="l00298"></a>00298 <span class="comment">// for(Int_t i=0; i&lt;(int)(q.size()); i++) cout&lt;&lt;q[i]&lt;&lt;endl;</span>
320     <a name="l00299"></a>00299 <span class="keywordflow">return</span> qaxis;
321     <a name="l00300"></a>00300 }
322     <a name="l00301"></a>00301
323     <a name="l00302"></a><a class="code" href="classCaloAxis.html#8bb7108aa0e18c85b471618cfcbba66a">00302</a> <span class="keywordtype">float</span> <a class="code" href="classCaloAxis.html#8bb7108aa0e18c85b471618cfcbba66a">GetQaxis</a>(<span class="keywordtype">float</span> toll){
324     <a name="l00303"></a>00303 <span class="keywordtype">float</span> qaxis=0;
325     <a name="l00304"></a>00304 <span class="keywordflow">for</span>(<span class="keywordtype">int</span> ih=0;ih&lt;<a class="code" href="classCaloAxis.html#fb48c08633030431b7e509ff1b34558d">cevent</a>-&gt;<a class="code" href="classCaloEvent.html#2b2597fe0263fe128101a6ab02520130">GetN</a>(); ih++){
326     <a name="l00305"></a>00305 <span class="keywordtype">float</span> <a class="code" href="classCaloAxis.html#bea775e33c9f4ea5c308eace6ac4a513">x</a> = <a class="code" href="classCaloAxis.html#fb48c08633030431b7e509ff1b34558d">cevent</a>-&gt;<a class="code" href="classCaloEvent.html#b48403ffdb4114d57e4ea4058438e814">xycoord</a>[ih];
327     <a name="l00306"></a>00306 <span class="keywordtype">float</span> z = <a class="code" href="classCaloAxis.html#fb48c08633030431b7e509ff1b34558d">cevent</a>-&gt;<a class="code" href="classCaloEvent.html#340ff2ae2167b37d0b6acc56b1475daf">zcoord</a>[ih];
328     <a name="l00307"></a>00307 <span class="keywordtype">float</span> <a class="code" href="classCaloAxis.html#06da2809e33c9fa67e151d210fdab8fe">q</a> = <a class="code" href="classCaloAxis.html#fb48c08633030431b7e509ff1b34558d">cevent</a>-&gt;<a class="code" href="classCaloEvent.html#688d5d2803d4688d4f1ef72cbaf42dcd">q</a>[ih];
329     <a name="l00308"></a>00308 <span class="comment">//int ip = cevent-&gt;idp[ih];</span>
330     <a name="l00309"></a>00309 <span class="keywordtype">float</span> d = fabs(x-<a class="code" href="classCaloAxis.html#e4b2ce253db4df8576ccc1f122a7e432">par</a>[0]-<a class="code" href="classCaloAxis.html#e4b2ce253db4df8576ccc1f122a7e432">par</a>[1]*z)/sqrt(<a class="code" href="classCaloAxis.html#e4b2ce253db4df8576ccc1f122a7e432">par</a>[1]*<a class="code" href="classCaloAxis.html#e4b2ce253db4df8576ccc1f122a7e432">par</a>[1]+1);
331     <a name="l00310"></a>00310 <span class="keywordflow">if</span>( d &lt; toll )qaxis+=q;
332     <a name="l00311"></a>00311 }
333     <a name="l00312"></a>00312 <span class="keywordflow">return</span> qaxis;
334     <a name="l00313"></a>00313 }
335     <a name="l00314"></a>00314
336     <a name="l00315"></a><a class="code" href="classCaloAxis.html#c59062096fd2f8addafd4c79c4f64382">00315</a> <span class="keywordtype">float</span> <a class="code" href="classCaloAxis.html#c59062096fd2f8addafd4c79c4f64382">GetCOG</a>(Int_t ip){
337     <a name="l00316"></a>00316 <span class="keywordflow">if</span>(ip&gt;=0 &amp;&amp; ip&lt;22)<span class="keywordflow">return</span> <a class="code" href="classCaloAxis.html#575a6d4f2ab14d0039902b0818395222">cog</a>[ip];
338     <a name="l00317"></a>00317 <span class="keywordflow">else</span> <span class="keywordflow">return</span> 0;
339     <a name="l00318"></a>00318 }
340     <a name="l00319"></a>00319
341     <a name="l00320"></a><a class="code" href="classCaloAxis.html#8c4cd0cb64825ff08e0f91e1759d4fa2">00320</a> <span class="keywordtype">float</span> <a class="code" href="classCaloAxis.html#8c4cd0cb64825ff08e0f91e1759d4fa2">GetQ</a>(Int_t ip){
342     <a name="l00321"></a>00321 <span class="keywordflow">if</span>(ip&gt;=0 &amp;&amp; ip&lt;22)<span class="keywordflow">return</span> <a class="code" href="classCaloAxis.html#2ff6400f0393e7ab05627a438c6eeca5">qpl</a>[ip];
343     <a name="l00322"></a>00322 <span class="keywordflow">else</span> <span class="keywordflow">return</span> 0;
344     <a name="l00323"></a>00323 }
345     <a name="l00324"></a>00324
346     <a name="l00325"></a><a class="code" href="classCaloAxis.html#79e084034e3a504494b41b8d14a210f1">00325</a> <span class="keywordtype">float</span> <a class="code" href="classCaloAxis.html#79e084034e3a504494b41b8d14a210f1">GetQstrip</a>(){
347     <a name="l00326"></a>00326
348     <a name="l00327"></a>00327 <span class="keywordtype">float</span> qaxis=0;
349     <a name="l00328"></a>00328 <span class="keywordflow">for</span>(Int_t i=0; i&lt;(int)(<a class="code" href="classCaloAxis.html#06da2809e33c9fa67e151d210fdab8fe">q</a>.size()); i++) qaxis +=<a class="code" href="classCaloAxis.html#06da2809e33c9fa67e151d210fdab8fe">q</a>[i];
350     <a name="l00329"></a>00329 <span class="keywordflow">if</span>(<a class="code" href="classCaloAxis.html#06da2809e33c9fa67e151d210fdab8fe">q</a>.size()==0) <span class="keywordflow">return</span> 0.;
351     <a name="l00330"></a>00330 <span class="keywordflow">return</span> qaxis/<a class="code" href="classCaloAxis.html#06da2809e33c9fa67e151d210fdab8fe">q</a>.size();
352     <a name="l00331"></a>00331 }
353     <a name="l00332"></a>00332
354     <a name="l00333"></a><a class="code" href="classCaloAxis.html#83364a0c2e1fa2a8cc1a1ab94cf9c05b">00333</a> <span class="keywordtype">float</span> <a class="code" href="classCaloAxis.html#83364a0c2e1fa2a8cc1a1ab94cf9c05b">GetYfit</a>(<span class="keywordtype">float</span> <a class="code" href="classCaloAxis.html#bea775e33c9f4ea5c308eace6ac4a513">x</a>){
355     <a name="l00334"></a>00334 <span class="keywordflow">return</span> <a class="code" href="classCaloAxis.html#e4b2ce253db4df8576ccc1f122a7e432">par</a>[0]+x*<a class="code" href="classCaloAxis.html#e4b2ce253db4df8576ccc1f122a7e432">par</a>[1];
356     <a name="l00335"></a>00335 }
357     <a name="l00336"></a>00336
358     <a name="l00337"></a><a class="code" href="classCaloAxis.html#9cb7da77d46cacded3ecb698d8402c75">00337</a> <span class="keywordtype">float</span> <a class="code" href="classCaloAxis.html#9cb7da77d46cacded3ecb698d8402c75">GetXmin</a>(){
359     <a name="l00338"></a>00338 <span class="keywordtype">float</span> zmin=9999;
360     <a name="l00339"></a>00339 <span class="keywordflow">for</span>(Int_t i=0; i&lt;(int)(<a class="code" href="classCaloAxis.html#bea775e33c9f4ea5c308eace6ac4a513">x</a>.size()); i++) <span class="keywordflow">if</span>(<a class="code" href="classCaloAxis.html#bea775e33c9f4ea5c308eace6ac4a513">x</a>[i]&lt;zmin) zmin = <a class="code" href="classCaloAxis.html#bea775e33c9f4ea5c308eace6ac4a513">x</a>[i];
361     <a name="l00340"></a>00340 <span class="keywordflow">return</span> zmin;
362     <a name="l00341"></a>00341 }
363     <a name="l00342"></a><a class="code" href="classCaloAxis.html#2b80915852cae6f6119b0b66ddf6da18">00342</a> <span class="keywordtype">float</span> <a class="code" href="classCaloAxis.html#2b80915852cae6f6119b0b66ddf6da18">GetXmax</a>(){
364     <a name="l00343"></a>00343 <span class="keywordtype">float</span> zmax=-9999;
365     <a name="l00344"></a>00344 <span class="keywordflow">for</span>(Int_t i=0; i&lt;(int)(<a class="code" href="classCaloAxis.html#bea775e33c9f4ea5c308eace6ac4a513">x</a>.size()); i++) <span class="keywordflow">if</span>(<a class="code" href="classCaloAxis.html#bea775e33c9f4ea5c308eace6ac4a513">x</a>[i]&gt;zmax) zmax = <a class="code" href="classCaloAxis.html#bea775e33c9f4ea5c308eace6ac4a513">x</a>[i];
366     <a name="l00345"></a>00345 <span class="keywordflow">return</span> zmax;
367     <a name="l00346"></a>00346
368     <a name="l00347"></a>00347 }
369     <a name="l00348"></a>00348
370     <a name="l00349"></a>00349
371     <a name="l00350"></a><a class="code" href="classCaloAxis.html#42f57a3a58bfeda6e03b25a6ba875fb6">00350</a> <span class="keywordtype">void</span> <a class="code" href="classCaloAxis.html#42f57a3a58bfeda6e03b25a6ba875fb6">Draw</a>(){
372     <a name="l00351"></a>00351
373     <a name="l00352"></a>00352 <span class="comment">// cout &lt;&lt; " CaloAxis::Draw()"&lt;&lt;endl;</span>
374     <a name="l00353"></a>00353 <span class="keywordflow">if</span>(<a class="code" href="classCaloAxis.html#be448c5906fb90329316b50796d2c721">GetN</a>()&lt;3)<span class="keywordflow">return</span>;
375     <a name="l00354"></a>00354
376     <a name="l00355"></a>00355 <a class="code" href="classCaloAxis.html#9e4372d0b9d9a3ec17fc2e1570039273">sel</a> = <span class="keyword">new</span> TPolyMarker(<a class="code" href="classCaloAxis.html#be448c5906fb90329316b50796d2c721">GetN</a>(),<a class="code" href="classCaloAxis.html#c3172d5d58795cc3de1522559a231066">GetY</a>(),<a class="code" href="classCaloAxis.html#2973ca2a8afc807f56805c9420b47ed0">GetX</a>());
377     <a name="l00356"></a>00356 <a class="code" href="classCaloAxis.html#9e4372d0b9d9a3ec17fc2e1570039273">sel</a>-&gt;SetMarkerSize(0.3);
378     <a name="l00357"></a>00357 <a class="code" href="classCaloAxis.html#9e4372d0b9d9a3ec17fc2e1570039273">sel</a>-&gt;SetMarkerColor(5);
379     <a name="l00358"></a>00358
380     <a name="l00359"></a>00359 <a class="code" href="classCaloAxis.html#7ca3199fb8575c566821219253d02ccf">axis</a> = <span class="keyword">new</span> TLine( <a class="code" href="classCaloAxis.html#e4b2ce253db4df8576ccc1f122a7e432">par</a>[0]+<a class="code" href="classCaloAxis.html#9cb7da77d46cacded3ecb698d8402c75">GetXmin</a>()*<a class="code" href="classCaloAxis.html#e4b2ce253db4df8576ccc1f122a7e432">par</a>[1], <a class="code" href="classCaloAxis.html#9cb7da77d46cacded3ecb698d8402c75">GetXmin</a>(), <a class="code" href="classCaloAxis.html#e4b2ce253db4df8576ccc1f122a7e432">par</a>[0]+<a class="code" href="classCaloAxis.html#2b80915852cae6f6119b0b66ddf6da18">GetXmax</a>()*<a class="code" href="classCaloAxis.html#e4b2ce253db4df8576ccc1f122a7e432">par</a>[1], <a class="code" href="classCaloAxis.html#2b80915852cae6f6119b0b66ddf6da18">GetXmax</a>() );
381     <a name="l00360"></a>00360 <a class="code" href="classCaloAxis.html#7ca3199fb8575c566821219253d02ccf">axis</a>-&gt;SetLineWidth(1);
382     <a name="l00361"></a>00361 <a class="code" href="classCaloAxis.html#7ca3199fb8575c566821219253d02ccf">axis</a>-&gt;SetLineColor(5);
383     <a name="l00362"></a>00362
384     <a name="l00363"></a>00363 <span class="comment">// cout &lt;&lt; sel &lt;&lt; endl;</span>
385     <a name="l00364"></a>00364 <span class="comment">// cout &lt;&lt; axis &lt;&lt; endl;</span>
386     <a name="l00365"></a>00365
387     <a name="l00366"></a>00366 <a class="code" href="classCaloAxis.html#9e4372d0b9d9a3ec17fc2e1570039273">sel</a>-&gt;Draw();
388     <a name="l00367"></a>00367 <a class="code" href="classCaloAxis.html#7ca3199fb8575c566821219253d02ccf">axis</a>-&gt;Draw();
389     <a name="l00368"></a>00368 };
390     <a name="l00369"></a>00369
391     <a name="l00370"></a>00370
392     <a name="l00371"></a>00371 <span class="comment">//-----------------------------------------------------------------</span>
393     <a name="l00372"></a>00372 <span class="comment">// fit the axis (optimized for non-interacting patterns)</span>
394     <a name="l00373"></a>00373 <span class="comment">//-----------------------------------------------------------------</span>
395     <a name="l00374"></a><a class="code" href="classCaloAxis.html#768db65112c0f2ecc7b0561772004e0f">00374</a> <span class="keywordtype">int</span> <a class="code" href="classCaloAxis.html#768db65112c0f2ecc7b0561772004e0f">FitAxis</a>( CaloLevel1* calo , Int_t view , Float_t toll , Bool_t usemechanicalalignement ){
396     <a name="l00375"></a>00375
397     <a name="l00376"></a>00376 <span class="comment">// cout &lt;&lt; "CaloAxis::FitAxis(...)"&lt;&lt;endl;</span>
398     <a name="l00377"></a>00377
399     <a name="l00378"></a>00378 <a class="code" href="classCaloAxis.html#00e32ad98be21c697b8a0671b90e8dd1">Set</a>(calo,view,usemechanicalalignement);
400     <a name="l00379"></a>00379
401     <a name="l00380"></a>00380 <span class="comment">// ------------------------------</span>
402     <a name="l00381"></a>00381 <span class="comment">// first : all hits included, w=1</span>
403     <a name="l00382"></a>00382 <span class="comment">// ------------------------------</span>
404     <a name="l00383"></a>00383 <span class="keywordflow">for</span>(<span class="keywordtype">int</span> ih=0;ih&lt;<a class="code" href="classCaloAxis.html#fb48c08633030431b7e509ff1b34558d">cevent</a>-&gt;<a class="code" href="classCaloEvent.html#2b2597fe0263fe128101a6ab02520130">GetN</a>(); ih++){
405     <a name="l00384"></a>00384 <span class="keywordtype">float</span> <a class="code" href="classCaloAxis.html#bea775e33c9f4ea5c308eace6ac4a513">x</a> = <a class="code" href="classCaloAxis.html#fb48c08633030431b7e509ff1b34558d">cevent</a>-&gt;<a class="code" href="classCaloEvent.html#b48403ffdb4114d57e4ea4058438e814">xycoord</a>[ih];
406     <a name="l00385"></a>00385 <span class="keywordtype">float</span> z = <a class="code" href="classCaloAxis.html#fb48c08633030431b7e509ff1b34558d">cevent</a>-&gt;<a class="code" href="classCaloEvent.html#340ff2ae2167b37d0b6acc56b1475daf">zcoord</a>[ih];
407     <a name="l00386"></a>00386 <span class="keywordtype">float</span> <a class="code" href="classCaloAxis.html#06da2809e33c9fa67e151d210fdab8fe">q</a> = <a class="code" href="classCaloAxis.html#fb48c08633030431b7e509ff1b34558d">cevent</a>-&gt;<a class="code" href="classCaloEvent.html#688d5d2803d4688d4f1ef72cbaf42dcd">q</a>[ih];
408     <a name="l00387"></a>00387 <a class="code" href="classCaloAxis.html#4953d4ace9da1118672a451d66afd01b">Add</a>(z,x,q);
409     <a name="l00388"></a>00388 }
410     <a name="l00389"></a>00389
411     <a name="l00390"></a>00390 <span class="keywordflow">if</span>(<a class="code" href="classCaloAxis.html#be448c5906fb90329316b50796d2c721">GetN</a>()&lt;3)<span class="keywordflow">return</span> 0;
412     <a name="l00391"></a>00391 <a class="code" href="classCaloAxis.html#58d8a61ebffbef5318b9539cd7364633">Fit</a>();
413     <a name="l00392"></a>00392 <span class="comment">// useful for debug cout &lt;&lt; "(1)"&lt;&lt;endl;</span>
414     <a name="l00393"></a>00393 <span class="comment">// cout &lt;&lt; " n. hits :"&lt;&lt;GetN()&lt;&lt;endl;</span>
415     <a name="l00394"></a>00394 <span class="comment">// cout &lt;&lt; " P0 P1 :"&lt;&lt;par[0]&lt;&lt; " "&lt;&lt;par[1]&lt;&lt;endl;</span>
416     <a name="l00395"></a>00395 <span class="comment">// cout &lt;&lt; " Chi2 :"&lt;&lt;GetChi2()&lt;&lt;endl;</span>
417     <a name="l00396"></a>00396 <span class="comment">// cout &lt;&lt; " Q axis :"&lt;&lt;GetQaxis()&lt;&lt;endl;</span>
418     <a name="l00397"></a>00397 <span class="comment">// cout &lt;&lt; " Q strip :"&lt;&lt;GetQstrip()&lt;&lt;endl;</span>
419     <a name="l00398"></a>00398
420     <a name="l00399"></a>00399 <span class="keywordtype">float</span> P0;
421     <a name="l00400"></a>00400 <span class="keywordtype">float</span> P1;
422     <a name="l00401"></a>00401 <span class="keywordtype">int</span> dof;
423     <a name="l00402"></a>00402
424     <a name="l00403"></a>00403 <span class="comment">// ---------------------------</span>
425     <a name="l00404"></a>00404 <span class="comment">// second : iteration</span>
426     <a name="l00405"></a>00405 <span class="comment">// ---------------------------</span>
427     <a name="l00406"></a>00406 <span class="keywordtype">int</span> niter=0; <span class="comment">// number of iterations</span>
428     <a name="l00407"></a>00407 <span class="keywordtype">int</span> pok[22];
429     <a name="l00408"></a>00408
430     <a name="l00409"></a>00409 <span class="comment">// prova</span>
431     <a name="l00410"></a>00410 <span class="comment">// pesi lineari</span>
432     <a name="l00411"></a>00411 <span class="keywordtype">float</span> W22=1.;
433     <a name="l00412"></a>00412 <span class="keywordtype">float</span> W1=10.;
434     <a name="l00413"></a>00413 <span class="keywordtype">float</span> b=(W22-W1)/(22.-1.);<span class="comment">// (w22 - w1)/(22-1)</span>
435     <a name="l00414"></a>00414 <span class="keywordtype">float</span> a=W1-b;
436     <a name="l00415"></a>00415
437     <a name="l00416"></a>00416 <span class="keywordflow">for</span>(<span class="keywordtype">int</span> i=0; i&lt;22; i++)pok[i]=0;
438     <a name="l00417"></a>00417 <span class="keywordflow">do</span>{
439     <a name="l00418"></a>00418 niter++;
440     <a name="l00419"></a>00419 <span class="keywordtype">float</span> ttoll = toll; <span class="comment">//tolerance (cm)</span>
441     <a name="l00420"></a>00420 <span class="keywordflow">if</span>(niter==1)ttoll = 10.*toll;
442     <a name="l00421"></a>00421 dof = <a class="code" href="classCaloAxis.html#be448c5906fb90329316b50796d2c721">GetN</a>(); <span class="comment">//previous fit</span>
443     <a name="l00422"></a>00422 P0 = <a class="code" href="classCaloAxis.html#e4b2ce253db4df8576ccc1f122a7e432">par</a>[0];
444     <a name="l00423"></a>00423 P1 = <a class="code" href="classCaloAxis.html#e4b2ce253db4df8576ccc1f122a7e432">par</a>[1];
445     <a name="l00424"></a>00424 <a class="code" href="classCaloAxis.html#a1b225fccc498fdf49c415497d405eac">Reset</a>(); <span class="comment">//fit reset</span>
446     <a name="l00425"></a>00425 <span class="keywordflow">for</span>(<span class="keywordtype">int</span> ih=0;ih&lt;<a class="code" href="classCaloAxis.html#fb48c08633030431b7e509ff1b34558d">cevent</a>-&gt;<a class="code" href="classCaloEvent.html#2b2597fe0263fe128101a6ab02520130">GetN</a>(); ih++){<span class="comment">//loop over selected hits</span>
447     <a name="l00426"></a>00426 <span class="keywordtype">float</span> <a class="code" href="classCaloAxis.html#bea775e33c9f4ea5c308eace6ac4a513">x</a> = <a class="code" href="classCaloAxis.html#fb48c08633030431b7e509ff1b34558d">cevent</a>-&gt;<a class="code" href="classCaloEvent.html#b48403ffdb4114d57e4ea4058438e814">xycoord</a>[ih];
448     <a name="l00427"></a>00427 <span class="keywordtype">float</span> z = <a class="code" href="classCaloAxis.html#fb48c08633030431b7e509ff1b34558d">cevent</a>-&gt;<a class="code" href="classCaloEvent.html#340ff2ae2167b37d0b6acc56b1475daf">zcoord</a>[ih];
449     <a name="l00428"></a>00428 <span class="keywordtype">float</span> <a class="code" href="classCaloAxis.html#06da2809e33c9fa67e151d210fdab8fe">q</a> = <a class="code" href="classCaloAxis.html#fb48c08633030431b7e509ff1b34558d">cevent</a>-&gt;<a class="code" href="classCaloEvent.html#688d5d2803d4688d4f1ef72cbaf42dcd">q</a>[ih];
450     <a name="l00429"></a>00429 <span class="keywordtype">int</span> ip = <a class="code" href="classCaloAxis.html#fb48c08633030431b7e509ff1b34558d">cevent</a>-&gt;<a class="code" href="classCaloEvent.html#3ad199ffd33f6300bff64e854c6ad652">idp</a>[ih];
451     <a name="l00430"></a>00430 <span class="keywordtype">float</span> d = fabs(x-P0-P1*z)/sqrt(P1*P1+1);<span class="comment">//distance to axis</span>
452     <a name="l00431"></a>00431 <span class="comment">// cout&lt;&lt;d&lt;&lt;endl;</span>
453     <a name="l00432"></a>00432 <span class="keywordflow">if</span>( d &lt; ttoll &amp;&amp; (niter==1 || (niter&gt;1 &amp;&amp; pok[ip]==1)) ){
454     <a name="l00433"></a>00433 <span class="comment">// if( d &lt; ttoll ){</span>
455     <a name="l00434"></a>00434 <span class="keywordtype">float</span> W=a+b*(ip+1);
456     <a name="l00435"></a>00435 <span class="comment">// cout &lt;&lt; ip &lt;&lt; " "&lt;&lt;W&lt;&lt;endl;</span>
457     <a name="l00436"></a>00436 <a class="code" href="classCaloAxis.html#4953d4ace9da1118672a451d66afd01b">Add</a>(z,x,q,W);
458     <a name="l00437"></a>00437 pok[ip]=1;
459     <a name="l00438"></a>00438 }
460     <a name="l00439"></a>00439 }
461     <a name="l00440"></a>00440 <span class="comment">// break the track if more than 3 planes are missing</span>
462     <a name="l00441"></a>00441 <span class="keywordtype">int</span> nmissing = 0;
463     <a name="l00442"></a>00442 <span class="keywordflow">for</span>(<span class="keywordtype">int</span> ip=0; ip&lt;22; ip++){
464     <a name="l00443"></a>00443 <span class="keywordflow">if</span>(pok[ip]==0){
465     <a name="l00444"></a>00444 nmissing++;
466     <a name="l00445"></a>00445 }<span class="keywordflow">else</span>{
467     <a name="l00446"></a>00446 nmissing = 0;
468     <a name="l00447"></a>00447 }
469     <a name="l00448"></a>00448 <span class="keywordflow">if</span>(nmissing==3){
470     <a name="l00449"></a>00449 <span class="keywordflow">for</span>(<span class="keywordtype">int</span> ipp=ip+1; ipp&lt;22; ipp++)pok[ipp]=0;
471     <a name="l00450"></a>00450 <span class="keywordflow">break</span>;
472     <a name="l00451"></a>00451 }
473     <a name="l00452"></a>00452 }
474     <a name="l00453"></a>00453
475     <a name="l00454"></a>00454
476     <a name="l00455"></a>00455 <span class="keywordflow">if</span>(niter==100)<span class="keywordflow">break</span>;
477     <a name="l00456"></a>00456 <span class="keywordflow">if</span>(<a class="code" href="classCaloAxis.html#be448c5906fb90329316b50796d2c721">GetN</a>()&lt;3)<span class="keywordflow">break</span>;
478     <a name="l00457"></a>00457 <a class="code" href="classCaloAxis.html#58d8a61ebffbef5318b9539cd7364633">Fit</a>();
479     <a name="l00458"></a>00458 }<span class="keywordflow">while</span>(niter==1 || <a class="code" href="classCaloAxis.html#be448c5906fb90329316b50796d2c721">GetN</a>() != dof);
480     <a name="l00459"></a>00459 <span class="comment">// cout &lt;&lt; " &gt;&gt;&gt; "&lt;&lt;GetN()&lt;&lt;endl;</span>
481     <a name="l00460"></a>00460 <span class="keywordflow">if</span>(<a class="code" href="classCaloAxis.html#be448c5906fb90329316b50796d2c721">GetN</a>()&lt;3)<span class="keywordflow">return</span> 0;
482     <a name="l00461"></a>00461 <a class="code" href="classCaloAxis.html#58d8a61ebffbef5318b9539cd7364633">Fit</a>();
483     <a name="l00462"></a>00462 <span class="comment">// cout &lt;&lt; "(2)"&lt;&lt;endl;</span>
484     <a name="l00463"></a>00463 <span class="comment">// cout &lt;&lt; " n. hits :"&lt;&lt;GetN()&lt;&lt;endl;</span>
485     <a name="l00464"></a>00464 <span class="comment">// cout &lt;&lt; " P0 P1 :"&lt;&lt;par[0]&lt;&lt; " "&lt;&lt;par[1]&lt;&lt;endl;</span>
486     <a name="l00465"></a>00465 <span class="comment">// cout &lt;&lt; " Chi2 :"&lt;&lt;GetChi2()&lt;&lt;endl;</span>
487     <a name="l00466"></a>00466 <span class="comment">// cout &lt;&lt; " Q axis :"&lt;&lt;GetQaxis()&lt;&lt;endl;</span>
488     <a name="l00467"></a>00467 <span class="comment">// cout &lt;&lt; " Q strip :"&lt;&lt;GetQstrip()&lt;&lt;endl;</span>
489     <a name="l00468"></a>00468
490     <a name="l00469"></a>00469 <span class="comment">// ---------------------------------------------</span>
491     <a name="l00470"></a>00470 <span class="comment">// third : selecting only closest strip-clusters</span>
492     <a name="l00471"></a>00471 <span class="comment">// and fit baricenters</span>
493     <a name="l00472"></a>00472 <span class="comment">// ---------------------------------------------</span>
494     <a name="l00473"></a>00473 P0 = <a class="code" href="classCaloAxis.html#e4b2ce253db4df8576ccc1f122a7e432">par</a>[0];
495     <a name="l00474"></a>00474 P1 = <a class="code" href="classCaloAxis.html#e4b2ce253db4df8576ccc1f122a7e432">par</a>[1];
496     <a name="l00475"></a>00475 <a class="code" href="classCaloAxis.html#a1b225fccc498fdf49c415497d405eac">Reset</a>();
497     <a name="l00476"></a>00476
498     <a name="l00477"></a>00477 <span class="keywordtype">float</span> dmin[22];
499     <a name="l00478"></a>00478 <span class="keywordtype">int</span> closest[22];
500     <a name="l00479"></a>00479 <span class="keywordflow">for</span>(<span class="keywordtype">int</span> ip=0; ip&lt;22; ip++){
501     <a name="l00480"></a>00480 dmin[ip]=999;
502     <a name="l00481"></a>00481 closest[ip]=-1;
503     <a name="l00482"></a>00482 }
504     <a name="l00483"></a>00483 <span class="keywordflow">for</span>(<span class="keywordtype">int</span> ih=0;ih&lt;<a class="code" href="classCaloAxis.html#fb48c08633030431b7e509ff1b34558d">cevent</a>-&gt;<a class="code" href="classCaloEvent.html#2b2597fe0263fe128101a6ab02520130">GetN</a>(); ih++){
505     <a name="l00484"></a>00484 <span class="keywordtype">float</span> <a class="code" href="classCaloAxis.html#bea775e33c9f4ea5c308eace6ac4a513">x</a> = <a class="code" href="classCaloAxis.html#fb48c08633030431b7e509ff1b34558d">cevent</a>-&gt;<a class="code" href="classCaloEvent.html#b48403ffdb4114d57e4ea4058438e814">xycoord</a>[ih];
506     <a name="l00485"></a>00485 <span class="keywordtype">float</span> z = <a class="code" href="classCaloAxis.html#fb48c08633030431b7e509ff1b34558d">cevent</a>-&gt;<a class="code" href="classCaloEvent.html#340ff2ae2167b37d0b6acc56b1475daf">zcoord</a>[ih];
507     <a name="l00486"></a>00486 <span class="comment">//float q = cevent-&gt;q[ih];</span>
508     <a name="l00487"></a>00487 <span class="keywordtype">int</span> ip = <a class="code" href="classCaloAxis.html#fb48c08633030431b7e509ff1b34558d">cevent</a>-&gt;<a class="code" href="classCaloEvent.html#3ad199ffd33f6300bff64e854c6ad652">idp</a>[ih];
509     <a name="l00488"></a>00488 <span class="comment">//int is = cevent-&gt;ids[ih];</span>
510     <a name="l00489"></a>00489 <span class="keywordtype">float</span> d = fabs(x-P0-P1*z)/sqrt(P1*P1+1);
511     <a name="l00490"></a>00490 <span class="keywordflow">if</span>( d &lt; toll &amp;&amp; d&lt;dmin[ip] &amp;&amp; pok[ip]==1 ){
512     <a name="l00491"></a>00491 <span class="comment">// Add(z,x,q,q);</span>
513     <a name="l00492"></a>00492 closest[ip]=ih;
514     <a name="l00493"></a>00493 <span class="comment">// cog[ip] += x*q;</span>
515     <a name="l00494"></a>00494 <span class="comment">// qpl[ip] += q;</span>
516     <a name="l00495"></a>00495 }
517     <a name="l00496"></a>00496 }
518     <a name="l00497"></a>00497 <span class="keywordflow">for</span>(Int_t ip=0; ip&lt;22; ip++){
519     <a name="l00498"></a>00498 <span class="keywordflow">if</span>(closest[ip]!=-1){
520     <a name="l00499"></a>00499 <span class="keywordtype">float</span> <a class="code" href="classCaloAxis.html#bea775e33c9f4ea5c308eace6ac4a513">x</a> = <a class="code" href="classCaloAxis.html#fb48c08633030431b7e509ff1b34558d">cevent</a>-&gt;<a class="code" href="classCaloEvent.html#b48403ffdb4114d57e4ea4058438e814">xycoord</a>[closest[ip]];
521     <a name="l00500"></a>00500 <span class="keywordtype">float</span> z = <a class="code" href="classCaloAxis.html#fb48c08633030431b7e509ff1b34558d">cevent</a>-&gt;<a class="code" href="classCaloEvent.html#340ff2ae2167b37d0b6acc56b1475daf">zcoord</a>[closest[ip]];
522     <a name="l00501"></a>00501 <span class="keywordtype">float</span> <a class="code" href="classCaloAxis.html#06da2809e33c9fa67e151d210fdab8fe">q</a> = <a class="code" href="classCaloAxis.html#fb48c08633030431b7e509ff1b34558d">cevent</a>-&gt;<a class="code" href="classCaloEvent.html#688d5d2803d4688d4f1ef72cbaf42dcd">q</a>[closest[ip]];
523     <a name="l00502"></a>00502 <span class="comment">//int is = cevent-&gt;ids[closest[ip]];</span>
524     <a name="l00503"></a>00503 <a class="code" href="classCaloAxis.html#4953d4ace9da1118672a451d66afd01b">Add</a>(z,x,q,q);
525     <a name="l00504"></a>00504 <a class="code" href="classCaloAxis.html#575a6d4f2ab14d0039902b0818395222">cog</a>[ip] += x*q;
526     <a name="l00505"></a>00505 <a class="code" href="classCaloAxis.html#2ff6400f0393e7ab05627a438c6eeca5">qpl</a>[ip] += q;
527     <a name="l00506"></a>00506 <span class="comment">// cout &lt;&lt; ip &lt;&lt; " -o- "&lt;&lt;is&lt;&lt;endl;</span>
528     <a name="l00507"></a>00507 }
529     <a name="l00508"></a>00508 }
530     <a name="l00509"></a>00509 <span class="comment">// -----------------</span>
531     <a name="l00510"></a>00510 <span class="comment">// add +/- one strip</span>
532     <a name="l00511"></a>00511 <span class="comment">// -----------------</span>
533     <a name="l00512"></a>00512 <span class="keywordflow">for</span>(<span class="keywordtype">int</span> ih=0;ih&lt;<a class="code" href="classCaloAxis.html#fb48c08633030431b7e509ff1b34558d">cevent</a>-&gt;<a class="code" href="classCaloEvent.html#2b2597fe0263fe128101a6ab02520130">GetN</a>(); ih++){
534     <a name="l00513"></a>00513 <span class="keywordtype">float</span> <a class="code" href="classCaloAxis.html#bea775e33c9f4ea5c308eace6ac4a513">x</a> = <a class="code" href="classCaloAxis.html#fb48c08633030431b7e509ff1b34558d">cevent</a>-&gt;<a class="code" href="classCaloEvent.html#b48403ffdb4114d57e4ea4058438e814">xycoord</a>[ih];
535     <a name="l00514"></a>00514 <span class="keywordtype">float</span> z = <a class="code" href="classCaloAxis.html#fb48c08633030431b7e509ff1b34558d">cevent</a>-&gt;<a class="code" href="classCaloEvent.html#340ff2ae2167b37d0b6acc56b1475daf">zcoord</a>[ih];
536     <a name="l00515"></a>00515 <span class="keywordtype">float</span> <a class="code" href="classCaloAxis.html#06da2809e33c9fa67e151d210fdab8fe">q</a> = <a class="code" href="classCaloAxis.html#fb48c08633030431b7e509ff1b34558d">cevent</a>-&gt;<a class="code" href="classCaloEvent.html#688d5d2803d4688d4f1ef72cbaf42dcd">q</a>[ih];
537     <a name="l00516"></a>00516 <span class="keywordtype">int</span> ip = <a class="code" href="classCaloAxis.html#fb48c08633030431b7e509ff1b34558d">cevent</a>-&gt;<a class="code" href="classCaloEvent.html#3ad199ffd33f6300bff64e854c6ad652">idp</a>[ih];
538     <a name="l00517"></a>00517 <span class="keywordtype">int</span> is = <a class="code" href="classCaloAxis.html#fb48c08633030431b7e509ff1b34558d">cevent</a>-&gt;<a class="code" href="classCaloEvent.html#593e99c95482f5e63e43eaa8c843bba2">ids</a>[ih];
539     <a name="l00518"></a>00518 <span class="keywordflow">if</span>( closest[ip]!=-1 ){
540     <a name="l00519"></a>00519 <span class="keywordtype">int</span> isc = <a class="code" href="classCaloAxis.html#fb48c08633030431b7e509ff1b34558d">cevent</a>-&gt;<a class="code" href="classCaloEvent.html#593e99c95482f5e63e43eaa8c843bba2">ids</a>[closest[ip]];
541     <a name="l00520"></a>00520 <span class="keywordflow">if</span>( is == isc+1 || is == isc-1 ){
542     <a name="l00521"></a>00521 <a class="code" href="classCaloAxis.html#4953d4ace9da1118672a451d66afd01b">Add</a>(z,x,q,q);
543     <a name="l00522"></a>00522 <a class="code" href="classCaloAxis.html#575a6d4f2ab14d0039902b0818395222">cog</a>[ip] += x*q;
544     <a name="l00523"></a>00523 <a class="code" href="classCaloAxis.html#2ff6400f0393e7ab05627a438c6eeca5">qpl</a>[ip] += q;
545     <a name="l00524"></a>00524 <span class="comment">// cout &lt;&lt; ip &lt;&lt; " -+- "&lt;&lt;is&lt;&lt;endl;</span>
546     <a name="l00525"></a>00525 }
547     <a name="l00526"></a>00526 }
548     <a name="l00527"></a>00527 }
549     <a name="l00528"></a>00528 <a class="code" href="classCaloAxis.html#58d8a61ebffbef5318b9539cd7364633">Fit</a>();
550     <a name="l00529"></a>00529
551     <a name="l00530"></a>00530 <span class="keywordflow">if</span>(<a class="code" href="classCaloAxis.html#be448c5906fb90329316b50796d2c721">GetN</a>()&lt;3)<span class="keywordflow">return</span> 0;
552     <a name="l00531"></a>00531 <span class="keywordflow">for</span>(<span class="keywordtype">int</span> ip=0; ip&lt;22; ip++){
553     <a name="l00532"></a>00532 <span class="keywordflow">if</span>(<a class="code" href="classCaloAxis.html#2ff6400f0393e7ab05627a438c6eeca5">qpl</a>[ip]!=0){
554     <a name="l00533"></a>00533 <a class="code" href="classCaloAxis.html#575a6d4f2ab14d0039902b0818395222">cog</a>[ip]=<a class="code" href="classCaloAxis.html#575a6d4f2ab14d0039902b0818395222">cog</a>[ip]/<a class="code" href="classCaloAxis.html#2ff6400f0393e7ab05627a438c6eeca5">qpl</a>[ip];
555     <a name="l00534"></a>00534 <span class="comment">// Add(z,cog[ip],cog[ip],0.7);</span>
556     <a name="l00535"></a>00535 }
557     <a name="l00536"></a>00536 }
558     <a name="l00537"></a>00537
559     <a name="l00538"></a>00538 <span class="comment">// cout &lt;&lt; "(3)"&lt;&lt;endl;</span>
560     <a name="l00539"></a>00539 <span class="comment">// cout &lt;&lt; " n. hits :"&lt;&lt;GetN()&lt;&lt;endl;</span>
561     <a name="l00540"></a>00540 <span class="comment">// cout &lt;&lt; " P0 P1 :"&lt;&lt;par[0]&lt;&lt; " "&lt;&lt;par[1]&lt;&lt;endl;</span>
562     <a name="l00541"></a>00541 <span class="comment">// cout &lt;&lt; " Chi2 :"&lt;&lt;GetChi2()&lt;&lt;endl;</span>
563     <a name="l00542"></a>00542 <span class="comment">// cout &lt;&lt; " Q axis :"&lt;&lt;GetQaxis()&lt;&lt;endl;</span>
564     <a name="l00543"></a>00543 <span class="comment">// cout &lt;&lt; " Q strip :"&lt;&lt;GetQstrip()&lt;&lt;endl;</span>
565     <a name="l00544"></a>00544 <span class="comment">// cout &lt;&lt; " ---------------------------"&lt;&lt;endl;</span>
566     <a name="l00545"></a>00545 <span class="keywordflow">return</span> 1;
567     <a name="l00546"></a>00546 };
568     <a name="l00547"></a>00547
569     <a name="l00548"></a>00548
570     <a name="l00549"></a><a class="code" href="classCaloAxis.html#02b1f08ab710247112097a662a24cc5a">00549</a> <span class="keywordtype">int</span> <a class="code" href="classCaloAxis.html#02b1f08ab710247112097a662a24cc5a">FitAxis</a>( CaloLevel1* calo , Int_t view , Float_t toll ){ <span class="keywordflow">return</span> <a class="code" href="classCaloAxis.html#02b1f08ab710247112097a662a24cc5a">FitAxis</a>(calo,view,toll,0); };
571     <a name="l00550"></a>00550
572     <a name="l00551"></a>00551
573     <a name="l00552"></a>00552
574     <a name="l00553"></a>00553 <span class="comment">//-----------------------------------------------------------------</span>
575     <a name="l00554"></a>00554 <span class="comment">// fit the axis (optimized for interacting patterns?? ...forse)</span>
576     <a name="l00555"></a>00555 <span class="comment">//-----------------------------------------------------------------</span>
577     <a name="l00556"></a><a class="code" href="classCaloAxis.html#bf6b255e057ff7e1f49970ecec2145ad">00556</a> <span class="keywordtype">int</span> <a class="code" href="classCaloAxis.html#bf6b255e057ff7e1f49970ecec2145ad">FitShower</a>( CaloLevel1* calo , Int_t view , Float_t toll , Bool_t usemechanicalalignement ){
578     <a name="l00557"></a>00557
579     <a name="l00558"></a>00558 <span class="comment">// cout &lt;&lt; "CaloAxis::FitAxis(...)"&lt;&lt;endl;</span>
580     <a name="l00559"></a>00559
581     <a name="l00560"></a>00560 <a class="code" href="classCaloAxis.html#00e32ad98be21c697b8a0671b90e8dd1">Set</a>(calo,view,usemechanicalalignement);
582     <a name="l00561"></a>00561
583     <a name="l00562"></a>00562 <span class="comment">// ------------------------------</span>
584     <a name="l00563"></a>00563 <span class="comment">// first :</span>
585     <a name="l00564"></a>00564 <span class="comment">// ------------------------------</span>
586     <a name="l00565"></a>00565 <span class="keywordflow">for</span>(<span class="keywordtype">int</span> ih=0;ih&lt;<a class="code" href="classCaloAxis.html#fb48c08633030431b7e509ff1b34558d">cevent</a>-&gt;<a class="code" href="classCaloEvent.html#2b2597fe0263fe128101a6ab02520130">GetN</a>(); ih++){
587     <a name="l00566"></a>00566 <span class="keywordtype">float</span> <a class="code" href="classCaloAxis.html#bea775e33c9f4ea5c308eace6ac4a513">x</a> = <a class="code" href="classCaloAxis.html#fb48c08633030431b7e509ff1b34558d">cevent</a>-&gt;<a class="code" href="classCaloEvent.html#b48403ffdb4114d57e4ea4058438e814">xycoord</a>[ih];
588     <a name="l00567"></a>00567 <span class="comment">//float z = cevent-&gt;zcoord[ih];</span>
589     <a name="l00568"></a>00568 <span class="keywordtype">float</span> <a class="code" href="classCaloAxis.html#06da2809e33c9fa67e151d210fdab8fe">q</a> = <a class="code" href="classCaloAxis.html#fb48c08633030431b7e509ff1b34558d">cevent</a>-&gt;<a class="code" href="classCaloEvent.html#688d5d2803d4688d4f1ef72cbaf42dcd">q</a>[ih];
590     <a name="l00569"></a>00569 <span class="keywordtype">int</span> ip = <a class="code" href="classCaloAxis.html#fb48c08633030431b7e509ff1b34558d">cevent</a>-&gt;<a class="code" href="classCaloEvent.html#3ad199ffd33f6300bff64e854c6ad652">idp</a>[ih];
591     <a name="l00570"></a>00570 <a class="code" href="classCaloAxis.html#575a6d4f2ab14d0039902b0818395222">cog</a>[ip] += x*q;
592     <a name="l00571"></a>00571 <a class="code" href="classCaloAxis.html#2ff6400f0393e7ab05627a438c6eeca5">qpl</a>[ip] += q;
593     <a name="l00572"></a>00572 <span class="comment">// Add(z,x,q);</span>
594     <a name="l00573"></a>00573 }
595     <a name="l00574"></a>00574 <span class="keywordflow">for</span>(<span class="keywordtype">int</span> ip=0; ip&lt;22; ip++){
596     <a name="l00575"></a>00575 <span class="keywordflow">if</span>(<a class="code" href="classCaloAxis.html#2ff6400f0393e7ab05627a438c6eeca5">qpl</a>[ip]!=0)<a class="code" href="classCaloAxis.html#575a6d4f2ab14d0039902b0818395222">cog</a>[ip]=<a class="code" href="classCaloAxis.html#575a6d4f2ab14d0039902b0818395222">cog</a>[ip]/<a class="code" href="classCaloAxis.html#2ff6400f0393e7ab05627a438c6eeca5">qpl</a>[ip];
597     <a name="l00576"></a>00576 <span class="comment">// cout &lt;&lt; ip &lt;&lt; " &gt; "&lt;&lt;qpl[ip]&lt;&lt;" "&lt;&lt;cog[ip]&lt;&lt;endl;</span>
598     <a name="l00577"></a>00577 }
599     <a name="l00578"></a>00578 <span class="comment">// ----------------</span>
600     <a name="l00579"></a>00579 <span class="comment">// plane of maximum</span>
601     <a name="l00580"></a>00580 <span class="comment">// ----------------</span>
602     <a name="l00581"></a>00581 <span class="keywordtype">float</span> qplmax=1;
603     <a name="l00582"></a>00582 <span class="keywordtype">int</span> ipmax = -1;
604     <a name="l00583"></a>00583 <span class="keywordtype">float</span> dmin[22];
605     <a name="l00584"></a>00584 <span class="keywordtype">int</span> closest[22];
606     <a name="l00585"></a>00585 <span class="keywordflow">for</span>(<span class="keywordtype">int</span> ip=0; ip&lt;22; ip++)if(qpl[ip]&gt;qplmax){
607     <a name="l00586"></a>00586 ipmax=ip;
608     <a name="l00587"></a>00587 qplmax=<a class="code" href="classCaloAxis.html#2ff6400f0393e7ab05627a438c6eeca5">qpl</a>[ip];
609     <a name="l00588"></a>00588 dmin[ip]=1000.;<span class="comment">//init</span>
610     <a name="l00589"></a>00589 closest[ip]=-1;<span class="comment">//init</span>
611     <a name="l00590"></a>00590 }
612     <a name="l00591"></a>00591 <span class="keywordflow">for</span>(<span class="keywordtype">int</span> ih=0;ih&lt;<a class="code" href="classCaloAxis.html#fb48c08633030431b7e509ff1b34558d">cevent</a>-&gt;<a class="code" href="classCaloEvent.html#2b2597fe0263fe128101a6ab02520130">GetN</a>(); ih++){
613     <a name="l00592"></a>00592 <span class="keywordtype">float</span> <a class="code" href="classCaloAxis.html#bea775e33c9f4ea5c308eace6ac4a513">x</a> = <a class="code" href="classCaloAxis.html#fb48c08633030431b7e509ff1b34558d">cevent</a>-&gt;<a class="code" href="classCaloEvent.html#b48403ffdb4114d57e4ea4058438e814">xycoord</a>[ih];
614     <a name="l00593"></a>00593 <span class="keywordtype">float</span> z = <a class="code" href="classCaloAxis.html#fb48c08633030431b7e509ff1b34558d">cevent</a>-&gt;<a class="code" href="classCaloEvent.html#340ff2ae2167b37d0b6acc56b1475daf">zcoord</a>[ih];
615     <a name="l00594"></a>00594 <span class="keywordtype">float</span> <a class="code" href="classCaloAxis.html#06da2809e33c9fa67e151d210fdab8fe">q</a> = <a class="code" href="classCaloAxis.html#fb48c08633030431b7e509ff1b34558d">cevent</a>-&gt;<a class="code" href="classCaloEvent.html#688d5d2803d4688d4f1ef72cbaf42dcd">q</a>[ih];
616     <a name="l00595"></a>00595 <span class="keywordtype">int</span> ip = <a class="code" href="classCaloAxis.html#fb48c08633030431b7e509ff1b34558d">cevent</a>-&gt;<a class="code" href="classCaloEvent.html#3ad199ffd33f6300bff64e854c6ad652">idp</a>[ih];
617     <a name="l00596"></a>00596 <span class="comment">// if( (ipmax&gt;3 &amp;&amp; ip&lt;=ipmax) || ipmax&lt;=3 )Add(z,x,q);</span>
618     <a name="l00597"></a>00597 <span class="keywordflow">if</span>( ip&lt;=ipmax || ip&lt;=3 ){
619     <a name="l00598"></a>00598 <a class="code" href="classCaloAxis.html#4953d4ace9da1118672a451d66afd01b">Add</a>(z,x,q);
620     <a name="l00599"></a>00599 }
621     <a name="l00600"></a>00600 }
622     <a name="l00601"></a>00601 <span class="keywordflow">if</span>(<a class="code" href="classCaloAxis.html#be448c5906fb90329316b50796d2c721">GetN</a>()&lt;3)<span class="keywordflow">return</span> 0;
623     <a name="l00602"></a>00602 <a class="code" href="classCaloAxis.html#58d8a61ebffbef5318b9539cd7364633">Fit</a>();
624     <a name="l00603"></a>00603 <span class="comment">// cout &lt;&lt; " n. hits :"&lt;&lt;GetN()&lt;&lt;endl;</span>
625     <a name="l00604"></a>00604 <span class="comment">// cout &lt;&lt; " P0 P1 :"&lt;&lt;par[0]&lt;&lt; " "&lt;&lt;par[1]&lt;&lt;endl;</span>
626     <a name="l00605"></a>00605 <span class="comment">// cout &lt;&lt; " Chi2 :"&lt;&lt;GetChi2()&lt;&lt;endl;</span>
627     <a name="l00606"></a>00606 <span class="comment">// cout &lt;&lt; " Q axis :"&lt;&lt;GetQaxis()&lt;&lt;endl;</span>
628     <a name="l00607"></a>00607 <span class="comment">// cout &lt;&lt; " Q strip :"&lt;&lt;GetQstrip()&lt;&lt;endl;</span>
629     <a name="l00608"></a>00608
630     <a name="l00609"></a>00609 <span class="keywordtype">float</span> P0;
631     <a name="l00610"></a>00610 <span class="keywordtype">float</span> P1;
632     <a name="l00611"></a>00611 <span class="keywordtype">int</span> dof;
633     <a name="l00612"></a>00612
634     <a name="l00613"></a>00613 <span class="comment">// ---------------------------</span>
635     <a name="l00614"></a>00614 <span class="comment">// second : iteration</span>
636     <a name="l00615"></a>00615 <span class="comment">// ---------------------------</span>
637     <a name="l00616"></a>00616 <span class="keywordtype">int</span> niter=0; <span class="comment">// number of iterations</span>
638     <a name="l00617"></a>00617 <span class="keywordtype">int</span> pok[22];
639     <a name="l00618"></a>00618 <span class="comment">// prova</span>
640     <a name="l00619"></a>00619 <span class="comment">// pesi lineari</span>
641     <a name="l00620"></a>00620 <span class="keywordtype">float</span> W22=1.;
642     <a name="l00621"></a>00621 <span class="keywordtype">float</span> W1=10.;
643     <a name="l00622"></a>00622 <span class="keywordtype">float</span> b=(W22-W1)/(22.-1.);<span class="comment">// (w22 - w1)/(22-1)</span>
644     <a name="l00623"></a>00623 <span class="keywordtype">float</span> a=W1-b;
645     <a name="l00624"></a>00624
646     <a name="l00625"></a>00625 <span class="keywordflow">for</span>(<span class="keywordtype">int</span> i=0; i&lt;22; i++){
647     <a name="l00626"></a>00626 pok[i]=0;
648     <a name="l00627"></a>00627 }
649     <a name="l00628"></a>00628 <span class="keywordflow">do</span>{
650     <a name="l00629"></a>00629 niter++;
651     <a name="l00630"></a>00630 <span class="keywordtype">float</span> ttoll = toll; <span class="comment">//tolerance (cm)</span>
652     <a name="l00631"></a>00631 <span class="keywordflow">if</span>(niter==1)ttoll = 10.*toll;
653     <a name="l00632"></a>00632 dof = <a class="code" href="classCaloAxis.html#be448c5906fb90329316b50796d2c721">GetN</a>(); <span class="comment">//previous fit</span>
654     <a name="l00633"></a>00633 P0 = <a class="code" href="classCaloAxis.html#e4b2ce253db4df8576ccc1f122a7e432">par</a>[0];
655     <a name="l00634"></a>00634 P1 = <a class="code" href="classCaloAxis.html#e4b2ce253db4df8576ccc1f122a7e432">par</a>[1];
656     <a name="l00635"></a>00635 <a class="code" href="classCaloAxis.html#a1b225fccc498fdf49c415497d405eac">Reset</a>(); <span class="comment">//fit reset</span>
657     <a name="l00636"></a>00636 <span class="keywordflow">for</span>(<span class="keywordtype">int</span> i=0; i&lt;22; i++){
658     <a name="l00637"></a>00637 <a class="code" href="classCaloAxis.html#575a6d4f2ab14d0039902b0818395222">cog</a>[i]=0;
659     <a name="l00638"></a>00638 <a class="code" href="classCaloAxis.html#2ff6400f0393e7ab05627a438c6eeca5">qpl</a>[i]=0;
660     <a name="l00639"></a>00639 }
661     <a name="l00640"></a>00640 <span class="keywordflow">for</span>(<span class="keywordtype">int</span> ih=0;ih&lt;<a class="code" href="classCaloAxis.html#fb48c08633030431b7e509ff1b34558d">cevent</a>-&gt;<a class="code" href="classCaloEvent.html#2b2597fe0263fe128101a6ab02520130">GetN</a>(); ih++){
662     <a name="l00641"></a>00641 <span class="keywordtype">float</span> <a class="code" href="classCaloAxis.html#bea775e33c9f4ea5c308eace6ac4a513">x</a> = <a class="code" href="classCaloAxis.html#fb48c08633030431b7e509ff1b34558d">cevent</a>-&gt;<a class="code" href="classCaloEvent.html#b48403ffdb4114d57e4ea4058438e814">xycoord</a>[ih];
663     <a name="l00642"></a>00642 <span class="keywordtype">float</span> z = <a class="code" href="classCaloAxis.html#fb48c08633030431b7e509ff1b34558d">cevent</a>-&gt;<a class="code" href="classCaloEvent.html#340ff2ae2167b37d0b6acc56b1475daf">zcoord</a>[ih];
664     <a name="l00643"></a>00643 <span class="keywordtype">float</span> <a class="code" href="classCaloAxis.html#06da2809e33c9fa67e151d210fdab8fe">q</a> = <a class="code" href="classCaloAxis.html#fb48c08633030431b7e509ff1b34558d">cevent</a>-&gt;<a class="code" href="classCaloEvent.html#688d5d2803d4688d4f1ef72cbaf42dcd">q</a>[ih];
665     <a name="l00644"></a>00644 <span class="keywordtype">int</span> ip = <a class="code" href="classCaloAxis.html#fb48c08633030431b7e509ff1b34558d">cevent</a>-&gt;<a class="code" href="classCaloEvent.html#3ad199ffd33f6300bff64e854c6ad652">idp</a>[ih];
666     <a name="l00645"></a>00645 <span class="keywordtype">float</span> d = fabs(x-P0-P1*z)/sqrt(P1*P1+1);
667     <a name="l00646"></a>00646 <span class="comment">// cout&lt;&lt;d&lt;&lt;endl;</span>
668     <a name="l00647"></a>00647 <span class="keywordflow">if</span>(
669     <a name="l00648"></a>00648 (niter==1 || (niter&gt;1 &amp;&amp; pok[ip]==1)) &amp;&amp;
670     <a name="l00649"></a>00649 (ip&lt;=ipmax || ip&lt;=3) &amp;&amp;
671     <a name="l00650"></a>00650 <span class="comment">// ((ipmax&gt;3 &amp;&amp; ip&lt;=ipmax) || ipmax&lt;=3) &amp;&amp;</span>
672     <a name="l00651"></a>00651 <span class="keyword">true</span>){
673     <a name="l00652"></a>00652
674     <a name="l00653"></a>00653 <span class="keywordflow">if</span>(niter==1 &amp;&amp; d&lt;dmin[ip]){
675     <a name="l00654"></a>00654 dmin[ip]=d;
676     <a name="l00655"></a>00655 closest[ip]=ih;
677     <a name="l00656"></a>00656 }
678     <a name="l00657"></a>00657
679     <a name="l00658"></a>00658 <span class="keywordtype">float</span> W=a+b*(ip+1);
680     <a name="l00659"></a>00659 <span class="comment">// cout &lt;&lt; ip &lt;&lt; " "&lt;&lt;W&lt;&lt;endl;</span>
681     <a name="l00660"></a>00660 <span class="keywordflow">if</span>( d &lt; ttoll || ( niter==2 &amp;&amp; ih==closest[ip] )){
682     <a name="l00661"></a>00661 <span class="comment">/* if(ip==ipmax &amp;&amp; ipmax&gt;1)Add(z,x,q,W*q); */</span>
683     <a name="l00662"></a>00662 <span class="comment">/* else Add(z,x,q,W); */</span>
684     <a name="l00663"></a>00663 <a class="code" href="classCaloAxis.html#575a6d4f2ab14d0039902b0818395222">cog</a>[ip] += x*q;
685     <a name="l00664"></a>00664 <a class="code" href="classCaloAxis.html#2ff6400f0393e7ab05627a438c6eeca5">qpl</a>[ip] += q;
686     <a name="l00665"></a>00665 <a class="code" href="classCaloAxis.html#4953d4ace9da1118672a451d66afd01b">Add</a>(z,x,q,W);
687     <a name="l00666"></a>00666 pok[ip]=1;
688     <a name="l00667"></a>00667 cout &lt;&lt; ip &lt;&lt; <span class="stringliteral">" "</span>&lt;&lt;<a class="code" href="classCaloAxis.html#2ff6400f0393e7ab05627a438c6eeca5">qpl</a>[ip]&lt;&lt;endl;
689     <a name="l00668"></a>00668 }
690     <a name="l00669"></a>00669 }
691     <a name="l00670"></a>00670 }
692     <a name="l00671"></a>00671 <span class="comment">// break the track if more than 3 planes are missing</span>
693     <a name="l00672"></a>00672 <span class="keywordtype">int</span> nmissing = 0;
694     <a name="l00673"></a>00673 <span class="keywordflow">for</span>(<span class="keywordtype">int</span> ip=0; ip&lt;22; ip++){
695     <a name="l00674"></a>00674 <span class="keywordflow">if</span>(pok[ip]==0){
696     <a name="l00675"></a>00675 nmissing++;
697     <a name="l00676"></a>00676 }<span class="keywordflow">else</span>{
698     <a name="l00677"></a>00677 nmissing = 0;
699     <a name="l00678"></a>00678 }
700     <a name="l00679"></a>00679 <span class="keywordflow">if</span>(nmissing==6){
701     <a name="l00680"></a>00680 <span class="keywordflow">for</span>(<span class="keywordtype">int</span> ipp=ip+1; ipp&lt;22; ipp++)pok[ipp]=0;
702     <a name="l00681"></a>00681 <span class="keywordflow">break</span>;
703     <a name="l00682"></a>00682 }
704     <a name="l00683"></a>00683 }
705     <a name="l00684"></a>00684
706     <a name="l00685"></a>00685 <span class="keywordflow">for</span>(<span class="keywordtype">int</span> ip=0; ip&lt;22; ip++)cout &lt;&lt; <a class="code" href="classCaloAxis.html#2ff6400f0393e7ab05627a438c6eeca5">qpl</a>[ip] &lt;&lt; <span class="stringliteral">" "</span>;
707     <a name="l00686"></a>00686 cout &lt;&lt; endl;
708     <a name="l00687"></a>00687
709     <a name="l00688"></a>00688
710     <a name="l00689"></a>00689 <span class="keywordflow">for</span>(<span class="keywordtype">int</span> ip=0; ip&lt;22; ip++)<span class="keywordflow">if</span>(<a class="code" href="classCaloAxis.html#2ff6400f0393e7ab05627a438c6eeca5">qpl</a>[ip]!=0)<a class="code" href="classCaloAxis.html#575a6d4f2ab14d0039902b0818395222">cog</a>[ip]=<a class="code" href="classCaloAxis.html#575a6d4f2ab14d0039902b0818395222">cog</a>[ip]/<a class="code" href="classCaloAxis.html#2ff6400f0393e7ab05627a438c6eeca5">qpl</a>[ip];
711     <a name="l00690"></a>00690
712     <a name="l00691"></a>00691
713     <a name="l00692"></a>00692 <span class="keywordflow">if</span>(niter==100)<span class="keywordflow">break</span>;
714     <a name="l00693"></a>00693 <span class="keywordflow">if</span>(<a class="code" href="classCaloAxis.html#be448c5906fb90329316b50796d2c721">GetN</a>()&lt;3)<span class="keywordflow">break</span>;
715     <a name="l00694"></a>00694
716     <a name="l00695"></a>00695 <a class="code" href="classCaloAxis.html#58d8a61ebffbef5318b9539cd7364633">Fit</a>();
717     <a name="l00696"></a>00696
718     <a name="l00697"></a>00697 }<span class="keywordflow">while</span>(niter==1 || <a class="code" href="classCaloAxis.html#be448c5906fb90329316b50796d2c721">GetN</a>() != dof);
719     <a name="l00698"></a>00698 <span class="comment">// cout &lt;&lt; " &gt;&gt;&gt; "&lt;&lt;GetN()&lt;&lt;endl;</span>
720     <a name="l00699"></a>00699 <span class="keywordflow">if</span>(<a class="code" href="classCaloAxis.html#be448c5906fb90329316b50796d2c721">GetN</a>()&lt;3)<span class="keywordflow">return</span> 0;
721     <a name="l00700"></a>00700 <a class="code" href="classCaloAxis.html#58d8a61ebffbef5318b9539cd7364633">Fit</a>();
722     <a name="l00701"></a>00701 <span class="comment">/* // cout &lt;&lt; " n. hits :"&lt;&lt;GetN()&lt;&lt;endl; */</span>
723     <a name="l00702"></a>00702 <span class="comment">/* // cout &lt;&lt; " P0 P1 :"&lt;&lt;par[0]&lt;&lt; " "&lt;&lt;par[1]&lt;&lt;endl; */</span>
724     <a name="l00703"></a>00703 <span class="comment">/* // cout &lt;&lt; " Chi2 :"&lt;&lt;GetChi2()&lt;&lt;endl; */</span>
725     <a name="l00704"></a>00704 <span class="comment">/* // cout &lt;&lt; " Q axis :"&lt;&lt;GetQaxis()&lt;&lt;endl; */</span>
726     <a name="l00705"></a>00705 <span class="comment">/* // cout &lt;&lt; " Q strip :"&lt;&lt;GetQstrip()&lt;&lt;endl; */</span>
727     <a name="l00706"></a>00706
728     <a name="l00707"></a>00707 <span class="comment">/* // --------------------------------------------- */</span>
729     <a name="l00708"></a>00708 <span class="comment">/* // third : selecting only closest strip-clusters */</span>
730     <a name="l00709"></a>00709 <span class="comment">/* // and fit baricenters */</span>
731     <a name="l00710"></a>00710 <span class="comment">/* // --------------------------------------------- */</span>
732     <a name="l00711"></a>00711 <span class="comment">/* P0 = par[0]; */</span>
733     <a name="l00712"></a>00712 <span class="comment">/* P1 = par[1]; */</span>
734     <a name="l00713"></a>00713 <span class="comment">/* Reset(); */</span>
735     <a name="l00714"></a>00714
736     <a name="l00715"></a>00715 <span class="comment">/* float dmin[22]; */</span>
737     <a name="l00716"></a>00716 <span class="comment">/* int closest[22]; */</span>
738     <a name="l00717"></a>00717 <span class="comment">/* for(int ip=0; ip&lt;22; ip++){ */</span>
739     <a name="l00718"></a>00718 <span class="comment">/* dmin[ip]=999; */</span>
740     <a name="l00719"></a>00719 <span class="comment">/* closest[ip]=-1; */</span>
741     <a name="l00720"></a>00720 <span class="comment">/* } */</span>
742     <a name="l00721"></a>00721 <span class="comment">/* for(int ih=0;ih&lt;cevent-&gt;GetN(); ih++){ */</span>
743     <a name="l00722"></a>00722 <span class="comment">/* float x = cevent-&gt;xycoord[ih]; */</span>
744     <a name="l00723"></a>00723 <span class="comment">/* float z = cevent-&gt;zcoord[ih]; */</span>
745     <a name="l00724"></a>00724 <span class="comment">/* float q = cevent-&gt;q[ih]; */</span>
746     <a name="l00725"></a>00725 <span class="comment">/* int ip = cevent-&gt;idp[ih]; */</span>
747     <a name="l00726"></a>00726 <span class="comment">/* int is = cevent-&gt;ids[ih]; */</span>
748     <a name="l00727"></a>00727 <span class="comment">/* float d = fabs(x-P0-P1*z)/sqrt(P1*P1+1); */</span>
749     <a name="l00728"></a>00728 <span class="comment">/* if( d &lt; toll &amp;&amp; d&lt;dmin[ip] &amp;&amp; pok[ip]==1 ){ */</span>
750     <a name="l00729"></a>00729 <span class="comment">/* // Add(z,x,q,q); */</span>
751     <a name="l00730"></a>00730 <span class="comment">/* closest[ip]=ih; */</span>
752     <a name="l00731"></a>00731 <span class="comment">/* // cog[ip] += x*q; */</span>
753     <a name="l00732"></a>00732 <span class="comment">/* // qpl[ip] += q; */</span>
754     <a name="l00733"></a>00733 <span class="comment">/* } */</span>
755     <a name="l00734"></a>00734 <span class="comment">/* } */</span>
756     <a name="l00735"></a>00735 <span class="comment">/* for(Int_t ip=0; ip&lt;22; ip++){ */</span>
757     <a name="l00736"></a>00736 <span class="comment">/* if(closest[ip]!=-1){ */</span>
758     <a name="l00737"></a>00737 <span class="comment">/* float x = cevent-&gt;xycoord[closest[ip]]; */</span>
759     <a name="l00738"></a>00738 <span class="comment">/* float z = cevent-&gt;zcoord[closest[ip]]; */</span>
760     <a name="l00739"></a>00739 <span class="comment">/* float q = cevent-&gt;q[closest[ip]]; */</span>
761     <a name="l00740"></a>00740 <span class="comment">/* int is = cevent-&gt;ids[closest[ip]]; */</span>
762     <a name="l00741"></a>00741 <span class="comment">/* Add(z,x,q,q); */</span>
763     <a name="l00742"></a>00742 <span class="comment">/* cog[ip] += x*q; */</span>
764     <a name="l00743"></a>00743 <span class="comment">/* qpl[ip] += q; */</span>
765     <a name="l00744"></a>00744 <span class="comment">/* // cout &lt;&lt; ip &lt;&lt; " -o- "&lt;&lt;is&lt;&lt;endl; */</span>
766     <a name="l00745"></a>00745 <span class="comment">/* } */</span>
767     <a name="l00746"></a>00746 <span class="comment">/* } */</span>
768     <a name="l00747"></a>00747 <span class="comment">/* // add +/- one strip */</span>
769     <a name="l00748"></a>00748 <span class="comment">/* for(int ih=0;ih&lt;cevent-&gt;GetN(); ih++){ */</span>
770     <a name="l00749"></a>00749 <span class="comment">/* float x = cevent-&gt;xycoord[ih]; */</span>
771     <a name="l00750"></a>00750 <span class="comment">/* float z = cevent-&gt;zcoord[ih]; */</span>
772     <a name="l00751"></a>00751 <span class="comment">/* float q = cevent-&gt;q[ih]; */</span>
773     <a name="l00752"></a>00752 <span class="comment">/* int ip = cevent-&gt;idp[ih]; */</span>
774     <a name="l00753"></a>00753 <span class="comment">/* int is = cevent-&gt;ids[ih]; */</span>
775     <a name="l00754"></a>00754 <span class="comment">/* if( closest[ip]!=-1 ){ */</span>
776     <a name="l00755"></a>00755 <span class="comment">/* int isc = cevent-&gt;ids[closest[ip]]; */</span>
777     <a name="l00756"></a>00756 <span class="comment">/* if( is == isc+1 || is == isc-1 ){ */</span>
778     <a name="l00757"></a>00757 <span class="comment">/* Add(z,x,q,q); */</span>
779     <a name="l00758"></a>00758 <span class="comment">/* cog[ip] += x*q; */</span>
780     <a name="l00759"></a>00759 <span class="comment">/* qpl[ip] += q; */</span>
781     <a name="l00760"></a>00760 <span class="comment">/* // cout &lt;&lt; ip &lt;&lt; " -+- "&lt;&lt;is&lt;&lt;endl; */</span>
782     <a name="l00761"></a>00761 <span class="comment">/* } */</span>
783     <a name="l00762"></a>00762 <span class="comment">/* } */</span>
784     <a name="l00763"></a>00763 <span class="comment">/* } */</span>
785     <a name="l00764"></a>00764 <span class="comment">/* Fit(); */</span>
786     <a name="l00765"></a>00765
787     <a name="l00766"></a>00766 <span class="comment">/* if(GetN()&lt;3)return 0; */</span>
788     <a name="l00767"></a>00767 <span class="comment">/* for(int ip=0; ip&lt;22; ip++){ */</span>
789     <a name="l00768"></a>00768 <span class="comment">/* if(qpl[ip]!=0){ */</span>
790     <a name="l00769"></a>00769 <span class="comment">/* cog[ip]=cog[ip]/qpl[ip]; */</span>
791     <a name="l00770"></a>00770 <span class="comment">/* // Add(z,cog[ip],cog[ip],0.7); */</span>
792     <a name="l00771"></a>00771 <span class="comment">/* } */</span>
793     <a name="l00772"></a>00772 <span class="comment">/* } */</span>
794     <a name="l00773"></a>00773
795     <a name="l00774"></a>00774 <span class="comment">// cout &lt;&lt; " n. hits :"&lt;&lt;GetN()&lt;&lt;endl;</span>
796     <a name="l00775"></a>00775 <span class="comment">// cout &lt;&lt; " P0 P1 :"&lt;&lt;par[0]&lt;&lt; " "&lt;&lt;par[1]&lt;&lt;endl;</span>
797     <a name="l00776"></a>00776 <span class="comment">// cout &lt;&lt; " Chi2 :"&lt;&lt;GetChi2()&lt;&lt;endl;</span>
798     <a name="l00777"></a>00777 <span class="comment">// cout &lt;&lt; " Q axis :"&lt;&lt;GetQaxis()&lt;&lt;endl;</span>
799     <a name="l00778"></a>00778 <span class="comment">// cout &lt;&lt; " Q strip :"&lt;&lt;GetQstrip()&lt;&lt;endl;</span>
800     <a name="l00779"></a>00779 <span class="comment">// cout &lt;&lt; " ---------------------------"&lt;&lt;endl;</span>
801     <a name="l00780"></a>00780 <span class="keywordflow">return</span> 1;
802     <a name="l00781"></a>00781 };
803     <a name="l00782"></a>00782
804     <a name="l00783"></a>00783
805     <a name="l00784"></a><a class="code" href="classCaloAxis.html#eb0e8858c1297b4b3e3f5f66f09f508e">00784</a> <span class="keywordtype">int</span> <a class="code" href="classCaloAxis.html#eb0e8858c1297b4b3e3f5f66f09f508e">FitShower</a>( CaloLevel1* calo , Int_t view , Float_t toll ){ <span class="keywordflow">return</span> <a class="code" href="classCaloAxis.html#768db65112c0f2ecc7b0561772004e0f">FitAxis</a>(calo,view,toll,0); };
806     <a name="l00785"></a>00785
807     <a name="l00786"></a>00786
808     <a name="l00787"></a>00787 };
809     <a name="l00788"></a>00788
810     <a name="l00789"></a>00789 <span class="comment">//===============================================================================</span>
811     <a name="l00790"></a>00790 <span class="comment">//</span>
812     <a name="l00791"></a>00791 <span class="comment">//</span>
813     <a name="l00792"></a>00792 <span class="comment">//</span>
814     <a name="l00793"></a>00793 <span class="comment">//</span>
815     <a name="l00794"></a>00794 <span class="comment">//===============================================================================</span>
816     <a name="l00795"></a>00795
817     <a name="l00796"></a>00796 <span class="preprocessor">#endif // CALOAXIS2_H_</span>
818     </pre></div></div>
819     <hr size="1"><address style="text-align: right;"><small>Generated on Wed May 27 15:21:27 2009 for PamCut by&nbsp;
820     <a href="http://www.doxygen.org/index.html">
821     <img src="doxygen.png" alt="doxygen" align="middle" border="0"></a> 1.5.9 </small></address>
822     </body>
823     </html>
  ViewVC Help
Powered by ViewVC 1.1.23