AliLog is implemented
authorbnandi <bnandi@f7af4fe6-9843-0410-8265-dc069ae4e863>
Thu, 27 Jan 2005 06:26:27 +0000 (06:26 +0000)
committerbnandi <bnandi@f7af4fe6-9843-0410-8265-dc069ae4e863>
Thu, 27 Jan 2005 06:26:27 +0000 (06:26 +0000)
PMD/AliPMD.cxx
PMD/AliPMDDDLRawData.cxx
PMD/AliPMDDigitizer.cxx
PMD/AliPMDDigitizer.h
PMD/AliPMDDiscriminator.cxx
PMD/AliPMDDiscriminator.h
PMD/AliPMDRawStream.cxx
PMD/AliPMDtracker.cxx
PMD/AliPMDtracker.h

index c37f422..e2b12dc 100644 (file)
@@ -45,6 +45,7 @@
 #include <TTree.h>
 #include <TVirtualMC.h>
 
+#include "AliLog.h"
 #include "AliLoader.h" 
 #include "AliPMDLoader.h" 
 #include "AliPMD.h"
@@ -108,11 +109,11 @@ AliLoader* AliPMD::MakeLoader(const char* topfoldername)
  
   if (fLoader)
     {
-      cout<<"Success"<<endl;
+      AliInfo("Success");
     }
   else
     {
-      cout<<"Failure"<<endl;
+      AliError("Failure");
     }
 
   return fLoader;
@@ -319,7 +320,7 @@ void AliPMD::Digits2Raw()
   fLoader->LoadDigits();
   TTree* digits = fLoader->TreeD();
   if (!digits) {
-    Error("Digits2Raw", "no digits tree");
+    AliError("No digits tree");
     return;
   }
 
index d8e290f..847347b 100644 (file)
@@ -19,6 +19,7 @@
 #include <TBranch.h>
 #include <TMath.h>
 
+#include "AliLog.h"
 #include "AliRawDataHeader.h"
 #include "AliBitPacking.h"
 #include "AliPMDdigit.h"
@@ -52,11 +53,15 @@ void AliPMDDDLRawData::WritePMDRawData(TTree *treeD)
   ofstream outfile;
 
   TBranch *branch = treeD->GetBranch("PMDDigit");
-  if (!branch) return;
+  if (!branch)
+    {
+      AliError("PMD Digit branch not found");
+      return;
+    }
   branch->SetAddress(&fDigits);  
   
-  //  Int_t   nmodules = (Int_t) treeD->GetEntries();
-  //  cout << " nmodules = " << nmodules << endl;
+  Int_t   nmodules = (Int_t) treeD->GetEntries();
+  AliDebug(1,Form("Number of modules inside treeD = %d",nmodules));
 
   const Int_t kSize         = 4608;
   const Int_t kDDL          = 6;
index c5fd143..780d72e 100644 (file)
@@ -31,6 +31,7 @@
 #include <TNtuple.h>
 #include <TParticle.h>
 
+#include "AliLog.h"
 #include "AliRun.h"
 #include "AliPMD.h"
 #include "AliPMDhit.h"
@@ -60,7 +61,6 @@ AliPMDDigitizer::AliPMDDigitizer() :
   fSDigits(0),
   fDigits(0),
   fCell(0),
-  fDebug(0),
   fNsdigit(0),
   fNdigit(0),
   fDetNo(0),
@@ -94,7 +94,6 @@ AliPMDDigitizer::AliPMDDigitizer(AliRunDigitizer* manager)
   fSDigits(new TClonesArray("AliPMDsdigit", 1000)),
   fDigits(new TClonesArray("AliPMDdigit", 1000)),
   fCell(0),
-  fDebug(0),
   fNsdigit(0),
   fNdigit(0),
   fDetNo(0),
@@ -152,7 +151,7 @@ void AliPMDDigitizer::OpengAliceFile(const char *file, Option_t *option)
   
   if (!fRunLoader)
    {
-     Error("Open","Can not open session for file %s.",file);
+     AliError(Form("Can not open session for file %s.",file));
    }
   
   if (!fRunLoader->GetAliRun()) fRunLoader->LoadgAlice();
@@ -161,24 +160,20 @@ void AliPMDDigitizer::OpengAliceFile(const char *file, Option_t *option)
 
   gAlice = fRunLoader->GetAliRun();
   
