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/N21-02-25/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::IOscCalcAdjustable *osc, const SystShifts &syst=SystShifts::Nominal()) const override
virtual void	SaveTo (TDirectory *dir, const std::string &name) const override
virtual stan::math::var	LogLikelihood (osc::IOscCalcAdjustableStan *osc, const SystShifts &syst=kNoShift) const

Static Public Member Functions
static std::unique_ptr< CountingExperiment >	LoadFrom (TDirectory *dir, const std::string &name)

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 20 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(), string, and ~CountingExperiment().

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

ana::CountingExperiment::~CountingExperiment

(

)

Definition at line 29 of file CountingExperiment.cxx.

References fCosmic.

Referenced by CountingExperiment().

  {
    delete fCosmic;
  }

Member Function Documentation

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

overridevirtual

Reimplemented from ana::IExperiment.

Definition at line 35 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);
  }

std::unique_ptr< CountingExperiment > ana::CountingExperiment::LoadFrom ( TDirectory *  dir, const std::string &  name  )

static

Definition at line 68 of file CountingExperiment.cxx.

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

Referenced by CountingExperiment().

  {
    dir = dir->GetDirectory(name.c_str()); // switch to subdir
    assert(dir);

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

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

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

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

    delete dir;

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

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

inlinevirtualinherited

Reimplemented in test::GaussQuadExperiment, ana::SingleSampleExperiment, and ana::MultiExperiment.

Definition at line 25 of file IExperiment.h.

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

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

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

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

overridevirtual

Reimplemented from ana::IExperiment.

Definition at line 47 of file CountingExperiment.cxx.

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

Referenced by CountingExperiment().

  {
    TDirectory* tmp = dir;

    dir = dir->mkdir(name.c_str()); // switch to subdir
    dir->cd();

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

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

    dir->Write();
    delete dir;

    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/N21-02-25/CAFAna/Experiment/CountingExperiment.h
• /cvmfs/nova-development.opensciencegrid.org/novasoft/releases/N21-02-25/CAFAna/Experiment/CountingExperiment.cxx