gXSecComp.cxx File Reference

#include <cassert>
#include <sstream>
#include <string>
#include <vector>
#include <TSystem.h>
#include <TFile.h>
#include <TDirectory.h>
#include <TGraph.h>
#include <TPostScript.h>
#include <TH1D.h>
#include <TMath.h>
#include <TCanvas.h>
#include <TPavesText.h>
#include <TText.h>
#include <TStyle.h>
#include <TLegend.h>
#include <TObjString.h>
#include "Framework/Conventions/GBuild.h"
#include "Framework/Messenger/Messenger.h"
#include "Framework/Numerical/RandomGen.h"
#include "Framework/ParticleData/PDGUtils.h"
#include "Framework/ParticleData/PDGCodes.h"
#include "Framework/Utils/CmdLnArgParser.h"
#include "Framework/Utils/StringUtils.h"
#include "Framework/Utils/Style.h"

Include dependency graph for gXSecComp.cxx:

Go to the source code of this file.

Functions
void	Init (void)
void	End (void)
void	OpenDir (void)
void	DirNameToProbe (void)
void	DirNameToTarget (void)
void	MakePlots (void)
void	MakePlotsCurrDir (void)
TH1F *	DrawFrame (TGraph *gr0, TGraph *gr1, TPad *p, const char *xt, const char *yt, double yminsc, double ymaxsc)
TH1F *	DrawFrame (double xmin, double xmax, double ymin, double ymax, TPad *p, const char *xt, const char *yt)
void	Draw (const char *plot, const char *title)
void	Draw (TGraph *gr, const char *opt)
TGraph *	TrimGraph (TGraph *gr, int max_np_per_decade)
TGraph *	DrawRatio (TGraph *gr0, TGraph *gr1)
void	GetCommandLineArgs (int argc, char **argv)
void	PrintSyntax (void)
bool	CheckRootFilename (string filename)
string	OutputFileName (string input_file_name)
int	main (int argc, char **argv)

Variables
string	gOptXSecFilename_curr = ""
string	gOptXSecFilename_ref0 = ""
string	gOptOutputFilename = ""
bool	gOptHaveRef
TFile *	gXSecFile_curr = 0
TFile *	gXSecFile_ref0 = 0
TDirectory *	gDirCurr = 0
TDirectory *	gDirRef0 = 0
string	gLabelCurr = ""
string	gLabelRef0 = ""
string	gDirName = ""
TPostScript *	gPS = 0
TCanvas *	gC = 0
TPad *	gPadTitle = 0
TPad *	gPadXSecs = 0
TPad *	gPadRatio = 0
TLegend *	gLS = 0
string	gCurrProbeLbl = ""
int	gCurrProbePdg = 0
bool	gCurrProbeIsNu = false
bool	gCurrProbeIsNuBar = false
string	gCurrTargetLbl = ""
bool	gCurrTargetHasP = false
bool	gCurrTargetHasN = false
bool	gCurrTargetIsFreeNuc = false

Function Documentation

CheckRootFilename()

bool CheckRootFilename

(

string

filename

)

Definition at line 807 of file gXSecComp.cxx.

{
  if(filename.size() == 0) return false;
 
  bool is_accessible = ! (gSystem->AccessPathName(filename.c_str()));
  if (!is_accessible) {
   LOG("gxscomp", pERROR)
       << "The input ROOT file [" << filename << "] is not accessible";
   return false;
  }
  return true;
}

References LOG, and pERROR.

Referenced by GetCommandLineArgs().

DirNameToProbe()

void DirNameToProbe

(

void

)

Definition at line 232 of file gXSecComp.cxx.

{
// Figure out the probe type from the input directory name
//
 
  string dirname = gDirName;
 
  int    pdg     = 0;
  string label   = "";
 
  if( dirname.find ("nu_e_bar") != string::npos ) 
  { 
    label = "#bar{#nu_{e}}";    
    pdg   = kPdgAntiNuE;   
  }
  else 
  if ( dirname.find ("nu_e") != string::npos ) 
  { 
    label = "#nu_{e}";          
    pdg   = kPdgNuE;       
  }
  else 
  if ( dirname.find ("nu_mu_bar") != string::npos ) 
  { 
    label = "#bar{#nu_{#mu}}";  
    pdg   = kPdgAntiNuMu;  
  }
  else 
  if ( dirname.find ("nu_mu") != string::npos ) 
  { 
    label ="#nu_{#mu}";       
    pdg   = kPdgNuMu;      
  }
  else 
  if ( dirname.find ("nu_tau_bar") != string::npos ) 
  { 
    label = "#bar{#nu_{#tau}}"; 
    pdg   = kPdgAntiNuTau; 
  }
  else 
  if ( dirname.find ("nu_tau") != string::npos ) 
  { 
    label = "#nu_{#tau}";      
    pdg   = kPdgNuTau;     
  }
 
  gCurrProbeLbl     = label;
  gCurrProbePdg     = pdg;
  gCurrProbeIsNu    = pdg::IsNeutrino     (pdg);
  gCurrProbeIsNuBar = pdg::IsAntiNeutrino (pdg);
}

References gCurrProbeIsNu, gCurrProbeIsNuBar, gCurrProbeLbl, gCurrProbePdg, gDirName, genie::pdg::IsAntiNeutrino(), genie::pdg::IsNeutrino(), genie::kPdgAntiNuE, genie::kPdgAntiNuMu, genie::kPdgAntiNuTau, genie::kPdgNuE, genie::kPdgNuMu, and genie::kPdgNuTau.

Referenced by MakePlotsCurrDir().

DirNameToTarget()

void DirNameToTarget

(

void

)

Definition at line 284 of file gXSecComp.cxx.

