[u/mrichter/AliRoot.git] / EVE / EveDet / AliEveTOFSector.cxx

// $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.                                                 *
 **************************************************************************/
#include "AliEveTOFSector.h"

#include <AliTOFdigit.h>
#include <AliTOFGeometry.h>

#include <TEveTrans.h>
#include <TEveManager.h>
#include <TStyle.h>

Bool_t           AliEveTOFSector::fgStaticInitDone    = kFALSE;
TEveFrameBox*    AliEveTOFSector::fgTOFsectorFrameBox = 0;
TEveRGBAPalette* AliEveTOFSector::fgTOFsectorPalette  = 0;

//_______________________________________________________
ClassImp(AliEveTOFSector)

/* ************************************************************************ */

AliEveTOFSector::AliEveTOFSector(const Text_t* n, const Text_t* t) :
  TEveQuadSet(n, t),
  fTOFgeometry(new AliTOFGeometry()),
  fTOFarray(0x0),
  fTOFtree(0x0),
  fSector(-1),
  fDx(0), fDy(0), fDz(0),
  fAutoTrans (kTRUE),
  //fMinTime   (0), fMaxTime (0),
  fThreshold (5), fMaxVal (80),
  fSectorID  (0),
  //fPlateFlag0(kTRUE), fPlateFlag1(kTRUE), fPlateFlag2(kTRUE), fPlateFlag3(kTRUE), fPlateFlag4(kTRUE),
  fPlateFlag (0x0),
  //fFrameColor(4),
  //fRnrFrame  (kFALSE),
  fGeoManager(0)
{

  fPlateFlag = new Bool_t[5];
  for (Int_t ii=0; ii<5; ii++) fPlateFlag[ii]=kTRUE;


  //fGeoManager = (TGeoManager*)gEve->GetGeometry("$REVESYS/alice-data/alice_fullgeo.root");
  if (!fGeoManager) {
    printf("ERROR: no TGeo\n");
  }

}
/* ************************************************************************ */

AliEveTOFSector::AliEveTOFSector(TGeoManager *localGeoManager,
		     Int_t nSector)
  :
  TEveQuadSet(Form("Sector%i",nSector)),
  fTOFgeometry(new AliTOFGeometry()),
  fTOFarray(0x0),
  fTOFtree(0x0),
  fSector(nSector),
  fDx(0), fDy(0), fDz(0),
  fAutoTrans (kTRUE),
  //fMinTime   (0), fMaxTime (0),
  fThreshold (5), fMaxVal  (80),
  fSectorID  (nSector),
  //fPlateFlag0(kTRUE), fPlateFlag1(kTRUE), fPlateFlag2(kTRUE), fPlateFlag3(kTRUE), fPlateFlag4(kTRUE),
  fPlateFlag (0x0),
  //fFrameColor(4),
  //fRnrFrame  (kFALSE),
  fGeoManager(localGeoManager)
{

  fPlateFlag = new Bool_t[5];
  for (Int_t ii=0; ii<5; ii++) fPlateFlag[ii]=kTRUE;

  /*
  if (!fGeoManager) {
    printf("ERROR: no TGeo\n");
  }
  */

  InitModule();

}
/* ************************************************************************ */

AliEveTOFSector::AliEveTOFSector(TGeoManager *localGeoManager,
		     Int_t nSector,
		     TClonesArray *tofArray)
  :
  TEveQuadSet(Form("Sector%i",nSector)),
  fTOFgeometry(new AliTOFGeometry()),
  fTOFarray(tofArray),
  fTOFtree(0x0),
  fSector(nSector),
  fDx(0), fDy(0), fDz(0),
  fAutoTrans (kTRUE),
  //fMinTime   (0), fMaxTime (0),
  fThreshold (5), fMaxVal    (80),
  fSectorID  (nSector),
  //fPlateFlag0(kTRUE), fPlateFlag1(kTRUE), fPlateFlag2(kTRUE), fPlateFlag3(kTRUE), fPlateFlag4(kTRUE),
  fPlateFlag (0x0),
  //fFrameColor(4),
  //fRnrFrame  (kFALSE),
  fGeoManager(localGeoManager)
{

  fPlateFlag = new Bool_t[5];
  for (Int_t ii=0; ii<5; ii++) fPlateFlag[ii]=kTRUE;

  InitModule();

}
/* ************************************************************************ */

