PHOS/Reconstruction.C

   1 #include "AliReconstruction.h"
   2 #include "TString.h"
   3 #include "Riostream.h"
   4 #include "AliPHOSGetter.h"
   5 #include "AliEMCALGetter.h"
   6
   7 void reco(TString opt, TString name)
   8 {
   9   AliReconstruction rec ;
  10   if ( !opt.Contains("T") )
  11     rec.SetRunTracking(kFALSE) ;
  12   if ( opt.Contains("R") )
  13     rec.SetRunReconstruction(name.Data()) ;
  14   if ( !opt.Contains("E") )
  15     rec.SetFillESD("") ;
  16   else
  17     rec.SetFillESD(name.Data()) ;
  18   rec.Run() ;
  19
  20   if ( name.Contains("PHOS") ) {
  21     cout << ">>>>>>>>>>>> PHOS " << endl ;
  22     AliPHOSGetter * gime = AliPHOSGetter::Instance("galice.root") ;
  23     Int_t event ;
  24     for (event = 0; event < gime->MaxEvent(); event++) {
  25       cout << "event # " << event << endl ;
  26       gime->Event(event, "RP") ;
  27       cout << "   EMC RecPoints  # " << gime->EmcRecPoints()->GetEntries() << endl ;
  28       cout << "   CPV RecPoints  # " << gime->CpvRecPoints()->GetEntries() << endl ;
  29       cout << "   Track Segments # " << gime->TrackSegments()->GetEntries() << endl ;
  30       cout << "   Rec Particles  # " << gime->RecParticles()->GetEntries() << endl ;
  31     }
  32   }
  33  if ( name.Contains("EMCAL") ) {
  34     cout << ">>>>>>>>>>>> EMCAL " << endl ;
  35     AliEMCALGetter * gime = AliEMCALGetter::Instance("galice.root") ;
  36     Int_t event ;
  37     for (event = 0; event < gime->MaxEvent(); event++) {
  38       cout << "event # " << event << endl ;
  39       gime->Event(event, "RP") ;
  40       cout << "       RecPoints  # " << gime->ECARecPoints()->GetEntries() << endl ;
  41       cout << "   Rec Particles  # " << gime->RecParticles()->GetEntries() << endl ;
  42     }
  43  }
  44 }