{
// Figure out the target type from the input directory name
//
 
  string label;
  string dirname = gDirName;
 
  if ( gCurrProbePdg == kPdgAntiNuE )  
  { 
    label = dirname.substr ( dirname.find("nu_e_bar")+9, dirname.length() ); 
  }
  else 
  if ( gCurrProbePdg == kPdgNuE )  
  { 
    label = dirname.substr ( dirname.find("nu_e")+5, dirname.length() ); 
  }
  else 
  if ( gCurrProbePdg == kPdgAntiNuMu )  
  { 
    label = dirname.substr ( dirname.find("nu_mu_bar")+10, dirname.length() ); 
  }
  else 
  if ( gCurrProbePdg == kPdgNuMu )  
  { 
    label = dirname.substr ( dirname.find("nu_mu")+6, dirname.length() ); 
  }
  else 
  if ( gCurrProbePdg == kPdgAntiNuTau )  
  { 
    label = dirname.substr ( dirname.find("nu_tau_bar")+11, dirname.length() ); 
  }
  else 
  if ( gCurrProbePdg == kPdgNuTau )  
  { 
    label = dirname.substr ( dirname.find("nu_tau")+7, dirname.length() ); 
  }
 
  bool free_n   = (label.size()==1 && label.find("n") !=string::npos);
  bool free_p   = (label.size()==2 && label.find("H1")!=string::npos);
  bool free_nuc = free_p || free_n;
 
  gCurrTargetHasP = (free_n) ? false : true;
  gCurrTargetHasN = (free_p) ? false : true;
 
  gCurrTargetLbl  = (free_nuc) ? "" : ("(" + label + ")");
 
  gCurrTargetIsFreeNuc  = free_nuc;
}

References gCurrProbePdg, gCurrTargetHasN, gCurrTargetHasP, gCurrTargetIsFreeNuc, gCurrTargetLbl, gDirName, genie::kPdgAntiNuE, genie::kPdgAntiNuMu, genie::kPdgAntiNuTau, genie::kPdgNuE, genie::kPdgNuMu, and genie::kPdgNuTau.

Referenced by MakePlotsCurrDir().

Draw() [1/2]

void Draw	(	const char *	plot,
		const char *	title )

Definition at line 580 of file gXSecComp.cxx.

{
  gPadTitle->cd();
  TPavesText hdr(10,10,90,90,1,"tr");
  hdr.AddText(title);
  hdr.SetFillColor(kWhite);
  hdr.Draw();
 
  gPadXSecs->cd();
  TGraph * gr_curr = (TGraph*) gDirCurr->Get(plot);
  TGraph * gr_ref0 = 0;
  if(gDirRef0) {
      gr_ref0 = (TGraph*) gDirRef0->Get(plot);
  }
 
  if(gr_curr == 0 && gr_ref0 ==0) return;
 
  if(!title) return;
 
  // trim points in reference plot (shown with markers) so that the markers
  // don't hide the current prediction (shown with a line).
  // Keep, at most, 20 points per decade.
  TGraph * gr_ref0_trim = TrimGraph(gr_ref0, 20);
 
  utils::style::Format(gr_curr,      1,1,1,1,1,1.0);
  utils::style::Format(gr_ref0_trim, 1,1,1,2,4,0.7);
 
  DrawFrame (gr_curr, gr_ref0, gPadXSecs, "E_{#nu} (GeV)", "#sigma (10^{-38} cm^{2})", 0.5, 1.5);
  Draw (gr_curr,      "L");
  Draw (gr_ref0_trim, "P");
 
  gLS->Clear();
  gLS->AddEntry(gr_curr, gLabelCurr.c_str(), "L");
  if(gOptHaveRef) {
    gLS->AddEntry(gr_ref0_trim, gLabelRef0.c_str(), "P");
  }
  gLS->Draw();
 
  // plot ratio of current and reference models
  if(gOptHaveRef) 
  {
    gPadRatio->cd();
    TGraph * gr_ratio = DrawRatio(gr_curr, gr_ref0);
    DrawFrame (gr_ratio, 0, gPadRatio, "E_{#nu} (GeV)", Form("%s / %s", gLabelCurr.c_str(), gLabelRef0.c_str()), 0.9, 1.1);
    Draw      (gr_ratio, "L");
  }
 
  gC ->Update();
 
  if(gr_ref0_trim) {
    delete gr_ref0_trim;
  }
 
  gPS->NewPage();
}

References Draw(), DrawFrame(), DrawRatio(), genie::utils::style::Format(), gC, gDirCurr, gDirRef0, gLabelCurr, gLabelRef0, gLS, gOptHaveRef, gPadRatio, gPadTitle, gPadXSecs, gPS, and TrimGraph().

Referenced by Draw(), MakePlots(), and MakePlotsCurrDir().

Draw() [2/2]

void Draw	(	TGraph *	gr,
		const char *	opt )

Definition at line 574 of file gXSecComp.cxx.

{
  if(!gr) return;
  gr->Draw(opt);
}

DrawFrame() [1/2]

TH1F * DrawFrame	(	double	xmin,
		double	xmax,
		double	ymin,
		double	ymax,
		TPad *	p,
		const char *	xt,
		const char *	yt )

Definition at line 561 of file gXSecComp.cxx.

{
  if(!p) return 0;
  TH1F * hf = (TH1F*) p->DrawFrame(xmin, ymin, xmax, ymax);
  hf->GetXaxis()->SetTitle(xt);
  hf->GetYaxis()->SetTitle(yt);
  hf->GetYaxis()->SetTitleSize(0.03);
  hf->GetYaxis()->SetTitleOffset(1.5);
  hf->GetXaxis()->SetLabelSize(0.03);
  hf->GetYaxis()->SetLabelSize(0.03);
  return hf;
}

DrawFrame() [2/2]

TH1F * DrawFrame	(	TGraph *	gr0,
		TGraph *	gr1,
		TPad *	p,
		const char *	xt,
		const char *	yt,
		double	yminsc,
		double	ymaxsc )

Definition at line 529 of file gXSecComp.cxx.

