[u/mrichter/AliRoot.git] / HMPID / HplotDA.C

#if !defined(__CINT__) || defined(__MAKECINT__)
#include <Riostream.h>     

#include <TSystem.h>
#include <TCanvas.h>
#include <TH2F.h>
#include <TH1F.h>
#include <TH1I.h>

#include "AliHMPIDRawStream.h"
#endif


TH2F hgPedMapMean[6]; 
TH2F hgPedMapSigma[6];
TH1F hgPedMapMean1D[6]; 
TH1F hgPedMapSigma1D[6]; 
TH2F hgPedMapMeanSigma[6];
TH1I hgDdlErr;

//++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
void Convert(Int_t ddl,Int_t r,Int_t d,Int_t a,Int_t &ch, Int_t &pc, Int_t &px, Int_t &py)
{
  
  Int_t a2y[6]={3,2,4,1,5,0};//pady for a given address (for single DILOGIC chip)
  
  Int_t ch=ddl/2;
  Int_t tmp=(r-1)/8;              Int_t pc=(ddl%2)? 5-2*tmp:2*tmp; 
                                  Int_t px=(d-1)*8+a/6;
        tmp=(ddl%2)?(24-r):r-1;   Int_t py=6*(tmp%8)+a2y[a%6];
}

//++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++      
void ProcPed(Int_t nDDL)
{

    if(gSystem->IsFileInIncludePath(Form("./HmpidPedDdl%02i.txt",nDDL)))
          ifstream infile(Form("./HmpidPedDdl%02i.txt",nDDL));
    else return;
    
    for(Int_t ni=0;ni<6;ni++)
      {
        hgPedMapMean[ni]  = new TH2F(Form("hPedMapMean_DDL%d_PC_%d" ,nDDL,ni),         Form("Pedestal: Mean, DDL(0-13): %d PC(0-5): %d;padx;pady" ,nDDL,ni),80,0,80,48,0,48);
        hgPedMapSigma[ni] = new TH2F(Form("hPedMapSigma_DDL%d_PC_%d",nDDL,ni),         Form("Pedestal: Sigma, DDL(0-13): %d PC(0-5): %d;padx;pady",nDDL,ni),80,0,80,48,0,48);
        hgPedMapMean1D[ni]  = new TH1F(Form("hPedMapMean1D_DDL%d_PC_%d" ,nDDL,ni),     Form("Pedestal: Mean, DDL(0-13): %d PC(0-5): %d;pad;Mean" ,nDDL,ni),3841,-0.5,3840.5);
        hgPedMapSigma1D[ni]  = new TH1F(Form("hPedMapSigma1D_DDL%d_PC_%d" ,nDDL,ni),   Form("Pedestal: Sigma, DDL(0-13): %d PC(0-5): %d;pad;Sigma" ,nDDL,ni),3841,-0.5,3840.5);
        hgPedMapMeanSigma[ni] = new TH2F(Form("hPedMapMeanSigma_DDL%d_PC_%d" ,nDDL,ni),Form("Pedestal, DDL(0-13): %d PC(0-5): %d;Mean;Sigma"),300,0,300,50,0,5);
       }
      
    Int_t nSigCut,r,d,a,hard;  Float_t mean,sigma;
    Int_t ch=0,pc=0,px=0,py=0;
    Int_t cnt=0;
    Int_t runNumber=-99999;
    Char_t tName[10];
    Int_t  ldcId;
    Int_t  timeStamp;
    Int_t  nEv;
    Printf("Start reading DDL: %d ...",nDDL);  
    infile>>tName>>runNumber;
    infile>>tName>>ldcId;
    infile>>tName>>timeStamp;
    infile>>tName>>nEv;  
    infile>>tName>>nSigCut;
    Printf("RunNumber: %d, LdcId: %d TimeStamp: %d nEv: %d",runNumber,ldcId,timeStamp,nEv);
    while(!infile.eof()){
      infile>>dec>>r>>d>>a>>mean>>sigma>>hex>>hard;
      Convert(nDDL,r,d,a,ch,pc,px,py); 
      hgPedMapMean[pc].Fill(px,py,mean);  
      hgPedMapSigma[pc].Fill(px,py,sigma);
      hgPedMapMean1D[pc].Fill(cnt%3840,mean);
      hgPedMapSigma1D[pc].Fill(cnt%3840,sigma);
      hgPedMapMeanSigma[pc].Fill(mean,sigma);
      cnt++;
      
    }
  infile.close();
  Printf("Stop reading DDL: %d ...",nDDL);
      
  gStyle->SetPalette(1);    
  
  TCanvas *cped = new TCanvas("cped","cped",800,800);     cped->Divide(2,2);
  TCanvas *cped2 = new TCanvas("cped2","cped2",800,800);  cped2->Divide(2,2);
  
  for(Int_t npc=0;npc<6;npc++) {
    if(hgPedMapMean[npc].GetEntries()==0) continue;
    cped->cd(1); hgPedMapMean[npc].Draw("surf1");
    cped->cd(2); hgPedMapSigma[npc].Draw("surf1");
    
    cped->cd(3); hgPedMapMean[npc].SetStats(0);hgPedMapMean[npc].Draw("colz");
    cped->cd(4); hgPedMapSigma[npc].SetStats(0);hgPedMapSigma[npc].Draw("colz");
    cped->SaveAs(Form("PedMap1_DDL%d_PC_%d.eps",nDDL,npc));
    cped->SaveAs(Form("PedMap1_DDL%d_PC_%d.gif",nDDL,npc));
    //cped->SaveAs(Form("PedMap1_DDL%d_PC_%d.pdf",nDDL,npc));
        
    cped2->cd(1); hgPedMapMean1D[npc].Draw();
    cped2->cd(2); hgPedMapSigma1D[npc].Draw();
    cped2->cd(3); hgPedMapMeanSigma[npc].Draw("colz");
    cped2->SaveAs(Form("PedMap2_DDL%d_PC_%d.eps",nDDL,npc));           
    cped2->SaveAs(Form("PedMap2_DDL%d_PC_%d.gif",nDDL,npc));           
    //cped2->SaveAs(Form("PedMap2_DDL%d_PC_%d.pdf",nDDL,npc));           
    
  }
    
       
}//ProcPed()
//++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++      
void ProcErr(Int_t nDDL)
{
  if(gSystem->IsFileInIncludePath(Form("./HmpidErrorsDdl%02i.txt",nDDL)))
          ifstream infile(Form("./HmpidErrorsDdl%02i.txt",nDDL));
    else return;
    
    hgDdlErr = new TH1I(Form("hPedErr_DDL%d",nDDL),Form("DDL(%d) Decoding Errors",nDDL),AliHMPIDRawStream::kSumErr+1,0,AliHMPIDRawStream::kSumErr);
     for(Int_t ilabel=0; ilabel< AliHMPIDRawStream::kSumErr; ilabel++) {
      hgDdlErr->SetStats(0);
      hgDdlErr->GetXaxis()->CenterLabels(kTRUE);
      hgDdlErr->GetXaxis()->SetBinLabel((ilabel+1),Form("%i  %s",ilabel+1,AliHMPIDRawStream::GetErrName(ilabel)));
      }
    hgDdlErr->SetFillColor(5);

    
    Int_t runNumber=-99999;
    Char_t tName[10];
    Int_t  ldcId;
    Int_t  timeStamp;
    Int_t  nEv;Int_t rerr;
     Printf("Start reading Error File DDL: %d ...",nDDL);  
    infile>>tName>>runNumber;
    infile>>tName>>ldcId;
    infile>>tName>>timeStamp;
    infile>>tName>>nEv;  
    Printf("RunNumber: %d, LdcId: %d TimeStamp: %d nEv: %d",runNumber,ldcId,timeStamp,nEv);
    for(Int_t ierr=0;ierr<AliHMPIDRawStream::kSumErr;ierr++)
    {
       infile>>rerr;hgDdlErr->SetBinContent(ierr+1,rerr);
     }
    
    
    TCanvas *cerr = new TCanvas("cped","cped");
    cerr->cd();
    hgDdlErr->Draw();
    cerr->SaveAs(Form("PedError_DDL%d.eps",nDDL)); 
    cerr->SaveAs(Form("PedError_DDL%d.gif",nDDL)); 
      
}
//++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++      
void HplotDA()
{
  for(Int_t i=0;i<14;i++) 
  {
   ProcPed(i);  
   ProcErr(i);
  }
}
//++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
Commit	Line	Data
20071be0	1	#if !defined(__CINT__) \|\| defined(__MAKECINT__)
	2	#include <Riostream.h>
	3
	4	#include <TSystem.h>
	5	#include <TCanvas.h>
	6	#include <TH2F.h>
	7	#include <TH1F.h>
	8	#include <TH1I.h>
	9
	10	#include "AliHMPIDRawStream.h"
	11	#endif
	12
	13
	14
	15	TH2F hgPedMapMean[6];
	16	TH2F hgPedMapSigma[6];
	17	TH1F hgPedMapMean1D[6];
	18	TH1F hgPedMapSigma1D[6];
	19	TH2F hgPedMapMeanSigma[6];
	20	TH1I hgDdlErr;
	21
	22	//++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
	23	void Convert(Int_t ddl,Int_t r,Int_t d,Int_t a,Int_t &ch, Int_t &pc, Int_t &px, Int_t &py)
	24	{
	25
	26	Int_t a2y[6]={3,2,4,1,5,0};//pady for a given address (for single DILOGIC chip)
	27
	28	Int_t ch=ddl/2;
	29	Int_t tmp=(r-1)/8; Int_t pc=(ddl%2)? 5-2tmp:2tmp;
	30	Int_t px=(d-1)*8+a/6;
	31	tmp=(ddl%2)?(24-r):r-1; Int_t py=6*(tmp%8)+a2y[a%6];
	32	}
	33
	34	//++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
	35	void ProcPed(Int_t nDDL)
	36	{
	37
	38	if(gSystem->IsFileInIncludePath(Form("./HmpidPedDdl%02i.txt",nDDL)))
	39	ifstream infile(Form("./HmpidPedDdl%02i.txt",nDDL));
	40	else return;
	41
	42	for(Int_t ni=0;ni<6;ni++)
	43	{
	44	hgPedMapMean[ni] = new TH2F(Form("hPedMapMean_DDL%d_PC_%d" ,nDDL,ni), Form("Pedestal: Mean, DDL(0-13): %d PC(0-5): %d;padx;pady" ,nDDL,ni),80,0,80,48,0,48);
	45	hgPedMapSigma[ni] = new TH2F(Form("hPedMapSigma_DDL%d_PC_%d",nDDL,ni), Form("Pedestal: Sigma, DDL(0-13): %d PC(0-5): %d;padx;pady",nDDL,ni),80,0,80,48,0,48);
	46	hgPedMapMean1D[ni] = new TH1F(Form("hPedMapMean1D_DDL%d_PC_%d" ,nDDL,ni), Form("Pedestal: Mean, DDL(0-13): %d PC(0-5): %d;pad;Mean" ,nDDL,ni),3841,-0.5,3840.5);
	47	hgPedMapSigma1D[ni] = new TH1F(Form("hPedMapSigma1D_DDL%d_PC_%d" ,nDDL,ni), Form("Pedestal: Sigma, DDL(0-13): %d PC(0-5): %d;pad;Sigma" ,nDDL,ni),3841,-0.5,3840.5);
	48	hgPedMapMeanSigma[ni] = new TH2F(Form("hPedMapMeanSigma_DDL%d_PC_%d" ,nDDL,ni),Form("Pedestal, DDL(0-13): %d PC(0-5): %d;Mean;Sigma"),300,0,300,50,0,5);
	49	}
	50
	51	Int_t nSigCut,r,d,a,hard; Float_t mean,sigma;
	52	Int_t ch=0,pc=0,px=0,py=0;
	53	Int_t cnt=0;
	54	Int_t runNumber=-99999;
	55	Char_t tName[10];
	56	Int_t ldcId;
	57	Int_t timeStamp;
	58	Int_t nEv;
	59	Printf("Start reading DDL: %d ...",nDDL);
	60	infile>>tName>>runNumber;
	61	infile>>tName>>ldcId;
	62	infile>>tName>>timeStamp;
	63	infile>>tName>>nEv;
	64	infile>>tName>>nSigCut;
65	Printf("RunNumber: %d, LdcId: %d TimeStamp: %d nEv: %d",runNumber,ldcId,timeStamp,nEv);
66	while(!infile.eof()){
67	infile>>dec>>r>>d>>a>>mean>>sigma>>hex>>hard;
68	Convert(nDDL,r,d,a,ch,pc,px,py);
69	hgPedMapMean[pc].Fill(px,py,mean);
70	hgPedMapSigma[pc].Fill(px,py,sigma);
71	hgPedMapMean1D[pc].Fill(cnt%3840,mean);
72	hgPedMapSigma1D[pc].Fill(cnt%3840,sigma);
73	hgPedMapMeanSigma[pc].Fill(mean,sigma);
74	cnt++;
75
76	}
77	infile.close();
78	Printf("Stop reading DDL: %d ...",nDDL);
79
80	gStyle->SetPalette(1);
81
82	TCanvas *cped = new TCanvas("cped","cped",800,800); cped->Divide(2,2);
83	TCanvas *cped2 = new TCanvas("cped2","cped2",800,800); cped2->Divide(2,2);
84
85	for(Int_t npc=0;npc<6;npc++) {
86	if(hgPedMapMean[npc].GetEntries()==0) continue;
87	cped->cd(1); hgPedMapMean[npc].Draw("surf1");
88	cped->cd(2); hgPedMapSigma[npc].Draw("surf1");
89
90	cped->cd(3); hgPedMapMean[npc].SetStats(0);hgPedMapMean[npc].Draw("colz");
91	cped->cd(4); hgPedMapSigma[npc].SetStats(0);hgPedMapSigma[npc].Draw("colz");
92	cped->SaveAs(Form("PedMap1_DDL%d_PC_%d.eps",nDDL,npc));
93	cped->SaveAs(Form("PedMap1_DDL%d_PC_%d.gif",nDDL,npc));
94	//cped->SaveAs(Form("PedMap1_DDL%d_PC_%d.pdf",nDDL,npc));
95
96	cped2->cd(1); hgPedMapMean1D[npc].Draw();
97	cped2->cd(2); hgPedMapSigma1D[npc].Draw();
98	cped2->cd(3); hgPedMapMeanSigma[npc].Draw("colz");
99	cped2->SaveAs(Form("PedMap2_DDL%d_PC_%d.eps",nDDL,npc));
100	cped2->SaveAs(Form("PedMap2_DDL%d_PC_%d.gif",nDDL,npc));
101	//cped2->SaveAs(Form("PedMap2_DDL%d_PC_%d.pdf",nDDL,npc));
102
103	}
104
105
106	}//ProcPed()
107	//++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
108	void ProcErr(Int_t nDDL)
109	{
110	if(gSystem->IsFileInIncludePath(Form("./HmpidErrorsDdl%02i.txt",nDDL)))
111	ifstream infile(Form("./HmpidErrorsDdl%02i.txt",nDDL));
112	else return;
113
114	hgDdlErr = new TH1I(Form("hPedErr_DDL%d",nDDL),Form("DDL(%d) Decoding Errors",nDDL),AliHMPIDRawStream::kSumErr+1,0,AliHMPIDRawStream::kSumErr);
115	for(Int_t ilabel=0; ilabel< AliHMPIDRawStream::kSumErr; ilabel++) {
116	hgDdlErr->SetStats(0);
117	hgDdlErr->GetXaxis()->CenterLabels(kTRUE);
118	hgDdlErr->GetXaxis()->SetBinLabel((ilabel+1),Form("%i %s",ilabel+1,AliHMPIDRawStream::GetErrName(ilabel)));
119	}
120	hgDdlErr->SetFillColor(5);
121
122
123	Int_t runNumber=-99999;
124	Char_t tName[10];
125	Int_t ldcId;
126	Int_t timeStamp;
127	Int_t nEv;Int_t rerr;
128	Printf("Start reading Error File DDL: %d ...",nDDL);
129	infile>>tName>>runNumber;
130	infile>>tName>>ldcId;
131	infile>>tName>>timeStamp;
132	infile>>tName>>nEv;
133	Printf("RunNumber: %d, LdcId: %d TimeStamp: %d nEv: %d",runNumber,ldcId,timeStamp,nEv);
134	for(Int_t ierr=0;ierr<AliHMPIDRawStream::kSumErr;ierr++)
135	{
136	infile>>rerr;hgDdlErr->SetBinContent(ierr+1,rerr);
137	}
138
139
140	TCanvas *cerr = new TCanvas("cped","cped");
141	cerr->cd();
142	hgDdlErr->Draw();
143	cerr->SaveAs(Form("PedError_DDL%d.eps",nDDL));
144	cerr->SaveAs(Form("PedError_DDL%d.gif",nDDL));
145
146	}
147	//++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
148	void HplotDA()
149	{
150	for(Int_t i=0;i<14;i++)
151	{
152	ProcPed(i);
153	ProcErr(i);
154	}
155	}
156	//++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++