AliEveTOFSector::AliEveTOFSector(TGeoManager *localGeoManager,
		     Int_t nSector,
		     TTree *tofTree)
  :
  TEveQuadSet(Form("Sector%i",nSector)),
  fTOFgeometry(new AliTOFGeometry()),
  fTOFarray(0x0),
  fTOFtree(tofTree),
  fSector(nSector),
  fDx(0), fDy(0), fDz(0),
  fAutoTrans (kTRUE),
  //fMinTime   (0), fMaxTime (0),
  fThreshold (5), fMaxVal    (80),
  fSectorID  (nSector),
  //fPlateFlag0(kTRUE), fPlateFlag1(kTRUE), fPlateFlag2(kTRUE), fPlateFlag3(kTRUE), fPlateFlag4(kTRUE),
  fPlateFlag (0x0),
  //fFrameColor(4),
  //fRnrFrame  (kFALSE),
  fGeoManager(localGeoManager)
{

  fPlateFlag = new Bool_t[5];
  for (Int_t ii=0; ii<5; ii++) fPlateFlag[ii]=kTRUE;

  InitModule();

}
/* ************************************************************************ */

AliEveTOFSector::~AliEveTOFSector()
{
  /*
  fGeoManager = 0x0;
  delete fGeoManager;

  delete fTOFarray;
  fTOFarray = 0x0;
  */
  delete fPlateFlag;

}

/* ************************************************************************ */
/*
void AliEveTOFSector::SetDigitsInfo(AliEveTOFDigitsInfo* info)
{
  if(fInfo) fInfo->DecRefCount();
  fInfo = info;
  if(fInfo) fInfo->IncRefCount();

}
*/
/* ************************************************************************ */
void AliEveTOFSector::InitStatics()
{
  if (fgStaticInitDone) return;

  Float_t dx = 124.5;
  Float_t dz =  29.;
  Float_t dy = 370.6*2.;
  fgTOFsectorFrameBox = new TEveFrameBox();

  fgTOFsectorFrameBox->SetAABox(-dx*0.5, -dy*0.5, -dz*0.5, dx, dy, dz);
  fgTOFsectorFrameBox->SetFrameColor((Color_t) 32);//31);

  //fgTOFsectorPalette  = new TEveRGBAPalette(0, 2048); // TOT
  //fgTOFsectorPalette->SetLimits(0, 2048); 
  //fgTOFsectorPalette  = new TEveRGBAPalette(0, 8192); // TDC
  fgTOFsectorPalette  = new TEveRGBAPalette(0, 100000); // TDC
  fgTOFsectorPalette->SetOverflowAction(2);
  fgTOFsectorPalette->SetUnderflowAction(2);

  fgStaticInitDone = kTRUE;
}

/* ************************************************************************ */
void AliEveTOFSector::InitModule()
{

  fDx = fTOFgeometry->XPad()*fTOFgeometry->NpadX();
  //fDy = fTOFgeometry->XPad()*fTOFgeometry->NpadX();
  fDz = fTOFgeometry->ZPad()*fTOFgeometry->NpadZ();

  if (!fgStaticInitDone) InitStatics();

  SetFrame(fgTOFsectorFrameBox);
  SetPalette(fgTOFsectorPalette);
  //fFrame   = fgTOFsectorFrameBox;
  //fPalette = fgTOFsectorPalette;

  LoadQuads();
  ComputeBBox();
  SetTrans();

}