{
  double xmin = 1E-5;
  double xmax = 1;
  double ymin = 1E-5;
  double ymax = 1;
 
  if(gr0) {  
    TAxis * x0 = gr0 -> GetXaxis();
    TAxis * y0 = gr0 -> GetYaxis();
    xmin = x0 -> GetXmin();
    xmax = x0 -> GetXmax();
    ymin = y0 -> GetXmin();
    ymax = y0 -> GetXmax();
  }
  if(gr1) {
     TAxis * x1 = gr1 -> GetXaxis();
     TAxis * y1 = gr1 -> GetYaxis();
     xmin = TMath::Min(xmin, x1 -> GetXmin());
     xmax = TMath::Max(xmax, x1 -> GetXmax());
     ymin = TMath::Min(ymin, y1 -> GetXmin());
     ymax = TMath::Max(ymax, y1 -> GetXmax());
  }
  xmin *= 0.5;
  xmax *= 1.5;
  ymin *= yminsc;
  ymax *= ymaxsc;
  xmin = TMath::Max(0.1, xmin);
  
  return DrawFrame(xmin, xmax, ymin, ymax, p, xt, yt);
}

References DrawFrame().

Referenced by Draw(), and DrawFrame().

DrawRatio()

TGraph * DrawRatio	(	TGraph *	gr0,
		TGraph *	gr1 )

Definition at line 698 of file gXSecComp.cxx.

{
// gr0 / gr1
 
  if(!gr0) return 0;
  if(!gr1) return 0;
 
  LOG("gxscomp", pDEBUG) << "Drawing ratio...";
 
  const int np = gr0->GetN();
 
  double * x = new double[np];
  double * y = new double[np];
 
  for(int i=0; i<np; i++) {
     x[i] = gr0->GetX()[i];
     y[i] = 0;
 
     if(gr0->Eval(x[i]) != 0. && gr1->Eval(x[i]) != 0.) 
     {
       y[i] = gr0->Eval(x[i]) / gr1->Eval(x[i]);
     }
     else {
       if(gr0->Eval(x[i]) != 0.) 
       {
         y[i] = -1;
       }
       else 
       {
         y[i] =  1;
       }
     }
  }
 
  TGraph * ratio = new TGraph(np,x,y);
 
  delete [] x;
  delete [] y;
 
  return ratio;
}

References LOG, and pDEBUG.

Referenced by Draw().

End()

void End

(

void

)

Definition at line 210 of file gXSecComp.cxx.

{
  gPS->Close();
 
  delete gC;
  delete gPS;
  delete gLS;
}

References gC, gLS, and gPS.

Referenced by main().

GetCommandLineArgs()

void GetCommandLineArgs	(	int	argc,
		char **	argv )

Definition at line 740 of file gXSecComp.cxx.

{
  LOG("gxscomp", pINFO) << "*** Parsing command line arguments";
 
  CmdLnArgParser parser(argc,argv);
 
  // get input GENIE cross section file
  if( parser.OptionExists('f') ) {
    string inp = parser.ArgAsString('f');
    if(inp.find(",") != string::npos) {
       vector<string> inpv = utils::str::Split(inp,",");
       assert(inpv.size()==2);
       gOptXSecFilename_curr = inpv[0];
       gLabelCurr = inpv[1];
    } else {
       gOptXSecFilename_curr = inp;
       gLabelCurr = "current";
    }
    bool ok = CheckRootFilename(gOptXSecFilename_curr.c_str());
    if(!ok) {
      PrintSyntax();
      exit(1);
    }
  } else {
    PrintSyntax();
    exit(1);
  }
 
  // get [reference] input GENIE cross section file
  if( parser.OptionExists('r') ) {
    string inp = parser.ArgAsString('r');
    if(inp.find(",") != string::npos) {
       vector<string> inpv = utils::str::Split(inp,",");
       assert(inpv.size()==2);
       gOptXSecFilename_ref0 = inpv[0];
       gLabelRef0 = inpv[1];
    } else {
       gOptXSecFilename_ref0 = inp;
       gLabelRef0 = "reference";
    }
    bool ok = CheckRootFilename(gOptXSecFilename_ref0.c_str());
    if(!ok) {
      PrintSyntax();
      exit(1);
    }
    gOptHaveRef = true;
  } else {
    LOG("gxscomp", pNOTICE) << "No reference cross section file";
    gOptHaveRef = false;
  }
 
  // get output filename
  if( parser.OptionExists('o') ) {
    gOptOutputFilename = parser.ArgAsString('o');
  } else {
    gOptOutputFilename = "xsec.ps";
  }
 
}

References genie::CmdLnArgParser::ArgAsString(), CheckRootFilename(), gLabelCurr, gLabelRef0, gOptHaveRef, gOptOutputFilename, gOptXSecFilename_curr, gOptXSecFilename_ref0, LOG, genie::CmdLnArgParser::OptionExists(), pINFO, pNOTICE, PrintSyntax(), and genie::utils::str::Split().

Referenced by main().

Init()

void Init

(

void

)

Definition at line 138 of file gXSecComp.cxx.

