[u/mrichter/AliRoot.git] / EVE / EveDet / AliEveITSModule.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 "AliEveITSModule.h"

#include <AliITSgeomTGeo.h>
#include <AliITSsegmentationSPD.h>
#include <AliITSsegmentationSDD.h>
#include <AliITSsegmentationSSD.h>

#include <AliITSdigitSPD.h>
#include <AliITSdigitSDD.h>
#include <AliITSdigitSSD.h>

#include <TStyle.h>


//______________________________________________________________________________
//
// Visualization of an ITS module.

ClassImp(AliEveITSModule)

Bool_t AliEveITSModule::fgStaticInitDone = kFALSE;

TEveFrameBox*    AliEveITSModule::fgSPDFrameBox = 0;
TEveFrameBox*    AliEveITSModule::fgSDDFrameBox = 0;
TEveFrameBox*    AliEveITSModule::fgSSDFrameBox = 0;

TEveRGBAPalette* AliEveITSModule::fgSPDPalette  = 0;
TEveRGBAPalette* AliEveITSModule::fgSDDPalette  = 0;
TEveRGBAPalette* AliEveITSModule::fgSSDPalette  = 0;

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

AliEveITSModule::AliEveITSModule(const Text_t* n, const Text_t* t) :
  TEveQuadSet(n, t),
  fInfo(0),
  fID(-1), fDetID(-1),
  fLayer(-1), fLadder(-1), fDet(-1),
  fDx(0), fDz(0), fDy(0)
{
  // Constructor.
}

AliEveITSModule::AliEveITSModule(Int_t gid, AliEveITSDigitsInfo* info) :
  TEveQuadSet(Form("ITS module %d", gid)),
  fInfo  (0),
  fID(-1), fDetID(-1),
  fLayer(-1), fLadder(-1), fDet(-1),
  fDx(0), fDz(0), fDy(0)
{
  // Constructor with module id and data-source.

  SetDigitsInfo(info);
  SetID(gid);
}

AliEveITSModule::~AliEveITSModule()
{
  // Destructor.

  if(fInfo) fInfo->DecRefCount();
}

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

void AliEveITSModule::InitStatics(AliEveITSDigitsInfo* info)
{
  // Initialize static variables,

  if (fgStaticInitDone) return;
  fgStaticInitDone = kTRUE;

  {
    Float_t dx = info->fSegSPD->Dx()*0.00005;
    Float_t dz = 3.50;

    fgSPDFrameBox = new TEveFrameBox();
    fgSPDFrameBox->SetAAQuadXZ(-dx, 0, -dz, 2*dx, 2*dz);
    fgSPDFrameBox->SetFrameColor((Color_t) 31);
    fgSPDFrameBox->SetFrameFill(kTRUE);
    fgSPDFrameBox->IncRefCount();
    fgSPDPalette  = new TEveRGBAPalette(info->fSPDMinVal,info->fSPDMaxVal);
    fgSPDPalette->IncRefCount();
  }

  {
    Float_t dx = info->fSegSDD->Dx()*0.0001;
    Float_t dz = info->fSegSDD->Dz()*0.00005;

    fgSDDFrameBox = new TEveFrameBox();
    fgSDDFrameBox->SetAAQuadXZ(-dx, 0, -dz, 2*dx, 2*dz);
    fgSDDFrameBox->SetFrameColor((Color_t) 32);
    fgSDDFrameBox->SetFrameFill(kTRUE);
    fgSDDFrameBox->IncRefCount();
    fgSDDPalette  = new TEveRGBAPalette(info->fSDDMinVal,info->fSDDMaxVal);
    fgSDDPalette->SetLimits(0, info->fSDDHighLim); // Set proper ADC range.
    fgSDDPalette->IncRefCount();
  }

  {
    Float_t dx = info->fSegSSD->Dx()*0.00005;
    Float_t dz = info->fSegSSD->Dz()*0.00005;

    fgSSDFrameBox = new TEveFrameBox();
    fgSSDFrameBox->SetAAQuadXZ(-dx, 0, -dz, 2*dx, 2*dz);
    fgSSDFrameBox->SetFrameColor((Color_t) 33);
    fgSSDFrameBox->SetFrameFill(kTRUE);
    fgSSDFrameBox->IncRefCount();
    fgSSDPalette  = new TEveRGBAPalette(info->fSSDMinVal,info->fSSDMaxVal);
    fgSSDPalette->SetLimits(0, info->fSSDHighLim); // Set proper ADC range.
    fgSSDPalette->IncRefCount();
  }

}

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

