[u/mrichter/AliRoot.git] / ITS / AnalyzeSDDGainAllMod.C

#if !defined(__CINT__) || defined(__MAKECINT__)
#include <TH2F.h>
#include <TCanvas.h>
#include <TStopwatch.h>
#include <TStyle.h>
#include <TLatex.h>
#include <TFile.h>
#include <TGrid.h>
#include "AliRawReader.h"
#include "AliRawReaderDate.h"
#include "AliRawReaderRoot.h"
#include "AliITSOnlineSDDBase.h"
#include "AliITSOnlineSDDCMN.h"
#include "AliITSOnlineSDDTP.h"
#include "AliITSRawStreamSDD.h"
#include "AliITSRawStreamSDDCompressed.h"
#endif

// Macro for the analysis of PULSER runs (equivalent to ITSSDDGAINda.cxx)
// Two functions named AnalyzeSDDGainAllModules: 
// The first is for analyzing a local raw data file and takes as agrument the file name.
// The second is for running on ALIEN
// All DDLs are analyzed, the argument nDDL selects the DDL to be plotted
// Origin: F. Prino (prino@to.infn.it)


void AnalyzeSDDGainAllMod(Char_t *datafil, 
			  Int_t adcfreq=20, 
			  Int_t nDDL=0, 
			  Int_t firstEv=18, 
			  Int_t lastEv=22, 
			  Float_t pascalDAC=100){

  const Int_t kTotDDL=24;
  const Int_t kModPerDDL=12;
  const Int_t kSides=2;
  Bool_t writtenoutput=kFALSE;


  TH2F** histo = new TH2F*[kTotDDL*kModPerDDL*kSides];
  AliITSOnlineSDDTP **anal=new AliITSOnlineSDDTP*[kTotDDL*kModPerDDL*kSides];
  Bool_t isFilled[kTotDDL*kModPerDDL*kSides];

  Char_t hisnam[20];
  for(Int_t iddl=0; iddl<kTotDDL;iddl++){
    for(Int_t imod=0; imod<kModPerDDL;imod++){
      for(Int_t isid=0;isid<kSides;isid++){
	Int_t index=kSides*(kModPerDDL*iddl+imod)+isid;
	anal[index]=new AliITSOnlineSDDTP(iddl,imod,isid,pascalDAC);
	if(adcfreq==40) anal[index]->SetLastGoodTB(254);
	else anal[index]->SetLastGoodTB(126);
	sprintf(hisnam,"h%02dc%02ds%d",iddl,imod,isid);
	histo[index]=new TH2F(hisnam,"",256,-0.5,255.5,256,-0.5,255.5);
	isFilled[index]=0;
      }
    }
  }

  TCanvas* c0 = new TCanvas("c0","Ev Display",900,900);
  gStyle->SetPalette(1);
  Char_t text[50];

  Int_t iev=firstEv;
  AliRawReader *rd; 
  if(strstr(datafil,".root")!=0){
    rd=new AliRawReaderRoot(datafil,iev);
  }else{
    rd=new AliRawReaderDate(datafil,iev);
  }
  TLatex *t0=new TLatex();
  t0->SetNDC();
  t0->SetTextSize(0.06);
  t0->SetTextColor(4);

  do{
    c0->Clear();
    c0->Divide(4,6,0.001,0.001);
    printf("Event # %d\n",iev);
    rd->Reset();
    for(Int_t iddl=0; iddl<kTotDDL;iddl++){
      for(Int_t imod=0; imod<kModPerDDL;imod++){
	for(Int_t isid=0;isid<kSides;isid++){
	  Int_t index=kSides*(kModPerDDL*iddl+imod)+isid;
	  histo[index]->Reset();
	}
      }
    }

    UChar_t cdhAttr=AliITSRawStreamSDD::ReadBlockAttributes(rd);
    UInt_t amSamplFreq=AliITSRawStreamSDD::ReadAMSamplFreqFromCDH(cdhAttr);
    AliITSRawStream* s=AliITSRawStreamSDD::CreateRawStreamSDD(rd,cdhAttr);
    if(!writtenoutput){
      printf("Use %s raw stream, sampling frequency %d MHz\n",s->ClassName(),amSamplFreq);
      writtenoutput=kTRUE;
    }
    while(s->Next()){
      Int_t iDDL=rd->GetDDLID();
      Int_t iCarlos=s->GetCarlosId();
      if(s->IsCompletedModule()) continue;
      if(s->IsCompletedDDL()) continue;
      if(iDDL>=0 && iDDL<kTotDDL){ 
	Int_t index=kSides*(kModPerDDL*iDDL+iCarlos)+s->GetChannel(); 
	histo[index]->Fill(s->GetCoord2(),s->GetCoord1(),s->GetSignal());
	isFilled[index]=1;
      }
    }
    delete s;
    iev++;
    for(Int_t iddl=0; iddl<kTotDDL;iddl++){
      for(Int_t imod=0; imod<kModPerDDL;imod++){
	for(Int_t isid=0;isid<kSides;isid++){
	  Int_t index=kSides*(kModPerDDL*iddl+imod)+isid;
	  anal[index]->AddEvent(histo[index]);
	  if(iddl==nDDL){
	    Int_t index2=kSides*imod+isid;
	    c0->cd(index2+1);
	    histo[index]->DrawCopy("colz");
	    sprintf(text,"DDL %d channel %d Side %d",nDDL,imod,isid);
	    t0->DrawLatex(0.15,0.92,text);
	    c0->Update();
	  }
	}
      }
    }
    printf(" --- OK\n");
  }while(rd->NextEvent()&&iev<=lastEv);

  TH1F *htotgain=new TH1F("htotgain","",100,0.2,4.2);
  TH1F *htotpeakpos=new TH1F("htotpeakpos","",256,-0.5,255.5);
  TH1F *hstatus=new TH1F("hstatus","",2,-0.5,1.5);

  TFile *outfil=new TFile("SDDgain-results.root","recreate");
  for(Int_t iddl=0; iddl<kTotDDL;iddl++){
    for(Int_t imod=0; imod<kModPerDDL;imod++){
      for(Int_t isid=0;isid<kSides;isid++){
	Int_t index=kSides*(kModPerDDL*iddl+imod)+isid;
	if(isFilled[index]){
	  anal[index]->ValidateAnodes();
	  anal[index]->WriteToASCII();
	  anal[index]->WriteToROOT(outfil);
	  for(Int_t ian=0; ian<256;ian++){
	    Float_t gain=anal[index]->GetChannelGain(ian);
	    Float_t ppos=anal[index]->GetTimeBinTPPeak(ian);
	    Int_t anstatus=anal[index]->IsAnodeGood(ian);
	    hstatus->Fill(anstatus);
	    htotgain->Fill(gain);
	    htotpeakpos->Fill(ppos);
	  }
	}
      }
    }
  }
  outfil->Close();

  // Draw Statistics of baselines and noise
  TCanvas *call=new TCanvas("call","General stats",700,700);
  call->Divide(2,2);
  call->cd(1);
  htotpeakpos->Draw();
  htotpeakpos->GetXaxis()->SetTitle("TP peak position (Time Bin)");
  htotpeakpos->GetXaxis()->SetTitleSize(0.07);
  htotpeakpos->GetXaxis()->SetTitleOffset(0.6);
  call->cd(2);
  htotgain->Draw();
  htotgain->GetXaxis()->SetTitle("Gain (ADC/DAC)");
  htotgain->GetXaxis()->SetTitleSize(0.07);
  htotgain->GetXaxis()->SetTitleOffset(0.6);
  call->cd(3);
  hstatus->Draw();
  hstatus->GetXaxis()->SetTitle("Anode Status (0=bad 1=good)");
  hstatus->GetXaxis()->SetTitleSize(0.07);
  hstatus->GetXaxis()->SetTitleOffset(0.6);
  call->Update();
  call->SaveAs("GenStatsPulser.gif");

  // Draw baselines and noisegain and TestPulse time bin for all modules

  TH1F** hgain = new TH1F*[kModPerDDL*kSides];
  TH1F** htptb = new TH1F*[kModPerDDL*kSides];

  TCanvas *c1=new TCanvas("c1","DDL: TP position",900,900);
  c1->SetBottomMargin(0.14);
  c1->Divide(4,6,0.001,0.001);
  TCanvas *c2=new TCanvas("c2","DDL: gain",900,900);
  c2->SetBottomMargin(0.14);
  c2->Divide(4,6,0.001,0.001);

  for(Int_t imod=0; imod<kModPerDDL;imod++){
    for(Int_t isid=0;isid<kSides;isid++){
      Int_t index1=kSides*(kModPerDDL*nDDL+imod)+isid;
      Int_t index2=kSides*imod+isid;
      sprintf(text,"DDL %d channel %d Side %d",nDDL,imod,isid);

      TLatex *t3=new TLatex(0.15,0.92,text);
      t3->SetNDC();
      t3->SetTextSize(0.06);
      t3->SetTextColor(4);
      sprintf(hisnam,"hgain%ds%d",imod,isid);
      hgain[index2]=new TH1F(hisnam,"",256,-0.5,255.5);
      sprintf(hisnam,"htptb%ds%d",imod,isid);
      htptb[index2]=new TH1F(hisnam,"",256,-0.5,255.5);
      for(Int_t ian=0;ian<256;ian++){
	hgain[index2]->SetBinContent(ian+1,anal[index1]->GetChannelGain(ian));
	htptb[index2]->SetBinContent(ian+1,anal[index1]->GetTimeBinTPPeak(ian));
      }

      c1->cd(index2+1);
      htptb[index2]->Draw();
    //    htptb[imod]->SetMinimum(0);
    //    htptb[imod]->SetMaximum(75);
      htptb[index2]->GetXaxis()->SetTitle("Anode");
      htptb[index2]->GetYaxis()->SetTitle("TP position (Time Bin)");
      htptb[index2]->GetXaxis()->SetTitleSize(0.07);
      htptb[index2]->GetYaxis()->SetTitleSize(0.07);
      htptb[index2]->GetXaxis()->SetTitleOffset(0.6);
      htptb[index2]->GetYaxis()->SetTitleOffset(0.7);
      t3->Draw();
      c1->Update();


      c2->cd(index2+1); 
      hgain[index2]->SetMinimum(0.);
      hgain[index2]->SetMaximum(4.);
      hgain[index2]->Draw();
      hgain[index2]->GetXaxis()->SetTitle("Anode");
      hgain[index2]->GetYaxis()->SetTitle("Gain");
      hgain[index2]->GetXaxis()->SetTitleSize(0.07);
      hgain[index2]->GetYaxis()->SetTitleSize(0.07);
      hgain[index2]->GetXaxis()->SetTitleOffset(0.6);
      hgain[index2]->GetYaxis()->SetTitleOffset(0.7);
      hgain[index2]->SetStats(0);
      t3->Draw();
      c2->Update();
    }
  }

  c1->SaveAs("TPtimebin.gif");
  c2->SaveAs("Gain.gif");

}

