]> git.uio.no Git - u/mrichter/AliRoot.git/blobdiff - STEER/AliAODRedCov.cxx
Implementation of BitPattern for ACORDE in AliACORDEQADataMakerSim
[u/mrichter/AliRoot.git] / STEER / AliAODRedCov.cxx
index 75494ffaab129859901f6771d44fb3cac357fa31..95b4142b27fdd80093e5962da039820e9b18fdec 100644 (file)
 
 /* $Id$ */
 
-#include <TMath.h>
-
 #include "AliAODRedCov.h"
 
 templateClassImp(AliAODRedCov)
 
-//______________________________________________________________________________
-template <Int_t N> template <class T> void AliAODRedCov<N>::GetCovMatrix(T *cmat) const
-{
-  //
-  // Returns the external cov matrix
-  //
-
-  for(Int_t i=0; i<N; ++i) {
-    // Off diagonal elements
-    for(Int_t j=0; j<i; ++j) {
-      printf("cmat[%2d] = fODia[%2d]*fDiag[%2d]*fDiag[%2d];\n",
-            i*(i+1)/2+j,(i-1)*i/2+j,j,i);
-      cmat[i*(i+1)/2+j] = fODia[(i-1)*i/2+j]*fDiag[j]*fDiag[i];}
-
-    // Diagonal elements
-    printf("cmat[%2d] = fDiag[%2d]*fDiag[%2d];\n",
-          i*(i+1)/2+i,i,i);
-    cmat[i*(i+1)/2+i] = fDiag[i]*fDiag[i];
-  }
-}
-
-
-//______________________________________________________________________________
-template <Int_t N> template <class T> void AliAODRedCov<N>::SetCovMatrix(T *cmat)
-{
-  //
-  // Sets the external cov matrix
-  //
-
-  if(cmat) {
-    // Diagonal elements first
-    for(Int_t i=0; i<N; ++i) {
-      printf("fDiag[%2d] = TMath::Sqrt(cmat[%2d]);\n",
-            i,i*(i+1)/2+i);
-      fDiag[i] = TMath::Sqrt(cmat[i*(i+1)/2+i]);}
-
-  // ... then the ones off diagonal
-  for(Int_t i=0; i<N; ++i) 
-    // Off diagonal elements
-    for(Int_t j=0; j<i; ++j) {
-      printf("fODia[%2d] = cmat[%2d]/(fDiag[%2d]*fDiag[%2d]);\n",
-            (i-1)*i/2+j,i*(i+1)/2+j,j,i);
-      fODia[(i-1)*i/2+j] = cmat[i*(i+1)/2+j]/(fDiag[j]*fDiag[i]);
-    }
-  } else {
-    for(Int_t i=0; i< N; ++i) fDiag[i]=-999.;
-    for(Int_t i=0; i< N*(N-1)/2; ++i) fODia[i]=0.;
-  }
-}