ana::ModularExtrapComponent Class Reference

Base class for component extrapolation. More...

#include "/cvmfs/nova-development.opensciencegrid.org/novasoft/releases/N21-04-19/CAFAna/Extrap/ModularExtrapComponent.h"

Inheritance diagram for ana::ModularExtrapComponent:

Classes
class	DivByZeroCounter

Public Member Functions
	ModularExtrapComponent ()
virtual	~ModularExtrapComponent ()
virtual void	SaveTo (TDirectory *dir, const std::string &name) const =0
virtual void	SavePlots (TDirectory *dir, double potFD) const =0
const OscillatableSpectrum &	Return () const
	Interface so that result of Eval() is called only once and cached. More...
virtual const IDecomp *	GetDecomp () const =0

Static Public Member Functions
static void	SetQuiet (bool quiet=true)

Protected Member Functions
virtual OscillatableSpectrum	Eval () const =0
	Core extrapolation math. More...

Static Protected Member Functions
static Spectrum	GetDecompResult (const IDecomp &, DecompResult)
	Helper function to pick out single Spectrum from a decomposition. More...
static std::string	DRToString (DecompResult)
	Helper function to turn a DecompResult into a string (for storage). More...
static DecompResult	StringToDR (std::string)
	Helper function to turn a string into a DecompResult (for loading). More...
static Ratio	FormSmartRatio (const Spectrum &num, const Spectrum &denom, std::string component, std::string location, const Spectrum &mult)
	Form Ratio, but be aware of zero division. More...
static void	ComparisonPlot (Spectrum mc, Spectrum notMC, double pot, std::string notMCLabel, int notMCColor, std::string latex, std::string title, std::string saveAs, bool restrictRange=false)

Private Attributes
bool	fEvaluated
OscillatableSpectrum	fCache

Static Private Attributes
static bool	fQuiet = false

Detailed Description

Base class for component extrapolation.

Definition at line 21 of file ModularExtrapComponent.h.

Constructor & Destructor Documentation

ana::ModularExtrapComponent::ModularExtrapComponent ( )

inline

Definition at line 25 of file ModularExtrapComponent.h.

: fEvaluated(false), fCache(OscillatableSpectrum::Uninitialized()) {}

virtual ana::ModularExtrapComponent::~ModularExtrapComponent ( )

inlinevirtual

Definition at line 26 of file ModularExtrapComponent.h.

References dir, ana::SaveTo(), and string.

{};

Member Function Documentation

void ana::ModularExtrapComponent::ComparisonPlot ( Spectrum  mc, Spectrum  notMC, double  pot, std::string  notMCLabel, int  notMCColor, std::string  latex, std::string  title, std::string  saveAs, bool  restrictRange = false  )

staticprotected

Definition at line 17 of file ModularExtrapComponentPlot.cxx.

References plot_validation_datamc::c, kRed, make_mec_shifts_plots::legend, and ana::Spectrum::ToTH1().

Referenced by ana::FluxReweight::SavePlots(), ana::TruthReweight::SavePlots(), and ana::RecoReweight::SavePlots().

   {
    TH1* histMC( mc.ToTH1(pot) );
    TH1* histNotMC( notMC.ToTH1(pot) );
    histMC->SetLineColor( kRed );
    histNotMC->SetLineColor( notMCColor );
    TCanvas c;
    histMC->SetMaximum( 1.5 * histMC->GetMaximum() );
    if (restrictRange) histMC->GetXaxis()->SetRangeUser( 0.5, 3.5);
    histMC->SetTitle( (latex+" "+title).c_str() );
    histMC->Draw("hist");
    histNotMC->Draw("same");
    TLegend legend( .15, .7, .4, .85 );
    legend.SetTextSize(0.05);
    legend.AddEntry( histNotMC, notMCLabel.c_str() );
    legend.AddEntry( histMC, "MC" );
    legend.SetFillStyle(0);
    legend.Draw();
    histMC->Write( (TString)"MC"+saveAs.c_str());
    histNotMC->Write( (TString)notMCLabel.c_str() + saveAs.c_str());

    c.Write( saveAs.c_str() );
    delete histMC;
    delete histNotMC;
  }

