X-Git-Url: http://git.uio.no/git/?a=blobdiff_plain;f=PHOS%2FAliPHOSPIDv0.cxx;h=2516937dbeb8f21e3901fe19831de2f1ec17ceb1;hb=75ab5a703d0b8475072def28f94a9462bc48d01f;hp=9479caa6aeb70b40583d6036c56301397e05a65d;hpb=88cb7938ca21d4a80991d4e7aa564008c29340f7;p=u%2Fmrichter%2FAliRoot.git diff --git a/PHOS/AliPHOSPIDv0.cxx b/PHOS/AliPHOSPIDv0.cxx index 9479caa6aeb..2516937dbeb 100644 --- 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(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(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(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(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(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())) ; } } }