void MakeDaikonHelper(Int_t whichPID) 
{
  gROOT->LoadMacro("NueAna/macros/LoadLibs.C"); LoadLibs();
  gSystem->Load("libNueAnaExtrapolation.so");

  MCReweight *mcr = &MCReweight::Instance();
  NeugenWeightCalculator *n = new NeugenWeightCalculator();
  mcr->AddWeightCalculator(n);

  TChain *nearChain = new TChain("ana_nue");
  TChain *farChain = new TChain("ana_nue");

  nearChain->Add("AnaNue_Daikon/AnaNue-Birch-n130*.root");
  farChain->Add("AnaNue_Daikon/AnaNue-Birch-f21*.root");
  Double_t farPOT  = 10.*1.02716e20.;
  Double_t nearPOT = 24*400*2.568e13;
  
  Selection::Selection_t sel = whichPID;
  Systematic::Systematic_t sys = Systematic::kShwDev;

  //make plot maker helper:
  NueFNHelper help(200,0.,100.);
  help.SetChains(nearChain,farChain,nearPOT,farPOT);
  
  //add a standard entry:  (numu CC osc pars, no nue oscillations; SKZP out to 30GeV)
  NueSystematic *nueSys = new NueSystematic("Standard");
  nueSys->AddSystematic(Systematic::kSKZP,Systematic::GetDefaultValue(Systematic::kSKZP));
  nueSys->AddSystematic(Systematic::kOscProb,Systematic::GetDefaultValue(Systematic::kOscProb));
  nueSys->SetOscParams(0.554,0.7854,0.,8.2e-5,2.7e-3,0,1);
  if(sys!=Systematic::kSKZP && sys!=Systematic::kOscProb) 
    nueSys->AddSystematic(sys,Systematic::GetDefaultValue(sys));  
  help.AddNueSystematic(nueSys);
  //run all the plot making code:
  help.MakeHelpers(sel);
  //write out the file:
  string outfiletag = string(Selection::AsString(sel)) + 
    "_Alter" + string(Systematic::AsString(sys));
  help.WriteFile(outfiletag.c_str());
}