| 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 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 List</span></a></li> |
| 19 |
|
|
<li><a href="globals.html"><span>File Members</span></a></li> |
| 20 |
|
|
</ul> |
| 21 |
|
|
</div> |
| 22 |
pam-fi |
1.2 |
<h1>CaloAxis2.h</h1><div class="fragment"><pre class="fragment"><a name="l00001"></a>00001 <span class="preprocessor">#ifndef CALOAXIS2_H_</span> |
| 23 |
pam-fi |
1.1 |
<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 <TString.h></span> |
| 29 |
|
|
<a name="l00008"></a>00008 <span class="preprocessor">#include <TH1F.h></span> |
| 30 |
|
|
<a name="l00009"></a>00009 <span class="preprocessor">#include <TH2F.h></span> |
| 31 |
|
|
<a name="l00010"></a>00010 <span class="preprocessor">#include <TMath.h></span> |
| 32 |
|
|
<a name="l00011"></a>00011 <span class="preprocessor">#include <TLine.h></span> |
| 33 |
|
|
<a name="l00012"></a>00012 <span class="preprocessor">#include <TPolyMarker.h></span> |
| 34 |
|
|
<a name="l00013"></a>00013 <span class="preprocessor">#include <TSelector.h></span> |
| 35 |
|
|
<a name="l00014"></a>00014 <span class="preprocessor">#include <TFile.h></span> |
| 36 |
|
|
<a name="l00015"></a>00015 |
| 37 |
|
|
<a name="l00016"></a>00016 <span class="preprocessor">#include <stdlib.h></span> |
| 38 |
|
|
<a name="l00017"></a>00017 <span class="preprocessor">#include <iostream></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 <PamLevel2.h></span> |
| 42 |
|
|
<a name="l00021"></a>00021 <span class="comment">//#include <CaloEvent.h></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 |
pam-fi |
1.2 |
<a name="l00030"></a>00030 <span class="keyword">class </span>CaloEvent{ |
| 52 |
pam-fi |
1.1 |
<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 |
pam-fi |
1.2 |
<a name="l00037"></a>00037 vector<int> ids; <span class="comment">//strip 0-95</span> |
| 59 |
|
|
<a name="l00038"></a>00038 vector<int> idp; <span class="comment">//plane 0-21</span> |
| 60 |
|
|
<a name="l00039"></a>00039 vector<float> q; <span class="comment">// signal</span> |
| 61 |
|
|
<a name="l00040"></a>00040 vector<float> zcoord; |
| 62 |
|
|
<a name="l00041"></a>00041 vector<float> xycoord; |
| 63 |
pam-fi |
1.1 |
<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 |
pam-fi |
1.2 |
<a name="l00047"></a>00047 <span class="keywordtype">void</span> Add(<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 ids.push_back(iis); |
| 70 |
|
|
<a name="l00049"></a>00049 idp.push_back(iip); |
| 71 |
|
|
<a name="l00050"></a>00050 q.push_back(fq); |
| 72 |
|
|
<a name="l00051"></a>00051 zcoord.push_back(fz); |
| 73 |
|
|
<a name="l00052"></a>00052 xycoord.push_back(fxy); |
| 74 |
pam-fi |
1.1 |
<a name="l00053"></a>00053 }; |
| 75 |
|
|
<a name="l00054"></a>00054 |
| 76 |
|
|
<a name="l00055"></a>00055 |
| 77 |
|
|
<a name="l00056"></a>00056 |
| 78 |
pam-fi |
1.2 |
<a name="l00057"></a>00057 <span class="keywordtype">int</span> GetN(){<span class="keywordflow">return</span> (<span class="keywordtype">int</span>)(q.size());}; |
| 79 |
pam-fi |
1.1 |
<a name="l00058"></a>00058 |
| 80 |
pam-fi |
1.2 |
<a name="l00059"></a>00059 <span class="keywordtype">void</span> Reset(){ |
| 81 |
|
|
<a name="l00060"></a>00060 ids.clear(); |
| 82 |
|
|
<a name="l00061"></a>00061 idp.clear(); |
| 83 |
|
|
<a name="l00062"></a>00062 q.clear(); |
| 84 |
|
|
<a name="l00063"></a>00063 zcoord.clear(); |
| 85 |
|
|
<a name="l00064"></a>00064 xycoord.clear(); |
| 86 |
pam-fi |
1.1 |
<a name="l00065"></a>00065 }; |
| 87 |
|
|
<a name="l00066"></a>00066 |
| 88 |
pam-fi |
1.2 |
<a name="l00067"></a>00067 <span class="keywordtype">void</span> Delete(){ Reset(); }; |
| 89 |
pam-fi |
1.1 |
<a name="l00068"></a>00068 |
| 90 |
pam-fi |
1.2 |
<a name="l00069"></a>00069 ~CaloEvent(){ Delete(); }; |
| 91 |
pam-fi |
1.1 |
<a name="l00070"></a>00070 |
| 92 |
pam-fi |
1.2 |
<a name="l00071"></a>00071 CaloEvent(){ }; |
| 93 |
pam-fi |
1.1 |
<a name="l00072"></a>00072 |
| 94 |
pam-fi |
1.2 |
<a name="l00073"></a>00073 CaloEvent( CaloLevel1* calo, <span class="keywordtype">int</span> view, Bool_t usemechanicalalignement){ Set(calo,view,usemechanicalalignement); }; |
| 95 |
pam-fi |
1.1 |
<a name="l00074"></a>00074 |
| 96 |
pam-fi |
1.2 |
<a name="l00075"></a>00075 CaloEvent( CaloLevel1* calo, <span class="keywordtype">int</span> view){ Set(calo,view); }; |
| 97 |
pam-fi |
1.1 |
<a name="l00076"></a>00076 |
| 98 |
pam-fi |
1.2 |
<a name="l00077"></a>00077 <span class="keywordtype">void</span> Set(CaloLevel1* calo, <span class="keywordtype">int</span> view, Bool_t usemechanicalalignement){ |
| 99 |
pam-fi |
1.1 |
<a name="l00078"></a>00078 |
| 100 |
pam-fi |
1.2 |
<a name="l00079"></a>00079 Reset(); |
| 101 |
pam-fi |
1.1 |
<a name="l00080"></a>00080 <span class="comment">// cout << " CaloEvent::Set()"<<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<calo->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->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 |
pam-fi |
1.2 |
<a name="l00096"></a>00096 <span class="keywordflow">if</span>(fq>0)Add(is,ip,fq,fz,fxy); |
| 118 |
pam-fi |
1.1 |
<a name="l00097"></a>00097 } |
| 119 |
|
|
<a name="l00098"></a>00098 } |
| 120 |
|
|
<a name="l00099"></a>00099 |
| 121 |
|
|
<a name="l00100"></a>00100 }; |
| 122 |
pam-fi |
1.2 |
<a name="l00101"></a>00101 <span class="keywordtype">void</span> Set(CaloLevel1* calo, <span class="keywordtype">int</span> view){ Set(calo,view,0); }; |
| 123 |
pam-fi |
1.1 |
<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 |
pam-fi |
1.2 |
<a name="l00111"></a>00111 <span class="keyword">class </span>CaloAxis { |
| 133 |
pam-fi |
1.1 |
<a name="l00112"></a>00112 |
| 134 |
|
|
<a name="l00113"></a>00113 <span class="keyword">public</span>: |
| 135 |
|
|
<a name="l00114"></a>00114 |
| 136 |
pam-fi |
1.2 |
<a name="l00115"></a>00115 CaloEvent *cevent; |
| 137 |
pam-fi |
1.1 |
<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 |
pam-fi |
1.2 |
<a name="l00121"></a>00121 vector<float> x;<span class="comment">// independent coordinates (z)</span> |
| 143 |
|
|
<a name="l00122"></a>00122 vector<float> y;<span class="comment">// dependente coordinate (x/y)</span> |
| 144 |
|
|
<a name="l00123"></a>00123 vector<float> w;<span class="comment">// fit weight</span> |
| 145 |
|
|
<a name="l00124"></a>00124 vector<float> q;<span class="comment">// signal</span> |
| 146 |
pam-fi |
1.1 |
<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 |
pam-fi |
1.2 |
<a name="l00129"></a>00129 Float_t par[2]; |
| 151 |
|
|
<a name="l00130"></a>00130 Float_t covpar[2][2]; |
| 152 |
pam-fi |
1.1 |
<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 |
pam-fi |
1.2 |
<a name="l00135"></a>00135 <span class="keywordtype">float</span> cog[22];<span class="comment">// baricenter</span> |
| 157 |
|
|
<a name="l00136"></a>00136 <span class="keywordtype">float</span> qpl[22]; <span class="comment">// charge in each plane</span> |
| 158 |
pam-fi |
1.1 |
<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 |
pam-fi |
1.2 |
<a name="l00141"></a>00141 TPolyMarker *sel; |
| 163 |
|
|
<a name="l00142"></a>00142 TLine *axis; |
| 164 |
pam-fi |
1.1 |
<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 |
pam-fi |
1.2 |
<a name="l00148"></a>00148 <span class="keywordtype">void</span> Init(){ |
| 170 |
pam-fi |
1.1 |
<a name="l00149"></a>00149 <span class="keywordflow">for</span>(Int_t i=0; i<2;i++){ |
| 171 |
pam-fi |
1.2 |
<a name="l00150"></a>00150 par[i]=0; |
| 172 |
|
|
<a name="l00151"></a>00151 <span class="keywordflow">for</span>(Int_t j=0; j<2;j++)covpar[i][j]=0; |
| 173 |
pam-fi |
1.1 |
<a name="l00152"></a>00152 }; |
| 174 |
pam-fi |
1.2 |
<a name="l00153"></a>00153 sel = 0; |
| 175 |
|
|
<a name="l00154"></a>00154 axis = 0; |
| 176 |
pam-fi |
1.1 |
<a name="l00155"></a>00155 <span class="keywordflow">for</span>(Int_t i=0; i<22; i++){ |
| 177 |
pam-fi |
1.2 |
<a name="l00156"></a>00156 cog[i]=0; |
| 178 |
|
|
<a name="l00157"></a>00157 qpl[i]=0; |
| 179 |
pam-fi |
1.1 |
<a name="l00158"></a>00158 } |
| 180 |
|
|
<a name="l00159"></a>00159 <span class="comment">// if(cevent)cevent->Reset();</span> |
| 181 |
|
|
<a name="l00160"></a>00160 } |
| 182 |
|
|
<a name="l00161"></a>00161 |
| 183 |
pam-fi |
1.2 |
<a name="l00162"></a>00162 <span class="keywordtype">void</span> Set(CaloLevel1* calo, Int_t view){ cevent->Set(calo,view,0); }; |
| 184 |
|
|
<a name="l00163"></a>00163 <span class="keywordtype">void</span> Set(CaloLevel1* calo, Int_t view, Bool_t usemechanicalalignement){ cevent->Set(calo,view,usemechanicalalignement); }; |
| 185 |
pam-fi |
1.1 |
<a name="l00164"></a>00164 |
| 186 |
pam-fi |
1.2 |
<a name="l00165"></a>00165 CaloAxis(){ |
| 187 |
|
|
<a name="l00166"></a>00166 cevent = <span class="keyword">new</span> CaloEvent(); |
| 188 |
|
|
<a name="l00167"></a>00167 Init(); |
| 189 |
pam-fi |
1.1 |
<a name="l00168"></a>00168 }; |
| 190 |
|
|
<a name="l00169"></a>00169 |
| 191 |
pam-fi |
1.2 |
<a name="l00170"></a>00170 CaloAxis(CaloLevel1* calo, Int_t view){ |
| 192 |
|
|
<a name="l00171"></a>00171 cevent = <span class="keyword">new</span> CaloEvent(); |
| 193 |
|
|
<a name="l00172"></a>00172 Init(); |
| 194 |
|
|
<a name="l00173"></a>00173 Set(calo,view); |
| 195 |
pam-fi |
1.1 |
<a name="l00174"></a>00174 }; |
| 196 |
pam-fi |
1.2 |
<a name="l00175"></a>00175 CaloAxis(CaloLevel1* calo, Int_t view, Bool_t usemechanicalalignement){ |
| 197 |
|
|
<a name="l00176"></a>00176 cevent = <span class="keyword">new</span> CaloEvent(); |
| 198 |
|
|
<a name="l00177"></a>00177 Init(); |
| 199 |
|
|
<a name="l00178"></a>00178 Set(calo,view,usemechanicalalignement); |
| 200 |
pam-fi |
1.1 |
<a name="l00179"></a>00179 }; |
| 201 |
|
|
<a name="l00180"></a>00180 |
| 202 |
|
|
<a name="l00181"></a>00181 |
| 203 |
pam-fi |
1.2 |
<a name="l00182"></a>00182 <span class="keywordtype">void</span> Reset(){ |
| 204 |
pam-fi |
1.1 |
<a name="l00183"></a>00183 <span class="comment">// cout << " CaloAxis::Reset()"<<endl;</span> |
| 205 |
pam-fi |
1.2 |
<a name="l00184"></a>00184 x.clear(); |
| 206 |
|
|
<a name="l00185"></a>00185 y.clear(); |
| 207 |
|
|
<a name="l00186"></a>00186 w.clear(); |
| 208 |
|
|
<a name="l00187"></a>00187 q.clear(); |
| 209 |
pam-fi |
1.1 |
<a name="l00188"></a>00188 <span class="comment">// cout << sel << endl;</span> |
| 210 |
pam-fi |
1.2 |
<a name="l00189"></a>00189 <span class="keywordflow">if</span>(sel)sel->Delete(); |
| 211 |
pam-fi |
1.1 |
<a name="l00190"></a>00190 <span class="comment">// cout << axis << endl;</span> |
| 212 |
pam-fi |
1.2 |
<a name="l00191"></a>00191 <span class="keywordflow">if</span>(axis)axis->Delete(); |
| 213 |
pam-fi |
1.1 |
<a name="l00192"></a>00192 <span class="comment">// if(cevent)cevent->Reset();</span> |
| 214 |
pam-fi |
1.2 |
<a name="l00193"></a>00193 Init(); |
| 215 |
pam-fi |
1.1 |
<a name="l00194"></a>00194 } |
| 216 |
|
|
<a name="l00195"></a>00195 |
| 217 |
pam-fi |
1.2 |
<a name="l00196"></a>00196 <span class="keywordtype">void</span> Delete(){ Reset(); <span class="keyword">delete</span> cevent; }; |
| 218 |
pam-fi |
1.1 |
<a name="l00197"></a>00197 |
| 219 |
pam-fi |
1.2 |
<a name="l00198"></a>00198 ~CaloAxis(){ Delete(); }; |
| 220 |
pam-fi |
1.1 |
<a name="l00199"></a>00199 |
| 221 |
|
|
<a name="l00200"></a>00200 |
| 222 |
pam-fi |
1.2 |
<a name="l00201"></a>00201 <span class="keywordtype">void</span> Add(<span class="keywordtype">float</span> xin, <span class="keywordtype">float</span> yin){ |
| 223 |
|
|
<a name="l00202"></a>00202 x.push_back(xin); |
| 224 |
|
|
<a name="l00203"></a>00203 y.push_back(yin); |
| 225 |
|
|
<a name="l00204"></a>00204 w.push_back(1.); |
| 226 |
|
|
<a name="l00205"></a>00205 q.push_back(0.); |
| 227 |
pam-fi |
1.1 |
<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>2)Fit(); */</span> |
| 229 |
|
|
<a name="l00208"></a>00208 }; |
| 230 |
pam-fi |
1.2 |
<a name="l00209"></a>00209 <span class="keywordtype">void</span> Add(<span class="keywordtype">float</span> xin, <span class="keywordtype">float</span> yin, <span class="keywordtype">float</span> qin){ |
| 231 |
|
|
<a name="l00210"></a>00210 x.push_back(xin); |
| 232 |
|
|
<a name="l00211"></a>00211 y.push_back(yin); |
| 233 |
|
|
<a name="l00212"></a>00212 q.push_back(qin); |
| 234 |
|
|
<a name="l00213"></a>00213 w.push_back(1.); |
| 235 |
pam-fi |
1.1 |
<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>2)Fit(); */</span> |
| 237 |
|
|
<a name="l00216"></a>00216 }; |
| 238 |
pam-fi |
1.2 |
<a name="l00217"></a>00217 <span class="keywordtype">void</span> Add(<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 x.push_back(xin); |
| 240 |
|
|
<a name="l00219"></a>00219 y.push_back(yin); |
| 241 |
|
|
<a name="l00220"></a>00220 w.push_back(win); |
| 242 |
|
|
<a name="l00221"></a>00221 q.push_back(qin); |
| 243 |
pam-fi |
1.1 |
<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>2)Fit(); */</span> |
| 245 |
|
|
<a name="l00224"></a>00224 }; |
| 246 |
|
|
<a name="l00225"></a>00225 |
| 247 |
pam-fi |
1.2 |
<a name="l00226"></a>00226 <span class="keywordtype">int</span> Fit(){ |
| 248 |
pam-fi |
1.1 |
<a name="l00227"></a>00227 |
| 249 |
|
|
<a name="l00228"></a>00228 <span class="comment">// cout << " CaloAxis::Fit()"<<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 |
pam-fi |
1.2 |
<a name="l00237"></a>00237 <span class="keywordflow">for</span>(Int_t i=0; i<(int)(x.size()); i++){ |
| 259 |
|
|
<a name="l00238"></a>00238 SSS += w[i]*w[i]; |
| 260 |
|
|
<a name="l00239"></a>00239 SSX += x[i]*w[i]*w[i]; |
| 261 |
|
|
<a name="l00240"></a>00240 SXX += x[i]*x[i]*w[i]*w[i]; |
| 262 |
|
|
<a name="l00241"></a>00241 SSY += y[i]*w[i]*w[i]; |
| 263 |
|
|
<a name="l00242"></a>00242 SXY += x[i]*y[i]*w[i]*w[i]; |
| 264 |
pam-fi |
1.1 |
<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<3)<span class="keywordflow">return</span> 0; |
| 269 |
|
|
<a name="l00248"></a>00248 <span class="comment">// cout << np << " points fitted -- ok "<<endl;</span> |
| 270 |
pam-fi |
1.2 |
<a name="l00249"></a>00249 par[0] = (SSY*SXX-SXY*SSX)/det; |
| 271 |
|
|
<a name="l00250"></a>00250 par[1] = (SXY*SSS-SSY*SSX)/det; |
| 272 |
pam-fi |
1.1 |
<a name="l00251"></a>00251 <span class="keywordflow">return</span> 1; |
| 273 |
|
|
<a name="l00252"></a>00252 }; |
| 274 |
|
|
<a name="l00253"></a>00253 |
| 275 |
pam-fi |
1.2 |
<a name="l00254"></a>00254 <span class="keywordtype">void</span> Print(){ |
| 276 |
pam-fi |
1.1 |
<a name="l00255"></a>00255 <span class="comment">// useful for debug for(Int_t i=0; i<(int)(x.size()); i++)if(w[i])cout << x[i] << " - "<<y[i]<<endl;</span> |
| 277 |
|
|
<a name="l00256"></a>00256 } |
| 278 |
|
|
<a name="l00257"></a>00257 |
| 279 |
pam-fi |
1.2 |
<a name="l00258"></a>00258 <span class="keywordtype">float</span> GetPar(Int_t i){<span class="keywordflow">return</span> par[i];}; |
| 280 |
pam-fi |
1.1 |
<a name="l00259"></a>00259 |
| 281 |
pam-fi |
1.2 |
<a name="l00260"></a>00260 <span class="keywordtype">int</span> GetN(){<span class="keywordflow">return</span> (<span class="keywordtype">int</span>)(x.size());}; |
| 282 |
pam-fi |
1.1 |
<a name="l00261"></a>00261 |
| 283 |
pam-fi |
1.2 |
<a name="l00262"></a>00262 <span class="keywordtype">float</span>* GetX(){ |
| 284 |
|
|
<a name="l00263"></a>00263 <span class="keywordtype">float</span> *xout = <span class="keyword">new</span> <span class="keywordtype">float</span>[(int)(x.size())]; |
| 285 |
|
|
<a name="l00264"></a>00264 <span class="keywordflow">for</span>(Int_t i=0; i<(int)(x.size()); i++)xout[i]=x[i]; |
| 286 |
pam-fi |
1.1 |
<a name="l00265"></a>00265 <span class="keywordflow">return</span> xout; |
| 287 |
|
|
<a name="l00266"></a>00266 } |
| 288 |
|
|
<a name="l00267"></a>00267 |
| 289 |
pam-fi |
1.2 |
<a name="l00268"></a>00268 <span class="keywordtype">float</span>* GetY(){ |
| 290 |
|
|
<a name="l00269"></a>00269 <span class="keywordtype">float</span> *yout = <span class="keyword">new</span> <span class="keywordtype">float</span>[(int)(y.size())]; |
| 291 |
|
|
<a name="l00270"></a>00270 <span class="keywordflow">for</span>(Int_t i=0; i<(int)(y.size()); i++) yout[i]=y[i]; |
| 292 |
pam-fi |
1.1 |
<a name="l00271"></a>00271 <span class="keywordflow">return</span> yout; |
| 293 |
|
|
<a name="l00272"></a>00272 } |
| 294 |
|
|
<a name="l00273"></a>00273 |
| 295 |
pam-fi |
1.2 |
<a name="l00274"></a>00274 <span class="keywordtype">float</span>* GetQ(){ |
| 296 |
|
|
<a name="l00275"></a>00275 <span class="keywordtype">float</span> *qout = <span class="keyword">new</span> <span class="keywordtype">float</span>[(int)(q.size())]; |
| 297 |
|
|
<a name="l00276"></a>00276 <span class="keywordflow">for</span>(<span class="keywordtype">int</span> i=0; i<(int)(q.size()); i++) qout[i]=q[i]; |
| 298 |
pam-fi |
1.1 |
<a name="l00277"></a>00277 <span class="keywordflow">return</span> qout; |
| 299 |
|
|
<a name="l00278"></a>00278 } |
| 300 |
|
|
<a name="l00279"></a>00279 |
| 301 |
pam-fi |
1.2 |
<a name="l00280"></a>00280 <span class="keywordtype">float</span> GetChi2(){ |
| 302 |
pam-fi |
1.1 |
<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 |
pam-fi |
1.2 |
<a name="l00283"></a>00283 <span class="keywordflow">for</span>(Int_t i=0; i<(int)(y.size()); i++){ |
| 305 |
|
|
<a name="l00284"></a>00284 chi2 += w[i]*w[i]*(y[i]-par[0]-par[1]*x[i])*(y[i]-par[0]-par[1]*x[i]); |
| 306 |
pam-fi |
1.1 |
<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<3) <span class="keywordflow">return</span> -999; <span class="comment">// original</span> |
| 310 |
|
|
<a name="l00289"></a>00289 <span class="comment">// if(nchi<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 |
pam-fi |
1.2 |
<a name="l00295"></a>00295 <span class="keywordtype">float</span> GetQaxis(){ |
| 317 |
pam-fi |
1.1 |
<a name="l00296"></a>00296 <span class="keywordtype">float</span> qaxis=0; |
| 318 |
pam-fi |
1.2 |
<a name="l00297"></a>00297 <span class="keywordflow">for</span>(Int_t i=0; i<(int)(q.size()); i++) qaxis +=q[i]; |
| 319 |
pam-fi |
1.1 |
<a name="l00298"></a>00298 <span class="comment">// for(Int_t i=0; i<(int)(q.size()); i++) cout<<q[i]<<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 |
pam-fi |
1.2 |
<a name="l00302"></a>00302 <span class="keywordtype">float</span> GetQaxis(<span class="keywordtype">float</span> toll){ |
| 324 |
pam-fi |
1.1 |
<a name="l00303"></a>00303 <span class="keywordtype">float</span> qaxis=0; |
| 325 |
pam-fi |
1.2 |
<a name="l00304"></a>00304 <span class="keywordflow">for</span>(<span class="keywordtype">int</span> ih=0;ih<cevent->GetN(); ih++){ |
| 326 |
|
|
<a name="l00305"></a>00305 <span class="keywordtype">float</span> x = cevent->xycoord[ih]; |
| 327 |
|
|
<a name="l00306"></a>00306 <span class="keywordtype">float</span> z = cevent->zcoord[ih]; |
| 328 |
|
|
<a name="l00307"></a>00307 <span class="keywordtype">float</span> q = cevent->q[ih]; |
| 329 |
pam-fi |
1.1 |
<a name="l00308"></a>00308 <span class="comment">//int ip = cevent->idp[ih];</span> |
| 330 |
pam-fi |
1.2 |
<a name="l00309"></a>00309 <span class="keywordtype">float</span> d = fabs(x-par[0]-par[1]*z)/sqrt(par[1]*par[1]+1); |
| 331 |
pam-fi |
1.1 |
<a name="l00310"></a>00310 <span class="keywordflow">if</span>( d < 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 |
pam-fi |
1.2 |
<a name="l00315"></a>00315 <span class="keywordtype">float</span> GetCOG(Int_t ip){ |
| 337 |
|
|
<a name="l00316"></a>00316 <span class="keywordflow">if</span>(ip>=0 && ip<22)<span class="keywordflow">return</span> cog[ip]; |
| 338 |
pam-fi |
1.1 |
<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 |
pam-fi |
1.2 |
<a name="l00320"></a>00320 <span class="keywordtype">float</span> GetQ(Int_t ip){ |
| 342 |
|
|
<a name="l00321"></a>00321 <span class="keywordflow">if</span>(ip>=0 && ip<22)<span class="keywordflow">return</span> qpl[ip]; |
| 343 |
pam-fi |
1.1 |
<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 |
pam-fi |
1.2 |
<a name="l00325"></a>00325 <span class="keywordtype">float</span> GetQstrip(){ |
| 347 |
pam-fi |
1.1 |
<a name="l00326"></a>00326 |
| 348 |
|
|
<a name="l00327"></a>00327 <span class="keywordtype">float</span> qaxis=0; |
| 349 |
pam-fi |
1.2 |
<a name="l00328"></a>00328 <span class="keywordflow">for</span>(Int_t i=0; i<(int)(q.size()); i++) qaxis +=q[i]; |
| 350 |
|
|
<a name="l00329"></a>00329 <span class="keywordflow">if</span>(q.size()==0) <span class="keywordflow">return</span> 0.; |
| 351 |
|
|
<a name="l00330"></a>00330 <span class="keywordflow">return</span> qaxis/q.size(); |
| 352 |
pam-fi |
1.1 |
<a name="l00331"></a>00331 } |
| 353 |
|
|
<a name="l00332"></a>00332 |
| 354 |
pam-fi |
1.2 |
<a name="l00333"></a>00333 <span class="keywordtype">float</span> GetYfit(<span class="keywordtype">float</span> x){ |
| 355 |
|
|
<a name="l00334"></a>00334 <span class="keywordflow">return</span> par[0]+x*par[1]; |
| 356 |
pam-fi |
1.1 |
<a name="l00335"></a>00335 } |
| 357 |
|
|
<a name="l00336"></a>00336 |
| 358 |
pam-fi |
1.2 |
<a name="l00337"></a>00337 <span class="keywordtype">float</span> GetXmin(){ |
| 359 |
pam-fi |
1.1 |
<a name="l00338"></a>00338 <span class="keywordtype">float</span> zmin=9999; |
| 360 |
pam-fi |
1.2 |
<a name="l00339"></a>00339 <span class="keywordflow">for</span>(Int_t i=0; i<(int)(x.size()); i++) <span class="keywordflow">if</span>(x[i]<zmin) zmin = x[i]; |
| 361 |
pam-fi |
1.1 |
<a name="l00340"></a>00340 <span class="keywordflow">return</span> zmin; |
| 362 |
|
|
<a name="l00341"></a>00341 } |
| 363 |
pam-fi |
1.2 |
<a name="l00342"></a>00342 <span class="keywordtype">float</span> GetXmax(){ |
| 364 |
pam-fi |
1.1 |
<a name="l00343"></a>00343 <span class="keywordtype">float</span> zmax=-9999; |
| 365 |
pam-fi |
1.2 |
<a name="l00344"></a>00344 <span class="keywordflow">for</span>(Int_t i=0; i<(int)(x.size()); i++) <span class="keywordflow">if</span>(x[i]>zmax) zmax = x[i]; |
| 366 |
pam-fi |
1.1 |
<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 |
pam-fi |
1.2 |
<a name="l00350"></a>00350 <span class="keywordtype">void</span> Draw(){ |
| 372 |
pam-fi |
1.1 |
<a name="l00351"></a>00351 |
| 373 |
|
|
<a name="l00352"></a>00352 <span class="comment">// cout << " CaloAxis::Draw()"<<endl;</span> |
| 374 |
pam-fi |
1.2 |
<a name="l00353"></a>00353 <span class="keywordflow">if</span>(GetN()<3)<span class="keywordflow">return</span>; |
| 375 |
pam-fi |
1.1 |
<a name="l00354"></a>00354 |
| 376 |
pam-fi |
1.2 |
<a name="l00355"></a>00355 sel = <span class="keyword">new</span> TPolyMarker(GetN(),GetY(),GetX()); |
| 377 |
|
|
<a name="l00356"></a>00356 sel->SetMarkerSize(0.3); |
| 378 |
|
|
<a name="l00357"></a>00357 sel->SetMarkerColor(5); |
| 379 |
pam-fi |
1.1 |
<a name="l00358"></a>00358 |
| 380 |
pam-fi |
1.2 |
<a name="l00359"></a>00359 axis = <span class="keyword">new</span> TLine( par[0]+GetXmin()*par[1], GetXmin(), par[0]+GetXmax()*par[1], GetXmax() ); |
| 381 |
|
|
<a name="l00360"></a>00360 axis->SetLineWidth(1); |
| 382 |
|
|
<a name="l00361"></a>00361 axis->SetLineColor(5); |
| 383 |
pam-fi |
1.1 |
<a name="l00362"></a>00362 |
| 384 |
|
|
<a name="l00363"></a>00363 <span class="comment">// cout << sel << endl;</span> |
| 385 |
|
|
<a name="l00364"></a>00364 <span class="comment">// cout << axis << endl;</span> |
| 386 |
|
|
<a name="l00365"></a>00365 |
| 387 |
pam-fi |
1.2 |
<a name="l00366"></a>00366 sel->Draw(); |
| 388 |
|
|
<a name="l00367"></a>00367 axis->Draw(); |
| 389 |
pam-fi |
1.1 |
<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 |
pam-fi |
1.2 |
<a name="l00374"></a>00374 <span class="keywordtype">int</span> FitAxis( CaloLevel1* calo , Int_t view , Float_t toll , Bool_t usemechanicalalignement ){ |
| 396 |
pam-fi |
1.1 |
<a name="l00375"></a>00375 |
| 397 |
|
|
<a name="l00376"></a>00376 <span class="comment">// cout << "CaloAxis::FitAxis(...)"<<endl;</span> |
| 398 |
|
|
<a name="l00377"></a>00377 |
| 399 |
pam-fi |
1.2 |
<a name="l00378"></a>00378 Set(calo,view,usemechanicalalignement); |
| 400 |
pam-fi |
1.1 |
<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 |
pam-fi |
1.2 |
<a name="l00383"></a>00383 <span class="keywordflow">for</span>(<span class="keywordtype">int</span> ih=0;ih<cevent->GetN(); ih++){ |
| 405 |
|
|
<a name="l00384"></a>00384 <span class="keywordtype">float</span> x = cevent->xycoord[ih]; |
| 406 |
|
|
<a name="l00385"></a>00385 <span class="keywordtype">float</span> z = cevent->zcoord[ih]; |
| 407 |
|
|
<a name="l00386"></a>00386 <span class="keywordtype">float</span> q = cevent->q[ih]; |
| 408 |
|
|
<a name="l00387"></a>00387 Add(z,x,q); |
| 409 |
pam-fi |
1.1 |
<a name="l00388"></a>00388 } |
| 410 |
|
|
<a name="l00389"></a>00389 |
| 411 |
pam-fi |
1.2 |
<a name="l00390"></a>00390 <span class="keywordflow">if</span>(GetN()<3)<span class="keywordflow">return</span> 0; |
| 412 |
|
|
<a name="l00391"></a>00391 Fit(); |
| 413 |
pam-fi |
1.1 |
<a name="l00392"></a>00392 <span class="comment">// useful for debug cout << "(1)"<<endl;</span> |
| 414 |
|
|
<a name="l00393"></a>00393 <span class="comment">// cout << " n. hits :"<<GetN()<<endl;</span> |
| 415 |
|
|
<a name="l00394"></a>00394 <span class="comment">// cout << " P0 P1 :"<<par[0]<< " "<<par[1]<<endl;</span> |
| 416 |
|
|
<a name="l00395"></a>00395 <span class="comment">// cout << " Chi2 :"<<GetChi2()<<endl;</span> |
| 417 |
|
|
<a name="l00396"></a>00396 <span class="comment">// cout << " Q axis :"<<GetQaxis()<<endl;</span> |
| 418 |
|
|
<a name="l00397"></a>00397 <span class="comment">// cout << " Q strip :"<<GetQstrip()<<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<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 |
pam-fi |
1.2 |
<a name="l00421"></a>00421 dof = GetN(); <span class="comment">//previous fit</span> |
| 443 |
|
|
<a name="l00422"></a>00422 P0 = par[0]; |
| 444 |
|
|
<a name="l00423"></a>00423 P1 = par[1]; |
| 445 |
|
|
<a name="l00424"></a>00424 Reset(); <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<cevent->GetN(); ih++){<span class="comment">//loop over selected hits</span> |
| 447 |
|
|
<a name="l00426"></a>00426 <span class="keywordtype">float</span> x = cevent->xycoord[ih]; |
| 448 |
|
|
<a name="l00427"></a>00427 <span class="keywordtype">float</span> z = cevent->zcoord[ih]; |
| 449 |
|
|
<a name="l00428"></a>00428 <span class="keywordtype">float</span> q = cevent->q[ih]; |
| 450 |
|
|
<a name="l00429"></a>00429 <span class="keywordtype">int</span> ip = cevent->idp[ih]; |
| 451 |
pam-fi |
1.1 |
<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<<d<<endl;</span> |
| 453 |
|
|
<a name="l00432"></a>00432 <span class="keywordflow">if</span>( d < ttoll && (niter==1 || (niter>1 && pok[ip]==1)) ){ |
| 454 |
|
|
<a name="l00433"></a>00433 <span class="comment">// if( d < 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 << ip << " "<<W<<endl;</span> |
| 457 |
pam-fi |
1.2 |
<a name="l00436"></a>00436 Add(z,x,q,W); |
| 458 |
pam-fi |
1.1 |
<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<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<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 |
pam-fi |
1.2 |
<a name="l00456"></a>00456 <span class="keywordflow">if</span>(GetN()<3)<span class="keywordflow">break</span>; |
| 478 |
|
|
<a name="l00457"></a>00457 Fit(); |
| 479 |
|
|
<a name="l00458"></a>00458 }<span class="keywordflow">while</span>(niter==1 || GetN() != dof); |
| 480 |
pam-fi |
1.1 |
<a name="l00459"></a>00459 <span class="comment">// cout << " >>> "<<GetN()<<endl;</span> |
| 481 |
pam-fi |
1.2 |
<a name="l00460"></a>00460 <span class="keywordflow">if</span>(GetN()<3)<span class="keywordflow">return</span> 0; |
| 482 |
|
|
<a name="l00461"></a>00461 Fit(); |
| 483 |
pam-fi |
1.1 |
<a name="l00462"></a>00462 <span class="comment">// cout << "(2)"<<endl;</span> |
| 484 |
|
|
<a name="l00463"></a>00463 <span class="comment">// cout << " n. hits :"<<GetN()<<endl;</span> |
| 485 |
|
|
<a name="l00464"></a>00464 <span class="comment">// cout << " P0 P1 :"<<par[0]<< " "<<par[1]<<endl;</span> |
| 486 |
|
|
<a name="l00465"></a>00465 <span class="comment">// cout << " Chi2 :"<<GetChi2()<<endl;</span> |
| 487 |
|
|
<a name="l00466"></a>00466 <span class="comment">// cout << " Q axis :"<<GetQaxis()<<endl;</span> |
| 488 |
|
|
<a name="l00467"></a>00467 <span class="comment">// cout << " Q strip :"<<GetQstrip()<<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 |
pam-fi |
1.2 |
<a name="l00473"></a>00473 P0 = par[0]; |
| 495 |
|
|
<a name="l00474"></a>00474 P1 = par[1]; |
| 496 |
|
|
<a name="l00475"></a>00475 Reset(); |
| 497 |
pam-fi |
1.1 |
<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<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 |
pam-fi |
1.2 |
<a name="l00483"></a>00483 <span class="keywordflow">for</span>(<span class="keywordtype">int</span> ih=0;ih<cevent->GetN(); ih++){ |
| 505 |
|
|
<a name="l00484"></a>00484 <span class="keywordtype">float</span> x = cevent->xycoord[ih]; |
| 506 |
|
|
<a name="l00485"></a>00485 <span class="keywordtype">float</span> z = cevent->zcoord[ih]; |
| 507 |
pam-fi |
1.1 |
<a name="l00486"></a>00486 <span class="comment">//float q = cevent->q[ih];</span> |
| 508 |
pam-fi |
1.2 |
<a name="l00487"></a>00487 <span class="keywordtype">int</span> ip = cevent->idp[ih]; |
| 509 |
pam-fi |
1.1 |
<a name="l00488"></a>00488 <span class="comment">//int is = cevent->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 < toll && d<dmin[ip] && 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<22; ip++){ |
| 519 |
|
|
<a name="l00498"></a>00498 <span class="keywordflow">if</span>(closest[ip]!=-1){ |
| 520 |
pam-fi |
1.2 |
<a name="l00499"></a>00499 <span class="keywordtype">float</span> x = cevent->xycoord[closest[ip]]; |
| 521 |
|
|
<a name="l00500"></a>00500 <span class="keywordtype">float</span> z = cevent->zcoord[closest[ip]]; |
| 522 |
|
|
<a name="l00501"></a>00501 <span class="keywordtype">float</span> q = cevent->q[closest[ip]]; |
| 523 |
pam-fi |
1.1 |
<a name="l00502"></a>00502 <span class="comment">//int is = cevent->ids[closest[ip]];</span> |
| 524 |
pam-fi |
1.2 |
<a name="l00503"></a>00503 Add(z,x,q,q); |
| 525 |
|
|
<a name="l00504"></a>00504 cog[ip] += x*q; |
| 526 |
|
|
<a name="l00505"></a>00505 qpl[ip] += q; |
| 527 |
pam-fi |
1.1 |
<a name="l00506"></a>00506 <span class="comment">// cout << ip << " -o- "<<is<<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 |
pam-fi |
1.2 |
<a name="l00512"></a>00512 <span class="keywordflow">for</span>(<span class="keywordtype">int</span> ih=0;ih<cevent->GetN(); ih++){ |
| 534 |
|
|
<a name="l00513"></a>00513 <span class="keywordtype">float</span> x = cevent->xycoord[ih]; |
| 535 |
|
|
<a name="l00514"></a>00514 <span class="keywordtype">float</span> z = cevent->zcoord[ih]; |
| 536 |
|
|
<a name="l00515"></a>00515 <span class="keywordtype">float</span> q = cevent->q[ih]; |
| 537 |
|
|
<a name="l00516"></a>00516 <span class="keywordtype">int</span> ip = cevent->idp[ih]; |
| 538 |
|
|
<a name="l00517"></a>00517 <span class="keywordtype">int</span> is = cevent->ids[ih]; |
| 539 |
pam-fi |
1.1 |
<a name="l00518"></a>00518 <span class="keywordflow">if</span>( closest[ip]!=-1 ){ |
| 540 |
pam-fi |
1.2 |
<a name="l00519"></a>00519 <span class="keywordtype">int</span> isc = cevent->ids[closest[ip]]; |
| 541 |
pam-fi |
1.1 |
<a name="l00520"></a>00520 <span class="keywordflow">if</span>( is == isc+1 || is == isc-1 ){ |
| 542 |
pam-fi |
1.2 |
<a name="l00521"></a>00521 Add(z,x,q,q); |
| 543 |
|
|
<a name="l00522"></a>00522 cog[ip] += x*q; |
| 544 |
|
|
<a name="l00523"></a>00523 qpl[ip] += q; |
| 545 |
pam-fi |
1.1 |
<a name="l00524"></a>00524 <span class="comment">// cout << ip << " -+- "<<is<<endl;</span> |
| 546 |
|
|
<a name="l00525"></a>00525 } |
| 547 |
|
|
<a name="l00526"></a>00526 } |
| 548 |
|
|
<a name="l00527"></a>00527 } |
| 549 |
pam-fi |
1.2 |
<a name="l00528"></a>00528 Fit(); |
| 550 |
pam-fi |
1.1 |
<a name="l00529"></a>00529 |
| 551 |
pam-fi |
1.2 |
<a name="l00530"></a>00530 <span class="keywordflow">if</span>(GetN()<3)<span class="keywordflow">return</span> 0; |
| 552 |
pam-fi |
1.1 |
<a name="l00531"></a>00531 <span class="keywordflow">for</span>(<span class="keywordtype">int</span> ip=0; ip<22; ip++){ |
| 553 |
pam-fi |
1.2 |
<a name="l00532"></a>00532 <span class="keywordflow">if</span>(qpl[ip]!=0){ |
| 554 |
|
|
<a name="l00533"></a>00533 cog[ip]=cog[ip]/qpl[ip]; |
| 555 |
pam-fi |
1.1 |
<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 << "(3)"<<endl;</span> |
| 560 |
|
|
<a name="l00539"></a>00539 <span class="comment">// cout << " n. hits :"<<GetN()<<endl;</span> |
| 561 |
|
|
<a name="l00540"></a>00540 <span class="comment">// cout << " P0 P1 :"<<par[0]<< " "<<par[1]<<endl;</span> |
| 562 |
|
|
<a name="l00541"></a>00541 <span class="comment">// cout << " Chi2 :"<<GetChi2()<<endl;</span> |
| 563 |
|
|
<a name="l00542"></a>00542 <span class="comment">// cout << " Q axis :"<<GetQaxis()<<endl;</span> |
| 564 |
|
|
<a name="l00543"></a>00543 <span class="comment">// cout << " Q strip :"<<GetQstrip()<<endl;</span> |
| 565 |
|
|
<a name="l00544"></a>00544 <span class="comment">// cout << " ---------------------------"<<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 |
pam-fi |
1.2 |
<a name="l00549"></a>00549 <span class="keywordtype">int</span> FitAxis( CaloLevel1* calo , Int_t view , Float_t toll ){ <span class="keywordflow">return</span> FitAxis(calo,view,toll,0); }; |
| 571 |
pam-fi |
1.1 |
<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 |
pam-fi |
1.2 |
<a name="l00556"></a>00556 <span class="keywordtype">int</span> FitShower( CaloLevel1* calo , Int_t view , Float_t toll , Bool_t usemechanicalalignement ){ |
| 578 |
pam-fi |
1.1 |
<a name="l00557"></a>00557 |
| 579 |
|
|
<a name="l00558"></a>00558 <span class="comment">// cout << "CaloAxis::FitAxis(...)"<<endl;</span> |
| 580 |
|
|
<a name="l00559"></a>00559 |
| 581 |
pam-fi |
1.2 |
<a name="l00560"></a>00560 Set(calo,view,usemechanicalalignement); |
| 582 |
pam-fi |
1.1 |
<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 |
pam-fi |
1.2 |
<a name="l00565"></a>00565 <span class="keywordflow">for</span>(<span class="keywordtype">int</span> ih=0;ih<cevent->GetN(); ih++){ |
| 587 |
|
|
<a name="l00566"></a>00566 <span class="keywordtype">float</span> x = cevent->xycoord[ih]; |
| 588 |
pam-fi |
1.1 |
<a name="l00567"></a>00567 <span class="comment">//float z = cevent->zcoord[ih];</span> |
| 589 |
pam-fi |
1.2 |
<a name="l00568"></a>00568 <span class="keywordtype">float</span> q = cevent->q[ih]; |
| 590 |
|
|
<a name="l00569"></a>00569 <span class="keywordtype">int</span> ip = cevent->idp[ih]; |
| 591 |
|
|
<a name="l00570"></a>00570 cog[ip] += x*q; |
| 592 |
|
|
<a name="l00571"></a>00571 qpl[ip] += q; |
| 593 |
pam-fi |
1.1 |
<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<22; ip++){ |
| 596 |
pam-fi |
1.2 |
<a name="l00575"></a>00575 <span class="keywordflow">if</span>(qpl[ip]!=0)cog[ip]=cog[ip]/qpl[ip]; |
| 597 |
pam-fi |
1.1 |
<a name="l00576"></a>00576 <span class="comment">// cout << ip << " > "<<qpl[ip]<<" "<<cog[ip]<<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<22; ip++)if(qpl[ip]>qplmax){ |
| 607 |
|
|
<a name="l00586"></a>00586 ipmax=ip; |
| 608 |
pam-fi |
1.2 |
<a name="l00587"></a>00587 qplmax=qpl[ip]; |
| 609 |
pam-fi |
1.1 |
<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 |
pam-fi |
1.2 |
<a name="l00591"></a>00591 <span class="keywordflow">for</span>(<span class="keywordtype">int</span> ih=0;ih<cevent->GetN(); ih++){ |
| 613 |
|
|
<a name="l00592"></a>00592 <span class="keywordtype">float</span> x = cevent->xycoord[ih]; |
| 614 |
|
|
<a name="l00593"></a>00593 <span class="keywordtype">float</span> z = cevent->zcoord[ih]; |
| 615 |
|
|
<a name="l00594"></a>00594 <span class="keywordtype">float</span> q = cevent->q[ih]; |
| 616 |
|
|
<a name="l00595"></a>00595 <span class="keywordtype">int</span> ip = cevent->idp[ih]; |
| 617 |
pam-fi |
1.1 |
<a name="l00596"></a>00596 <span class="comment">// if( (ipmax>3 && ip<=ipmax) || ipmax<=3 )Add(z,x,q);</span> |
| 618 |
|
|
<a name="l00597"></a>00597 <span class="keywordflow">if</span>( ip<=ipmax || ip<=3 ){ |
| 619 |
pam-fi |
1.2 |
<a name="l00598"></a>00598 Add(z,x,q); |
| 620 |
pam-fi |
1.1 |
<a name="l00599"></a>00599 } |
| 621 |
|
|
<a name="l00600"></a>00600 } |
| 622 |
pam-fi |
1.2 |
<a name="l00601"></a>00601 <span class="keywordflow">if</span>(GetN()<3)<span class="keywordflow">return</span> 0; |
| 623 |
|
|
<a name="l00602"></a>00602 Fit(); |
| 624 |
pam-fi |
1.1 |
<a name="l00603"></a>00603 <span class="comment">// cout << " n. hits :"<<GetN()<<endl;</span> |
| 625 |
|
|
<a name="l00604"></a>00604 <span class="comment">// cout << " P0 P1 :"<<par[0]<< " "<<par[1]<<endl;</span> |
| 626 |
|
|
<a name="l00605"></a>00605 <span class="comment">// cout << " Chi2 :"<<GetChi2()<<endl;</span> |
| 627 |
|
|
<a name="l00606"></a>00606 <span class="comment">// cout << " Q axis :"<<GetQaxis()<<endl;</span> |
| 628 |
|
|
<a name="l00607"></a>00607 <span class="comment">// cout << " Q strip :"<<GetQstrip()<<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<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 |
pam-fi |
1.2 |
<a name="l00632"></a>00632 dof = GetN(); <span class="comment">//previous fit</span> |
| 654 |
|
|
<a name="l00633"></a>00633 P0 = par[0]; |
| 655 |
|
|
<a name="l00634"></a>00634 P1 = par[1]; |
| 656 |
|
|
<a name="l00635"></a>00635 Reset(); <span class="comment">//fit reset</span> |
| 657 |
pam-fi |
1.1 |
<a name="l00636"></a>00636 <span class="keywordflow">for</span>(<span class="keywordtype">int</span> i=0; i<22; i++){ |
| 658 |
pam-fi |
1.2 |
<a name="l00637"></a>00637 cog[i]=0; |
| 659 |
|
|
<a name="l00638"></a>00638 qpl[i]=0; |
| 660 |
pam-fi |
1.1 |
<a name="l00639"></a>00639 } |
| 661 |
pam-fi |
1.2 |
<a name="l00640"></a>00640 <span class="keywordflow">for</span>(<span class="keywordtype">int</span> ih=0;ih<cevent->GetN(); ih++){ |
| 662 |
|
|
<a name="l00641"></a>00641 <span class="keywordtype">float</span> x = cevent->xycoord[ih]; |
| 663 |
|
|
<a name="l00642"></a>00642 <span class="keywordtype">float</span> z = cevent->zcoord[ih]; |
| 664 |
|
|
<a name="l00643"></a>00643 <span class="keywordtype">float</span> q = cevent->q[ih]; |
| 665 |
|
|
<a name="l00644"></a>00644 <span class="keywordtype">int</span> ip = cevent->idp[ih]; |
| 666 |
pam-fi |
1.1 |
<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<<d<<endl;</span> |
| 668 |
|
|
<a name="l00647"></a>00647 <span class="keywordflow">if</span>( |
| 669 |
|
|
<a name="l00648"></a>00648 (niter==1 || (niter>1 && pok[ip]==1)) && |
| 670 |
|
|
<a name="l00649"></a>00649 (ip<=ipmax || ip<=3) && |
| 671 |
|
|
<a name="l00650"></a>00650 <span class="comment">// ((ipmax>3 && ip<=ipmax) || ipmax<=3) &&</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 && d<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 << ip << " "<<W<<endl;</span> |
| 681 |
|
|
<a name="l00660"></a>00660 <span class="keywordflow">if</span>( d < ttoll || ( niter==2 && ih==closest[ip] )){ |
| 682 |
|
|
<a name="l00661"></a>00661 <span class="comment">/* if(ip==ipmax && ipmax>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 |
pam-fi |
1.2 |
<a name="l00663"></a>00663 cog[ip] += x*q; |
| 685 |
|
|
<a name="l00664"></a>00664 qpl[ip] += q; |
| 686 |
|
|
<a name="l00665"></a>00665 Add(z,x,q,W); |
| 687 |
pam-fi |
1.1 |
<a name="l00666"></a>00666 pok[ip]=1; |
| 688 |
pam-fi |
1.2 |
<a name="l00667"></a>00667 cout << ip << <span class="stringliteral">" "</span><<qpl[ip]<<endl; |
| 689 |
pam-fi |
1.1 |
<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<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<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 |
pam-fi |
1.2 |
<a name="l00685"></a>00685 <span class="keywordflow">for</span>(<span class="keywordtype">int</span> ip=0; ip<22; ip++)cout << qpl[ip] << <span class="stringliteral">" "</span>; |
| 707 |
pam-fi |
1.1 |
<a name="l00686"></a>00686 cout << endl; |
| 708 |
|
|
<a name="l00687"></a>00687 |
| 709 |
|
|
<a name="l00688"></a>00688 |
| 710 |
pam-fi |
1.2 |
<a name="l00689"></a>00689 <span class="keywordflow">for</span>(<span class="keywordtype">int</span> ip=0; ip<22; ip++)<span class="keywordflow">if</span>(qpl[ip]!=0)cog[ip]=cog[ip]/qpl[ip]; |
| 711 |
pam-fi |
1.1 |
<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 |
pam-fi |
1.2 |
<a name="l00693"></a>00693 <span class="keywordflow">if</span>(GetN()<3)<span class="keywordflow">break</span>; |
| 715 |
pam-fi |
1.1 |
<a name="l00694"></a>00694 |
| 716 |
pam-fi |
1.2 |
<a name="l00695"></a>00695 Fit(); |
| 717 |
pam-fi |
1.1 |
<a name="l00696"></a>00696 |
| 718 |
pam-fi |
1.2 |
<a name="l00697"></a>00697 }<span class="keywordflow">while</span>(niter==1 || GetN() != dof); |
| 719 |
pam-fi |
1.1 |
<a name="l00698"></a>00698 <span class="comment">// cout << " >>> "<<GetN()<<endl;</span> |
| 720 |
pam-fi |
1.2 |
<a name="l00699"></a>00699 <span class="keywordflow">if</span>(GetN()<3)<span class="keywordflow">return</span> 0; |
| 721 |
|
|
<a name="l00700"></a>00700 Fit(); |
| 722 |
pam-fi |
1.1 |
<a name="l00701"></a>00701 <span class="comment">/* // cout << " n. hits :"<<GetN()<<endl; */</span> |
| 723 |
|
|
<a name="l00702"></a>00702 <span class="comment">/* // cout << " P0 P1 :"<<par[0]<< " "<<par[1]<<endl; */</span> |
| 724 |
|
|
<a name="l00703"></a>00703 <span class="comment">/* // cout << " Chi2 :"<<GetChi2()<<endl; */</span> |
| 725 |
|
|
<a name="l00704"></a>00704 <span class="comment">/* // cout << " Q axis :"<<GetQaxis()<<endl; */</span> |
| 726 |
|
|
<a name="l00705"></a>00705 <span class="comment">/* // cout << " Q strip :"<<GetQstrip()<<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<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<cevent->GetN(); ih++){ */</span> |
| 743 |
|
|
<a name="l00722"></a>00722 <span class="comment">/* float x = cevent->xycoord[ih]; */</span> |
| 744 |
|
|
<a name="l00723"></a>00723 <span class="comment">/* float z = cevent->zcoord[ih]; */</span> |
| 745 |
|
|
<a name="l00724"></a>00724 <span class="comment">/* float q = cevent->q[ih]; */</span> |
| 746 |
|
|
<a name="l00725"></a>00725 <span class="comment">/* int ip = cevent->idp[ih]; */</span> |
| 747 |
|
|
<a name="l00726"></a>00726 <span class="comment">/* int is = cevent->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 < toll && d<dmin[ip] && 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<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->xycoord[closest[ip]]; */</span> |
| 759 |
|
|
<a name="l00738"></a>00738 <span class="comment">/* float z = cevent->zcoord[closest[ip]]; */</span> |
| 760 |
|
|
<a name="l00739"></a>00739 <span class="comment">/* float q = cevent->q[closest[ip]]; */</span> |
| 761 |
|
|
<a name="l00740"></a>00740 <span class="comment">/* int is = cevent->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 << ip << " -o- "<<is<<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<cevent->GetN(); ih++){ */</span> |
| 770 |
|
|
<a name="l00749"></a>00749 <span class="comment">/* float x = cevent->xycoord[ih]; */</span> |
| 771 |
|
|
<a name="l00750"></a>00750 <span class="comment">/* float z = cevent->zcoord[ih]; */</span> |
| 772 |
|
|
<a name="l00751"></a>00751 <span class="comment">/* float q = cevent->q[ih]; */</span> |
| 773 |
|
|
<a name="l00752"></a>00752 <span class="comment">/* int ip = cevent->idp[ih]; */</span> |
| 774 |
|
|
<a name="l00753"></a>00753 <span class="comment">/* int is = cevent->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->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 << ip << " -+- "<<is<<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()<3)return 0; */</span> |
| 788 |
|
|
<a name="l00767"></a>00767 <span class="comment">/* for(int ip=0; ip<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 << " n. hits :"<<GetN()<<endl;</span> |
| 796 |
|
|
<a name="l00775"></a>00775 <span class="comment">// cout << " P0 P1 :"<<par[0]<< " "<<par[1]<<endl;</span> |
| 797 |
|
|
<a name="l00776"></a>00776 <span class="comment">// cout << " Chi2 :"<<GetChi2()<<endl;</span> |
| 798 |
|
|
<a name="l00777"></a>00777 <span class="comment">// cout << " Q axis :"<<GetQaxis()<<endl;</span> |
| 799 |
|
|
<a name="l00778"></a>00778 <span class="comment">// cout << " Q strip :"<<GetQstrip()<<endl;</span> |
| 800 |
|
|
<a name="l00779"></a>00779 <span class="comment">// cout << " ---------------------------"<<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 |
pam-fi |
1.2 |
<a name="l00784"></a>00784 <span class="keywordtype">int</span> FitShower( CaloLevel1* calo , Int_t view , Float_t toll ){ <span class="keywordflow">return</span> FitAxis(calo,view,toll,0); }; |
| 806 |
pam-fi |
1.1 |
<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 |
pam-fi |
1.2 |
<hr size="1"><address style="text-align: right;"><small>Generated on Fri Jun 5 15:08:58 2009 for PamCut by |
| 820 |
pam-fi |
1.1 |
<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> |
Parent Directory
| 
Revision Log