Typo corrected

[u/mrichter/AliRoot.git] / ITS / AliITSetfSDD.cxx

diff --git a/ITS/AliITSetfSDD.cxx b/ITS/AliITSetfSDD.cxx
index 2b87db07930f4c9a33ff41fb89b793edf3d331e3..817bbb728d9c365a9a0cf1487ea7e51fd076cd0f 100644 (file)
--- a/ITS/AliITSetfSDD.cxx
+++ b/ITS/AliITSetfSDD.cxx
@@ -12,9 +12,23 @@
  * about the suitability of this software for any purpose. It is          *
  * provided "as is" without express or implied warranty.                  *
  **************************************************************************/
+ 
+/*
+$Log$
+Revision 1.8  2002/10/14 14:57:00  hristov
+Merging the VirtualMC branch to the main development branch (HEAD)
 
-#include <iostream.h>
+Revision 1.6.6.1  2002/06/10 17:51:15  hristov
+Merged with v3-08-02
+
+Revision 1.7  2002/04/24 22:02:31  nilsen
+New SDigits and Digits routines, and related changes,  (including new
+noise values).
+
+*/
+#include <Riostream.h>
 #include <TMath.h>
+#include <TString.h>
 #include "AliITSetfSDD.h"
 
 ////////////////////////////////////////////////////////////////////////
@@ -33,15 +47,16 @@ Int_t ppower(Int_t b, Int_t e) {
   return power;
 }
 
-AliITSetfSDD::AliITSetfSDD(Double_t timestep)
+AliITSetfSDD::AliITSetfSDD(Double_t timestep, Int_t amplif)
 {
   // sampling time in ns
 
+  fTimeDelay = 53.5;
+  if(amplif == 2) fTimeDelay = 35.5;
   fSamplingTime = timestep;
 
   fT0 = 0.;
   fDf = ppower(10,9)/(kMaxNofSamples*fSamplingTime);
-  fA0 = 9000.;
 
   Int_t i,j;
   for(i=0; i<kMaxNofPoles; i++) {
@@ -52,14 +67,38 @@ AliITSetfSDD::AliITSetfSDD(Double_t timestep)
     fPoleR[i] = 0.;
     fPoleI[i] = 0.;
   }
-  fPoleM[0] = 1.;
-  fPoleR[0] = -2100000.;
-  fPoleI[0] = fPoleR[0];
-  fPoleM[1] = 1.;
-  fPoleR[1] = -2100000.;
-  fPoleI[1] = -fPoleR[1];
+  // Alice
+
+  // PASCAL amplif
+  fA0 = 5.53269815e+11; 
+  fPoleM[0] = 3.;
+  fPoleR[0] = -8280000.; 
+  fPoleI[0] = 0.; 
+
+  if(amplif == 2) { // OLA amplif.
+    fA0 = 24000.;
+    fPoleM[0] = 1.;
+    fPoleR[0] = -3000000.;
+    fPoleI[0] = 4000000.;
+    fPoleM[1] = 1.;
+    fPoleR[1] = fPoleR[0];
+    fPoleI[1] = -fPoleI[0]; 
+  }
 
-   // Compute Transfer Function
+  if( amplif == 3 ) { // old PASCAL
+    fA0 = 16500.; // AL: 16500.;  // TB: 24000.; // 26000.; // 24000.; // 18000.; 
+    fPoleM[0] = 1.;
+    fPoleR[0] = -4140000.; // AL: -4140000.; // TB: -3000000.; // -3750000.; // -3500000; // -3000000.; 
+    fPoleI[0] = 0.; // AL: 0.; // TB: 4000000.; // 3750000.; // 3500000.; // 3000000.; 
+    fPoleM[1] = 1.;
+    fPoleR[1] = fPoleR[0];
+    fPoleI[1] = -fPoleI[0]; 
+  }
+
+  //cout << "fA0: " << fA0 << endl;
+  //cout << "fTimeDelay: " << fTimeDelay << endl;
+  
+  // Compute Transfer Function
 
   Double_t PI = acos(-1.);
   for(i=0; i<=kMaxNofSamples/2; i++) {
@@ -128,6 +167,7 @@ AliITSetfSDD::AliITSetfSDD(Double_t timestep)
 
 void AliITSetfSDD::PrintElectronics()
 {
+  cout << "Time Delay " << fTimeDelay << endl;
   cout << "Sampling Time " << fSamplingTime << endl;
   cout << "Number of Time Samples " << kMaxNofSamples << endl;
   cout << "fT0 " << fT0 << endl;