Adding statistical function (Marian)
[u/mrichter/AliRoot.git] / PHOS / AliPHOSSDigitizer.cxx
index a5c2d58..bcffe97 100644 (file)
 
 /* $Id$ */
 
+/* History of cvs commits:
+ *
+ * $Log$
+ * Revision 1.49  2006/05/10 06:42:53  kharlov
+ * Remove redundant loop over primaries
+ *
+ * Revision 1.48  2006/04/22 10:30:17  hristov
+ * Add fEnergy to AliPHOSDigit and operate with EMC amplitude in energy units (Yu.Kharlov)
+ *
+ * Revision 1.47  2005/05/28 14:19:04  schutz
+ * Compilation warnings fixed by T.P.
+ *
+ */
+
 //_________________________________________________________________________
 // This is a TTask that makes SDigits out of Hits
 // The name of the TTask is also the title of the branch that will contain 
 
 
 // --- ROOT system ---
-#include "TFile.h"
-#include "TROOT.h"
 #include "TBenchmark.h"
+               //#include "TObjectTable.h"
 
 // --- Standard library ---
 
 // --- AliRoot header files ---
-#include "AliRun.h"
+#include "AliLog.h"
+#include "AliPHOSGeometry.h" 
 #include "AliPHOSDigit.h"
 #include "AliPHOSGetter.h"
 #include "AliPHOSHit.h"
 #include "AliPHOSSDigitizer.h"
+//#include "AliMemoryWatcher.h"
 
 ClassImp(AliPHOSSDigitizer)
 
            
 //____________________________________________________________________________ 
-  AliPHOSSDigitizer::AliPHOSSDigitizer():TTask("","")
+AliPHOSSDigitizer::AliPHOSSDigitizer() : 
+  TTask("",""),
+  fA(0.f), fB(0.f),
+  fPrimThreshold(0.f),
+  fDefaultInit(kTRUE),
+  fEventFolderName(""),
+  fInit(kFALSE),
+  fSDigitsInRun(0),
+  fFirstEvent(0),
+  fLastEvent(0)
 {
   // ctor
-  InitParameters() ;
-  fDefaultInit = kTRUE ; 
 }
 
 //____________________________________________________________________________ 
