PHOS/AliPHOSReconstruct.C

   1 // YS Subatech Mai 2002
   2 // YK Subatech 6 Aug 2002
   3
   4 // 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 Hits2SDigits( 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  //  AliEMCALSDigitizer * sde = new AliEMCALSDigitizer(fileName) ;
  42 //   if (split)
  43 //     sde->SetSplitFile() ;
  44 //   sde->ExecuteTask("deb") ;
  45
  46 //   delete sde ;
  47
  48 }
  49
  50 //________________________________________________________________________
  51 void SDigits2Digits( Bool_t split=kFALSE, TString fileName = "galice.root") {
  52
  53  // usage :
  54   // 1. write SDigits in the same file as SDigits --------------- (OK)
  55   //root [0] .L Reconstruct.C++
  56   //root [1] SDigits2Digits()
  57   // 2. write SDigits in a separate file, one per detector, from SDigits --------------- (OK)
  58   //root [0] .L Reconstruct.C++
  59   //root [1] SDigitsDigits(kTRUE) // Digits saved in [DET}.Digits.root (DET=PHOS, EMCAL)
  60
  61   delete gAlice ;
  62   gAlice = 0 ;
  63
  64   // PHOS
  65   AliPHOSDigitizer * dp = 0 ;
  66
  67   if (split) {
  68     dp = new AliPHOSDigitizer("PHOS.SDigits.root") ;
  69     dp->SetSplitFile() ; }
  70   else
  71     dp = new AliPHOSDigitizer(fileName) ;
  72
  73   dp->ExecuteTask("deb") ;
  74
  75   delete dp ;
  76
  77   //EMCAL
  78 //   AliEMCALDigitizer * de = 0 ;
  79
  80 //   if (split) {
  81 //     de = new AliEMCALDigitizer("EMCAL.SDigits.root") ;
  82 //     de->SetSplitFile() ;
  83 //   } else
  84 //     de = new AliEMCALDigitizer(fileName) ;
  85
  86 //   de->ExecuteTask("deb") ;
  87
  88 //   delete de ;
  89 }
  90
  91 //________________________________________________________________________
  92 void Digits2RecPoints( Bool_t split=kFALSE, TString fileName = "galice.root") {
  93
  94  // usage :
  95   // 1. write RecPoints in the same file as Digits --------------- OK
  96   //root [0] .L Reconstruct.C++
  97   //root [1] Digits2RecPoints()
  98   // 2. write RecPoints in a separate file, one per detector, from Digits --------------- OK
  99   //root [0] .L Reconstruct.C++
 100   //root [1] Digits2RecPoints(kTRUE) // RecPoints saved in [DET}.RecPoints.root (DET=PHOS, EMCAL)
 101
 102   delete gAlice ;
 103   gAlice = 0 ;
 104
 105 // PHOS
 106   AliPHOSClusterizer * cp = 0 ;
 107
 108   if (split) {
 109     cp = new AliPHOSClusterizerv1("PHOS.Digits.root") ;
 110     cp->SetSplitFile() ; }
 111   else
 112     cp = new AliPHOSClusterizerv1(fileName) ;
 113
 114   cp->ExecuteTask("deb") ;
 115
 116   delete cp ;
 117
 118 //   //EMCAL
 119 //   AliEMCALClusterizerv1 * ce = 0 ;
 120
 121 //   if (split) {
 122 //     ce = new AliEMCALClusterizerv1("EMCAL.Digits.root") ;
 123 //     ce->SetSplitFile() ;
 124 //   } else
 125 //     ce = new AliEMCALClusterizerv1(fileName) ;
 126
 127 //   ce->ExecuteTask("deb") ;
 128
 129 //   delete ce ;
 130 }
 131
 132 //________________________________________________________________________
 133 void RecPoints2TrackSegments( Bool_t split=kFALSE, TString fileName = "galice.root") {
 134
 135  // usage :
 136   // 1. write TrackSegments in the same file as RecPoints --------------- (OK)
 137   //root [0] .L Reconstruct.C++
 138   //root [1] RecPoints2TrackSegments()
 139   // 2. write TrackSegments in a separate file, one per detector, from RecPoints --------------- (Not needed)
 140   //root [0] .L Reconstruct.C++
 141   //root [1] RecPoints2TrackSegments(kTRUE) // TrackSegments saved in [DET}.RecData.root (DET=PHOS, EMCAL)
 142
 143   delete gAlice ;
 144   gAlice = 0 ;
 145
 146   AliPHOSTrackSegmentMaker * tmp = 0 ;
 147
 148   if (split)
 149     tmp = new AliPHOSTrackSegmentMakerv1("PHOS.RecData.root") ;
 150   else
 151     tmp = new AliPHOSTrackSegmentMakerv1(fileName) ;
 152
 153   tmp->ExecuteTask("deb") ;
 154
 155   delete tmp ;
 156 }
 157
 158 //________________________________________________________________________
 159 void TrackSegments2RecParticles( Bool_t split=kFALSE, TString fileName = "galice.root") {
 160
 161  // usage :
 162   // 1. write RecParticles in the same file as TrackSegments ---------------  (OK)
 163   //root [0] .L Reconstruct.C++
 164   //root [1] TrackSegments2RecParticles()
 165   // 2. write RecParticles in a separate file, one per detector, from TrackSegments --------------- (Not needed)
 166   //root [0] .L Reconstruct.C++
 167   //root [1] TrackSegments2RecParticles(kTRUE) // RecParticles saved in [DET}.RecData.root (DET=PHOS, EMCAL)
 168
 169   delete gAlice ;
 170   gAlice = 0 ;
 171
 172   AliPHOSPID * pp = 0 ;
 173
 174   if (split)
 175     pp = new AliPHOSPIDv1("PHOS.RecData.root") ;
 176   else
 177     pp = new AliPHOSPIDv1(fileName) ;
 178
 179   pp->ExecuteTask("deb") ;
 180
 181   delete pp ;
 182 }
 183
 184 //________________________________________________________________________
 185 void Digits2RecParticles( Bool_t split=kFALSE, TString fileName = "galice.root") {
 186
 187  // usage :
 188   // 1. write RecPoints, TrackSegments and RecParticles in the same file as Digits --------------- (OK)
 189   //root [0] .L Reconstruct.C++
 190   //root [1] Digits2RecParticles()
 191   // 2. write RecPoints , TrackSegments and RecParticles in a separate file, one per detector, from Digits --------------- (OK)
 192   //root [0] .L Reconstruct.C++
 193   //root [1] Digits2RecParticles(kTRUE) // TrackSegments saved in [DET}.RecData.root (DET=PHOS, EMCAL)
 194
 195
 196   delete gAlice ;
 197   gAlice = 0 ;
 198
 199   // PHOS
 200   AliPHOSClusterizer * cp = 0 ;
 201
 202   if (split) {
 203     cp = new AliPHOSClusterizerv1("PHOS.Digits.root") ;
 204     cp->SetSplitFile() ; }
 205   else
 206     cp = new AliPHOSClusterizerv1(fileName) ;
 207
 208   cp->ExecuteTask("deb") ;
 209
 210   if (split)
 211     delete cp ;
 212
 213   AliPHOSTrackSegmentMaker * tmp = 0 ;
 214
 215   if (split)
 216     tmp = new AliPHOSTrackSegmentMakerv1("PHOS.RecData.root") ;
 217   else
 218     tmp = new AliPHOSTrackSegmentMakerv1(fileName) ;
 219
 220   tmp->ExecuteTask("deb") ;
 221
 222   AliPHOSPID * pp = 0 ;
 223
 224   if (split)
 225     pp = new AliPHOSPIDv1("PHOS.RecData.root") ;
 226   else
 227     pp = new AliPHOSPIDv1(fileName) ;
 228
 229   pp->ExecuteTask("deb") ;
 230
 231   delete tmp;
 232   delete pp ;
 233 }
 234
 235 //________________________________________________________________________
 236 void Hits2Digits (Bool_t split=kFALSE, TString fileName = "galice.root") {
 237   // usage :
 238   // 1. write (S)Digits in the same file as Hits --------------- (OK)
 239   //root [0] .L Reconstruct.C++
 240   //root [1] Hits2Digits()
 241   // 2. write (S)Digits in a separate file, one per detector, from Hits --------------- (OK)
 242   //root [0] .L Reconstruct.C++
 243   //root [1] Hits2Digits(kTRUE) // SDigits saved in [DET}.SDigits.root (DET=PHOS, EMCAL)
 244                                 // Digits  saved in [DET}.Digits.root  (DET=PHOS, EMCAL)
 245
 246   delete gAlice ;
 247   gAlice = 0 ;
 248
 249   //PHOS
 250   AliPHOSSDigitizer * sdp = new AliPHOSSDigitizer(fileName) ;
 251   if (split)
 252     sdp->SetSplitFile() ;
 253   sdp->ExecuteTask("deb") ;
 254
 255   if (split)
 256     delete sdp ;
 257
 258   AliPHOSDigitizer * dp = 0 ;
 259
 260   if (split) {
 261     dp = new AliPHOSDigitizer("PHOS.SDigits.root") ;
 262     dp->SetSplitFile() ; }
 263   else
 264     dp = new AliPHOSDigitizer(fileName) ;
 265
 266   dp->ExecuteTask("deb") ;
 267
 268   if (split)
 269     delete dp ;
 270
 271   if (!split) {
 272     delete sdp ;
 273     delete dp ;
 274   }
 275
 276 //   //EMCAL
 277 //   AliEMCALSDigitizer * sde = new AliEMCALSDigitizer(fileName) ;
 278 //   if (split)
 279 //     sde->SetSplitFile() ;
 280 //   sde->ExecuteTask("deb") ;
 281
 282 //   delete sde ;
 283
 284 //   AliEMCALDigitizer * de = 0 ;
 285 //   if (split) {
 286 //     de = new AliEMCALDigitizer("EMCAL.SDigits.root") ;
 287 //     de->SetSplitFile() ;
 288 //   } else
 289 //     de = new AliEMCALDigitizer(fileName) ;
 290
 291 //   de->ExecuteTask("deb") ;
 292
 293 //   delete de ;
 294
 295 }
 296
 297
 298