ana::TwoSampleDecomp Class Reference

Class for data driven decomposition of signal and background using two different samples with slightly different spectra. More...

#include "/cvmfs/nova-development.opensciencegrid.org/novasoft/releases/N21-02-25/CAFAna/Decomp/TwoSampleDecomp.h"

Inheritance diagram for ana::TwoSampleDecomp:

Public Member Functions
	TwoSampleDecomp (SpectrumLoaderBase &loaderMC, SpectrumLoaderBase &loaderData, std::string label, const Binning &bins, const Var &var, const Cut &firstSample, const Cut &isSignal=kNoCut, const Cut &anaCut=kNoCut, const SystShifts &shift=kNoShift, const Var &weight=kUnweighted)
	Constructor for creating spectra from a SpectrumLoader. Must supply cuts to define samples, first or second. For instance first: left side, second: right side. Also must define cut to specefiy signal/background. Analysis cuts are optional, but likely necessary. More...
	TwoSampleDecomp (SpectrumLoaderBase &loaderMC, SpectrumLoaderBase &loaderData, const HistAxis &axis, const Cut &firstSample, const Cut &isSignal=kNoCut, const Cut &anaCut=kNoCut, const SystShifts &shift=kNoShift, const Var &weight=kUnweighted)
virtual	~TwoSampleDecomp ()
virtual Spectrum	NCTotalComponent () const override
	Returns the decomposed NC component. More...
virtual Spectrum	NCComponent () const override
virtual Spectrum	NCAntiComponent () const override
virtual Spectrum	NumuComponent () const override
	Returns the decomposed numu CC component. More...
virtual Spectrum	NueComponent () const override
	Could return the nue CC component, but not implemented for two sample decomposition. More...
virtual Spectrum	AntiNumuComponent () const override
virtual Spectrum	AntiNueComponent () const override
void	Decomp () const
	Preform the decomposition using the ratios in the MC sample. More...
void	RatioCalc () const
	Calculate the ratios for the signal and background, second/first. More...
void	SaveTo (TDirectory *dir, const std::string &name) const override
void	DrawSigBkgOverlay ()
void	DrawTwoSamplesWithRatios ()
void	DrawRatios ()
virtual Spectrum	Data_Component () const
virtual Spectrum	MC_NCTotalComponent () const
virtual Spectrum	MC_NCComponent () const
virtual Spectrum	MC_NCAntiComponent () const
virtual Spectrum	MC_NumuComponent () const
virtual Spectrum	MC_AntiNumuComponent () const
virtual Spectrum	MC_NueComponent () const
virtual Spectrum	MC_AntiNueComponent () const

Static Public Member Functions
static std::unique_ptr< TwoSampleDecomp >	LoadFrom (TDirectory *dir, const std::string &name)

Public Attributes
const SpectrumLoaderBase &	fLoaderMC
const SpectrumLoaderBase &	fLoaderData
const HistAxis *	fAxis
	Axis used for spectra. More...
Spectrum	fFirstSig
	Spectrum for first sample, signal. More...
Spectrum	fSecondSig
	Spectrum for second sample, signal. More...
Spectrum	fFirstBkg
	Spectrum for first sample, background. More...
Spectrum	fSecondBkg
	Spectrum for second sample, background. More...
Spectrum	fTrueSig
	True signal spectrum. More...
Spectrum	fTrueBkg
	True backgroundh spectrum. More...
Spectrum	fFirst
	Spectrum of events that pass the firstSample cut. More...
Spectrum	fSecond
	Spectrum of events that fail the firstSample cut. More...
Spectrum	fSig
	Decomposed spectrum of signal events. More...
Spectrum	fBkg
	Decomposed spectrum of signal events. More...
TH1D *	fSigRatio
	Ratio of second/first spectra for signal. More...
TH1D *	fBkgRatio
	Ratio of second/first spectra for bkg. More...
bool	fDecomposed
	Indicates whether decomposition was done. More...

Protected Member Functions
	TwoSampleDecomp ()

Detailed Description

Class for data driven decomposition of signal and background using two different samples with slightly different spectra.

Definition at line 21 of file TwoSampleDecomp.h.

Constructor & Destructor Documentation

ana::TwoSampleDecomp::TwoSampleDecomp	(	SpectrumLoaderBase &	loaderMC,
		SpectrumLoaderBase &	loaderData,
		std::string	label,
		const Binning &	bins,
		const Var &	var,
		const Cut &	firstSample,
		const Cut &	isSignal = kNoCut,
		const Cut &	anaCut = kNoCut,
		const SystShifts &	shift = kNoShift,
		const Var &	weight = kUnweighted
	)