/* ************************************************************************ */
void AliEveTOFSector::LoadQuads()
{

  //Int_t n_col = gStyle->GetNumberOfColors();

  Int_t vol[5] = {fSectorID, -1, -1, -1, -1};
  Int_t informations[4] = {-1, -1, -1, -1};
  Int_t dummy[3] = {-1, -1, -1};
  Int_t tdc = -1;
  Int_t tot = -1;

  Double_t **coord = new Double_t*[4];
  for (Int_t ii=0; ii<4; ii++) coord[ii] = new Double_t[3];

  Reset(kQT_FreeQuad, kFALSE, 32);

  //printf(" fTOFarray->GetEntries() = %4i \n",fTOFarray->GetEntries());

  if (fTOFtree && !fTOFarray) {
    //printf("Hello world\n");
    TClonesArray* digitsTOFnew = new TClonesArray("AliTOFdigit",  300);

    fTOFarray = new TClonesArray("AliTOFdigit",  300);
    TClonesArray &ldigits = *fTOFarray;
    Int_t newCounter = 0;

    AliTOFdigit *digs;

    fTOFtree->SetBranchAddress("TOF",&digitsTOFnew);
    fTOFtree->GetEntry(0);
    for (Int_t digitNumber=0; digitNumber<digitsTOFnew->GetEntries(); digitNumber++) {

      //if (digitNumber==digitsTOF->GetEntries()-1) printf(" Hello  4 -> %3i digit of %i \n", digitNumber+1, digitsTOF->GetEntries());

      digs = (AliTOFdigit*)digitsTOFnew->UncheckedAt(digitNumber);

      if (digs->GetSector()!=fSectorID) continue;

      vol[1] = digs->GetPlate();  // Plate Number (0-4)
      vol[2] = digs->GetStrip();  // Strip Number (0-14/18)
      vol[3] = digs->GetPadx();   // Pad Number in x direction (0-47)
      vol[4] = digs->GetPadz();   // Pad Number in z direction (0-1)

      informations[0] = digs->GetTdc();
      informations[1] = digs->GetAdc();
      informations[2] = digs->GetToT();
      informations[3] = digs->GetTdcND();
      new (ldigits[newCounter++]) AliTOFdigit(dummy, vol, informations);
    }
  }

  AliTOFdigit *tofDigit;
  //printf("  0x%lx\n",fTOFarray);
  for (Int_t ii=0; ii<fTOFarray->GetEntries(); ii++) {

    tofDigit = (AliTOFdigit*)fTOFarray->UncheckedAt(ii);

    if (fPlateFlag[tofDigit->GetPlate()]) {

      vol[1] = tofDigit->GetPlate();
      vol[2] = tofDigit->GetStrip();
      vol[3] = tofDigit->GetPadz();
      vol[4] = tofDigit->GetPadx();

      tot = tofDigit->GetToT();
      tdc = tofDigit->GetTdc();

      //if (fSector==4 && fPlate==2 && fStrip==0) printf(" %2i  %1i\n", iPadX, iPadZ);
      //if (iPadX==23 || iPadX==24) printf(" %2i  %1i %2i \n", fSector, fPlate, fStrip);
      //if (vol[0]==4 && vol[1]==2 && vol[2]==7) {

      for (Int_t kk=0; kk<4; kk++) for (Int_t jj=0; jj<3; jj++) coord[kk][jj]=0.;

      fTOFgeometry->DetToSectorRF(vol, coord);
      /*
	printf("\n");
	printf("  %1i   %2i,   %f  %f  %f \n", vol[3], vol[4], coord[0][0], coord[0][1], coord[0][2]);
	printf("  %1i   %2i,   %f  %f  %f \n", vol[3], vol[4], coord[1][0], coord[1][1], coord[1][2]);
	printf("  %1i   %2i,   %f  %f  %f \n", vol[3], vol[4], coord[2][0], coord[2][1], coord[2][2]);
	printf("  %1i   %2i,   %f  %f  %f \n", vol[3], vol[4], coord[3][0], coord[3][1], coord[3][2]);
      */
      Float_t vertices[12]={(Float_t)coord[0][0], (Float_t)coord[0][1], (Float_t)coord[0][2],
			    (Float_t)coord[1][0], (Float_t)coord[1][1], (Float_t)coord[1][2],
			    (Float_t)coord[2][0], (Float_t)coord[2][1], (Float_t)coord[2][2],
			    (Float_t)coord[3][0], (Float_t)coord[3][1], (Float_t)coord[3][2]};

      AddQuad(vertices);
      //AddQuad((Float_t*)coord);
      //AddQuad(coord[0], coord[1], coord[2], 2.5, 3.5);
      //AddQuad(-2.5*0.5, -3.5*0.5, 2.5, 3.5);

      // In principle could have color based on number of neigbours. We
      // can insert the time-of-flight value for each pad

      //QuadValue((Int_t)tot);
      QuadValue((Int_t)tdc);
      QuadId(tofDigit);

    //}
    } // closed if control on plates switched on
  } // closed loop on TOF sector digits

  RefitPlex();

  fTOFarray = 0x0;

}

/* ************************************************************ */
void AliEveTOFSector::SetTrans()
{
  InitMainTrans();
  TEveTrans& t = RefMainTrans();

  //Int_t det[5] = {fSector, -1, -1, -1, -1};
  Char_t path[100];

  Int_t localSector = fSector;
  if (!fAutoTrans) localSector = 4;

  //fTOFgeometry->GetVolumePath(det,path);
  fTOFgeometry->GetVolumePath(localSector,path);
  fGeoManager->cd(path);
  TGeoHMatrix global = *fGeoManager->GetCurrentMatrix();
  Double_t *rotMat = global.GetRotationMatrix();
  Double_t *tr = global.GetTranslation();

  t.SetBaseVec(1, rotMat[0], rotMat[3], rotMat[6]);
  t.SetBaseVec(2, rotMat[1], rotMat[4], rotMat[7]);
  t.SetBaseVec(3, rotMat[2], rotMat[5], rotMat[8]);
  t.SetBaseVec(4, tr);
}