{
  gC = new TCanvas("c","",20,20,500,650);
  gC->SetBorderMode(0);
  gC->SetFillColor(0);
 
  // create output file
  gPS = new TPostScript(gOptOutputFilename.c_str(), 111);
 
  // front page
  gPS->NewPage();
  gC->Range(0,0,100,100);
  TPavesText hdr(10,40,90,70,3,"tr");
  hdr.AddText("GENIE cross sections");
  hdr.AddText(" ");
  hdr.AddText(" ");
  hdr.AddText("Plotting data from: ");
  hdr.AddText(gOptXSecFilename_curr.c_str());
  if(gOptHaveRef) {
     hdr.AddText(" ");
     hdr.AddText("Comparing with reference data (red circles) from: ");
     hdr.AddText(gOptXSecFilename_ref0.c_str());
  } else {
     hdr.AddText(" ");
  }
  hdr.AddText(" ");
  hdr.AddText(" ");
  hdr.Draw();
  gC->Update();
 
  gPS->NewPage();
 
  //
  if(gOptHaveRef) {
    gPadTitle = new TPad("gPadTitle","",0.05,0.90,0.95,0.97);
    gPadXSecs = new TPad("gPadXSecs","",0.05,0.35,0.95,0.88);
    gPadRatio = new TPad("gPadRatio","",0.05,0.03,0.95,0.32);
  } else {
    gPadTitle = new TPad("gPadTitle","",0.05,0.90,0.95,0.97);
    gPadXSecs = new TPad("gPadXSecs","",0.05,0.05,0.95,0.88);
    gPadRatio = new TPad("gPadRatio","",0.05,0.03,0.95,0.04);
  }
 
  gPadTitle->Range(0,0,100,100);
 
  gPadTitle->SetBorderMode(0);
  gPadTitle->SetFillColor(0);
  gPadXSecs->SetBorderMode(0);
  gPadXSecs->SetFillColor(0);
  gPadRatio->SetBorderMode(0);
  gPadRatio->SetFillColor(0);
 
  gPadXSecs->SetGridx();
  gPadXSecs->SetGridy();
  gPadXSecs->SetLogx();
  gPadXSecs->SetLogy();
  gPadRatio->SetGridx();
  gPadRatio->SetGridy();
  gPadRatio->SetLogx();
 
  gPadTitle->Draw();
  gPadXSecs->Draw();
  gPadRatio->Draw();
 
  gPadXSecs->cd();
 
  gLS = new TLegend(0.80,0.25,0.90,0.45);
  gLS -> SetFillColor(0);
//gLS -> SetFillStyle(0);
  gLS -> SetBorderSize(0);
}

References gC, gLS, gOptHaveRef, gOptOutputFilename, gOptXSecFilename_curr, gOptXSecFilename_ref0, gPadRatio, gPadTitle, gPadXSecs, and gPS.

Referenced by genie::BLI2DNonUnifGrid::BLI2DNonUnifGrid(), genie::BLI2DNonUnifGrid::BLI2DNonUnifGrid(), genie::BLI2DUnifGrid::BLI2DUnifGrid(), genie::BLI2DUnifGrid::BLI2DUnifGrid(), ClassImp(), ClassImp(), ClassImp(), genie::GSimFiles::GSimFiles(), gst::gst(), and main().

main()

int main	(	int	argc,
		char **	argv )

Definition at line 125 of file gXSecComp.cxx.

{
  GetCommandLineArgs (argc,argv);
  utils::style::SetDefaultStyle();
  Init();
  MakePlots();
  End();
  
  LOG("gxscomp", pINFO)  << "Done!";
 
  return 0;
}

References End(), GetCommandLineArgs(), Init(), LOG, MakePlots(), pINFO, and genie::utils::style::SetDefaultStyle().

MakePlots()

void MakePlots

(

void

)

Definition at line 334 of file gXSecComp.cxx.

{
  // Open files
  gXSecFile_curr = new TFile(gOptXSecFilename_curr.c_str(), "read");
  gXSecFile_ref0 = 0;
  if (gOptHaveRef) {
    gXSecFile_ref0 = new TFile(gOptXSecFilename_ref0.c_str(), "read");
  }
 
  // Get list of directories in the input file
  TList * directories = gXSecFile_curr->GetListOfKeys();
 
  // Loop over directories & generate plots for each one
  int ndir = directories->GetEntries();
  for(int idir=0; idir<ndir; idir++) {
    TObjString * dir = (TObjString*) directories->At(idir);
    gDirName = dir->GetString().Data();
    MakePlotsCurrDir();
  }
}

References dir, gDirName, gOptHaveRef, gOptXSecFilename_curr, gOptXSecFilename_ref0, gXSecFile_curr, gXSecFile_ref0, and MakePlotsCurrDir().

Referenced by main().

MakePlotsCurrDir()

void MakePlotsCurrDir

(

void

)

Definition at line 355 of file gXSecComp.cxx.