void AliEveITSModule::SetDigitsInfo(AliEveITSDigitsInfo* info)
{
  // Set data and geometry source.

  if (fInfo == info) return;
  if (fInfo) fInfo->DecRefCount();
  fInfo = info;
  if (fInfo) fInfo->IncRefCount();
}

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

void AliEveITSModule::SetID(Int_t gid, Bool_t trans)
{
  // Set detector id.

  static const TEveException eH("AliEveITSModule::SetID ");

  if (fInfo == 0)
    throw(eH + "AliEveITSDigitsInfo not set.");

  Int_t firstSPD = AliITSgeomTGeo::GetModuleIndex(1,1,1);
  Int_t lastSSD  = AliITSgeomTGeo::GetNModules() - 1;
  if (gid < firstSPD || gid > lastSSD)
  {
    throw(eH + Form("%d is not valid. ID range from %d to %d", gid,
		    firstSPD, lastSSD ));
  }

  fID = gid;

  if (!fgStaticInitDone) {
    InitStatics(fInfo);

    fgSPDFrameBox->IncRefCount(this);
    fgSPDPalette->IncRefCount();

    fgSDDFrameBox->IncRefCount(this);
    fgSDDPalette->IncRefCount();

    fgSSDFrameBox->IncRefCount(this);
    fgSSDPalette->IncRefCount();
  }

  AliITSgeomTGeo::GetModuleId(fID, fLayer, fLadder, fDet);
  TString strLadder = "Ladder";
  TString strSensor = "Sensor";
  TString symname;
  Int_t   id, nsector, nstave, nladder, rest;

  if (fID <= (AliITSgeomTGeo::GetModuleIndex(3,1,1) - 1))
  {
    // SPD

    SetFrame(fgSPDFrameBox);
    SetPalette(fgSPDPalette);

    symname += strLadder;
    if (fID < 80)
    {
      nsector = fID/8;
      rest    = fID - 8*nsector;
      nstave  = 1;
    }
    else
    {
      id      = fID - 80;
      nsector = id/8;
      rest    = id - 8*nsector;
      nstave  = 1;
    }
    if (rest < 4) nstave = 0;
    rest    -= 4*nstave;
    symname += rest;
    SetName(symname);
    fDetID = 0;
    fDx = fInfo->fSegSPD->Dx()*0.00005;
    fDz = 3.50;
    fDy = fInfo->fSegSPD->Dy()*0.00005;

  }
  else if (fID <= (AliITSgeomTGeo::GetModuleIndex(5,1,1) - 1))
  {
    // SDD

    SetFrame(fgSDDFrameBox);
    SetPalette(fgSDDPalette);

    symname += strSensor;
    if (fID < 324)
    {
      id      = fID - 240;
      nladder = id/6;
      rest    = id - 6*nladder;
    }
    else
    {
      id      = fID - 324;
      nladder = id/8;
      rest    = id - 8*nladder;
    }
    symname += rest;
    SetName(symname);
    fDetID = 1;
    fDx = fInfo->fSegSDD->Dx()*0.0001;
    fDz = fInfo->fSegSDD->Dz()*0.00005;
    fDy = fInfo->fSegSDD->Dy()*0.00005;

  }
  else
  {
    // SSD

    SetFrame(fgSSDFrameBox);
    SetPalette(fgSSDPalette);

    symname += strSensor;
    if (fID < 1248)
    {
      id      = fID - 500;
      nladder = id/22;
      rest    = id - nladder*22;
    }
    else
    {
      id      = fID - 1248;
      nladder = id/25;
      rest    = id - nladder*25;
    }
    symname += rest;
    SetName(symname);
    fDetID = 2;
    fInfo->fSegSSD->SetLayer(fLayer);
    fDx = fInfo->fSegSSD->Dx()*0.00005;
    fDz = fInfo->fSegSSD->Dz()*0.00005;
    fDy = fInfo->fSegSSD->Dy()*0.00005;

  }

  LoadQuads();
  ComputeBBox();
  if (trans)
    SetTrans();
}