//-----------------------------------------------------

void AliEveTOFSector::SetSectorID(Int_t id)
{
  fSectorID = id;
  fSector   = id;
  if (fAutoTrans)
    SetTrans(); // Force repositioning.

  LoadQuads();
}

//-----------------------------------------------------

void AliEveTOFSector::SetPlate(Int_t nPlate, Bool_t r)
{

  fPlateFlag[nPlate] = r;

  //printf("   HELLO World ! %i %i %i \n", nPlate, r, fPlateFlag[nPlate]);
}

/******************************************************************************/

void AliEveTOFSector::SetThreshold(Short_t t)
{
  fThreshold = TMath::Min(t, (Short_t)(fMaxVal - 1));
  // ClearColorArray();
}

/******************************************************************************/

void AliEveTOFSector::SetMaxVal(Int_t mv)
{
  fMaxVal = TMath::Max(mv, (Int_t)(fThreshold + 1));
  //ClearColorArray();
}

/******************************************************************************/

void AliEveTOFSector::DigitSelected(Int_t idx)
{
  // Override control-click from TEveQuadSet

  DigitBase_t* qb   = GetDigit(idx);
  TObject* obj   = qb->fId.GetObject();
  AliTOFdigit* digs = dynamic_cast<AliTOFdigit*>(obj);
  // printf("AliEveTOFSector::QuadSelected "); Print();
  /*
  printf("  idx = %5i, value = %5d, obj = 0x%lx, digit = 0x%lx  ",
	 idx, qb->fValue, (ULong_t)obj, (ULong_t)digs);
  */
  if (digs)
    printf("\n Sector = %2i  Plate = %1i  Strip = %2i  PadZ = %1i PadX = %2i  ToT = %3i  Tof = %5i\n",
	   fSector , digs->GetPlate(), digs->GetStrip(), digs->GetPadz(), digs->GetPadx(), digs->GetToT(), digs->GetTdc());
  else printf("\n");

}

/******************************************************************************/
Commit	Line	Data
d810d0de	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 *
51346b82	7	* full copyright notice. *
d810d0de	8	**************************************************************************/
d810d0de	9	#include "AliEveTOFSector.h"
8616d353	10
8616d353	11	#include <AliTOFdigit.h>
28a21457	12	#include <AliTOFGeometry.h>
8616d353	13
a15e6d7d	14	#include <TEveTrans.h>
a15e6d7d	15	#include <TEveManager.h>
8616d353	16	#include <TStyle.h>
d810d0de	17
fd31e9de	18	Bool_t AliEveTOFSector::fgStaticInitDone = kFALSE;
d810d0de	19	TEveFrameBox* AliEveTOFSector::fgTOFsectorFrameBox = 0;
d810d0de	20	TEveRGBAPalette* AliEveTOFSector::fgTOFsectorPalette = 0;
8616d353	21
8616d353	22	//_______________________________________________________
d810d0de	23	ClassImp(AliEveTOFSector)
8616d353	24
	25	/* ************************************************************************ */
	26
d810d0de	27	AliEveTOFSector::AliEveTOFSector(const Text_t* n, const Text_t* t) :
84aff7a4	28	TEveQuadSet(n, t),
28a21457	29	fTOFgeometry(new AliTOFGeometry()),
8616d353	30	fTOFarray(0x0),
	31	fTOFtree(0x0),
	32	fSector(-1),
	33	fDx(0), fDy(0), fDz(0),
	34	fAutoTrans (kTRUE),
ee9d8eca	35	//fMinTime (0), fMaxTime (0),
ee9d8eca	36	fThreshold (5), fMaxVal (80),
8616d353	37	fSectorID (0),
ee9d8eca	38	//fPlateFlag0(kTRUE), fPlateFlag1(kTRUE), fPlateFlag2(kTRUE), fPlateFlag3(kTRUE), fPlateFlag4(kTRUE),
	39	fPlateFlag (0x0),
	40	//fFrameColor(4),
	41	//fRnrFrame (kFALSE),