Constructor for creating spectra from a SpectrumLoader. Must supply cuts to define samples, first or second. For instance first: left side, second: right side. Also must define cut to specefiy signal/background. Analysis cuts are optional, but likely necessary.

Detailed description starts here.

Parameters

loaderMC	Spectrum loader for MC
loaderData	Spectrum loader for data
label	Name of the thing
bins	Desired binning for spectra
var	Variable to be plotted in spectra
firstSample	Cut to define samples, true: first, false: second
isSignal	Signal or not? true: sig, false: bkg
anaCut	Cut to define selection(e.g. numu CC)
weight	Weight for events in spectrum

Definition at line 24 of file TwoSampleDecomp.cxx.

     :TwoSampleDecomp(loaderMC, loaderData, HistAxis(label, bins, var), 
                      firstSample, isSignal, anaCut, shift, weight)
    {}

ana::TwoSampleDecomp::TwoSampleDecomp	(	SpectrumLoaderBase &	loaderMC,
		SpectrumLoaderBase &	loaderData,
		const HistAxis &	axis,
		const Cut &	firstSample,
		const Cut &	isSignal = kNoCut,
		const Cut &	anaCut = kNoCut,
		const SystShifts &	shift = kNoShift,
		const Var &	weight = kUnweighted
	)

Parameters

loaderMC	Spectrum loader for MC
loaderData	Spectrum loader for data
axis	Variable to be plotted in spectra
firstSample	Cut to define samples, true: first, false: second
isSignal	Signal or not? true: sig, false: bkg
anaCut	Cut to define selection(e.g. numu CC)
weight	Weight for events in spectrum

Definition at line 42 of file TwoSampleDecomp.cxx.

References om::cout, allTimeWatchdog::endl, and ana::LabelsAndBins::NDimensions().

     :
  fLoaderMC  (loaderMC),
  fLoaderData(loaderData),
  fAxis(&axis),
  fFirstSig (loaderMC, axis,  firstSample &&  isSignal && anaCut, 
             shift, weight),
  fSecondSig(loaderMC, axis, !firstSample &&  isSignal && anaCut, 
             shift, weight),
  fFirstBkg (loaderMC, axis,  firstSample && !isSignal && anaCut, 
             shift, weight),
  fSecondBkg(loaderMC, axis, !firstSample && !isSignal && anaCut, 
             shift, weight),
  fTrueSig  (loaderMC, axis,                  isSignal && anaCut, 
             shift, weight),
  fTrueBkg  (loaderMC, axis,                 !isSignal && anaCut, 
             shift, weight),
  fFirst    (loaderData, axis,  firstSample            && anaCut, 
             shift, weight),
  fSecond   (loaderData, axis, !firstSample            && anaCut, 
             shift, weight),
  fSig      (Spectrum::Uninitialized()),
  fBkg      (Spectrum::Uninitialized()),
  fSigRatio (0),
  fBkgRatio (0),
  fDecomposed(false)
  {
    if(axis.NDimensions() != 1){
      std::cout << "Error: TwoSampleDecomp not equipped to handle multi-dimensional axis." << std::endl;
      abort();
    }
  }

virtual ana::TwoSampleDecomp::~TwoSampleDecomp ( )

inlinevirtual

Definition at line 57 of file TwoSampleDecomp.h.

{};

ana::TwoSampleDecomp::TwoSampleDecomp ( )

protected

Definition at line 86 of file TwoSampleDecomp.cxx.

                                  :
  fLoaderMC  (kNullLoader),
  fLoaderData(kNullLoader),
  fFirstSig (Spectrum::Uninitialized()),
  fSecondSig(Spectrum::Uninitialized()),
  fFirstBkg (Spectrum::Uninitialized()),
  fSecondBkg(Spectrum::Uninitialized()),
  fTrueSig  (Spectrum::Uninitialized()),
  fTrueBkg  (Spectrum::Uninitialized()),
  fFirst    (Spectrum::Uninitialized()),
  fSecond   (Spectrum::Uninitialized()),
  fSig      (Spectrum::Uninitialized()),
  fBkg      (Spectrum::Uninitialized()),
  fSigRatio (0),
  fBkgRatio (0),
  fDecomposed(false)
  {  }

Member Function Documentation

virtual Spectrum ana::TwoSampleDecomp::AntiNueComponent ( ) const

inlineoverridevirtual

Implements ana::IDecomp.

Definition at line 75 of file TwoSampleDecomp.h.