-  if (fDebug) {
-    if (gAlice)
-      {
-       printf("<AliPMDdigitizer::Open> ");
-       printf("AliRun object found on file.\n");
-      }
-    else
-      {
-       printf("<AliPMDdigitizer::Open> ");
-       printf("Could not find AliRun object.\n");
-      }
-  }
+  if (gAlice)
+    {
+      AliDebug(1,"Alirun object found");
+    }
+  else
+    {
+      AliError("Could not found Alirun object");
+    }
   
   fPMD  = (AliPMD*)gAlice->GetDetector("PMD");
   fPMDLoader = fRunLoader->GetLoader("PMDLoader");
   if (fPMDLoader == 0x0)
     {
-      cerr<<"Hits2Digits : Can not find PMD or PMDLoader\n";
+      AliError("Can not find PMDLoader");
     }
 
   const char *cHS = strstr(option,"HS");
@@ -207,7 +202,6 @@ void AliPMDDigitizer::Hits2SDigits(Int_t ievt)
   // This reads the PMD Hits tree and assigns the right track number
   // to a cell and stores in the summable digits tree
   //
-  // cout << " -------- Beginning of Hits2SDigits ----------- " << endl;
 
   const Int_t kPi0 = 111;
   const Int_t kGamma = 22;
@@ -227,9 +221,9 @@ void AliPMDDigitizer::Hits2SDigits(Int_t ievt)
   if (!fSDigits) fSDigits = new TClonesArray("AliPMDsdigit", 1000);
   ResetSDigit();
 
-  if (fDebug) printf("Event Number =  %d \n",ievt); 
+  AliDebug(1,Form("Event Number = %d",ievt));
   Int_t nparticles = fRunLoader->GetHeader()->GetNtrack();
-  if (fDebug) printf("Number of Particles = %d \n", nparticles);
+  AliDebug(1,Form("Number of Particles = %d",nparticles));
   fRunLoader->GetEvent(ievt);
   // ------------------------------------------------------- //
   // Pointer to specific detector hits.
@@ -238,8 +232,7 @@ void AliPMDDigitizer::Hits2SDigits(Int_t ievt)
   TTree* treeH = fPMDLoader->TreeH();
   
   Int_t ntracks    = (Int_t) treeH->GetEntries();
-  if (fDebug) printf("Number of Tracks in the TreeH = %d \n", ntracks);
-
+  AliDebug(1,Form("Number of Tracks in the TreeH = %d", ntracks));
   TTree* treeS = fPMDLoader->TreeS();
   if (treeS == 0x0)
     {
@@ -360,8 +353,7 @@ void AliPMDDigitizer::Hits2SDigits(Int_t ievt)
                  ypad = vol1;
                }
 
-             //cout << "zpos = " << zPos << " edep = " << edep << endl;
-
+             AliDebug(2,Form("Zposition = %f Edeposition = %f",zPos,edep));
              Float_t zposition = TMath::Abs(zPos);
              if (zposition < fZPos)
                {
@@ -434,7 +426,6 @@ void AliPMDDigitizer::Hits2SDigits(Int_t ievt)
   fPMDLoader->WriteSDigits("OVERWRITE");
   ResetCellADC();
 
-  //  cout << " -------- End of Hits2SDigit ----------- " << endl;
 }
 //____________________________________________________________________________
 
@@ -460,10 +451,9 @@ void AliPMDDigitizer::Hits2Digits(Int_t ievt)
   if (!fDigits) fDigits = new TClonesArray("AliPMDdigit", 1000);
   ResetDigit();
 
-  if (fDebug) printf("Event Number =  %d \n",ievt); 
-
+  AliDebug(1,Form("Event Number =  %d",ievt)); 
   Int_t nparticles = fRunLoader->GetHeader()->GetNtrack();
-  if (fDebug) printf("Number of Particles = %d \n", nparticles);
+  AliDebug(1,Form("Number of Particles = %d", nparticles));
   fRunLoader->GetEvent(ievt);
   // ------------------------------------------------------- //
   // Pointer to specific detector hits.
@@ -474,11 +464,11 @@ void AliPMDDigitizer::Hits2Digits(Int_t ievt)
 
   if (fPMDLoader == 0x0)
     {
-      cerr<<"Hits2Digits method : Can not find PMD or PMDLoader\n";
+      AliError("Can not find PMD or PMDLoader");
     }
   TTree* treeH = fPMDLoader->TreeH();
   Int_t ntracks    = (Int_t) treeH->GetEntries();
-  if (fDebug) printf("Number of Tracks in the TreeH = %d \n", ntracks);
+  AliDebug(1,Form("Number of Tracks in the TreeH = %d", ntracks));
   fPMDLoader->LoadDigits("recreate");
   TTree* treeD = fPMDLoader->TreeD();
   if (treeD == 0x0)
@@ -604,8 +594,7 @@ void AliPMDDigitizer::Hits2Digits(Int_t ievt)
                  ypad = vol1;
                }
 