8616d353	42	fGeoManager(0)
	43	{
	44
	45	fPlateFlag = new Bool_t[5];
	46	for (Int_t ii=0; ii<5; ii++) fPlateFlag[ii]=kTRUE;
	47
ee9d8eca	48
ee9d8eca	49	//fGeoManager = (TGeoManager*)gEve->GetGeometry("$REVESYS/alice-data/alice_fullgeo.root");
8616d353	50	if (!fGeoManager) {
	51	printf("ERROR: no TGeo\n");
	52	}
fd31e9de	53
ee9d8eca	54	}
8616d353	55	/* ************************************************************************ */
8616d353	56
ee9d8eca	57	AliEveTOFSector::AliEveTOFSector(TGeoManager *localGeoManager,
	58	Int_t nSector)
	59	:
	60	TEveQuadSet(Form("Sector%i",nSector)),
28a21457	61	fTOFgeometry(new AliTOFGeometry()),
8616d353	62	fTOFarray(0x0),
	63	fTOFtree(0x0),
	64	fSector(nSector),
	65	fDx(0), fDy(0), fDz(0),
	66	fAutoTrans (kTRUE),
ee9d8eca	67	//fMinTime (0), fMaxTime (0),
fd31e9de	68	fThreshold (5), fMaxVal (80),
8616d353	69	fSectorID (nSector),
ee9d8eca	70	//fPlateFlag0(kTRUE), fPlateFlag1(kTRUE), fPlateFlag2(kTRUE), fPlateFlag3(kTRUE), fPlateFlag4(kTRUE),
	71	fPlateFlag (0x0),
	72	//fFrameColor(4),
	73	//fRnrFrame (kFALSE),
8616d353	74	fGeoManager(localGeoManager)
	75	{
	76
	77	fPlateFlag = new Bool_t[5];
	78	for (Int_t ii=0; ii<5; ii++) fPlateFlag[ii]=kTRUE;
	79
	80	/*
	81	if (!fGeoManager) {
	82	printf("ERROR: no TGeo\n");
	83	}
	84	*/
	85
	86	InitModule();
fd31e9de	87
ee9d8eca	88	}
8616d353	89	/* ************************************************************************ */
8616d353	90
d810d0de	91	AliEveTOFSector::AliEveTOFSector(TGeoManager *localGeoManager,
8616d353	92	Int_t nSector,
	93	TClonesArray *tofArray)
	94	:
84aff7a4	95	TEveQuadSet(Form("Sector%i",nSector)),
28a21457	96	fTOFgeometry(new AliTOFGeometry()),
8616d353	97	fTOFarray(tofArray),
	98	fTOFtree(0x0),
	99	fSector(nSector),
	100	fDx(0), fDy(0), fDz(0),
	101	fAutoTrans (kTRUE),
ee9d8eca	102	//fMinTime (0), fMaxTime (0),
ee9d8eca	103	fThreshold (5), fMaxVal (80),
8616d353	104	fSectorID (nSector),
ee9d8eca	105	//fPlateFlag0(kTRUE), fPlateFlag1(kTRUE), fPlateFlag2(kTRUE), fPlateFlag3(kTRUE), fPlateFlag4(kTRUE),
	106	fPlateFlag (0x0),
	107	//fFrameColor(4),
	108	//fRnrFrame (kFALSE),
8616d353	109	fGeoManager(localGeoManager)
	110	{
	111
	112	fPlateFlag = new Bool_t[5];
	113	for (Int_t ii=0; ii<5; ii++) fPlateFlag[ii]=kTRUE;
	114
	115	InitModule();
	116
	117	}
	118	/* ************************************************************************ */
	119
ee9d8eca	120	AliEveTOFSector::AliEveTOFSector(TGeoManager *localGeoManager,
	121	Int_t nSector,
	122	TTree *tofTree)
	123	:
	124	TEveQuadSet(Form("Sector%i",nSector)),
28a21457	125	fTOFgeometry(new AliTOFGeometry()),
8616d353	126	fTOFarray(0x0),
	127	fTOFtree(tofTree),
	128	fSector(nSector),
	129	fDx(0), fDy(0), fDz(0),
	130	fAutoTrans (kTRUE),
