X-Git-Url: http://git.uio.no/git/?a=blobdiff_plain;f=STEER%2FAliGausCorr.cxx;h=24df66a38db54e4aa98640e8856b7b3f5fb71123;hb=4cfb400721deb33b5bad5aed8516f5045eadfa6b;hp=15f8f2663e7f75d4ba1fb350e16a507bb556b891;hpb=740ebff31b6584a78bc7e49341261dc254ba3e71;p=u%2Fmrichter%2FAliRoot.git

diff --git a/STEER/AliGausCorr.cxx b/STEER/AliGausCorr.cxx
index 15f8f2663e7..24df66a38db 100644
--- a/STEER/AliGausCorr.cxx
+++ b/STEER/AliGausCorr.cxx
@@ -13,6 +13,8 @@
  * provided "as is" without express or implied warranty.                  *
  **************************************************************************/
 
+/* $Id$ */
+
 ////////////////////////////////////////////////////////////////////////
 // Class used to generate correlated gaussian numbers with mean
 // zero and known covariance matrix.
@@ -21,28 +23,34 @@
 // Rep. Prog. Phys. 43 (1980) 1145-1189. 
 // M.Masera 15.03.2001 9:30 - modified on 26.02.2002 17:40
 ////////////////////////////////////////////////////////////////////////
-#include <iostream.h>
+
+#include <Riostream.h>
 #include <TArrayD.h>
-#include <TMatrixD.h>
+#include <TMath.h>
 #include <TRandom.h>
+
 #include "AliGausCorr.h"
 
 ClassImp(AliGausCorr)
 
-//________________________________________________________
-AliGausCorr::AliGausCorr()
+//_______________________________________________________________________
+AliGausCorr::AliGausCorr():
+  fSize(0),
+  fCv(0)
 {
+  //
   // Default constructor
-  fSize = 0;
-  fCv = 0;
+  //
 }
 
-//________________________________________________________
-AliGausCorr::AliGausCorr(const TMatrixD & vec, Int_t size) 
+//_______________________________________________________________________
+AliGausCorr::AliGausCorr(const TMatrixD & vec, Int_t size):
+  fSize(size),
+  fCv(new TMatrixD(fSize,fSize))
 {
+  //
   // Standard constructor
-  fSize = size;
-  fCv = new TMatrixD(fSize,fSize);
+  //
   for(Int_t j=0;j<fSize;j++){
     double accum = 0;
     for(Int_t k=0;k<j;k++){
@@ -59,26 +67,28 @@ AliGausCorr::AliGausCorr(const TMatrixD & vec, Int_t size)
   }
 }
 
-//________________________________________________________
-AliGausCorr::AliGausCorr(const AliGausCorr & tgcorr)
+//_______________________________________________________________________
+AliGausCorr::AliGausCorr(const AliGausCorr & tgcorr):
+  TObject(tgcorr),
+  fSize(tgcorr.fSize),
+  fCv(new TMatrixD(fSize,fSize))
 {
+  //
   // Copy contructor
-
-  fSize = tgcorr.fSize;
-  fCv = new TMatrixD(fSize,fSize);
+  //
   for(Int_t i=0;i<fSize;i++){
     for(Int_t j=0;j<fSize;j++)(*fCv)(i,j)=(*tgcorr.fCv)(i,j);
   }
 }
 
-//________________________________________________________
+//_______________________________________________________________________
 AliGausCorr::~AliGausCorr()
 {
   // Destructor
   delete fCv;
 }
 
-//________________________________________________________
+//_______________________________________________________________________
 void AliGausCorr::GetGaussN(TArrayD &vec) const
 {
   // return fSize correlated gaussian numbers
@@ -102,7 +112,7 @@ void AliGausCorr::GetGaussN(TArrayD &vec) const
 
 }
 
-//________________________________________________________
+//_______________________________________________________________________
 void AliGausCorr::PrintCv() const
 {
   // Printout of the "square root" cov. matrix 
@@ -121,9 +131,10 @@ void AliGausCorr::PrintCv() const
   printf("\n");
 }
 
-//________________________________________________________
+//_______________________________________________________________________
 AliGausCorr & AliGausCorr::operator=(const AliGausCorr & tgcorr)
 {
+  // Assignment operator
   if(&tgcorr != this && tgcorr.fSize!=fSize){
     if(fCv)delete fCv;
     fSize = tgcorr.fSize;
@@ -137,8 +148,3 @@ AliGausCorr & AliGausCorr::operator=(const AliGausCorr & tgcorr)
 
   return *this;
 }
-
-
-
-
-