ana::OscCurve Class Reference

Transition probability for any one channel as a function of energy. More...

#include <OscCurve.h>

Public Member Functions
	OscCurve (osc::IOscCalc *calc, int from, int to, bool LoverE)
	OscCurve (TH1 *h)
virtual	~OscCurve ()
	OscCurve (const OscCurve &rhs)
OscCurve &	operator= (const OscCurve &rhs)
TH1D *	ToTH1 (bool title=false) const

Protected Attributes
int	fFrom
int	fTo
TH1D *	fHist

Detailed Description

Transition probability for any one channel as a function of energy.

Definition at line 18 of file OscCurve.h.

Constructor & Destructor Documentation

ana::OscCurve::OscCurve	(	osc::IOscCalc *	calc,
		int	from,
		int	to,
		bool	LoverE
	)

Definition at line 28 of file OscCurve.cxx.

    : fFrom(from), fTo(to)
  {
    DontAddDirectory guard;

    fHist = LoverE ? HistCache::New("True L / E (km / GeV);Probability", kTrueLOverEBins) : HistCache::New(";True Energy (GeV);Probability", kTrueEnergyBins);

    // We have extra knowledge that calculators of this type have special
    // modes allowing calculation in L/E and an intrinsic energy smearing.
    OscCalcSterileApprox* approx = DowncastToSterileApprox(calc, true);

    if(approx){
      for(int i = 0; i < fHist->GetNbinsX()+2; ++i){
        if(LoverE){
          const double LElo = fHist->GetXaxis()->GetBinLowEdge(i);
          const double LEhi = fHist->GetXaxis()->GetBinUpEdge(i);
          fHist->SetBinContent(i, approx->P_LoverE(from, to, LElo, LEhi));
        }
        else{
          const double E = fHist->GetBinCenter(i);
          const double Elo = fHist->GetXaxis()->GetBinLowEdge(i);
          const double Ehi = fHist->GetXaxis()->GetBinUpEdge(i);
          // Use 2% resolution (intended to be << the resolution of any actual
          // event) or the bin width, whichever is larger
          fHist->SetBinContent(i, approx->P_range(from, to,
                                                  std::min(Elo, 0.98*E),
                                                  std::max(Ehi, 1.02*E)));
        }
        fHist->SetBinError(i, 0);
      }
    }
    else{
      if(LoverE && !IsNoOscillations(calc)){
        std::cout << "Trying to use a calculator which is not OscCalcSterileApprox with an L/E axis. Will have to code up additional hacks for this to work" << std::endl;
        abort();
      }

      for(int i = 0; i < fHist->GetNbinsX()+2; ++i){
        const double E = fHist->GetBinCenter(i);
        fHist->SetBinContent(i, E > 0 ? calc->P(from, to, E) : 0);
        fHist->SetBinError(i, 0);
      }
    }
  }

ana::OscCurve::OscCurve

(

TH1 *

h

)

Definition at line 74 of file OscCurve.cxx.

  {
    DontAddDirectory guard;

    const TString className = h->ClassName();

    if(className == "TH1D"){
      // Shortcut if types match
      fHist = HistCache::Copy((TH1D*)h);
    }
    else{
      fHist = HistCache::New("", h->GetXaxis());
      fHist->Add(h);
    }

    fHist->SetTitle(";True Energy (GeV);Probability");
  }

ana::OscCurve::~OscCurve ( )

virtual

Definition at line 93 of file OscCurve.cxx.

  {
    if(fHist && fHist->GetDirectory()){
      static bool once = true;
      if(once){
        once = false;
        std::cerr << "OscCurve's fHist is associated with a directory. How did that happen?" << std::endl;
      }
    }

    HistCache::Delete(fHist);
  }

ana::OscCurve::OscCurve

(

const OscCurve &

rhs

)

Definition at line 107 of file OscCurve.cxx.

  {
    DontAddDirectory guard;

    assert(rhs.fHist);
    fHist = HistCache::Copy(rhs.fHist);
  }

Member Function Documentation

OscCurve & ana::OscCurve::operator=

(

const OscCurve &

rhs

)

Definition at line 116 of file OscCurve.cxx.

  {
    if(&rhs == this) return *this;

    DontAddDirectory guard;

    HistCache::Delete(fHist);
    assert(rhs.fHist);
    fHist = HistCache::Copy(rhs.fHist);
    return *this;
  }

TH1D * ana::OscCurve::ToTH1

(

bool

title = false

)

const

Definition at line 129 of file OscCurve.cxx.

  {
    // Could have a file temporarily open
    DontAddDirectory guard;

    TH1D* ret = HistCache::Copy(fHist);

    if(title){
      // Don't do this work unless it's explicitly requested
      std::map<int, std::string> nus;
      nus[12] = nus[-12] = "e";
      nus[14] = nus[-14] = "#mu";
      nus[16] = nus[-16] = "#tau";
      nus[0] = "active";
      const std::string nu = (fFrom > 0) ? "#nu" : "#bar#nu";

      ret->SetTitle((nu+"_{"+nus[fFrom]+"}#rightarrow"+nu+"_{"+nus[fTo]+"}").c_str());
    }

    return ret;
  }

Member Data Documentation

int ana::OscCurve::fFrom

protected

Definition at line 30 of file OscCurve.h.

TH1D* ana::OscCurve::fHist

protected

Definition at line 31 of file OscCurve.h.

int ana::OscCurve::fTo

protected

Definition at line 30 of file OscCurve.h.

---

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

• OscCurve.h
• OscCurve.cxx