{
  LOG("gxscomp", pINFO)  << "Plotting graphs from directory: " << gDirName;
 
  OpenDir         ();
  DirNameToProbe  ();
  DirNameToTarget ();
 
  const char * probestr  = gCurrProbeLbl.c_str();
  const char * targetstr = gCurrTargetLbl.c_str();
 
  LOG("gxscomp", pINFO)  << "Probe  : " << probestr;
  LOG("gxscomp", pINFO)  << "Target : " << targetstr;
 
  //
  // Start plotting...
  //
 
  if(!gCurrTargetIsFreeNuc) {
    Draw("tot_cc", Form("%s + %s, TOT CC",probestr,targetstr));
    Draw("tot_nc", Form("%s + %s, TOT NC",probestr,targetstr));
  }
 
  if(gCurrTargetHasN) {
      Draw("tot_cc_n",               Form("%s + n %s, TOT CC"                                , probestr, targetstr));
      Draw("tot_nc_n",               Form("%s + n %s, TOT NC"                                , probestr, targetstr));
    if(gCurrProbeIsNu) {
      Draw("qel_cc_n",               Form("%s + n %s, QEL CC"                                , probestr, targetstr));
    }
      Draw("qel_nc_n",               Form("%s + n %s, NCEL"                                  , probestr, targetstr));
      Draw("res_cc_n",               Form("%s + n %s, RES CC"                                , probestr, targetstr));
      Draw("res_nc_n",               Form("%s + n %s, RES NC"                                , probestr, targetstr));
      Draw("res_cc_n_1232P33",       Form("%s + n %s, RES CC, P33(1232)"                     , probestr, targetstr));
      Draw("res_cc_n_1535S11",       Form("%s + n %s, RES CC, S11(1535)"                     , probestr, targetstr));
      Draw("res_cc_n_1520D13",       Form("%s + n %s, RES CC, D13(1520)"                     , probestr, targetstr));
      Draw("res_cc_n_1650S11",       Form("%s + n %s, RES CC, S11(1650)"                     , probestr, targetstr));
      Draw("res_cc_n_1700D13",       Form("%s + n %s, RES CC, D13(1700)"                     , probestr, targetstr));
      Draw("res_cc_n_1675D15",       Form("%s + n %s, RES CC, D15(1675)"                     , probestr, targetstr));
      Draw("res_cc_n_1620S31",       Form("%s + n %s, RES CC, S31(1620)"                     , probestr, targetstr));
      Draw("res_cc_n_1700D33",       Form("%s + n %s, RES CC, D33(1700)"                     , probestr, targetstr));
      Draw("res_cc_n_1440P11",       Form("%s + n %s, RES CC, P11(1440)"                     , probestr, targetstr));
      Draw("res_cc_n_1720P13",       Form("%s + n %s, RES CC, P13(1720)"                     , probestr, targetstr));
      Draw("res_cc_n_1680F15",       Form("%s + n %s, RES CC, F15(1680)"                     , probestr, targetstr));
      Draw("res_cc_n_1910P31",       Form("%s + n %s, RES CC, P31(1910)"                     , probestr, targetstr));
      Draw("res_cc_n_1920P33",       Form("%s + n %s, RES CC, P33(1920)"                     , probestr, targetstr));
      Draw("res_cc_n_1905F35",       Form("%s + n %s, RES CC, F35(1905)"                     , probestr, targetstr));
      Draw("res_cc_n_1950F37",       Form("%s + n %s, RES CC, F37(1950)"                     , probestr, targetstr));
      Draw("res_cc_n_1710P11",       Form("%s + n %s, RES CC, P11(1710)"                     , probestr, targetstr));
      Draw("res_nc_n_1232P33",       Form("%s + n %s, RES NC, P33(1232)"                     , probestr, targetstr));
      Draw("res_nc_n_1535S11",       Form("%s + n %s, RES NC, S11(1535)"                     , probestr, targetstr));
      Draw("res_nc_n_1520D13",       Form("%s + n %s, RES NC, D13(1520)"                     , probestr, targetstr));
      Draw("res_nc_n_1650S11",       Form("%s + n %s, RES NC, S11(1650)"                     , probestr, targetstr));
      Draw("res_nc_n_1700D13",       Form("%s + n %s, RES NC, D13(1700)"                     , probestr, targetstr));
      Draw("res_nc_n_1675D15",       Form("%s + n %s, RES NC, D15(1675)"                     , probestr, targetstr));
      Draw("res_nc_n_1620S31",       Form("%s + n %s, RES NC, S31(1620)"                     , probestr, targetstr));
      Draw("res_nc_n_1700D33",       Form("%s + n %s, RES NC, D33(1700)"                     , probestr, targetstr));
      Draw("res_nc_n_1440P11",       Form("%s + n %s, RES NC, P11(1440)"                     , probestr, targetstr));
      Draw("res_nc_n_1720P13",       Form("%s + n %s, RES NC, P13(1720)"                     , probestr, targetstr));
      Draw("res_nc_n_1680F15",       Form("%s + n %s, RES NC, F15(1680)"                     , probestr, targetstr));
      Draw("res_nc_n_1910P31",       Form("%s + n %s, RES NC, P31(1910)"                     , probestr, targetstr));
      Draw("res_nc_n_1920P33",       Form("%s + n %s, RES NC, P33(1920)"                     , probestr, targetstr));
      Draw("res_nc_n_1905F35",       Form("%s + n %s, RES NC, F35(1905)"                     , probestr, targetstr));
      Draw("res_nc_n_1950F37",       Form("%s + n %s, RES NC, F37(1950)"                     , probestr, targetstr));
      Draw("res_nc_n_1710P11",       Form("%s + n %s, RES NC, P11(1710)"                     , probestr, targetstr));
      Draw("dis_cc_n",               Form("%s + n %s, DIS CC"                                , probestr, targetstr));
      Draw("dis_nc_n",               Form("%s + n %s, DIS NC"                                , probestr, targetstr));
    if(gCurrProbeIsNu) {    
      Draw("dis_cc_n_ubarsea",       Form("%s + n %s, DIS CC (#bar{u}_{sea})"                , probestr, targetstr));
      Draw("dis_cc_n_dval",          Form("%s + n %s, DIS CC (d_{val})"                      , probestr, targetstr));
      Draw("dis_cc_n_dsea",          Form("%s + n %s, DIS CC (d_{sea})"                      , probestr, targetstr));
      Draw("dis_cc_n_ssea",          Form("%s + n %s, DIS CC (s_{sea})"                      , probestr, targetstr));
    }
    if(gCurrProbeIsNuBar) {    
      Draw("dis_cc_n_sbarsea",       Form("%s + n %s, DIS CC (#bar{s}_{sea})"                , probestr, targetstr));
      Draw("dis_cc_n_dbarsea",       Form("%s + n %s, DIS CC (#bar{d}_{val})"                , probestr, targetstr));
      Draw("dis_cc_n_uval",          Form("%s + n %s, DIS CC (u_{val})"                      , probestr, targetstr));
      Draw("dis_cc_n_usea",          Form("%s + n %s, DIS CC (u_{sea})"                      , probestr, targetstr));
    }
      Draw("dis_nc_n_sbarsea",       Form("%s + n %s, DIS NC (#bar{s}_{sea})"                , probestr, targetstr));
      Draw("dis_nc_n_ubarsea",       Form("%s + n %s, DIS NC (#bar{u}_{sea})"                , probestr, targetstr));
      Draw("dis_nc_n_dbarsea",       Form("%s + n %s, DIS NC (#bar{d}_{sea})"                , probestr, targetstr));
      Draw("dis_nc_n_dval",          Form("%s + n %s, DIS NC (d_{val})"                      , probestr, targetstr));
      Draw("dis_nc_n_dsea",          Form("%s + n %s, DIS NC (d_{sea})"                      , probestr, targetstr));
      Draw("dis_nc_n_uval",          Form("%s + n %s, DIS NC (u_{val})"                      , probestr, targetstr));
      Draw("dis_nc_n_usea",          Form("%s + n %s, DIS NC (u_{sea})"                      , probestr, targetstr));
      Draw("dis_nc_n_ssea",          Form("%s + n %s, DIS NC (s_{sea})"                      , probestr, targetstr));
    if(gCurrProbeIsNu) {    
      Draw("dis_cc_n_dval_charm",    Form("%s + n %s, DIS CC (d_{val} -> c)"                 , probestr, targetstr));
      Draw("dis_cc_n_dsea_charm",    Form("%s + n %s, DIS CC (d_{sea} -> c)"                 , probestr, targetstr));
      Draw("dis_cc_n_ssea_charm",    Form("%s + n %s, DIS CC (s_{sea} -> c)"                 , probestr, targetstr));
    }
    if(gCurrProbeIsNuBar) {    
      Draw("dis_cc_n_dbarsea_charm", Form("%s + n %s, DIS CC (#bar{d}_{sea} -> #bar{c})"     , probestr, targetstr));
      Draw("dis_cc_n_sbarsea_charm", Form("%s + n %s, DIS CC (#bar{s}_{sea} -> #bar{c})"     , probestr, targetstr));
    }
  }//N>0?
 
 
  if(gCurrTargetHasP) {
      Draw("tot_cc_p",               Form("%s + p %s, TOT CC"                                , probestr, targetstr));
      Draw("tot_nc_p",               Form("%s + p %s, TOT NC"                                , probestr, targetstr));
    if(gCurrProbeIsNuBar) {
      Draw("qel_cc_p",               Form("%s + p %s, QEL CC"                                , probestr, targetstr));
    }
      Draw("qel_nc_p",               Form("%s + p %s, NCEL"                                  , probestr, targetstr));
      Draw("res_cc_p",               Form("%s + p %s, RES CC"                                , probestr, targetstr));
      Draw("res_nc_p",               Form("%s + p %s, RES NC"                                , probestr, targetstr));
      Draw("res_cc_p_1232P33",       Form("%s + p %s, RES CC, P33(1232)"                     , probestr, targetstr));
      Draw("res_cc_p_1535S11",       Form("%s + p %s, RES CC, S11(1535)"                     , probestr, targetstr));
      Draw("res_cc_p_1520D13",       Form("%s + p %s, RES CC, D13(1520)"                     , probestr, targetstr));
      Draw("res_cc_p_1650S11",       Form("%s + p %s, RES CC, S11(1650)"                     , probestr, targetstr));
      Draw("res_cc_p_1700D13",       Form("%s + p %s, RES CC, D13(1700)"                     , probestr, targetstr));
      Draw("res_cc_p_1675D15",       Form("%s + p %s, RES CC, D15(1675)"                     , probestr, targetstr));
      Draw("res_cc_p_1620S31",       Form("%s + p %s, RES CC, S31(1620)"                     , probestr, targetstr));
      Draw("res_cc_p_1700D33",       Form("%s + p %s, RES CC, D33(1700)"                     , probestr, targetstr));
      Draw("res_cc_p_1440P11",       Form("%s + p %s, RES CC, P11(1440)"                     , probestr, targetstr));
      Draw("res_cc_p_1720P13",       Form("%s + p %s, RES CC, P13(1720)"                     , probestr, targetstr));
      Draw("res_cc_p_1680F15",       Form("%s + p %s, RES CC, F15(1680)"                     , probestr, targetstr));
      Draw("res_cc_p_1910P31",       Form("%s + p %s, RES CC, P31(1910)"                     , probestr, targetstr));
      Draw("res_cc_p_1920P33",       Form("%s + p %s, RES CC, P33(1920)"                     , probestr, targetstr));
      Draw("res_cc_p_1905F35",       Form("%s + p %s, RES CC, F35(1905)"                     , probestr, targetstr));
      Draw("res_cc_p_1950F37",       Form("%s + p %s, RES CC, F37(1950)"                     , probestr, targetstr));
      Draw("res_cc_p_1710P11",       Form("%s + p %s, RES CC, P11(1710)"                     , probestr, targetstr));
      Draw("res_nc_p_1232P33",       Form("%s + p %s, RES NC, P33(1232)"                     , probestr, targetstr));
      Draw("res_nc_p_1535S11",       Form("%s + p %s, RES NC, S11(1535)"                     , probestr, targetstr));
      Draw("res_nc_p_1520D13",       Form("%s + p %s, RES NC, D13(1520)"                     , probestr, targetstr));
      Draw("res_nc_p_1650S11",       Form("%s + p %s, RES NC, S11(1650)"                     , probestr, targetstr));
      Draw("res_nc_p_1700D13",       Form("%s + p %s, RES NC, D13(1700)"                     , probestr, targetstr));
      Draw("res_nc_p_1675D15",       Form("%s + p %s, RES NC, D15(1675)"                     , probestr, targetstr));
      Draw("res_nc_p_1620S31",       Form("%s + p %s, RES NC, S31(1620)"                     , probestr, targetstr));
      Draw("res_nc_p_1700D33",       Form("%s + p %s, RES NC, D33(1700)"                     , probestr, targetstr));
      Draw("res_nc_p_1440P11",       Form("%s + p %s, RES NC, P11(1440)"                     , probestr, targetstr));
      Draw("res_nc_p_1720P13",       Form("%s + p %s, RES NC, P13(1720)"                     , probestr, targetstr));
      Draw("res_nc_p_1680F15",       Form("%s + p %s, RES NC, F15(1680)"                     , probestr, targetstr));
      Draw("res_nc_p_1910P31",       Form("%s + p %s, RES NC, P31(1910)"                     , probestr, targetstr));
      Draw("res_nc_p_1920P33",       Form("%s + p %s, RES NC, P33(1920)"                     , probestr, targetstr));
      Draw("res_nc_p_1905F35",       Form("%s + p %s, RES NC, F35(1905)"                     , probestr, targetstr));
      Draw("res_nc_p_1950F37",       Form("%s + p %s, RES NC, F37(1950)"                     , probestr, targetstr));
      Draw("res_nc_p_1710P11",       Form("%s + p %s, RES NC, P11(1710)"                     , probestr, targetstr));
      Draw("dis_cc_p",               Form("%s + p %s, DIS CC"                                , probestr, targetstr));
      Draw("dis_nc_p",               Form("%s + p %s, DIS NC"                                , probestr, targetstr));
    if(gCurrProbeIsNu) {    
      Draw("dis_cc_p_ubarsea",       Form("%s + p %s, DIS CC (#bar{u}_{sea})"                , probestr, targetstr));
      Draw("dis_cc_p_dval",          Form("%s + p %s, DIS CC (d_{val})"                      , probestr, targetstr));
      Draw("dis_cc_p_dsea",          Form("%s + p %s, DIS CC (d_{sea})"                      , probestr, targetstr));
      Draw("dis_cc_p_ssea",          Form("%s + p %s, DIS CC (s_{sea})"                      , probestr, targetstr));
    }
    if(gCurrProbeIsNuBar) {    
      Draw("dis_cc_n_sbarsea",       Form("%s + n %s, DIS CC (#bar{s}_{sea})"                , probestr, targetstr));
      Draw("dis_cc_n_dbarsea",       Form("%s + n %s, DIS CC (#bar{d}_{val})"                , probestr, targetstr));
      Draw("dis_cc_n_uval",          Form("%s + n %s, DIS CC (u_{val})"                      , probestr, targetstr));
      Draw("dis_cc_n_usea",          Form("%s + n %s, DIS CC (u_{sea})"                      , probestr, targetstr));
    }
      Draw("dis_nc_p_sbarsea",       Form("%s + p %s, DIS NC (#bar{s}_{sea})"                , probestr, targetstr));
      Draw("dis_nc_p_ubarsea",       Form("%s + p %s, DIS NC (#bar{u}_{sea})"                , probestr, targetstr));
      Draw("dis_nc_p_dbarsea",       Form("%s + p %s, DIS NC (#bar{d}_{sea})"                , probestr, targetstr));
      Draw("dis_nc_p_dval",          Form("%s + p %s, DIS NC (d_{val})"                      , probestr, targetstr));
      Draw("dis_nc_p_dsea",          Form("%s + p %s, DIS NC (d_{sea})"                      , probestr, targetstr));
      Draw("dis_nc_p_uval",          Form("%s + p %s, DIS NC (u_{val})"                      , probestr, targetstr));
      Draw("dis_nc_p_usea",          Form("%s + p %s, DIS NC (u_{sea})"                      , probestr, targetstr));
      Draw("dis_nc_p_ssea",          Form("%s + p %s, DIS NC (s_{sea})"                      , probestr, targetstr));
    if(gCurrProbeIsNu) {    
      Draw("dis_cc_p_dval_charm",    Form("%s + p %s, DIS CC (d_{val} -> charm)"             , probestr, targetstr));
      Draw("dis_cc_p_dsea_charm",    Form("%s + p %s, DIS CC (d_{sea} -> charm)"             , probestr, targetstr));
      Draw("dis_cc_p_ssea_charm",    Form("%s + p %s, DIS CC (s_{sea} -> charm)"             , probestr, targetstr));
    }
    if(gCurrProbeIsNuBar) {    
      Draw("dis_cc_p_dbarsea_charm", Form("%s + p %s, DIS CC (#bar{d}_{sea} -> #bar{charm})" , probestr, targetstr));
      Draw("dis_cc_p_sbarsea_charm", Form("%s + p %s, DIS CC (#bar{s}_{sea} -> #bar{charm})" , probestr, targetstr));
    }
  }//Z>0
 
}