-             //cout << "-zpos = " << -zPos << endl;
-
+             AliDebug(2,Form("ZPosition = %f Edeposition = %d",zPos,edep));
              Float_t zposition = TMath::Abs(zPos);
 
              if (zposition < fZPos)
@@ -682,8 +671,7 @@ void AliPMDDigitizer::Hits2Digits(Int_t ievt)
   
   fPMDLoader->WriteDigits("OVERWRITE");
   ResetCellADC();
-  
-  //  cout << " -------- End of Hits2Digit ----------- " << endl;
+
 }
 //____________________________________________________________________________
 
@@ -693,12 +681,17 @@ void AliPMDDigitizer::SDigits2Digits(Int_t ievt)
   // This reads the PMD sdigits tree and converts energy deposition
   // in a cell to ADC and stores in the digits tree
   //
-  //  cout << " -------- Beginning of SDigits2Digit ----------- " << endl;
+
   fRunLoader->GetEvent(ievt);
 
   TTree* treeS = fPMDLoader->TreeS();
   AliPMDsdigit  *pmdsdigit;
   TBranch *branch = treeS->GetBranch("PMDSDigit");
+  if(!branch)
+    {
+      AliError("PMD Sdigit branch does not exist");
+      return;
+    }
   if (!fSDigits) fSDigits = new TClonesArray("AliPMDsdigit", 1000);
   branch->SetAddress(&fSDigits);
 
@@ -717,12 +710,13 @@ void AliPMDDigitizer::SDigits2Digits(Int_t ievt)
   Float_t edep, adc;
 
   Int_t nmodules = (Int_t) treeS->GetEntries();
