{
   gROOT->Reset();
gStyle->SetOptStat(0); 

      Float_t tanb[6000];
      Float_t mh0[6000];
      Float_t mA0[6000];    
      Float_t gtot1[6000];
      Float_t ww2[6000];
      Float_t ww3[6000];
      Float_t gg1[6000];
      Float_t gg3[6000];
      Float_t zz3[6000];
      Float_t sigma_kmr_mh0[6000];


   //open the ASCII file
    FILE *fp = fopen("h0bin2-sa-msusy1000-full.txt","r");

     for(Int_t i=0; i<1701; i++){
    fscanf(fp,"%g %g %g %g %*g %*g %g %*g %g %*g %g %g %*g %*g %*g %*g %*g %*g %*g %*g %g",&tanb[i],&mA0[i],&mh0[i],&gtot1[i],&gg1[i],&gg3[i],&ww2[i],&ww3[i],&zz3[i]);

      if (mh0[i]>0) {
	sigma_kmr_mh0[i]=13333*pow((136/(16+mh0[i])),3.3)*pow((120/mh0[i]),3)*gg1[i]*ww2[i];}
      else {sigma_kmr_mh0[i]=0;}

  }

   fclose(fp);

   TCanvas *c = new TCanvas("c","h->WW",0,0,500,300);
   c->Range(0,0,25,18);
   c->SetFillColor(10);

   
   gStyle->SetHistLineStyle(2);
   TGraph2D *w6 = new TGraph2D(i,mA0,tanb,sigma_kmr_mh0);
   w6->SetNpx(200);
   w6->SetNpy(200);
   TH2D* h6 = w6->GetHistogram("h6");
   double contours[3]={0,0.2,1.0};  
   h6->SetContour(3,contours);
   TH2D* f6 = new TH2D("f6","",100,90,250,100,10,50);
   f6->GetXaxis()->SetTitle("m_{A} [GeV]");
   f6->GetYaxis()->SetTitle("tan #beta");
   f6->SetTitle("Cross section");
   f6->Draw();
   h6->Draw("cont3");

gStyle->SetLineStyleString(21,"20 80");
gStyle->SetHistLineStyle(21);
      TGraph2D *w7 = new TGraph2D(i,mA0,tanb,mh0);
      w7->SetNpx(200);
      w7->SetNpy(200);
      TH2D* h7 = w7->GetHistogram("h7");
      double contoursm[4]={100,120,122,122.8};//sa
      h7->SetContour(4,contoursm);
      h7->Draw("cont3same");

c->Print("marek.ps");

}