ana::CountingExperiment Class Reference

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

#include <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 override
Public Member Functions inherited from ana::IExperiment
virtual	~IExperiment ()

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.

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

ana::CountingExperiment::~CountingExperiment

(

)

Definition at line 26 of file CountingExperiment.cxx.

  {
    delete fCosmic;
  }

Member Function Documentation

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

overridevirtual

Implements ana::IExperiment.

Definition at line 32 of file CountingExperiment.cxx.

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

    return LogLikelihood(exp, obs);
  }

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

static

Definition at line 59 of file CountingExperiment.cxx.

  {
    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;
  }

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

overridevirtual

Reimplemented from ana::IExperiment.

Definition at line 43 of file CountingExperiment.cxx.

  {
    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.

Spectrum ana::CountingExperiment::fData

protected

Definition at line 26 of file CountingExperiment.h.

const IPrediction* ana::CountingExperiment::fMC

protected

Definition at line 25 of file CountingExperiment.h.

---

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

• CountingExperiment.h
• CountingExperiment.cxx