void DrawMatterEffects(Double_t x,Int_t hier=1){

  gROOT->LoadMacro("NueAna/NueAnaTools/OscProb.cxx");
  TF1 *form1 = new TF1("form1",ElecAppear,1,10,9);
  TF1 *form2 = new TF1("form2",ElecAppear,1,10,9);
  TF1 *form3 = new TF1("form3",ElecAppear,1,10,9);
  
  form1->SetLineColor(1);
  form2->SetLineColor(4);
  form3->SetLineColor(2);
  
  form1->SetParameters(735,0.7854,0.554,0.161,hier*2.7e-3,8.2e-5,
		       2.65,0,1);
  form2->SetParameters(735,0.7854,0.554,0.161,hier*2.7e-3,8.2e-5,
		       2.65,TMath::Pi()/x,1);
  form3->SetParameters(735,0.7854,0.554,0.161,hier*2.7e-3,8.2e-5,
		       2.65,-TMath::Pi()/x,1);
  
  TCanvas *can1 = new TCanvas();
  form2->Draw();
  form1->Draw("same");
  form3->Draw("same");
  
  TH1* h1 = form1->GetHistogram();  
  TH1* h2 = form2->GetHistogram();
  TH1* h3 = form3->GetHistogram();
  TCanvas *can2 = new TCanvas();
  h2->Draw();
  h1->Draw("same");
  h3->Draw("same");


  TCanvas *can3 = new TCanvas();
  can3->SetGridx();
  can3->SetGridy();
  TH1 *a = (TH1*) h2->Clone("a");
  TH1 *b = (TH1*) h3->Clone("b");
  a->Divide(h1);
  b->Divide(h1);
  a->Draw();
  a->SetMaximum(1.5);
  a->SetMinimum(0.5);
  b->Draw("same");
}