{






int NP100=41;


double ex[41];
double ey[41];

for(int i=0;i<41;i++){

ex[i]=0.;
ey[i]=0.;


}






  const double zRF[41 ]={
         0,
      0.05,     0.1,    0.15,     0.2,    0.25,     0.3,    0.35,     0.4,
      0.45,     0.5,    0.55,     0.6,    0.65,     0.7,    0.75,     0.8,
      0.85,     0.9,    0.95,       1,    1.05,     1.1,    1.15,     1.2,
      1.25,     1.3,    1.35,     1.4,    1.45,     1.5,    1.55,     1.6,
      1.65,     1.7,    1.75,     1.8,    1.85,     1.9,    1.95,       2,
  };   

const double RF_0[41 ]={
         0,
         0,       0,       0,       0,       0,       0,       0,       0,
         0,       0,       0,       0,0.00443727,0.00726099,0.0133118,0.014522,
  0.0173457,0.0137152,0.019766,0.0229931,0.0165389,0.0181525,0.014522,0.0133118,
  0.0104881,0.00927793,0.0129084,0.00806777,0.00968132,0.00847116,0.00282372,0.00443727,
  0.00605083,0.0068576,0.00524405,0.0068576,0.00242033,0.00403388,0.00403388,       0,
  };   


TGraphErrors  *gr1 = new TGraphErrors(NP100,zRF,RF_0,ex,ey); 
gr1->SetMarkerColor(2);gr1->SetMarkerSize(.7);gr1->SetMarkerStyle(2);

TH2F *hpx = new TH2F ("hpx","test",100,0,2,100,0,0.1);
TCanvas  *c1 = new TCanvas("c1","test",10,10,700,500);
c1->Divide(1,1,0.00,0.00);




c1_1->SetTickx();c1_1->SetTicky();c1_1->SetBorderMode(0);
c1->cd(1); hpx->Draw();


TF1 *myFit1T1 = new TF1("myFit1T1","TMath::Erf((x-[0])/[1])",0.3,1.5);
myFit1T1->SetParameter(0,1);
myFit1T1->SetLineColor(kBlack);
gr1->Fit("myFit1T1");
gr1->Draw("P");












}