std::string ana::ModularExtrapComponent::DRToString ( DecompResult  dr )

staticprotected

Helper function to turn a DecompResult into a string (for storage).

Definition at line 98 of file ModularExtrapComponent.cxx.

References ana::assert(), ana::NC, ana::NCbar, ana::NCtot, ana::nue, ana::nuebar, ana::numu, and ana::numubar.

Referenced by ana::TruthReweight::SaveTo(), and ana::RecoReweight::SaveTo().

   {
    switch (dr){
      case DecompResult::nue     : return "Nue";
      case DecompResult::numu    : return "Numu";
      case DecompResult::nuebar  : return "NueBar";
      case DecompResult::numubar : return "NumuBar";
      case DecompResult::NC      : return "NC";
      case DecompResult::NCbar      : return "NCbar";
      case DecompResult::NCtot      : return "NCtot";
    }
    assert( 0 && "Bad DecompResult" );
  }

virtual OscillatableSpectrum ana::ModularExtrapComponent::Eval ( ) const

protectedpure virtual

Core extrapolation math.

Implemented in ana::RecoReweight, ana::TruthReweight, ana::NoReweight, and ana::FluxReweight.

Referenced by Return().

Ratio ana::ModularExtrapComponent::FormSmartRatio ( const Spectrum &  num, const Spectrum &  denom, std::string  component, std::string  location, const Spectrum &  mult  )

staticprotected

Form Ratio, but be aware of zero division.

Necessary because root thinks n/0==0, so we lose events when low stats cause empty bins. If zero division occurs, ratio is set to 1 and warning is issued. If numerator is 0 and ratio will be multiplied by 0 later, the warning is suppressed.

Definition at line 130 of file ModularExtrapComponent.cxx.

References ana::assert(), bin, hadd_many_files::counter, ana::ModularExtrapComponent::DivByZeroCounter::fBins, fQuiet, ana::Spectrum::GetBinnings(), ana::Spectrum::GetEigen(), and ana::Spectrum::GetLabels().

Referenced by ana::FluxReweight::Eval(), ana::TruthReweight::Eval(), and ana::RecoReweight::Eval().

   {
    const Eigen::ArrayXd numa = num.GetEigen(1e20);
    const Eigen::ArrayXd denoma = denom.GetEigen(1e20);
    const Eigen::ArrayXd multa = mult.GetEigen(1e20);

    Eigen::ArrayXd ratioa = numa;
    ratioa.setZero();

    assert( (numa.size() == denoma.size()) && "Bin Mismatch" );
    assert( (numa.size() == multa.size()) && "Bin Mismatch" );

    static DivByZeroCounter counter;
    for (int bin = 0; bin < numa.size(); ++bin)
    {
      if ( denoma[bin] != 0 ){
        ratioa[bin] = numa[bin] / denoma[bin];
      } else {
        ratioa[bin] = 1;
        if (    numa[bin] != 0
             || multa[bin] != 0 )
          counter.fBins.insert(std::make_tuple(component, location, bin));
      }
    }

    return Ratio( std::move(ratioa), num.GetLabels(), num.GetBinnings() );
  }

virtual const IDecomp* ana::ModularExtrapComponent::GetDecomp ( ) const

pure virtual

Implemented in ana::RecoReweight, ana::TruthReweight, ana::NoReweight, and ana::FluxReweight.

Spectrum ana::ModularExtrapComponent::GetDecompResult ( const IDecomp &  decomp, DecompResult  dr  )

staticprotected

Helper function to pick out single Spectrum from a decomposition.

Definition at line 80 of file ModularExtrapComponent.cxx.

References ana::IDecomp::AntiNueComponent(), ana::IDecomp::AntiNumuComponent(), ana::assert(), ana::NC, ana::IDecomp::NCAntiComponent(), ana::NCbar, ana::IDecomp::NCComponent(), ana::NCtot, ana::IDecomp::NCTotalComponent(), ana::nue, ana::nuebar, ana::IDecomp::NueComponent(), ana::numu, ana::numubar, and ana::IDecomp::NumuComponent().

