ana::PredictionExtrap Class Reference

Take the output of an extrapolation and oscillate it as required. More...

#include <PredictionExtrap.h>

Inheritance diagram for ana::PredictionExtrap:

Public Member Functions
	PredictionExtrap (IExtrap *extrap)
	Takes ownership of extrap. More...
virtual	~PredictionExtrap ()
virtual Spectrum	Predict (osc::IOscCalc *calc) const override
virtual Spectrum	PredictComponent (osc::IOscCalc *calc, Flavors::Flavors_t flav, Current::Current_t curr, Sign::Sign_t sign) const override
OscillatableSpectrum	ComponentCC (int from, int to) const override
virtual void	SaveTo (TDirectory *dir) const override
	PredictionExtrap ()=delete
IExtrap *	GetExtrap () const
Public Member Functions inherited from ana::IPrediction
virtual	~IPrediction ()
virtual Spectrum	PredictUnoscillated () const
virtual Spectrum	PredictSyst (osc::IOscCalc *calc, const SystShifts &syst) const
virtual Spectrum	PredictComponentSyst (osc::IOscCalc *calc, const SystShifts &syst, Flavors::Flavors_t flav, Current::Current_t curr, Sign::Sign_t sign) const

Static Public Member Functions
static std::unique_ptr < PredictionExtrap >	LoadFrom (TDirectory *dir)

Protected Attributes
IExtrap *	fExtrap

Detailed Description

Take the output of an extrapolation and oscillate it as required.

Definition at line 10 of file PredictionExtrap.h.

Constructor & Destructor Documentation

ana::PredictionExtrap::PredictionExtrap

(

IExtrap *

extrap

)

Takes ownership of extrap.

Definition at line 13 of file PredictionExtrap.cxx.

    : fExtrap(extrap)
  {
  }

ana::PredictionExtrap::~PredictionExtrap ( )

virtual

Definition at line 19 of file PredictionExtrap.cxx.

  {
    //    delete fExtrap;
  }

ana::PredictionExtrap::PredictionExtrap ( )

delete

Member Function Documentation

OscillatableSpectrum ana::PredictionExtrap::ComponentCC ( int  from, int  to  ) const

overridevirtual

Reimplemented from ana::IPrediction.

Definition at line 76 of file PredictionExtrap.cxx.

  {
    if(from == +12 && to == +12) return fExtrap->NueSurvComponent();
    if(from == -12 && to == -12) return fExtrap->AntiNueSurvComponent();

    if(from == +12 && to == +14) return fExtrap->NumuAppComponent();
    if(from == -12 && to == -14) return fExtrap->AntiNumuAppComponent();

    if(from == +12 && to == +16) return fExtrap->TauFromEComponent();
    if(from == -12 && to == -16) return fExtrap->AntiTauFromEComponent();

    if(from == +14 && to == +12) return fExtrap->NueAppComponent();
    if(from == -14 && to == -12) return fExtrap->AntiNueAppComponent();

    if(from == +14 && to == +14) return fExtrap->NumuSurvComponent();
    if(from == -14 && to == -14) return fExtrap->AntiNumuSurvComponent();

    if(from == +14 && to == +16) return fExtrap->TauFromMuComponent();
    if(from == -14 && to == -16) return fExtrap->AntiTauFromMuComponent();

    assert(0 && "Not reached");
  }

IExtrap* ana::PredictionExtrap::GetExtrap ( ) const

inline

Definition at line 32 of file PredictionExtrap.h.

{return fExtrap;}

std::unique_ptr< PredictionExtrap > ana::PredictionExtrap::LoadFrom ( TDirectory *  dir )

static

Definition at line 120 of file PredictionExtrap.cxx.

  {
    assert(dir->GetDirectory("extrap"));
    IExtrap* extrap = ana::LoadFrom<IExtrap>(dir->GetDirectory("extrap")).release();

    return std::unique_ptr<PredictionExtrap>(new PredictionExtrap(extrap));
  }

Spectrum ana::PredictionExtrap::Predict ( osc::IOscCalc *  calc ) const

overridevirtual

