| 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 |
|
|
<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 <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 |
|
|
<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<int> <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<int> <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<float> <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<float> <a class="code" href="classCaloEvent.html#340ff2ae2167b37d0b6acc56b1475daf">zcoord</a>; |
| 62 |
|
|
<a name="l00041"></a><a class="code" href="classCaloEvent.html#b48403ffdb4114d57e4ea4058438e814">00041</a> vector<float> <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 << " 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 |
|
|
<a name="l00096"></a>00096 <span class="keywordflow">if</span>(fq>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<float> <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<float> <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<float> <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<float> <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<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<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<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->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>-><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>-><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 << " CaloAxis::Reset()"<<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 << sel << 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>->Delete(); |
| 211 |
|
|
<a name="l00190"></a>00190 <span class="comment">// cout << axis << 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>->Delete(); |
| 213 |
|
|
<a name="l00192"></a>00192 <span class="comment">// if(cevent)cevent->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>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>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>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 << " 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 |
|
|
<a name="l00237"></a>00237 <span class="keywordflow">for</span>(Int_t i=0; i<(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<3)<span class="keywordflow">return</span> 0; |
| 269 |
|
|
<a name="l00248"></a>00248 <span class="comment">// cout << np << " points fitted -- ok "<<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<(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 |
|
|
<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<(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<(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<(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<(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<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 |
|
|
<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<(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<(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 |
|
|
<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<<a class="code" href="classCaloAxis.html#fb48c08633030431b7e509ff1b34558d">cevent</a>-><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>-><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>-><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>-><a class="code" href="classCaloEvent.html#688d5d2803d4688d4f1ef72cbaf42dcd">q</a>[ih]; |
| 329 |
|
|
<a name="l00308"></a>00308 <span class="comment">//int ip = cevent->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 < 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>=0 && ip<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>=0 && ip<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<(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<(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]<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<(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]>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 << " CaloAxis::Draw()"<<endl;</span> |
| 374 |
|
|
<a name="l00353"></a>00353 <span class="keywordflow">if</span>(<a class="code" href="classCaloAxis.html#be448c5906fb90329316b50796d2c721">GetN</a>()<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>->SetMarkerSize(0.3); |
| 378 |
|
|
<a name="l00357"></a>00357 <a class="code" href="classCaloAxis.html#9e4372d0b9d9a3ec17fc2e1570039273">sel</a>->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>->SetLineWidth(1); |
| 382 |
|
|
<a name="l00361"></a>00361 <a class="code" href="classCaloAxis.html#7ca3199fb8575c566821219253d02ccf">axis</a>->SetLineColor(5); |
| 383 |
|
|
<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 |
|
|
<a name="l00366"></a>00366 <a class="code" href="classCaloAxis.html#9e4372d0b9d9a3ec17fc2e1570039273">sel</a>->Draw(); |
| 388 |
|
|
<a name="l00367"></a>00367 <a class="code" href="classCaloAxis.html#7ca3199fb8575c566821219253d02ccf">axis</a>->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 << "CaloAxis::FitAxis(...)"<<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<<a class="code" href="classCaloAxis.html#fb48c08633030431b7e509ff1b34558d">cevent</a>-><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>-><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>-><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>-><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>()<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 << "(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 |
|
|
<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<<a class="code" href="classCaloAxis.html#fb48c08633030431b7e509ff1b34558d">cevent</a>-><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>-><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>-><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>-><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>-><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<<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 |
|
|
<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<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 |
|
|
<a name="l00456"></a>00456 <span class="keywordflow">if</span>(<a class="code" href="classCaloAxis.html#be448c5906fb90329316b50796d2c721">GetN</a>()<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 << " >>> "<<GetN()<<endl;</span> |
| 481 |
|
|
<a name="l00460"></a>00460 <span class="keywordflow">if</span>(<a class="code" href="classCaloAxis.html#be448c5906fb90329316b50796d2c721">GetN</a>()<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 << "(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 |
|
|
<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<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<<a class="code" href="classCaloAxis.html#fb48c08633030431b7e509ff1b34558d">cevent</a>-><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>-><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>-><a class="code" href="classCaloEvent.html#340ff2ae2167b37d0b6acc56b1475daf">zcoord</a>[ih]; |
| 507 |
|
|
<a name="l00486"></a>00486 <span class="comment">//float q = cevent->q[ih];</span> |
| 508 |
|
|
<a name="l00487"></a>00487 <span class="keywordtype">int</span> ip = <a class="code" href="classCaloAxis.html#fb48c08633030431b7e509ff1b34558d">cevent</a>-><a class="code" href="classCaloEvent.html#3ad199ffd33f6300bff64e854c6ad652">idp</a>[ih]; |
| 509 |
|
|
<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 |
|
|
<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>-><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>-><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>-><a class="code" href="classCaloEvent.html#688d5d2803d4688d4f1ef72cbaf42dcd">q</a>[closest[ip]]; |
| 523 |
|
|
<a name="l00502"></a>00502 <span class="comment">//int is = cevent->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 << 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 |
|
|
<a name="l00512"></a>00512 <span class="keywordflow">for</span>(<span class="keywordtype">int</span> ih=0;ih<<a class="code" href="classCaloAxis.html#fb48c08633030431b7e509ff1b34558d">cevent</a>-><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>-><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>-><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>-><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>-><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>-><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>-><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 << ip << " -+- "<<is<<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>()<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<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 << "(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 |
|
|
<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 << "CaloAxis::FitAxis(...)"<<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<<a class="code" href="classCaloAxis.html#fb48c08633030431b7e509ff1b34558d">cevent</a>-><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>-><a class="code" href="classCaloEvent.html#b48403ffdb4114d57e4ea4058438e814">xycoord</a>[ih]; |
| 588 |
|
|
<a name="l00567"></a>00567 <span class="comment">//float z = cevent->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>-><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>-><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<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 << 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 |
|
|
<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<<a class="code" href="classCaloAxis.html#fb48c08633030431b7e509ff1b34558d">cevent</a>-><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>-><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>-><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>-><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>-><a class="code" href="classCaloEvent.html#3ad199ffd33f6300bff64e854c6ad652">idp</a>[ih]; |
| 617 |
|
|
<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 |
|
|
<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>()<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 << " 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 |
|
|
<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<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<<a class="code" href="classCaloAxis.html#fb48c08633030431b7e509ff1b34558d">cevent</a>-><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>-><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>-><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>-><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>-><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<<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 |
|
|
<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 << ip << <span class="stringliteral">" "</span><<<a class="code" href="classCaloAxis.html#2ff6400f0393e7ab05627a438c6eeca5">qpl</a>[ip]<<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<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 |
|
|
<a name="l00685"></a>00685 <span class="keywordflow">for</span>(<span class="keywordtype">int</span> ip=0; ip<22; ip++)cout << <a class="code" href="classCaloAxis.html#2ff6400f0393e7ab05627a438c6eeca5">qpl</a>[ip] << <span class="stringliteral">" "</span>; |
| 707 |
|
|
<a name="l00686"></a>00686 cout << 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<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>()<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 << " >>> "<<GetN()<<endl;</span> |
| 720 |
|
|
<a name="l00699"></a>00699 <span class="keywordflow">if</span>(<a class="code" href="classCaloAxis.html#be448c5906fb90329316b50796d2c721">GetN</a>()<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 << " 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 |
|
|
<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 |
| 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> |
Parent Directory
| 
Revision Log