| 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: TrkCuts/RigFillCut/RigFillCut.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>TrkCuts/RigFillCut/RigFillCut.cpp</h1><a href="RigFillCut_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"> * RigFillCut.cpp</span> |
| 24 |
<a name="l00003"></a>00003 <span class="comment"> *</span> |
| 25 |
<a name="l00004"></a>00004 <span class="comment"> * Created on: 27-mar-2009</span> |
| 26 |
<a name="l00005"></a>00005 <span class="comment"> * Author: Nicola Mori, S. Ricciarini</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="RigFillCut_8h.html">RigFillCut.h</a>"</span> |
| 30 |
<a name="l00011"></a>00011 |
| 31 |
<a name="l00012"></a><a class="code" href="classRigFillCut.html#5441206a7a24f1f5d19dc1a64685fcd4">00012</a> <span class="keywordtype">int</span> <a class="code" href="classRigFillCut.html#5441206a7a24f1f5d19dc1a64685fcd4" title="The rigidity and cutoff rigidity check.">RigFillCut::Check</a>(PamLevel2 *event) { |
| 32 |
<a name="l00013"></a>00013 |
| 33 |
<a name="l00014"></a>00014 <span class="keywordtype">float</span> rigThreshold = <a class="code" href="classRigFillCut.html#5634e95eb75291fd866f995a324a3a0a">_thresholdCoeff</a> * <span class="keyword">event</span>->GetOrbitalInfo()->GetCutoffSVL(); |
| 34 |
<a name="l00015"></a>00015 |
| 35 |
<a name="l00016"></a>00016 <span class="keywordflow">if</span> (rigThreshold > <a class="code" href="classRigFillCut.html#24c7fece1b02c4222dc2fd9d0d13bcc0">_binning</a>[<a class="code" href="classRigFillCut.html#24c7fece1b02c4222dc2fd9d0d13bcc0">_binning</a>.size() - 1] || rigThreshold < <a class="code" href="classRigFillCut.html#24c7fece1b02c4222dc2fd9d0d13bcc0">_binning</a>[0]) { |
| 36 |
<a name="l00017"></a>00017 <span class="keywordflow">return</span> <a class="code" href="RigFillCut_8h.html#7ea9d9857d7f7dc01a217d60fc93ebc42732da7410302fac27f6e86017ece272" title="Event is discarded because its threshold rigidity is outside the threshold rigidity...">THRESHRIGOUT</a>; |
| 37 |
<a name="l00018"></a>00018 } |
| 38 |
<a name="l00019"></a>00019 |
| 39 |
<a name="l00020"></a>00020 <span class="keywordtype">float</span> rig = <span class="keyword">event</span>->GetTrack(0)->GetTrkTrack()->GetRigidity(); |
| 40 |
<a name="l00021"></a>00021 |
| 41 |
<a name="l00022"></a>00022 <span class="keywordflow">if</span> (rig > <a class="code" href="classRigFillCut.html#24c7fece1b02c4222dc2fd9d0d13bcc0">_binning</a>[<a class="code" href="classRigFillCut.html#24c7fece1b02c4222dc2fd9d0d13bcc0">_binning</a>.size() - 1] || rig < <a class="code" href="classRigFillCut.html#24c7fece1b02c4222dc2fd9d0d13bcc0">_binning</a>[0]) { |
| 42 |
<a name="l00023"></a>00023 <span class="keywordflow">return</span> <a class="code" href="RigFillCut_8h.html#7ea9d9857d7f7dc01a217d60fc93ebc4277656d13d84afccc03957b79674e573" title="Event is discarded because its rigidity is outside the rigidity binning limits.">RIGOUT</a>; |
| 43 |
<a name="l00024"></a>00024 } |
| 44 |
<a name="l00025"></a>00025 |
| 45 |
<a name="l00026"></a>00026 <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>; |
| 46 |
<a name="l00027"></a>00027 } |
| 47 |
<a name="l00028"></a>00028 |
| 48 |
<a name="l00029"></a><a class="code" href="classRigFillCut.html#537f1b31f208ebe8f190f254be1fdc6c">00029</a> <span class="keywordtype">void</span> <a class="code" href="classRigFillCut.html#537f1b31f208ebe8f190f254be1fdc6c" title="The histogram filling.">RigFillCut::OnGood</a>(PamLevel2 *event) { |
| 49 |
<a name="l00030"></a>00030 |
| 50 |
<a name="l00031"></a>00031 <span class="keywordtype">float</span> rigThreshold = <a class="code" href="classRigFillCut.html#5634e95eb75291fd866f995a324a3a0a">_thresholdCoeff</a> * <span class="keyword">event</span>->GetOrbitalInfo()->GetCutoffSVL(); |
| 51 |
<a name="l00032"></a>00032 <span class="keywordtype">float</span> rig = <span class="keyword">event</span>->GetTrack(0)->GetTrkTrack()->GetRigidity(); |
| 52 |
<a name="l00033"></a>00033 |
| 53 |
<a name="l00034"></a>00034 <span class="comment">//Bin identification: value must be less than (<) bin maximum</span> |
| 54 |
<a name="l00035"></a>00035 |
| 55 |
<a name="l00036"></a>00036 <span class="keywordtype">int</span> i = 1; |
| 56 |
<a name="l00037"></a>00037 <span class="keywordflow">while</span> (rig >= <a class="code" href="classRigFillCut.html#24c7fece1b02c4222dc2fd9d0d13bcc0">_binning</a>[i]) { |
| 57 |
<a name="l00038"></a>00038 i++; |
| 58 |
<a name="l00039"></a>00039 } |
| 59 |
<a name="l00040"></a>00040 i--; |
| 60 |
<a name="l00041"></a>00041 |
| 61 |
<a name="l00042"></a>00042 <span class="comment">// cout << _binning[i] << " < " << rig << " < " << _binning[i+1] << endl;</span> |
| 62 |
<a name="l00043"></a>00043 |
| 63 |
<a name="l00044"></a>00044 <span class="keywordtype">int</span> j = 1; |
| 64 |
<a name="l00045"></a>00045 <span class="keywordflow">while</span> (rigThreshold >= <a class="code" href="classRigFillCut.html#24c7fece1b02c4222dc2fd9d0d13bcc0">_binning</a>[j]) { |
| 65 |
<a name="l00046"></a>00046 j++; |
| 66 |
<a name="l00047"></a>00047 } |
| 67 |
<a name="l00048"></a>00048 j--; |
| 68 |
<a name="l00049"></a>00049 |
| 69 |
<a name="l00050"></a>00050 <span class="comment">// cout << _binning[j] << " < " << rigCutoff << " < " << _binning[j+1] << endl;</span> |
| 70 |
<a name="l00051"></a>00051 <span class="comment">// cout << endl;</span> |
| 71 |
<a name="l00052"></a>00052 |
| 72 |
<a name="l00053"></a>00053 <a class="code" href="classRigFillCut.html#9c0b18d019325388ab4d99cd304a5909">_histogram</a>[i][j]++; |
| 73 |
<a name="l00054"></a>00054 } |
| 74 |
</pre></div></div> |
| 75 |
<hr size="1"><address style="text-align: right;"><small>Generated on Wed May 27 15:21:27 2009 for PamCut by |
| 76 |
<a href="http://www.doxygen.org/index.html"> |
| 77 |
<img src="doxygen.png" alt="doxygen" align="middle" border="0"></a> 1.5.9 </small></address> |
| 78 |
</body> |
| 79 |
</html> |
Parent Directory
| 
Revision Log