ToFdEdx_patch/macros/plot_desatBB.cxx

{

TString name1 = "../GoodParam/desaturation_beta.txt";
ifstream fin1(name1);


gROOT->SetStyle("Plain");
gStyle->SetOptStat(0);
gStyle->SetOptFit(0);
gStyle->SetFuncColor( kRed );

Float_t  tp[48][3];

char *s = new char[1];

Int_t idummy;
Float_t  p0,p1,p2;


for (int i=0; i<48; i++) {
fin1>>idummy>>p0>>p1>>p2;
cout<<" "<<idummy<<" "<<p0<<" "<<p1<<" "<<p2<<endl;
tp[i][0]=p0;
tp[i][1]=p1;
tp[i][2]=p2;
}
fin1.close();


TCanvas *c = new TCanvas("c","histograms",800,600);

TF1 *f1[48];

for (Int_t i=0; i<48; i++) {

cout<<"PMT "<<i<<endl;

/*
double ToFNaNuclei::f_desatBB( TArrayF &p, float x )
{
  return 
    p[0] +
    p[1]*x +
    p[2]*x*x;
}

*/

f1[i] = new TF1("f1","[0]+[1]*x+[2]*x*x",0,1000);
f1[i].SetParameters(tp[i][0],tp[i][1],tp[i][2]);
f1[i]->SetLineColor(kRed);


for (Int_t ii=1; ii<20; ii++) {
Float_t xb = ii*20.;
Float_t yhelp = f1[i]->Eval(xb);
cout<<xb<<" "<<yhelp<<endl;
                          }


//TH2F *hr = new TH2F("hr","",2,0,100,2,0,100);
TH2F *hr = new TH2F("hr","",2,0,10,2,0,10);
hr->SetXTitle("adclin");
hr->SetYTitle("f_desatBB");
hr->Draw();

f1[i]->Draw("same");

//gPad->SetLogy();
gPad->SetGrid();

c->Update();
gets(s);

}

//==========================================================================
//==========================================================================
//==========================================================================

TCanvas *c1 = new TCanvas("c1","histograms1",1100,800);
c1->Divide(8,6);

TF1 *f1[48];

for (Int_t i=0; i<48; i++) {

cout<<"PMT "<<i<<endl;

c1->cd(i+1);


hr->Draw();
f1[i]->Draw("same");

//gPad->SetLogy();
gPad->SetGrid();

}


}


}
1	{
2
3	TString name1 = "../GoodParam/desaturation_beta.txt";
4	ifstream fin1(name1);
5
6
7	gROOT->SetStyle("Plain");
8	gStyle->SetOptStat(0);
9	gStyle->SetOptFit(0);
10	gStyle->SetFuncColor( kRed );
11
12	Float_t tp[48][3];
13
14	char *s = new char[1];
15
16	Int_t idummy;
17	Float_t p0,p1,p2;
18
19
20	for (int i=0; i<48; i++) {
21	fin1>>idummy>>p0>>p1>>p2;
22	cout<<" "<<idummy<<" "<<p0<<" "<<p1<<" "<<p2<<endl;
23	tp[i][0]=p0;
24	tp[i][1]=p1;
25	tp[i][2]=p2;
26	}
27	fin1.close();
28
29
30	TCanvas *c = new TCanvas("c","histograms",800,600);
31
32	TF1 *f1[48];
33
34	for (Int_t i=0; i<48; i++) {
35
36	cout<<"PMT "<<i<<endl;
37
38	/*
39	double ToFNaNuclei::f_desatBB( TArrayF &p, float x )
40	{
41	return
42	p[0] +
43	p[1]*x +
44	p[2]xx;
45	}
46
47	*/
48
49	f1[i] = new TF1("f1","[0]+[1]x+[2]x*x",0,1000);
50	f1[i].SetParameters(tp[i][0],tp[i][1],tp[i][2]);
51	f1[i]->SetLineColor(kRed);
52
53
54	for (Int_t ii=1; ii<20; ii++) {
55	Float_t xb = ii*20.;
56	Float_t yhelp = f1[i]->Eval(xb);
57	cout<<xb<<" "<<yhelp<<endl;
58	}
59
60
61
62	//TH2F *hr = new TH2F("hr","",2,0,100,2,0,100);
63	TH2F *hr = new TH2F("hr","",2,0,10,2,0,10);
64	hr->SetXTitle("adclin");
65	hr->SetYTitle("f_desatBB");
66	hr->Draw();
67
68	f1[i]->Draw("same");
69
70	//gPad->SetLogy();
71	gPad->SetGrid();
72
73	c->Update();
74	gets(s);
75
76	}
77
78	//==========================================================================
79	//==========================================================================
80	//==========================================================================
81
82	TCanvas *c1 = new TCanvas("c1","histograms1",1100,800);
83	c1->Divide(8,6);
84
85	TF1 *f1[48];
86
87	for (Int_t i=0; i<48; i++) {
88
89	cout<<"PMT "<<i<<endl;
90
91	c1->cd(i+1);
92
93
94	hr->Draw();
95	f1[i]->Draw("same");
96
97	//gPad->SetLogy();
98	gPad->SetGrid();
99
100	}
101
102
103	}
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120	}