{
gStyle->SetOptStat(kFALSE);
int entries=4;
float cut[entries];
double Response[entries];
double Et[entries];
double Error[entries]
cut[0] = 25;
cut[1] = 50;
cut[2] = 75;
cut[3] = 100;
float range[2];


TFile g( "40GeVTest.2.root", "read" );
g.Get("RbRespEtEta06Prots");
hist2d= (TH2F *)g.Get("RbRespEtEta06Prots");
double start=0;

range[0]=hist2d->GetYaxis()->GetBinUpEdge(0);
range[1]=hist2d->GetYaxis()->GetBinUpEdge(hist2d->GetNbinsY());

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

hist1d = hist2d->ProjectionX("hist1d",start,cut[i]);
Responce[i]= hist1d->GetMean(1);
Error[i] = hist1d->GetRMS(1)/sqrt(hist1d->GetEntries());

hist2d->SetAxisRange(hist2d->GetYaxis()->GetBinUpEdge(start),cut[i],"Y");

Et[i]= hist2d->GetMean(2);
hist2d->SetAxisRange(range[0],range[1],"Y");


start = cut[i];
}

hist2d->Draw();

}