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

Contents of /PamCut/doc/Doxygen/html/TofNucleiZCut_8cpp_source.html

Parent Directory Parent Directory | Revision Log Revision Log


Revision 1.2 - (show annotations) (download) (as text)
Fri Sep 25 15:41:25 2009 UTC (15 years, 6 months ago) by pam-fi
Branch: MAIN
CVS Tags: HEAD
Changes since 1.1: +0 -0 lines
File MIME type: text/html
FILE REMOVED
Removed from repository.

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/TofNucleiZCut/TofNucleiZCut.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&nbsp;Page</span></a></li>
12 <li><a href="annotated.html"><span>Classes</span></a></li>
13 <li class="current"><a href="files.html"><span>Files</span></a></li>
14 </ul>
15 </div>
16 <div class="tabs">
17 <ul>
18 <li><a href="files.html"><span>File&nbsp;List</span></a></li>
19 <li><a href="globals.html"><span>File&nbsp;Members</span></a></li>
20 </ul>
21 </div>
22 <h1>TofCuts/TofNucleiZCut/TofNucleiZCut.cpp</h1><a href="TofNucleiZCut_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"> * TofNucleiZCut.cpp</span>
24 <a name="l00003"></a>00003 <span class="comment"> *</span>
25 <a name="l00004"></a>00004 <span class="comment"> * Created on: 21-apr-2009</span>
26 <a name="l00005"></a>00005 <span class="comment"> * Author: Nicola Mori</span>
27 <a name="l00006"></a>00006 <span class="comment"> */</span>
28 <a name="l00007"></a>00007
29 <a name="l00008"></a>00008 <span class="preprocessor">#ifndef NO_TOFNUCLEI</span>
30 <a name="l00009"></a>00009 <span class="preprocessor"></span>
31 <a name="l00012"></a>00012 <span class="preprocessor">#include "<a class="code" href="TofNucleiZCut_8h.html">TofNucleiZCut.h</a>"</span>
32 <a name="l00013"></a>00013
33 <a name="l00014"></a>00014 <span class="keyword">const</span> <span class="keywordtype">float</span> <a class="code" href="classTofNucleiZCut.html#e9e277a7761ca24ec76c5ba04fe93ed5">TofNucleiZCut::_sigmaZdefault</a>[] = { 0.12, 0.15, 0.20, 0.25, 0.32, 0.40, 0.45, 0.50 }; <span class="comment">// the charge-width in a layer for Z=1..8</span>
34 <a name="l00015"></a>00015 <span class="keyword">const</span> <span class="keywordtype">unsigned</span> <span class="keywordtype">int</span> <a class="code" href="classTofNucleiZCut.html#bd6900afd20211e94489644410e1367f">TofNucleiZCut::_maskArray</a>[] = { <a class="code" href="TofNucleiZCut_8h.html#2d58b1b799df7d79374ed3ab6cbd204a9f207c4c51713cc11b8a8040caa150c5" title="S11.">TOFNUCLEIZ_S11</a>, <a class="code" href="TofNucleiZCut_8h.html#2d58b1b799df7d79374ed3ab6cbd204a51c80f879c6c303b029c1198bd36dfa9" title="S12.">TOFNUCLEIZ_S12</a>, <a class="code" href="TofNucleiZCut_8h.html#2d58b1b799df7d79374ed3ab6cbd204a7502fcae3146992720651315bc8c477f" title="S21.">TOFNUCLEIZ_S21</a>, <a class="code" href="TofNucleiZCut_8h.html#2d58b1b799df7d79374ed3ab6cbd204aa19e694ccfcf5d719221b3c182db460f" title="S22.">TOFNUCLEIZ_S22</a>,
35 <a name="l00016"></a>00016 <a class="code" href="TofNucleiZCut_8h.html#2d58b1b799df7d79374ed3ab6cbd204a5f2d7e3b31610f8c9265aac5603803a5" title="S31.">TOFNUCLEIZ_S31</a>, <a class="code" href="TofNucleiZCut_8h.html#2d58b1b799df7d79374ed3ab6cbd204a6009fe0dd60a524c817c0ce5dfc8aeeb" title="S32.">TOFNUCLEIZ_S32</a> }; <span class="comment">//The layer - mask code association</span>
36 <a name="l00017"></a>00017
37 <a name="l00018"></a>00018
38 <a name="l00019"></a><a class="code" href="classTofNucleiZCut.html#0da6346342abf76e1dd481360ccccde8">00019</a> <a class="code" href="classTofNucleiZCut.html#0da6346342abf76e1dd481360ccccde8" title="Constructor.">TofNucleiZCut::TofNucleiZCut</a>(<span class="keyword">const</span> <span class="keywordtype">char</span> *cutName, <span class="keywordtype">unsigned</span> <span class="keywordtype">int</span> Z, <span class="keywordtype">float</span> lowerLimit, <span class="keywordtype">float</span> upperLimit, <span class="keywordtype">unsigned</span> <span class="keywordtype">int</span> minLayers,
39 <a name="l00020"></a>00020 <span class="keywordtype">unsigned</span> <span class="keywordtype">int</span> layersMask) :
40 <a name="l00021"></a>00021 <a class="code" href="classPamCut.html" title="An abstract class to apply cuts to Pamela data.">PamCut</a>(cutName), _Z(Z), _mean(8, 6), _sigma(8, 6), _lowerLimit(lowerLimit), _upperLimit(upperLimit), _minLayers(
41 <a name="l00022"></a>00022 minLayers), _layersMask(layersMask), _tofNuclei(NULL) {
42 <a name="l00023"></a>00023
43 <a name="l00024"></a>00024 <span class="comment">//TODO Put calibration parameters on a file; load only those parameters needed for the current Z.</span>
44 <a name="l00025"></a>00025
45 <a name="l00026"></a>00026 <span class="comment">// Hard-coded calibrations</span>
46 <a name="l00027"></a>00027 <span class="comment">// H</span>
47 <a name="l00028"></a>00028 <a class="code" href="classTofNucleiZCut.html#cf873efeadca326eedc7fe647c44bd89">_mean</a>[0][0] = 0.987; <span class="comment">// S11 peak</span>
48 <a name="l00029"></a>00029 <a class="code" href="classTofNucleiZCut.html#eea4796a4958fef682803939034579cf">_sigma</a>[0][0] = 0.11; <span class="comment">// S11 sigma</span>
49 <a name="l00030"></a>00030 <a class="code" href="classTofNucleiZCut.html#cf873efeadca326eedc7fe647c44bd89">_mean</a>[0][1] = 0.988; <span class="comment">// S12 peak</span>
50 <a name="l00031"></a>00031 <a class="code" href="classTofNucleiZCut.html#eea4796a4958fef682803939034579cf">_sigma</a>[0][1] = 0.12; <span class="comment">// S12 sigma</span>
51 <a name="l00032"></a>00032 <a class="code" href="classTofNucleiZCut.html#cf873efeadca326eedc7fe647c44bd89">_mean</a>[0][2] = 0.961; <span class="comment">// S21 peak</span>
52 <a name="l00033"></a>00033 <a class="code" href="classTofNucleiZCut.html#eea4796a4958fef682803939034579cf">_sigma</a>[0][2] = 0.083;<span class="comment">// S21 sigma</span>
53 <a name="l00034"></a>00034 <a class="code" href="classTofNucleiZCut.html#cf873efeadca326eedc7fe647c44bd89">_mean</a>[0][3] = 0.966; <span class="comment">// S22 peak</span>
54 <a name="l00035"></a>00035 <a class="code" href="classTofNucleiZCut.html#eea4796a4958fef682803939034579cf">_sigma</a>[0][3] = 0.10; <span class="comment">// S22 sigma</span>
55 <a name="l00036"></a>00036 <a class="code" href="classTofNucleiZCut.html#cf873efeadca326eedc7fe647c44bd89">_mean</a>[0][4] = 0.981; <span class="comment">// S31 peak</span>
56 <a name="l00037"></a>00037 <a class="code" href="classTofNucleiZCut.html#eea4796a4958fef682803939034579cf">_sigma</a>[0][4] = 0.092;<span class="comment">// S31 sigma</span>
57 <a name="l00038"></a>00038 <a class="code" href="classTofNucleiZCut.html#cf873efeadca326eedc7fe647c44bd89">_mean</a>[0][5] = 0.979; <span class="comment">// S32 peak</span>
58 <a name="l00039"></a>00039 <a class="code" href="classTofNucleiZCut.html#eea4796a4958fef682803939034579cf">_sigma</a>[0][5] = 0.095;<span class="comment">// S32 sigma</span>
59 <a name="l00040"></a>00040
60 <a name="l00041"></a>00041 <span class="comment">// He</span>
61 <a name="l00042"></a>00042 <a class="code" href="classTofNucleiZCut.html#cf873efeadca326eedc7fe647c44bd89">_mean</a>[1][0] = 1.96; <span class="comment">// S11 peak</span>
62 <a name="l00043"></a>00043 <a class="code" href="classTofNucleiZCut.html#eea4796a4958fef682803939034579cf">_sigma</a>[1][0] = 0.12;<span class="comment">// S11 sigma</span>
63 <a name="l00044"></a>00044 <a class="code" href="classTofNucleiZCut.html#cf873efeadca326eedc7fe647c44bd89">_mean</a>[1][1] = 1.95; <span class="comment">// S12 peak</span>
64 <a name="l00045"></a>00045 <a class="code" href="classTofNucleiZCut.html#eea4796a4958fef682803939034579cf">_sigma</a>[1][1] = 0.14;<span class="comment">// S12 sigma</span>
65 <a name="l00046"></a>00046 <a class="code" href="classTofNucleiZCut.html#cf873efeadca326eedc7fe647c44bd89">_mean</a>[1][2] = 1.95; <span class="comment">// S21 peak</span>
66 <a name="l00047"></a>00047 <a class="code" href="classTofNucleiZCut.html#eea4796a4958fef682803939034579cf">_sigma</a>[1][2] = 0.13;<span class="comment">// S21 sigma</span>
67 <a name="l00048"></a>00048 <a class="code" href="classTofNucleiZCut.html#cf873efeadca326eedc7fe647c44bd89">_mean</a>[1][3] = 1.96; <span class="comment">// S22 peak</span>
68 <a name="l00049"></a>00049 <a class="code" href="classTofNucleiZCut.html#eea4796a4958fef682803939034579cf">_sigma</a>[1][3] = 0.14;<span class="comment">// S22 sigma</span>
69 <a name="l00050"></a>00050 <a class="code" href="classTofNucleiZCut.html#cf873efeadca326eedc7fe647c44bd89">_mean</a>[1][4] = 1.99; <span class="comment">// S31 peak</span>
70 <a name="l00051"></a>00051 <a class="code" href="classTofNucleiZCut.html#eea4796a4958fef682803939034579cf">_sigma</a>[1][4] = 0.14;<span class="comment">// S31 sigma</span>
71 <a name="l00052"></a>00052 <a class="code" href="classTofNucleiZCut.html#cf873efeadca326eedc7fe647c44bd89">_mean</a>[1][5] = 2.00; <span class="comment">// S32 peak</span>
72 <a name="l00053"></a>00053 <a class="code" href="classTofNucleiZCut.html#eea4796a4958fef682803939034579cf">_sigma</a>[1][5] = 0.15;<span class="comment">// S32 sigma</span>
73 <a name="l00054"></a>00054
74 <a name="l00055"></a>00055 <span class="comment">// Heavier nuclei are not calibrated yet; using default values</span>
75 <a name="l00056"></a>00056 <span class="keywordflow">for</span> (<span class="keywordtype">unsigned</span> <span class="keywordtype">int</span> z = 3; z &lt; 9; z++) {
76 <a name="l00057"></a>00057 <span class="keywordflow">for</span> (<span class="keywordtype">int</span> i = 0; i &lt; 6; i++) {
77 <a name="l00058"></a>00058 <a class="code" href="classTofNucleiZCut.html#cf873efeadca326eedc7fe647c44bd89">_mean</a>[z - 1][i] = z;
78 <a name="l00059"></a>00059 <a class="code" href="classTofNucleiZCut.html#eea4796a4958fef682803939034579cf">_sigma</a>[z - 1][i] = <a class="code" href="classTofNucleiZCut.html#e9e277a7761ca24ec76c5ba04fe93ed5">_sigmaZdefault</a>[z - 1];
79 <a name="l00060"></a>00060 }
80 <a name="l00061"></a>00061 }
81 <a name="l00062"></a>00062
82 <a name="l00063"></a>00063 <span class="preprocessor">#ifdef DEBUGPAMCUT</span>
83 <a name="l00064"></a>00064 <span class="preprocessor"></span>
84 <a name="l00065"></a>00065 TString hId;
85 <a name="l00066"></a>00066 TString hTitle;
86 <a name="l00067"></a>00067
87 <a name="l00068"></a>00068 <span class="keywordflow">for</span> (UInt_t j = 0; j &lt; 12; j++) {
88 <a name="l00069"></a>00069 hId.Form(<span class="stringliteral">"h_tof_z_vs_beta_%i"</span>, j);
89 <a name="l00070"></a>00070 hTitle.Form(<span class="stringliteral">"TOF Z vs beta (%i)"</span>, j);
90 <a name="l00071"></a>00071 h_tof_z_beta[j] = <span class="keyword">new</span> TH2F(hId.Data(), hTitle.Data(), 50, 0, 1, 50, 0, 10);
91 <a name="l00072"></a>00072 }
92 <a name="l00073"></a>00073
93 <a name="l00074"></a>00074 <span class="preprocessor">#endif</span>
94 <a name="l00075"></a>00075 <span class="preprocessor"></span>
95 <a name="l00076"></a>00076 }
96 <a name="l00077"></a>00077
97 <a name="l00078"></a><a class="code" href="classTofNucleiZCut.html#86aa8211e4afee39ff38080347bf651d">00078</a> <span class="keywordtype">int</span> <a class="code" href="classTofNucleiZCut.html#86aa8211e4afee39ff38080347bf651d" title="The TofNucleiZ check.">TofNucleiZCut::Check</a>(PamLevel2 *event) {
98 <a name="l00079"></a>00079
99 <a name="l00080"></a>00080 <span class="keywordflow">if</span> (_Z &lt; 1 || _Z &gt; 8)
100 <a name="l00081"></a>00081 <span class="keywordflow">return</span> <a class="code" href="TofNucleiZCut_8h.html#7cbf6ce6906f1524de86f3d84f23aa54f09f54efe798d0a77429385071fc64fe" title="Discarded because charge is not between 1 and 8.">TOFNUCLEIZ_ILLEGALZ</a>;
101 <a name="l00082"></a>00082
102 <a name="l00083"></a>00083 <span class="comment">// Check if ToFNuclei has already been initialized</span>
103 <a name="l00084"></a>00084 <span class="keywordflow">if</span> (<a class="code" href="classTofNucleiZCut.html#6e52d79381364353da0780db39b98492">_tofNuclei</a> == NULL)
104 <a name="l00085"></a>00085 <a class="code" href="classTofNucleiZCut.html#6e52d79381364353da0780db39b98492">_tofNuclei</a> = <span class="keyword">new</span> ToFNuclei(event);
105 <a name="l00086"></a>00086
106 <a name="l00087"></a>00087 Float_t *charge = <a class="code" href="classTofNucleiZCut.html#6e52d79381364353da0780db39b98492">_tofNuclei</a>-&gt;Get_Charge_ToF_std_layer(); <span class="comment">//Returns a 6 elements array</span>
107 <a name="l00088"></a>00088 <span class="keywordtype">unsigned</span> <span class="keywordtype">int</span> nGood = 0;
108 <a name="l00089"></a>00089 <span class="keywordtype">bool</span> good[6];
109 <a name="l00090"></a>00090
110 <a name="l00091"></a>00091 <span class="preprocessor">#ifdef DEBUGPAMCUT</span>
111 <a name="l00092"></a>00092 <span class="preprocessor"></span> Int_t seqno = <span class="keyword">event</span>-&gt;GetTrack(0)-&gt;GetTrkTrack()-&gt;GetSeqNo();
112 <a name="l00093"></a>00093 Int_t ntrack_tof = 0;
113 <a name="l00094"></a>00094 <span class="keywordflow">for</span> (Int_t i=0; i&lt;<span class="keyword">event</span>-&gt;GetToFLevel2()-&gt;ntrk(); i++) {
114 <a name="l00095"></a>00095 <span class="keywordflow">if</span> ( event-&gt;GetToFLevel2()-&gt;GetToFTrkVar(i)-&gt;trkseqno == seqno ) ntrack_tof = i;
115 <a name="l00096"></a>00096 }
116 <a name="l00097"></a>00097 Float_t beta=<span class="keyword">event</span>-&gt;GetToFLevel2()-&gt;CalcBeta(ntrack_tof,10.,10.,20.); <span class="comment">// beta[12]</span>
117 <a name="l00098"></a>00098 <span class="preprocessor">#endif</span>
118 <a name="l00099"></a>00099 <span class="preprocessor"></span>
119 <a name="l00100"></a>00100 <span class="preprocessor">#ifdef DEBUGPAMCUT</span>
120 <a name="l00101"></a>00101 <span class="preprocessor"></span> <span class="keywordflow">for</span> (<span class="keywordtype">int</span> i = 0; i &lt; 6; i++) {
121 <a name="l00102"></a>00102 h_tof_z_beta[i]-&gt;Fill(beta, charge[i]);
122 <a name="l00103"></a>00103 }
123 <a name="l00104"></a>00104 <span class="preprocessor">#endif</span>
124 <a name="l00105"></a>00105 <span class="preprocessor"></span>
125 <a name="l00106"></a>00106 <span class="keywordflow">for</span> (<span class="keywordtype">int</span> i = 0; i &lt; 6; i++) { <span class="comment">// determine good, non-masked layers</span>
126 <a name="l00107"></a>00107 <span class="keywordflow">if</span> ((charge[i] &lt; 1000) &amp;&amp; (charge[i] &gt; 0) &amp;&amp; ((<a class="code" href="classTofNucleiZCut.html#bd6900afd20211e94489644410e1367f">_maskArray</a>[i] &amp; <a class="code" href="classTofNucleiZCut.html#59b391889d5a9e2c4a127adb92f0277e">_layersMask</a>) != <a class="code" href="classTofNucleiZCut.html#bd6900afd20211e94489644410e1367f">_maskArray</a>[i])) {
127 <a name="l00108"></a>00108 good[i] = <span class="keyword">true</span>;
128 <a name="l00109"></a>00109 nGood++;
129 <a name="l00110"></a>00110 }
130 <a name="l00111"></a>00111 <span class="keywordflow">else</span> {
131 <a name="l00112"></a>00112 good[i] = <span class="keyword">false</span>;
132 <a name="l00113"></a>00113 }
133 <a name="l00114"></a>00114 }
134 <a name="l00115"></a>00115
135 <a name="l00116"></a>00116 <span class="keywordflow">if</span> (nGood &lt; <a class="code" href="classTofNucleiZCut.html#4c74f3475f72039dea6c87ec92276c1f">_minLayers</a>) {
136 <a name="l00117"></a>00117 <span class="keywordflow">return</span> <a class="code" href="TofNucleiZCut_8h.html#7cbf6ce6906f1524de86f3d84f23aa54179c0935201f5d6ac6d7f00ebd1e959c" title="Discarded because charge can be measured in too few layers.">TOFNUCLEIZ_TOOFEWLAYERS</a>;
137 <a name="l00118"></a>00118 }
138 <a name="l00119"></a>00119
139 <a name="l00120"></a>00120 <span class="keywordflow">for</span> (<span class="keywordtype">int</span> i = 0; i &lt; 6; i++) { <span class="comment">// check done only for good, non-masked layers</span>
140 <a name="l00121"></a>00121 <span class="keywordflow">if</span> (good[i]) {
141 <a name="l00122"></a>00122 <span class="keywordflow">if</span> ((charge[i] &gt; <a class="code" href="classTofNucleiZCut.html#cf873efeadca326eedc7fe647c44bd89">_mean</a>[<a class="code" href="classTofNucleiZCut.html#356fa811a65d5194dbbc4651d83530ef">_Z</a> - 1][i] + <a class="code" href="classTofNucleiZCut.html#884c0edaccd86f4b121572386c42a290">_upperLimit</a> * <a class="code" href="classTofNucleiZCut.html#eea4796a4958fef682803939034579cf">_sigma</a>[<a class="code" href="classTofNucleiZCut.html#356fa811a65d5194dbbc4651d83530ef">_Z</a> - 1][i]) || ((charge[i] &lt; <a class="code" href="classTofNucleiZCut.html#cf873efeadca326eedc7fe647c44bd89">_mean</a>[<a class="code" href="classTofNucleiZCut.html#356fa811a65d5194dbbc4651d83530ef">_Z</a> - 1][i]
142 <a name="l00123"></a>00123 - <a class="code" href="classTofNucleiZCut.html#d4436c3e8efcd48f8befda3c6395f4ed">_lowerLimit</a> * <a class="code" href="classTofNucleiZCut.html#eea4796a4958fef682803939034579cf">_sigma</a>[<a class="code" href="classTofNucleiZCut.html#356fa811a65d5194dbbc4651d83530ef">_Z</a> - 1][i]))) {
143 <a name="l00124"></a>00124 <span class="keywordflow">return</span> <a class="code" href="TofNucleiZCut_8h.html#7cbf6ce6906f1524de86f3d84f23aa54ef9ad5021a1b83dd6a53eb225044d7c6" title="Discarded because charge is out of bounds in a layer.">TOFNUCLEIZ_OUTOFBOUNDS</a>;
144 <a name="l00125"></a>00125 }
145 <a name="l00126"></a>00126 }
146 <a name="l00127"></a>00127 }
147 <a name="l00128"></a>00128
148 <a name="l00129"></a>00129 <span class="preprocessor">#ifdef DEBUGPAMCUT</span>
149 <a name="l00130"></a>00130 <span class="preprocessor"></span> <span class="keywordflow">for</span> (<span class="keywordtype">int</span> i = 0; i &lt; 6; i++) {
150 <a name="l00131"></a>00131 h_tof_z_beta[6+i]-&gt;Fill(beta, charge[i]);
151 <a name="l00132"></a>00132 }
152 <a name="l00133"></a>00133 <span class="preprocessor">#endif</span>
153 <a name="l00134"></a>00134 <span class="preprocessor"></span>
154 <a name="l00135"></a>00135 <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>;
155 <a name="l00136"></a>00136
156 <a name="l00137"></a>00137 }
157 <a name="l00138"></a>00138
158 <a name="l00139"></a>00139 <span class="preprocessor">#endif </span><span class="comment">/* NO_TOFNUCLEI */</span>
159 </pre></div></div>
160 <hr size="1"><address style="text-align: right;"><small>Generated on Wed May 27 15:21:27 2009 for PamCut by&nbsp;
161 <a href="http://www.doxygen.org/index.html">
162 <img src="doxygen.png" alt="doxygen" align="middle" border="0"></a> 1.5.9 </small></address>
163 </body>
164 </html>

  ViewVC Help
Powered by ViewVC 1.1.23