| 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: TofCuts/TofPatternCut/TofPatternCut.cpp 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>TofCuts/TofPatternCut/TofPatternCut.cpp</h1><a href="TofPatternCut_8cpp.html">Go to the documentation of this file.</a><div class="fragment"><pre class="fragment"><a name="l00001"></a>00001 <span class="comment">/*</span> |
| 23 |
<a name="l00002"></a>00002 <span class="comment"> * TofPatternCut.cpp</span> |
| 24 |
<a name="l00003"></a>00003 <span class="comment"> *</span> |
| 25 |
<a name="l00004"></a>00004 <span class="comment"> * Created on: 12-mar-2009</span> |
| 26 |
<a name="l00005"></a>00005 <span class="comment"> * Author: Elena Vannuccini</span> |
| 27 |
<a name="l00006"></a>00006 <span class="comment"> */</span> |
| 28 |
<a name="l00007"></a>00007 |
| 29 |
<a name="l00010"></a>00010 <span class="preprocessor">#include "<a class="code" href="TofPatternCut_8h.html">TofPatternCut.h</a>"</span> |
| 30 |
<a name="l00011"></a>00011 |
| 31 |
<a name="l00012"></a><a class="code" href="classTofPatternCut.html#b960628bc87a72be3c8aedab6b036dce">00012</a> <span class="keywordtype">int</span> <a class="code" href="classTofPatternCut.html#b960628bc87a72be3c8aedab6b036dce" title="The ToF-pattern check.">TofPatternCut::Check</a>(PamLevel2 *event) { |
| 32 |
<a name="l00013"></a>00013 |
| 33 |
<a name="l00014"></a>00014 <span class="keywordflow">if</span> (!event) |
| 34 |
<a name="l00015"></a>00015 <span class="keywordflow">return</span> (0); |
| 35 |
<a name="l00016"></a>00016 <span class="comment">//</span> |
| 36 |
<a name="l00017"></a>00017 ToFLevel2 *tofl2 = <span class="keyword">event</span>->GetToFLevel2(); |
| 37 |
<a name="l00018"></a>00018 <span class="keywordflow">if</span> (!tofl2) |
| 38 |
<a name="l00019"></a>00019 <span class="keywordflow">return</span> (0); |
| 39 |
<a name="l00020"></a>00020 ToFTrkVar *tof = 0; |
| 40 |
<a name="l00021"></a>00021 <span class="keywordflow">if</span> (!<a class="code" href="classTofPatternCut.html#28f5094cd7898788d4e490e75384b538">_standalone</a> && <a class="code" href="classTofPatternCut.html#9c8d4ac314f4f3a0b62fae04ecc32fd7">_notrk</a> >= 0 && _notrk < event->GetTrkLevel2()->GetNTracks()) |
| 41 |
<a name="l00022"></a>00022 tof = <span class="keyword">event</span>->GetTrack(<a class="code" href="classTofPatternCut.html#9c8d4ac314f4f3a0b62fae04ecc32fd7">_notrk</a>)->GetToFTrack(); |
| 42 |
<a name="l00023"></a>00023 <span class="keywordflow">if</span> (<a class="code" href="classTofPatternCut.html#28f5094cd7898788d4e490e75384b538">_standalone</a>) |
| 43 |
<a name="l00024"></a>00024 tof = <span class="keyword">event</span>->GetToFStoredTrack(-1); |
| 44 |
<a name="l00025"></a>00025 <span class="keywordflow">if</span> (!tof) |
| 45 |
<a name="l00026"></a>00026 <span class="keywordflow">return</span> (0); |
| 46 |
<a name="l00027"></a>00027 |
| 47 |
<a name="l00028"></a>00028 <span class="preprocessor">#ifdef DEBUGPAMCUT</span> |
| 48 |
<a name="l00029"></a>00029 <span class="preprocessor"></span> cout << <span class="stringliteral">">>> TofPatternCut CHECK <<<"</span> << endl; |
| 49 |
<a name="l00030"></a>00030 <span class="preprocessor">#endif</span> |
| 50 |
<a name="l00031"></a>00031 <span class="preprocessor"></span> |
| 51 |
<a name="l00032"></a>00032 <span class="keywordtype">bool</span> TOF__OK = <span class="keyword">false</span>; |
| 52 |
<a name="l00033"></a>00033 |
| 53 |
<a name="l00034"></a>00034 <span class="preprocessor">#ifdef DEBUGPAMCUT</span> |
| 54 |
<a name="l00035"></a>00035 <span class="preprocessor"></span> cout << <span class="stringliteral">"Check topology:"</span> << endl; |
| 55 |
<a name="l00036"></a>00036 <span class="keywordflow">for</span> (<span class="keywordtype">int</span> ip = 0; ip < 6; ip++) |
| 56 |
<a name="l00037"></a>00037 cout << <span class="stringliteral">"plane "</span> << ip << <span class="stringliteral">" n.hit paddles "</span> << tofl2->GetNHitPaddles(ip) << endl; |
| 57 |
<a name="l00038"></a>00038 <span class="keywordflow">for</span> (<span class="keywordtype">int</span> ip = 0; ip < 6; ip++) |
| 58 |
<a name="l00039"></a>00039 cout << <span class="stringliteral">"plane "</span> << ip << <span class="stringliteral">" n.hit PMTs outside track "</span> << <a class="code" href="classTofPatternCut.html#2f9755367321e7f86f862506865f260c" title="Counts the number of hit PMTs (with tdc signal) outside track.">GetNHitPMTsOutsideTrack</a>(event, ip) << endl; |
| 59 |
<a name="l00040"></a>00040 <span class="preprocessor">#endif</span> |
| 60 |
<a name="l00041"></a>00041 <span class="preprocessor"></span> |
| 61 |
<a name="l00042"></a>00042 <span class="comment">// -----------------------------------------</span> |
| 62 |
<a name="l00043"></a>00043 <span class="comment">// at least one hit paddle on both S1 and S2</span> |
| 63 |
<a name="l00044"></a>00044 <span class="comment">// -----------------------------------------</span> |
| 64 |
<a name="l00045"></a>00045 <span class="keywordflow">if</span> ( |
| 65 |
<a name="l00046"></a>00046 <span class="comment">//-----------------------------------</span> |
| 66 |
<a name="l00047"></a>00047 <span class="comment">// not more than 1 paddle per layer</span> |
| 67 |
<a name="l00048"></a>00048 <span class="comment">//-----------------------------------</span> |
| 68 |
<a name="l00049"></a>00049 tofl2 && tofl2->GetNHitPaddles(0) <= 1 && <span class="comment">//S11</span> |
| 69 |
<a name="l00050"></a>00050 tofl2->GetNHitPaddles(1) <= 1 && <span class="comment">//S12</span> |
| 70 |
<a name="l00051"></a>00051 tofl2->GetNHitPaddles(2) <= 1 && <span class="comment">//S21</span> |
| 71 |
<a name="l00052"></a>00052 tofl2->GetNHitPaddles(3) <= 1 && <span class="comment">//S22</span> |
| 72 |
<a name="l00053"></a>00053 <span class="comment">// tofl2->GetNHitPaddles(4) <= 1 && //S31</span> |
| 73 |
<a name="l00054"></a>00054 <span class="comment">// tofl2->GetNHitPaddles(5) <= 1 && //S32</span> |
| 74 |
<a name="l00055"></a>00055 <span class="comment">//-----------------------------------</span> |
| 75 |
<a name="l00056"></a>00056 <span class="comment">// at least 1 paddle per plane</span> |
| 76 |
<a name="l00057"></a>00057 <span class="comment">//-----------------------------------</span> |
| 77 |
<a name="l00058"></a>00058 (tofl2->GetNHitPaddles(0) == 1 || tofl2->GetNHitPaddles(1) == 1) && <span class="comment">//S1</span> |
| 78 |
<a name="l00059"></a>00059 (tofl2->GetNHitPaddles(2) == 1 || tofl2->GetNHitPaddles(3) == 1) && <span class="comment">//S2</span> |
| 79 |
<a name="l00060"></a>00060 <span class="comment">// (tofl2->GetNHitPaddles(4)==1 || tofl2->GetNHitPaddles(5)==1) && //S3</span> |
| 80 |
<a name="l00061"></a>00061 <span class="comment">//-----------------------------------</span> |
| 81 |
<a name="l00062"></a>00062 <span class="comment">// n.hit pmts outside track</span> |
| 82 |
<a name="l00063"></a>00063 <span class="comment">//-----------------------------------</span> |
| 83 |
<a name="l00064"></a>00064 <span class="comment">// GetHitPMTsOutsideTrack(event,0)<2 && //S11 //troppo (puo` succedere )</span> |
| 84 |
<a name="l00065"></a>00065 <span class="comment">// GetHitPMTsOutsideTrack(event,1)<2 && //S12</span> |
| 85 |
<a name="l00066"></a>00066 <span class="comment">// GetHitPMTsOutsideTrack(event,2)<2 && //S21</span> |
| 86 |
<a name="l00067"></a>00067 <span class="comment">// GetHitPMTsOutsideTrack(event,3)<2 && //S22</span> |
| 87 |
<a name="l00068"></a>00068 <a class="code" href="classTofPatternCut.html#2f9755367321e7f86f862506865f260c" title="Counts the number of hit PMTs (with tdc signal) outside track.">GetNHitPMTsOutsideTrack</a>(event, 0) < 3 && <span class="comment">//S11</span> |
| 88 |
<a name="l00069"></a>00069 <a class="code" href="classTofPatternCut.html#2f9755367321e7f86f862506865f260c" title="Counts the number of hit PMTs (with tdc signal) outside track.">GetNHitPMTsOutsideTrack</a>(event, 1) < 3 && <span class="comment">//S12</span> |
| 89 |
<a name="l00070"></a>00070 <span class="comment">// GetHitPMTsOutsideTrack(event,2)<3 && //S21</span> |
| 90 |
<a name="l00071"></a>00071 <span class="comment">// GetHitPMTsOutsideTrack(event,3)<3 && //S22</span> |
| 91 |
<a name="l00072"></a>00072 <span class="keyword">true</span>) |
| 92 |
<a name="l00073"></a>00073 TOF__OK = <span class="keyword">true</span>; |
| 93 |
<a name="l00074"></a>00074 |
| 94 |
<a name="l00075"></a>00075 <span class="keywordflow">if</span> (!TOF__OK) { |
| 95 |
<a name="l00076"></a>00076 <span class="preprocessor">#ifdef DEBUGPAMCUT</span> |
| 96 |
<a name="l00077"></a>00077 <span class="preprocessor"></span> cout << <span class="stringliteral">"(((DISCARDED)))"</span> << endl; |
| 97 |
<a name="l00078"></a>00078 <span class="preprocessor">#endif</span> |
| 98 |
<a name="l00079"></a>00079 <span class="preprocessor"></span> <span class="keywordflow">return</span> (0); |
| 99 |
<a name="l00080"></a>00080 }; |
| 100 |
<a name="l00081"></a>00081 |
| 101 |
<a name="l00082"></a>00082 <span class="preprocessor">#ifdef DEBUGPAMCUT</span> |
| 102 |
<a name="l00083"></a>00083 <span class="preprocessor"></span> cout << <span class="stringliteral">"--> OK"</span> << endl; |
| 103 |
<a name="l00084"></a>00084 cout << <span class="stringliteral">"Check track concistency:"</span> << endl; |
| 104 |
<a name="l00085"></a>00085 <span class="preprocessor">#endif</span> |
| 105 |
<a name="l00086"></a>00086 <span class="preprocessor"></span> |
| 106 |
<a name="l00087"></a>00087 <span class="comment">// -----------------------------------------------</span> |
| 107 |
<a name="l00088"></a>00088 <span class="comment">// the hit paddles MUST be associated to the track</span> |
| 108 |
<a name="l00089"></a>00089 <span class="comment">// (in collaborazione con wolfgang)</span> |
| 109 |
<a name="l00090"></a>00090 <span class="comment">// -----------------------------------------------</span> |
| 110 |
<a name="l00091"></a>00091 <span class="comment">//</span> |
| 111 |
<a name="l00092"></a>00092 <span class="comment">// first build the map of hit PMTs along the track</span> |
| 112 |
<a name="l00093"></a>00093 <span class="comment">// and discard the event if there is an hit paddle,</span> |
| 113 |
<a name="l00094"></a>00094 <span class="comment">// but it is not the one associated to the track</span> |
| 114 |
<a name="l00095"></a>00095 <span class="comment">//</span> |
| 115 |
<a name="l00096"></a>00096 |
| 116 |
<a name="l00097"></a>00097 <span class="keywordtype">int</span> hitplane[] = { 0, 0, 0, 0, 0, 0 }; |
| 117 |
<a name="l00098"></a>00098 <span class="keywordtype">int</span> hitpaddle[] = { -1, -1, -1, -1, -1, -1 }; |
| 118 |
<a name="l00099"></a>00099 <span class="keywordflow">for</span> (Int_t ip = 0; ip < tof->npmttdc; ip++) {<span class="comment">//loop over tdc hits associated to the track</span> |
| 119 |
<a name="l00100"></a>00100 Int_t iplane = -1; |
| 120 |
<a name="l00101"></a>00101 Int_t ipaddle = -1; |
| 121 |
<a name="l00102"></a>00102 Int_t ipmt = -1; |
| 122 |
<a name="l00103"></a>00103 Int_t <span class="keywordtype">id</span> = tof->pmttdc[ip];<span class="comment">//pmt id</span> |
| 123 |
<a name="l00104"></a>00104 tofl2->GetPMTName(<span class="keywordtype">id</span>, iplane, ipaddle, ipmt); |
| 124 |
<a name="l00105"></a>00105 <span class="keywordflow">if</span> (tof->tdcflag[ip] == 0) { |
| 125 |
<a name="l00106"></a>00106 hitplane[iplane] = 1; <span class="comment">//there is a true tdc signal associated to the track</span> |
| 126 |
<a name="l00107"></a>00107 <span class="preprocessor">#ifdef DEBUGPAMCUT</span> |
| 127 |
<a name="l00108"></a>00108 <span class="preprocessor"></span> <span class="keywordflow">if</span> (hitpaddle[iplane] >= 0 && hitpaddle[iplane] != ipaddle) |
| 128 |
<a name="l00109"></a>00109 cout << <span class="stringliteral">"ORRORE!!!!"</span> << endl; |
| 129 |
<a name="l00110"></a>00110 <span class="preprocessor">#endif</span> |
| 130 |
<a name="l00111"></a>00111 <span class="preprocessor"></span> hitpaddle[iplane] = ipaddle;<span class="comment">//store the id of the paddle associated to the track</span> |
| 131 |
<a name="l00112"></a>00112 <span class="comment">//</span> |
| 132 |
<a name="l00113"></a>00113 <span class="comment">//</span> |
| 133 |
<a name="l00114"></a>00114 <span class="comment">//</span> |
| 134 |
<a name="l00115"></a>00115 <span class="keywordflow">if</span> (tofl2->GetNHitPaddles(iplane) > 0 && <span class="comment">//if there is a hit paddle in this plane...</span> |
| 135 |
<a name="l00116"></a>00116 !tofl2->HitPaddle(iplane, hitpaddle[iplane]) && <span class="comment">//... and the paddle associated to the track is not hit</span> |
| 136 |
<a name="l00117"></a>00117 <span class="keyword">true</span>) { |
| 137 |
<a name="l00118"></a>00118 <span class="preprocessor">#ifdef DEBUGPAMCUT</span> |
| 138 |
<a name="l00119"></a>00119 <span class="preprocessor"></span> cout << <span class="stringliteral">"hit plane "</span> << iplane << <span class="stringliteral">" -- track PMT "</span> << hitpaddle[iplane] << <span class="stringliteral">" is not hit "</span> << endl; |
| 139 |
<a name="l00120"></a>00120 cout << <span class="stringliteral">"(((DISCARDED)))"</span> << endl; |
| 140 |
<a name="l00121"></a>00121 <span class="preprocessor">#endif</span> |
| 141 |
<a name="l00122"></a>00122 <span class="preprocessor"></span> <span class="keywordflow">return</span> (0);<span class="comment">//discard the event</span> |
| 142 |
<a name="l00123"></a>00123 } |
| 143 |
<a name="l00124"></a>00124 } |
| 144 |
<a name="l00125"></a>00125 }; |
| 145 |
<a name="l00126"></a>00126 <span class="comment">//</span> |
| 146 |
<a name="l00127"></a>00127 <span class="comment">// hence check the planes with no PMTs associated to the track:</span> |
| 147 |
<a name="l00128"></a>00128 <span class="comment">// the event is discarded if there is a hit paddle</span> |
| 148 |
<a name="l00129"></a>00129 <span class="comment">//</span> |
| 149 |
<a name="l00130"></a>00130 <span class="keywordflow">for</span> (<span class="keywordtype">int</span> iplane = 0; iplane < 4; iplane++) {<span class="comment">//loop over S11/S12/S21/S22</span> |
| 150 |
<a name="l00131"></a>00131 <span class="keywordflow">if</span> (hitplane[iplane] == 0 && <span class="comment">//if there are no PMTs associated to the track...</span> |
| 151 |
<a name="l00132"></a>00132 tofl2->GetNHitPaddles(iplane) > 0 && <span class="comment">//...and the plane has hit paddles</span> |
| 152 |
<a name="l00133"></a>00133 <span class="keyword">true</span>) { |
| 153 |
<a name="l00134"></a>00134 <span class="preprocessor">#ifdef DEBUGPAMCUT</span> |
| 154 |
<a name="l00135"></a>00135 <span class="preprocessor"></span> cout << <span class="stringliteral">"hit plane "</span> << iplane << <span class="stringliteral">" has no PMTs associated to the track "</span> << endl; |
| 155 |
<a name="l00136"></a>00136 cout << <span class="stringliteral">"(((DISCARDED)))"</span> << endl; |
| 156 |
<a name="l00137"></a>00137 <span class="preprocessor">#endif</span> |
| 157 |
<a name="l00138"></a>00138 <span class="preprocessor"></span> <span class="keywordflow">return</span> (0);<span class="comment">//discard the event</span> |
| 158 |
<a name="l00139"></a>00139 } |
| 159 |
<a name="l00140"></a>00140 } |
| 160 |
<a name="l00141"></a>00141 |
| 161 |
<a name="l00142"></a>00142 <span class="comment">//</span> |
| 162 |
<a name="l00143"></a>00143 <span class="comment">// forse non c'e` bisogno di tutto cio`?</span> |
| 163 |
<a name="l00144"></a>00144 <span class="comment">//</span> |
| 164 |
<a name="l00145"></a>00145 <span class="comment">// for(int iplane=0; iplane<4; iplane++){//loop over S11/S12/S21/S22</span> |
| 165 |
<a name="l00146"></a>00146 <span class="comment">// //retrieve the id of the paddle traversed by the track</span> |
| 166 |
<a name="l00147"></a>00147 <span class="comment">// // (...there is some tolerance...)</span> |
| 167 |
<a name="l00148"></a>00148 <span class="comment">// int ipaddle = tofl2->GetPaddleIdOfTrack(tof->xtr_tof[iplane],tof->ytr_tof[iplane],iplane);</span> |
| 168 |
<a name="l00149"></a>00149 <span class="comment">// //check if the traversed paddle is hit</span> |
| 169 |
<a name="l00150"></a>00150 <span class="comment">// bool OK = true;</span> |
| 170 |
<a name="l00151"></a>00151 <span class="comment">// if(</span> |
| 171 |
<a name="l00152"></a>00152 <span class="comment">// tofl2->GetNHitPaddles(iplane)>0 &&//if there is a hit paddle in this plane...</span> |
| 172 |
<a name="l00153"></a>00153 <span class="comment">// !tofl2->HitPaddle(iplane,ipaddle) && //...and the paddle traversed by the track is not hit</span> |
| 173 |
<a name="l00154"></a>00154 <span class="comment">// true)OK = false;//..discard the event</span> |
| 174 |
<a name="l00155"></a>00155 <span class="comment">// /// hence try to recover some events...</span> |
| 175 |
<a name="l00156"></a>00156 <span class="comment">// if(</span> |
| 176 |
<a name="l00157"></a>00157 <span class="comment">// tofl2->GetNHitPaddles(iplane)>0 &&//if there is a hit paddle...</span> |
| 177 |
<a name="l00158"></a>00158 <span class="comment">// ipaddle==-1 &&//...and the track does not traverse any paddle...</span> |
| 178 |
<a name="l00159"></a>00159 <span class="comment">// tofl2->HitPaddle(iplane,hitpaddle[iplane]) &&//... BUT there are tdc signals belonging to a hit paddle</span> |
| 179 |
<a name="l00160"></a>00160 <span class="comment">// true)OK=true;//recover</span> |
| 180 |
<a name="l00161"></a>00161 |
| 181 |
<a name="l00162"></a>00162 <span class="comment">// if(!OK)return (0);</span> |
| 182 |
<a name="l00163"></a>00163 |
| 183 |
<a name="l00164"></a>00164 <span class="comment">// }</span> |
| 184 |
<a name="l00165"></a>00165 |
| 185 |
<a name="l00166"></a>00166 <span class="comment">// if( !TOF__OK )return (0);</span> |
| 186 |
<a name="l00167"></a>00167 |
| 187 |
<a name="l00168"></a>00168 <span class="preprocessor">#ifdef DEBUGPAMCUT</span> |
| 188 |
<a name="l00169"></a>00169 <span class="preprocessor"></span> cout << <span class="stringliteral">"--> OK"</span> << endl; |
| 189 |
<a name="l00170"></a>00170 <span class="preprocessor">#endif</span> |
| 190 |
<a name="l00171"></a>00171 <span class="preprocessor"></span> |
| 191 |
<a name="l00172"></a>00172 <span class="keywordflow">return</span> <a class="code" href="CommonDefs_8h.html#c857f4b76ae431837cbf61fd74863b6c" title="The implementations of PamCut::Check() must return this value if the event satisfy...">CUTOK</a>; <span class="comment">//good pattern</span> |
| 192 |
<a name="l00173"></a>00173 |
| 193 |
<a name="l00174"></a>00174 } |
| 194 |
<a name="l00175"></a>00175 <span class="comment">//===========================================================</span> |
| 195 |
<a name="l00176"></a>00176 <span class="comment">//</span> |
| 196 |
<a name="l00177"></a>00177 <span class="comment">//</span> |
| 197 |
<a name="l00178"></a>00178 <span class="comment">//</span> |
| 198 |
<a name="l00179"></a>00179 <span class="comment">//</span> |
| 199 |
<a name="l00180"></a>00180 <span class="comment">//</span> |
| 200 |
<a name="l00181"></a>00181 <span class="comment">//===========================================================</span> |
| 201 |
<a name="l00182"></a><a class="code" href="classTofPatternCut.html#2f9755367321e7f86f862506865f260c">00182</a> <span class="keywordtype">int</span> <a class="code" href="classTofPatternCut.html#2f9755367321e7f86f862506865f260c" title="Counts the number of hit PMTs (with tdc signal) outside track.">TofPatternCut::GetNHitPMTsOutsideTrack</a>(PamLevel2 *event, <span class="keywordtype">int</span> plane) { |
| 202 |
<a name="l00183"></a>00183 |
| 203 |
<a name="l00184"></a>00184 <span class="keywordflow">if</span> (!event) |
| 204 |
<a name="l00185"></a>00185 <span class="keywordflow">return</span> (0); |
| 205 |
<a name="l00186"></a>00186 <span class="comment">//</span> |
| 206 |
<a name="l00187"></a>00187 ToFLevel2 *tofl2 = <span class="keyword">event</span>->GetToFLevel2(); |
| 207 |
<a name="l00188"></a>00188 <span class="keywordflow">if</span> (!tofl2) |
| 208 |
<a name="l00189"></a>00189 <span class="keywordflow">return</span> (0); |
| 209 |
<a name="l00190"></a>00190 ToFTrkVar *tof = 0; |
| 210 |
<a name="l00191"></a>00191 <span class="keywordflow">if</span> (!<a class="code" href="classTofPatternCut.html#28f5094cd7898788d4e490e75384b538">_standalone</a> && <a class="code" href="classTofPatternCut.html#9c8d4ac314f4f3a0b62fae04ecc32fd7">_notrk</a> >= 0 && _notrk < event->GetTrkLevel2()->GetNTracks()) |
| 211 |
<a name="l00192"></a>00192 tof = <span class="keyword">event</span>->GetTrack(<a class="code" href="classTofPatternCut.html#9c8d4ac314f4f3a0b62fae04ecc32fd7">_notrk</a>)->GetToFTrack(); |
| 212 |
<a name="l00193"></a>00193 <span class="keywordflow">if</span> (<a class="code" href="classTofPatternCut.html#28f5094cd7898788d4e490e75384b538">_standalone</a>) |
| 213 |
<a name="l00194"></a>00194 tof = <span class="keyword">event</span>->GetToFStoredTrack(-1); |
| 214 |
<a name="l00195"></a>00195 <span class="keywordflow">if</span> (!tof) |
| 215 |
<a name="l00196"></a>00196 <span class="keywordflow">return</span> (0); |
| 216 |
<a name="l00197"></a>00197 |
| 217 |
<a name="l00198"></a>00198 <span class="keywordtype">int</span> nn = 0; |
| 218 |
<a name="l00199"></a>00199 <span class="comment">// cout << "%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%"<<endl;</span> |
| 219 |
<a name="l00200"></a>00200 <span class="keywordflow">for</span> (<span class="keywordtype">int</span> ip = 0; ip < tofl2->npmt(); ip++) { <span class="comment">//loop over ALL PMTs</span> |
| 220 |
<a name="l00201"></a>00201 Int_t iplane = -1; |
| 221 |
<a name="l00202"></a>00202 Int_t ipaddle = -1; |
| 222 |
<a name="l00203"></a>00203 Int_t ipmt = -1; |
| 223 |
<a name="l00204"></a>00204 Int_t <span class="keywordtype">id</span> = tofl2->GetToFPMT(ip)->pmt_id; <span class="comment">//PMT id</span> |
| 224 |
<a name="l00205"></a>00205 tofl2->GetPMTName(<span class="keywordtype">id</span>, iplane, ipaddle, ipmt); |
| 225 |
<a name="l00206"></a>00206 |
| 226 |
<a name="l00207"></a>00207 <span class="keywordflow">if</span> (iplane != plane) |
| 227 |
<a name="l00208"></a>00208 <span class="keywordflow">continue</span>; |
| 228 |
<a name="l00209"></a>00209 <span class="comment">// cout << endl<<plane << " -- "<<id;</span> |
| 229 |
<a name="l00210"></a>00210 <span class="keywordtype">float</span> tdc = tofl2->GetToFPMT(ip)->tdc; |
| 230 |
<a name="l00211"></a>00211 <span class="keywordflow">if</span> (tdc >= 4095.) |
| 231 |
<a name="l00212"></a>00212 <span class="keywordflow">continue</span>; |
| 232 |
<a name="l00213"></a>00213 <span class="comment">// spiegazione:</span> |
| 233 |
<a name="l00214"></a>00214 <span class="comment">// quando c'e` adc pile-up, succede che un pmt abbia un segnale adc (adc<4095) ma NON tdc (tdc==4095),</span> |
| 234 |
<a name="l00215"></a>00215 <span class="comment">// che e` il contrario di quanto avviene di solito. quindi, se nel livello 2 c'e` un pmt colpito</span> |
| 235 |
<a name="l00216"></a>00216 <span class="comment">// e questo ha tdc==4095, implica che c'e` solo il segnale adc e che probabilmente non e` un vero segnale,</span> |
| 236 |
<a name="l00217"></a>00217 <span class="comment">// ma e` il pile-up dell'evento precedente</span> |
| 237 |
<a name="l00218"></a>00218 |
| 238 |
<a name="l00219"></a>00219 <span class="comment">// loop over TDC hits along the track</span> |
| 239 |
<a name="l00220"></a>00220 <span class="keywordtype">int</span> iht = -1; |
| 240 |
<a name="l00221"></a>00221 <span class="keywordflow">for</span> (iht = 0; iht < tof->npmttdc; iht++) |
| 241 |
<a name="l00222"></a>00222 <span class="keywordflow">if</span> (<span class="keywordtype">id</span> == tof->pmttdc[iht]) |
| 242 |
<a name="l00223"></a>00223 <span class="keywordflow">break</span>; <span class="comment">//match</span> |
| 243 |
<a name="l00224"></a>00224 <span class="keywordflow">if</span> (iht >= 0 && iht < tof->npmttdc) |
| 244 |
<a name="l00225"></a>00225 <span class="keywordflow">continue</span>; |
| 245 |
<a name="l00226"></a>00226 <span class="comment">// loop over ADC hits along the track (credo sia inutile...)</span> |
| 246 |
<a name="l00227"></a>00227 <span class="comment">// int iha = -1;</span> |
| 247 |
<a name="l00228"></a>00228 <span class="comment">// for(iha=0; iha<tof->npmtadc; iha++)if(id==tof->pmtadc[iha])break; //match</span> |
| 248 |
<a name="l00229"></a>00229 <span class="comment">// if( iha>=0 && iha <tof->npmtadc )continue;</span> |
| 249 |
<a name="l00230"></a>00230 <span class="comment">// get here if outside track</span> |
| 250 |
<a name="l00231"></a>00231 nn++; |
| 251 |
<a name="l00232"></a>00232 |
| 252 |
<a name="l00233"></a>00233 <span class="comment">// cout << " OUTSIDE ";</span> |
| 253 |
<a name="l00234"></a>00234 } |
| 254 |
<a name="l00235"></a>00235 <span class="keywordflow">return</span> nn; |
| 255 |
<a name="l00236"></a>00236 |
| 256 |
<a name="l00237"></a>00237 } |
| 257 |
<a name="l00238"></a>00238 |
| 258 |
</pre></div></div> |
| 259 |
<hr size="1"><address style="text-align: right;"><small>Generated on Wed May 27 15:21:27 2009 for PamCut by |
| 260 |
<a href="http://www.doxygen.org/index.html"> |
| 261 |
<img src="doxygen.png" alt="doxygen" align="middle" border="0"></a> 1.5.9 </small></address> |
| 262 |
</body> |
| 263 |
</html> |
Parent Directory
| 
Revision Log