References Decomp(), dir, DrawRatios(), DrawSigBkgOverlay(), DrawTwoSamplesWithRatios(), fBkg, LoadFrom(), RatioCalc(), SaveTo(), and string.

{Decomp(); return fBkg;}

virtual Spectrum ana::TwoSampleDecomp::AntiNumuComponent ( ) const

inlineoverridevirtual

Implements ana::IDecomp.

Definition at line 74 of file TwoSampleDecomp.h.

References Decomp(), and fBkg.

{Decomp(); return fBkg;}

Spectrum ana::IDecomp::Data_Component ( ) const

virtualinherited

Reimplemented in ana::BENDecomp, ana::MichelDecomp, ana::FluxDecomp, ana::ProportionalDecomp, ana::FakeDecomp, ana::NumuDecomp, ana::NueDecomp, and ana::NCDecomp.

Definition at line 16 of file IDecomp.cxx.

References om::cout, and allTimeWatchdog::endl.

Referenced by ana::GetNDComponents(), ana::GetNDCompsFromDecomp(), ana::GetNDDecompsFromDecomp(), and ana::IDecomp::NCTotalComponent().

  { 
    std::cout<<"IDecomp::Data_Componet() is not implemented for this class, check your Decomp"<<std::endl;
    abort();
  }

void ana::TwoSampleDecomp::Decomp

(

)

const

Preform the decomposition using the ratios in the MC sample.

Definition at line 160 of file TwoSampleDecomp.cxx.

References ana::assert(), fBkg, fBkgRatio, fDecomposed, fFirst, fLoaderData, fLoaderMC, fSecond, fSig, fSigRatio, ana::Spectrum::GetBinnings(), ana::Spectrum::GetLabels(), ana::SpectrumLoaderBase::Gone(), MECModelEnuComparisons::i, ana::Spectrum::Livetime(), plotROC::nBins, ana::Spectrum::POT(), RatioCalc(), util::sqr(), std::sqrt(), string, plotROC::title, and ana::Spectrum::ToTH1().

