[u/mrichter/AliRoot.git] / RICH / AliRICHDisplFast.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 "AliRICHDisplFast.h"
#include "AliRICH.h"
#include "AliRICHChamber.h"
#include "AliRICHParam.h"
#include <AliLoader.h>
#include <TCanvas.h>
#include <TPolyLine.h>
#include <TParticle.h>
#include <TStyle.h>
#include <TMath.h>
#include <TLatex.h>

ClassImp(AliRICHDisplFast)
//__________________________________________________________________________________________________
void AliRICHDisplFast::ShowEvent(Int_t iEvtNmin,Int_t iEvtNmax)
{
  TH2F *pDigitsH2[8];
  char titobj[11],titdisp[20];

  AliRICH *pRich = (AliRICH*)gAlice->GetDetector("RICH");
  Bool_t isDigits  =!pRich->GetLoader()->LoadDigits();
  if(!isDigits){Error("ShoEvent","No digits. Nothing to display.");return;}
  
  TCanvas *canvas = new TCanvas("RICHDisplay","RICH Display",0,0,1226,900);
//  TCanvas *canvas = (TCanvas*)gROOT->FindObjectAny("RICHDisplay");
  
  gStyle->SetPalette(1);

  canvas->Divide(3,3);
  
  for(Int_t iChamber=1;iChamber<=7;iChamber++) {
    sprintf(titobj,"pDigitsH2_%i",iChamber);
    sprintf(titdisp,"Chamber  %i",iChamber);
    pDigitsH2[iChamber] = new TH2F(titobj,titdisp,165,0,AliRICHParam::PcSizeX(),144,0,AliRICHParam::PcSizeY());
    pDigitsH2[iChamber]->SetMarkerColor(kGreen); 
    pDigitsH2[iChamber]->SetMarkerStyle(29); 
    pDigitsH2[iChamber]->SetMarkerSize(0.4);
    pDigitsH2[iChamber]->SetStats(kFALSE);
  }
  
  if(iEvtNmax>gAlice->GetEventsPerRun()) iEvtNmax=gAlice->GetEventsPerRun();

  TText *t = new TText();
  t->SetTextSize(0.10);

  TText *tit = new TText(0.1,0.6,"RICH Display");
  tit->SetTextSize(0.10);

  for(Int_t iEventN=iEvtNmin;iEventN<=iEvtNmax;iEventN++) {
        
    canvas->cd(1);
    sprintf(titdisp,"Event Number %i",iEventN);
    t->SetText(0.2,0.4,titdisp);
    t->Draw();
    tit->Draw();
    pRich->GetLoader()->GetRunLoader()->GetEvent(iEventN);
    pRich->GetLoader()->TreeD()->GetEntry(0);
    for(Int_t iChamber=1;iChamber<=7;iChamber++) {
      pDigitsH2[iChamber]->Reset();    
      for(Int_t j=0;j<pRich->Digits(iChamber)->GetEntries();j++) {//digits loop
        AliRICHdigit *pDig = (AliRICHdigit*)pRich->Digits(iChamber)->At(j);
        TVector2 x2=AliRICHParam::Pad2Loc(pDig->Pad());
        pDigitsH2[iChamber]->Fill(x2.X(),x2.Y());
      }
      if(iChamber==1) canvas->cd(7);
      if(iChamber==2) canvas->cd(8);
      if(iChamber==3) canvas->cd(4);
      if(iChamber==4) canvas->cd(5);
      if(iChamber==5) canvas->cd(6);
      if(iChamber==6) canvas->cd(2);
      if(iChamber==7) canvas->cd(3);
      pDigitsH2[iChamber]->Draw();
    }
    canvas->Update();
    canvas->Modified();
    
    if(iEvtNmin<iEvtNmax) gPad->WaitPrimitive();
//    canvas->cd(3);
//    gPad->Clear();
  }
}
//__________________________________________________________________________________________________
void AliRICHDisplFast::Exec(Option_t *)
{
  AliRICH *pRich = (AliRICH*)gAlice->GetDetector("RICH");
  Bool_t isHits    =!pRich->GetLoader()->LoadHits();
  Bool_t isDigits  =!pRich->GetLoader()->LoadDigits();
  Bool_t isClusters=!pRich->GetLoader()->LoadRecPoints();
  
  if(!isHits && !isDigits && !isClusters){Error("Exec","No hits digits and clusters. Nothing to display.");return;}
  
  TCanvas *Display = new TCanvas("Display","RICH Display",0,0,600,600);
  gStyle->SetPalette(1);
  
  TH2F *pHitsH2=0,*pDigitsH2=0,*pClustersH2=0;
  
  if(isHits)     pHitsH2     = new TH2F("pHitsH2"  ,  "Event Display;x,cm;y,cm",165,0,AliRICHParam::PcSizeX(),
                                                                      144,0,AliRICHParam::PcSizeY());
  if(pHitsH2)    pHitsH2->SetStats(kFALSE);
  
  if(isDigits)   pDigitsH2   = new TH2F("pDigitsH2"  ,"Event Display",165,0,AliRICHParam::PcSizeX(),
                                                                      144,0,AliRICHParam::PcSizeY());
  if(isClusters) pClustersH2 = new TH2F("pClustersH2","Event Display",165,0,AliRICHParam::PcSizeX(),
                                                                      144,0,AliRICHParam::PcSizeY());
  
  for(Int_t iEventN=0;iEventN<gAlice->GetEventsPerRun();iEventN++){//events Loop
    pRich->GetLoader()->GetRunLoader()->GetEvent(iEventN);

  
    Int_t nPrimaries = (Int_t)pRich->GetLoader()->TreeH()->GetEntries();
    TObjArray * Hits = new TObjArray[nPrimaries];
  
    for(Int_t i=0;i<nPrimaries;i++) {
      pRich->GetLoader()->TreeH()->GetEntry(i);
      Int_t nHits = pRich->Hits()->GetEntries();
      for(Int_t k=0;k<nHits;k++)         Hits[i].Add(pRich->Hits()->At(k));
    
    }
//display all the staff on chamber by chamber basis           
    for(Int_t iChamber=1;iChamber<=7;iChamber++){//chambers loop       
      if(isHits)     pHitsH2    ->Reset();     
      if(isDigits)   pDigitsH2  ->Reset();     
      if(isClusters) pClustersH2->Reset();
//deals with hits
      for(Int_t i=0;i<nPrimaries;i++){//prims loop
        pRich->GetLoader()->TreeH()->GetEntry(i);
        Int_t nHits = pRich->Hits()->GetEntries();
        for(Int_t j=0;j<nHits;j++){//hits loop
          AliRICHhit *pHit = (AliRICHhit*)Hits[i].At(j);
          if(pHit->C()==iChamber){
            TVector3 hitGlobX3= pHit->OutX3();
            TVector2 hitLocX2 = pRich->C(iChamber)->Mrs2Pc(hitGlobX3);
            pHitsH2->Fill(hitLocX2.X(),hitLocX2.Y(),200);
          }//if
        }//hits loop         
      }//prims loop
      pHitsH2->SetTitle(Form("event %i chamber %2i",iEventN,iChamber));
      pHitsH2->SetMarkerColor(kRed); pHitsH2->SetMarkerStyle(29); pHitsH2->SetMarkerSize(0.4);
      pHitsH2->Draw();
      DrawSectors();
      TLatex l; l.SetNDC(); l.SetTextSize(0.02);
      if(!isHits)     {l.SetTextColor(kRed)  ;l.DrawLatex(0.1,0.01,"No Hits"    );}
      if(!isDigits)   {l.SetTextColor(kGreen);l.DrawLatex(0.4,0.01,"No DIGITS"  );}
      if(!isClusters) {l.SetTextColor(kBlue) ;l.DrawLatex(0.8,0.01,"No CLUSTERS");}
      Display->Update();
      Display->Modified();
      gPad->WaitPrimitive();
//deals with digits      
      if(isDigits){
        pRich->GetLoader()->TreeD()->GetEntry(0);
        for(Int_t j=0;j<pRich->Digits(iChamber)->GetEntries();j++){//digits loop
          AliRICHdigit *pDig = (AliRICHdigit*)pRich->Digits(iChamber)->At(j);
	  TVector2 x2=AliRICHParam::Pad2Loc(pDig->Pad());
	  pDigitsH2->Fill(x2.X(),x2.Y(),100);
        }//digits loop
        pDigitsH2->SetMarkerColor(kGreen); pDigitsH2->SetMarkerStyle(29); pDigitsH2->SetMarkerSize(0.4);
        pDigitsH2->Draw("same");
        Display->Update();
        Display->Modified();       
        gPad->WaitPrimitive();
      }//if(isDigits)      
//deals with clusters      
      if(isClusters){
        pRich->GetLoader()->TreeR()->GetEntry(0);
        for(Int_t j=0;j<pRich->Clusters(iChamber)->GetEntries();j++){//clusters loop
          AliRICHcluster *pClus = (AliRICHcluster*)pRich->Clusters(iChamber)->At(j);
          pClustersH2->Fill(pClus->X(),pClus->Y(),50);
        }//clusters loop
        pClustersH2->SetMarkerColor(kBlue); pClustersH2->SetMarkerStyle(29);  pClustersH2->SetMarkerSize(0.4);
        pClustersH2->Draw("same");
        Display->Update();
        Display->Modified();
        gPad->WaitPrimitive();
      }//if(isClusters)
    }//chambers loop
    delete [] Hits;
  }//events Loop
  pRich->GetLoader()->UnloadHits();
  if(isDigits)   pRich->GetLoader()->UnloadDigits();
  if(isClusters) pRich->GetLoader()->UnloadRecPoints();
}//Exec()
//__________________________________________________________________________________________________
void AliRICHDisplFast::DrawSectors() 
{ 
  AliRICHParam p;
  Double_t xLeft[5]  = {0,0,p.SectorSizeX(),p.SectorSizeX(),0};
  Double_t xRight[5] = {p.SectorSizeX()+p.DeadZone(),p.SectorSizeX()+p.DeadZone(),p.PcSizeX(),p.PcSizeX(),p.SectorSizeX()+p.DeadZone()};
  
  Double_t yDown[5]   = {0,p.SectorSizeY(),p.SectorSizeY(),0,0};
  Double_t yCenter[5] = {  p.SectorSizeY()+p.DeadZone(),2*p.SectorSizeY()+p.DeadZone(),2*p.SectorSizeY()+p.DeadZone(),
                           p.SectorSizeY()+p.DeadZone(),p.SectorSizeY()+p.DeadZone()};  
  Double_t yUp[5]     = {2*p.SectorSizeY()+2*p.DeadZone(),p.PcSizeY(),p.PcSizeY(),2*p.SectorSizeY()+2*p.DeadZone(),2*p.SectorSizeY()+2*p.DeadZone()};
  
  TPolyLine *sec1 = new TPolyLine(5,xLeft ,yDown);    sec1->SetLineColor(21);  sec1->Draw();
  TPolyLine *sec2 = new TPolyLine(5,xRight,yDown);    sec2->SetLineColor(21);  sec2->Draw();
  TPolyLine *sec3 = new TPolyLine(5,xLeft ,yCenter);  sec3->SetLineColor(21);  sec3->Draw();
  TPolyLine *sec4 = new TPolyLine(5,xRight,yCenter);  sec4->SetLineColor(21);  sec4->Draw();
  TPolyLine *sec5 = new TPolyLine(5,xLeft, yUp);      sec5->SetLineColor(21);  sec5->Draw();
  TPolyLine *sec6 = new TPolyLine(5,xRight,yUp);      sec6->SetLineColor(21);  sec6->Draw();
}//DrawSectors()
Commit	Line	Data
6b17b320	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	**************************************************************************/
6b17b320	15	#include "AliRICHDisplFast.h"
53fd478b	16	#include "AliRICH.h"
53fd478b	17	#include "AliRICHChamber.h"
6b17b320	18	#include "AliRICHParam.h"
6b17b320	19	#include <AliLoader.h>
6b17b320	20	#include <TCanvas.h>
af3d25a6	21	#include <TPolyLine.h>
6b17b320	22	#include <TParticle.h>
6b17b320	23	#include <TStyle.h>
6b17b320	24	#include <TMath.h>
d4c94996	25	#include <TLatex.h>
6b17b320	26
aed240d4	27	ClassImp(AliRICHDisplFast)
1cb5a97c	28	//__________________________________________________________________________________________________
	29	void AliRICHDisplFast::ShowEvent(Int_t iEvtNmin,Int_t iEvtNmax)
	30	{
	31	TH2F *pDigitsH2[8];
	32	char titobj[11],titdisp[20];
	33
	34	AliRICH pRich = (AliRICH)gAlice->GetDetector("RICH");
	35	Bool_t isDigits =!pRich->GetLoader()->LoadDigits();
	36	if(!isDigits){Error("ShoEvent","No digits. Nothing to display.");return;}
	37
	38	TCanvas *canvas = new TCanvas("RICHDisplay","RICH Display",0,0,1226,900);
	39	// TCanvas canvas = (TCanvas)gROOT->FindObjectAny("RICHDisplay");
	40
	41	gStyle->SetPalette(1);
6b17b320	42
1cb5a97c	43	canvas->Divide(3,3);
	44
	45	for(Int_t iChamber=1;iChamber<=7;iChamber++) {
	46	sprintf(titobj,"pDigitsH2_%i",iChamber);
	47	sprintf(titdisp,"Chamber %i",iChamber);
	48	pDigitsH2[iChamber] = new TH2F(titobj,titdisp,165,0,AliRICHParam::PcSizeX(),144,0,AliRICHParam::PcSizeY());
	49	pDigitsH2[iChamber]->SetMarkerColor(kGreen);
	50	pDigitsH2[iChamber]->SetMarkerStyle(29);
	51	pDigitsH2[iChamber]->SetMarkerSize(0.4);
	52	pDigitsH2[iChamber]->SetStats(kFALSE);
	53	}
	54
	55	if(iEvtNmax>gAlice->GetEventsPerRun()) iEvtNmax=gAlice->GetEventsPerRun();
	56
	57	TText *t = new TText();
	58	t->SetTextSize(0.10);
	59
	60	TText *tit = new TText(0.1,0.6,"RICH Display");
	61	tit->SetTextSize(0.10);
	62
	63	for(Int_t iEventN=iEvtNmin;iEventN<=iEvtNmax;iEventN++) {
	64
	65	canvas->cd(1);
	66	sprintf(titdisp,"Event Number %i",iEventN);
	67	t->SetText(0.2,0.4,titdisp);
	68	t->Draw();
	69	tit->Draw();
	70	pRich->GetLoader()->GetRunLoader()->GetEvent(iEventN);
	71	pRich->GetLoader()->TreeD()->GetEntry(0);
	72	for(Int_t iChamber=1;iChamber<=7;iChamber++) {
	73	pDigitsH2[iChamber]->Reset();
	74	for(Int_t j=0;j<pRich->Digits(iChamber)->GetEntries();j++) {//digits loop
	75	AliRICHdigit pDig = (AliRICHdigit)pRich->Digits(iChamber)->At(j);
	76	TVector2 x2=AliRICHParam::Pad2Loc(pDig->Pad());
	77	pDigitsH2[iChamber]->Fill(x2.X(),x2.Y());
	78	}
	79	if(iChamber==1) canvas->cd(7);
	80	if(iChamber==2) canvas->cd(8);
	81	if(iChamber==3) canvas->cd(4);
	82	if(iChamber==4) canvas->cd(5);
	83	if(iChamber==5) canvas->cd(6);
	84	if(iChamber==6) canvas->cd(2);
	85	if(iChamber==7) canvas->cd(3);
	86	pDigitsH2[iChamber]->Draw();
	87	}
	88	canvas->Update();
	89	canvas->Modified();
	90
	91	if(iEvtNmin<iEvtNmax) gPad->WaitPrimitive();
	92	// canvas->cd(3);
	93	// gPad->Clear();
	94	}
	95	}