void AnalyzeSDDGainAllMod(Int_t nrun, Int_t n2, Int_t year=2009, Char_t* dir="LHC09b_SDD",
			  Int_t adcfreq=20, 
			  Int_t nDDL=0, 
			  Int_t firstEv=18, 
			  Int_t lastEv=22, 
			  Float_t pascalDAC=100){


  TGrid::Connect("alien:",0,0,"t");
  Char_t filnam[200];
  sprintf(filnam,"alien:///alice/data/%d/%s/%09d/raw/%02d%09d%03d.10.root",year,dir,nrun,year-2000,nrun,n2);
  printf("Open file %s\n",filnam);
  AnalyzeSDDGainAllMod(filnam,adcfreq,nDDL,firstEv,lastEv,pascalDAC);
}
Commit	Line	Data
91282711	1	#if !defined(__CINT__) \|\| defined(__MAKECINT__)
	2	#include <TH2F.h>
	3	#include <TCanvas.h>
	4	#include <TStopwatch.h>
	5	#include <TStyle.h>
	6	#include <TLatex.h>
	7	#include <TFile.h>
	8	#include <TGrid.h>
	9	#include "AliRawReader.h"
	10	#include "AliRawReaderDate.h"
	11	#include "AliRawReaderRoot.h"
	12	#include "AliITSOnlineSDDBase.h"
	13	#include "AliITSOnlineSDDCMN.h"
	14	#include "AliITSOnlineSDDTP.h"
	15	#include "AliITSRawStreamSDD.h"
