[u/mrichter/AliRoot.git] / PMD / AliPMDRootDataRead.C

// To read PMD raw root data and fetch the adc value for each cell
void AliPMDRootDataRead(Char_t *file="rawfile.root",Int_t runNr = 0,
			Int_t NEVT = 10)
{
  TObjArray pmdddlcont;

  AliCDBManager *man = AliCDBManager::Instance();
  man->SetDefaultStorage("local://$ALICE_ROOT/OCDB");
  man->SetRun(runNr);

  gBenchmark->Start("");
  gStyle->SetOptStat(0);

  Int_t   xpad, ypad;
  Float_t xx, yy;

  AliPMDUtility cc;
  TH2F *h2 = new TH2F("h2","Y vs. X",200,-100.,100.,200,-100.,100.);

  TH1F *h1 = new TH1F("h1","",200,0.,200.);


  for(Int_t ievt=495; ievt < NEVT; ievt++)
    {
  
      //AliRawReaderRoot reader("09000095029017.10.root",ievt);
      AliRawReaderRoot reader(file,ievt);

      cout<<" Processing Event No  : "<<ievt<<endl;

      AliPMDRawStream stream(&reader);
      
      Int_t iddl = -1;
      
      while((iddl = stream.DdlData(&pmdddlcont)) >= 0)
	{

	  Int_t ientries = pmdddlcont.GetEntries();
	  
	  cout << "iddl = " << iddl << " ientries = " << ientries << endl;
	  
	  for (Int_t ient = 0; ient < ientries; ient++)
	    {
	      AliPMDddldata *pmdddl = (AliPMDddldata*)pmdddlcont.UncheckedAt(ient);
	  
	      Int_t det = pmdddl->GetDetector();
	      Int_t smn = pmdddl->GetSMN();
	      Int_t mcm = pmdddl->GetMCM();
	      Int_t pbus = pmdddl->GetPatchBusId();
	      //Int_t chno = pmdddl->GetChannel();
	      
	      Int_t row = pmdddl->GetRow();
	      Int_t col = pmdddl->GetColumn();
	      Int_t sig = pmdddl->GetSignal();

	      //if (iddl == 5)
	      //printf("%d %d %d %d %d %d\n",iddl,det,smn,pbus,row,col);
	      
	      if(mcm == 0) continue;
	      

	      if(smn <12)
		{
		  xpad = col;
		  ypad = row;
		}
	      else if(smn >=12 && smn < 24)
		{
		  xpad = row;
		  ypad = col;
		}
	      
	      if (det == 0)
		{
		  cc.RectGeomCellPos(smn,xpad,ypad,xx,yy);
		  h2->Fill(xx,yy,sig); 
		  h1->Fill(sig);
		}
	      
	    }
	  
      pmdddlcont.Delete();
	}
      
    }//event loop
  h2->Draw();
  cc.SetWriteModule(1);  // set here 0 not to print the module no
  cc.DrawPMDModule(0);   // 0 : for preshower, 1:for CPV plane
  gBenchmark->Show("");
}
Commit	Line	Data
	1	// To read PMD raw root data and fetch the adc value for each cell
	2	void AliPMDRootDataRead(Char_t *file="rawfile.root",Int_t runNr = 0,
	3	Int_t NEVT = 10)
	4	{
	5	TObjArray pmdddlcont;
	6
	7	AliCDBManager *man = AliCDBManager::Instance();
	8	man->SetDefaultStorage("local://$ALICE_ROOT/OCDB");
	9	man->SetRun(runNr);
	10
	11	gBenchmark->Start("");
	12	gStyle->SetOptStat(0);
	13
	14	Int_t xpad, ypad;
	15	Float_t xx, yy;
	16
	17	AliPMDUtility cc;
	18	TH2F *h2 = new TH2F("h2","Y vs. X",200,-100.,100.,200,-100.,100.);
	19
	20	TH1F *h1 = new TH1F("h1","",200,0.,200.);
	21
	22
	23	for(Int_t ievt=495; ievt < NEVT; ievt++)
	24	{
	25
	26	//AliRawReaderRoot reader("09000095029017.10.root",ievt);
	27	AliRawReaderRoot reader(file,ievt);
	28
	29	cout<<" Processing Event No : "<<ievt<<endl;
	30
	31	AliPMDRawStream stream(&reader);
	32
	33	Int_t iddl = -1;
	34
	35	while((iddl = stream.DdlData(&pmdddlcont)) >= 0)
	36	{
	37
	38	Int_t ientries = pmdddlcont.GetEntries();
	39
	40	cout << "iddl = " << iddl << " ientries = " << ientries << endl;
	41
	42	for (Int_t ient = 0; ient < ientries; ient++)
	43	{
	44	AliPMDddldata pmdddl = (AliPMDddldata)pmdddlcont.UncheckedAt(ient);
	45
	46	Int_t det = pmdddl->GetDetector();
	47	Int_t smn = pmdddl->GetSMN();
	48	Int_t mcm = pmdddl->GetMCM();
	49	Int_t pbus = pmdddl->GetPatchBusId();
	50	//Int_t chno = pmdddl->GetChannel();
	51
	52	Int_t row = pmdddl->GetRow();
	53	Int_t col = pmdddl->GetColumn();
	54	Int_t sig = pmdddl->GetSignal();
	55
	56	//if (iddl == 5)
	57	//printf("%d %d %d %d %d %d\n",iddl,det,smn,pbus,row,col);
	58
	59	if(mcm == 0) continue;
	60
	61
	62	if(smn <12)
	63	{
	64	xpad = col;
	65	ypad = row;
	66	}
	67	else if(smn >=12 && smn < 24)
	68	{
	69	xpad = row;
	70	ypad = col;
	71	}
	72
	73	if (det == 0)
	74	{
	75	cc.RectGeomCellPos(smn,xpad,ypad,xx,yy);
	76	h2->Fill(xx,yy,sig);
	77	h1->Fill(sig);
	78	}
	79
	80	}
	81
	82	pmdddlcont.Delete();
	83	}
	84
	85	}//event loop
	86	h2->Draw();
	87	cc.SetWriteModule(1); // set here 0 not to print the module no
	88	cc.DrawPMDModule(0); // 0 : for preshower, 1:for CPV plane
	89	gBenchmark->Show("");
	90	}