ToFdEdx_patch/macros/plot_parpos.cxx

{

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


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

Float_t  tp[48][4];

char *s = new char[1];

Int_t idummy;
Float_t  p0,p1,p2,p3;


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


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

TF1 *f1[48];

TGraph *gr[48];


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

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

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

f1[i] = new TF1("f1","[0]+[1]*x+[2]*x*x+[3]*x*x*x",0,1000);
f1[i].SetParameters(tp[i][0],tp[i][1],tp[i][2],tp[i][3]);
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,1000,2,0,1000);
hr->SetXTitle("adccorr");
hr->SetYTitle("f_pos");
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);


TH2F *hr = new TH2F("hr","",2,0,1000,2,0,1000);
hr->SetXTitle("adccorr");
hr->SetYTitle("f_pos");
hr->Draw();


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

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

}


}


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