PHOS/AliPHOSReconstruct.C

   1 // YS Subatech Mai 2002
   2 // YK Subatech 6 Aug 2002
   3
   4 // PHOS Reconstruction chain:
   5 // Hits -> SDigits -> Digits -> RecPoints -> TrackSegments -> RecParticles
   6
   7 //Root
   8 #include "TString.h"
   9
  10 //AliRoot
  11 #include "STEER/AliRun.h"
  12 #include "PHOS/AliPHOSSDigitizer.h"
  13 #include "PHOS/AliPHOSDigitizer.h"
  14 #include "PHOS/AliPHOSClusterizerv1.h"
  15 #include "PHOS/AliPHOSTrackSegmentMakerv1.h"
  16 #include "PHOS/AliPHOSPIDv1.h"
  17 #include "EMCAL/AliEMCALSDigitizer.h"
  18 #include "EMCAL/AliEMCALDigitizer.h"
  19 #include "EMCAL/AliEMCALClusterizerv1.h"
  20
  21 void PHOSHits2SDigits( Bool_t split=kFALSE, TString fileName = "galice.root") {
  22
  23   // usage :
  24   // 1. write SDigits in the same file as Hits --------------- (OK)
  25   //root [0] .L Reconstruct.C++
  26   //root [1] SDigits2Digits()
  27   // 2. write SDigits in a separate file, one per detector, from Hits --------------- (OK)
  28   //root [0] .L Reconstruct.C++
  29   //root [1] SDigits2Digits(kTRUE) // SDigits saved in [DET}.SDigits.root (DET=PHOS, EMCAL)
  30
  31   delete gAlice ;
  32   gAlice = 0 ;
  33
  34   AliPHOSSDigitizer * sdp = new AliPHOSSDigitizer(fileName) ;
  35   if (split)
  36     sdp->SetSplitFile() ;
  37   sdp->ExecuteTask("deb") ;
  38
  39   delete sdp ;
  40 }
  41
  42 //________________________________________________________________________
  43 void PHOSSDigits2Digits( Bool_t split=kFALSE, TString fileName = "galice.root") {
  44
  45  // usage :
  46   // 1. write SDigits in the same file as SDigits --------------- (OK)
  47   //root [0] .L Reconstruct.C++
  48   //root [1] SDigits2Digits()
  49   // 2. write SDigits in a separate file, one per detector, from SDigits --------------- (OK)
  50   //root [0] .L Reconstruct.C++
  51   //root [1] SDigitsDigits(kTRUE) // Digits saved in [DET}.Digits.root (DET=PHOS, EMCAL)
  52
  53   delete gAlice ;
  54   gAlice = 0 ;
  55
  56   // PHOS
  57   AliPHOSDigitizer * dp = 0 ;
  58
  59   if (split) {
  60     dp = new AliPHOSDigitizer("PHOS.SDigits.root") ;
  61     dp->SetSplitFile() ; }
  62   else
  63     dp = new AliPHOSDigitizer(fileName) ;
  64
  65   dp->ExecuteTask("deb") ;
  66
  67   delete dp ;
  68 }
  69
  70 //________________________________________________________________________
  71 void PHOSDigits2RecPoints( Bool_t split=kFALSE, TString fileName = "galice.root") {
  72
  73  // usage :
  74   // 1. write RecPoints in the same file as Digits --------------- OK
  75   //root [0] .L Reconstruct.C++
  76   //root [1] Digits2RecPoints()
  77   // 2. write RecPoints in a separate file, one per detector, from Digits --------------- OK
  78   //root [0] .L Reconstruct.C++
  79   //root [1] Digits2RecPoints(kTRUE) // RecPoints saved in [DET}.RecPoints.root (DET=PHOS, EMCAL)
  80
  81   delete gAlice ;
  82   gAlice = 0 ;
  83
  84   AliPHOSClusterizer * cp = 0 ;
  85
  86   if (split) {
  87     cp = new AliPHOSClusterizerv1("PHOS.Digits.root") ;
  88     cp->SetSplitFile() ; }
  89   else
  90     cp = new AliPHOSClusterizerv1(fileName) ;
  91
  92   cp->ExecuteTask("deb") ;
  93
  94   delete cp ;
  95 }
  96
  97 //________________________________________________________________________
  98 void PHOSRecPoints2TrackSegments( Bool_t split=kFALSE, TString fileName = "galice.root") {
  99
 100  // usage :
 101   // 1. write TrackSegments in the same file as RecPoints --------------- (OK)
 102   //root [0] .L Reconstruct.C++
 103   //root [1] RecPoints2TrackSegments()
 104   // 2. write TrackSegments in a separate file, one per detector, from RecPoints --------------- (Not needed)
 105   //root [0] .L Reconstruct.C++
 106   //root [1] RecPoints2TrackSegments(kTRUE) // TrackSegments saved in [DET}.RecData.root (DET=PHOS, EMCAL)
 107
 108   delete gAlice ;
 109   gAlice = 0 ;
 110
 111   AliPHOSTrackSegmentMaker * tmp = 0 ;
 112
 113   if (split)
 114     tmp = new AliPHOSTrackSegmentMakerv1("PHOS.RecData.root") ;
 115   else
 116     tmp = new AliPHOSTrackSegmentMakerv1(fileName) ;
 117
 118   tmp->ExecuteTask("deb") ;
 119
 120   delete tmp ;
 121 }
 122
 123 //________________________________________________________________________
 124 void PHOSTrackSegments2RecParticles( Bool_t split=kFALSE, TString fileName = "galice.root") {
 125
 126  // usage :
 127   // 1. write RecParticles in the same file as TrackSegments ---------------  (OK)
 128   //root [0] .L Reconstruct.C++
 129   //root [1] TrackSegments2RecParticles()
 130   // 2. write RecParticles in a separate file, one per detector, from TrackSegments --------------- (Not needed)
 131   //root [0] .L Reconstruct.C++
 132   //root [1] TrackSegments2RecParticles(kTRUE) // RecParticles saved in [DET}.RecData.root (DET=PHOS, EMCAL)
 133
 134   delete gAlice ;
 135   gAlice = 0 ;
 136
 137   AliPHOSPID * pp = 0 ;
 138
 139   if (split)
 140     pp = new AliPHOSPIDv1("PHOS.RecData.root") ;
 141   else
 142     pp = new AliPHOSPIDv1(fileName) ;
 143
 144   pp->ExecuteTask("deb") ;
 145
 146   delete pp ;
 147 }
 148
 149 //________________________________________________________________________
 150 void PHOSDigits2RecParticles( Bool_t split=kFALSE, TString fileName = "galice.root") {
 151
 152  // usage :
 153   // 1. write RecPoints, TrackSegments and RecParticles in the same file as Digits --------------- (OK)
 154   //root [0] .L Reconstruct.C++
 155   //root [1] Digits2RecParticles()
 156   // 2. write RecPoints , TrackSegments and RecParticles in a separate file, one per detector, from Digits --------------- (OK)
 157   //root [0] .L Reconstruct.C++
 158   //root [1] Digits2RecParticles(kTRUE) // TrackSegments saved in [DET}.RecData.root (DET=PHOS, EMCAL)
 159
 160
 161   delete gAlice ;
 162   gAlice = 0 ;
 163
 164   // PHOS
 165   AliPHOSClusterizer * cp = 0 ;
 166
 167   if (split) {
 168     cp = new AliPHOSClusterizerv1("PHOS.Digits.root") ;
 169     cp->SetSplitFile() ; }
 170   else
 171     cp = new AliPHOSClusterizerv1(fileName) ;
 172
 173   cp->ExecuteTask("deb") ;
 174
 175   if (split)
 176     delete cp ;
 177
 178   AliPHOSTrackSegmentMaker * tmp = 0 ;
 179
 180   if (split)
 181     tmp = new AliPHOSTrackSegmentMakerv1("PHOS.RecData.root") ;
 182   else
 183     tmp = new AliPHOSTrackSegmentMakerv1(fileName) ;
 184
 185   tmp->ExecuteTask("deb") ;
 186
 187   AliPHOSPID * pp = 0 ;
 188
 189   if (split)
 190     pp = new AliPHOSPIDv1("PHOS.RecData.root") ;
 191   else
 192     pp = new AliPHOSPIDv1(fileName) ;
 193
 194   pp->ExecuteTask("deb") ;
 195
 196   delete tmp;
 197   delete pp ;
 198 }
 199
 200 //________________________________________________________________________
 201 void PHOSHits2Digits (Bool_t split=kFALSE, TString fileName = "galice.root") {
 202   // usage :
 203   // 1. write (S)Digits in the same file as Hits --------------- (OK)
 204   //root [0] .L Reconstruct.C++
 205   //root [1] Hits2Digits()
 206   // 2. write (S)Digits in a separate file, one per detector, from Hits --------------- (OK)
 207   //root [0] .L Reconstruct.C++
 208   //root [1] Hits2Digits(kTRUE) // SDigits saved in [DET}.SDigits.root (DET=PHOS, EMCAL)
 209                                 // Digits  saved in [DET}.Digits.root  (DET=PHOS, EMCAL)
 210
 211   delete gAlice ;
 212   gAlice = 0 ;
 213
 214   //PHOS
 215   AliPHOSSDigitizer * sdp = new AliPHOSSDigitizer(fileName) ;
 216   if (split)
 217     sdp->SetSplitFile() ;
 218   sdp->ExecuteTask("deb") ;
 219
 220   if (split)
 221     delete sdp ;
 222
 223   AliPHOSDigitizer * dp = 0 ;
 224
 225   if (split) {
 226     dp = new AliPHOSDigitizer("PHOS.SDigits.root") ;
 227     dp->SetSplitFile() ; }
 228   else
 229     dp = new AliPHOSDigitizer(fileName) ;
 230
 231   dp->ExecuteTask("deb") ;
 232
 233   if (split)
 234     delete dp ;
 235
 236   if (!split) {
 237     delete sdp ;
 238     delete dp ;
 239   }
 240 }
 241
 242
 243