EMCAL/AliEMCALReconstruct.C

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