[u/mrichter/AliRoot.git] / RICH / AliRICHDigit.cxx

//  **************************************************************************
//  * Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved. *
//  *                                                                        *
//  * Author: The ALICE Off-line Project.                                    *
//  * Contributors are mentioned in the code where appropriate.              *
//  *                                                                        *
//  * Permission to use, copy, modify and distribute this software and its   *
//  * documentation strictly for non-commercial purposes is hereby granted   *
//  * without fee, provided that the above copyright notice appears in all   *
//  * copies and that both the copyright notice and this permission notice   *
//  * appear in the supporting documentation. The authors make no claims     *
//  * about the suitability of this software for any purpose. It is          *
//  * provided "as is" without express or implied warranty.                  *
//  **************************************************************************

#include "AliRICHDigit.h" //class header
#include "AliRICHHit.h"   //Hit2Sdi() 
#include <TClonesArray.h> //Hit2Sdi() 
#include <TCanvas.h>      //TestSeg()
#include <TLine.h>        //TestSeg()
#include <TLatex.h>       //TestSeg()
#include <TPolyLine.h>    //DrawPc() 
#include <TGStatusBar.h>  //Zoom()
#include <TRootCanvas.h>  //Zoom() 
ClassImp(AliRICHDigit)

/*
Any given LDC collects data from a number of D-RORC cards connected to this same LDC by separate DDLs. This data is stored in corresponding number
of DDL buffers.

Each DDL buffer corresponds to a single D-RORC. The data this buffer contains are hardware generated by D-RORC. The buffer starts with common header 
which size and structure is standartized and mandatory for all detectors.
The header contains among other words, so called Equipment ID word. This unique value for each D-RORC is calculated as detector ID << 8 + DDL index. 
For RICH the detector ID is 6 (reffered in the code as kRichRawId) while DDL indexes are from 0 to 13.

Common header might be followed by the private one although  RICH has no any private header, just uses the common one.

Single RICH D-RORC serves one half of a chamber i.e. 3 photocathodes even LDC for left part( PCs 0-2-4) and odd LDC for right part(1-3-5) 
as it's seen from electronics side.

So the LDC -chamber-ddl map is:
DDL index  0 -> ch 1 left -> DDL ID 0x600          DDL index  1 -> ch 1 right -> DDL ID 0x601 
DDL index  2 -> ch 2 left -> DDL ID 0x602          DDL index  3 -> ch 2 right -> DDL ID 0x603 
DDL index  4 -> ch 3 left -> DDL ID 0x604          DDL index  5 -> ch 3 right -> DDL ID 0x605 
DDL index  6 -> ch 4 left -> DDL ID 0x606          DDL index  7 -> ch 4 right -> DDL ID 0x607 
DDL index  8 -> ch 5 left -> DDL ID 0x608          DDL index  9 -> ch 5 right -> DDL ID 0x609 
DDL index 10 -> ch 6 left -> DDL ID 0x60a          DDL index 11 -> ch 6 right -> DDL ID 0x60b 
DDL index 12 -> ch 7 left -> DDL ID 0x60c          DDL index 13 -> ch 7 right -> DDL ID 0x60d 

RICH FEE as seen by single D-RORC is composed from a number of DILOGIC chips organized in vertical stack of rows. 
Each DILOGIC chip serves 48 channels for the 8x6 pads (reffered in the code as kDilX,kDilY). Channels counted from 0 to 47.

The mapping inside DILOGIC chip has the following structure (see from electronics side):

5  04 10 16 22 28 34 40 46
4  02 08 14 20 26 32 38 44
3  00 06 12 18 24 30 36 42
2  01 07 13 19 25 31 37 43
1  03 09 15 21 27 33 39 45
0  05 11 17 23 29 35 41 47
    
    0  1  2  3  4  5  6  7  padx

10 DILOGIC chips composes so called "row" in horizontal direction (reffered in the code as kNdil), so the row is 80x6 pads structure. 
DILOGIC chips in the row are counted from right to left as seen from electronics side, from 1 to 10.
24 rows are piled up forming the whole FEE served by single D-RORC, so one DDL sees 80x144 pads separated in 3 photocathodes.
Rows are counted from 1 to 24 from top    to bottom for right  half of the chamber (PCs 1-3-5) as seen from electronics side, meaning even LDC number
                          and from bottom to top    for left   half of the chamber (PCs 0-2-4) as seen from electronics side, meaning odd LDC number.

RICH raw word is 32 bits with the structure:   
   00000             rrrrr                      dddd                               aaaaaa                          qqqqqqqqqqqq 
 5 bits zero  5 bits row number (1..24)  4 bits DILOGIC chip number (1..10) 6 bits DILOGIC address (0..47)  12 bits QDC value (0..4095)
*/
//++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
Float_t AliRICHDigit::Hit2Sdi(AliRICHHit *pHit,TClonesArray *pSdiLst)
{
// Creates a list of sdigits out of provided hit
// Arguments: pHit- hit
//   Returns: none
  
  Float_t  x=(pHit->LorsX() > SizePcX())? pHit->LorsX()-SizePcX()-SizeDead():pHit->LorsX();          //sagita is for PC (0-64) and not for chamber   
  Float_t  qdcEle=34.06311+0.2337070*x+5.807476e-3*x*x-2.956471e-04*x*x*x+2.310001e-06*x*x*x*x;      //reparametrised from DiMauro
  
  Int_t iNele=Int_t(pHit->E()/26e-9);  if(iNele<1) iNele=1;                                          //number of electrons created by hit
  Float_t qdcTot=0; for(Int_t i=1;i<=iNele;i++) qdcTot-=qdcEle*TMath::Log(gRandom->Rndm()+1e-6);     //total qdc fro hit, 1e-6 is a protection against 0 from rndm
  
  AliRICHDigit dd(1,1,1,pHit->LorsX(),pHit->LorsY());                                                //tmp digit to shift hit y to the nearest anod wire  
  Float_t y=  (pHit->LorsY() > dd.LorsY()) ? dd.LorsY()+0.21 :  dd.LorsY()-0.21;
  
  Int_t iSdiCnt=pSdiLst->GetEntries();
  for(Int_t i=0;i<9;i++){                                      //affected pads loop
    AliRICHDigit dig(pHit->Ch(),qdcTot,pHit->Tid(),pHit->LorsX(),y,i); //c,q,tid,x,y   create tmp sdigit for pad i around hit position
    if(dig.PadX()==-1) continue;
    new((*pSdiLst)[iSdiCnt++]) AliRICHDigit(dig);
  }
  return qdcTot;
}//Hit2Sdi()
//++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
void AliRICHDigit::Print(Option_t*)const
{
// Print current digit  
// Arguments: option string not used
//   Returns: none    
  Printf("pad=(%2i,%2i,%3i,%3i),pos=(%7.2f,%7.2f) QDC=%8.3f, TID=(%5i,%5i,%5i) raw r=%2i d=%2i a=%2i",
               Ch(),Pc(),PadX(),PadY(),LorsX(),LorsY(), Q(),  fTracks[0],fTracks[1],fTracks[2], Row(), Dilogic(), Addr());
}
//++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
void AliRICHDigit::PrintSize()
{
// Print all segmentaion related sizes
// Arguments: none
//   Returns: none    
  Printf("-->pad    =(%6.2f,%6.2f) cm dead zone %.2f cm\n"
         "-->PC     =(%6.2f,%6.2f) cm\n"
         "-->all PCs=(%6.2f,%6.2f) cm",
               SizePadX(),SizePadY(),SizeDead(),SizePcX(),SizePcY(),SizeAllX(),SizeAllY());
}
//++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
void AliRICHDigit::TestSeg()
{
// Draws the picture of segmentation   
// Arguments: none
//   Returns: none     
  TCanvas *pC=new TCanvas("pads","View from electronics side, IP is behind the picture.");pC->ToggleEventStatus();
  gPad->AddExec("test","AliRICHDigit::Zoom()");
  DrawPc();  
}//TestSeg()
//++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
void AliRICHDigit::Zoom()
{
// Show info about current cursur position in status bar of the canvas
// Arguments: none
//   Returns: none     
  TCanvas *pC=(TCanvas*)gPad; 
  TRootCanvas *pRC= (TRootCanvas*)pC->GetCanvasImp();
  TGStatusBar *pBar=pRC->GetStatusBar();
  pBar->SetParts(5);
  Float_t x=gPad->AbsPixeltoX(gPad->GetEventX());
  Float_t y=gPad->AbsPixeltoY(gPad->GetEventY());
  AliRICHDigit dig(1,100,1,x,y);  
  if(IsInDead(x,y))
    pBar->SetText("Out of sensitive area",4);    
  else
    pBar->SetText(Form("p%i x%i y%i r%i d%i a%i (%.2f,%.2f)",dig.Pc(),dig.PadX(),dig.PadY(),dig.Row(),dig.Dilogic(),dig.Addr(),dig.LorsX(),dig.LorsY()),4);
  if(gPad->GetEvent()==1){
    new TCanvas("zoom",Form("Row %i DILOGIC %i",dig.Row(),dig.Dilogic()));  
  }
}//Zoom()
//++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
void AliRICHDigit::DrawPc(Bool_t isFill) 
{ 
// Utility methode draws RICH chamber PCs on event display.
// Arguments: none
//   Returns: none      
//  y6  ----------  ----------
//      |        |  |        |
//      |    4   |  |    5   |
//  y5  ----------  ----------
//
//  y4  ----------  ----------
//      |        |  |        |
//      |    2   |  |    3   |   view from electronics side
//  y3  ----------  ----------
//          
//  y2  ----------  ----------
//      |        |  |        |
//      |    0   |  |    1   |
//  y1  ----------  ----------
//      x1      x2  x3       x4
  gPad->Range(-5,-5,SizeAllX()+5,SizeAllY()+5); 
  Float_t x1=0,x2=SizePcX(),x3=SizePcX()+SizeDead(),                                                  x4=SizeAllX();
  Float_t y1=0,y2=SizePcY(),y3=SizePcY()+SizeDead(),y4=2*SizePcY()+SizeDead(),y5=SizeAllY()-SizePcY(),y6=SizeAllY();

  Float_t xL[5]={x1,x1,x2,x2,x1}; //clockwise
  Float_t xR[5]={x3,x3,x4,x4,x3};  
  Float_t yD[5]={y1,y2,y2,y1,y1};
  Float_t yC[5]={y3,y4,y4,y3,y3};  
  Float_t yU[5]={y5,y6,y6,y5,y5};
    
  TLatex t; t.SetTextSize(0.01); t.SetTextAlign(22);
  Int_t iColLeft=29,iColRight=41;
  TPolyLine *pc=0;  TLine *pL; Float_t x0=0,y0=0,wRow=kDilY*SizePadY(),wDil=kDilX*SizePadX();
  for(Int_t iPc=0;iPc<kPcAll;iPc++){
    if(iPc==4) {pc=new TPolyLine(5,xL,yU);t.DrawText(x1-3.5,y6-20,"PC4");x0=x1;y0=y5;} if(iPc==5) {pc=new TPolyLine(5,xR,yU);t.DrawText(x4+3.5,y6-20,"PC5");x0=x3;y0=y5;}
    if(iPc==2) {pc=new TPolyLine(5,xL,yC);t.DrawText(x1-3.5,y4-20,"PC2");x0=x1;y0=y3;} if(iPc==3) {pc=new TPolyLine(5,xR,yC);t.DrawText(x4+3.5,y4-20,"PC3");x0=x3;y0=y3;}
    if(iPc==0) {pc=new TPolyLine(5,xL,yD);t.DrawText(x1-3.5,y2-20,"PC0");x0=x1;y0=y1;} if(iPc==1) {pc=new TPolyLine(5,xR,yD);t.DrawText(x4+3.5,y2-20,"PC1");x0=x3;y0=y1;}    
    (iPc%2)? pc->SetFillColor(iColLeft): pc->SetFillColor(iColRight);
    if(isFill) pc->Draw("f"); else pc->Draw();
    for(Int_t i=1;i<=8 ;i++){//draw row lines (horizontal)
      Float_t y=y0+i*wRow;
      Int_t row=i+iPc/2*8; if(iPc%2!=0) row=25-row; t.DrawText(x0-1,y -3,Form("r%i",row)); 
      if(i==8) break;                                      //do not draw the last line of PC
      pL=new TLine(x0,y,x0+SizePcX(),y);   pL->Draw(); 
    }    
    for(Int_t iDil=1;iDil<=10;iDil++){Float_t x=x0+iDil*wDil;t.DrawText(x -3,y0-1,Form("d%i",11-iDil));   if(iDil==10) break;    pL=new TLine(x,y0,x,y0+SizePcY());  pL->Draw();}  
  }      
}//DrawPc()
//++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
Commit	Line	Data
0fe8fa07	1	// **************************************************************************
	2	// * Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved. *
	3	// * *
	4	// * Author: The ALICE Off-line Project. *
	5	// * Contributors are mentioned in the code where appropriate. *
	6	// * *
	7	// * Permission to use, copy, modify and distribute this software and its *
	8	// * documentation strictly for non-commercial purposes is hereby granted *
	9	// * without fee, provided that the above copyright notice appears in all *
	10	// * copies and that both the copyright notice and this permission notice *
	11	// * appear in the supporting documentation. The authors make no claims *
	12	// * about the suitability of this software for any purpose. It is *
	13	// * provided "as is" without express or implied warranty. *
	14	// **************************************************************************
	15