void AliEveITSModule::LoadQuads()
{
  // Read module data from source and create low-level objects for
  // visualization - called quads.

  TClonesArray *digits  = fInfo->GetDigits(fID, fDetID);
  Int_t         ndigits = digits ? digits->GetEntriesFast() : 0;

  Float_t       x, z, dpx, dpz;
  Int_t         i, j;

  switch(fDetID)
  {

    case 0: { // SPD
      AliITSsegmentationSPD* seg =  fInfo->fSegSPD;

      Reset(kQT_RectangleXZFixedY, kFALSE, 32);

      for (Int_t k=0; k<ndigits; ++k)
      {
	AliITSdigit *d = (AliITSdigit*) digits->UncheckedAt(k);
	j = d->GetCoord1();
	i = d->GetCoord2();
	x  = -0.5*seg->Dx() + i*seg->Dpx(0);
	x *=  0.0001;
	fInfo->GetSPDLocalZ(j, z);
	dpx = seg->Dpx(i)*0.0001;
	dpz = seg->Dpz(j)*0.0001;

	AddQuad(x, z, dpx, dpz);
	QuadValue(1); // In principle could have color based on number of neigbours
	QuadId(d);
      }
      break;
    }

    case 1: { // SDD
      AliITSsegmentationSDD *seg =  fInfo->fSegSDD;

      Reset(kQT_RectangleXZFixedY, kFALSE, 32);

      for (Int_t k=0; k<ndigits; ++k)
      {
	AliITSdigit* d = (AliITSdigit*) digits->UncheckedAt(k);

	// if (d->GetSignal() > fgSDDThreshold)
	{
	  j = d->GetCoord1();
	  i = d->GetCoord2();
	  seg->DetToLocal(i, j, x, z);
	  dpx = seg->Dpx(i)*0.0001;
	  dpz = seg->Dpz(j)*0.0001;

	  AddQuad(x-2*dpx, z - dpz*0.5, 4*dpx, dpz);
	  QuadValue(d->GetSignal());
	  QuadId(d);
	}
      }
      break;
    }

    case 2: { // SSD
      AliITSsegmentationSSD* seg = fInfo->fSegSSD;

      Reset(kQT_LineXZFixedY, kFALSE, 32);

      Float_t ap, an; // positive/negative angles -> offsets
      seg->Angles(ap, an);
      ap =   TMath::Tan(ap) * fDz;
      an = - TMath::Tan(an) * fDz;

      for (Int_t k=0; k<ndigits; ++k)
      {
	AliITSdigit *d = (AliITSdigit*) digits->UncheckedAt(k);
	// if(d->GetSignal() > fgSSDThreshold)
	{
	  j = d->GetCoord1();
	  i = d->GetCoord2();
	  seg->DetToLocal(i,j,x,z);

	  Float_t a = ( d->GetCoord1() == 1) ? ap : an;

	  AddLine(x-a, -fDz, 2*a, 2*fDz);
	  QuadValue(d->GetSignal());
	  QuadId(d);
	  // printf("%3d -> %3d -> %8x\n", d->GetSignal(), ci, fQuads.back().color);
	}
      }
      break;
    }

  } // end switch

  RefitPlex();
}

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