References DirNameToProbe(), DirNameToTarget(), Draw(), gCurrProbeIsNu, gCurrProbeIsNuBar, gCurrProbeLbl, gCurrTargetHasN, gCurrTargetHasP, gCurrTargetIsFreeNuc, gCurrTargetLbl, gDirName, LOG, OpenDir(), and pINFO.

Referenced by MakePlots().

OpenDir()

void OpenDir

(

void

)

Definition at line 219 of file gXSecComp.cxx.

{
  gDirCurr  = (TDirectory *) gXSecFile_curr->Get(gDirName.c_str());
  gDirRef0  = 0;
  if(gXSecFile_ref0) {
    gDirRef0 = (TDirectory *) gXSecFile_ref0->Get(gDirName.c_str());
  }
 
  if(!gDirRef0) {
    LOG("gxscomp", pINFO)  << "No reference plots will be shown.";
  }
}

References gDirCurr, gDirName, gDirRef0, gXSecFile_curr, gXSecFile_ref0, LOG, and pINFO.

Referenced by MakePlotsCurrDir().

OutputFileName()

string OutputFileName

(

string

input_file_name

)

Definition at line 820 of file gXSecComp.cxx.

{
// Builds the output filename based on the name of the input filename
// Perfors the following conversion: name.root -> name.nuxsec_test.ps
 
  unsigned int L = inpname.length();
 
  // if the last 4 characters are "root" (ROOT file extension) then
  // remove them
  if(inpname.substr(L-4, L).find("root") != string::npos) {
    inpname.erase(L-4, L);
  }
 
  ostringstream name;
  name << inpname << "nuxsec_test.ps";
 
  return gSystem->BaseName(name.str().c_str());
}