ee9d8eca	131	//fMinTime (0), fMaxTime (0),
ee9d8eca	132	fThreshold (5), fMaxVal (80),
8616d353	133	fSectorID (nSector),
ee9d8eca	134	//fPlateFlag0(kTRUE), fPlateFlag1(kTRUE), fPlateFlag2(kTRUE), fPlateFlag3(kTRUE), fPlateFlag4(kTRUE),
	135	fPlateFlag (0x0),
	136	//fFrameColor(4),
	137	//fRnrFrame (kFALSE),
8616d353	138	fGeoManager(localGeoManager)
	139	{
	140
	141	fPlateFlag = new Bool_t[5];
	142	for (Int_t ii=0; ii<5; ii++) fPlateFlag[ii]=kTRUE;
	143
	144	InitModule();
	145
	146	}
	147	/* ************************************************************************ */
	148
d810d0de	149	AliEveTOFSector::~AliEveTOFSector()
8616d353	150	{
	151	/*
	152	fGeoManager = 0x0;
	153	delete fGeoManager;
	154
	155	delete fTOFarray;
	156	fTOFarray = 0x0;
	157	*/
	158	delete fPlateFlag;
	159
	160	}
	161
	162	/* ************************************************************************ */
	163	/*
d810d0de	164	void AliEveTOFSector::SetDigitsInfo(AliEveTOFDigitsInfo* info)
8616d353	165	{
	166	if(fInfo) fInfo->DecRefCount();
	167	fInfo = info;
	168	if(fInfo) fInfo->IncRefCount();
	169
	170	}
	171	*/
	172	/* ************************************************************************ */
d810d0de	173	void AliEveTOFSector::InitStatics()
8616d353	174	{
	175	if (fgStaticInitDone) return;
	176
	177	Float_t dx = 124.5;
	178	Float_t dz = 29.;
	179	Float_t dy = 370.6*2.;
84aff7a4	180	fgTOFsectorFrameBox = new TEveFrameBox();
8616d353	181
	182	fgTOFsectorFrameBox->SetAABox(-dx0.5, -dy0.5, -dz*0.5, dx, dy, dz);
	183	fgTOFsectorFrameBox->SetFrameColor((Color_t) 32);//31);
	184
84aff7a4	185	//fgTOFsectorPalette = new TEveRGBAPalette(0, 2048); // TOT
ee9d8eca	186	//fgTOFsectorPalette->SetLimits(0, 2048);
82dabae7	187	//fgTOFsectorPalette = new TEveRGBAPalette(0, 8192); // TDC
	188	fgTOFsectorPalette = new TEveRGBAPalette(0, 100000); // TDC
	189	fgTOFsectorPalette->SetOverflowAction(2);
	190	fgTOFsectorPalette->SetUnderflowAction(2);
8616d353	191
	192	fgStaticInitDone = kTRUE;
	193	}
	194
	195	/* ************************************************************************ */
d810d0de	196	void AliEveTOFSector::InitModule()
8616d353	197	{
	198
	199	fDx = fTOFgeometry->XPad()*fTOFgeometry->NpadX();
	200	//fDy = fTOFgeometry->XPad()*fTOFgeometry->NpadX();
	201	fDz = fTOFgeometry->ZPad()*fTOFgeometry->NpadZ();
	202
	203	if (!fgStaticInitDone) InitStatics();
	204
	205	SetFrame(fgTOFsectorFrameBox);
	206	SetPalette(fgTOFsectorPalette);
	207	//fFrame = fgTOFsectorFrameBox;
	208	//fPalette = fgTOFsectorPalette;
	209
51346b82	210	LoadQuads();
8616d353	211	ComputeBBox();
	212	SetTrans();
	213
	214	}
	215
	216	/* ************************************************************************ */