void AliEveITSModule::SetTrans()
{
  // Set transformation matrix based on module id (use geometry to
  // retrieve this information).

  fHMTrans.SetFrom(*AliITSgeomTGeo::GetMatrix(fID));
}

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

void AliEveITSModule::DigitSelected(Int_t idx)
{
  // Override secondary select (alt-click) from TEveQuadSet.

  DigitBase_t *qb  = GetDigit(idx);
  TObject     *obj = qb->fId.GetObject();
  AliITSdigit *d   = dynamic_cast<AliITSdigit*>(obj);
  printf("AliEveITSModule::QuadSelected "); Print();
  printf("  idx=%d, value=%d, obj=0x%lx, digit=0x%lx\n",
	 idx, qb->fValue, (ULong_t)obj, (ULong_t)d);
  if (d)
    printf("  coord1=%3d coord2=%3d signal=%d\n",
	   d->GetCoord1(), d->GetCoord2(), d->GetSignal());

}

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

void AliEveITSModule::Print(Option_t* ) const
{
  // Print object summary information.

  printf("AliEveITSModule: ID %d, layer %d, ladder %d, det %d\n",
         fID, fLayer, fLadder, fDetID);
}
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	**************************************************************************/
0879c50b	9
d810d0de	10	#include "AliEveITSModule.h"
5a5a1232	11
12f6a2f3	12	#include <AliITSgeomTGeo.h>
0879c50b	13	#include <AliITSsegmentationSPD.h>
	14	#include <AliITSsegmentationSDD.h>
	15	#include <AliITSsegmentationSSD.h>
	16
5a5a1232	17	#include <AliITSdigitSPD.h>
	18	#include <AliITSdigitSDD.h>
	19	#include <AliITSdigitSSD.h>
	20
	21	#include <TStyle.h>
5a5a1232	22
9efc978d	23
698e2c9b	24	//______________________________________________________________________________
	25	//
	26	// Visualization of an ITS module.
	27
	28	ClassImp(AliEveITSModule)
	29
0879c50b	30	Bool_t AliEveITSModule::fgStaticInitDone = kFALSE;
d810d0de	31
	32	TEveFrameBox* AliEveITSModule::fgSPDFrameBox = 0;
	33	TEveFrameBox* AliEveITSModule::fgSDDFrameBox = 0;
	34	TEveFrameBox* AliEveITSModule::fgSSDFrameBox = 0;
9efc978d	35
d810d0de	36	TEveRGBAPalette* AliEveITSModule::fgSPDPalette = 0;
	37	TEveRGBAPalette* AliEveITSModule::fgSDDPalette = 0;
	38	TEveRGBAPalette* AliEveITSModule::fgSSDPalette = 0;
9efc978d	39
57ffa5fb	40	/******************************************************************************/
5a5a1232	41
d810d0de	42	AliEveITSModule::AliEveITSModule(const Text_t* n, const Text_t* t) :
84aff7a4	43	TEveQuadSet(n, t),
265ecb21	44	fInfo(0),
	45	fID(-1), fDetID(-1),
	46	fLayer(-1), fLadder(-1), fDet(-1),
9efc978d	47	fDx(0), fDz(0), fDy(0)
698e2c9b	48	{
	49	// Constructor.
	50	}
5a5a1232	51
d810d0de	52	AliEveITSModule::AliEveITSModule(Int_t gid, AliEveITSDigitsInfo* info) :
84aff7a4	53	TEveQuadSet(Form("ITS module %d", gid)),
27db2029	54	fInfo (0),
265ecb21	55	fID(-1), fDetID(-1),
265ecb21	56	fLayer(-1), fLadder(-1), fDet(-1),
9efc978d	57	fDx(0), fDz(0), fDy(0)
5a5a1232	58	{
698e2c9b	59	// Constructor with module id and data-source.
698e2c9b	60
27db2029	61	SetDigitsInfo(info);
9c8221e9	62	SetID(gid);
5a5a1232	63	}
5a5a1232	64
d810d0de	65	AliEveITSModule::~AliEveITSModule()
5a5a1232	66	{
698e2c9b	67	// Destructor.
698e2c9b	68
5a5a1232	69	if(fInfo) fInfo->DecRefCount();
	70	}
	71