Referenced by AntiNueComponent(), AntiNumuComponent(), DrawSigBkgOverlay(), DrawTwoSamplesWithRatios(), NCTotalComponent(), NueComponent(), and NumuComponent().

                                    {
    if(fDecomposed) return;

    assert(fLoaderMC.Gone() 
          && "Call Go() on MC SpectrumLoader before TwoSampleDecomp::Decomp");
    assert(fLoaderData.Gone() 
          && "Call Go() on Data SpectrumLoader before TwoSampleDecomp::Decomp");

    RatioCalc();

    TH1D* hFirst = fFirst.ToTH1(fFirst.POT());
    TH1D* hSecond = fSecond.ToTH1(fSecond.POT());
    
    // Get the number of bins from the axis 
    int nBins = hFirst->GetNbinsX();

    // Make sure the histogram binning is sane 
    assert(hFirst->GetNbinsX() == nBins && 
       "Inconsistent binning across spectra.");
    assert(hSecond->GetNbinsX() == nBins && 
       "Inconsistent binning across spectra.");
    assert(fSigRatio->GetNbinsX() == nBins && 
       "MC and data have different numbers of bins.");
    assert(fBkgRatio->GetNbinsX() == nBins && 
       "MC and data have different numbers of bins.");

    std::stringstream title;
    title  << std::string(";") << hFirst->GetXaxis()->GetTitle() 
           << std::string(";Events");


    float axisMin = hFirst->GetBinLowEdge(1);
    float axisMax = hFirst->GetBinLowEdge(nBins+1);

    TH1D* hSig = new TH1D("Sig", title.str().c_str(), nBins, axisMin, axisMax);
    TH1D* hBkg = new TH1D("Bkg", title.str().c_str(), nBins, axisMin, axisMax);

    for(int i = 1; i <= nBins; ++i)
    {
      //These guys are what we're after, signal and background. Initialize them.  
      double nSig    = 0;
      double nBkg    = 0;
      double errSig  = 0;
      double errBkg  = 0;


      // Grab all of the things out of the histograms and store them so that 
      // the math looks a bit nicer.
      double rSig    = fSigRatio->GetBinContent(i);
      double rSigErr = fSigRatio->GetBinError(i);
      double rBkg    = fBkgRatio->GetBinContent(i);
      double rBkgErr = fBkgRatio->GetBinError(i);
      double nFirst  = hFirst->GetBinContent(i);
      double nSecond = hSecond->GetBinContent(i);

      // d2 is the square of the determinant of the matrix we're inverting.  
      // It shows up in the math a lot, so we'll store it for neatness.  
      double d2  = util::sqr(rSig - rBkg);


      //We're about to divide by d2 until the cows come home, so make sure it's
      //not zero.  When it is zero, it's probably because those bins were empty.  
      if  (d2 != 0) 
      {

        //Do the math in a way that stores the common terms in the error formula
        double sSig  = (1+rSig); 
        double sBkg  = (1+rBkg);
        nSig      = sSig * (nSecond - rBkg*nFirst) / (rSig - rBkg);
        nBkg      = sBkg * (nFirst*rSig - nSecond) / (rSig - rBkg);

        // Here's the error formula 
        errSig = sqrt((nSecond * util::sqr(sSig))/d2 + 
                    (nFirst*util::sqr(sSig) *util::sqr(rBkg))/d2 + 
                    util::sqr(rSigErr) * (-((sSig*(nSecond - nFirst*rBkg))/d2) + 
                    util::sqr(nSecond - nFirst*rBkg)/(rSig - rBkg)) + 
                    util::sqr(-((nFirst*(sSig))/(rSig - rBkg)) + 
                    (sSig*(nSecond - nFirst*rBkg))/d2) * util::sqr(rBkgErr));
        errBkg = sqrt((nSecond * util::sqr(sBkg))/d2 + 
                (nFirst*util::sqr(rSig) * util::sqr(sBkg))/d2 + 
                util::sqr(rSigErr)*util::sqr(-(((-nSecond + nFirst*rSig)*sBkg)/
                (rSig - rBkg)*2) + 
                (nFirst*sBkg)/(rSig - rBkg)) + 
                util::sqr((-nSecond + nFirst*rSig)/(rSig - rBkg) + 
                ((-nSecond + nFirst*rSig) *sBkg)/d2)*util::sqr(rBkgErr));
      }
      
      // Set the values in the resultant histograms, signal and background.  
      hSig->SetBinContent(i, nSig);
      hBkg->SetBinContent(i, nBkg);
      hSig->SetBinError(i, errSig);
      hBkg->SetBinError(i, errBkg);
    }

    fSig = Spectrum(Eigen::ArrayXd(Eigen::Map<Eigen::ArrayXd>(hSig->GetArray(), hSig->GetNbinsX()+2)),
                    HistAxis(fFirst.GetLabels(), fFirst.GetBinnings()),
                    fFirst.POT(), fFirst.Livetime());
    fBkg = Spectrum(Eigen::ArrayXd(Eigen::Map<Eigen::ArrayXd>(hBkg->GetArray(), hBkg->GetNbinsX()+2)),
                    HistAxis(fFirst.GetLabels(), fFirst.GetBinnings()),
                    fFirst.POT(), fFirst.Livetime());

    fDecomposed = true;
  }

void ana::TwoSampleDecomp::DrawRatios

(

)

Definition at line 408 of file TwoSampleDecomp.cxx.

References fBkgRatio, and fSigRatio.

Referenced by AntiNueComponent(), and runTwoSampleDecomp().

  {
    fSigRatio->Draw("");
    fBkgRatio->SetLineStyle(2);
    fBkgRatio->Draw("same");
  }

void ana::TwoSampleDecomp::DrawSigBkgOverlay

(

)

Definition at line 333 of file TwoSampleDecomp.cxx.

References Decomp(), fBkg, fFirst, fSecond, fSig, fTrueBkg, fTrueSig, kBlue, kGreen, kRed, ana::Spectrum::POT(), and ana::Spectrum::ToTH1().

Referenced by AntiNueComponent(), and runTwoSampleDecomp().

  {
    Decomp();

    Spectrum data = fFirst + fSecond;

    TH1D* hData = data.ToTH1(data.POT());
    hData->Draw();
  
    TH1D* hSig = fSig.ToTH1(data.POT());
    TH1D* hBkg = fBkg.ToTH1(data.POT());
  
    hSig->SetLineColor(kGreen+2);
    hBkg->SetLineColor(kBlue-4);
    hSig->SetMarkerColor(kGreen+2);
    hBkg->SetMarkerColor(kBlue-4);
  
    hSig->Draw("same");
    hBkg->Draw("same");
  
  
    TH1D* hTrueSig = fTrueSig.ToTH1(data.POT());
    TH1D* hTrueBkg = fTrueBkg.ToTH1(data.POT());
  
    hTrueSig->SetLineColor(kGreen+2);
    hTrueBkg->SetLineColor(kBlue-4);
  
    hTrueSig->Draw("hist same");
    hTrueBkg->Draw("hist same");
  
    TH1D* hTrueTotal = (fTrueSig + fTrueBkg).ToTH1(data.POT());
  
    hTrueTotal->SetLineColor(kRed);
    hTrueTotal->Draw("hist same");
    
      
  }

