[u/mrichter/AliRoot.git] / ACORDE / AcoReco.C

//_____________________________________________________//
//                                                     //
//     This macro reads ACORDE Raw-Rootfied files      //
//     and draws Histograms of Acorde's Bitpattern     //    
//                                                     //
//     Author:                                         //
//                                                     //
//   Mario Rodriguez Cahuantzi <mrodrigu@mail.cern.ch> //
//                                                     //
//     Also comments to:                               //
//                                                     //
//   Mario Sitta <sitta@to.infn.it>                    //
//   Arturo Fernandez <afernan@mail.cern.ch>           //
//   Luciano Diaz <luciano.diaz@nucleares.unam.mx>     //
//   Eleazar Cuautle <ecuautle@nucleares.unam.mx>      //
//____________________________________________________ //


#include <TClonesArray.h>
#include <stdio.h>

//void AcoReco(char* fileName="08000020614001.20.root")
void AcoReco(char* fileName)
{
  TStopwatch timer;
  timer.Start();
  AliRawReader* rawReader = new AliRawReaderRoot(fileName); 
  AliRawReader* rCount = new AliRawReaderRoot(fileName);
  AliACORDERawStream* rawStream  = new AliACORDERawStream(rawReader);    
  TH1D *h1 = new TH1D("h1","ACORDE - Single Muon Hits",60,1,60);
  TH1D *h2 = new TH1D("h2","ACORDE - Hit Multiplicity",60,1,60);
  size_t contSingle=0;
  size_t contMulti=0;
  UInt_t dy[4];
  bool kroSingle[60],kroMulti[60];
  UInt_t tmpDy;
  Int_t dyma[60],DyM,nEvents=0;
  for(Int_t m=0;m<60;m++) {kroSingle[m]=0;kroMulti[m]=0;dyma[m]=0;}
  
 
  Int_t nEvents = rawStream->GetNEvents(fileName);
  
  printf("Numero de eventos: %d \n",nEvents);

for (Int_t i=1; i<=nEvents; i++) 
{
	if (!rawReader->NextEvent()) break;
	rawStream->Reset();
	if (!rawStream->Next()) break;
	dy[0]=rawStream->GetWord(0);
	dy[1]=rawStream->GetWord(1);
	dy[2]=rawStream->GetWord(2);
	dy[3]=rawStream->GetWord(3);
	tmpDy=dy[0];
	for(Int_t r=0;r<30;++r)
	{
		kroSingle[r] = tmpDy & 1;
		tmpDy>>=1;
	}
	tmpDy=dy[1];
	for(Int_t r=30;r<60;++r)
	{
		kroSingle[r] = tmpDy & 1;
		tmpDy>>=1;
	}
	tmpDy=dy[2];
	for(Int_t r=0;r<30;++r)
	{
		kroMulti[r] = tmpDy & 1;
		tmpDy>>=1;
	}
	tmpDy=dy[3];
	for(Int_t r=30;r<60;++r)
	{
		kroMulti[r] = tmpDy & 1;
		tmpDy>>=1;
	}
	contSingle=0;	
	for(Int_t r=0;r<60;++r) 
	{
		if(kroSingle[r]==1) 
		{
			dyma[r]=dyma[r]+1;
			h1->Fill(r+1);
			contSingle=contSingle+1;
		}
		
	}h2->Fill(contSingle);
}

TCanvas *acorde = new TCanvas("ACORDE","ACORDE-Hist-Real Data",1);
acorde->Divide(2,1);
acorde->cd(1);
h1->GetXaxis()->SetTitle("No. of Channel");
h1->GetYaxis()->SetTitle("No. of Hits");
h1->SetFillColor(kRed);
h1->Draw();
acorde->cd(2);
gPad->SetLogy();
h2->GetXaxis()->SetTitle("No. of Modules");
h2->GetYaxis()->SetTitle("Multiplicity");
h2->SetFillColor(kBlue);
h2->Draw();

delete rawReader;
delete rawStream;
timer.Stop();
timer.Print();
}
Commit	Line	Data
bd438c48	1	//_____________________________________________________//
	2	// //
	3	// This macro reads ACORDE Raw-Rootfied files //
	4	// and draws Histograms of Acorde's Bitpattern //
	5	// //
	6	// Author: //
	7	// //
	8	// Mario Rodriguez Cahuantzi <mrodrigu@mail.cern.ch> //
	9	// //
	10	// Also comments to: //
	11	// //
	12	// Mario Sitta <sitta@to.infn.it> //
	13	// Arturo Fernandez <afernan@mail.cern.ch> //
	14	// Luciano Diaz <luciano.diaz@nucleares.unam.mx> //
	15	// Eleazar Cuautle <ecuautle@nucleares.unam.mx> //
	16	//____________________________________________________ //
	17
	18
	19	#include <TClonesArray.h>
	20	#include <stdio.h>
	21
	22	//void AcoReco(char* fileName="08000020614001.20.root")
	23	void AcoReco(char* fileName)
	24	{
	25	TStopwatch timer;
	26	timer.Start();
	27	AliRawReader* rawReader = new AliRawReaderRoot(fileName);
	28	AliRawReader* rCount = new AliRawReaderRoot(fileName);
	29	AliACORDERawStream* rawStream = new AliACORDERawStream(rawReader);
	30	TH1D *h1 = new TH1D("h1","ACORDE - Single Muon Hits",60,1,60);
	31	TH1D *h2 = new TH1D("h2","ACORDE - Hit Multiplicity",60,1,60);
	32	size_t contSingle=0;
	33	size_t contMulti=0;
	34	UInt_t dy[4];
	35	bool kroSingle[60],kroMulti[60];
	36	UInt_t tmpDy;
	37	Int_t dyma[60],DyM,nEvents=0;
	38	for(Int_t m=0;m<60;m++) {kroSingle[m]=0;kroMulti[m]=0;dyma[m]=0;}
	39
	40
	41	Int_t nEvents = rawStream->GetNEvents(fileName);
	42
	43	printf("Numero de eventos: %d \n",nEvents);
	44
	45	for (Int_t i=1; i<=nEvents; i++)
	46	{
	47	if (!rawReader->NextEvent()) break;
	48	rawStream->Reset();
	49	if (!rawStream->Next()) break;
	50	dy[0]=rawStream->GetWord(0);
	51	dy[1]=rawStream->GetWord(1);
	52	dy[2]=rawStream->GetWord(2);
	53	dy[3]=rawStream->GetWord(3);
	54	tmpDy=dy[0];
	55	for(Int_t r=0;r<30;++r)
	56	{
	57	kroSingle[r] = tmpDy & 1;
	58	tmpDy>>=1;
	59	}
	60	tmpDy=dy[1];
	61	for(Int_t r=30;r<60;++r)
	62	{
	63	kroSingle[r] = tmpDy & 1;
	64	tmpDy>>=1;
65	}
66	tmpDy=dy[2];
67	for(Int_t r=0;r<30;++r)
68	{
69	kroMulti[r] = tmpDy & 1;
70	tmpDy>>=1;
71	}
72	tmpDy=dy[3];
73	for(Int_t r=30;r<60;++r)
74	{
75	kroMulti[r] = tmpDy & 1;
76	tmpDy>>=1;
77	}
78	contSingle=0;
79	for(Int_t r=0;r<60;++r)
80	{
81	if(kroSingle[r]==1)
82	{
83	dyma[r]=dyma[r]+1;
84	h1->Fill(r+1);
85	contSingle=contSingle+1;
86	}
87
88	}h2->Fill(contSingle);
89	}
90
91	TCanvas *acorde = new TCanvas("ACORDE","ACORDE-Hist-Real Data",1);
92	acorde->Divide(2,1);
93	acorde->cd(1);
94	h1->GetXaxis()->SetTitle("No. of Channel");
95	h1->GetYaxis()->SetTitle("No. of Hits");
96	h1->SetFillColor(kRed);
97	h1->Draw();
98	acorde->cd(2);
99	gPad->SetLogy();
100	h2->GetXaxis()->SetTitle("No. of Modules");
101	h2->GetYaxis()->SetTitle("Multiplicity");
102	h2->SetFillColor(kBlue);
103	h2->Draw();
104
105	delete rawReader;
106	delete rawStream;
107	timer.Stop();
108	timer.Print();
109	}