57ffa5fb	72	/******************************************************************************/
5a5a1232	73
d810d0de	74	void AliEveITSModule::InitStatics(AliEveITSDigitsInfo* info)
9efc978d	75	{
698e2c9b	76	// Initialize static variables,
698e2c9b	77
9efc978d	78	if (fgStaticInitDone) return;
9c8221e9	79	fgStaticInitDone = kTRUE;
9efc978d	80
	81	{
	82	Float_t dx = info->fSegSPD->Dx()*0.00005;
51346b82	83	Float_t dz = 3.50;
9efc978d	84
84aff7a4	85	fgSPDFrameBox = new TEveFrameBox();
9efc978d	86	fgSPDFrameBox->SetAAQuadXZ(-dx, 0, -dz, 2dx, 2dz);
9efc978d	87	fgSPDFrameBox->SetFrameColor((Color_t) 31);
32e219c2	88	fgSPDFrameBox->SetFrameFill(kTRUE);
32e219c2	89	fgSPDFrameBox->IncRefCount();
84aff7a4	90	fgSPDPalette = new TEveRGBAPalette(info->fSPDMinVal,info->fSPDMaxVal);
32e219c2	91	fgSPDPalette->IncRefCount();
9efc978d	92	}
	93
	94	{
	95	Float_t dx = info->fSegSDD->Dx()*0.0001;
	96	Float_t dz = info->fSegSDD->Dz()*0.00005;
	97
84aff7a4	98	fgSDDFrameBox = new TEveFrameBox();
9efc978d	99	fgSDDFrameBox->SetAAQuadXZ(-dx, 0, -dz, 2dx, 2dz);
9efc978d	100	fgSDDFrameBox->SetFrameColor((Color_t) 32);
32e219c2	101	fgSDDFrameBox->SetFrameFill(kTRUE);
32e219c2	102	fgSDDFrameBox->IncRefCount();
84aff7a4	103	fgSDDPalette = new TEveRGBAPalette(info->fSDDMinVal,info->fSDDMaxVal);
32e219c2	104	fgSDDPalette->SetLimits(0, info->fSDDHighLim); // Set proper ADC range.
32e219c2	105	fgSDDPalette->IncRefCount();
9efc978d	106	}
	107
	108	{
	109	Float_t dx = info->fSegSSD->Dx()*0.00005;
	110	Float_t dz = info->fSegSSD->Dz()*0.00005;
	111
84aff7a4	112	fgSSDFrameBox = new TEveFrameBox();
9efc978d	113	fgSSDFrameBox->SetAAQuadXZ(-dx, 0, -dz, 2dx, 2dz);
9efc978d	114	fgSSDFrameBox->SetFrameColor((Color_t) 33);
32e219c2	115	fgSSDFrameBox->SetFrameFill(kTRUE);
32e219c2	116	fgSSDFrameBox->IncRefCount();
84aff7a4	117	fgSSDPalette = new TEveRGBAPalette(info->fSSDMinVal,info->fSSDMaxVal);
32e219c2	118	fgSSDPalette->SetLimits(0, info->fSSDHighLim); // Set proper ADC range.
32e219c2	119	fgSSDPalette->IncRefCount();
9efc978d	120	}
	121
	122	}
	123
57ffa5fb	124	/******************************************************************************/
9efc978d	125
d810d0de	126	void AliEveITSModule::SetDigitsInfo(AliEveITSDigitsInfo* info)
5a5a1232	127	{
698e2c9b	128	// Set data and geometry source.
698e2c9b	129
9c8221e9	130	if (fInfo == info) return;
	131	if (fInfo) fInfo->DecRefCount();
	132	fInfo = info;
	133	if (fInfo) fInfo->IncRefCount();
	134	}
	135