void ana::TwoSampleDecomp::DrawTwoSamplesWithRatios

(

)

Definition at line 372 of file TwoSampleDecomp.cxx.

References Decomp(), fFirst, fFirstBkg, fFirstSig, fSecondBkg, fSecondSig, make_syst_table_plots::h, hists, kBlue, kRed, m, cet::sqlite::max(), ana::Spectrum::POT(), and ana::Spectrum::ToTH1().

Referenced by AntiNueComponent(), and runTwoSampleDecomp().

  {
    Decomp();

    TH1D* hFirstSig = fFirstSig.ToTH1(fFirst.POT());
    TH1D* hSecondSig = fSecondSig.ToTH1(fFirst.POT());
    TH1D* hFirstBkg = fFirstBkg.ToTH1(fFirst.POT());
    TH1D* hSecondBkg = fSecondBkg.ToTH1(fFirst.POT());

    
    std::vector<TH1D*> hists = {hFirstSig, hSecondSig, hFirstBkg, hSecondBkg};

    float max = 0;
    for(const auto& h:hists){
      float m = h->GetBinContent(h->GetMaximumBin());
      if(m > max) max = m;
    }

    for(const auto& h:hists) h->SetMaximum(1.2*max);

    hFirstSig->SetLineColor(kRed);
    hSecondSig->SetLineColor(kBlue-4);
    hFirstBkg->SetLineColor(kRed);
    hSecondBkg->SetLineColor(kBlue-4);

    hFirstBkg->SetLineStyle(2);
    hSecondBkg->SetLineStyle(2);

    hFirstSig->Draw("hist");
    hSecondSig->Draw("hist same");
    hFirstBkg->Draw("hist same");
    hSecondBkg->Draw("hist same");
        
    
  }

std::unique_ptr< TwoSampleDecomp > ana::TwoSampleDecomp::LoadFrom ( TDirectory *  dir, const std::string &  name  )

static

Definition at line 133 of file TwoSampleDecomp.cxx.

References ana::assert(), dir, ana::Spectrum::LoadFrom(), and runNovaSAM::ret.

Referenced by AntiNueComponent().

  {
    dir = dir->GetDirectory(name.c_str()); // switch to subdir
    assert(dir);

    std::unique_ptr<TwoSampleDecomp> ret(new TwoSampleDecomp);

    ret->fFirstSig = *Spectrum::LoadFrom(dir, "FirstSig");
    ret->fSecondSig = *Spectrum::LoadFrom(dir, "SecondSig");
    ret->fFirstBkg = *Spectrum::LoadFrom(dir, "FirstBkg");
    ret->fSecondBkg = *Spectrum::LoadFrom(dir, "SecondBkg");
    ret->fTrueSig = *Spectrum::LoadFrom(dir, "TrueSig");
    ret->fTrueBkg = *Spectrum::LoadFrom(dir, "TrueBkg");
    ret->fFirst = *Spectrum::LoadFrom(dir, "First");
    ret->fSecond = *Spectrum::LoadFrom(dir, "Second");
    ret->fSig = *Spectrum::LoadFrom(dir, "Sig");
    ret->fBkg = *Spectrum::LoadFrom(dir, "Bkg");    

    delete dir;

    return ret;
  }

Spectrum ana::IDecomp::MC_AntiNueComponent ( ) const

virtualinherited

Reimplemented in ana::BENDecomp, ana::MichelDecomp, ana::FluxDecomp, ana::ProportionalDecomp, and ana::NumuDecomp.

Definition at line 46 of file IDecomp.cxx.

References om::cout, and allTimeWatchdog::endl.

Referenced by ana::GetNDCompsFromDecomp(), and ana::IDecomp::NCTotalComponent().

  {
    std::cout << "IDecomp: not implemented" << std::endl; abort();
  }

Spectrum ana::IDecomp::MC_AntiNumuComponent ( ) const

virtualinherited

Reimplemented in ana::BENDecomp, ana::MichelDecomp, ana::FluxDecomp, ana::ProportionalDecomp, and ana::NumuDecomp.

Definition at line 38 of file IDecomp.cxx.

References om::cout, and allTimeWatchdog::endl.

