[u/mrichter/AliRoot.git] / T0 / AliT0CalibViewer.C

/* $Id$ */


#if !defined(__CINT__)
#include "TControlBar.h"
#include "TString.h"
#include "TRandom.h"
#include "TH1F.h"
#include "TF1.h"
#include "TCanvas.h"

#include "AliRun.h"
#include "AliT0CalibData.h"
#include "AliT0AlignData.h"
#include "AliCDBMetaData.h"
#include "AliCDBId.h"
#include "AliCDBEntry.h"
#include "AliCDBManager.h"
#include "AliCDBStorage.h"
#endif
TFile *gFile=0;

void AliT0CalibViewer()
{
  TControlBar *menu = new TControlBar("vertical","T0 CalibViewer");
  menu->AddButton("Open file","OpenFile()",
		  "Open file");
  menu->AddButton("Draw CFD","DrawCFD()",
		  "Draw all CFD");
  menu->AddButton("Draw LED","DrawLED()",
		  "Draw all LED");
  menu->AddButton("Draw LED-CFD","DrawLEDminCFD()",
		  "Draw LED-CFD");
  menu->AddButton("Draw QTC","DrawQTC()",
                  "Draw all QTC");
  menu->AddButton("Draw CFDvsQTC","DrawCFDvsQTC()",
		  "Draw CFD vs QTC");
  menu->AddButton("Draw CFDvsLED","DrawCFDvsLED()",
		  "Draw CFD vs LED-CFD");
  menu->Show();
}
void OpenFile()
{

const char *ft[]={"T0 raw files","t0tree*.root","All files","*",0,0};
  TString dir(".");
  TGFileInfo fi; fi.fFileTypes=ft; fi.fIniDir=StrDup(dir);
  new TGFileDialog(gClient->GetRoot(), 0x0, kFDOpen, &fi);

  if(!fi.fFilename) return;
  if(gFile){ gFile->Close(); gFile=0;}
  
  gFile=TFile::Open(fi.fFilename); 

}

//------------------------------------------------------------------------
void DrawCFD()
{
  Int_t npeaks = 20;
  Int_t sigma=3.;
  Bool_t down=false;
  Int_t index[20];
   
  TCanvas *c1 = new TCanvas("c1", "c1",0,48,1280,951);
  c1->Divide(4,3);
  gStyle->SetOptFit(1111);
  //c1->Divide(2,2);
  Char_t buf1[10];
  for (Int_t i=0; i<12; i++)
    {
      c1->cd(i+1);
      sprintf(buf1,"T0_C_%i_CFD",i+1);
      TH1F *cfd = (TH1F*) gFile->Get(buf1);
      //     cfd->Draw();
      TSpectrum *s = new TSpectrum(2*npeaks,1);
      Int_t nfound = s->Search(cfd,sigma," ",0.05);
      cout<<"Found "<<nfound<<" peaks sigma "<<sigma<<endl;;
      if(nfound!=0){
	Float_t *xpeak = s->GetPositionX();
  	TMath::Sort(nfound, xpeak, index,down);
  	Float_t xp = xpeak[index[0]];
	Int_t xbin = cfd->GetXaxis()->FindBin(xp);
	Float_t yp = cfd->GetBinContent(xbin);
	cout<<"xbin = "<<xbin<<"\txpeak = "<<xpeak[1]<<"\typeak = "<<yp<<endl;
	Float_t hmax = xp+3*sigma;
	Float_t hmin = xp-3*sigma;
	cout<<hmin<< " "<<hmax<<endl;
	cfd->GetXaxis()->SetRange(hmin-20,hmax+20);
	//	cout<<" cfd range "<<mean<<" rms "<<rms<<" "<<hmin<<" "<<hmax<<endl;
	TF1 *g1 = new TF1("g1", "gaus", hmin, hmax);
	//	cfd->Fit("g1","R");
	cfd->Draw();
      }
    }
 
}
//------------------------------------------------------------------------
void DrawLED()
{
  Int_t npeaks = 10;
  Int_t sigma=10.;
  Bool_t down=false;
  Int_t index[20];
  
  TCanvas *c1 = new TCanvas("c1", "c1",0,48,1280,951);
  c1->Divide(4,3);
  
  Char_t buf1[20];
  for (Int_t i=0; i<12; i++)
    {
      c1->cd(i+1);
      sprintf(buf1,"T0_C_%i_LED",i+1);
      TH1F *cfd = (TH1F*) gFile->Get(buf1);
      cfd->Draw();
      TSpectrum *s = new TSpectrum(2*npeaks,1);
      Int_t nfound = s->Search(cfd,sigma," ",0.2);
      cout<<"Found "<<nfound<<" peaks sigma "<<sigma<<endl;;
      if(nfound!=0) {
	Float_t *xpeak = s->GetPositionX();
	TMath::Sort(nfound, xpeak, index,down);
	Float_t xp = xpeak[index[0]];
        Int_t xbin = cfd->GetXaxis()->FindBin(xp);
        Float_t yp = cfd->GetBinContent(xbin);
        cout<<"xbin = "<<xbin<<"\txpeak = "<<xpeak[index[0]]<<"\typeak = "<<yp<<endl;
        Float_t hmin=xp-3*sigma;
        Float_t hmax =xp+3*sigma;
        cfd->GetXaxis()->SetRange(hmin,hmax);
 	TF1 *g1 = new TF1("g1", "gaus", hmin, hmax);
	cfd->Fit("g1","R");
	//      cfd->Draw();
	
      }
      
    }
}

