| 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/TofBetaRigCut/TofBetaRigCut.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/TofBetaRigCut/TofBetaRigCut.cpp</h1><a href="TofBetaRigCut_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"> * TofBetaRigCut.cpp</span> |
| 24 |
<a name="l00003"></a>00003 <span class="comment"> *</span> |
| 25 |
<a name="l00004"></a>00004 <span class="comment"> * Created on: 13-mar-2009</span> |
| 26 |
<a name="l00005"></a>00005 <span class="comment"> * Author: Sergio Ricciarini, Nicola Mori</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="TofBetaRigCut_8h.html">TofBetaRigCut.h</a>"</span> |
| 30 |
<a name="l00011"></a>00011 |
| 31 |
<a name="l00012"></a><a class="code" href="classTofBetaRigCut.html#661c3a16f981d43c28ccea04242a2a17">00012</a> <span class="keywordtype">int</span> <a class="code" href="classTofBetaRigCut.html#661c3a16f981d43c28ccea04242a2a17" title="The rigidity check.">TofBetaRigCut::Check</a>(PamLevel2 *event) { |
| 32 |
<a name="l00013"></a>00013 |
| 33 |
<a name="l00014"></a>00014 TrkTrack *trkTrack = <span class="keyword">event</span>->GetTrack(0)->GetTrkTrack(); |
| 34 |
<a name="l00015"></a>00015 |
| 35 |
<a name="l00016"></a>00016 Int_t seqno = trkTrack->GetSeqNo(); |
| 36 |
<a name="l00017"></a>00017 Int_t ntrack_tof = 0; |
| 37 |
<a name="l00018"></a>00018 <span class="keywordflow">for</span> (Int_t i = 0; i < <span class="keyword">event</span>->GetToFLevel2()->ntrk(); i++) { |
| 38 |
<a name="l00019"></a>00019 <span class="keywordflow">if</span> (event->GetToFLevel2()->GetToFTrkVar(i)->trkseqno == seqno) |
| 39 |
<a name="l00020"></a>00020 ntrack_tof = i; |
| 40 |
<a name="l00021"></a>00021 } |
| 41 |
<a name="l00022"></a>00022 Float_t betaModTof = fabs(event->GetToFLevel2()->CalcBeta(ntrack_tof, 10., 10., 20.)); <span class="comment">// beta[12]</span> |
| 42 |
<a name="l00023"></a>00023 |
| 43 |
<a name="l00024"></a>00024 <span class="comment">// Float_t rigMod = _mass*fabs(beta)*pow((1-pow(beta,2)),-0.5)/abs(_charge);</span> |
| 44 |
<a name="l00025"></a>00025 |
| 45 |
<a name="l00026"></a>00026 Float_t betaModTrk = (trkTrack->GetRigidity() * abs(<a class="code" href="classTofBetaRigCut.html#f96c5c502466dce216acef86b40f90c5">_charge</a>)) / pow(pow(trkTrack->GetRigidity() * <a class="code" href="classTofBetaRigCut.html#f96c5c502466dce216acef86b40f90c5">_charge</a>, 2) + pow( |
| 46 |
<a name="l00027"></a>00027 <a class="code" href="classTofBetaRigCut.html#6f5740800a492eda945ee8be637b70aa">_mass</a>, 2), 0.5); |
| 47 |
<a name="l00028"></a>00028 |
| 48 |
<a name="l00029"></a>00029 <span class="preprocessor">#ifdef DEBUGPAMCUT</span> |
| 49 |
<a name="l00030"></a>00030 <span class="preprocessor"></span> |
| 50 |
<a name="l00031"></a>00031 h_trk_vs_tof_betamod[0]->Fill(betaModTof,betaModTrk); |
| 51 |
<a name="l00032"></a>00032 |
| 52 |
<a name="l00033"></a>00033 <span class="preprocessor">#endif</span> |
| 53 |
<a name="l00034"></a>00034 <span class="preprocessor"></span> |
| 54 |
<a name="l00035"></a>00035 <span class="keywordflow">if</span> (fabs(betaModTof - betaModTrk) > <a class="code" href="classTofBetaRigCut.html#64d4ca7f488b148a47ec57d63a21bd7a">_threshold</a> * betaModTof) { |
| 55 |
<a name="l00036"></a>00036 <span class="keywordflow">return</span> 0; |
| 56 |
<a name="l00037"></a>00037 } |
| 57 |
<a name="l00038"></a>00038 |
| 58 |
<a name="l00039"></a>00039 <span class="preprocessor">#ifdef DEBUGPAMCUT</span> |
| 59 |
<a name="l00040"></a>00040 <span class="preprocessor"></span> |
| 60 |
<a name="l00041"></a>00041 h_trk_vs_tof_betamod[1]->Fill(betaModTof,betaModTrk); |
| 61 |
<a name="l00042"></a>00042 |
| 62 |
<a name="l00043"></a>00043 <span class="preprocessor">#endif</span> |
| 63 |
<a name="l00044"></a>00044 <span class="preprocessor"></span> |
| 64 |
<a name="l00045"></a>00045 <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>; |
| 65 |
<a name="l00046"></a>00046 } |
| 66 |
</pre></div></div> |
| 67 |
<hr size="1"><address style="text-align: right;"><small>Generated on Wed May 27 15:21:27 2009 for PamCut by |
| 68 |
<a href="http://www.doxygen.org/index.html"> |
| 69 |
<img src="doxygen.png" alt="doxygen" align="middle" border="0"></a> 1.5.9 </small></address> |
| 70 |
</body> |
| 71 |
</html> |
Parent Directory
| 
Revision Log