PrintSyntax()

void PrintSyntax

(

void

)

Definition at line 800 of file gXSecComp.cxx.

{
  LOG("gxscomp", pNOTICE)
    << "\n\n" << "Syntax:" << "\n"
    << " gxscomp  -f xsec_file [-r reference_xsec_file] [-o output]\n";
}

References LOG, and pNOTICE.

Referenced by GetCommandLineArgs().

TrimGraph()

TGraph * TrimGraph	(	TGraph *	gr,
		int	max_np_per_decade )

Definition at line 636 of file gXSecComp.cxx.

{
  if(!gr) return 0;
 
  const int np = gr->GetN();
  vector<bool> remv(np);
 
  int    fp   = 0; 
  int    lp   = 0; 
 
  while(1) {
    double xmin = gr->GetX()[fp];
    double xmax = 10 * xmin;
    int ndec = 0;
    for(int i=fp; i<np; i++) {
       remv[i] = false;
       lp = i;
       double x = gr->GetX()[i];
       if(x > xmin && x <= xmax) ndec++;
       if(x > xmax) break;
    }
    if(ndec > max_np_per_decade) {
       double keep_prob = (double)max_np_per_decade/ (double) ndec;
       int keep_rate = TMath::FloorNint(1./keep_prob);
       int ithrow = 0;
       for(int i=fp; i<=lp; i++) {
         if(ithrow % keep_rate) {
           remv[i] = true;
         }
         ithrow++;
       } 
    }
    if(lp == np-1) break;
    fp = lp+1;
  }
 
  int np2 = 0;
  for(int i=0; i<np; i++) {
   if(!remv[i]) {
     np2++;
   }
  }
 
  double * x = new double[np2];
  double * y = new double[np2];
 
  int j=0;
  for(int i=0; i<np; i++) {
   if(!remv[i]) {
     x[j] = gr->GetX()[i];
     y[j] = gr->GetY()[i];
     j++;
   }
  }
  TGraph * gr2 = new TGraph(np2,x,y);
 
  delete [] x;
  delete [] y;
 
  return gr2;
}

