#include "constants.txt"

void forum()
{// begin
  
  char text[64];  
  //------------- Style --------------
  gROOT->SetStyle("Plain");
  gStyle->SetOptStat(0);
  double ft = 10 * 10 + 2;
  double sz = 0.055*0.6;

  //-------------------------------------------------------------------------//

  //-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-
  gROOT->ProcessLine(".L functions_new.C"); 
  //-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-

  //low      energy regime: reg = 0
  int reg = 0;

  //- Zenith angle (degrees)
  double za1  = 0.;
  double za   = za1;
  
  int npx = 100;

  //-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-

  double Enulo = 0.080;   // GeV
  double Enuhi = 0.600;   // GeV

  //-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-

  //- Array of crossing times
  first = true;
  GetCrossingTimes(za);

  printf("Total L: %f km \n",(times[2*Nshells]-times[0])*c_light);

  printf("[INFO] New OscP_JC, start! \n");

  TF1 *fOscPmaao_2;
  fOscPmaao_2 = new TF1("fOscPmaao_2",OscP_JC,1.25*Enulo,Enuhi,5);
  fOscPmaao_2->SetParameter(0,za);
  fOscPmaao_2->SetParameter(1,tlohi[0]);
  fOscPmaao_2->SetParameter(2,tlohi[1]);
  fOscPmaao_2->SetParameter(3,1);
  fOscPmaao_2->SetParameter(4,1);  //2nd order Magnus
  fOscPmaao_2->SetLineColor(4);
  fOscPmaao_2->SetLineWidth(1);
  fOscPmaao_2->SetLineStyle(1);
  fOscPmaao_2->SetNpx(npx);

  printf("[INFO] New OscP_JC, done! \n");

  //-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-

  //-- Drawing section
  
  TLatex *lat = new TLatex();
  lat->SetNDC();
  lat->SetTextFont(ft);
  lat->SetTextSize(sz);

  double yto = 1.2;
  double lm  = 0.18;  
  double rm  = 0.05;  
  double bm  = 0.15;

  TH1F *frame_OscP = new TH1F("frame_OscP","",1,Enulo,Enuhi);
  frame_OscP->SetMinimum(0);
  frame_OscP->SetMaximum(0.76);
  frame_OscP->GetXaxis()->SetLabelFont(ft);
  frame_OscP->GetXaxis()->SetLabelSize(sz);
  frame_OscP->GetXaxis()->SetTitleFont(ft);
  frame_OscP->GetXaxis()->SetTitleSize(sz);
  frame_OscP->GetXaxis()->SetTitle("E (GeV)              ");
  frame_OscP->GetYaxis()->SetLabelFont(ft);
  frame_OscP->GetYaxis()->SetLabelSize(sz);
  frame_OscP->GetYaxis()->SetTitleOffset(yto);
  frame_OscP->GetYaxis()->SetTitleFont(ft);
  frame_OscP->GetYaxis()->SetTitleSize(sz);
  frame_OscP->GetYaxis()->SetTitle("P(#nu_{#mu}#rightarrow#nu_{#mu})        "); 

  TLegend *leg1 = new TLegend(0.65,0.65,0.95,0.75);
  leg1->SetTextFont(ft);
  leg1->SetTextSize(sz);
  leg1->SetFillColor(0);
  leg1->SetLineColor(0);

  int ny = 1;
  
  TCanvas *canvOscP = new TCanvas("canvOscP","",1.3*500,ny*500);
  
  canvOscP->Divide(1,1);
  
  leg1->AddEntry(fOscPmaao_2 ,"Magnus O(2)","l");

  canvOscP->cd(1);
  gPad->SetBottomMargin(bm);
  gPad->SetRightMargin(rm);
  gPad->SetLeftMargin(lm);
  gPad->SetTicks(1,1);

  frame_OscP->Draw();
  fOscPmaao_2->Draw("same");

  leg1->Draw();
  sprintf(text,"cos#Theta=%3.1f",cos(za*c_pi/180.));
  lat->DrawLatex(0.65,0.35,text);
  
  canvOscP->Print("test_23Apr12.eps");

}//end