ana::CountingExperiment Class Reference

Compare a data spectrum to MC expectation using only the event count. More...

#include "/cvmfs/nova-development.opensciencegrid.org/novasoft/releases/N20-07-07/CAFAna/Experiment/CountingExperiment.h"

Inheritance diagram for ana::CountingExperiment:

Public Member Functions
	CountingExperiment (const IPrediction *p, const Spectrum &d, const Spectrum &cosmic)
	CountingExperiment (const IPrediction *p, const Spectrum &d)
	Version without cosmics may be wanted for MC studies. More...
	~CountingExperiment ()
virtual double	ChiSq (osc::IOscCalculatorAdjustable *osc, const SystShifts &syst=SystShifts::Nominal()) const override
virtual void	SaveTo (TDirectory *dir) const override
virtual void	Derivative (osc::IOscCalculator *calc, const SystShifts &shift, std::unordered_map< const ISyst *, double > &dchi) const
virtual stan::math::var	LogLikelihood (osc::IOscCalculatorAdjustableStan *osc, const SystShifts &syst=kNoShift) const

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

Protected Attributes
const IPrediction *	fMC
Spectrum	fData
TH1 *	fCosmic

Detailed Description

Compare a data spectrum to MC expectation using only the event count.

Definition at line 12 of file CountingExperiment.h.

Constructor & Destructor Documentation

ana::CountingExperiment::CountingExperiment	(	const IPrediction *	p,
		const Spectrum &	d,
		const Spectrum &	cosmic
	)

Definition at line 17 of file CountingExperiment.cxx.

    : fMC(p), fData(d),
      fCosmic(cosmic.ToTH1(d.Livetime(), kLivetime))
  {
  }

ana::CountingExperiment::CountingExperiment ( const IPrediction *  p, const Spectrum &  d  )

inline

Version without cosmics may be wanted for MC studies.

Definition at line 17 of file CountingExperiment.h.

References ChiSq(), dir, LoadFrom(), ana::SystShifts::Nominal(), SaveTo(), plotSysts::syst, and ~CountingExperiment().

: fMC(p), fData(d), fCosmic(0) {}

ana::CountingExperiment::~CountingExperiment

(

)

Definition at line 26 of file CountingExperiment.cxx.

References fCosmic.

Referenced by CountingExperiment().

  {
    delete fCosmic;
  }

Member Function Documentation

double ana::CountingExperiment::ChiSq ( osc::IOscCalculatorAdjustable *  osc, const SystShifts &  syst = SystShifts::Nominal()  ) const

overridevirtual

Reimplemented from ana::IExperiment.

Definition at line 32 of file CountingExperiment.cxx.

References stan::math::exp(), fCosmic, fData, fMC, ana::Spectrum::Integral(), ana::LogLikelihood(), ana::Spectrum::POT(), and ana::IPrediction::PredictSyst().

Referenced by CountingExperiment().

  {
    double exp = fMC->PredictSyst(osc, syst).Integral(fData.POT());
    if (fCosmic) exp += fCosmic->Integral(0,-1);
    double obs = fData.Integral(fData.POT());

    // use namespace to disabiguate from method inherited from IExperiment
    return ana::LogLikelihood(exp, obs);
  }

virtual void ana::IExperiment::Derivative ( osc::IOscCalculator *  calc, const SystShifts &  shift, std::unordered_map< const ISyst *, double > &  dchi  ) const

inlinevirtualinherited

Reimplemented in ana::SingleSampleExperiment, ana::Ssth23Constraint, ana::AtmConstraint, ana::MultiExperiment, ana::T2KToyExperiment, ana::GaussianConstraint, ana::SolarConstraints, ana::Dmsq32Constraint, and ana::ReactorExperiment.

Definition at line 25 of file IExperiment.h.

Referenced by ana::MinuitFitter::Gradient(), and ana::MinuitFitter::SupportsDerivatives().

    {
      // Optional to implement
      //
      // NB the convention is to *add* your contribution to the dchi values.
      //
      // If unimplemented, this default will be called, signaling no result to
      // the caller
      dchi.clear();
    }

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

static

Definition at line 60 of file CountingExperiment.cxx.

References ana::assert(), ana::Spectrum::LoadFrom(), ana::LoadFrom< IPrediction >(), mc, runNovaSAM::release, and runNovaSAM::ret.

Referenced by CountingExperiment(), and ana::LoadFrom< IExperiment >().

  {
    TObjString* ptag = (TObjString*)dir->Get("type");
    assert(ptag);
    assert(ptag->GetString() == "CountingExperiment");

    assert(dir->GetDirectory("mc"));
    assert(dir->GetDirectory("data"));
    

    const IPrediction* mc = ana::LoadFrom<IPrediction>(dir->GetDirectory("mc")).release();
    const std::unique_ptr<Spectrum> data = Spectrum::LoadFrom(dir->GetDirectory("data"));

    TH1* cosmic = 0;
    if(dir->Get("cosmic")) cosmic = (TH1*)dir->Get("cosmic");

    auto ret = std::make_unique<CountingExperiment>(mc, *data);
    if(cosmic) ret->fCosmic = cosmic;
    return ret;
  }

virtual stan::math::var ana::IExperiment::LogLikelihood ( osc::IOscCalculatorAdjustableStan *  osc, const SystShifts &  syst = kNoShift  ) const

inlinevirtualinherited

Reimplemented in test::GaussQuadExperiment.

Definition at line 38 of file IExperiment.h.

References ana::assert(), dir, and ana::IExperiment::SaveTo().

Referenced by ana::StanFitter::log_prob().

      {
        assert(false && "unimplemented");
        return 0;
      };

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

overridevirtual

Reimplemented from ana::IExperiment.

Definition at line 44 of file CountingExperiment.cxx.

References dir, fCosmic, fData, fMC, ana::IPrediction::SaveTo(), ana::Spectrum::SaveTo(), and tmp.

Referenced by CountingExperiment().

  {
    TDirectory* tmp = dir;

    dir->cd();
    TObjString("CountingExperiment").Write("type");

    fMC->SaveTo(dir->mkdir("mc"));
    fData.SaveTo(dir->mkdir("data"));
    
    if(fCosmic) fCosmic->Write("cosmic");

    tmp->cd();
  }

Member Data Documentation

TH1* ana::CountingExperiment::fCosmic

protected

Definition at line 27 of file CountingExperiment.h.

Referenced by ChiSq(), SaveTo(), and ~CountingExperiment().

Spectrum ana::CountingExperiment::fData

protected

Definition at line 26 of file CountingExperiment.h.

Referenced by ChiSq(), and SaveTo().

const IPrediction* ana::CountingExperiment::fMC

protected

Definition at line 25 of file CountingExperiment.h.

Referenced by ChiSq(), and SaveTo().

---

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

• /cvmfs/nova-development.opensciencegrid.org/novasoft/releases/N20-07-07/CAFAna/Experiment/CountingExperiment.h
• /cvmfs/nova-development.opensciencegrid.org/novasoft/releases/N20-07-07/CAFAna/Experiment/CountingExperiment.cxx