Move contents of EVE/Alieve to EVE/EveDet as most code will remain there.

[u/mrichter/AliRoot.git] / EVE / test-macros / tpc_gui_hack.C

// $Id$
// Main authors: Matevz Tadel & Alja Mrak-Tadel: 2006, 2007

/**************************************************************************
 * Copyright(c) 1998-2008, ALICE Experiment at CERN, all rights reserved. *
 * See http://aliceinfo.cern.ch/Offline/AliRoot/License.html for          *
 * full copyright notice.                                                 *
 **************************************************************************/

// Function to spawn a gui for reading rootified raw-data from TPC sector test.

#ifdef __CINT__

class AliRawReaderRoot;

class AliEveTPCData;
class AliEveTPCSector2D;
class AliEveTPCSector3D;

#else

#include <TEve.h>
#include <TEveManager.h>
#include <EveDet/AliEveTPCData.h>
#include <EveDet/AliEveTPCSector2D.h>
#include <EveDet/AliEveTPCSector3D.h>

#include <RAW/AliRawReaderRoot.h>
#include <TPC/AliTPCRawStream.h>

#include <TSystem.h>
#include <TStyle.h>

#endif


AliEveTPCSectorData* su = 0;
AliEveTPCSectorData* sl = 0;

AliEveTPCLoader* loader = 0;

void tpc_gui_hack(const char *file=0, Int_t ievent=0)
{
  gStyle->SetPalette(1, 0);

  AliEveTPCLoader* l = new AliEveTPCLoader; loader = l;
  AliEveTPCData*   d = new AliEveTPCData;
  // d->SetLoadPedestal(5);
  d->SetLoadThreshold(5);
  d->SetAutoPedestal(kTRUE);
  l->SetData(d);
  l->SetDoubleSR(kTRUE);
  l->SetInitParams(50, 980, 15);
  // l->SetTPCEquipementMap("EquipmentIdMap.data");

  su = d->GetSectorData( 4, kTRUE);
  sl = d->GetSectorData(13, kTRUE);

  gEve->AddElement(l);
  gEve->NotifyBrowser(l);
  gEve->DrawElement(l);

  if(file != 0) {
    l->SetFile(file);
    l->OpenFile();
    l->GotoEvent(ievent);
  }
}

void disable_pad(Int_t row, Int_t pad, AliEveTPCSectorData* sd,
                 Int_t thrExt=10, Float_t thrFac=2)
{
  if(row < 0 || row >= AliEveTPCSectorData::GetNAllRows())
    { printf("row off, %d\n", row); return; }

  Int_t npads = AliEveTPCSectorData::GetNPadsInRow(row);
  if(pad < 0) pad = npads + pad;
  if(pad < 0 || pad >= npads) { printf("pad off\n"); return; }

  sd->AddPadRowHack(row, pad, thrExt, thrFac);
}

void disable_std()
{
  // The noisy pads in lower inner seg, left edge
  disable_pad(30, 1, sl, 200, 5);
  disable_pad(30, 2, sl, 200, 5);
  disable_pad(31, 1, sl, 200, 5);

  // The noisy pads in lower inner seg, middle
  disable_pad(31, 44, sl);
  disable_pad(31, 43, sl);
  disable_pad(30, 43, sl);

  for(Int_t r=16; r<32; ++r) {
    disable_pad(r, 0, su, 20, 3);
    disable_pad(r, 1, su, 20, 3);
    disable_pad(r, 2, su, 20, 3);
    disable_pad(r, 3, su, 20, 3);
    disable_pad(r, 4, su, 20, 3);
  }

  { // Top 12, 4 pads on the negative side.
    AliEveTPCSectorData::SegmentInfo& o1si = AliEveTPCSectorData::GetOut1Seg();
    Int_t last = o1si.GetLastRow();
    for(Int_t r=last - 11; r<=last; ++r) {
      disable_pad(r, -1, su, 20, 4);
      disable_pad(r, -2, su, 10, 3);
      disable_pad(r, -3, su, 10, 3);
      disable_pad(r, -4, su,  5, 2.5);
      disable_pad(r, -5, su,  5, 2.5);
      disable_pad(r,  0, su,  5, 2.5);
      disable_pad(r,  1, su, 20, 4);
      disable_pad(r,  2, su, 10, 3);
      disable_pad(r,  3, su, 10, 3);
      disable_pad(r,  4, su,  5, 2.5);

      disable_pad(r, -1, sl, 20, 4);
      disable_pad(r, -2, sl, 10, 3);
      disable_pad(r, -3, sl, 10, 3);
      disable_pad(r, -4, sl,  5, 2.5);
      disable_pad(r, -5, sl,  5, 2.5);
    }

    disable_pad(last-8, -4, sl, 5, 2.5);

    disable_pad(last-9, -2, sl, 5, 2);
    disable_pad(last-9, -3, sl, 5, 2.5);
    disable_pad(last-9, -4, sl, 5, 2.5);
    disable_pad(last-9, -5, sl, 5, 2.5);
    disable_pad(last-8, -5, sl, 5, 2.5);
  }

  {
    AliEveTPCSectorData::SegmentInfo& o2si = AliEveTPCSectorData::GetOut2Seg();
    Int_t first = o2si.GetFirstRow();
    Int_t last  = o2si.GetLastRow();

    for(Int_t r=first; r<=last; ++r) {
      disable_pad(r, -1, su, 30, 3.5);
      disable_pad(r, -2, su, 25, 3);
      disable_pad(r, -3, su, 20, 3);
      disable_pad(r,  0, su, 30, 3.5);
      disable_pad(r,  1, su, 25, 3);

      disable_pad(r,  0, sl, 30, 3.5);
      disable_pad(r,  1, sl, 25, 3);
      disable_pad(r,  2, sl, 20, 3);
      disable_pad(r,  3, sl, 20, 3);
      disable_pad(r,  4, sl, 20, 3);
      disable_pad(r,  5, sl, 20, 3);
      disable_pad(r, -1, sl, 30, 3.5);
      disable_pad(r, -2, sl, 25, 3);
      disable_pad(r, -3, sl, 25, 3);
      disable_pad(r, -4, sl, 25, 3);
      disable_pad(r, -5, sl, 25, 3);
    }

    for(Int_t pad=3; pad<30; ++pad) {
      disable_pad(last, -pad-1, sl, 20, 3);
      disable_pad(last,    pad, su, 20, 3);
    }
    for(Int_t pad=3; pad<15; ++pad) {
      disable_pad(last-1, -pad-1, sl, 20, 3);
      disable_pad(last-2, -pad-1, sl, 20, 3);
    }
    disable_pad(last, 2, su, 20, 3);
  }

  loader->LoadEvent(); loader->UpdateSectors();
}