Referenced by ana::TruthReweight::Eval(), ana::RecoReweight::Eval(), ana::TruthReweight::SavePlots(), and ana::RecoReweight::SavePlots().

   {
    switch (dr){
      case DecompResult::nue     : return decomp.NueComponent();
      case DecompResult::numu    : return decomp.NumuComponent();
      case DecompResult::nuebar  : return decomp.AntiNueComponent();
      case DecompResult::numubar : return decomp.AntiNumuComponent();
      case DecompResult::NCtot      : return decomp.NCTotalComponent();
      case DecompResult::NCbar      : return decomp.NCAntiComponent();
      case DecompResult::NC      : return decomp.NCComponent();
    }
    assert( 0 && "Bad DecompResult" );
  }

const OscillatableSpectrum & ana::ModularExtrapComponent::Return

(

)

const

Interface so that result of Eval() is called only once and cached.

Definition at line 68 of file ModularExtrapComponent.cxx.

References Eval(), fCache, and fEvaluated.

Referenced by ana::FluxReweight::SavePlots(), ana::TruthReweight::SavePlots(), and ana::RecoReweight::SavePlots().

  {
    if (!fEvaluated)
    {
      fCache = Eval();
      fEvaluated = true;
    }
    return fCache;
  }

virtual void ana::ModularExtrapComponent::SavePlots ( TDirectory *  dir, double  potFD  ) const

pure virtual

Implemented in ana::RecoReweight, ana::TruthReweight, ana::NoReweight, and ana::FluxReweight.

virtual void ana::ModularExtrapComponent::SaveTo ( TDirectory *  dir, const std::string &  name  ) const

pure virtual

Implemented in ana::RecoReweight, ana::TruthReweight, and ana::NoReweight.

static void ana::ModularExtrapComponent::SetQuiet ( bool  quiet = true )

inlinestatic

Definition at line 29 of file ModularExtrapComponent.h.

References ana::GetDecomp(), latex(), location, mc, num, pot, string, and plotROC::title.

{ fQuiet = quiet; }

DecompResult ana::ModularExtrapComponent::StringToDR ( std::string  str )

staticprotected

Helper function to turn a string into a DecompResult (for loading).

Definition at line 115 of file ModularExtrapComponent.cxx.

References ana::assert(), ana::NC, ana::NCbar, ana::NCtot, ana::nue, ana::nuebar, ana::numu, and ana::numubar.

Referenced by ana::TruthReweight::LoadFrom(), and ana::RecoReweight::LoadFrom().

   {
    if      (str=="Nue")     return DecompResult::nue;
    else if (str=="Numu")    return DecompResult::numu;
    else if (str=="NueBar")  return DecompResult::nuebar;
    else if (str=="NumuBar") return DecompResult::numubar;
    else if (str=="NC")      return DecompResult::NC;
    else if (str=="NCbar")      return DecompResult::NCbar;
    else if (str=="NCtot")      return DecompResult::NCtot;
    else assert( 0 && "Bad DecompResult String" );
  }

Member Data Documentation

OscillatableSpectrum ana::ModularExtrapComponent::fCache

mutableprivate

Definition at line 77 of file ModularExtrapComponent.h.

Referenced by Return().

bool ana::ModularExtrapComponent::fEvaluated

mutableprivate

Definition at line 76 of file ModularExtrapComponent.h.

Referenced by Return().

bool ana::ModularExtrapComponent::fQuiet = false

staticprivate

Definition at line 78 of file ModularExtrapComponent.h.

Referenced by FormSmartRatio().

---

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

• /cvmfs/nova-development.opensciencegrid.org/novasoft/releases/N21-04-19/CAFAna/Extrap/ModularExtrapComponent.h
• /cvmfs/nova-development.opensciencegrid.org/novasoft/releases/N21-04-19/CAFAna/Extrap/ModularExtrapComponent.cxx
• /cvmfs/nova-development.opensciencegrid.org/novasoft/releases/N21-04-19/CAFAna/Extrap/ModularExtrapComponentPlot.cxx