[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++) 
{

//       printf("No Event %d",i);	
       if (!rawReader->NextEvent()) break;
	rawStream->Reset();
      // printf("No Event %d",i);
        
	if (!rawStream->Next()) continue;
	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	{
b7f06c9d	47
	48	// printf("No Event %d",i);
	49	if (!rawReader->NextEvent()) break;
bd438c48	50	rawStream->Reset();
b7f06c9d	51	// printf("No Event %d",i);
	52
	53	if (!rawStream->Next()) continue;
bd438c48	54	dy[0]=rawStream->GetWord(0);
	55	dy[1]=rawStream->GetWord(1);
	56	dy[2]=rawStream->GetWord(2);
	57	dy[3]=rawStream->GetWord(3);
	58	tmpDy=dy[0];
	59	for(Int_t r=0;r<30;++r)
	60	{
	61	kroSingle[r] = tmpDy & 1;
	62	tmpDy>>=1;
	63	}
	64	tmpDy=dy[1];
	65	for(Int_t r=30;r<60;++r)
	66	{
	67	kroSingle[r] = tmpDy & 1;
	68	tmpDy>>=1;
	69	}
	70	tmpDy=dy[2];
	71	for(Int_t r=0;r<30;++r)
	72	{
	73	kroMulti[r] = tmpDy & 1;
	74	tmpDy>>=1;
	75	}
	76	tmpDy=dy[3];
	77	for(Int_t r=30;r<60;++r)
	78	{
	79	kroMulti[r] = tmpDy & 1;
	80	tmpDy>>=1;
	81	}
	82	contSingle=0;
	83	for(Int_t r=0;r<60;++r)
	84	{
	85	if(kroSingle[r]==1)
	86	{
	87	dyma[r]=dyma[r]+1;
	88	h1->Fill(r+1);
	89	contSingle=contSingle+1;
	90	}
	91
	92	}h2->Fill(contSingle);
	93	}
	94
	95	TCanvas *acorde = new TCanvas("ACORDE","ACORDE-Hist-Real Data",1);
	96	acorde->Divide(2,1);
	97	acorde->cd(1);
	98	h1->GetXaxis()->SetTitle("No. of Channel");
	99	h1->GetYaxis()->SetTitle("No. of Hits");
	100	h1->SetFillColor(kRed);
	101	h1->Draw();
	102	acorde->cd(2);
	103	gPad->SetLogy();
	104	h2->GetXaxis()->SetTitle("No. of Modules");
	105	h2->GetYaxis()->SetTitle("Multiplicity");
	106	h2->SetFillColor(kBlue);
	107	h2->Draw();
	108
	109	delete rawReader;
	110	delete rawStream;
	111	timer.Stop();
	112	timer.Print();
	113	}