NuMMRunNuBar.cxx

Go to the documentation of this file.

#include <cmath>
#include <cassert>


#include "Rtypes.h"
#include "TH1D.h"

#include "Conventions/Detector.h"
#include "Conventions/SimFlag.h"
#include "MessageService/MsgService.h"
#include "NtupleUtils/NuCuts.h"
#include "NtupleUtils/NuHistInterpolator.h"
#include "NtupleUtils/NuMMHelperPRL.h"
#include "NtupleUtils/NuMMRunNuBar.h"
#include "NtupleUtils/NuMatrixSpectrum.h"

ClassImp(NuMMRunNuBar)

CVSID("$Id: NuMMRunNuBar.cxx,v 1.5 2008/10/28 22:59:29 ahimmel Exp $");

int NuMMRunNuBar::counter = 0;

//____________________________________________________________________72
NuMMRunNuBar::NuMMRunNuBar()
  : NuMMRun()
{
  fCached = false;
  
  fndNuData_orig = 0;
  fndBarData_orig = 0;
}

//____________________________________________________________________72
NuMMRunNuBar::NuMMRunNuBar(NuMMHelperPRL* helper,
                 NuHistInterpolator* ndData,
                 NuHistInterpolator* fdData)
  : NuMMRun(helper, ndData, fdData)
{
  fCached = false;
  
  fndNuData_orig = 0;
  fndBarData_orig = 0;
}


void NuMMRunNuBar::MakeFDPredNoPair(const NuMMParameters& pars, NuMatrixSpectrum **nu, NuMatrixSpectrum **nubar)
{
  pair<NuMatrixSpectrum, NuMatrixSpectrum> fits = MakeFDPred(pars);
  if (nu) *nu = new NuMatrixSpectrum(fits.first);
  if (nubar) *nubar = new NuMatrixSpectrum(fits.second);
}


NuMatrixSpectrum* NuMMRunNuBar::MakeFDPredNuMu(const NuMMParameters& pars)
{
  pair<NuMatrixSpectrum, NuMatrixSpectrum> fits = MakeFDPred(pars);
  return new NuMatrixSpectrum(fits.first);
}

NuMatrixSpectrum* NuMMRunNuBar::MakeFDPredNuBar(const NuMMParameters& pars)
{
  pair<NuMatrixSpectrum, NuMatrixSpectrum> fits = MakeFDPred(pars);
  return new NuMatrixSpectrum(fits.second);
}


NuMatrixSpectrum* NuMMRunNuBar::ModifyNDNuData()
{
  if (fndNuData_orig) {
    MAXMSG("NuMMRunNuBar",Msg::kWarning, 5) << "You did not reset the Nu data after modifying it.  Data should always be reset.  Doing so now." << endl;
    ResetNDNuData();
  }

  if (!fQuietMode) {MAXMSG("NuMMRunNuBar",Msg::kInfo, 5) << "Making a Nu backup copy." << endl;}
  
  // Backup the original object, create a new one for modifying
  fndNuData_orig = fndNuData;
  fndNuData = new NuMatrixSpectrum(*fndNuData);
  
  // Redo caching since the ND Data has changed
  fCached = false;
  
  return fndNuData;
}


NuMatrixSpectrum* NuMMRunNuBar::ModifyNDBarData()
{
  if (fndBarData_orig) {
    MAXMSG("NuMMRunNuBar",Msg::kWarning, 5) << "You did not reset the Bar data after modifying it.  Data should always be reset.  Doing so now." << endl;
    ResetNDBarData();
  }

  if (!fQuietMode) {MAXMSG("NuMMRunNuBar",Msg::kInfo, 5) << "Making a Bar backup copy." << endl;}
  
  // Backup the original object, create a new one for modifying
  fndBarData_orig = fndBarData;
  fndBarData = new NuMatrixSpectrum(*fndBarData);
  
  // Redo caching since the ND Data has changed
  fCached = false;
  return fndBarData;
}


void NuMMRunNuBar::ResetNDNuData()
{
  if (fndNuData_orig) {
    if (!fQuietMode) {MAXMSG("NuMMRunNuBar",Msg::kInfo, 5) << "Reverting Nu to backup copy." << endl;}
    
    delete fndNuData;
    fndNuData = fndNuData_orig;
    fndNuData_orig = 0;
    
    // Redo caching since the ND Data has changed
    fCached = false;
  }
  else {
    MAXMSG("NuMMRunNuBar",Msg::kWarning, 5) << "Resetting unmodified Nu data.  Not harmful, just superfluous." << endl;
  }
}

void NuMMRunNuBar::ResetNDBarData()
{
  if (fndBarData_orig) {
    if (!fQuietMode) {MAXMSG("NuMMRunNuBar",Msg::kInfo, 5) << "Reverting Bar to backup copy." << endl;}
    
    delete fndBarData;
    fndBarData = fndBarData_orig;
    fndBarData_orig = 0;
    
    // Redo caching since the ND Data has changed
    fCached = false;
  }
  else {
    MAXMSG("NuMMRunNuBar",Msg::kWarning, 5) << "Resetting unmodified Bars data.  Not harmful, just superfluous." << endl;
  }
}

NuMMHelperCPT* NuMMRunNuBar::GetHelper()
{
  fCached = false;
  return fHelper;
}

---