+  AliDebug(1,Form("Number of modules = %d",nmodules));
 
   for (Int_t imodule = 0; imodule < nmodules; imodule++)
     {
       treeS->GetEntry(imodule); 
       Int_t nentries = fSDigits->GetLast();
-      //cout << " nentries = " << nentries << endl;
+      AliDebug(2,Form("Number of entries per module = %d",nentries+1));
       for (Int_t ient = 0; ient < nentries+1; ient++)
        {
          pmdsdigit = (AliPMDsdigit*)fSDigits->UncheckedAt(ient);
@@ -740,26 +734,20 @@ void AliPMDDigitizer::SDigits2Digits(Int_t ievt)
       ResetDigit();
     }
   fPMDLoader->WriteDigits("OVERWRITE");
-  //  cout << " -------- End of SDigits2Digit ----------- " << endl;
+
 }
 //____________________________________________________________________________
 void AliPMDDigitizer::Exec(Option_t *option)
 {
   // Does the event merging and digitization
-
-  fDebug = 0;
   const char *cdeb = strstr(option,"deb");
   if(cdeb)
     {
-      cout << "**************** PMD Exec *************** " << endl;
-      fDebug = 1;
+      AliDebug(100," *** PMD Exec is called ***");
     }
-  
+
   Int_t ninputs = fManager->GetNinputs();
-  if(fDebug)
-    {
-      cout << " Number of files = " << ninputs << endl;
-    }
+  AliDebug(1,Form("Number of files to be processed = %d",ninputs));
   ResetCellADC();
 
   for (Int_t i = 0; i < ninputs; i++)
@@ -773,7 +761,7 @@ void AliPMDDigitizer::Exec(Option_t *option)
   fPMDLoader = fRunLoader->GetLoader("PMDLoader");
   if (fPMDLoader == 0x0)
     {
-      cerr<<"AliPMDDigitizer::Exec : Can not find PMD or PMDLoader\n";
+      AliError("Can not find PMD or PMDLoader");
     }
   fPMDLoader->LoadDigits("update");
   TTree* treeD = fPMDLoader->TreeD();
@@ -845,19 +833,13 @@ void AliPMDDigitizer::MergeSDigits(Int_t filenumber, Int_t troffset)
   Float_t edep;
   
   Int_t nmodules = (Int_t) treeS->GetEntries();
-  if(fDebug)
-    {
-      cout << " nmodules = " << nmodules << endl;
-      cout << " tr offset = " << troffset << endl;
-    }
+  AliDebug(1,Form("Number of Modules in the treeS = %d",nmodules));
+  AliDebug(1,Form("Track Offset = %d",troffset));
   for (Int_t imodule = 0; imodule < nmodules; imodule++)
     {
       treeS->GetEntry(imodule); 
       Int_t nentries = fSDigits->GetLast();
-      if(fDebug)
-       {
-         cout << " nentries = " << nentries << endl;
-       }
+      AliDebug(2,Form("Number of Entries per Module = %d",nentries));
       for (Int_t ient = 0; ient < nentries+1; ient++)
        {
          pmdsdigit = (AliPMDsdigit*)fSDigits->UncheckedAt(ient);
@@ -964,7 +946,6 @@ void AliPMDDigitizer::TrackAssignment2Cell()
       iyp      = cell->GetY();
       edep     = cell->GetEdep();
       Int_t nn = fPRECounter[ism][ixp][iyp];
-      //      cout << " nn = " << nn << endl;
       pmdTrack[ism][ixp][iyp][nn] = (Int_t) mtrackno;
       pmdEdep[ism][ixp][iyp][nn] = edep;
       fPRECounter[ism][ixp][iyp]++;
@@ -1099,10 +1080,6 @@ void AliPMDDigitizer::MeV2ADC(Float_t mev, Float_t & adc) const
   // To be done
   //
 
-  //  adc = mev*1.;
-
-  
-  // %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% //
   // PS Test in September 2003
   // MeV - ADC conversion for 10bit ADC
 
@@ -1131,8 +1108,6 @@ void AliPMDDigitizer::MeV2ADC(Float_t mev, Float_t & adc) const
       adc = 3000.0;
     }
 
-  // %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% //
-
 }
 //____________________________________________________________________________
 void AliPMDDigitizer::AddSDigit(Int_t trnumber, Int_t det, Int_t smnumber, 
index 3205a38..39b17a8 100644 (file)
@@ -75,7 +75,6 @@ class AliPMDDigitizer:public AliDigitizer
 
   TObjArray     fCell;       //! List of pmd cells
 
-  Int_t   fDebug;            // Debug switch
   Int_t   fNsdigit;          // Summable digits counter
   Int_t   fNdigit;           // Digits counter
   Int_t   fDetNo;            // Detector Number (0:PRE, 1:CPV)
@@ -92,7 +91,7 @@ class AliPMDDigitizer:public AliDigitizer
   Int_t   fCPVTrackNo[fgkTotUM][fgkRow][fgkCol]; //! CPV Array containing track number
 
 
-  ClassDef(AliPMDDigitizer,4)    // To digitize PMD Hits
+  ClassDef(AliPMDDigitizer,5)    // To digitize PMD Hits
 };
 #endif
 
index 25717d5..5f16505 100644 (file)
@@ -34,6 +34,7 @@
 #include <TNtuple.h>
 #include <TParticle.h>
 
+#include "AliLog.h"
 #include "AliPMDcluster.h"
 #include "AliPMDclupid.h"
 #include "AliPMDDiscriminator.h"
@@ -41,8 +42,7 @@
 ClassImp(AliPMDDiscriminator)
 
 AliPMDDiscriminator::AliPMDDiscriminator() :
-  fDiscrim(0),
-  fDebug(0)
+  fDiscrim(0)
 {
 //
 // Default Constructor
@@ -94,7 +94,8 @@ void AliPMDDiscriminator::EmpDiscrimination(TObjArray *pmdcontin, TObjArray *pmd
   AliPMDclupid   *pmdclout = 0;
 
   Int_t nentries1 = pmdcontin->GetEntries();
-  cout << " nentries1 = " << nentries1 << endl;
+
+  AliDebug(1,Form("Number of total clusters from CPV PRE = %d",nentries1));
   for (Int_t ient1 = 0; ient1 < nentries1; ient1++)
     {
       pmdcl = (AliPMDcluster*)pmdcontin->UncheckedAt(ient1);
index 214bee1..0f03271 100644 (file)
@@ -18,7 +18,7 @@ class TNtuple;
 class AliPMDcluster;
 class AliPMDclupid;
 
-class AliPMDDiscriminator
+class AliPMDDiscriminator : public TObject
 {
 
  public:
@@ -35,9 +35,8 @@ class AliPMDDiscriminator
  protected:
 
   Int_t   fDiscrim;       // To switch on different discrimination method
-  Int_t   fDebug;         // Debugging switch (0:NO, 1:YES)
 
-  ClassDef(AliPMDDiscriminator,1) // To run PMD discrimination
+  ClassDef(AliPMDDiscriminator,2) // To run PMD discrimination
 };
 #endif
 
index 8f157e8..0f16529 100644 (file)
@@ -26,6 +26,7 @@
 ///
 ///////////////////////////////////////////////////////////////////////////////
 
+#include "AliLog.h"
 #include "AliPMDRawStream.h"
 #include "AliRawReader.h"
 
@@ -66,7 +67,7 @@ AliPMDRawStream::AliPMDRawStream(const AliPMDRawStream& stream) :
 {
 // copy constructor
 
-  Fatal("AliPMDRawStream", "copy constructor not implemented");
+  AliFatal("Copy constructor not implemented");
 }
 
 //_____________________________________________________________________________
@@ -75,7 +76,7 @@ AliPMDRawStream& AliPMDRawStream::operator = (const AliPMDRawStream&
 {
 // assignment operator
 
-  Fatal("operator =", "assignment operator not implemented");
+  AliFatal("operator = assignment operator not implemented");
   return *this;
 }
 
index 4f264a6..8e128b8 100644 (file)
@@ -44,7 +44,7 @@
 
 #include "AliESDPmdTrack.h"
 #include "AliESD.h"
-
+#include "AliLog.h"
 
 ClassImp(AliPMDtracker)
 
@@ -58,7 +58,6 @@ AliPMDtracker::AliPMDtracker():
   fPMDrecpoint(0),
   fPMDclin(0),
   fPMDclout(0),
-  fDebug(0),
   fXvertex(0.),
   fYvertex(0.),
   fZvertex(0.),
@@ -111,16 +110,21 @@ void AliPMDtracker::Clusters2Tracks(AliESD *event)
   Float_t clusdata[5];
 
   TBranch *branch = fTreeR->GetBranch("PMDRecpoint");
-  if (!branch) return;
+  if (!branch)
+    {
+      AliError("PMDRecpoint branch not found");
+      return;
+    }
   branch->SetAddress(&fRecpoints);  
   
   Int_t   nmodules = (Int_t) fTreeR->GetEntries();
-  cout << " nmodules = " << nmodules << endl;
+  AliDebug(1,Form("Number of modules filled in treeR = %d",nmodules));
   for (Int_t imodule = 0; imodule < nmodules; imodule++)
     {
       fTreeR->GetEntry(imodule); 
       Int_t nentries = fRecpoints->GetLast();
-      //      cout << " nentries = " << nentries << endl;
+      AliDebug(2,Form("Number of clusters per modules filled in treeR = %d"
+                     ,nentries));
       for(Int_t ient = 0; ient < nentries+1; ient++)
        {
          fPMDrecpoint = (AliPMDrecpoint1*)fRecpoints->UncheckedAt(ient);
@@ -151,7 +155,8 @@ void AliPMDtracker::Clusters2Tracks(AliESD *event)
   Float_t zglobal = kzpos + (Float_t) fZvertex;
 
   Int_t nentries2 = fPMDcontout->GetEntries();
-  cout << " nentries2 = " << nentries2 << endl;
+  AliDebug(1,Form("Number of clusters coming after discrimination = %d"
+                 ,nentries2));
   for (Int_t ient1 = 0; ient1 < nentries2; ient1++)
     {
       fPMDclout = (AliPMDclupid*)fPMDcontout->UncheckedAt(ient1);
@@ -236,11 +241,6 @@ void AliPMDtracker::SetVertex(Double_t vtx[3], Double_t evtx[3])
   fSigmaZ  = evtx[2];
 }
 //--------------------------------------------------------------------//
-void AliPMDtracker::SetDebug(Int_t idebug)
-{
-  fDebug = idebug;
-}
-//--------------------------------------------------------------------//
 void AliPMDtracker::ResetClusters()
 {
   if (fRecpoints) fRecpoints->Clear();
index c0acb15..9e9eaf2 100644 (file)
@@ -23,7 +23,7 @@ class AliPMDDiscriminator;
 
 class AliESD;
 
-class AliPMDtracker
+class AliPMDtracker:public TObject
 {
 
  public:
@@ -34,7 +34,6 @@ class AliPMDtracker
   void LoadClusters(TTree *treein);
   void Clusters2Tracks(AliESD *event);
   void SetVertex(Double_t vtx[3], Double_t evtx[3]);
-  void SetDebug(Int_t idebug);
   void ResetClusters();
 
  protected:
@@ -50,7 +49,6 @@ class AliPMDtracker
   AliPMDcluster       *fPMDclin;
   AliPMDclupid        *fPMDclout;
 
-  Int_t    fDebug;          // Debugging switch (0:NO, 1:YES)
   Double_t fXvertex;        // X-vertex position
   Double_t fYvertex;        // Y-vertex position
   Double_t fZvertex;        // Z-vertex position
@@ -58,7 +56,7 @@ class AliPMDtracker
   Double_t fSigmaY;         // Y-vertex error
   Double_t fSigmaZ;         // Z-vertex error
 
-  ClassDef(AliPMDtracker,1) // To run PMD clustering
+  ClassDef(AliPMDtracker,2) // To run PMD clustering
 };
 #endif