void ErrorsFromProfile(){
  //  TFile *hfile = new TFile("fvtx_mc_eval.root");
m1=new TH1F("m1","DCA resolution vs. Momentum",15,1,12);
e1=new TProfile("e1","DCA resolution vs. Momentum",15,1,12,"s");

TCut rres="abs(x0reco*(px0mc/sqrt(px0mc**2+py0mc**2))+y0reco*(py0mc/sqrt(px0mc**2+py0mc**2)))<.5";
mc_trk_eval->Draw("x0reco*(px0mc/sqrt(px0mc**2+py0mc**2))+y0reco*(py0mc/sqrt(px0mc**2+py0mc**2)):sqrt(px0mc**2+py0mc**2+pz0mc**2)>>e1",rres&&"abs(x0reco)>0&&abs(y0reco)>0&&abs(z0mc)<10");
//sngmuons->Draw("x0reco*(px0mc/sqrt(px0mc**2+py0mc**2))+y0reco*(py0mc/sqrt(px0mc**2+py0mc**2)):sqrt(px0mc**2+py0mc**2+pz0mc**2)>>e1",rres&&"abs(x0reco)>0&&abs(y0reco)>0&&abs(z0mc)<10");
m1->Sumw2();
m1->SetBinContent(0,e1->GetBinError(0));
m1->SetBinContent(1,e1->GetBinError(1));
m1->SetBinContent(2,e1->GetBinError(2));
m1->SetBinContent(3,e1->GetBinError(3));
m1->SetBinContent(4,e1->GetBinError(4));
m1->SetBinContent(5,e1->GetBinError(5));
m1->SetBinContent(6,e1->GetBinError(6));
m1->SetBinContent(7,e1->GetBinError(7));
m1->SetBinContent(8,e1->GetBinError(8));
m1->SetBinContent(9,e1->GetBinError(9));
m1->SetBinContent(10,e1->GetBinError(10));
m1->SetBinContent(11,e1->GetBinError(11));
m1->SetBinContent(12,e1->GetBinError(12));
m1->SetBinContent(13,e1->GetBinError(13));
m1->SetBinContent(14,e1->GetBinError(14));
m1->SetMarkerStyle(2);
m1->Draw("p");
}