{

Int_t np = 100;

TFile *fdata;
fdata = new TFile("/twist/tw00y/systematics/data/set2/anal9/set2anal9_nxc.root");

Double_t p[np];
Double_t costh[np];
Double_t negcosth[np];
for(ip=0;ip<np;ip++){
  p[ip] = ip*(50.0 - 44.0)/(np - 1) + 44.0;
  costh[ip] = sqrt(1 - pow(38.1,2)/pow(p[ip],2));
  negcosth[ip] = -costh[ip];
//   cout<<"ip="<<p<<"\n";
//   cout<<"p="<<p[ip]<<"\n";
//   cout<<"costh="<<costh[ip]<<"\n";
//   cout<<"\n";
}

TGraph *gtrpos = new TGraph(np,p,costh);
TGraph *gtrneg = new TGraph(np,p,negcosth);

Double_t pz[2] = {16.1, 27.4};
Double_t pzcosth[2] = {0.85, 0.5};
Double_t negpzcosth[2] = {-0.85, -0.5};
TGraph *gzpos = new TGraph(2,pz,pzcosth);
TGraph *gzneg = new TGraph(2,pz,negpzcosth);

Double_t p5[2] = {27.4, 44.0};
Double_t p85[2] = {16.1, 50};
Double_t costh5[2] = {0.5, 0.5};
Double_t costhneg5[2] = {-0.5, -0.5};
Double_t costh85[2] = {0.85, 0.85};
Double_t costhneg85[2] = {-0.85, -0.85};
TGraph *gcos5 = new TGraph(2,p5,costh5);
TGraph *gcos5neg = new TGraph(2,p5,costhneg5);
TGraph *gcos85 = new TGraph(2,p85,costh85);
TGraph *gcos85neg = new TGraph(2,p85,costhneg85);

Double_t highp[2] = {50.0, 50.0};
Double_t highpcosth[2] = {0.65, 0.85};
Double_t highpnegcosth[2] = {-0.65, -0.85};
TGraph *ghighp = new TGraph(2,highp,highpcosth);
TGraph *ghighpneg = new TGraph(2,highp,highpnegcosth);

gtrpos->SetLineWidth(3);
gtrneg->SetLineWidth(3);
gzpos->SetLineWidth(3);
gzneg->SetLineWidth(3);
gcos5->SetLineWidth(3);
gcos5neg->SetLineWidth(3);
gcos85->SetLineWidth(3);
gcos85neg->SetLineWidth(3);
ghighp->SetLineWidth(3);
ghighpneg->SetLineWidth(3);

h0225->SetStats(kFALSE);
h0225->SetTitle("");
h0225->GetXaxis()->SetTitle("Momentum (MeV/c)");
h0225->GetXaxis()->CenterTitle();
h0225->GetYaxis()->SetTitle("cos#theta");
h0225->GetYaxis()->CenterTitle();
h0225->GetYaxis()->SetTitleOffset(1.1);

gStyle->SetPalette(1);
TCanvas *cfid;
cfid = new TCanvas("cfid","fiducial",400,400);
cfid->SetFillColor(0);
cfid->SetBorderSize(0.0);
cfid->SetBorderMode(0.0);
cfid->SetTopMargin(0.02);
cfid->SetRightMargin(0.02);
h0225->Draw("col");
gtrpos->Draw();
gtrneg->Draw();
gzpos->Draw();
gzneg->Draw();
gcos5->Draw();
gcos5neg->Draw();
gcos85->Draw();
gcos85neg->Draw();
ghighp->Draw();
ghighpneg->Draw();

}