e30ca504	16	#include "AliRICHDigit.h" //class header
	17	#include "AliRICHHit.h" //Hit2Sdi()
	18	#include <TClonesArray.h> //Hit2Sdi()
	19	#include <TCanvas.h> //TestSeg()
	20	#include <TLine.h> //TestSeg()
	21	#include <TLatex.h> //TestSeg()
	22	#include <TPolyLine.h> //DrawPc()
	23	#include <TGStatusBar.h> //Zoom()
	24	#include <TRootCanvas.h> //Zoom()
0fe8fa07	25	ClassImp(AliRICHDigit)
0fe8fa07	26
9de8d990	27	/*
	28	Any given LDC collects data from a number of D-RORC cards connected to this same LDC by separate DDLs. This data is stored in corresponding number
	29	of DDL buffers.
	30
	31	Each DDL buffer corresponds to a single D-RORC. The data this buffer contains are hardware generated by D-RORC. The buffer starts with common header
	32	which size and structure is standartized and mandatory for all detectors.
	33	The header contains among other words, so called Equipment ID word. This unique value for each D-RORC is calculated as detector ID << 8 + DDL index.
	34	For RICH the detector ID is 6 (reffered in the code as kRichRawId) while DDL indexes are from 0 to 13.
	35
	36	Common header might be followed by the private one although RICH has no any private header, just uses the common one.
	37
e30ca504	38	Single RICH D-RORC serves one half of a chamber i.e. 3 photocathodes even LDC for left part( PCs 0-2-4) and odd LDC for right part(1-3-5)
9de8d990	39	as it's seen from electronics side.
	40
	41	So the LDC -chamber-ddl map is:
	42	DDL index 0 -> ch 1 left -> DDL ID 0x600 DDL index 1 -> ch 1 right -> DDL ID 0x601
	43	DDL index 2 -> ch 2 left -> DDL ID 0x602 DDL index 3 -> ch 2 right -> DDL ID 0x603
	44	DDL index 4 -> ch 3 left -> DDL ID 0x604 DDL index 5 -> ch 3 right -> DDL ID 0x605
	45	DDL index 6 -> ch 4 left -> DDL ID 0x606 DDL index 7 -> ch 4 right -> DDL ID 0x607
	46	DDL index 8 -> ch 5 left -> DDL ID 0x608 DDL index 9 -> ch 5 right -> DDL ID 0x609
	47	DDL index 10 -> ch 6 left -> DDL ID 0x60a DDL index 11 -> ch 6 right -> DDL ID 0x60b
	48	DDL index 12 -> ch 7 left -> DDL ID 0x60c DDL index 13 -> ch 7 right -> DDL ID 0x60d
	49
	50	RICH FEE as seen by single D-RORC is composed from a number of DILOGIC chips organized in vertical stack of rows.
e30ca504	51	Each DILOGIC chip serves 48 channels for the 8x6 pads (reffered in the code as kDilX,kDilY). Channels counted from 0 to 47.
	52
	53	The mapping inside DILOGIC chip has the following structure (see from electronics side):
9de8d990	54
e30ca504	55	5 04 10 16 22 28 34 40 46
	56	4 02 08 14 20 26 32 38 44
	57	3 00 06 12 18 24 30 36 42
	58	2 01 07 13 19 25 31 37 43
	59	1 03 09 15 21 27 33 39 45
	60	0 05 11 17 23 29 35 41 47
	61
	62	0 1 2 3 4 5 6 7 padx
9de8d990	63
e30ca504	64	10 DILOGIC chips composes so called "row" in horizontal direction (reffered in the code as kNdil), so the row is 80x6 pads structure.
9de8d990	65	DILOGIC chips in the row are counted from right to left as seen from electronics side, from 1 to 10.
e30ca504	66	24 rows are piled up forming the whole FEE served by single D-RORC, so one DDL sees 80x144 pads separated in 3 photocathodes.
	67	Rows are counted from 1 to 24 from top to bottom for right half of the chamber (PCs 1-3-5) as seen from electronics side, meaning even LDC number
	68	and from bottom to top for left half of the chamber (PCs 0-2-4) as seen from electronics side, meaning odd LDC number.
9de8d990	69
	70	RICH raw word is 32 bits with the structure:
	71	00000 rrrrr dddd aaaaaa qqqqqqqqqqqq
	72	5 bits zero 5 bits row number (1..24) 4 bits DILOGIC chip number (1..10) 6 bits DILOGIC address (0..47) 12 bits QDC value (0..4095)
	73	*/