Referenced by ana::GetNDCompsFromDecomp(), and ana::IDecomp::NCTotalComponent().

  {
    std::cout << "IDecomp: not implemented" << std::endl; abort();
  }

Spectrum ana::IDecomp::MC_NCAntiComponent ( ) const

virtualinherited

Reimplemented in ana::MichelDecomp, ana::FluxDecomp, and ana::ProportionalDecomp.

Definition at line 30 of file IDecomp.cxx.

References om::cout, and allTimeWatchdog::endl.

Referenced by ana::IDecomp::NCTotalComponent().

  {
    std::cout << "IDecomp: not implemented" << std::endl; abort();
  }

Spectrum ana::IDecomp::MC_NCComponent ( ) const

virtualinherited

Reimplemented in ana::MichelDecomp, ana::FluxDecomp, and ana::ProportionalDecomp.

Definition at line 26 of file IDecomp.cxx.

References om::cout, and allTimeWatchdog::endl.

Referenced by ana::IDecomp::NCTotalComponent().

  {
    std::cout << "IDecomp: not implemented" << std::endl; abort();
  }

Spectrum ana::IDecomp::MC_NCTotalComponent ( ) const

virtualinherited

Reimplemented in ana::BENDecomp, ana::MichelDecomp, ana::FluxDecomp, ana::ProportionalDecomp, and ana::NumuDecomp.

Definition at line 22 of file IDecomp.cxx.

References om::cout, and allTimeWatchdog::endl.

Referenced by ana::GetNDCompsFromDecomp(), and ana::IDecomp::NCTotalComponent().

  {
    std::cout << "IDecomp: not implemented" << std::endl; abort();
  }

Spectrum ana::IDecomp::MC_NueComponent ( ) const

virtualinherited

Reimplemented in ana::BENDecomp, ana::MichelDecomp, ana::FluxDecomp, ana::ProportionalDecomp, and ana::NumuDecomp.

Definition at line 42 of file IDecomp.cxx.

References om::cout, and allTimeWatchdog::endl.

Referenced by ana::GetNDCompsFromDecomp(), and ana::IDecomp::NCTotalComponent().

  {
    std::cout << "IDecomp: not implemented" << std::endl; abort();
  }

Spectrum ana::IDecomp::MC_NumuComponent ( ) const

virtualinherited

Reimplemented in ana::BENDecomp, ana::MichelDecomp, ana::FluxDecomp, ana::ProportionalDecomp, and ana::NumuDecomp.

Definition at line 34 of file IDecomp.cxx.

References om::cout, and allTimeWatchdog::endl.

Referenced by ana::GetNDCompsFromDecomp(), and ana::IDecomp::NCTotalComponent().

  {
    std::cout << "IDecomp: not implemented" << std::endl; abort();
  }

Spectrum ana::TwoSampleDecomp::NCAntiComponent ( ) const

overridevirtual

Implements ana::IDecomp.

Definition at line 158 of file TwoSampleDecomp.cxx.

References om::cout, and allTimeWatchdog::endl.

Referenced by NCTotalComponent().

{std::cout << "TwoSampleDecomp::NCAntiComponent() not implemented" << std::endl; abort();}

Spectrum ana::TwoSampleDecomp::NCComponent ( ) const

overridevirtual

Implements ana::IDecomp.

Definition at line 156 of file TwoSampleDecomp.cxx.

References om::cout, and allTimeWatchdog::endl.

Referenced by NCTotalComponent().

{std::cout << "TwoSampleDecomp::NCComponent() not implemented" << std::endl; abort();}

virtual Spectrum ana::TwoSampleDecomp::NCTotalComponent ( ) const

inlineoverridevirtual

Returns the decomposed NC component.

Reimplemented from ana::IDecomp.

Definition at line 62 of file TwoSampleDecomp.h.

References Decomp(), fBkg, NCAntiComponent(), and NCComponent().

{Decomp(); return fBkg;};

virtual Spectrum ana::TwoSampleDecomp::NueComponent ( ) const

inlineoverridevirtual

Could return the nue CC component, but not implemented for two sample decomposition.

Implements ana::IDecomp.

Definition at line 72 of file TwoSampleDecomp.h.

References Decomp(), and fBkg.

{Decomp(); return fBkg;};

virtual Spectrum ana::TwoSampleDecomp::NumuComponent ( ) const

inlineoverridevirtual

Returns the decomposed numu CC component.

Implements ana::IDecomp.

Definition at line 68 of file TwoSampleDecomp.h.

References Decomp(), and fSig.

{Decomp(); return fSig;};

