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 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/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 < 9; z++) { |
76 |
<a name="l00057"></a>00057 <span class="keywordflow">for</span> (<span class="keywordtype">int</span> i = 0; i < 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 < 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 < 1 || _Z > 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>->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>->GetTrack(0)->GetTrkTrack()->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<<span class="keyword">event</span>->GetToFLevel2()->ntrk(); i++) { |
114 |
<a name="l00095"></a>00095 <span class="keywordflow">if</span> ( event->GetToFLevel2()->GetToFTrkVar(i)->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>->GetToFLevel2()->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 < 6; i++) { |
121 |
<a name="l00102"></a>00102 h_tof_z_beta[i]->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 < 6; i++) { <span class="comment">// determine good, non-masked layers</span> |
126 |
<a name="l00107"></a>00107 <span class="keywordflow">if</span> ((charge[i] < 1000) && (charge[i] > 0) && ((<a class="code" href="classTofNucleiZCut.html#bd6900afd20211e94489644410e1367f">_maskArray</a>[i] & <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 < <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 < 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] > <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] < <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 < 6; i++) { |
150 |
<a name="l00131"></a>00131 h_tof_z_beta[6+i]->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 |
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> |