e30ca504	74	//++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
	75	Float_t AliRICHDigit::Hit2Sdi(AliRICHHit pHit,TClonesArray pSdiLst)
	76	{
	77	// Creates a list of sdigits out of provided hit
	78	// Arguments: pHit- hit
	79	// Returns: none
	80
	81	Float_t x=(pHit->LorsX() > SizePcX())? pHit->LorsX()-SizePcX()-SizeDead():pHit->LorsX(); //sagita is for PC (0-64) and not for chamber
	82	Float_t qdcEle=34.06311+0.2337070x+5.807476e-3xx-2.956471e-04xxx+2.310001e-06xxxx; //reparametrised from DiMauro
	83
	84	Int_t iNele=Int_t(pHit->E()/26e-9); if(iNele<1) iNele=1; //number of electrons created by hit
	85	Float_t qdcTot=0; for(Int_t i=1;i<=iNele;i++) qdcTot-=qdcEle*TMath::Log(gRandom->Rndm()+1e-6); //total qdc fro hit, 1e-6 is a protection against 0 from rndm
	86
	87	AliRICHDigit dd(1,1,1,pHit->LorsX(),pHit->LorsY()); //tmp digit to shift hit y to the nearest anod wire
	88	Float_t y= (pHit->LorsY() > dd.LorsY()) ? dd.LorsY()+0.21 : dd.LorsY()-0.21;
	89
	90	Int_t iSdiCnt=pSdiLst->GetEntries();
	91	for(Int_t i=0;i<9;i++){ //affected pads loop
	92	AliRICHDigit dig(pHit->Ch(),qdcTot,pHit->Tid(),pHit->LorsX(),y,i); //c,q,tid,x,y create tmp sdigit for pad i around hit position
	93	if(dig.PadX()==-1) continue;
	94	new((*pSdiLst)[iSdiCnt++]) AliRICHDigit(dig);
	95	}
	96	return qdcTot;
	97	}//Hit2Sdi()
