EVE/test-macros/tpc_gui_hack.C

   1 // $Id$
   2 // Main authors: Matevz Tadel & Alja Mrak-Tadel: 2006, 2007
   3
   4 /**************************************************************************
   5  * Copyright(c) 1998-2008, ALICE Experiment at CERN, all rights reserved. *
   6  * See http://aliceinfo.cern.ch/Offline/AliRoot/License.html for          *
   7  * full copyright notice.                                                 *
   8  **************************************************************************/
   9
  10 // Function to spawn a gui for reading rootified raw-data from TPC sector test.
  11
  12 #ifdef __CINT__
  13
  14 class AliRawReaderRoot;
  15
  16 class AliEveTPCData;
  17 class AliEveTPCSector2D;
  18 class AliEveTPCSector3D;
  19
  20 #else
  21
  22 #include <TEve.h>
  23 #include <TEveManager.h>
  24 #include <EveDet/AliEveTPCData.h>
  25 #include <EveDet/AliEveTPCSector2D.h>
  26 #include <EveDet/AliEveTPCSector3D.h>
  27
  28 #include <RAW/AliRawReaderRoot.h>
  29 #include <TPC/AliTPCRawStream.h>
  30
  31 #include <TSystem.h>
  32 #include <TStyle.h>
  33
  34 #endif
  35
  36
  37 AliEveTPCSectorData* su = 0;
  38 AliEveTPCSectorData* sl = 0;
  39
  40 AliEveTPCLoader* loader = 0;
  41
  42 void tpc_gui_hack(const char *file=0, Int_t ievent=0)
  43 {
  44   gStyle->SetPalette(1, 0);
  45
  46   AliEveTPCLoader* l = new AliEveTPCLoader; loader = l;
  47   AliEveTPCData*   d = new AliEveTPCData;
  48   // d->SetLoadPedestal(5);
  49   d->SetLoadThreshold(5);
  50   d->SetAutoPedestal(kTRUE);
  51   l->SetData(d);
  52   l->SetDoubleSR(kTRUE);
  53   l->SetInitParams(50, 980, 15);
  54   // l->SetTPCEquipementMap("EquipmentIdMap.data");
  55
  56   su = d->GetSectorData( 4, kTRUE);
  57   sl = d->GetSectorData(13, kTRUE);
  58
  59   gEve->AddElement(l);
  60   gEve->NotifyBrowser(l);
  61   gEve->DrawElement(l);
  62
  63   if(file != 0) {
  64     l->SetFile(file);
  65     l->OpenFile();
  66     l->GotoEvent(ievent);
  67   }
  68 }
  69
  70 void disable_pad(Int_t row, Int_t pad, AliEveTPCSectorData* sd,
  71                  Int_t thrExt=10, Float_t thrFac=2)
  72 {
  73   if(row < 0 || row >= AliEveTPCSectorData::GetNAllRows())
  74     { printf("row off, %d\n", row); return; }
  75
  76   Int_t npads = AliEveTPCSectorData::GetNPadsInRow(row);
  77   if(pad < 0) pad = npads + pad;
  78   if(pad < 0 || pad >= npads) { printf("pad off\n"); return; }
  79
  80   sd->AddPadRowHack(row, pad, thrExt, thrFac);
  81 }
  82
  83 void disable_std()
  84 {
  85   // The noisy pads in lower inner seg, left edge
  86   disable_pad(30, 1, sl, 200, 5);
  87   disable_pad(30, 2, sl, 200, 5);
  88   disable_pad(31, 1, sl, 200, 5);
  89
  90   // The noisy pads in lower inner seg, middle
  91   disable_pad(31, 44, sl);
  92   disable_pad(31, 43, sl);
  93   disable_pad(30, 43, sl);
  94
  95   for(Int_t r=16; r<32; ++r) {
  96     disable_pad(r, 0, su, 20, 3);
  97     disable_pad(r, 1, su, 20, 3);
  98     disable_pad(r, 2, su, 20, 3);
  99     disable_pad(r, 3, su, 20, 3);
 100     disable_pad(r, 4, su, 20, 3);
 101   }
 102
 103   { // Top 12, 4 pads on the negative side.
 104     AliEveTPCSectorData::SegmentInfo& o1si = AliEveTPCSectorData::GetOut1Seg();
 105     Int_t last = o1si.GetLastRow();
 106     for(Int_t r=last - 11; r<=last; ++r) {
 107       disable_pad(r, -1, su, 20, 4);
 108       disable_pad(r, -2, su, 10, 3);
 109       disable_pad(r, -3, su, 10, 3);
 110       disable_pad(r, -4, su,  5, 2.5);
 111       disable_pad(r, -5, su,  5, 2.5);
 112       disable_pad(r,  0, su,  5, 2.5);
 113       disable_pad(r,  1, su, 20, 4);
 114       disable_pad(r,  2, su, 10, 3);
 115       disable_pad(r,  3, su, 10, 3);
 116       disable_pad(r,  4, su,  5, 2.5);
 117
 118       disable_pad(r, -1, sl, 20, 4);
 119       disable_pad(r, -2, sl, 10, 3);
 120       disable_pad(r, -3, sl, 10, 3);
 121       disable_pad(r, -4, sl,  5, 2.5);
 122       disable_pad(r, -5, sl,  5, 2.5);
 123     }
 124
 125     disable_pad(last-8, -4, sl, 5, 2.5);
 126
 127     disable_pad(last-9, -2, sl, 5, 2);
 128     disable_pad(last-9, -3, sl, 5, 2.5);
 129     disable_pad(last-9, -4, sl, 5, 2.5);
 130     disable_pad(last-9, -5, sl, 5, 2.5);
 131     disable_pad(last-8, -5, sl, 5, 2.5);
 132   }
 133
 134   {
 135     AliEveTPCSectorData::SegmentInfo& o2si = AliEveTPCSectorData::GetOut2Seg();
 136     Int_t first = o2si.GetFirstRow();
 137     Int_t last  = o2si.GetLastRow();
 138
 139     for(Int_t r=first; r<=last; ++r) {
 140       disable_pad(r, -1, su, 30, 3.5);
 141       disable_pad(r, -2, su, 25, 3);
 142       disable_pad(r, -3, su, 20, 3);
 143       disable_pad(r,  0, su, 30, 3.5);
 144       disable_pad(r,  1, su, 25, 3);
 145
 146       disable_pad(r,  0, sl, 30, 3.5);
 147       disable_pad(r,  1, sl, 25, 3);
 148       disable_pad(r,  2, sl, 20, 3);
 149       disable_pad(r,  3, sl, 20, 3);
 150       disable_pad(r,  4, sl, 20, 3);
 151       disable_pad(r,  5, sl, 20, 3);
 152       disable_pad(r, -1, sl, 30, 3.5);
 153       disable_pad(r, -2, sl, 25, 3);
 154       disable_pad(r, -3, sl, 25, 3);
 155       disable_pad(r, -4, sl, 25, 3);
 156       disable_pad(r, -5, sl, 25, 3);
 157     }
 158
 159     for(Int_t pad=3; pad<30; ++pad) {
 160       disable_pad(last, -pad-1, sl, 20, 3);
 161       disable_pad(last,    pad, su, 20, 3);
 162     }
 163     for(Int_t pad=3; pad<15; ++pad) {
 164       disable_pad(last-1, -pad-1, sl, 20, 3);
 165       disable_pad(last-2, -pad-1, sl, 20, 3);
 166     }
 167     disable_pad(last, 2, su, 20, 3);
 168   }
 169
 170   loader->LoadEvent(); loader->UpdateSectors();
 171 }