-AliPHOSSDigitizer::AliPHOSSDigitizer(const char * alirunFileName, const char * eventFolderName):
-  TTask("PHOS"+AliConfig::fgkSDigitizerTaskName, alirunFileName),
-  fEventFolderName(eventFolderName)
+AliPHOSSDigitizer::AliPHOSSDigitizer(const char * alirunFileName, 
+                                    const char * eventFolderName):
+  TTask("PHOS"+AliConfig::Instance()->GetSDigitizerTaskName(), alirunFileName),
+  fA(0.f), fB(0.f),
+  fPrimThreshold(0.f),
+  fDefaultInit(kFALSE),
+  fEventFolderName(eventFolderName),
+  fInit(kFALSE),
+  fSDigitsInRun(0),
+  fFirstEvent(0),
+  fLastEvent(0)
 {
-
   // ctor
   InitParameters() ; 
   Init();
@@ -88,19 +117,14 @@ AliPHOSSDigitizer::AliPHOSSDigitizer(const char * alirunFileName, const char * e
 }
 
 //____________________________________________________________________________ 
-AliPHOSSDigitizer::AliPHOSSDigitizer(const AliPHOSSDigitizer & sd) {
-  //cpy ctor 
-
-  fA             = sd.fA ;
-  fB             = sd.fB ;
-  fPrimThreshold = sd.fPrimThreshold ;
-  fSDigitsInRun  = sd.fSDigitsInRun ;
-  SetName(sd.GetName()) ; 
-  SetTitle(sd.GetTitle()) ; 
-  fEventFolderName = sd.fEventFolderName;
+AliPHOSSDigitizer::~AliPHOSSDigitizer() {
+  //dtor
+  //  AliPHOSGetter * gime =
+  //  AliPHOSGetter::Instance(GetTitle(),fEventFolderName.Data());  
+  AliPHOSGetter * gime =
+    AliPHOSGetter::Instance();  
+  gime->PhosLoader()->CleanSDigitizer();
 }
-
-
 //____________________________________________________________________________ 
 void AliPHOSSDigitizer::Init()
 {
@@ -122,7 +146,7 @@ void AliPHOSSDigitizer::Init()
 
   gime->PostSDigitizer(this);
   gime->PhosLoader()->GetSDigitsDataLoader()->GetBaseTaskLoader()->SetDoNotReload(kTRUE);
-  
 }
 
 //____________________________________________________________________________ 
@@ -138,7 +162,13 @@ void AliPHOSSDigitizer::InitParameters()
 //____________________________________________________________________________
 void AliPHOSSDigitizer::Exec(Option_t *option) 
 { 
-  // Collects all hits in the same active volume into digit
+  // Steering method to produce summable digits for events
+  // in the range from fFirstEvent to fLastEvent.
+  // This range is optionally set by SetEventRange().
+  // if fLastEvent=-1 (by default), then process events until the end.
+  //
+  // Summable digit is a sum of all hits in the same active
+  // volume into digit
   
   if (strstr(option, "print") ) {
     Print() ; 
@@ -152,48 +182,42 @@ void AliPHOSSDigitizer::Exec(Option_t *option)
 
   //switch off reloading of this task while getting event
   if (!fInit) { // to prevent overwrite existing file
-    Error( "Exec", "Give a version name different from %s", fEventFolderName.Data() ) ;
+    AliError( Form("Give a version name different from %s", fEventFolderName.Data()) ) ;
     return ;
   }
 
+  if (fLastEvent == -1) 
+    fLastEvent = gime->MaxEvent() - 1 ;
+  else 
+    fLastEvent = TMath::Min(fFirstEvent, gime->MaxEvent()); // only ine event at the time 
+  Int_t nEvents   = fLastEvent - fFirstEvent + 1;
+  
+  Int_t ievent, i;
 
-  Int_t nevents = gime->MaxEvent() ; 
-  Int_t ievent ;
-  for(ievent = 0; ievent < nevents; ievent++){
-
-    gime->Event(ievent,"H") ;
+  //AliMemoryWatcher memwatcher;
 
+  for (ievent = fFirstEvent; ievent <= fLastEvent; ievent++) {    
+    gime->Event(ievent,"H") ;   
     TTree * treeS = gime->TreeS(); 
     TClonesArray * hits = gime->Hits() ;
     TClonesArray * sdigits = gime->SDigits() ;
     sdigits->Clear();
     Int_t nSdigits = 0 ;
     //Now make SDigits from hits, for PHOS it is the same, so just copy    
-    Int_t nPrim =  static_cast<Int_t>((gime->TreeH())->GetEntries()) ; 
-    // Attention nPrim is the number of primaries tracked by Geant 
-    // and this number could be different to the number of Primaries in TreeK;
-    Int_t iprim ;
-
-    for (iprim = 0 ; iprim < nPrim ; iprim ++) { 
-      //=========== Get the PHOS branch from Hits Tree for the Primary iprim
-      gime->Track(iprim) ;
-      Int_t i;
-      for ( i = 0 ; i < hits->GetEntries() ; i++ ) {
-       AliPHOSHit * hit = dynamic_cast<AliPHOSHit *>(hits->At(i)) ;
-       // Assign primary number only if contribution is significant
-       
-       if( hit->GetEnergy() > fPrimThreshold)
-         new((*sdigits)[nSdigits]) AliPHOSDigit(hit->GetPrimary(),hit->GetId(),
-                                                Digitize(hit->GetEnergy()), hit->GetTime()) ;
-       else
-         new((*sdigits)[nSdigits]) AliPHOSDigit( -1              , hit->GetId(), 
-                                                 Digitize(hit->GetEnergy()), hit->GetTime()) ;
-       nSdigits++ ;    
-       
-      }
+    for ( i = 0 ; i < hits->GetEntries() ; i++ ) {
+      AliPHOSHit * hit = dynamic_cast<AliPHOSHit *>(hits->At(i)) ;
+      // Assign primary number only if contribution is significant
+      
+      if( hit->GetEnergy() > fPrimThreshold)
+       new((*sdigits)[nSdigits]) AliPHOSDigit(hit->GetPrimary(),hit->GetId(),
+                                              hit->GetEnergy() ,hit->GetTime()) ;
+      else
+       new((*sdigits)[nSdigits]) AliPHOSDigit(-1               ,hit->GetId(), 
+                                              hit->GetEnergy() ,hit->GetTime()) ;
+      nSdigits++ ;     
+      
+    }
  
-    } // loop over iprim
-
     sdigits->Sort() ;
 
     nSdigits = sdigits->GetEntriesFast() ;
@@ -201,7 +225,6 @@ void AliPHOSSDigitizer::Exec(Option_t *option)
     fSDigitsInRun += nSdigits ;  
     sdigits->Expand(nSdigits) ;
 
-    Int_t i ;
     for (i = 0 ; i < nSdigits ; i++) { 
       AliPHOSDigit * digit = dynamic_cast<AliPHOSDigit *>(sdigits->At(i)) ; 
       digit->SetIndexInList(i) ;     
@@ -222,24 +245,31 @@ void AliPHOSSDigitizer::Exec(Option_t *option)
     //Next - SDigitizer
 
     gime->WriteSDigitizer("OVERWRITE");
-
+    //gObjectTable->Print() ; 
+  
     if(strstr(option,"deb"))
       PrintSDigits(option) ;
-  }
+      
+    //memwatcher.Watch(ievent); 
+  }// event loop
   
   Unload();
 
-  gime->PhosLoader()->GetSDigitsDataLoader()->GetBaseTaskLoader()->SetDoNotReload(kTRUE);
+  //  gime->PhosLoader()->GetSDigitsDataLoader()->GetBaseTaskLoader()->SetDoNotReload(kTRUE);
   
   if(strstr(option,"tim")){
     gBenchmark->Stop("PHOSSDigitizer");
     Info("Exec","   took %f seconds for SDigitizing  %f seconds per event",
-        gBenchmark->GetCpuTime("PHOSSDigitizer"), gBenchmark->GetCpuTime("PHOSSDigitizer")/nevents) ;
+        gBenchmark->GetCpuTime("PHOSSDigitizer"), gBenchmark->GetCpuTime("PHOSSDigitizer")/nEvents) ;
   }
+  
+  //TFile f("out.root","RECREATE");
+  //memwatcher.WriteToFile(); 
+  //f.Close();
 }
 
 //__________________________________________________________________
-void AliPHOSSDigitizer::Print()const
+void AliPHOSSDigitizer::Print(const Option_t *)const
 {
   // Prints parameters of SDigitizer
   Info("Print", "\n------------------- %s -------------", GetName() ) ; 
@@ -271,34 +301,27 @@ void AliPHOSSDigitizer::PrintSDigits(Option_t * option)
 
   AliPHOSGetter * gime = AliPHOSGetter::Instance() ; 
   const TClonesArray * sdigits = gime->SDigits() ;
-
-  TString message ; 
-  message  = "\nAliPHOSSDigitiser: event " ;
-  message += gAlice->GetEvNumber(); 
-  message += "\n      Number of entries in SDigits list " ;  
-  message += sdigits->GetEntriesFast() ; 
-  char * tempo = new char[8192]; 
   
+  Info( "\nPrintSDigits", "event # %d %d sdigits", gAlice->GetEvNumber(), sdigits->GetEntriesFast() ) ; 
+
   if(strstr(option,"all")||strstr(option,"EMC")){
     
     //loop over digits
     AliPHOSDigit * digit;
-    message += "\nEMC sdigits\n" ;
-    message += "\n   Id  Amplitude    Time          Index Nprim: Primaries list \n" ;    
+    printf("\nEMC sdigits\n") ; 
     Int_t maxEmc = gime->PHOSGeometry()->GetNModules()*gime->PHOSGeometry()->GetNCristalsInModule() ;
     Int_t index ;
     for (index = 0 ; (index < sdigits->GetEntriesFast()) && 
         ((dynamic_cast<AliPHOSDigit *> (sdigits->At(index)))->GetId() <= maxEmc) ; index++) {
       digit = dynamic_cast<AliPHOSDigit *>( sdigits->At(index) ) ;
-      if(digit->GetNprimary() == 0) 
-       continue;
-      sprintf(tempo, "\n%6d  %8d    %6.5e %4d      %2d :",
-             digit->GetId(), digit->GetAmp(), digit->GetTime(), digit->GetIndexInList(), digit->GetNprimary()) ;  
-      message += tempo ; 
+      //  if(digit->GetNprimary() == 0) 
+      //       continue;
+//       printf("%6d  %8d    %6.5e %4d      %2d :\n", // YVK
+      printf("%6d  %.4f    %6.5e %4d      %2d :\n",
+             digit->GetId(), digit->GetEnergy(), digit->GetTime(), digit->GetIndexInList(), digit->GetNprimary()) ;  
       Int_t iprimary;
       for (iprimary=0; iprimary<digit->GetNprimary(); iprimary++) {
-       sprintf(tempo, "%d ",digit->GetPrimary(iprimary+1) ) ; 
-       message += tempo ; 
+       printf("%d ",digit->GetPrimary(iprimary+1) ) ; 
       }  
     }    
   }
@@ -307,32 +330,27 @@ void AliPHOSSDigitizer::PrintSDigits(Option_t * option)
     
     //loop over CPV digits
     AliPHOSDigit * digit;
-    
-    message += "CPV sdigits\n" ;
-    message += "\n   Id  Amplitude    Index Nprim: Primaries list \n" ;    
+    printf("\nCPV sdigits\n") ; 
     Int_t maxEmc = gime->PHOSGeometry()->GetNModules()*gime->PHOSGeometry()->GetNCristalsInModule() ;
     Int_t index ;
     for (index = 0 ; index < sdigits->GetEntriesFast(); index++) {
       digit = dynamic_cast<AliPHOSDigit *>( sdigits->At(index) ) ;
       if(digit->GetId() > maxEmc){
-       sprintf(tempo, "\n%6d  %8d    %4d      %2d :",
+       printf("\n%6d  %8d    %4d      %2d :",
                digit->GetId(), digit->GetAmp(), digit->GetIndexInList(), digit->GetNprimary()) ;  
-       message += tempo ; 
        Int_t iprimary;
        for (iprimary=0; iprimary<digit->GetNprimary(); iprimary++) {
-         sprintf(tempo, "%d ",digit->GetPrimary(iprimary+1) ) ; 
-         message += tempo ; 
+         printf("%d ",digit->GetPrimary(iprimary+1) ) ; 
        }
       }    
     }
   }
-  delete []tempo ; 
-  Info("PrintSDigits", message.Data() ) ;
 }
 
 //____________________________________________________________________________ 
 void AliPHOSSDigitizer::Unload() const
 {
+  // Unloads the objects from the folder
   AliPHOSGetter * gime = AliPHOSGetter::Instance() ; 
   AliPHOSLoader * loader = gime->PhosLoader() ; 
   loader->UnloadHits() ;