Updating geometry for new scintillaton modules

[u/mrichter/AliRoot.git] / EMCAL / AliEMCALReconstruct.C

// YS Subatech Mai 2002
// YK Subatech 6 Aug 2002

// EMCAL Reconstruction chain:
// Hits -> SDigits -> Digits -> RecPoints -> TrackSegments -> RecParticles

//Root
#if !defined(__CINT__) || defined(__MAKECINT__)
#include "TString.h"

//AliRoot
#include "STEER/AliRun.h"
#include "EMCAL/AliEMCALSDigitizer.h"
#include "EMCAL/AliEMCALDigitizer.h"
#include "EMCAL/AliEMCALClusterizerv1.h"
// #include "EMCAL/AliEMCALTrackSegmentMakerv1.h"
// #include "EMCAL/AliEMCALPIDv1.h"
#endif

void EMCALHits2SDigits( Bool_t split=kFALSE, TString fileName = "galice.root") {

  // usage : 
  // 1. write SDigits in the same file as Hits --------------- (OK)
  //root [0] .L Reconstruct.C++
  //root [1] SDigits2Digits()
  // 2. write SDigits in a separate file, one per detector, from Hits --------------- (OK)
  //root [0] .L Reconstruct.C++
  //root [1] SDigits2Digits(kTRUE) // SDigits saved in [DET}.SDigits.root (DET=EMCAL, EMCAL)

  delete gAlice ; 
  gAlice = 0 ; 
  
  AliEMCALSDigitizer * sdp = new AliEMCALSDigitizer(fileName) ; 
  if (split) 
    sdp->SetSplitFile() ;
  sdp->ExecuteTask("deb") ; 

  delete sdp ;
}

//________________________________________________________________________
void EMCALSDigits2Digits( Bool_t split=kFALSE, TString fileName = "galice.root") {
  
 // usage : 
  // 1. write SDigits in the same file as SDigits --------------- (OK)
  //root [0] .L Reconstruct.C++
  //root [1] SDigits2Digits()
  // 2. write SDigits in a separate file, one per detector, from SDigits --------------- (OK)
  //root [0] .L Reconstruct.C++
  //root [1] SDigitsDigits(kTRUE) // Digits saved in [DET}.Digits.root (DET=EMCAL, EMCAL)

  delete gAlice ; 
  gAlice = 0 ; 
  
  // EMCAL
  AliEMCALDigitizer * dp = 0 ; 
 
  if (split) {
    dp = new AliEMCALDigitizer("EMCAL.SDigits.root") ; 
    dp->SetSplitFile() ; } 
  else 
    dp = new AliEMCALDigitizer(fileName) ; 
  
  dp->ExecuteTask("deb") ; 
  
  delete dp ;
}

//________________________________________________________________________
void EMCALDigits2RecPoints( Bool_t split=kFALSE, TString fileName = "galice.root") {
  
 // usage : 
  // 1. write RecPoints in the same file as Digits --------------- OK 
  //root [0] .L Reconstruct.C++
  //root [1] Digits2RecPoints()
  // 2. write RecPoints in a separate file, one per detector, from Digits --------------- OK 
  //root [0] .L Reconstruct.C++
  //root [1] Digits2RecPoints(kTRUE) // RecPoints saved in [DET}.RecPoints.root (DET=EMCAL, EMCAL)

  delete gAlice ; 
  gAlice = 0 ; 
 
  AliEMCALClusterizer * cp = 0 ; 
 
  if (split) {
    cp = new AliEMCALClusterizerv1("EMCAL.Digits.root") ; 
    cp->SetSplitFile() ; } 
  else 
    cp = new AliEMCALClusterizerv1(fileName) ; 
  
  cp->ExecuteTask("deb") ; 
  
  delete cp ;
}

// //________________________________________________________________________
// void EMCALRecPoints2TrackSegments( Bool_t split=kFALSE, TString fileName = "galice.root") {
  
//  // usage : 
//   // 1. write TrackSegments in the same file as RecPoints --------------- (OK) 
//   //root [0] .L Reconstruct.C++
//   //root [1] RecPoints2TrackSegments()
//   // 2. write TrackSegments in a separate file, one per detector, from RecPoints --------------- (Not needed) 
//   //root [0] .L Reconstruct.C++
//   //root [1] RecPoints2TrackSegments(kTRUE) // TrackSegments saved in [DET}.RecData.root (DET=EMCAL, EMCAL)

//   delete gAlice ; 
//   gAlice = 0 ; 
  
//   AliEMCALTrackSegmentMaker * tmp = 0 ; 
 
