HLT/PHOS/html/AliHLTPHOSPulseGenerator_8cxx-source.html

   1 <!DOCTYPE HTML PUBLIC "-//W3C//DTD HTML 4.01 Transitional//EN">
   2 <html><head><meta http-equiv="Content-Type" content="text/html;charset=iso-8859-1">
   3 <title>AliHLTPHOS: AliHLTPHOSPulseGenerator.cxx Source File</title>
   4 <link href="doxygen.css" rel="stylesheet" type="text/css">
   5 </head><body>
   6 <!-- Generated by Doxygen 1.3.9.1 -->
   7 <div class="qindex"><a class="qindex" href="main.html">Main&nbsp;Page</a> | <a class="qindex" href="hierarchy.html">Class&nbsp;Hierarchy</a> | <a class="qindex" href="classes.html">Alphabetical&nbsp;List</a> | <a class="qindex" href="annotated.html">Class&nbsp;List</a> | <a class="qindex" href="files.html">File&nbsp;List</a> | <a class="qindex" href="functions.html">Class&nbsp;Members</a> | <a class="qindex" href="globals.html">File&nbsp;Members</a></div>
   8 <h1>AliHLTPHOSPulseGenerator.cxx</h1><a href="AliHLTPHOSPulseGenerator_8cxx.html">Go to the documentation of this file.</a><div class="fragment"><pre class="fragment">00001 <span class="comment">/**************************************************************************</span>
   9 00002 <span class="comment"> * Copyright(c) 2006, ALICE Experiment at CERN, All rights reserved.      *</span>
  10 00003 <span class="comment"> *                                                                        *</span>
  11 00004 <span class="comment"> * Author: Per Thomas Hille for the ALICE HLT Project.                    *</span>
  12 00005 <span class="comment"> * Contributors are mentioned in the code where appropriate.              *</span>
  13 00006 <span class="comment"> *                                                                        *</span>
  14 00007 <span class="comment"> * Permission to use, copy, modify and distribute this software and its   *</span>
  15 00008 <span class="comment"> * documentation strictly for non-commercial purposes is hereby granted   *</span>
  16 00009 <span class="comment"> * without fee, provided that the above copyright notice appears in all   *</span>
  17 00010 <span class="comment"> * copies and that both the copyright notice and this permission notice   *</span>
  18 00011 <span class="comment"> * appear in the supporting documentation. The authors make no claims     *</span>
  19 00012 <span class="comment"> * about the suitability of this software for any purpose. It is          *</span>
  20 00013 <span class="comment"> * provided "as is" without express or implied warranty.                  *</span>
  21 00014 <span class="comment"> **************************************************************************/</span>
  22 00015
  23 00016
  24 00017 <span class="preprocessor">#include "<a class="code" href="AliHLTPHOSPulseGenerator_8h.html">AliHLTPHOSPulseGenerator.h</a>"</span>
  25 00018 <span class="comment">//#include &lt;stdio.h&gt;</span>
  26 00019 <span class="preprocessor">#include &lt;cmath&gt;</span>
  27 00020 <span class="preprocessor">#include &lt;iostream&gt;</span>
  28 00021
  29 00022
  30 00023
  31 00024 <span class="keyword">using</span> std::cout;
  32 00025 <span class="keyword">using</span> std::endl;
  33 00026
  34 00027
  35 <a name="l00028"></a><a class="code" href="AliHLTPHOSPulseGenerator_8cxx.html#a0">00028</a> <a class="code" href="AliHLTPHOSRawAnalyzerPeakFinder_8cxx.html#a0">ClassImp</a>(<a class="code" href="classAliHLTPHOSPulseGenerator.html">AliHLTPHOSPulseGenerator</a>)
  36 00029
  37 00030
  38 00033 <a class="code" href="classAliHLTPHOSPulseGenerator.html">AliHLTPHOSPulseGenerator</a>::<a class="code" href="classAliHLTPHOSPulseGenerator.html">AliHLTPHOSPulseGenerator</a>(): fAmplitude(0), fNSamples(0),fTau(0), fSampleFreq(0), fTZero(0), fDataPtr(0), fDT(0), fEvent(0)
  39 00034 {
  40 00035   cout &lt;&lt; <span class="stringliteral">"You cannot invoke the Pulsgenerator without parameters"</span> &lt;&lt; endl;
  41 00036 }
  42 00037
  43 <a name="l00038"></a><a class="code" href="classAliHLTPHOSPulseGenerator.html#AliHLTPHOSPulseGeneratora3">00038</a> <a class="code" href="classAliHLTPHOSPulseGenerator.html#AliHLTPHOSPulseGeneratora0">AliHLTPHOSPulseGenerator::AliHLTPHOSPulseGenerator</a>(<span class="keyword">const</span> <a class="code" href="classAliHLTPHOSPulseGenerator.html">AliHLTPHOSPulseGenerator</a> &amp;): fAmplitude(0), fNSamples(0),fTau(0), fSampleFreq(0), fTZero(0), fDataPtr(0), fDT(0), fEvent(0)
  44 00039 {
  45 00040
  46 00041 }
  47 00042
  48 00043
  49 <a name="l00044"></a><a class="code" href="classAliHLTPHOSPulseGenerator.html#AliHLTPHOSPulseGeneratora1">00044</a> <a class="code" href="classAliHLTPHOSPulseGenerator.html#AliHLTPHOSPulseGeneratora1">AliHLTPHOSPulseGenerator::~AliHLTPHOSPulseGenerator</a>()
  50 00045 {
  51 00046   <span class="keyword">delete</span> fDataPtr;
  52 00047   <a class="code" href="classAliHLTPHOSPulseGenerator.html#AliHLTPHOSPulseGeneratorr5">fDataPtr</a>=0;
  53 00048 }
  54 00049
  55 <a name="l00056"></a><a class="code" href="classAliHLTPHOSPulseGenerator.html#AliHLTPHOSPulseGeneratora2">00056</a> <a class="code" href="classAliHLTPHOSPulseGenerator.html#AliHLTPHOSPulseGeneratora0">AliHLTPHOSPulseGenerator::AliHLTPHOSPulseGenerator</a>(<span class="keywordtype">double</span> a, <span class="keywordtype">double</span> t0, <span class="keywordtype">int</span> N, <span class="keywordtype">double</span> t, <span class="keywordtype">double</span> fs): fAmplitude(a), fNSamples(N),fTau(0), fSampleFreq(fs), fTZero(0), fDataPtr(0), fDT(0), fEvent(0)
  56 00057 {
  57 00058
  58 00059   <a class="code" href="classAliHLTPHOSPulseGenerator.html#AliHLTPHOSPulseGeneratorr5">fDataPtr</a> = <span class="keyword">new</span> <span class="keywordtype">double</span>[100];
  59 00060
  60 00061
  61 00062
  62 00063   <a class="code" href="classAliHLTPHOSPulseGenerator.html#AliHLTPHOSPulseGeneratora12">SetAmplitude</a>(a);
  63 00064   <a class="code" href="classAliHLTPHOSPulseGenerator.html#AliHLTPHOSPulseGeneratora13">SetDT</a>(fs);
  64 00065   <a class="code" href="classAliHLTPHOSPulseGenerator.html#AliHLTPHOSPulseGeneratora14">SetTZero</a>(t0);
  65 00066   <a class="code" href="classAliHLTPHOSPulseGenerator.html#AliHLTPHOSPulseGeneratorr1">fNSamples</a>=N;
  66 00067   <a class="code" href="classAliHLTPHOSPulseGenerator.html#AliHLTPHOSPulseGeneratorr2">fTau</a>=t;
  67 00068   <a class="code" href="classAliHLTPHOSPulseGenerator.html#AliHLTPHOSPulseGeneratorr3">fSampleFreq</a>=fs;
  68 00069   <span class="comment">//  dT=tau/fs;   //Function values are calculated at intervals dT</span>
  69 00070   <span class="comment">//  fDT=1/fs;   //Function values are calculated at intervals dT</span>
  70 00071   <a class="code" href="classAliHLTPHOSPulseGenerator.html#AliHLTPHOSPulseGeneratord0">MakePulse</a>(<a class="code" href="classAliHLTPHOSPulseGenerator.html#AliHLTPHOSPulseGeneratorr5">fDataPtr</a>);
  71 00072 }
  72 00073
  73 00079 <span class="keywordtype">void</span>
  74 <a name="l00080"></a><a class="code" href="classAliHLTPHOSPulseGenerator.html#AliHLTPHOSPulseGeneratora5">00080</a> <a class="code" href="classAliHLTPHOSPulseGenerator.html#AliHLTPHOSPulseGeneratora5">AliHLTPHOSPulseGenerator::AddBaseline</a>(<span class="keywordtype">double</span> baselineLevel, <span class="keywordtype">double</span> *samples)
  75 00081 {
  76 00082   <span class="keywordtype">double</span> *tmpSamples;
  77 00083   tmpSamples = samples;
  78 00084   printf(<span class="stringliteral">"\nbaselineLevel = %f\n"</span>, baselineLevel);
  79 00085   cout &lt;&lt; <span class="stringliteral">"AddBaseline not implemented yet"</span> &lt;&lt; endl;
  80 00086 }
  81 00087
  82 00093 <span class="keywordtype">void</span>
  83 <a name="l00094"></a><a class="code" href="classAliHLTPHOSPulseGenerator.html#AliHLTPHOSPulseGeneratora6">00094</a> <a class="code" href="classAliHLTPHOSPulseGenerator.html#AliHLTPHOSPulseGeneratora6">AliHLTPHOSPulseGenerator::AddNoise</a>(<span class="keywordtype">double</span> *dataPtr, <span class="keywordtype">double</span> *sigma)
  84 00095 {
  85 00096   printf(<span class="stringliteral">"\ndataPtr = %f, sigma = %f\n"</span>, *dataPtr, *sigma);
  86 00097   cout &lt;&lt; <span class="stringliteral">"AddNoise is not implemented yet"</span> &lt;&lt; endl;
  87 00098 }
  88 00099
  89 00100
  90 00107 <span class="keywordtype">void</span>
  91 <a name="l00108"></a><a class="code" href="classAliHLTPHOSPulseGenerator.html#AliHLTPHOSPulseGeneratora7">00108</a> <a class="code" href="classAliHLTPHOSPulseGenerator.html#AliHLTPHOSPulseGeneratora6">AliHLTPHOSPulseGenerator::AddNoise</a>(<span class="keywordtype">double</span> *dataPtr, <span class="keywordtype">double</span> *sigma, <span class="keywordtype">double</span> cutoff)
  92 00109 {
  93 00110   printf(<span class="stringliteral">"\ndataPtr = %f, sigma = %f, cutoff = %f\n"</span>, *dataPtr, *sigma, cutoff);
  94 00111   cout &lt;&lt; <span class="stringliteral">"AddNoise is not implemeted yet"</span> &lt;&lt; endl;
  95 00112 }
  96 00113
  97 00121 <span class="keywordtype">double</span> *
  98 <a name="l00122"></a><a class="code" href="classAliHLTPHOSPulseGenerator.html#AliHLTPHOSPulseGeneratora8">00122</a> <a class="code" href="classAliHLTPHOSPulseGenerator.html#AliHLTPHOSPulseGeneratora8">AliHLTPHOSPulseGenerator::AddPretriggerSamples</a>(<span class="keywordtype">double</span> baselineLevel, <span class="keywordtype">double</span> *samples)
  99 00123 {
 100 00124   printf(<span class="stringliteral">"\nbaslinelevel = %f, samples = %f\n"</span>, baselineLevel, *samples);
 101 00125   cout &lt;&lt; <span class="stringliteral">"AddPretriggerSamples not implemented yet"</span> &lt;&lt; endl;
 102 00126   <span class="keywordflow">return</span> 0;
 103 00127 }
 104 00128
 105 00129
 106 00133 <span class="keywordtype">double</span> *
 107 <a name="l00134"></a><a class="code" href="classAliHLTPHOSPulseGenerator.html#AliHLTPHOSPulseGeneratora9">00134</a> <a class="code" href="classAliHLTPHOSPulseGenerator.html#AliHLTPHOSPulseGeneratora9">AliHLTPHOSPulseGenerator::GetPulse</a>()
 108 00135 {
 109 00136   <span class="keywordflow">return</span> fDataPtr;
 110 00137 }
 111 00138
 112 00139
 113 00144 <span class="keywordtype">double</span> *
 114 <a name="l00145"></a><a class="code" href="classAliHLTPHOSPulseGenerator.html#AliHLTPHOSPulseGeneratora10">00145</a> <a class="code" href="classAliHLTPHOSPulseGenerator.html#AliHLTPHOSPulseGeneratora9">AliHLTPHOSPulseGenerator::GetPulse</a>(<span class="keywordtype">double</span> a, <span class="keywordtype">double</span> t0)
 115 00146 {
 116 00147   <span class="keywordflow">return</span> fDataPtr;
 117 00148 }
 118 00149
 119 00154 <span class="keywordtype">void</span>
 120 <a name="l00155"></a><a class="code" href="classAliHLTPHOSPulseGenerator.html#AliHLTPHOSPulseGeneratora11">00155</a> <a class="code" href="classAliHLTPHOSPulseGenerator.html#AliHLTPHOSPulseGeneratora11">AliHLTPHOSPulseGenerator::Quantisize</a>(<span class="keywordtype">double</span> *dataPtr)
 121 00156 {
 122 00157   <span class="keywordtype">double</span> *dtaPtr;
 123 00158   dtaPtr = <span class="keyword">new</span> <span class="keywordtype">double</span>[100];
 124 00159   dtaPtr = dataPtr;
 125 00160   <span class="comment">//  cout &lt;&lt; "Quantisize is not implemented yet" &lt;&lt; endl;</span>
 126 00161 }
 127 00162
 128 00163 <span class="keywordtype">void</span>
 129 <a name="l00164"></a><a class="code" href="classAliHLTPHOSPulseGenerator.html#AliHLTPHOSPulseGeneratora12">00164</a> <a class="code" href="classAliHLTPHOSPulseGenerator.html#AliHLTPHOSPulseGeneratora12">AliHLTPHOSPulseGenerator::SetAmplitude</a>(<span class="keywordtype">double</span> a)
 130 00165 {
 131 00166   <a class="code" href="classAliHLTPHOSPulseGenerator.html#AliHLTPHOSPulseGeneratorr0">fAmplitude</a>=a;
 132 00167 }
 133 00168
 134 00169 <span class="keywordtype">void</span>
 135 <a name="l00170"></a><a class="code" href="classAliHLTPHOSPulseGenerator.html#AliHLTPHOSPulseGeneratora13">00170</a> <a class="code" href="classAliHLTPHOSPulseGenerator.html#AliHLTPHOSPulseGeneratora13">AliHLTPHOSPulseGenerator::SetDT</a>(<span class="keywordtype">double</span> fs)
 136 00171 {
 137 00172   <a class="code" href="classAliHLTPHOSPulseGenerator.html#AliHLTPHOSPulseGeneratorr6">fDT</a>=1/fs;
 138 00173 }
 139 00174
 140 00175 <span class="keywordtype">void</span>
 141 <a name="l00176"></a><a class="code" href="classAliHLTPHOSPulseGenerator.html#AliHLTPHOSPulseGeneratora14">00176</a> <a class="code" href="classAliHLTPHOSPulseGenerator.html#AliHLTPHOSPulseGeneratora14">AliHLTPHOSPulseGenerator::SetTZero</a>(<span class="keywordtype">double</span> t0)
 142 00177 {
 143 00178   <a class="code" href="classAliHLTPHOSPulseGenerator.html#AliHLTPHOSPulseGeneratorr4">fTZero</a> = -t0/1000; <span class="comment">// Since time is in nanoseconds and the samplingfrequency is in MHz -&gt; divide by 1000</span>
 144 00179 }
 145 00180
 146 00181
 147 00182 <span class="keywordtype">void</span>
 148 <a name="l00183"></a><a class="code" href="classAliHLTPHOSPulseGenerator.html#AliHLTPHOSPulseGeneratord0">00183</a> <a class="code" href="classAliHLTPHOSPulseGenerator.html#AliHLTPHOSPulseGeneratord0">AliHLTPHOSPulseGenerator::MakePulse</a>(<span class="keywordtype">double</span> *dtaPtr)
 149 00184 {
 150 00185 <span class="keywordflow">for</span>(<span class="keywordtype">int</span> i=0; i&lt;fNSamples; i++)
 151 00186   {
 152 00187     dtaPtr[i]=<a class="code" href="classAliHLTPHOSPulseGenerator.html#AliHLTPHOSPulseGeneratorr0">fAmplitude</a>*exp((Double_t)2)*pow((i*<a class="code" href="classAliHLTPHOSPulseGenerator.html#AliHLTPHOSPulseGeneratorr6">fDT</a>-<a class="code" href="classAliHLTPHOSPulseGenerator.html#AliHLTPHOSPulseGeneratorr4">fTZero</a>)/<a class="code" href="classAliHLTPHOSPulseGenerator.html#AliHLTPHOSPulseGeneratorr2">fTau</a>, 2)*exp(-2*(i*<a class="code" href="classAliHLTPHOSPulseGenerator.html#AliHLTPHOSPulseGeneratorr6">fDT</a>-fTZero)/fTau);
 153 00188   }
 154 00189 }
 155 </pre></div><hr size="1"><address style="align: right;"><small>Generated on Sun Feb 18 22:17:03 2007 for AliHLTPHOS by&nbsp;
 156 <a href="http://www.doxygen.org/index.html">
 157 <img src="doxygen.png" alt="doxygen" align="middle" border="0"></a> 1.3.9.1 </small></address>
 158 </body>
 159 </html>