a8cd2ab2	16	#include "AliITSRawStreamSDDCompressed.h"
91282711	17	#endif
	18
	19	// Macro for the analysis of PULSER runs (equivalent to ITSSDDGAINda.cxx)
	20	// Two functions named AnalyzeSDDGainAllModules:
	21	// The first is for analyzing a local raw data file and takes as agrument the file name.
	22	// The second is for running on ALIEN
	23	// All DDLs are analyzed, the argument nDDL selects the DDL to be plotted
	24	// Origin: F. Prino (prino@to.infn.it)
	25
	26
a8cd2ab2	27	void AnalyzeSDDGainAllMod(Char_t *datafil,
	28	Int_t adcfreq=20,
	29	Int_t nDDL=0,
e5b2f7f2	30	Int_t firstEv=18,
	31	Int_t lastEv=22,
	32	Float_t pascalDAC=100){
91282711	33
91282711	34	const Int_t kTotDDL=24;
	35	const Int_t kModPerDDL=12;
	36	const Int_t kSides=2;
e5b2f7f2	37	Bool_t writtenoutput=kFALSE;
91282711	38
	39
	40	TH2F** histo = new TH2F[kTotDDLkModPerDDL*kSides];
	41	AliITSOnlineSDDTP *anal=new AliITSOnlineSDDTP[kTotDDLkModPerDDLkSides];
	42	Bool_t isFilled[kTotDDLkModPerDDLkSides];
	43
	44	Char_t hisnam[20];
	45	for(Int_t iddl=0; iddl<kTotDDL;iddl++){
	46	for(Int_t imod=0; imod<kModPerDDL;imod++){
	47	for(Int_t isid=0;isid<kSides;isid++){
	48	Int_t index=kSides(kModPerDDLiddl+imod)+isid;
	49	anal[index]=new AliITSOnlineSDDTP(iddl,imod,isid,pascalDAC);
b255f56d	50	if(adcfreq==40) anal[index]->SetLastGoodTB(254);
b255f56d	51	else anal[index]->SetLastGoodTB(126);
91282711	52	sprintf(hisnam,"h%02dc%02ds%d",iddl,imod,isid);
	53	histo[index]=new TH2F(hisnam,"",256,-0.5,255.5,256,-0.5,255.5);
	54	isFilled[index]=0;
	55	}
	56	}
	57	}
	58
	59	TCanvas* c0 = new TCanvas("c0","Ev Display",900,900);
	60	gStyle->SetPalette(1);
	61	Char_t text[50];
	62
	63	Int_t iev=firstEv;
	64	AliRawReader *rd;
	65	if(strstr(datafil,".root")!=0){
	66	rd=new AliRawReaderRoot(datafil,iev);
	67	}else{
	68	rd=new AliRawReaderDate(datafil,iev);
	69	}
	70	TLatex *t0=new TLatex();
	71	t0->SetNDC();
	72	t0->SetTextSize(0.06);
	73	t0->SetTextColor(4);
	74
	75	do{
	76	c0->Clear();
	77	c0->Divide(4,6,0.001,0.001);
	78	printf("Event # %d\n",iev);
91282711	79	rd->Reset();
	80	for(Int_t iddl=0; iddl<kTotDDL;iddl++){
	81	for(Int_t imod=0; imod<kModPerDDL;imod++){
	82	for(Int_t isid=0;isid<kSides;isid++){
	83	Int_t index=kSides(kModPerDDLiddl+imod)+isid;
	84	histo[index]->Reset();
	85	}
	86	}
	87	}
e5b2f7f2	88
	89	UChar_t cdhAttr=AliITSRawStreamSDD::ReadBlockAttributes(rd);
	90	UInt_t amSamplFreq=AliITSRawStreamSDD::ReadAMSamplFreqFromCDH(cdhAttr);
	91	AliITSRawStream* s=AliITSRawStreamSDD::CreateRawStreamSDD(rd,cdhAttr);
	92	if(!writtenoutput){
	93	printf("Use %s raw stream, sampling frequency %d MHz\n",s->ClassName(),amSamplFreq);
	94	writtenoutput=kTRUE;
a8cd2ab2	95	}
a8cd2ab2	96	while(s->Next()){
91282711	97	Int_t iDDL=rd->GetDDLID();
a8cd2ab2	98	Int_t iCarlos=s->GetCarlosId();
	99	if(s->IsCompletedModule()) continue;
	100	if(s->IsCompletedDDL()) continue;
de075dae	101	if(iDDL>=0 && iDDL<kTotDDL){
a8cd2ab2	102	Int_t index=kSides(kModPerDDLiDDL+iCarlos)+s->GetChannel();
a8cd2ab2	103	histo[index]->Fill(s->GetCoord2(),s->GetCoord1(),s->GetSignal());
91282711	104	isFilled[index]=1;
	105	}
	106	}
a8cd2ab2	107	delete s;
91282711	108	iev++;
	109	for(Int_t iddl=0; iddl<kTotDDL;iddl++){
	110	for(Int_t imod=0; imod<kModPerDDL;imod++){
	111	for(Int_t isid=0;isid<kSides;isid++){
	112	Int_t index=kSides(kModPerDDLiddl+imod)+isid;
	113	anal[index]->AddEvent(histo[index]);
	114	if(iddl==nDDL){
	115	Int_t index2=kSides*imod+isid;
	116	c0->cd(index2+1);
	117	histo[index]->DrawCopy("colz");
	118	sprintf(text,"DDL %d channel %d Side %d",nDDL,imod,isid);
	119	t0->DrawLatex(0.15,0.92,text);
	120	c0->Update();
	121	}
	122	}
	123	}
	124	}
	125	printf(" --- OK\n");
	126	}while(rd->NextEvent()&&iev<=lastEv);
	127
de075dae	128	TH1F *htotgain=new TH1F("htotgain","",100,0.2,4.2);
91282711	129	TH1F *htotpeakpos=new TH1F("htotpeakpos","",256,-0.5,255.5);
	130	TH1F *hstatus=new TH1F("hstatus","",2,-0.5,1.5);
	131
	132	TFile *outfil=new TFile("SDDgain-results.root","recreate");
	133	for(Int_t iddl=0; iddl<kTotDDL;iddl++){
	134	for(Int_t imod=0; imod<kModPerDDL;imod++){
	135	for(Int_t isid=0;isid<kSides;isid++){
	136	Int_t index=kSides(kModPerDDLiddl+imod)+isid;
	137	if(isFilled[index]){
	138	anal[index]->ValidateAnodes();
	139	anal[index]->WriteToASCII();
	140	anal[index]->WriteToROOT(outfil);
	141	for(Int_t ian=0; ian<256;ian++){
	142	Float_t gain=anal[index]->GetChannelGain(ian);
	143	Float_t ppos=anal[index]->GetTimeBinTPPeak(ian);
	144	Int_t anstatus=anal[index]->IsAnodeGood(ian);
	145	hstatus->Fill(anstatus);
	146	htotgain->Fill(gain);
	147	htotpeakpos->Fill(ppos);
	148	}
	149	}
	150	}
	151	}
	152	}
	153	outfil->Close();
	154
	155	// Draw Statistics of baselines and noise
	156	TCanvas *call=new TCanvas("call","General stats",700,700);
	157	call->Divide(2,2);
	158	call->cd(1);
	159	htotpeakpos->Draw();
	160	htotpeakpos->GetXaxis()->SetTitle("TP peak position (Time Bin)");
	161	htotpeakpos->GetXaxis()->SetTitleSize(0.07);
	162	htotpeakpos->GetXaxis()->SetTitleOffset(0.6);
	163	call->cd(2);
	164	htotgain->Draw();
	165	htotgain->GetXaxis()->SetTitle("Gain (ADC/DAC)");
	166	htotgain->GetXaxis()->SetTitleSize(0.07);
	167	htotgain->GetXaxis()->SetTitleOffset(0.6);
	168	call->cd(3);
	169	hstatus->Draw();
	170	hstatus->GetXaxis()->SetTitle("Anode Status (0=bad 1=good)");
	171	hstatus->GetXaxis()->SetTitleSize(0.07);
	172	hstatus->GetXaxis()->SetTitleOffset(0.6);
	173	call->Update();
	174	call->SaveAs("GenStatsPulser.gif");
	175
	176	// Draw baselines and noisegain and TestPulse time bin for all modules
	177
	178	TH1F** hgain = new TH1F[kModPerDDLkSides];
	179	TH1F** htptb = new TH1F[kModPerDDLkSides];
	180
	181	TCanvas *c1=new TCanvas("c1","DDL: TP position",900,900);
	182	c1->SetBottomMargin(0.14);
	183	c1->Divide(4,6,0.001,0.001);
	184	TCanvas *c2=new TCanvas("c2","DDL: gain",900,900);
	185	c2->SetBottomMargin(0.14);
	186	c2->Divide(4,6,0.001,0.001);
	187
	188	for(Int_t imod=0; imod<kModPerDDL;imod++){
	189	for(Int_t isid=0;isid<kSides;isid++){
	190	Int_t index1=kSides(kModPerDDLnDDL+imod)+isid;
	191	Int_t index2=kSides*imod+isid;
	192	sprintf(text,"DDL %d channel %d Side %d",nDDL,imod,isid);
193
194	TLatex *t3=new TLatex(0.15,0.92,text);
195	t3->SetNDC();
196	t3->SetTextSize(0.06);
197	t3->SetTextColor(4);
198	sprintf(hisnam,"hgain%ds%d",imod,isid);
199	hgain[index2]=new TH1F(hisnam,"",256,-0.5,255.5);
200	sprintf(hisnam,"htptb%ds%d",imod,isid);
201	htptb[index2]=new TH1F(hisnam,"",256,-0.5,255.5);
202	for(Int_t ian=0;ian<256;ian++){
203	hgain[index2]->SetBinContent(ian+1,anal[index1]->GetChannelGain(ian));
204	htptb[index2]->SetBinContent(ian+1,anal[index1]->GetTimeBinTPPeak(ian));
205	}
206
207	c1->cd(index2+1);
208	htptb[index2]->Draw();
209	// htptb[imod]->SetMinimum(0);
210	// htptb[imod]->SetMaximum(75);
211	htptb[index2]->GetXaxis()->SetTitle("Anode");
212	htptb[index2]->GetYaxis()->SetTitle("TP position (Time Bin)");
213	htptb[index2]->GetXaxis()->SetTitleSize(0.07);
214	htptb[index2]->GetYaxis()->SetTitleSize(0.07);
215	htptb[index2]->GetXaxis()->SetTitleOffset(0.6);
216	htptb[index2]->GetYaxis()->SetTitleOffset(0.7);
217	t3->Draw();
218	c1->Update();
219
220
221	c2->cd(index2+1);
222	hgain[index2]->SetMinimum(0.);
223	hgain[index2]->SetMaximum(4.);
224	hgain[index2]->Draw();
225	hgain[index2]->GetXaxis()->SetTitle("Anode");
226	hgain[index2]->GetYaxis()->SetTitle("Gain");
227	hgain[index2]->GetXaxis()->SetTitleSize(0.07);
228	hgain[index2]->GetYaxis()->SetTitleSize(0.07);
229	hgain[index2]->GetXaxis()->SetTitleOffset(0.6);
230	hgain[index2]->GetYaxis()->SetTitleOffset(0.7);
231	hgain[index2]->SetStats(0);
232	t3->Draw();
233	c2->Update();
234	}
235	}
236
237	c1->SaveAs("TPtimebin.gif");
238	c2->SaveAs("Gain.gif");
239
240	}
241
e5b2f7f2	242	void AnalyzeSDDGainAllMod(Int_t nrun, Int_t n2, Int_t year=2009, Char_t* dir="LHC09b_SDD",
a8cd2ab2	243	Int_t adcfreq=20,
a8cd2ab2	244	Int_t nDDL=0,
e5b2f7f2	245	Int_t firstEv=18,
	246	Int_t lastEv=22,
	247	Float_t pascalDAC=100){
	248
a8cd2ab2	249
91282711	250	TGrid::Connect("alien:",0,0,"t");
91282711	251	Char_t filnam[200];
e5b2f7f2	252	sprintf(filnam,"alien:///alice/data/%d/%s/%09d/raw/%02d%09d%03d.10.root",year,dir,nrun,year-2000,nrun,n2);
91282711	253	printf("Open file %s\n",filnam);
e5b2f7f2	254	AnalyzeSDDGainAllMod(filnam,adcfreq,nDDL,firstEv,lastEv,pascalDAC);
91282711	255	}