//   if (split)
//     tmp = new AliEMCALTrackSegmentMakerv1("EMCAL.RecData.root") ; 
//   else 
//     tmp = new AliEMCALTrackSegmentMakerv1(fileName) ; 
  
//   tmp->ExecuteTask("deb") ; 
  
//   delete tmp ;
// }

// //________________________________________________________________________
// void EMCALTrackSegments2RecParticles( Bool_t split=kFALSE, TString fileName = "galice.root") {
  
//  // usage : 
//   // 1. write RecParticles in the same file as TrackSegments ---------------  (OK)
//   //root [0] .L Reconstruct.C++
//   //root [1] TrackSegments2RecParticles()
//   // 2. write RecParticles in a separate file, one per detector, from TrackSegments --------------- (Not needed) 
//   //root [0] .L Reconstruct.C++
//   //root [1] TrackSegments2RecParticles(kTRUE) // RecParticles saved in [DET}.RecData.root (DET=EMCAL, EMCAL)

//   delete gAlice ; 
//   gAlice = 0 ; 
  
//   AliEMCALPID * pp = 0 ; 
 
//   if (split) 
//     pp = new AliEMCALPIDv1("EMCAL.RecData.root") ; 
//   else 
//     pp = new AliEMCALPIDv1(fileName) ; 
  
//   pp->ExecuteTask("deb") ; 
  
//   delete pp ;
// }

// //________________________________________________________________________
// void EMCALDigits2RecParticles( Bool_t split=kFALSE, TString fileName = "galice.root") {
  
//  // usage : 
//   // 1. write RecPoints, TrackSegments and RecParticles in the same file as Digits --------------- (OK)
//   //root [0] .L Reconstruct.C++
//   //root [1] Digits2RecParticles()
//   // 2. write RecPoints , TrackSegments and RecParticles in a separate file, one per detector, from Digits --------------- (OK)
//   //root [0] .L Reconstruct.C++
//   //root [1] Digits2RecParticles(kTRUE) // TrackSegments saved in [DET}.RecData.root (DET=EMCAL, EMCAL)

 
//   delete gAlice ; 
//   gAlice = 0 ; 
 
//   // EMCAL
//   AliEMCALClusterizer * cp = 0 ; 
 
//   if (split) {
//     cp = new AliEMCALClusterizerv1("EMCAL.Digits.root") ; 
//     cp->SetSplitFile() ; } 
//   else 
//     cp = new AliEMCALClusterizerv1(fileName) ; 
  
//   cp->ExecuteTask("deb") ; 

//   if (split) 
//     delete cp ;
  
//   AliEMCALTrackSegmentMaker * tmp = 0 ; 
  
//   if (split) 
//     tmp = new AliEMCALTrackSegmentMakerv1("EMCAL.RecData.root") ; 
//   else 
//     tmp = new AliEMCALTrackSegmentMakerv1(fileName) ; 
  
//   tmp->ExecuteTask("deb") ; 
  
//   AliEMCALPID * pp = 0 ; 
 
//   if (split) 
//     pp = new AliEMCALPIDv1("EMCAL.RecData.root") ; 
//   else 
//     pp = new AliEMCALPIDv1(fileName) ; 
  
//   pp->ExecuteTask("deb") ; 
  
//   delete tmp; 
//   delete pp ; 
// }

//________________________________________________________________________
void EMCALHits2Digits (Bool_t split=kFALSE, TString fileName = "galice.root") {
  // usage : 
  // 1. write (S)Digits in the same file as Hits --------------- (OK)
  //root [0] .L Reconstruct.C++
  //root [1] Hits2Digits()
  // 2. write (S)Digits in a separate file, one per detector, from Hits --------------- (OK)
  //root [0] .L Reconstruct.C++
  //root [1] Hits2Digits(kTRUE) // SDigits saved in [DET}.SDigits.root (DET=EMCAL, EMCAL)
                                // Digits  saved in [DET}.Digits.root  (DET=EMCAL, EMCAL)

  delete gAlice ; 
  gAlice = 0 ; 
  
  //EMCAL
  AliEMCALSDigitizer * sdp = new AliEMCALSDigitizer(fileName) ; 
  if (split) 
    sdp->SetSplitFile() ;
  sdp->ExecuteTask("deb") ; 

  if (split) 
    delete sdp ; 

  AliEMCALDigitizer * dp = 0 ; 
 
  if (split) {
    dp = new AliEMCALDigitizer("EMCAL.SDigits.root") ; 
    dp->SetSplitFile() ; } 
  else 
    dp = new AliEMCALDigitizer(fileName) ; 
  
  dp->ExecuteTask("deb") ; 

  if (split) 
    delete dp ; 

  if (!split) { 
    delete sdp ; 
    delete dp ; 
  }
}