9de8d990	98	//++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
0fe8fa07	99	void AliRICHDigit::Print(Option_t*)const
0fe8fa07	100	{
e30ca504	101	// Print current digit
	102	// Arguments: option string not used
	103	// Returns: none
	104	Printf("pad=(%2i,%2i,%3i,%3i),pos=(%7.2f,%7.2f) QDC=%8.3f, TID=(%5i,%5i,%5i) raw r=%2i d=%2i a=%2i",
	105	Ch(),Pc(),PadX(),PadY(),LorsX(),LorsY(), Q(), fTracks[0],fTracks[1],fTracks[2], Row(), Dilogic(), Addr());
0fe8fa07	106	}
9de8d990	107	//++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
e30ca504	108	void AliRICHDigit::PrintSize()
0422a446	109	{
e30ca504	110	// Print all segmentaion related sizes
	111	// Arguments: none
	112	// Returns: none
	113	Printf("-->pad =(%6.2f,%6.2f) cm dead zone %.2f cm\n"
	114	"-->PC =(%6.2f,%6.2f) cm\n"
	115	"-->all PCs=(%6.2f,%6.2f) cm",
	116	SizePadX(),SizePadY(),SizeDead(),SizePcX(),SizePcY(),SizeAllX(),SizeAllY());
0422a446	117	}
9de8d990	118	//++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
e30ca504	119	void AliRICHDigit::TestSeg()
	120	{
	121	// Draws the picture of segmentation
	122	// Arguments: none
	123	// Returns: none
	124	TCanvas *pC=new TCanvas("pads","View from electronics side, IP is behind the picture.");pC->ToggleEventStatus();
	125	gPad->AddExec("test","AliRICHDigit::Zoom()");
	126	DrawPc();
	127	}//TestSeg()
	128	//++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
	129	void AliRICHDigit::Zoom()
	130	{
	131	// Show info about current cursur position in status bar of the canvas
	132	// Arguments: none
	133	// Returns: none
	134	TCanvas pC=(TCanvas)gPad;
	135	TRootCanvas pRC= (TRootCanvas)pC->GetCanvasImp();
	136	TGStatusBar *pBar=pRC->GetStatusBar();
	137	pBar->SetParts(5);
	138	Float_t x=gPad->AbsPixeltoX(gPad->GetEventX());
	139	Float_t y=gPad->AbsPixeltoY(gPad->GetEventY());
	140	AliRICHDigit dig(1,100,1,x,y);
	141	if(IsInDead(x,y))
	142	pBar->SetText("Out of sensitive area",4);
	143	else
	144	pBar->SetText(Form("p%i x%i y%i r%i d%i a%i (%.2f,%.2f)",dig.Pc(),dig.PadX(),dig.PadY(),dig.Row(),dig.Dilogic(),dig.Addr(),dig.LorsX(),dig.LorsY()),4);
	145	if(gPad->GetEvent()==1){
	146	new TCanvas("zoom",Form("Row %i DILOGIC %i",dig.Row(),dig.Dilogic()));
	147	}
	148	}//Zoom()
	149	//++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
	150	void AliRICHDigit::DrawPc(Bool_t isFill)
	151	{
	152	// Utility methode draws RICH chamber PCs on event display.
	153	// Arguments: none
	154	// Returns: none
	155	// y6 ---------- ----------
	156	// \| \| \| \|
	157	// \| 4 \| \| 5 \|
	158	// y5 ---------- ----------
	159	//
	160	// y4 ---------- ----------
	161	// \| \| \| \|
	162	// \| 2 \| \| 3 \| view from electronics side
	163	// y3 ---------- ----------
	164	//
	165	// y2 ---------- ----------
	166	// \| \| \| \|
	167	// \| 0 \| \| 1 \|
	168	// y1 ---------- ----------
	169	// x1 x2 x3 x4
	170	gPad->Range(-5,-5,SizeAllX()+5,SizeAllY()+5);
	171	Float_t x1=0,x2=SizePcX(),x3=SizePcX()+SizeDead(), x4=SizeAllX();
	172	Float_t y1=0,y2=SizePcY(),y3=SizePcY()+SizeDead(),y4=2*SizePcY()+SizeDead(),y5=SizeAllY()-SizePcY(),y6=SizeAllY();
	173
	174	Float_t xL[5]={x1,x1,x2,x2,x1}; //clockwise
	175	Float_t xR[5]={x3,x3,x4,x4,x3};
	176	Float_t yD[5]={y1,y2,y2,y1,y1};
	177	Float_t yC[5]={y3,y4,y4,y3,y3};
	178	Float_t yU[5]={y5,y6,y6,y5,y5};
	179
	180	TLatex t; t.SetTextSize(0.01); t.SetTextAlign(22);
	181	Int_t iColLeft=29,iColRight=41;
	182	TPolyLine pc=0; TLine pL; Float_t x0=0,y0=0,wRow=kDilYSizePadY(),wDil=kDilXSizePadX();
183	for(Int_t iPc=0;iPc<kPcAll;iPc++){
184	if(iPc==4) {pc=new TPolyLine(5,xL,yU);t.DrawText(x1-3.5,y6-20,"PC4");x0=x1;y0=y5;} if(iPc==5) {pc=new TPolyLine(5,xR,yU);t.DrawText(x4+3.5,y6-20,"PC5");x0=x3;y0=y5;}
185	if(iPc==2) {pc=new TPolyLine(5,xL,yC);t.DrawText(x1-3.5,y4-20,"PC2");x0=x1;y0=y3;} if(iPc==3) {pc=new TPolyLine(5,xR,yC);t.DrawText(x4+3.5,y4-20,"PC3");x0=x3;y0=y3;}
186	if(iPc==0) {pc=new TPolyLine(5,xL,yD);t.DrawText(x1-3.5,y2-20,"PC0");x0=x1;y0=y1;} if(iPc==1) {pc=new TPolyLine(5,xR,yD);t.DrawText(x4+3.5,y2-20,"PC1");x0=x3;y0=y1;}
187	(iPc%2)? pc->SetFillColor(iColLeft): pc->SetFillColor(iColRight);
188	if(isFill) pc->Draw("f"); else pc->Draw();
189	for(Int_t i=1;i<=8 ;i++){//draw row lines (horizontal)
190	Float_t y=y0+i*wRow;
191	Int_t row=i+iPc/2*8; if(iPc%2!=0) row=25-row; t.DrawText(x0-1,y -3,Form("r%i",row));
192	if(i==8) break; //do not draw the last line of PC
193	pL=new TLine(x0,y,x0+SizePcX(),y); pL->Draw();
194	}
195	for(Int_t iDil=1;iDil<=10;iDil++){Float_t x=x0+iDil*wDil;t.DrawText(x -3,y0-1,Form("d%i",11-iDil)); if(iDil==10) break; pL=new TLine(x,y0,x,y0+SizePcY()); pL->Draw();}
196	}
197	}//DrawPc()
198	//++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++