Spelling

[u/mrichter/AliRoot.git] / PHOS / AliPHOSPIDv1.cxx

diff --git a/PHOS/AliPHOSPIDv1.cxx b/PHOS/AliPHOSPIDv1.cxx
index 86b1598a3177754da67ba1c7311ad0342789800d..fda233fdc6f77570a66bae4f77f0114715a0d18c 100644 (file)
--- a/PHOS/AliPHOSPIDv1.cxx
+++ b/PHOS/AliPHOSPIDv1.cxx
@@ -160,7 +160,7 @@ AliPHOSPIDv1::AliPHOSPIDv1() :
   fPPhoton(0),
   fPPi0(0),
   fParameters(0),
-  fVtx(0.), 
+  fVtx(0.,0.,0.), 
   fTFphoton(0),
   fTFpiong(0),
   fTFkaong(0),
@@ -195,7 +195,7 @@ AliPHOSPIDv1::AliPHOSPIDv1(const AliPHOSPIDv1 & pid ) :
   fPPhoton(0),
   fPPi0(0),
   fParameters(0),
-  fVtx(0.), 
+  fVtx(0.,0.,0.), 
   fTFphoton(0),
   fTFpiong(0),
   fTFkaong(0),
@@ -230,7 +230,7 @@ AliPHOSPIDv1::AliPHOSPIDv1(AliPHOSGeometry *geom):
   fPPhoton(0),
   fPPi0(0),
   fParameters(0),
-  fVtx(0.), 
+  fVtx(0.,0.,0.), 
   fTFphoton(0),
   fTFpiong(0),
   fTFkaong(0),
@@ -473,7 +473,7 @@ void AliPHOSPIDv1::InitParameters()
   fERecWeight ->SetParameters(fERecWeightPar[0],fERecWeightPar[1] ,fERecWeightPar[2] ,fERecWeightPar[3]) ; 
 
 
-  for (Int_t i =0; i<  AliPID::kSPECIESN ; i++)
+  for (Int_t i =0; i<  AliPID::kSPECIESCN ; i++)
     fInitPID[i] = 1.;
  
 }
@@ -944,7 +944,7 @@ void  AliPHOSPIDv1::MakePID()
 {
   // construct the PID weight from a Bayesian Method
   
-  const Int_t kSPECIES = AliPID::kSPECIESN ;
+  const Int_t kSPECIES = AliPID::kSPECIESCN ;
  
   Int_t nparticles = fRecParticles->GetEntriesFast() ;
 
@@ -962,8 +962,8 @@ void  AliPHOSPIDv1::MakePID()
   Double_t * sw  [kSPECIES] ;
   //Info("MakePID","Begin MakePID"); 
   
-  for (Int_t i =0; i< kSPECIES; i++){
-    stof[i] = new Double_t[nparticles] ;
+         for (Int_t i =0; i< kSPECIES; i++){
+           stof[i] = new Double_t[nparticles] ;
     sdp [i] = new Double_t[nparticles] ;
     scpv[i] = new Double_t[nparticles] ;
     sw  [i] = new Double_t[nparticles] ;
@@ -995,10 +995,13 @@ void  AliPHOSPIDv1::MakePID()
     //    Info("MakePID", "TOF");
     Float_t  en   = emc->GetEnergy();    
     Double_t time = emc->GetTime() ;
-    //    cout<<">>>>>>>Energy "<<en<<"Time "<<time<<endl;
+    //      cout<<">>>>>>>Energy "<<en<<"Time "<<time<<endl;
    
     // now get the signals probability
     // s(pid) in the Bayesian formulation
+
+    //Initialize anused species
+    for(Int_t iii=0; iii<kSPECIES; iii++)stof[iii][index]=0. ;
     
     stof[AliPID::kPhoton][index]   = 1.; 
     stof[AliPID::kElectron][index] = 1.;
@@ -1057,6 +1060,8 @@ void  AliPHOSPIDv1::MakePID()
     //DP: still to be done 
 
     //dispersion is not well defined if the cluster is only in few crystals
+    //Initialize anused species
+    for(Int_t iii=0; iii<kSPECIES; iii++)sdp[iii][index]=0. ;
     
     sdp[AliPID::kPhoton][index]   = 1. ;
     sdp[AliPID::kElectron][index] = 1. ;
@@ -1121,6 +1126,8 @@ void  AliPHOSPIDv1::MakePID()
       }
     //    else
     //      cout<<">>>>>>>>>>>CHARGED>>>>>>>>>>>"<<endl;
+    //Initialize anused species
+    for(Int_t iii=0; iii<kSPECIES; iii++)scpv[iii][index]=0. ;
     
     scpv[AliPID::kPion][index]     =  pcpvcharged  ; 
     scpv[AliPID::kKaon][index]     =  pcpvcharged  ; 
@@ -1168,6 +1175,8 @@ void  AliPHOSPIDv1::MakePID()
     }
 
     //Weight to apply to hadrons due to energy reconstruction
+    //Initialize anused species
+    for(Int_t iii=0; iii<kSPECIES; iii++)sw[iii][index]=1. ;
 
     Float_t weight = fERecWeight ->Eval(en) ;
  
@@ -1632,10 +1641,10 @@ void AliPHOSPIDv1::GetVertex(void)
 //_______________________________________________________________________
 void AliPHOSPIDv1::SetInitPID(const Double_t *p) {
   // Sets values for the initial population of each particle type 
-  for (Int_t i=0; i<AliPID::kSPECIESN; i++) fInitPID[i] = p[i];
+  for (Int_t i=0; i<AliPID::kSPECIESCN; i++) fInitPID[i] = p[i];
 }
 //_______________________________________________________________________
 void AliPHOSPIDv1::GetInitPID(Double_t *p) const {
   // Gets values for the initial population of each particle type 
-  for (Int_t i=0; i<AliPID::kSPECIESN; i++) p[i] = fInitPID[i];
+  for (Int_t i=0; i<AliPID::kSPECIESCN; i++) p[i] = fInitPID[i];
 }