void ana::TwoSampleDecomp::RatioCalc

(

)

const

Calculate the ratios for the signal and background, second/first.

Definition at line 266 of file TwoSampleDecomp.cxx.

References ana::assert(), ana::bins, fAxis, fBkgRatio, fFirstBkg, fFirstSig, fSecondBkg, fSecondSig, fSigRatio, ana::LabelsAndBins::GetBinnings(), MECModelEnuComparisons::i, plotROC::nBins, ana::Spectrum::POT(), string, plotROC::title, and ana::Spectrum::ToTH1().

Referenced by AntiNueComponent(), and Decomp().

  {


    // Check to make sure we haven't already computed ratios 
    // histogram pointers should be null
    assert(!fSigRatio && "Ratio histograms already constructed.");
    assert(!fBkgRatio && "Ratio histograms already constructed.");

    // Pull histograms out of the required Spectrum object
    TH1D* hFirstSig  = fFirstSig .ToTH1(fFirstSig.POT());        
    TH1D* hSecondSig = fSecondSig.ToTH1(fSecondSig.POT());      
    TH1D* hFirstBkg  = fFirstBkg .ToTH1(fFirstBkg.POT());   
    TH1D* hSecondBkg = fSecondBkg.ToTH1(fSecondBkg.POT());

    // Make the axis labels for the ratio histograms out of the 
    std::stringstream title;
    title  << std::string(";") << hFirstSig->GetXaxis()->GetTitle() 
           << std::string(";Ratio");

    // Get the number of bins from the axis 
    int nBins = hFirstSig->GetNbinsX();

    // Make sure the histogram binning is sane 
    assert(hFirstSig->GetNbinsX() == nBins &&
                                       "Inconsistent binning across spectra.");
    assert(hSecondSig->GetNbinsX() == nBins &&
                                       "Inconsistent binning across spectra.");
    assert(hFirstBkg->GetNbinsX() == nBins &&
                                       "Inconsistent binning across spectra.");
    assert(hSecondBkg->GetNbinsX() == nBins &&
                                       "Inconsistent binning across spectra.");

    const Binning bins = fAxis->GetBinnings()[0];

    // Construct the ratio histogram objects 
    fSigRatio = new TH1D("SigRatio", title.str().c_str(), nBins, 
                         bins.Min(), bins.Max());
    fBkgRatio = new TH1D("BkgRatio", title.str().c_str(), nBins, 
                         bins.Min(), bins.Max());


    // Loop over the bins in the histograms to compute the ratios bin-by-bin
    for (int i = 1; i <= nBins; ++i)
    {
      // Ratios for signal 
      if(hFirstSig->GetBinContent(i) != 0 and hSecondSig->GetBinContent(i) != 0)
      {
        // if non zero, make an object that does the math
        RatioError ratioError(hSecondSig->GetBinContent(i),  
                              hFirstSig->GetBinContent(i));
        // Set the bins 
        fSigRatio->SetBinContent(i, ratioError.fRatio);
        fSigRatio->SetBinError  (i, ratioError.fError);
      }
      // Repeat for background 
      if(hFirstBkg->GetBinContent(i) != 0 and hSecondBkg->GetBinContent(i) != 0)
      {
        RatioError ratioError(hSecondBkg->GetBinContent(i),  
                              hFirstBkg->GetBinContent(i));
        fBkgRatio->SetBinContent(i, ratioError.fRatio);
        fBkgRatio->SetBinError  (i, ratioError.fError);
      }
    }
  }

void ana::TwoSampleDecomp::SaveTo ( TDirectory *  dir, const std::string &  name  ) const

overridevirtual

Implements ana::IDecomp.

Definition at line 106 of file TwoSampleDecomp.cxx.

References dir, fBkg, fFirst, fFirstBkg, fFirstSig, fSecond, fSecondBkg, fSecondSig, fSig, fTrueBkg, fTrueSig, ana::Spectrum::SaveTo(), and tmp.

Referenced by AntiNueComponent(), and runTwoSampleDecomp().

  {
    TDirectory* tmp = gDirectory;

    dir = dir->mkdir(name.c_str()); // switch to subdir
    dir->cd();

    TObjString("TwoSampleDecomp").Write("type");
                
    fFirstSig.SaveTo(dir, "FirstSig");
    fSecondSig.SaveTo(dir, "SecondSig");
    fFirstBkg.SaveTo(dir, "FirstBkg");
    fSecondBkg.SaveTo(dir, "SecondBkg");
    fTrueSig.SaveTo(dir, "TrueSig");
    fTrueBkg.SaveTo(dir, "TrueBkg");
    fFirst.SaveTo(dir, "First");
    fSecond.SaveTo(dir, "Second");
    fSig.SaveTo(dir, "Sig");
    fBkg.SaveTo(dir, "Bkg");

    dir->Write();
    delete dir;

    tmp->cd();
  }

