ana::SingleSampleExperiment Class Reference

Compare a single data spectrum to the MC + cosmics expectation. More...

#include <SingleSampleExperiment.h>

Inheritance diagram for ana::SingleSampleExperiment:

Public Member Functions
	SingleSampleExperiment (const IPrediction *pred, const Spectrum &data, const Spectrum &cosmicInTime, const Spectrum &cosmicOutOfTime)
	SingleSampleExperiment (const IPrediction *pred, const Spectrum &data)
	In MC studies you might not want to bother with cosmics. More...
	fCosmicOutOfTime (0,{},{}, 0, 0)
	fMask (0)
virtual	~SingleSampleExperiment ()
virtual double	ChiSq (osc::IOscCalcAdjustable *osc, const SystShifts &syst=SystShifts::Nominal()) const override
virtual void	SaveTo (TDirectory *dir) const override
	SingleSampleExperiment (const SingleSampleExperiment &)=delete
SingleSampleExperiment &	operator= (const SingleSampleExperiment &)=delete
	SingleSampleExperiment (SingleSampleExperiment &&s)
void	SetMaskHist (double xmin=0, double xmax=-1, double ymin=0, double ymax=-1)
Public Member Functions inherited from ana::IExperiment
virtual	~IExperiment ()

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

Protected Attributes
const IPrediction *	fMC
Spectrum	fData
Spectrum	fCosmicInTime
Spectrum	fCosmicOutOfTime
TH1 *	fMask

Detailed Description

Compare a single data spectrum to the MC + cosmics expectation.

Definition at line 10 of file SingleSampleExperiment.h.

Constructor & Destructor Documentation

ana::SingleSampleExperiment::SingleSampleExperiment	(	const IPrediction *	pred,
		const Spectrum &	data,
		const Spectrum &	cosmicInTime,
		const Spectrum &	cosmicOutOfTime
	)

Parameters

pred	Source of oscillated MC beam predictions
data	Data spectrum to compare to
cosmicInTime	In-time cosmic ray background component
cosmicOutOfTime	Out-of-time cosmic ray background component

Definition at line 16 of file SingleSampleExperiment.cxx.

    : fMC(pred), fData(data), fCosmicInTime(cosmicInTime), fCosmicOutOfTime(cosmicOutOfTime), fMask(0)
  {
    if(cosmicInTime.POT() > 0){
      std::cout << "SingleSampleExperiment: in-time cosmics have nonzero POT. "
                << "Did you confuse the two cosmics arguments?"
                << std::endl;
      abort();
    }
  }

ana::SingleSampleExperiment::SingleSampleExperiment ( const IPrediction *  pred, const Spectrum &  data  )

inline

In MC studies you might not want to bother with cosmics.

Definition at line 23 of file SingleSampleExperiment.h.

      : fMC(pred), fData(data), fCosmicInTime(0, {}, {}, 0, 0), fCosmicOutOfTime(0, {}, {}, 0, 0), fMask(0)

ana::SingleSampleExperiment::~SingleSampleExperiment ( )

virtual

Definition at line 31 of file SingleSampleExperiment.cxx.

  {
    delete fMask;
  }

ana::SingleSampleExperiment::SingleSampleExperiment ( const SingleSampleExperiment &  )

delete

ana::SingleSampleExperiment::SingleSampleExperiment ( SingleSampleExperiment &&  s )

inline

Definition at line 42 of file SingleSampleExperiment.h.

      : fMC(s.fMC),
        fData(std::move(s.fData)),
        fCosmicInTime(std::move(s.fCosmicInTime)),
        fCosmicOutOfTime(std::move(s.fCosmicOutOfTime))
    {
      s.fMC = nullptr;
    };

Member Function Documentation

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

overridevirtual

Implements ana::IExperiment.

