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	mayorov	1.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			}