57ffa5fb	136	/******************************************************************************/
9c8221e9	137
d810d0de	138	void AliEveITSModule::SetID(Int_t gid, Bool_t trans)
9c8221e9	139	{
698e2c9b	140	// Set detector id.
698e2c9b	141
d810d0de	142	static const TEveException eH("AliEveITSModule::SetID ");
9efc978d	143
12f6a2f3	144	if (fInfo == 0)
d810d0de	145	throw(eH + "AliEveITSDigitsInfo not set.");
9c8221e9	146
ff67a00a	147	Int_t firstSPD = AliITSgeomTGeo::GetModuleIndex(1,1,1);
	148	Int_t lastSSD = AliITSgeomTGeo::GetNModules() - 1;
	149	if (gid < firstSPD \|\| gid > lastSSD)
12f6a2f3	150	{
9c8221e9	151	throw(eH + Form("%d is not valid. ID range from %d to %d", gid,
12f6a2f3	152	firstSPD, lastSSD ));
12f6a2f3	153	}
9c8221e9	154
	155	fID = gid;
	156
1d3c2b09	157	if (!fgStaticInitDone) {
1d3c2b09	158	InitStatics(fInfo);
51346b82	159
1d3c2b09	160	fgSPDFrameBox->IncRefCount(this);
1d3c2b09	161	fgSPDPalette->IncRefCount();
51346b82	162
1d3c2b09	163	fgSDDFrameBox->IncRefCount(this);
	164	fgSDDPalette->IncRefCount();
	165
	166	fgSSDFrameBox->IncRefCount(this);
	167	fgSSDPalette->IncRefCount();
	168	}
9c8221e9	169
12f6a2f3	170	AliITSgeomTGeo::GetModuleId(fID, fLayer, fLadder, fDet);
ba6cfb13	171	TString strLadder = "Ladder";
	172	TString strSensor = "Sensor";
	173	TString symname;
9c8221e9	174	Int_t id, nsector, nstave, nladder, rest;
51346b82	175
12f6a2f3	176	if (fID <= (AliITSgeomTGeo::GetModuleIndex(3,1,1) - 1))
9efc978d	177	{
	178	// SPD
	179
9c8221e9	180	SetFrame(fgSPDFrameBox);
9c8221e9	181	SetPalette(fgSPDPalette);
51346b82	182
9efc978d	183	symname += strLadder;
	184	if (fID < 80)
	185	{
ba6cfb13	186	nsector = fID/8;
9efc978d	187	rest = fID - 8*nsector;
	188	nstave = 1;
	189	}
	190	else
	191	{
	192	id = fID - 80;
ba6cfb13	193	nsector = id/8;
9efc978d	194	rest = id - 8*nsector;
9efc978d	195	nstave = 1;
ba6cfb13	196	}
9efc978d	197	if (rest < 4) nstave = 0;
	198	rest -= 4*nstave;
	199	symname += rest;
	200	SetName(symname);
	201	fDetID = 0;
	202	fDx = fInfo->fSegSPD->Dx()*0.00005;
51346b82	203	fDz = 3.50;
9efc978d	204	fDy = fInfo->fSegSPD->Dy()*0.00005;
9efc978d	205
5a5a1232	206	}
12f6a2f3	207	else if (fID <= (AliITSgeomTGeo::GetModuleIndex(5,1,1) - 1))
9efc978d	208	{
	209	// SDD
	210
9c8221e9	211	SetFrame(fgSDDFrameBox);
9c8221e9	212	SetPalette(fgSDDPalette);
51346b82	213
9efc978d	214	symname += strSensor;
	215	if (fID < 324)
	216	{
	217	id = fID - 240;
ba6cfb13	218	nladder = id/6;
9efc978d	219	rest = id - 6*nladder;
	220	}
	221	else
	222	{
	223	id = fID - 324;
ba6cfb13	224	nladder = id/8;
9efc978d	225	rest = id - 8*nladder;
ba6cfb13	226	}
9efc978d	227	symname += rest;
	228	SetName(symname);
	229	fDetID = 1;
	230	fDx = fInfo->fSegSDD->Dx()*0.0001;
	231	fDz = fInfo->fSegSDD->Dz()*0.00005;
	232	fDy = fInfo->fSegSDD->Dy()*0.00005;
	233
5a5a1232	234	}
9efc978d	235	else
	236	{
	237	// SSD
	238
9c8221e9	239	SetFrame(fgSSDFrameBox);
	240	SetPalette(fgSSDPalette);
	241
9efc978d	242	symname += strSensor;
	243	if (fID < 1248)
	244	{
	245	id = fID - 500;
ba6cfb13	246	nladder = id/22;
9efc978d	247	rest = id - nladder*22;
	248	}
	249	else
	250	{
	251	id = fID - 1248;
ba6cfb13	252	nladder = id/25;
9efc978d	253	rest = id - nladder*25;
ba6cfb13	254	}
9efc978d	255	symname += rest;
	256	SetName(symname);
	257	fDetID = 2;
51346b82	258	fInfo->fSegSSD->SetLayer(fLayer);
9efc978d	259	fDx = fInfo->fSegSSD->Dx()*0.00005;
	260	fDz = fInfo->fSegSSD->Dz()*0.00005;
	261	fDy = fInfo->fSegSSD->Dy()*0.00005;
	262
5a5a1232	263	}
5a5a1232	264
51346b82	265	LoadQuads();
5a5a1232	266	ComputeBBox();
97e1e1bc	267	if (trans)
97e1e1bc	268	SetTrans();
5a5a1232	269	}
5a5a1232	270
d810d0de	271	void AliEveITSModule::LoadQuads()
5a5a1232	272	{
698e2c9b	273	// Read module data from source and create low-level objects for
698e2c9b	274	// visualization - called quads.
9efc978d	275
b3bc5327	276	TClonesArray *digits = fInfo->GetDigits(fID, fDetID);
	277	Int_t ndigits = digits ? digits->GetEntriesFast() : 0;
	278
51346b82	279	Float_t x, z, dpx, dpz;
b3bc5327	280	Int_t i, j;
5a5a1232	281
9efc978d	282	switch(fDetID)
	283	{
	284
	285	case 0: { // SPD
51346b82	286	AliITSsegmentationSPD* seg = fInfo->fSegSPD;
9efc978d	287
84aff7a4	288	Reset(kQT_RectangleXZFixedY, kFALSE, 32);
5a5a1232	289
9efc978d	290	for (Int_t k=0; k<ndigits; ++k)
9efc978d	291	{
b3bc5327	292	AliITSdigit d = (AliITSdigit) digits->UncheckedAt(k);
5a5a1232	293	j = d->GetCoord1();
5a5a1232	294	i = d->GetCoord2();
9efc978d	295	x = -0.5seg->Dx() + iseg->Dpx(0);
	296	x *= 0.0001;
	297	fInfo->GetSPDLocalZ(j, z);
5a5a1232	298	dpx = seg->Dpx(i)*0.0001;
	299	dpz = seg->Dpz(j)*0.0001;
	300
9efc978d	301	AddQuad(x, z, dpx, dpz);
9efc978d	302	QuadValue(1); // In principle could have color based on number of neigbours
ffa36681	303	QuadId(d);
5a5a1232	304	}
9efc978d	305	break;
5a5a1232	306	}
5a5a1232	307
9efc978d	308	case 1: { // SDD
51346b82	309	AliITSsegmentationSDD *seg = fInfo->fSegSDD;
5a5a1232	310
84aff7a4	311	Reset(kQT_RectangleXZFixedY, kFALSE, 32);
9efc978d	312
	313	for (Int_t k=0; k<ndigits; ++k)
	314	{
b3bc5327	315	AliITSdigit* d = (AliITSdigit*) digits->UncheckedAt(k);
9efc978d	316
	317	// if (d->GetSignal() > fgSDDThreshold)
	318	{
	319	j = d->GetCoord1();
	320	i = d->GetCoord2();
	321	seg->DetToLocal(i, j, x, z);
	322	dpx = seg->Dpx(i)*0.0001;
	323	dpz = seg->Dpz(j)*0.0001;
	324
97e1e1bc	325	AddQuad(x-2dpx, z - dpz0.5, 4*dpx, dpz);
9efc978d	326	QuadValue(d->GetSignal());
ffa36681	327	QuadId(d);
5a5a1232	328	}
9efc978d	329	}
	330	break;
	331	}
	332
	333	case 2: { // SSD
51346b82	334	AliITSsegmentationSSD* seg = fInfo->fSegSSD;
9efc978d	335
84aff7a4	336	Reset(kQT_LineXZFixedY, kFALSE, 32);
9efc978d	337
	338	Float_t ap, an; // positive/negative angles -> offsets
	339	seg->Angles(ap, an);
	340	ap = TMath::Tan(ap) * fDz;
	341	an = - TMath::Tan(an) * fDz;
	342
	343	for (Int_t k=0; k<ndigits; ++k)
	344	{
b3bc5327	345	AliITSdigit d = (AliITSdigit) digits->UncheckedAt(k);
9efc978d	346	// if(d->GetSignal() > fgSSDThreshold)
	347	{
	348	j = d->GetCoord1();
	349	i = d->GetCoord2();
	350	seg->DetToLocal(i,j,x,z);
	351
	352	Float_t a = ( d->GetCoord1() == 1) ? ap : an;
	353
e974c097	354	AddLine(x-a, -fDz, 2a, 2fDz);
9efc978d	355	QuadValue(d->GetSignal());
ffa36681	356	QuadId(d);
9efc978d	357	// printf("%3d -> %3d -> %8x\n", d->GetSignal(), ci, fQuads.back().color);
5a5a1232	358	}
5a5a1232	359	}
9efc978d	360	break;
5a5a1232	361	}
5a5a1232	362
9efc978d	363	} // end switch
	364
	365	RefitPlex();
5a5a1232	366	}
5a5a1232	367
57ffa5fb	368	/******************************************************************************/
5a5a1232	369
d810d0de	370	void AliEveITSModule::SetTrans()
5a5a1232	371	{
4852ff6f	372	// Set transformation matrix based on module id (use geometry to
	373	// retrieve this information).
	374
12f6a2f3	375	fHMTrans.SetFrom(*AliITSgeomTGeo::GetMatrix(fID));
5a5a1232	376	}
5a5a1232	377
57ffa5fb	378	/******************************************************************************/
5a5a1232	379
d810d0de	380	void AliEveITSModule::DigitSelected(Int_t idx)
ffa36681	381	{
12f6a2f3	382	// Override secondary select (alt-click) from TEveQuadSet.
ffa36681	383
12f6a2f3	384	DigitBase_t *qb = GetDigit(idx);
	385	TObject *obj = qb->fId.GetObject();
	386	AliITSdigit d = dynamic_cast<AliITSdigit>(obj);
d810d0de	387	printf("AliEveITSModule::QuadSelected "); Print();
ffa36681	388	printf(" idx=%d, value=%d, obj=0x%lx, digit=0x%lx\n",
	389	idx, qb->fValue, (ULong_t)obj, (ULong_t)d);
	390	if (d)
	391	printf(" coord1=%3d coord2=%3d signal=%d\n",
	392	d->GetCoord1(), d->GetCoord2(), d->GetSignal());
	393
	394	}
	395
57ffa5fb	396	/******************************************************************************/
ffa36681	397
d810d0de	398	void AliEveITSModule::Print(Option_t* ) const
5a5a1232	399	{
698e2c9b	400	// Print object summary information.
	401
	402	printf("AliEveITSModule: ID %d, layer %d, ladder %d, det %d\n",
	403	fID, fLayer, fLadder, fDetID);
5a5a1232	404	}