ana::MultiExperiment Class Reference

Combine multiple component experiments. More...

#include <MultiExperiment.h>

Inheritance diagram for ana::MultiExperiment:

Public Member Functions
	MultiExperiment (std::vector< const IExperiment * > expts={})
void	Add (const IExperiment *expt)
virtual double	ChiSq (osc::IOscCalcAdjustable *osc, const SystShifts &syst=SystShifts::Nominal()) const override
void	SetSystCorrelations (int idx, const std::vector< std::pair< const ISyst *, const ISyst * >> &corrs)
virtual void	SaveTo (TDirectory *dir) const override
Public Member Functions inherited from ana::IExperiment
virtual	~IExperiment ()

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

Protected Attributes
std::vector< const IExperiment * >	fExpts
std::vector< std::vector < std::pair< const ISyst *, const ISyst * > > >	fSystCorrelations

Detailed Description

Combine multiple component experiments.

Definition at line 11 of file MultiExperiment.h.

Constructor & Destructor Documentation

ana::MultiExperiment::MultiExperiment ( std::vector< const IExperiment * >  expts = {} )

inline

Definition at line 14 of file MultiExperiment.h.

                                                          {}) : fExpts(expts)
    {
      fSystCorrelations.resize(expts.size());
    }

Member Function Documentation

void ana::MultiExperiment::Add ( const IExperiment *  expt )

inline

Definition at line 19 of file MultiExperiment.h.

{fExpts.push_back(expt);}

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

overridevirtual

Implements ana::IExperiment.

Definition at line 13 of file MultiExperiment.cxx.

  {
    double ret = 0;
    for(unsigned int n = 0; n < fExpts.size(); ++n){
      // Don't waste time fiddling with the systematics if for sure there's
      // nothing to do.
      if(fSystCorrelations[n].empty()){
        ret += fExpts[n]->ChiSq(osc, syst);
      }
      else{
        // Make a local copy we're going to rewrite into the terms this
        // sub-experiment will accept.
        SystShifts localShifts = syst;
        for(auto it: fSystCorrelations[n]){
          // We're mapping prim -> sec
          const ISyst* prim = it.first;
          const ISyst* sec = it.second;
          if(syst.GetShift(prim) != 0){
            // sec can be unset, which means there's no representation needed
            // of prim in the sub-experiment.
            if(sec) localShifts.SetShift(sec, syst.GetShift(prim));
            // We've either translated or discarded prim, so drop it here.
            localShifts.SetShift(prim, 0);
          }
        }
        ret += fExpts[n]->ChiSq(osc, localShifts);
      }
    }
    return ret;
  }

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

static

Definition at line 97 of file MultiExperiment.cxx.

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

    std::vector<const IExperiment*> expts;

    for(int i = 0; ; ++i){
      TDirectory* subdir = dir->GetDirectory(TString::Format("expt%d", i));
      if(!subdir) break;

      expts.push_back(ana::LoadFrom<IExperiment>(subdir).release());
    }

    assert(!expts.empty());

    return std::unique_ptr<MultiExperiment>(new MultiExperiment(expts));
  }

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

overridevirtual

Reimplemented from ana::IExperiment.

Definition at line 75 of file MultiExperiment.cxx.

  {
    bool hasCorr = false;
    for(auto it: fSystCorrelations) if(!it.empty()) hasCorr = true;

    if(hasCorr){
      std::cerr << "Warning in MultiExperiment: systematic correlations are set and will not be serialized by this call to SaveTo(). You will have to re-set them once you load the experiment back in." << std::endl;
    }

    TDirectory* tmp = dir;

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

    for(unsigned int i = 0; i < fExpts.size(); ++i){
      fExpts[i]->SaveTo(dir->mkdir(TString::Format("expt%d", i)));
    }

    tmp->cd();
  }

void ana::MultiExperiment::SetSystCorrelations	(	int	idx,
		const std::vector< std::pair< const ISyst *, const ISyst * >> &	corrs
	)

For the subexperiment idx, set up a mapping between systematics

Each element in the vector is a pair from a "primary" systematic to a "secondary". When this MultiExperiment is called with a primary systematic shifted, the sub-experiment will be called with the secondary systematic set to the same value (and the primary unset).

You can pass NULL for a secondary to indicate that the systematic simply has no effect on the experiment in question and should be filtered out.

Definition at line 47 of file MultiExperiment.cxx.

  {
    // Sanity-check the mapping
    std::map<const ISyst*, const ISyst*> already;
    for(auto it: corrs){
      assert(it.first != it.second);

      // Don't worry if second element is null pointer
      if (!it.second) continue;
      if(already.find(it.second) == already.end()){
        already[it.second] = it.first;
      }
      else{
        std::cout << "MultiExperiment::SetSystCorrelations(): Warning!\n"
                  << "In experiment " << idx << " both "
                  << already[it.second]->ShortName() << " and "
                  << it.first->ShortName()
                  << " are configured to map to " << it.second->ShortName()
                  << ". That's probably not what you want." << std::endl;
      }
    }

    // Apply it
    fSystCorrelations[idx] = corrs;
  }

Member Data Documentation

std::vector<const IExperiment*> ana::MultiExperiment::fExpts

protected

Definition at line 42 of file MultiExperiment.h.

std::vector<std::vector<std::pair<const ISyst*, const ISyst*> > > ana::MultiExperiment::fSystCorrelations

protected

Definition at line 44 of file MultiExperiment.h.

---

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

• MultiExperiment.h
• MultiExperiment.cxx