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");
} |
Partager