d810d0de	217	void AliEveTOFSector::LoadQuads()
8616d353	218	{
8616d353	219
8616d353	220	//Int_t n_col = gStyle->GetNumberOfColors();
	221
	222	Int_t vol[5] = {fSectorID, -1, -1, -1, -1};
	223	Int_t informations[4] = {-1, -1, -1, -1};
	224	Int_t dummy[3] = {-1, -1, -1};
	225	Int_t tdc = -1;
	226	Int_t tot = -1;
	227
	228	Double_t *coord = new Double_t[4];
	229	for (Int_t ii=0; ii<4; ii++) coord[ii] = new Double_t[3];
	230
ee9d8eca	231	Reset(kQT_FreeQuad, kFALSE, 32);
ee9d8eca	232
8616d353	233	//printf(" fTOFarray->GetEntries() = %4i \n",fTOFarray->GetEntries());
	234
	235	if (fTOFtree && !fTOFarray) {
	236	//printf("Hello world\n");
	237	TClonesArray* digitsTOFnew = new TClonesArray("AliTOFdigit", 300);
	238
	239	fTOFarray = new TClonesArray("AliTOFdigit", 300);
	240	TClonesArray &ldigits = *fTOFarray;
	241	Int_t newCounter = 0;
	242
	243	AliTOFdigit *digs;
	244
	245	fTOFtree->SetBranchAddress("TOF",&digitsTOFnew);
	246	fTOFtree->GetEntry(0);
	247	for (Int_t digitNumber=0; digitNumber<digitsTOFnew->GetEntries(); digitNumber++) {
	248
	249	//if (digitNumber==digitsTOF->GetEntries()-1) printf(" Hello 4 -> %3i digit of %i \n", digitNumber+1, digitsTOF->GetEntries());
51346b82	250
8616d353	251	digs = (AliTOFdigit*)digitsTOFnew->UncheckedAt(digitNumber);
	252
	253	if (digs->GetSector()!=fSectorID) continue;
	254
	255	vol[1] = digs->GetPlate(); // Plate Number (0-4)
	256	vol[2] = digs->GetStrip(); // Strip Number (0-14/18)
ee9d8eca	257	vol[3] = digs->GetPadx(); // Pad Number in x direction (0-47)
ee9d8eca	258	vol[4] = digs->GetPadz(); // Pad Number in z direction (0-1)
8616d353	259
	260	informations[0] = digs->GetTdc();
	261	informations[1] = digs->GetAdc();
	262	informations[2] = digs->GetToT();
	263	informations[3] = digs->GetTdcND();
	264	new (ldigits[newCounter++]) AliTOFdigit(dummy, vol, informations);
	265	}
	266	}
	267
	268	AliTOFdigit *tofDigit;
	269	//printf(" 0x%lx\n",fTOFarray);
	270	for (Int_t ii=0; ii<fTOFarray->GetEntries(); ii++) {
	271
	272	tofDigit = (AliTOFdigit*)fTOFarray->UncheckedAt(ii);
51346b82	273
8616d353	274	if (fPlateFlag[tofDigit->GetPlate()]) {
	275
	276	vol[1] = tofDigit->GetPlate();
	277	vol[2] = tofDigit->GetStrip();
	278	vol[3] = tofDigit->GetPadz();
	279	vol[4] = tofDigit->GetPadx();
	280
	281	tot = tofDigit->GetToT();
	282	tdc = tofDigit->GetTdc();
	283
	284	//if (fSector==4 && fPlate==2 && fStrip==0) printf(" %2i %1i\n", iPadX, iPadZ);
	285	//if (iPadX==23 \|\| iPadX==24) printf(" %2i %1i %2i \n", fSector, fPlate, fStrip);
	286	//if (vol[0]==4 && vol[1]==2 && vol[2]==7) {
	287
	288	for (Int_t kk=0; kk<4; kk++) for (Int_t jj=0; jj<3; jj++) coord[kk][jj]=0.;
	289
	290	fTOFgeometry->DetToSectorRF(vol, coord);
	291	/*
	292	printf("\n");
	293	printf(" %1i %2i, %f %f %f \n", vol[3], vol[4], coord[0][0], coord[0][1], coord[0][2]);
	294	printf(" %1i %2i, %f %f %f \n", vol[3], vol[4], coord[1][0], coord[1][1], coord[1][2]);
	295	printf(" %1i %2i, %f %f %f \n", vol[3], vol[4], coord[2][0], coord[2][1], coord[2][2]);
	296	printf(" %1i %2i, %f %f %f \n", vol[3], vol[4], coord[3][0], coord[3][1], coord[3][2]);
	297	*/
	298	Float_t vertices[12]={(Float_t)coord[0][0], (Float_t)coord[0][1], (Float_t)coord[0][2],
	299	(Float_t)coord[1][0], (Float_t)coord[1][1], (Float_t)coord[1][2],
	300	(Float_t)coord[2][0], (Float_t)coord[2][1], (Float_t)coord[2][2],
	301	(Float_t)coord[3][0], (Float_t)coord[3][1], (Float_t)coord[3][2]};
51346b82	302
8616d353	303	AddQuad(vertices);
	304	//AddQuad((Float_t*)coord);
	305	//AddQuad(coord[0], coord[1], coord[2], 2.5, 3.5);
	306	//AddQuad(-2.50.5, -3.50.5, 2.5, 3.5);
	307
	308	// In principle could have color based on number of neigbours. We
	309	// can insert the time-of-flight value for each pad
ee9d8eca	310
8616d353	311	//QuadValue((Int_t)tot);
	312	QuadValue((Int_t)tdc);
	313	QuadId(tofDigit);
51346b82	314
8616d353	315	//}
	316	} // closed if control on plates switched on
	317	} // closed loop on TOF sector digits
	318
	319	RefitPlex();
	320
	321	fTOFarray = 0x0;
51346b82	322
8616d353	323	}
	324
	325	/* ************************************************************ */