//------------------------------------------------------------------------
void DrawQTC()
{
  
  TCanvas *c1 = new TCanvas("c1", "QTC",0,48,1280,951);
  // c1->Divide(4,3);
  c1->Divide(2,3);
  
  Char_t buf1[10];
  for (Int_t i=0; i<12; i++)
    {
      c1->cd(i+1);
      sprintf(buf1,"QTC%i",i+1);
      TH1F *qtc = (TH1F*) gFile->Get(buf1);
      Float_t mean = qtc->GetMean();
      Float_t rms = qtc->GetRMS();
      Float_t hminR=mean - 0.1*mean;
      Float_t hmaxR =mean + 0.1*mean;
      qtc->GetXaxis()->SetRange(hminR,hmaxR);
      Float_t hmin=mean - 3*rms;
      Float_t hmax =mean + 3*rms;
      qtc->GetXaxis()->SetRange(hmin,hmax);
      // TF1 *g2 = new TF1("g2", "gaus", hmin, hmax);
       //             qtc->Fit("g2","RQ");
       qtc->Draw();
    }
  
  
}

//------------------------------------------------------------------------
void DrawLEDminCFD()
{
  Int_t npeaks = 10;
  Int_t sigma=10.;
  Bool_t down=false;
  Int_t index[20];
  
  TCanvas *c1 = new TCanvas("c1", "c1",0,48,1280,951);
  c1->Divide(4,3);
  
  Char_t buf1[20];
  for (Int_t i=0; i<12; i++)
    {
      c1->cd(i+1);
      sprintf(buf1,"LED-CFD%i",i+1);
      TH1F *cfd = (TH1F*) gFile->Get(buf1);
      cfd->Draw();
      TSpectrum *s = new TSpectrum(2*npeaks,1);
      Int_t nfound = s->Search(cfd,sigma," ",0.2);
      cout<<"Found "<<nfound<<" peaks sigma "<<sigma<<endl;;
      if(nfound!=0) {
	Float_t *xpeak = s->GetPositionX();
	TMath::Sort(nfound, xpeak, index,down);
	Float_t xp = xpeak[index[0]];
        Int_t xbin = cfd->GetXaxis()->FindBin(xp);
        Float_t yp = cfd->GetBinContent(xbin);
        cout<<"xbin = "<<xbin<<"\txpeak = "<<xpeak[index[0]]<<"\typeak = "<<yp<<endl;
        Float_t hmin=xp-3*sigma;
        Float_t hmax =xp+3*sigma;
        cfd->GetXaxis()->SetRange(hmin,hmax);
	TF1 *g1 = new TF1("g1", "gaus", hmin, hmax);
        cfd->Fit("g1","RQ");

	
      }
      
    }
  /*
  TCanvas *c1 = new TCanvas("c1", "c1",0,48,1280,951);
  c1->Divide(2,2);
  Char_t buf1[10];
  for (Int_t i=0; i<4; i++)
    {
      c1->cd(i+1);
      sprintf(buf1,"LED-CFD%i",i+1);
      TH1F *cfd = (TH1F*) file->Get(buf1);
      //  cout<<buf1<<" "<<cfd<<endl;
      //     cfd->Draw();
      //   cfd->GetXaxis()->SetRange(0,100);
      Float_t mean = cfd->GetMean();
      Float_t rms = cfd->GetRMS();
      Float_t hmin=mean - 3*rms;
      Float_t hmax =mean + 3*rms;
      cfd->GetXaxis()->SetRange(hmin-10,hmax+10);
      cout<<" cfd range "<<mean<<" rms "<<rms<<" "<<hmin<<" "<<hmax<<endl;
      //     TF1 *g1 = new TF1("g1", "gaus", hmin, hmax);
      //  cfd->Fit("g1","RQ");
      cfd->Draw();
    }
  */

}
//------------------------------------------------------------------------
void DrawCFDvsQTC()
{
  Int_t npeaks = 20;
  Int_t sigma=3.;
  Bool_t down=false;

  Int_t index[20];
  Char_t buf1[10], buf2[10];
  
  TCanvas *c1 = new TCanvas("c1", "c1",0,48,1280,951);
  gStyle->SetOptStat(0);
  c1->Divide(3,2);
  
  for (Int_t i=0; i<5; i++)
    {
      c1->cd(i+1);
      sprintf(buf1,"T0_C_%i_CFD",i+1);
      sprintf(buf2,"CFDvsQTC%i",i+1);
      cout<<buf1<<" "<<buf2<<endl;
      TH2F *qtc_cfd = (TH2F*) gFile->Get(buf2);
      TH1F *cfd = (TH1F*) gFile->Get(buf1);
      //       cfd->Draw();
      TSpectrum *s = new TSpectrum(2*npeaks,1);
      Int_t nfound = s->Search(cfd,sigma," ",0.05);
      cout<<"Found "<<nfound<<" peaks sigma "<<sigma<<endl;;
      if(nfound!=0){
	Float_t *xpeak = s->GetPositionX();
  	TMath::Sort(nfound, xpeak, index,down);
  	Float_t xp = xpeak[index[0]];
	Int_t xbin = cfd->GetXaxis()->FindBin(xp);
	Float_t yp = cfd->GetBinContent(xbin);
	cout<<"xbin = "<<xbin<<"\txpeak = "<<xpeak[1]<<"\typeak = "<<yp<<endl;
	Float_t hmax = xp+10*sigma;
	Float_t hmin = xp-10*sigma;
	cout<<hmin<< " "<<hmax<<endl;
	//	    cfd->GetXaxis()->SetRange(hmin,hmax);
	//	    TF1 *g1 = new TF1("g1", "gaus", hmin, hmax);
	// cfd->Fit("g1","R");
	Int_t hminbin=  qtc_cfd->GetXaxis()->GetFirst();
	Int_t hmaxbin=  qtc_cfd->GetXaxis()->GetLast();
	Int_t nbins= qtc_cfd->GetXaxis()->GetNbins();
	cout<<"  qtc_cfd "<<hminbin<<" "<<hmaxbin<<" "<<nbins<<endl;
	//  qtc_cfd->Draw();
	TProfile *pr_y = qtc_cfd->ProfileX();
	
	pr_y->SetMaximum(hmax);
	pr_y->SetMinimum(hmin);
	Int_t np=nbins/20;
	Double_t *xx = new Double_t[np];
	Double_t *yy = new Double_t[np];
	Int_t ng=0;
	Double_t yg=0;
	for (Int_t ip=1; ip<nbins; ip++)
	  {
	    if(ip%20 != 0 ) {
	      if (pr_y->GetBinContent(ip) !=0)
		yg +=pr_y->GetBinContent(ip);
	      //	cout<<ng<<" "<<pr_y->GetBinContent(ip)<<" "<<yg<<endl;
	      ng++;}
	    else {
	      xx[ip/20] = Float_t (ip);
	      yy[ip/20] = yg/ng;
	      yg=0;
	      ng=0;
	      //	      cout<<ip<<" "<<ip/20<<" "<< xx[ip/20]<<" "<< yy[ip/20]<<endl;
	    }
	  }
	TH2F *hr = new TH2F("hr"," ",np,0,nbins, np, hmin, hmax);
	hr->Draw();
	TGraph *gr = new TGraph(np,xx,yy);
	gr->SetMinimum(hmin);
	gr->SetMaximum(hmax);
	gr->SetMarkerStyle(20);
	gr->Draw("P");
	//	delete [] xx;
	//	delete [] yy;
	
	// pr_y->Rebin(10);
		   //  pr_y->Draw();
	    
      }
      
    }
  

}
//------------------------------------------------------------------------
void DrawCFDvsLED()
{
  Int_t runNumber=1098;
  
  Int_t npeaks = 20;
  Int_t sigma=2.;
  Char_t buf1[10], buf2[10];
  
  for (Int_t i=0; i<1; i++)
    {
      //   c1->cd(i+1);
      sprintf(buf1,"T0_C_%i_CFD",i+1);
      sprintf(buf2,"CFD_LED%i",i+1);
      
      TH2F *qtc_cfd = (TH2F*) gFile->Get(buf2);
      TH1F *cfd = (TH1F*) gFile->Get(buf1);
      //       cfd->Draw();
      TSpectrum *s = new TSpectrum(2*npeaks,1);
      Int_t nfound = s->Search(cfd,sigma," ",0.05);
      cout<<"Found "<<nfound<<" peaks sigma "<<sigma<<endl;;
      if(nfound!=0){
	Double_t max=0.0; 
	Double_t tabmax[2] = {0.0, 0.0};
	Float_t *xpeak = s->GetPositionX();
	for(Int_t k=0; k<1 ;k++)
	  {
	    Float_t xp = xpeak[k];
	    Int_t xbin = cfd->GetXaxis()->FindBin(xp);
	    Float_t yp = cfd->GetBinContent(xbin);
	    cout<<"xbin = "<<xbin<<"\txpeak = "<<xpeak[k]<<"\typeak = "<<yp<<endl;
	  }
	Float_t hmin=xp-10*sigma;
	Float_t hmax =xp+10*sigma;
	cout<<hmin<< " "<<hmax<<endl;
	//	    cfd->GetXaxis()->SetRange(hmin,hmax);
	//	    TF1 *g1 = new TF1("g1", "gaus", hmin, hmax);
	// cfd->Fit("g1","R");
	Int_t hminbin=  qtc_cfd->GetXaxis()->GetFirst();
	Int_t hmaxbin=  qtc_cfd->GetXaxis()->GetLast();
	Int_t nbins= qtc_cfd->GetXaxis()->GetNbins();
	cout<<"  qtc_cfd "<<hminbin<<" "<<hmaxbin<<" "<<nbins<<endl;
	//  qtc_cfd->Draw();
	TProfile *pr_y = qtc_cfd->ProfileX();
	
	pr_y->SetMaximum(hmax);
	pr_y->SetMinimum(hmin);
	Int_t np=nbins/20;
	Double_t *xx = new Double_t[np];
	Double_t *yy = new Double_t[np];
	Int_t ng=0;
	Double_t yg=0;
	for (Int_t ip=1; ip<nbins; ip++)
	  {
	    if(ip%20 != 0 ) {
	      if (pr_y->GetBinContent(ip) !=0)
		yg +=pr_y->GetBinContent(ip);
	      //	cout<<ng<<" "<<pr_y->GetBinContent(ip)<<" "<<yg<<endl;
	      ng++;}
	    else {
	      xx[ip/20] = Float_t (ip);
	      yy[ip/20] = yg/ng;
	      yg=0;
	      ng=0;
	      cout<<ip<<" "<<ip/20<<" "<< xx[ip/20]<<" "<< yy[ip/20]<<endl;
	    }
	  }
	TH2F *hr = new TH2F("hr"," ",np,0,nbins, np, hmin, hmax);
	hr->Draw();
	TGraph *gr = new TGraph(np,xx,yy);
	gr->SetMinimum(hmin);
	gr->SetMaximum(hmax);
	gr->Draw("P");
	delete [] xx;
	delete [] yy;
	
	// pr_y->Rebin(10);
	//  pr_y->Draw();
	
      }
      
    }
  

}
Commit	Line	Data
79de2f45	1	/* $Id$ */
	2
	3
	4	#if !defined(__CINT__)
	5	#include "TControlBar.h"
	6	#include "TString.h"
	7	#include "TRandom.h"
	8	#include "TH1F.h"
	9	#include "TF1.h"
	10	#include "TCanvas.h"
	11
	12	#include "AliRun.h"
	13	#include "AliT0CalibData.h"
	14	#include "AliT0AlignData.h"
	15	#include "AliCDBMetaData.h"
	16	#include "AliCDBId.h"
	17	#include "AliCDBEntry.h"
	18	#include "AliCDBManager.h"
	19	#include "AliCDBStorage.h"
	20	#endif
	21	TFile *gFile=0;
	22
	23	void AliT0CalibViewer()
	24	{
	25	TControlBar *menu = new TControlBar("vertical","T0 CalibViewer");
	26	menu->AddButton("Open file","OpenFile()",
	27	"Open file");
	28	menu->AddButton("Draw CFD","DrawCFD()",
	29	"Draw all CFD");
	30	menu->AddButton("Draw LED","DrawLED()",
	31	"Draw all LED");
	32	menu->AddButton("Draw LED-CFD","DrawLEDminCFD()",
	33	"Draw LED-CFD");
	34	menu->AddButton("Draw QTC","DrawQTC()",
	35	"Draw all QTC");
	36	menu->AddButton("Draw CFDvsQTC","DrawCFDvsQTC()",
	37	"Draw CFD vs QTC");
	38	menu->AddButton("Draw CFDvsLED","DrawCFDvsLED()",
	39	"Draw CFD vs LED-CFD");
	40	menu->Show();
	41	}
	42	void OpenFile()
	43	{
	44
	45	const char ft[]={"T0 raw files","t0tree.root","All files","*",0,0};
	46	TString dir(".");
	47	TGFileInfo fi; fi.fFileTypes=ft; fi.fIniDir=StrDup(dir);
	48	new TGFileDialog(gClient->GetRoot(), 0x0, kFDOpen, &fi);
	49
	50	if(!fi.fFilename) return;
	51	if(gFile){ gFile->Close(); gFile=0;}
	52
	53	gFile=TFile::Open(fi.fFilename);
	54
	55	}
	56
	57	//------------------------------------------------------------------------
	58	void DrawCFD()
	59	{
	60	Int_t npeaks = 20;
	61	Int_t sigma=3.;
	62	Bool_t down=false;
	63	Int_t index[20];
	64
65	TCanvas *c1 = new TCanvas("c1", "c1",0,48,1280,951);
66	c1->Divide(4,3);
67	gStyle->SetOptFit(1111);
68	//c1->Divide(2,2);
69	Char_t buf1[10];
70	for (Int_t i=0; i<12; i++)
71	{
72	c1->cd(i+1);
73	sprintf(buf1,"T0_C_%i_CFD",i+1);
74	TH1F cfd = (TH1F) gFile->Get(buf1);
75	// cfd->Draw();
76	TSpectrum s = new TSpectrum(2npeaks,1);
77	Int_t nfound = s->Search(cfd,sigma," ",0.05);
78	cout<<"Found "<<nfound<<" peaks sigma "<<sigma<<endl;;
79	if(nfound!=0){
80	Float_t *xpeak = s->GetPositionX();
81	TMath::Sort(nfound, xpeak, index,down);
82	Float_t xp = xpeak[index[0]];
83	Int_t xbin = cfd->GetXaxis()->FindBin(xp);
84	Float_t yp = cfd->GetBinContent(xbin);
85	cout<<"xbin = "<<xbin<<"\txpeak = "<<xpeak[1]<<"\typeak = "<<yp<<endl;
86	Float_t hmax = xp+3*sigma;
87	Float_t hmin = xp-3*sigma;
88	cout<<hmin<< " "<<hmax<<endl;
89	cfd->GetXaxis()->SetRange(hmin-20,hmax+20);
90	// cout<<" cfd range "<<mean<<" rms "<<rms<<" "<<hmin<<" "<<hmax<<endl;
91	TF1 *g1 = new TF1("g1", "gaus", hmin, hmax);
92	// cfd->Fit("g1","R");
93	cfd->Draw();
94	}
95	}
96
97	}
98	//------------------------------------------------------------------------
99	void DrawLED()
100	{
101	Int_t npeaks = 10;
102	Int_t sigma=10.;
103	Bool_t down=false;
104	Int_t index[20];
105
106	TCanvas *c1 = new TCanvas("c1", "c1",0,48,1280,951);
107	c1->Divide(4,3);
108
109	Char_t buf1[20];
110	for (Int_t i=0; i<12; i++)
111	{
112	c1->cd(i+1);
113	sprintf(buf1,"T0_C_%i_LED",i+1);
114	TH1F cfd = (TH1F) gFile->Get(buf1);
115	cfd->Draw();
116	TSpectrum s = new TSpectrum(2npeaks,1);
117	Int_t nfound = s->Search(cfd,sigma," ",0.2);
118	cout<<"Found "<<nfound<<" peaks sigma "<<sigma<<endl;;
119	if(nfound!=0) {
120	Float_t *xpeak = s->GetPositionX();
121	TMath::Sort(nfound, xpeak, index,down);
122	Float_t xp = xpeak[index[0]];
123	Int_t xbin = cfd->GetXaxis()->FindBin(xp);
124	Float_t yp = cfd->GetBinContent(xbin);
125	cout<<"xbin = "<<xbin<<"\txpeak = "<<xpeak[index[0]]<<"\typeak = "<<yp<<endl;
126	Float_t hmin=xp-3*sigma;
127	Float_t hmax =xp+3*sigma;
128	cfd->GetXaxis()->SetRange(hmin,hmax);
129	TF1 *g1 = new TF1("g1", "gaus", hmin, hmax);
130	cfd->Fit("g1","R");
131	// cfd->Draw();
132
133	}
134
135	}
136	}
137
138	//------------------------------------------------------------------------
139	void DrawQTC()
140	{
141
142	TCanvas *c1 = new TCanvas("c1", "QTC",0,48,1280,951);
143	// c1->Divide(4,3);
144	c1->Divide(2,3);
145
146	Char_t buf1[10];
147	for (Int_t i=0; i<12; i++)
148	{
149	c1->cd(i+1);
150	sprintf(buf1,"QTC%i",i+1);
151	TH1F qtc = (TH1F) gFile->Get(buf1);
152	Float_t mean = qtc->GetMean();
153	Float_t rms = qtc->GetRMS();
154	Float_t hminR=mean - 0.1*mean;
155	Float_t hmaxR =mean + 0.1*mean;
156	qtc->GetXaxis()->SetRange(hminR,hmaxR);
157	Float_t hmin=mean - 3*rms;
158	Float_t hmax =mean + 3*rms;
159	qtc->GetXaxis()->SetRange(hmin,hmax);
160	// TF1 *g2 = new TF1("g2", "gaus", hmin, hmax);
161	// qtc->Fit("g2","RQ");
162	qtc->Draw();
163	}
164
165
166	}
167
168	//------------------------------------------------------------------------
169	void DrawLEDminCFD()
170	{
171	Int_t npeaks = 10;
172	Int_t sigma=10.;
173	Bool_t down=false;
174	Int_t index[20];
175
176	TCanvas *c1 = new TCanvas("c1", "c1",0,48,1280,951);
177	c1->Divide(4,3);
178
179	Char_t buf1[20];
180	for (Int_t i=0; i<12; i++)
181	{
182	c1->cd(i+1);
183	sprintf(buf1,"LED-CFD%i",i+1);
184	TH1F cfd = (TH1F) gFile->Get(buf1);
185	cfd->Draw();
186	TSpectrum s = new TSpectrum(2npeaks,1);
187	Int_t nfound = s->Search(cfd,sigma," ",0.2);
188	cout<<"Found "<<nfound<<" peaks sigma "<<sigma<<endl;;
189	if(nfound!=0) {
190	Float_t *xpeak = s->GetPositionX();
191	TMath::Sort(nfound, xpeak, index,down);
192	Float_t xp = xpeak[index[0]];
193	Int_t xbin = cfd->GetXaxis()->FindBin(xp);
194	Float_t yp = cfd->GetBinContent(xbin);
195	cout<<"xbin = "<<xbin<<"\txpeak = "<<xpeak[index[0]]<<"\typeak = "<<yp<<endl;
196	Float_t hmin=xp-3*sigma;
197	Float_t hmax =xp+3*sigma;
198	cfd->GetXaxis()->SetRange(hmin,hmax);
199	TF1 *g1 = new TF1("g1", "gaus", hmin, hmax);
200	cfd->Fit("g1","RQ");
201
202
203	}
204
205	}
206	/*
207	TCanvas *c1 = new TCanvas("c1", "c1",0,48,1280,951);
208	c1->Divide(2,2);
209	Char_t buf1[10];
210	for (Int_t i=0; i<4; i++)
211	{
212	c1->cd(i+1);
213	sprintf(buf1,"LED-CFD%i",i+1);
214	TH1F cfd = (TH1F) file->Get(buf1);
215	// cout<<buf1<<" "<<cfd<<endl;
216	// cfd->Draw();
217	// cfd->GetXaxis()->SetRange(0,100);
218	Float_t mean = cfd->GetMean();
219	Float_t rms = cfd->GetRMS();
220	Float_t hmin=mean - 3*rms;
221	Float_t hmax =mean + 3*rms;
222	cfd->GetXaxis()->SetRange(hmin-10,hmax+10);
223	cout<<" cfd range "<<mean<<" rms "<<rms<<" "<<hmin<<" "<<hmax<<endl;
224	// TF1 *g1 = new TF1("g1", "gaus", hmin, hmax);
225	// cfd->Fit("g1","RQ");
226	cfd->Draw();
227	}
228	*/
229
230	}
231	//------------------------------------------------------------------------
232	void DrawCFDvsQTC()
233	{
234	Int_t npeaks = 20;
235	Int_t sigma=3.;
236	Bool_t down=false;
237
238	Int_t index[20];
239	Char_t buf1[10], buf2[10];
240
241	TCanvas *c1 = new TCanvas("c1", "c1",0,48,1280,951);
242	gStyle->SetOptStat(0);
243	c1->Divide(3,2);
244
245	for (Int_t i=0; i<5; i++)
246	{
247	c1->cd(i+1);
248	sprintf(buf1,"T0_C_%i_CFD",i+1);
249	sprintf(buf2,"CFDvsQTC%i",i+1);
250	cout<<buf1<<" "<<buf2<<endl;
251	TH2F qtc_cfd = (TH2F) gFile->Get(buf2);
252	TH1F cfd = (TH1F) gFile->Get(buf1);
253	// cfd->Draw();
254	TSpectrum s = new TSpectrum(2npeaks,1);
255	Int_t nfound = s->Search(cfd,sigma," ",0.05);
256	cout<<"Found "<<nfound<<" peaks sigma "<<sigma<<endl;;
257	if(nfound!=0){
258	Float_t *xpeak = s->GetPositionX();
259	TMath::Sort(nfound, xpeak, index,down);
260	Float_t xp = xpeak[index[0]];
261	Int_t xbin = cfd->GetXaxis()->FindBin(xp);
262	Float_t yp = cfd->GetBinContent(xbin);
263	cout<<"xbin = "<<xbin<<"\txpeak = "<<xpeak[1]<<"\typeak = "<<yp<<endl;
264	Float_t hmax = xp+10*sigma;
265	Float_t hmin = xp-10*sigma;
266	cout<<hmin<< " "<<hmax<<endl;
267	// cfd->GetXaxis()->SetRange(hmin,hmax);
268	// TF1 *g1 = new TF1("g1", "gaus", hmin, hmax);
269	// cfd->Fit("g1","R");
270	Int_t hminbin= qtc_cfd->GetXaxis()->GetFirst();
271	Int_t hmaxbin= qtc_cfd->GetXaxis()->GetLast();
272	Int_t nbins= qtc_cfd->GetXaxis()->GetNbins();
273	cout<<" qtc_cfd "<<hminbin<<" "<<hmaxbin<<" "<<nbins<<endl;
274	// qtc_cfd->Draw();
275	TProfile *pr_y = qtc_cfd->ProfileX();
276
277	pr_y->SetMaximum(hmax);
278	pr_y->SetMinimum(hmin);
279	Int_t np=nbins/20;
280	Double_t *xx = new Double_t[np];
281	Double_t *yy = new Double_t[np];
282	Int_t ng=0;
283	Double_t yg=0;
284	for (Int_t ip=1; ip<nbins; ip++)
285	{
286	if(ip%20 != 0 ) {
287	if (pr_y->GetBinContent(ip) !=0)
288	yg +=pr_y->GetBinContent(ip);
289	// cout<<ng<<" "<<pr_y->GetBinContent(ip)<<" "<<yg<<endl;
290	ng++;}
291	else {
292	xx[ip/20] = Float_t (ip);
293	yy[ip/20] = yg/ng;
294	yg=0;
295	ng=0;
296	// cout<<ip<<" "<<ip/20<<" "<< xx[ip/20]<<" "<< yy[ip/20]<<endl;
297	}
298	}
299	TH2F *hr = new TH2F("hr"," ",np,0,nbins, np, hmin, hmax);
300	hr->Draw();
301	TGraph *gr = new TGraph(np,xx,yy);
302	gr->SetMinimum(hmin);
303	gr->SetMaximum(hmax);
304	gr->SetMarkerStyle(20);
305	gr->Draw("P");
306	// delete [] xx;
307	// delete [] yy;
308
309	// pr_y->Rebin(10);
310	// pr_y->Draw();
311
312	}
313
314	}
315
316
317	}
318	//------------------------------------------------------------------------
319	void DrawCFDvsLED()
320	{
321	Int_t runNumber=1098;
322
323	Int_t npeaks = 20;
324	Int_t sigma=2.;
325	Char_t buf1[10], buf2[10];
326
327	for (Int_t i=0; i<1; i++)
328	{
329	// c1->cd(i+1);
330	sprintf(buf1,"T0_C_%i_CFD",i+1);
331	sprintf(buf2,"CFD_LED%i",i+1);
332
333	TH2F qtc_cfd = (TH2F) gFile->Get(buf2);
334	TH1F cfd = (TH1F) gFile->Get(buf1);
335	// cfd->Draw();
336	TSpectrum s = new TSpectrum(2npeaks,1);
337	Int_t nfound = s->Search(cfd,sigma," ",0.05);
338	cout<<"Found "<<nfound<<" peaks sigma "<<sigma<<endl;;
339	if(nfound!=0){
340	Double_t max=0.0;
341	Double_t tabmax[2] = {0.0, 0.0};
342	Float_t *xpeak = s->GetPositionX();
343	for(Int_t k=0; k<1 ;k++)
344	{
345	Float_t xp = xpeak[k];
346	Int_t xbin = cfd->GetXaxis()->FindBin(xp);
347	Float_t yp = cfd->GetBinContent(xbin);
348	cout<<"xbin = "<<xbin<<"\txpeak = "<<xpeak[k]<<"\typeak = "<<yp<<endl;
349	}
350	Float_t hmin=xp-10*sigma;
351	Float_t hmax =xp+10*sigma;
352	cout<<hmin<< " "<<hmax<<endl;
353	// cfd->GetXaxis()->SetRange(hmin,hmax);
354	// TF1 *g1 = new TF1("g1", "gaus", hmin, hmax);
355	// cfd->Fit("g1","R");
356	Int_t hminbin= qtc_cfd->GetXaxis()->GetFirst();
357	Int_t hmaxbin= qtc_cfd->GetXaxis()->GetLast();
358	Int_t nbins= qtc_cfd->GetXaxis()->GetNbins();
359	cout<<" qtc_cfd "<<hminbin<<" "<<hmaxbin<<" "<<nbins<<endl;
360	// qtc_cfd->Draw();
361	TProfile *pr_y = qtc_cfd->ProfileX();
362
363	pr_y->SetMaximum(hmax);
364	pr_y->SetMinimum(hmin);
365	Int_t np=nbins/20;
366	Double_t *xx = new Double_t[np];
367	Double_t *yy = new Double_t[np];
368	Int_t ng=0;
369	Double_t yg=0;
370	for (Int_t ip=1; ip<nbins; ip++)
371	{
372	if(ip%20 != 0 ) {
373	if (pr_y->GetBinContent(ip) !=0)
374	yg +=pr_y->GetBinContent(ip);
375	// cout<<ng<<" "<<pr_y->GetBinContent(ip)<<" "<<yg<<endl;
376	ng++;}
377	else {
378	xx[ip/20] = Float_t (ip);
379	yy[ip/20] = yg/ng;
380	yg=0;
381	ng=0;
382	cout<<ip<<" "<<ip/20<<" "<< xx[ip/20]<<" "<< yy[ip/20]<<endl;
383	}
384	}
385	TH2F *hr = new TH2F("hr"," ",np,0,nbins, np, hmin, hmax);
386	hr->Draw();
387	TGraph *gr = new TGraph(np,xx,yy);
388	gr->SetMinimum(hmin);
389	gr->SetMaximum(hmax);
390	gr->Draw("P");
391	delete [] xx;
392	delete [] yy;
393
394	// pr_y->Rebin(10);
395	// pr_y->Draw();
396
397	}
398
399	}
400
401
402	}