1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
| {
int N=225;
double lamdaalpha[N];
double lamdabeta[N];
double khi2[N],k[N], I[N],J[N];
ifstream f(Form("../2-chi2/val.txt"));// des données a partir d'un fichier txt
TH2D *h = new TH2D("h","h",15,0.05,1.55,15,0.05,1.55);
for(int i=0; i<N; i++){
f>>k[i]>>I[i]>>J[i]>>lamdaalpha[i]>>lamdabeta[i]>>khi2[i];
cout<<k[i]<<I[i]<<J[i]<<lamdaalpha[i]<<" "<<lamdabeta[i]<<" "<<k[i]<<" "<<khi2[i]<<endl;
h->Fill(lamdaalpha[i],lamdabeta[i],khi2[i]);//remplissage de hitogramme
}
TCanvas *c1=new TCanvas("c1","c1",800,600);
c1->SetGridx();
c1->SetGridy();
Double_t contours[3];
contours[0] = 0.5;
contours[1] = 1.;
contours[2] = 2.;
h->GetXaxis()->SetTitle("#lambda_{#alpha}");
h->GetYaxis()->SetTitle("#lambda_{#beta}");
h ->SetContour(3, contours);
h->Draw("cont list");
} |