aed240d4	96	//__________________________________________________________________________________________________
e42a7b46	97	void AliRICHDisplFast::Exec(Option_t *)
6b17b320	98	{
d4c94996	99	AliRICH pRich = (AliRICH)gAlice->GetDetector("RICH");
	100	Bool_t isHits =!pRich->GetLoader()->LoadHits();
	101	Bool_t isDigits =!pRich->GetLoader()->LoadDigits();
	102	Bool_t isClusters=!pRich->GetLoader()->LoadRecPoints();
	103
	104	if(!isHits && !isDigits && !isClusters){Error("Exec","No hits digits and clusters. Nothing to display.");return;}
	105
	106	TCanvas *Display = new TCanvas("Display","RICH Display",0,0,600,600);
	107	gStyle->SetPalette(1);
	108
	109	TH2F pHitsH2=0,pDigitsH2=0,*pClustersH2=0;
	110
e42a7b46	111	if(isHits) pHitsH2 = new TH2F("pHitsH2" , "Event Display;x,cm;y,cm",165,0,AliRICHParam::PcSizeX(),
e42a7b46	112	144,0,AliRICHParam::PcSizeY());
d4c94996	113	if(pHitsH2) pHitsH2->SetStats(kFALSE);
d4c94996	114
e42a7b46	115	if(isDigits) pDigitsH2 = new TH2F("pDigitsH2" ,"Event Display",165,0,AliRICHParam::PcSizeX(),
	116	144,0,AliRICHParam::PcSizeY());
	117	if(isClusters) pClustersH2 = new TH2F("pClustersH2","Event Display",165,0,AliRICHParam::PcSizeX(),
	118	144,0,AliRICHParam::PcSizeY());
d4c94996	119
aed240d4	120	for(Int_t iEventN=0;iEventN<gAlice->GetEventsPerRun();iEventN++){//events Loop
aed240d4	121	pRich->GetLoader()->GetRunLoader()->GetEvent(iEventN);
6b17b320	122
c94cc386	123
aed240d4	124	Int_t nPrimaries = (Int_t)pRich->GetLoader()->TreeH()->GetEntries();
aed240d4	125	TObjArray * Hits = new TObjArray[nPrimaries];
6b17b320	126
aed240d4	127	for(Int_t i=0;i<nPrimaries;i++) {
	128	pRich->GetLoader()->TreeH()->GetEntry(i);
	129	Int_t nHits = pRich->Hits()->GetEntries();
	130	for(Int_t k=0;k<nHits;k++) Hits[i].Add(pRich->Hits()->At(k));
	131
6b17b320	132	}
d4c94996	133	//display all the staff on chamber by chamber basis
	134	for(Int_t iChamber=1;iChamber<=7;iChamber++){//chambers loop
	135	if(isHits) pHitsH2 ->Reset();
	136	if(isDigits) pDigitsH2 ->Reset();
	137	if(isClusters) pClustersH2->Reset();
	138	//deals with hits
aed240d4	139	for(Int_t i=0;i<nPrimaries;i++){//prims loop
6b17b320	140	pRich->GetLoader()->TreeH()->GetEntry(i);
6b17b320	141	Int_t nHits = pRich->Hits()->GetEntries();
aed240d4	142	for(Int_t j=0;j<nHits;j++){//hits loop
6b17b320	143	AliRICHhit pHit = (AliRICHhit)Hits[i].At(j);
aed240d4	144	if(pHit->C()==iChamber){
d4c94996	145	TVector3 hitGlobX3= pHit->OutX3();
998b831f	146	TVector2 hitLocX2 = pRich->C(iChamber)->Mrs2Pc(hitGlobX3);
d4c94996	147	pHitsH2->Fill(hitLocX2.X(),hitLocX2.Y(),200);
aed240d4	148	}//if
	149	}//hits loop
	150	}//prims loop
d4c94996	151	pHitsH2->SetTitle(Form("event %i chamber %2i",iEventN,iChamber));
aed240d4	152	pHitsH2->SetMarkerColor(kRed); pHitsH2->SetMarkerStyle(29); pHitsH2->SetMarkerSize(0.4);
aed240d4	153	pHitsH2->Draw();
af3d25a6	154	DrawSectors();
d4c94996	155	TLatex l; l.SetNDC(); l.SetTextSize(0.02);
	156	if(!isHits) {l.SetTextColor(kRed) ;l.DrawLatex(0.1,0.01,"No Hits" );}
	157	if(!isDigits) {l.SetTextColor(kGreen);l.DrawLatex(0.4,0.01,"No DIGITS" );}
	158	if(!isClusters) {l.SetTextColor(kBlue) ;l.DrawLatex(0.8,0.01,"No CLUSTERS");}
6b17b320	159	Display->Update();
6b17b320	160	Display->Modified();
998b831f	161	gPad->WaitPrimitive();
d4c94996	162	//deals with digits
	163	if(isDigits){
	164	pRich->GetLoader()->TreeD()->GetEntry(0);
	165	for(Int_t j=0;j<pRich->Digits(iChamber)->GetEntries();j++){//digits loop
	166	AliRICHdigit pDig = (AliRICHdigit)pRich->Digits(iChamber)->At(j);
e42a7b46	167	TVector2 x2=AliRICHParam::Pad2Loc(pDig->Pad());
d4c94996	168	pDigitsH2->Fill(x2.X(),x2.Y(),100);
	169	}//digits loop
	170	pDigitsH2->SetMarkerColor(kGreen); pDigitsH2->SetMarkerStyle(29); pDigitsH2->SetMarkerSize(0.4);
	171	pDigitsH2->Draw("same");
	172	Display->Update();
	173	Display->Modified();
998b831f	174	gPad->WaitPrimitive();
d4c94996	175	}//if(isDigits)
	176	//deals with clusters
	177	if(isClusters){
	178	pRich->GetLoader()->TreeR()->GetEntry(0);
	179	for(Int_t j=0;j<pRich->Clusters(iChamber)->GetEntries();j++){//clusters loop
	180	AliRICHcluster pClus = (AliRICHcluster)pRich->Clusters(iChamber)->At(j);
	181	pClustersH2->Fill(pClus->X(),pClus->Y(),50);
	182	}//clusters loop
	183	pClustersH2->SetMarkerColor(kBlue); pClustersH2->SetMarkerStyle(29); pClustersH2->SetMarkerSize(0.4);
	184	pClustersH2->Draw("same");
	185	Display->Update();
	186	Display->Modified();
998b831f	187	gPad->WaitPrimitive();
d4c94996	188	}//if(isClusters)
d4c94996	189	}//chambers loop
aed240d4	190	delete [] Hits;
d4c94996	191	}//events Loop
aed240d4	192	pRich->GetLoader()->UnloadHits();
d4c94996	193	if(isDigits) pRich->GetLoader()->UnloadDigits();
d4c94996	194	if(isClusters) pRich->GetLoader()->UnloadRecPoints();
aed240d4	195	}//Exec()
aed240d4	196	//__________________________________________________________________________________________________
af3d25a6	197	void AliRICHDisplFast::DrawSectors()
af3d25a6	198	{
e42a7b46	199	AliRICHParam p;
	200	Double_t xLeft[5] = {0,0,p.SectorSizeX(),p.SectorSizeX(),0};
	201	Double_t xRight[5] = {p.SectorSizeX()+p.DeadZone(),p.SectorSizeX()+p.DeadZone(),p.PcSizeX(),p.PcSizeX(),p.SectorSizeX()+p.DeadZone()};
	202
	203	Double_t yDown[5] = {0,p.SectorSizeY(),p.SectorSizeY(),0,0};
	204	Double_t yCenter[5] = { p.SectorSizeY()+p.DeadZone(),2p.SectorSizeY()+p.DeadZone(),2p.SectorSizeY()+p.DeadZone(),
	205	p.SectorSizeY()+p.DeadZone(),p.SectorSizeY()+p.DeadZone()};
	206	Double_t yUp[5] = {2p.SectorSizeY()+2p.DeadZone(),p.PcSizeY(),p.PcSizeY(),2p.SectorSizeY()+2p.DeadZone(),2p.SectorSizeY()+2p.DeadZone()};
	207
	208	TPolyLine *sec1 = new TPolyLine(5,xLeft ,yDown); sec1->SetLineColor(21); sec1->Draw();
	209	TPolyLine *sec2 = new TPolyLine(5,xRight,yDown); sec2->SetLineColor(21); sec2->Draw();
	210	TPolyLine *sec3 = new TPolyLine(5,xLeft ,yCenter); sec3->SetLineColor(21); sec3->Draw();
	211	TPolyLine *sec4 = new TPolyLine(5,xRight,yCenter); sec4->SetLineColor(21); sec4->Draw();
	212	TPolyLine *sec5 = new TPolyLine(5,xLeft, yUp); sec5->SetLineColor(21); sec5->Draw();
	213	TPolyLine *sec6 = new TPolyLine(5,xRight,yUp); sec6->SetLineColor(21); sec6->Draw();
d4c94996	214	}//DrawSectors()