Definition at line 37 of file SingleSampleExperiment.cxx.

  {
    Spectrum p = fMC->PredictSyst(calc, syst);

    if(fCosmicOutOfTime.POT() > 0){ // if out-of-time cosmics supplied
      p += fCosmicOutOfTime;
    }

    TH1D* hpred = p.ToTH1(fData.POT());

    // "Livetime" here means number of readout windows.
    if(fCosmicInTime.Livetime() > 0){ // if in-cosmics supplied
      // This many are already included in the beam MC
      const double beamLivetime = p.Livetime() * fData.POT()/p.POT();
      // So we need to take this many from the cosmics-only to match the data
      const double intimeLivetime = fData.Livetime() - beamLivetime;

      TH1D* hcosmic = fCosmicInTime.ToTH1(intimeLivetime, kLivetime);
      hpred->Add(hcosmic);
      HistCache::Delete(hcosmic);
    }

    TH1D* hdata = fData.ToTH1(fData.POT());

    // If a valid mask has been set, zero out the offending bins
    if (fMask){
      assert(hpred->GetNbinsX() == fMask->GetNbinsX());
      assert(hdata->GetNbinsX() == fMask->GetNbinsX());

      for(int i = 0; i < fMask->GetNbinsX()+2; ++i){
        if (fMask->GetBinContent(i+1) == 1) continue;
        hpred->SetBinContent(i+1, 0);
        hdata->SetBinContent(i+1, 0);
      }
    }

    const double ll = LogLikelihood(hpred, hdata);

    HistCache::Delete(hpred);
    HistCache::Delete(hdata);

    return ll;
  }

ana::SingleSampleExperiment::fCosmicOutOfTime	(	0	,
		{}	,
		{}	,
		0	,
		0
	)

ana::SingleSampleExperiment::fMask ( 0  )

inline

Definition at line 25 of file SingleSampleExperiment.h.

      : fMC(pred), fData(data), fCosmicInTime(0, {}, {}, 0, 0), fCosmicOutOfTime(0, {}, {}, 0, 0), fMask(0)
    {
    }

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

static

Definition at line 99 of file SingleSampleExperiment.cxx.

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

    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"));

    // Legacy format. This is the in-time cosmics
    if(dir->GetDirectory("cosmic")){
      const std::unique_ptr<Spectrum> cosmicInTime = Spectrum::LoadFrom(dir->GetDirectory("cosmic"));
      return std::make_unique<SingleSampleExperiment>(mc, *data, *cosmicInTime, Spectrum(0, {}, {}, 0, 0));
    }

    const std::unique_ptr<Spectrum> cosmicInTime = Spectrum::LoadFrom(dir->GetDirectory("cosmicInTime"));
    const std::unique_ptr<Spectrum> cosmicOutOfTime = Spectrum::LoadFrom(dir->GetDirectory("cosmicOutOfTime"));

    return std::make_unique<SingleSampleExperiment>(mc, *data, *cosmicInTime, *cosmicOutOfTime);
  }

SingleSampleExperiment& ana::SingleSampleExperiment::operator= ( const SingleSampleExperiment &  )

delete

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

overridevirtual

Reimplemented from ana::IExperiment.

Definition at line 83 of file SingleSampleExperiment.cxx.

  {
    TDirectory* tmp = dir;

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

    fMC->SaveTo(dir->mkdir("mc"));
    fData.SaveTo(dir->mkdir("data"));
    fCosmicInTime.SaveTo(dir->mkdir("cosmicInTime"));
    fCosmicOutOfTime.SaveTo(dir->mkdir("cosmicOutOfTime"));

    tmp->cd();
  }

void ana::SingleSampleExperiment::SetMaskHist	(	double	xmin = 0,
		double	xmax = -1,
		double	ymin = 0,
		double	ymax = -1
	)

Definition at line 125 of file SingleSampleExperiment.cxx.

  {
    fMask = GetMaskHist(fData, xmin, xmax, ymin, ymax);
  }

Member Data Documentation

Spectrum ana::SingleSampleExperiment::fCosmicInTime

protected

Definition at line 57 of file SingleSampleExperiment.h.

Spectrum ana::SingleSampleExperiment::fCosmicOutOfTime

protected

Definition at line 57 of file SingleSampleExperiment.h.

Spectrum ana::SingleSampleExperiment::fData

protected

Definition at line 56 of file SingleSampleExperiment.h.

TH1* ana::SingleSampleExperiment::fMask

protected

Definition at line 58 of file SingleSampleExperiment.h.

const IPrediction* ana::SingleSampleExperiment::fMC

protected

Definition at line 55 of file SingleSampleExperiment.h.

---

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

• SingleSampleExperiment.h
• SingleSampleExperiment.cxx