fADCValuesHigh4x4 is not persistent

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

diff --git a/PHOS/AliPHOSPIDv0.cxx b/PHOS/AliPHOSPIDv0.cxx
index 9479caa6aeb70b40583d6036c56301397e05a65d..2516937dbeb8f21e3901fe19831de2f1ec17ceb1 100644 (file)
--- a/PHOS/AliPHOSPIDv0.cxx
+++ b/PHOS/AliPHOSPIDv0.cxx
@@ -15,6 +15,11 @@
 
 /* $Id$ */
 
+/* History of cvs commits:
+ *
+ * $Log$
+ */
+
 //_________________________________________________________________________
 // Implementation version v0 of the PHOS particle identifier 
 // Particle identification based on the 
@@ -57,26 +62,21 @@
 //            Completely redesined by Dmitri Peressounko, March 2001
 
 // --- ROOT system ---
-#include "TROOT.h"
 #include "TTree.h"
-#include "TFile.h"
 #include "TF2.h"
 #include "TFormula.h"
 #include "TCanvas.h"
-#include "TFolder.h"
-#include "TSystem.h"
+
 #include "TBenchmark.h"
 // --- Standard library ---
 
 // --- AliRoot header files ---
-
+#include "AliLog.h"
 #include "AliRun.h"
 #include "AliGenerator.h"
-#include "AliPHOS.h"
 #include "AliPHOSPIDv0.h"
-#include "AliPHOSClusterizerv1.h"
+#include "AliPHOSEmcRecPoint.h"
 #include "AliPHOSTrackSegment.h"
-#include "AliPHOSTrackSegmentMakerv1.h"
 #include "AliPHOSRecParticle.h"
 #include "AliPHOSGeometry.h"
 #include "AliPHOSLoader.h"
@@ -175,14 +175,15 @@ void  AliPHOSPIDv0::Exec(Option_t * option)
   AliRunLoader* runget = AliRunLoader::GetRunLoader(GetTitle());
   if(runget == 0x0) 
    {
-     Error("Exec","Can not find run getter in event folder \"%s\"",GetTitle());
+     AliError(Form("Can not find run getter in event folder \"%s\"",
+                  GetTitle()));
      return;
    }
   
   AliPHOSLoader* gime = dynamic_cast<AliPHOSLoader*>(runget->GetLoader("PHOSLoader"));
   if ( gime == 0 ) 
    {
-     Error("Exec","Could not obtain the Loader object !"); 
+     AliError("Could not obtain the Loader object !"); 
      return ;
    } 
 
@@ -196,7 +197,9 @@ void  AliPHOSPIDv0::Exec(Option_t * option)
   
   for(ievent = 0; ievent < nevents; ievent++){
     runget->GetEvent(ievent);
-    Info("Exec", "event %d %d %d", ievent, gime->EmcRecPoints(), gime->TrackSegments()) ;
+    AliInfo(Form("event %d %d %d", 
+                ievent, gime->EmcRecPoints(), 
+                gime->TrackSegments())) ;
     MakeRecParticles() ;
     
     WriteRecParticles();
@@ -211,8 +214,9 @@ void  AliPHOSPIDv0::Exec(Option_t * option)
   
   if(strstr(option,"tim")){
     gBenchmark->Stop("PHOSPID");
-    Info("Exec", "took %f seconds for PID %f seconds per event", 
-        gBenchmark->GetCpuTime("PHOSPID"), gBenchmark->GetCpuTime("PHOSPID")/nevents) ; 
+    AliInfo(Form("took %f seconds for PID %f seconds per event", 
+                gBenchmark->GetCpuTime("PHOSPID"), 
+                gBenchmark->GetCpuTime("PHOSPID")/nevents)) ; 
   }
   
 }
@@ -225,14 +229,15 @@ void AliPHOSPIDv0::Init()
   AliRunLoader* runget = AliRunLoader::GetRunLoader(GetTitle());
   if(runget == 0x0) 
    {
-     Error("Exec","Can not find run getter in event folder \"%s\"",GetTitle());
+     AliError(Form("Can not find run getter in event folder \"%s\"",
+                  GetTitle()));
      return;
    }
   
   AliPHOSLoader* gime = dynamic_cast<AliPHOSLoader*>(runget->GetLoader("PHOSLoader"));
   if ( gime == 0 ) 
    {
-     Error("Exec","Could not obtain the Loader object !"); 
+     AliError("Could not obtain the Loader object !"); 
      return ;
    } 
    
@@ -242,7 +247,9 @@ void AliPHOSPIDv0::Init()
 }
 
 //____________________________________________________________________________