Member Data Documentation

const HistAxis* ana::TwoSampleDecomp::fAxis

Axis used for spectra.

Definition at line 95 of file TwoSampleDecomp.h.

Referenced by RatioCalc().

Spectrum ana::TwoSampleDecomp::fBkg

mutable

Decomposed spectrum of signal events.

Definition at line 108 of file TwoSampleDecomp.h.

Referenced by AntiNueComponent(), AntiNumuComponent(), Decomp(), DrawSigBkgOverlay(), NCTotalComponent(), NueComponent(), and SaveTo().

TH1D* ana::TwoSampleDecomp::fBkgRatio

mutable

Ratio of second/first spectra for bkg.

Definition at line 112 of file TwoSampleDecomp.h.

Referenced by Decomp(), DrawRatios(), and RatioCalc().

bool ana::TwoSampleDecomp::fDecomposed

mutable

Indicates whether decomposition was done.

Definition at line 113 of file TwoSampleDecomp.h.

Referenced by Decomp().

Spectrum ana::TwoSampleDecomp::fFirst

Spectrum of events that pass the firstSample cut.

Definition at line 104 of file TwoSampleDecomp.h.

Referenced by Decomp(), DrawSigBkgOverlay(), DrawTwoSamplesWithRatios(), and SaveTo().

Spectrum ana::TwoSampleDecomp::fFirstBkg

Spectrum for first sample, background.

Definition at line 98 of file TwoSampleDecomp.h.

Referenced by DrawTwoSamplesWithRatios(), RatioCalc(), and SaveTo().

Spectrum ana::TwoSampleDecomp::fFirstSig

Spectrum for first sample, signal.

Definition at line 96 of file TwoSampleDecomp.h.

Referenced by DrawTwoSamplesWithRatios(), RatioCalc(), and SaveTo().

const SpectrumLoaderBase& ana::TwoSampleDecomp::fLoaderData

Definition at line 93 of file TwoSampleDecomp.h.

Referenced by Decomp().

const SpectrumLoaderBase& ana::TwoSampleDecomp::fLoaderMC

Definition at line 92 of file TwoSampleDecomp.h.

Referenced by Decomp().

Spectrum ana::TwoSampleDecomp::fSecond

Spectrum of events that fail the firstSample cut.

Definition at line 105 of file TwoSampleDecomp.h.

Referenced by Decomp(), DrawSigBkgOverlay(), and SaveTo().

Spectrum ana::TwoSampleDecomp::fSecondBkg

Spectrum for second sample, background.

Definition at line 99 of file TwoSampleDecomp.h.

Referenced by DrawTwoSamplesWithRatios(), RatioCalc(), and SaveTo().

Spectrum ana::TwoSampleDecomp::fSecondSig

Spectrum for second sample, signal.

Definition at line 97 of file TwoSampleDecomp.h.

Referenced by DrawTwoSamplesWithRatios(), RatioCalc(), and SaveTo().

Spectrum ana::TwoSampleDecomp::fSig

mutable

Decomposed spectrum of signal events.

Definition at line 107 of file TwoSampleDecomp.h.

Referenced by Decomp(), DrawSigBkgOverlay(), NumuComponent(), and SaveTo().

TH1D* ana::TwoSampleDecomp::fSigRatio

mutable

Ratio of second/first spectra for signal.

Definition at line 111 of file TwoSampleDecomp.h.

Referenced by Decomp(), DrawRatios(), and RatioCalc().

Spectrum ana::TwoSampleDecomp::fTrueBkg

True backgroundh spectrum.

Definition at line 102 of file TwoSampleDecomp.h.

Referenced by DrawSigBkgOverlay(), and SaveTo().

Spectrum ana::TwoSampleDecomp::fTrueSig

True signal spectrum.

Definition at line 101 of file TwoSampleDecomp.h.

Referenced by DrawSigBkgOverlay(), and SaveTo().

---

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

• /cvmfs/nova-development.opensciencegrid.org/novasoft/releases/N21-02-25/CAFAna/Decomp/TwoSampleDecomp.h
• /cvmfs/nova-development.opensciencegrid.org/novasoft/releases/N21-02-25/CAFAna/Decomp/TwoSampleDecomp.cxx