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

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

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

//
// Class to visualize the TOF digit information
// in TOF sector frame
//
// Author: A. De Caro (email: decaro@sa.infn.t)
//

#include <TGeoManager.h>
#include <TClonesArray.h>
#include <TTree.h>

#include <TEveTrans.h>

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

#include "AliEveTOFSector.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)
{

  // default ctr

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


  //fGeoManager = AliEveEventManager::AssertGeometry();
  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)
{

  // ctr

  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)
{

  // ctr

  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)
{

  // ctr

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

  InitModule();

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

AliEveTOFSector::~AliEveTOFSector()
{

  // dtr

  /*
  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()
{
  //
  // To initialize statistic variables
  //

  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));
  fgTOFsectorFrameBox->IncRefCount();

  //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);
  fgTOFsectorPalette->IncRefCount();

  fgStaticInitDone = kTRUE;
}

/* ************************************************************************ */
void AliEveTOFSector::InitModule()
{
  //
  // To initialize TOF sector frame variables
  //

  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()
{
  //
  // Load TOF digit as illuminated pad
  //

  //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[4][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()
{
  //
  // Set the translation matrix for TOF sector
  //

  InitMainTrans();
  TEveTrans& t = RefMainTrans();

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

  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)
{
  //
  // Set TOF sector ID to be visualized
  //

  fSectorID = id;
  fSector   = id;
  if (fAutoTrans)
    SetTrans(); // Force repositioning.

  LoadQuads();
}

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

void AliEveTOFSector::SetPlate(Int_t nPlate, Bool_t r)
{
  //
  // Set visualization flag to visualize TOF modules
  //

  fPlateFlag[nPlate] = r;

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

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

void AliEveTOFSector::SetThreshold(Short_t t)
{
  //
  // Set visualization threshold
  //

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

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

void AliEveTOFSector::SetMaxVal(Int_t mv)
{
  //
  // Set visualization max value
  //

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

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

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

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