-void  AliPHOSPIDv0::MakeRecParticles(){
+void  AliPHOSPIDv0::MakeRecParticles()
+{
+  // Reconstructs the particles from the tracksegments
 
   TString taskName(GetName()) ; 
   taskName.Remove(taskName.Index(Version())-1) ;
@@ -250,14 +257,15 @@ void  AliPHOSPIDv0::MakeRecParticles(){
   AliRunLoader* runget = AliRunLoader::GetRunLoader(GetTitle());
   if(runget == 0x0) 
    {
-     Error("Exec","Can not find run getter in event folder \"%s\"",GetTitle());
+     AliError(Form("Can not find run getter in event folder \"%s\"",
+                  GetTitle()));
      return;
    }
   
   AliPHOSLoader* gime = dynamic_cast<AliPHOSLoader*>(runget->GetLoader("PHOSLoader"));
   if ( gime == 0 ) 
    {
-     Error("Exec","Could not obtain the Loader object !"); 
+     AliError("Could not obtain the Loader object !"); 
      return ;
    } 
 
@@ -339,7 +347,7 @@ void  AliPHOSPIDv0::MakeRecParticles(){
 }
 
 //____________________________________________________________________________
-void  AliPHOSPIDv0:: Print() const
+void  AliPHOSPIDv0:: Print(const Option_t *) const
 {
   // Print the parameters used for the particle type identification
   TString message ; 
@@ -351,23 +359,24 @@ void  AliPHOSPIDv0:: Print() const
   message += "    RecParticles Branch title   %s\n" ;  
   message += "with parameters:\n"  ;
   message += "    Maximal EMC - CPV  distance (cm) %f\n" ;
-  Info("Print", message.Data(),  
+  AliInfo(Form( message.Data(),  
        GetTitle(), 
        fRecPointsTitle.Data(), 
        fTrackSegmentsTitle.Data(), 
        fRecParticlesTitle.Data(), 
-       fCpvEmcDistance );
+       fCpvEmcDistance ));
 
   if(fIDOptions.Contains("dis",TString::kIgnoreCase ))
-    Info("Print", "                    dispersion cut %f",  fDispersion ) ;
+    AliInfo(Form("                    dispersion cut %f",  fDispersion )) ;
   if(fIDOptions.Contains("ell",TString::kIgnoreCase ))
-    Info("Print", "             Eliptic cuts function: %s",  fFormula->GetTitle() ) ;
+    AliInfo(Form("             Eliptic cuts function: %s",  
+                fFormula->GetTitle() )) ;
   if(fIDOptions.Contains("tim",TString::kIgnoreCase ))
-    Info("Print",  "             Time Gate used: %f",  fTimeGate) ;
+    AliInfo(Form("             Time Gate used: %f",  fTimeGate)) ;
 }
 
 //____________________________________________________________________________
-void  AliPHOSPIDv0::SetShowerProfileCut(char * formula)
+void  AliPHOSPIDv0::SetShowerProfileCut(const char * formula)
 {
   //set shape of the cut on the axis of ellipce, drown around shouer
   //shower considered "narrow" if Formula(lambda[0],lambda[1]) > 0.
@@ -378,18 +387,20 @@ void  AliPHOSPIDv0::SetShowerProfileCut(char * formula)
 //____________________________________________________________________________
 void  AliPHOSPIDv0::WriteRecParticles()
 {
+  // Saves the reconstructed particles too a file
  
   AliRunLoader* runget = AliRunLoader::GetRunLoader(GetTitle());
   if(runget == 0x0) 
    {
-     Error("Exec","Can not find run getter in event folder \"%s\"",GetTitle());
+     AliError(Form("Can not find run getter in event folder \"%s\"",
+                  GetTitle()));
      return;
    }
   
   AliPHOSLoader* gime = dynamic_cast<AliPHOSLoader*>(runget->GetLoader("PHOSLoader"));
   if ( gime == 0 ) 
    {
-     Error("Exec","Could not obtain the Loader object !"); 
+     AliError("Could not obtain the Loader object !"); 
      return ;
    } 
 
@@ -456,7 +467,7 @@ void  AliPHOSPIDv0::PlotDispersionCuts()const
 }
 
 //____________________________________________________________________________
-TVector3 AliPHOSPIDv0::GetMomentumDirection(AliPHOSEmcRecPoint * emc, AliPHOSRecPoint * cpv)const 
+TVector3 AliPHOSPIDv0::GetMomentumDirection(AliPHOSEmcRecPoint * emc, AliPHOSRecPoint * )const 
 { 
   // Calculates the momentum direction:
   //   1. if only a EMC RecPoint, direction is given by IP and this RecPoint
@@ -508,14 +519,15 @@ void AliPHOSPIDv0::PrintRecParticles(Option_t * option)
   AliRunLoader* runget = AliRunLoader::GetRunLoader(GetTitle());
   if(runget == 0x0) 
    {
-     Error("WriteRecParticles","Can not find run getter in event folder \"%s\"",GetTitle());
+     AliError(Form("Can not find run getter in event folder \"%s\"",
+                  GetTitle()));
      return;
    }
   
   AliPHOSLoader* gime = dynamic_cast<AliPHOSLoader*>(runget->GetLoader("PHOSLoader"));
   if ( gime == 0 ) 
    {
-     Error("WriteRecParticles","Could not obtain the Loader object !"); 
+     AliError("Could not obtain the Loader object !"); 
      return ;
    } 
 
@@ -526,10 +538,11 @@ void AliPHOSPIDv0::PrintRecParticles(Option_t * option)
   TString message ; 
   message  = "event %d\n" ; 
   message += "       found %d RecParticles\n" ; 
-  Info("PrintRecParticles", message.Data(), gAlice->GetEvNumber(), recParticles->GetEntriesFast() ) ;   
+  AliInfo(Form(message.Data(), 
+              gAlice->GetEvNumber(), recParticles->GetEntriesFast() )) ;   
 
   if(strstr(option,"all")) {  // printing found TS
-    Info("PrintRecParticles","  PARTICLE   Index    \n"  ) ; 
+    AliInfo("  PARTICLE   Index"  ) ; 
    
     Int_t index ;
     for (index = 0 ; index < recParticles->GetEntries() ; index++) {
@@ -567,7 +580,8 @@ void AliPHOSPIDv0::PrintRecParticles(Option_t * option)
       //    Int_t nprimaries;
       //    primaries = rp->GetPrimaries(nprimaries);
       
-      Info("PrintRecParticles", "          %s     %d\n",  particle, rp->GetIndexInList()) ;
+      AliInfo(Form("          %s     %d",  
+                  particle, rp->GetIndexInList())) ;
     }
   }  
 }