NuMultiRunXFitAnalysis.cxx

Go to the documentation of this file.

// $Id: NuMultiRunXFitAnalysis.cxx,v 1.1 2007/12/27 22:13:46 evans Exp $
//
// Wrapper to perform the cross section fitting extrapolation and
// analyses for multiple runs.
//
// Justin Evans
// j.evans2@physics.ox.ac.uk
//
//
// $Log: NuMultiRunXFitAnalysis.cxx,v $
// Revision 1.1  2007/12/27 22:13:46  evans
// New class to allow the XFit extrapolation to be performed on multiple
// runs (e.g. SKZP periods), along with necessary interface additions to
// NuXFitAnalysis.
//
//

#include "MessageService/MsgService.h"
#include "NtupleUtils/NuMultiRunXFitAnalysis.h"

ClassImp(NuMultiRunXFitAnalysis)

CVSID("$Id: NuMultiRunXFitAnalysis.cxx,v 1.1 2007/12/27 22:13:46 evans Exp $");

//____________________________________________________________________72
NuMultiRunXFitAnalysis::NuMultiRunXFitAnalysis
(vector<NuFDXFitConfig>& vFDConfig)
{
  fanalysisSetting = NuXFit::kUnknown;

  fvFDConfig = vFDConfig;

  fFitter = new TFitterMinuit();
  fFitter->CreateMinimizer();
  fFitter->SetMinuitFCN(this);
}

//____________________________________________________________________72
double NuMultiRunXFitAnalysis
::operator () (const vector<double>& pars) const
{
  static Int_t fitcounter = 0;
  if (!(fitcounter%50)){
    cout << "Minuit FD XFit call " << fitcounter << endl;
  }
  ++fitcounter;
  if (NuXFit::kCPT == fanalysisSetting){
    Double_t like = 0;
    for (vector<NuXFitAnalysis*>::const_iterator it = fvExtrapolators.begin();
         it != fvExtrapolators.end();
         ++it){
      (*it)->FillCPTFDPrediction(pars);
      like += (*it)->CalculateLikelihood();
    }
    return like;
  }
  else if(NuXFit::kTransition == fanalysisSetting){
    Double_t like = 0;
    for (vector<NuXFitAnalysis*>::const_iterator it = fvExtrapolators.begin();
         it != fvExtrapolators.end();
         ++it){
      (*it)->FillTransitionFDPrediction(pars);
      like += (*it)->CalculateLikelihood();
    }
    return like;
  }
  else{
    MSG("NuXFitAnalysis.cxx",Msg::kFatal)
      << "Incorrect analysis setting" << endl;
    return -1.0;
  }
}

//____________________________________________________________________72
NuMultiRunXFitAnalysis::~NuMultiRunXFitAnalysis()
{
}

//____________________________________________________________________72
void NuMultiRunXFitAnalysis::PerformFit()
{
  this->ConfigureExtrapolators();
  fFitter->Minimize();
  fFitter->PrintResults(0,0);
}

//____________________________________________________________________72
void NuMultiRunXFitAnalysis::ConfigureExtrapolators()
{
  for (vector<NuFDXFitConfig>::iterator it = fvFDConfig.begin();
       it != fvFDConfig.end();
       ++it){
    NuXFitAnalysis* xAnal = new NuXFitAnalysis((*it).xmlConfig);
    xAnal->WriteOutput(false);
    xAnal->NDFitResults((*it).ndFitResults);
    xAnal->NDFitBinningNuMuCC((*it).ndFitBinsNuMuCC);
    xAnal->NDFitBinningNuMuBarCC((*it).ndFitBinsNuMuBarCC);
    xAnal->FDNuMuCCData(*(*it).fdNuMuCCData);
    xAnal->FDNuMuBarCCData(*(*it).fdNuMuBarCCData);
    xAnal->FDDataPoT((*it).fdDataPoT);
    xAnal->FDMC((*it).fdMCFileName);
    xAnal->ConfigureForExternalFit((*it).analSetting);
    fvExtrapolators.push_back(xAnal);

    fanalysisSetting = (*it).analSetting;
  }
}

//____________________________________________________________________72
void NuMultiRunXFitAnalysis::SetupPars()
{
  if (NuXFit::kCPT == fanalysisSetting){
    this->SetupCPTPars();
    return;
  }
  else if (NuXFit::kTransition == fanalysisSetting){
    this->SetupTransitionPars();
    return;
  }
  else{
    MSG("NuXFitAnalysis.cxx",Msg::kFatal)
      << "Incorrect analysis setting" << endl;
    return;
  }
}

//____________________________________________________________________72
void NuMultiRunXFitAnalysis::SetupCPTPars()
{
  fFitter->SetParameter(0,
                        "Dm2",
                        0.003,0.2,1.0,0.0);
  fFitter->SetParameter(1,
                        "Sn2",
                        1.0,0.2,1.0,0.0);
  fFitter->SetParameter(2,
                        "Dm2Bar",
                        0.003,0.2,1.0,0.0);
  fFitter->SetParameter(3,
                        "Sn2Bar",
                        1.0,0.2,1.0,0.0);
}

//____________________________________________________________________72
void NuMultiRunXFitAnalysis::SetupTransitionPars()
{
  //Setup fit parameters
  fFitter->SetParameter(0,
                        "transitionProb",
                        0.0,0.2,1.0,0.0);
}

---