Referenced by Draw().

Variable Documentation

gC

TCanvas* gC = 0

Definition at line 110 of file gXSecComp.cxx.

Referenced by Draw(), End(), and Init().

gCurrProbeIsNu

bool gCurrProbeIsNu = false

Definition at line 117 of file gXSecComp.cxx.

Referenced by DirNameToProbe(), and MakePlotsCurrDir().

gCurrProbeIsNuBar

bool gCurrProbeIsNuBar = false

Definition at line 118 of file gXSecComp.cxx.

Referenced by DirNameToProbe(), and MakePlotsCurrDir().

gCurrProbeLbl

string gCurrProbeLbl = ""

Definition at line 115 of file gXSecComp.cxx.

Referenced by DirNameToProbe(), and MakePlotsCurrDir().

gCurrProbePdg

int gCurrProbePdg = 0

Definition at line 116 of file gXSecComp.cxx.

Referenced by DirNameToProbe(), and DirNameToTarget().

gCurrTargetHasN

bool gCurrTargetHasN = false

Definition at line 121 of file gXSecComp.cxx.

Referenced by DirNameToTarget(), and MakePlotsCurrDir().

gCurrTargetHasP

bool gCurrTargetHasP = false

Definition at line 120 of file gXSecComp.cxx.

Referenced by DirNameToTarget(), and MakePlotsCurrDir().

gCurrTargetIsFreeNuc

bool gCurrTargetIsFreeNuc = false

Definition at line 122 of file gXSecComp.cxx.

Referenced by DirNameToTarget(), and MakePlotsCurrDir().

gCurrTargetLbl

string gCurrTargetLbl = ""

Definition at line 119 of file gXSecComp.cxx.

Referenced by DirNameToTarget(), and MakePlotsCurrDir().

gDirCurr

TDirectory* gDirCurr = 0

Definition at line 104 of file gXSecComp.cxx.

Referenced by Draw(), and OpenDir().

gDirName

string gDirName = ""

Definition at line 108 of file gXSecComp.cxx.

Referenced by DirNameToProbe(), DirNameToTarget(), MakePlots(), MakePlotsCurrDir(), and OpenDir().

gDirRef0

TDirectory* gDirRef0 = 0

Definition at line 105 of file gXSecComp.cxx.

Referenced by Draw(), and OpenDir().

gLabelCurr

string gLabelCurr = ""

Definition at line 106 of file gXSecComp.cxx.

Referenced by Draw(), and GetCommandLineArgs().

gLabelRef0

string gLabelRef0 = ""

Definition at line 107 of file gXSecComp.cxx.

Referenced by Draw(), and GetCommandLineArgs().

gLS

TLegend* gLS = 0

Definition at line 114 of file gXSecComp.cxx.

Referenced by Draw(), End(), and Init().

gOptHaveRef

bool gOptHaveRef

Definition at line 99 of file gXSecComp.cxx.

Referenced by Draw(), GetCommandLineArgs(), Init(), and MakePlots().

gOptOutputFilename

string gOptOutputFilename = ""

Definition at line 98 of file gXSecComp.cxx.

Referenced by GetCommandLineArgs(), Init(), and main().

gOptXSecFilename_curr

string gOptXSecFilename_curr = ""

Definition at line 96 of file gXSecComp.cxx.

Referenced by GetCommandLineArgs(), Init(), and MakePlots().

gOptXSecFilename_ref0

string gOptXSecFilename_ref0 = ""

Definition at line 97 of file gXSecComp.cxx.

Referenced by GetCommandLineArgs(), Init(), and MakePlots().

gPadRatio

TPad* gPadRatio = 0

Definition at line 113 of file gXSecComp.cxx.

Referenced by Draw(), and Init().

gPadTitle

TPad* gPadTitle = 0

Definition at line 111 of file gXSecComp.cxx.

Referenced by Draw(), and Init().

gPadXSecs

TPad* gPadXSecs = 0

Definition at line 112 of file gXSecComp.cxx.

Referenced by Draw(), and Init().

gPS

TPostScript* gPS = 0

Definition at line 109 of file gXSecComp.cxx.

Referenced by Draw(), End(), and Init().

gXSecFile_curr

TFile* gXSecFile_curr = 0

Definition at line 102 of file gXSecComp.cxx.

Referenced by MakePlots(), and OpenDir().

gXSecFile_ref0

TFile* gXSecFile_ref0 = 0

Definition at line 103 of file gXSecComp.cxx.

Referenced by MakePlots(), and OpenDir().