d810d0de	326	void AliEveTOFSector::SetTrans()
8616d353	327	{
a15e6d7d	328	InitMainTrans();
a15e6d7d	329	TEveTrans& t = RefMainTrans();
8616d353	330
	331	//Int_t det[5] = {fSector, -1, -1, -1, -1};
	332	Char_t path[100];
	333
	334	Int_t localSector = fSector;
	335	if (!fAutoTrans) localSector = 4;
	336
	337	//fTOFgeometry->GetVolumePath(det,path);
	338	fTOFgeometry->GetVolumePath(localSector,path);
	339	fGeoManager->cd(path);
	340	TGeoHMatrix global = *fGeoManager->GetCurrentMatrix();
	341	Double_t *rotMat = global.GetRotationMatrix();
	342	Double_t *tr = global.GetTranslation();
	343
a15e6d7d	344	t.SetBaseVec(1, rotMat[0], rotMat[3], rotMat[6]);
	345	t.SetBaseVec(2, rotMat[1], rotMat[4], rotMat[7]);
	346	t.SetBaseVec(3, rotMat[2], rotMat[5], rotMat[8]);
	347	t.SetBaseVec(4, tr);
8616d353	348	}
	349
	350	//-----------------------------------------------------
	351
d810d0de	352	void AliEveTOFSector::SetSectorID(Int_t id)
8616d353	353	{
	354	fSectorID = id;
	355	fSector = id;
	356	if (fAutoTrans)
	357	SetTrans(); // Force repositioning.
	358
	359	LoadQuads();
	360	}
	361
	362	//-----------------------------------------------------
	363
d810d0de	364	void AliEveTOFSector::SetPlate(Int_t nPlate, Bool_t r)
8616d353	365	{
	366
	367	fPlateFlag[nPlate] = r;
	368
	369	//printf(" HELLO World ! %i %i %i \n", nPlate, r, fPlateFlag[nPlate]);
	370	}
	371
57ffa5fb	372	/******************************************************************************/
8616d353	373
d810d0de	374	void AliEveTOFSector::SetThreshold(Short_t t)
8616d353	375	{
	376	fThreshold = TMath::Min(t, (Short_t)(fMaxVal - 1));
	377	// ClearColorArray();
	378	}
	379
57ffa5fb	380	/******************************************************************************/
8616d353	381
d810d0de	382	void AliEveTOFSector::SetMaxVal(Int_t mv)
8616d353	383	{
	384	fMaxVal = TMath::Max(mv, (Int_t)(fThreshold + 1));
	385	//ClearColorArray();
	386	}
	387
57ffa5fb	388	/******************************************************************************/
8616d353	389
d810d0de	390	void AliEveTOFSector::DigitSelected(Int_t idx)
8616d353	391	{
ee9d8eca	392	// Override control-click from TEveQuadSet
8616d353	393
84aff7a4	394	DigitBase_t* qb = GetDigit(idx);
8616d353	395	TObject* obj = qb->fId.GetObject();
8616d353	396	AliTOFdigit* digs = dynamic_cast<AliTOFdigit*>(obj);
d810d0de	397	// printf("AliEveTOFSector::QuadSelected "); Print();
ee9d8eca	398	/*
8616d353	399	printf(" idx = %5i, value = %5d, obj = 0x%lx, digit = 0x%lx ",
8616d353	400	idx, qb->fValue, (ULong_t)obj, (ULong_t)digs);
ee9d8eca	401	*/
8616d353	402	if (digs)
ee9d8eca	403	printf("\n Sector = %2i Plate = %1i Strip = %2i PadZ = %1i PadX = %2i ToT = %3i Tof = %5i\n",
ee9d8eca	404	fSector , digs->GetPlate(), digs->GetStrip(), digs->GetPadz(), digs->GetPadx(), digs->GetToT(), digs->GetTdc());
8616d353	405	else printf("\n");
	406
	407	}
	408
57ffa5fb	409	/******************************************************************************/