[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 "AliRICH.h"
#include "AliRICHDisplFast.h"
#include "AliRICHParam.h"
#include <AliLoader.h>
#include <Riostream.h>
#include <TCanvas.h>
#include <TParticle.h>
#include <TStyle.h>
#include <TF1.h>
#include <TH2.h>
#include <TMath.h>
#include <TRandom.h>
#include <TText.h>

#define NPointsOfRing 201

// Geometry of the RICH at Star...

static const Int_t nPadX      = AliRICHParam::NpadsY();
static const Int_t nPadY      = AliRICHParam::NpadsX();
static const Float_t PadSizeX = AliRICHParam::PadSizeY();
static const Float_t PadSizeY = AliRICHParam::PadSizeX();
static const Float_t spacer   = AliRICHParam::DeadZone();
static const Float_t degree = 180/3.1415926535;

static const Float_t pi = TMath::Pi();

static const Float_t RadiatorWidth = AliRICHParam::FreonThickness();
static const Float_t QuartzWidth   = AliRICHParam::QuartzThickness();
static const Float_t GapWidth      = AliRICHParam::RadiatorToPads();

static const Float_t Xmin = -AliRICHParam::PcSizeY()/2.;
static const Float_t Xmax =  AliRICHParam::PcSizeY()/2.;
static const Float_t Ymin = -AliRICHParam::PcSizeX()/2.;
static const Float_t Ymax =  AliRICHParam::PcSizeX()/2.;


AliRICHDisplFast::AliRICHDisplFast()
{
}

void AliRICHDisplFast::Display()
{

  TText *text;

  TH2F *h2_disp   = new TH2F("h2_disp"  ,"STAR-RICH Event Display",165,Xmin,Xmax,144,Ymin,Ymax);
  TH2F *gHitsH2   = new TH2F("gHitsH2"  ,"STAR-RICH Event Display",165,Xmin,Xmax,144,Ymin,Ymax);
  TH2F *h2_dispad = new TH2F("h2_dispad","STAR-RICH Event Display",161,0.,161.,145,0.,145.);
  TH2F *gClustersH2 = new TH2F("gClusH2"  ,"STAR-RICH Event Display",165,Xmin,Xmax,144,Ymin,Ymax);

  TCanvas *Display = new TCanvas("Display","Star Display",0,0,800,800);
    
  // Display event...

  gStyle->SetPalette(1);

  AliRICH *pRich = (AliRICH*)gAlice->GetDetector("RICH");
  pRich->GetLoader()->LoadRecPoints();
  pRich->GetLoader()->LoadDigits();
  pRich->GetLoader()->LoadHits();
  pRich->GetLoader()->TreeD()->GetEntry(0);
  pRich->GetLoader()->TreeR()->GetEntry(0);

  cout << " start to display..." << endl;     
  
  Int_t nPrimaries = (Int_t)pRich->GetLoader()->TreeH()->GetEntries();
  TObjArray Hits[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));
    }
  }
           
  for(Int_t iChamber=1;iChamber<=7;iChamber++) {
    
     Int_t nDigits = pRich->Digits(iChamber)->GetEntries();
     Int_t nClusters = pRich->Clusters(iChamber)->GetEntries();
   
     Display->ToggleEventStatus();
     Display->Modified();
     
     text = new TText(0,0,"");
     text->SetTextFont(61);
     text->SetTextSize(0.03);
     text->SetTextAlign(22);                                                       
      
//     Display->Resize();

     h2_disp->Reset();
     gHitsH2->Reset();
     gClustersH2->Reset();
//     h2_dispad->Reset();

     Double_t xpad,ypad;

      for(Int_t i=0;i<nPrimaries;i++) {
        pRich->GetLoader()->TreeH()->GetEntry(i);
        Int_t nHits = pRich->Hits()->GetEntries();

        for(Int_t j=0;j<nHits;j++) {
        
          AliRICHhit *pHit = (AliRICHhit*)Hits[i].At(j);
          cout << " chamber " << iChamber << " hit n. " << j << " ch " << pHit->Chamber() << endl;
          if(pHit->C()==iChamber) {
            TVector3 xyzhit(pHit->X(),pHit->Y(),pHit->Z());
            TVector3 hitlocal = pRich->C(iChamber)->Glob2Loc(xyzhit);
            gHitsH2->Fill(hitlocal.X(),hitlocal.Y());
//            pHit->Print("");
          }
         }         
       }
     
     for(Int_t j=0;j<nDigits;j++)
	{
          AliRICHdigit *pDigit = (AliRICHdigit*)pRich->Digits(iChamber)->At(j);
	  AliRICHParam::Pad2Loc(pDigit->X(),pDigit->Y(),xpad,ypad);
          Float_t charge = (Float_t)pDigit->Q();
	  h2_disp->Fill(xpad,ypad,charge);
          h2_dispad->Fill(pDigit->X(),pDigit->Y(),pDigit->Q());
          pDigit->Print("");
	}
        
     for(Int_t j=0;j<nClusters;j++)
        {
          AliRICHcluster *pCluster = (AliRICHcluster*)pRich->Clusters(iChamber)->At(j);
//          gClustersH2->Fill(pCluster->X(),pCluster->Y(),pCluster->Q());
          gClustersH2->Fill(pCluster->X(),pCluster->Y(),100);
        }

        
      cout << " -----------------" << endl;
      
             
      // end loop
                    

      h2_disp->SetMaximum(200);
      h2_disp->SetStats(0);
      h2_disp->Draw("colz");
      h2_disp->SetTitle(Form("Chamber %2i",iChamber));

/*      
      h2_dispad->SetTitle(Form("Chamber %2i",iChamber));
      h2_dispad->SetMaximum(200);
      h2_dispad->SetStats(0);
      h2_dispad->Draw("colz");
*/
      Display->Update();
      Display->Modified();
     
      getchar();
      gHitsH2->SetTitle(Form("Chamber %2i",iChamber));
      gHitsH2->SetMarkerColor(kBlue);
      gHitsH2->SetMarkerStyle(29);
      gHitsH2->SetMarkerSize(0.4);
      gHitsH2->Draw("same");
      Display->Update();
      Display->Modified();
       
      getchar();
      gClustersH2->SetTitle(Form("Chamber %2i",iChamber));
      gClustersH2->SetMarkerColor(kBlue);
      gClustersH2->SetMarkerStyle(29);
      gClustersH2->SetMarkerSize(0.4);
      gClustersH2->Draw("same");
      Display->Update();
      Display->Modified();
      getchar();
     }
}
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	**************************************************************************/
	15
	16	#include "AliRICH.h"
	17	#include "AliRICHDisplFast.h"
	18	#include "AliRICHParam.h"
	19	#include <AliLoader.h>
	20	#include <Riostream.h>
	21	#include <TCanvas.h>
	22	#include <TParticle.h>
	23	#include <TStyle.h>
	24	#include <TF1.h>
	25	#include <TH2.h>
	26	#include <TMath.h>
	27	#include <TRandom.h>
	28	#include <TText.h>
	29
	30	#define NPointsOfRing 201
	31
	32	// Geometry of the RICH at Star...
	33
	34	static const Int_t nPadX = AliRICHParam::NpadsY();
	35	static const Int_t nPadY = AliRICHParam::NpadsX();
	36	static const Float_t PadSizeX = AliRICHParam::PadSizeY();
	37	static const Float_t PadSizeY = AliRICHParam::PadSizeX();
	38	static const Float_t spacer = AliRICHParam::DeadZone();
	39	static const Float_t degree = 180/3.1415926535;
	40
	41	static const Float_t pi = TMath::Pi();
	42
	43	static const Float_t RadiatorWidth = AliRICHParam::FreonThickness();
	44	static const Float_t QuartzWidth = AliRICHParam::QuartzThickness();
	45	static const Float_t GapWidth = AliRICHParam::RadiatorToPads();
	46
	47	static const Float_t Xmin = -AliRICHParam::PcSizeY()/2.;
	48	static const Float_t Xmax = AliRICHParam::PcSizeY()/2.;
	49	static const Float_t Ymin = -AliRICHParam::PcSizeX()/2.;
	50	static const Float_t Ymax = AliRICHParam::PcSizeX()/2.;
	51
	52
	53	AliRICHDisplFast::AliRICHDisplFast()
	54	{
	55	}
	56
	57	void AliRICHDisplFast::Display()
	58	{
	59
	60	TText *text;
	61
80ffc56e	62	TH2F *h2_disp = new TH2F("h2_disp" ,"STAR-RICH Event Display",165,Xmin,Xmax,144,Ymin,Ymax);
80ffc56e	63	TH2F *gHitsH2 = new TH2F("gHitsH2" ,"STAR-RICH Event Display",165,Xmin,Xmax,144,Ymin,Ymax);
6b17b320	64	TH2F *h2_dispad = new TH2F("h2_dispad","STAR-RICH Event Display",161,0.,161.,145,0.,145.);
80ffc56e	65	TH2F *gClustersH2 = new TH2F("gClusH2" ,"STAR-RICH Event Display",165,Xmin,Xmax,144,Ymin,Ymax);
6b17b320	66
	67	TCanvas *Display = new TCanvas("Display","Star Display",0,0,800,800);
	68
	69	// Display event...
	70
c94cc386	71	gStyle->SetPalette(1);
c94cc386	72
6b17b320	73	AliRICH pRich = (AliRICH)gAlice->GetDetector("RICH");
c94cc386	74	pRich->GetLoader()->LoadRecPoints();
6b17b320	75	pRich->GetLoader()->LoadDigits();
6b17b320	76	pRich->GetLoader()->LoadHits();
c94cc386	77	pRich->GetLoader()->TreeD()->GetEntry(0);
c94cc386	78	pRich->GetLoader()->TreeR()->GetEntry(0);
6b17b320	79
c94cc386	80	cout << " start to display..." << endl;
c94cc386	81
6b17b320	82	Int_t nPrimaries = (Int_t)pRich->GetLoader()->TreeH()->GetEntries();
	83	TObjArray Hits[nPrimaries];
	84
	85	for(Int_t i=0;i<nPrimaries;i++) {
	86	pRich->GetLoader()->TreeH()->GetEntry(i);
	87	Int_t nHits = pRich->Hits()->GetEntries();
	88	for(Int_t k=0;k<nHits;k++) {
	89	Hits[i].Add(pRich->Hits()->At(k));
	90	}
	91	}
	92
	93	for(Int_t iChamber=1;iChamber<=7;iChamber++) {
	94
c94cc386	95	Int_t nDigits = pRich->Digits(iChamber)->GetEntries();
c94cc386	96	Int_t nClusters = pRich->Clusters(iChamber)->GetEntries();
6b17b320	97
	98	Display->ToggleEventStatus();
	99	Display->Modified();
	100
	101	text = new TText(0,0,"");
	102	text->SetTextFont(61);
	103	text->SetTextSize(0.03);
	104	text->SetTextAlign(22);
	105
	106	// Display->Resize();
	107
	108	h2_disp->Reset();
	109	gHitsH2->Reset();
80ffc56e	110	gClustersH2->Reset();
6b17b320	111	// h2_dispad->Reset();
6b17b320	112
c712cb2f	113	Double_t xpad,ypad;
6b17b320	114
6b17b320	115	for(Int_t i=0;i<nPrimaries;i++) {
	116	pRich->GetLoader()->TreeH()->GetEntry(i);
	117	Int_t nHits = pRich->Hits()->GetEntries();
	118
	119	for(Int_t j=0;j<nHits;j++) {
	120
	121	AliRICHhit pHit = (AliRICHhit)Hits[i].At(j);
	122	cout << " chamber " << iChamber << " hit n. " << j << " ch " << pHit->Chamber() << endl;
	123	if(pHit->C()==iChamber) {
	124	TVector3 xyzhit(pHit->X(),pHit->Y(),pHit->Z());
c712cb2f	125	TVector3 hitlocal = pRich->C(iChamber)->Glob2Loc(xyzhit);
6b17b320	126	gHitsH2->Fill(hitlocal.X(),hitlocal.Y());
c94cc386	127	// pHit->Print("");
6b17b320	128	}
	129	}
	130	}
c94cc386	131
	132	for(Int_t j=0;j<nDigits;j++)
	133	{
	134	AliRICHdigit pDigit = (AliRICHdigit)pRich->Digits(iChamber)->At(j);
	135	AliRICHParam::Pad2Loc(pDigit->X(),pDigit->Y(),xpad,ypad);
	136	Float_t charge = (Float_t)pDigit->Q();
	137	h2_disp->Fill(xpad,ypad,charge);
	138	h2_dispad->Fill(pDigit->X(),pDigit->Y(),pDigit->Q());
	139	pDigit->Print("");
	140	}
	141
	142	for(Int_t j=0;j<nClusters;j++)
	143	{
	144	AliRICHcluster pCluster = (AliRICHcluster)pRich->Clusters(iChamber)->At(j);
	145	// gClustersH2->Fill(pCluster->X(),pCluster->Y(),pCluster->Q());
	146	gClustersH2->Fill(pCluster->X(),pCluster->Y(),100);
	147	}
	148
	149
	150
	151
	152	cout << " -----------------" << endl;
	153
6b17b320	154
	155	// end loop
	156
	157
	158	h2_disp->SetMaximum(200);
	159	h2_disp->SetStats(0);
	160	h2_disp->Draw("colz");
	161	h2_disp->SetTitle(Form("Chamber %2i",iChamber));
	162
	163	/*
	164	h2_dispad->SetTitle(Form("Chamber %2i",iChamber));
	165	h2_dispad->SetMaximum(200);
	166	h2_dispad->SetStats(0);
	167	h2_dispad->Draw("colz");
	168	*/
	169	Display->Update();
	170	Display->Modified();
	171
	172	getchar();
	173	gHitsH2->SetTitle(Form("Chamber %2i",iChamber));
	174	gHitsH2->SetMarkerColor(kBlue);
	175	gHitsH2->SetMarkerStyle(29);
	176	gHitsH2->SetMarkerSize(0.4);
	177	gHitsH2->Draw("same");
6b17b320	178	Display->Update();
	179	Display->Modified();
	180
	181	getchar();
c94cc386	182	gClustersH2->SetTitle(Form("Chamber %2i",iChamber));
	183	gClustersH2->SetMarkerColor(kBlue);
	184	gClustersH2->SetMarkerStyle(29);
	185	gClustersH2->SetMarkerSize(0.4);
	186	gClustersH2->Draw("same");
	187	Display->Update();
	188	Display->Modified();
	189	getchar();
	190	}
6b17b320	191	}
6b17b320	192