ana::ExposureSource Class Reference

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

Inheritance diagram for ana::ExposureSource:

Classes
struct	CutDetails

Public Member Functions
	~ExposureSource ()
void	AddSink (const Cut &cut, IExposureSink *s)
void	HandleSpill (const caf::SRSpillProxy *spill) override
void	Finalize ()
	POT/livetime is not filled into sinks until this is called! More...
void	ReportExposures (double denom=0) const

Protected Attributes
std::vector< std::pair< Cut, CutDetails > >	fCuts
bool	fFinalized = false

Detailed Description

Definition at line 15 of file ExposureSource.h.

Constructor & Destructor Documentation

ana::ExposureSource::~ExposureSource

(

)

Definition at line 13 of file ExposureSource.cxx.

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

  {
    if(!fFinalized && !fCuts.empty()){
      std::cout << "ExposureSource destroyed without Finalize() being called" << std::endl;
      abort();
    }
  }

Member Function Documentation

void ana::ExposureSource::AddSink	(	const Cut &	cut,
		IExposureSink *	s
	)

Definition at line 22 of file ExposureSource.cxx.

References ana::assert(), fCuts, fFinalized, ana::_Cut< T >::ID(), it, and make_pair().

Referenced by ana::SpectrumLoaderBase::AddReweightableSpectrum(), and ana::SpectrumLoaderBase::AddSpectrum().

  {
    assert(!fFinalized);

    // It's more important that we can iterate fast in HandleSpill() than that
    // this function is fast.

    for(auto& it: fCuts){
      if(it.first.ID() == cut.ID()){
        it.second.sinks.push_back(s);
        return;
      }
    }

    fCuts.push_back(std::make_pair(cut, CutDetails()));
    fCuts.back().second.sinks.push_back(s);
  }

void ana::ExposureSource::Finalize

(

)

POT/livetime is not filled into sinks until this is called!

Definition at line 79 of file ExposureSource.cxx.

References om::cout, allTimeWatchdog::endl, fCuts, fFinalized, ana::IExposureSink::HandleLivetime(), ana::IExposureSink::HandlePOT(), and it.

Referenced by ana::SpectrumLoader::Go().

  {
    if(fFinalized){
      std::cout << "ExposureSource::Finalize() called twice!" << std::endl;
      abort();
    }

    for(auto& it: fCuts){
      for(IExposureSink* s: it.second.sinks){
        s->HandlePOT(it.second.pot);
        s->HandleLivetime(it.second.livetime);
      }
    }

    fFinalized = true;
  }

void ana::ExposureSource::HandleSpill ( const caf::SRSpillProxy *  spill )

overridevirtual

Implements ana::ISpillSink.

Definition at line 41 of file ExposureSource.cxx.

References ana::assert(), om::cout, cut, caf::Proxy< caf::SRSpill >::det, allTimeWatchdog::endl, fCuts, fFinalized, caf::Proxy< caf::SRSpill >::ismc, it, caf::kFARDET, livetime, ana::_Cut< T >::Livetime(), ana::_Cut< T >::POT(), pot, and caf::Proxy< caf::SRSpill >::spillpot.

Referenced by ana::SpectrumLoader::AccumulateExposures().

  {
    assert(!fFinalized);

    for(auto& it: fCuts){
      const Cut& cut = it.first;
      double pot = cut.POT(spill);
      // Not implemented -> fall back to usual counting
      if(pot < 0) pot = spill->spillpot;


      // FD MC is 1 neutrino-per-event, so the livetime count will be all
      // screwed up there. Just go with 0, ie "unknown".
      double livetime = 0;
      if(spill->det != caf::kFARDET || !spill->ismc) // TODO || fSource == kCosmic)
        livetime = cut.Livetime(spill);

      // Not implemented -> no good way to count
      if(livetime < 0){
        if(spill->det == caf::kFARDET && !spill->ismc){
          // People might well expect FD data to have a livetime. Warn.
          static bool once = true;
          if(once){
            once = false;
            std::cout << "WARNING: no way to compute livetime for FD data spectrum. If you want a livetime you need to be applying one of the cuts from TimingCuts.h or similar. You probably should be anyway to remove bad data near the spill ends." << std::endl;
          }
        }
        // Otherwise the user probably doesn't care. Just don't sum it into the
        // livetime
        livetime = 0;
      }

      it.second.pot += pot;
      it.second.livetime += livetime;
    }
  }

void ana::ExposureSource::ReportExposures

(

double

denom = 0

)

const

Definition at line 104 of file ExposureSource.cxx.

References om::cout, allTimeWatchdog::endl, fCuts, it, ana::POTPlusPercent(), and pots.

Referenced by ana::SpectrumLoader::Go().

  {
    std::cout << "Total exposure:" << std::endl;

    // Normally using a floating-point number as a key is a bad idea. But this
    // seems to work, and if it goes wrong all it breaks is this printout.
    std::map<double, int> livetimes;
    std::map<double, int> pots;
    for(auto it: fCuts){
      livetimes[it.second.livetime] += it.second.sinks.size();
      pots[it.second.pot] += it.second.sinks.size();
    }

    if(pots.size() == 1){
      std::cout << "  ";
      POTPlusPercent(std::cout, pots.begin()->first, denom);
      std::cout << std::endl;
    }
    else{
      for(auto it: pots){
        std::cout << "  ";
        POTPlusPercent(std::cout, pots.begin()->first, denom);
        std::cout << " for " << it.second << " spectra" << std::endl;
      }
    }

    if(livetimes.size() == 1){
      if(livetimes.begin()->first == 0){
        std::cout << "  Unknown livetime" << std::endl;
      }
      else{
        std::cout << "  " << livetimes.begin()->first << " seconds" << std::endl;
      }
    }
    else{
      for(auto it: livetimes){
        std::cout << "  " << it.first << " seconds for " << it.second << " spectra" << std::endl;
      }
    }
  }

Member Data Documentation

std::vector<std::pair<Cut, CutDetails> > ana::ExposureSource::fCuts

protected

Definition at line 39 of file ExposureSource.h.

Referenced by AddSink(), Finalize(), HandleSpill(), ReportExposures(), and ~ExposureSource().

bool ana::ExposureSource::fFinalized = false

protected

Definition at line 41 of file ExposureSource.h.

Referenced by AddSink(), Finalize(), HandleSpill(), and ~ExposureSource().

---

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

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