ana::MultiExperimentSBN Class Reference

Combine multiple component experiments. More...

#include <MultiExperimentSBN.h>

Inheritance diagram for ana::MultiExperimentSBN:

Public Member Functions
	MultiExperimentSBN (std::vector< const IExperiment * > expts={}, std::vector< int > exptnames={})
void	Add (const IExperiment *expt, float l)
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 < MultiExperimentSBN >	LoadFrom (TDirectory *dir)

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

Detailed Description

Combine multiple component experiments.

Definition at line 13 of file MultiExperimentSBN.h.

Constructor & Destructor Documentation

ana::MultiExperimentSBN::MultiExperimentSBN ( std::vector< const IExperiment * >  expts = {}, std::vector< int >  exptnames = {}  )

inline

Definition at line 19 of file MultiExperimentSBN.h.

                                                             {},
                       std::vector<int> exptnames = {}) 
        : fExpts(expts), fExptNames(exptnames)
    {
      fSystCorrelations.resize(expts.size());
    }

Member Function Documentation

void ana::MultiExperimentSBN::Add ( const IExperiment *  expt, float  l  )

inline

Definition at line 26 of file MultiExperimentSBN.h.

                                              {
                  fExpts.push_back(expt);
                  fExptNames.push_back(l);
                  }

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

overridevirtual

Implements ana::IExperiment.

Definition at line 13 of file MultiExperimentSBN.cxx.

  {
    assert(fExpts.size() == fExptNames.size());
    
    double ret = 0;

    for(unsigned int n = 0; n < fExpts.size(); ++n){
      // Set baseline for each of the experiments
      osc->SetL(kBLs[fExptNames[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< MultiExperimentSBN > ana::MultiExperimentSBN::LoadFrom ( TDirectory *  dir )

static

Definition at line 103 of file MultiExperimentSBN.cxx.

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

    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<MultiExperimentSBN>(new MultiExperimentSBN(expts));
  }

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

overridevirtual

Reimplemented from ana::IExperiment.

Definition at line 81 of file MultiExperimentSBN.cxx.

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

    if(hasCorr){
      std::cerr << "Warning in MultiExperimentSBN: 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("MultiExperimentSBN").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::MultiExperimentSBN::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 53 of file MultiExperimentSBN.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 << "MultiExperimentSBN::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<int> ana::MultiExperimentSBN::fExptNames

protected

Definition at line 53 of file MultiExperimentSBN.h.

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

protected

Definition at line 52 of file MultiExperimentSBN.h.

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

protected

Definition at line 54 of file MultiExperimentSBN.h.

---

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

• MultiExperimentSBN.h
• MultiExperimentSBN.cxx