Implements ana::IPrediction.

Definition at line 25 of file PredictionExtrap.cxx.

  {
    return PredictComponent(calc,
                            Flavors::kAll,
                            Current::kBoth,
                            Sign::kBoth);
  }

Spectrum ana::PredictionExtrap::PredictComponent ( osc::IOscCalc *  calc, Flavors::Flavors_t  flav, Current::Current_t  curr, Sign::Sign_t  sign  ) const

overridevirtual

Implements ana::IPrediction.

Definition at line 34 of file PredictionExtrap.cxx.

  {
    // Get binning
    Spectrum ret = fExtrap->NCComponentFromNumu().Oscillated(calc, 14, 0);
    ret.Clear();

    if(curr & Current::kCC){
      if(flav & Flavors::kNuEToNuE    && sign & Sign::kNu)     {
        ret += fExtrap->NueSurvComponent().    Oscillated(calc, +12, +12);
      }
      if(flav & Flavors::kNuEToNuE    && sign & Sign::kAntiNu) ret += fExtrap->AntiNueSurvComponent().Oscillated(calc, -12, -12);

      if(flav & Flavors::kNuEToNuMu   && sign & Sign::kNu)     ret += fExtrap->NumuAppComponent().    Oscillated(calc, +12, +14);
      if(flav & Flavors::kNuEToNuMu   && sign & Sign::kAntiNu) ret += fExtrap->AntiNumuAppComponent().Oscillated(calc, -12, -14);

      if(flav & Flavors::kNuEToNuTau  && sign & Sign::kNu)     ret += fExtrap->TauFromEComponent().    Oscillated(calc, +12, +16);
      if(flav & Flavors::kNuEToNuTau  && sign & Sign::kAntiNu) ret += fExtrap->AntiTauFromEComponent().Oscillated(calc, -12, -16);

      if(flav & Flavors::kNuMuToNuE   && sign & Sign::kNu)     ret += fExtrap->NueAppComponent().    Oscillated(calc, +14, +12);
      if(flav & Flavors::kNuMuToNuE   && sign & Sign::kAntiNu) ret += fExtrap->AntiNueAppComponent().Oscillated(calc, -14, -12);

      if(flav & Flavors::kNuMuToNuMu  && sign & Sign::kNu)     ret += fExtrap->NumuSurvComponent().    Oscillated(calc, +14, +14);
      if(flav & Flavors::kNuMuToNuMu  && sign & Sign::kAntiNu) ret += fExtrap->AntiNumuSurvComponent().Oscillated(calc, -14, -14);

      if(flav & Flavors::kNuMuToNuTau && sign & Sign::kNu)     ret += fExtrap->TauFromMuComponent().    Oscillated(calc, +14, +16);
      if(flav & Flavors::kNuMuToNuTau && sign & Sign::kAntiNu) ret += fExtrap->AntiTauFromMuComponent().Oscillated(calc, -14, -16);
    }
    if(curr & Current::kNC){
      assert(flav == Flavors::kAll); // Don't know how to calculate anything else
      assert(sign == Sign::kBoth);   // Why would you want to split NCs out by sign?

      ret += fExtrap->NCComponentFromNumu().Oscillated(calc, 14, 0);
      ret += fExtrap->NCComponentFromNue().Oscillated(calc, 12, 0);
    }

    return ret;
  }

void ana::PredictionExtrap::SaveTo ( TDirectory *  dir ) const

overridevirtual

Reimplemented from ana::IPrediction.

Reimplemented in ana::PredictionNoExtrap.

Definition at line 106 of file PredictionExtrap.cxx.

  {
    TDirectory* tmp = gDirectory;

    dir->cd();

    TObjString("PredictionExtrap").Write("type");

    fExtrap->SaveTo(dir->mkdir("extrap"));

    tmp->cd();
  }

Member Data Documentation

IExtrap* ana::PredictionExtrap::fExtrap

protected

Definition at line 35 of file PredictionExtrap.h.

---

The documentation for this class was generated from the following files:

• PredictionExtrap.h
• PredictionExtrap.cxx