cmf::CovarianceBinUtility Class Reference

#include "/cvmfs/nova-development.opensciencegrid.org/novasoft/releases/N21-02-25/CovarianceMatrixFit/utilities/CovarianceBinUtility.h"

Public Member Functions
void	Initialize (fhicl::ParameterSet const &pset)
double const &	NuMuLowEdge ()
double const &	NuELowEdge (cmf::MetaData const &md)
double const &	NCLowEdge (cmf::DetType_t const &det)
double const &	SelectionLowEdge (cmf::MetaData const &md)
std::map< double, int > const &	NuMuHighEdges ()
std::map< double, int > const &	NuEHighEdges (cmf::DetType_t const &det)
std::map< double, int > const &	NCHighEdges (cmf::DetType_t const &det)
std::map< double, int > const &	SelectionHighEdges (cmf::MetaData const &md)
std::map< double, int > const &	SelectionHighEdges (long const &key)
void	SelectionHistBinning (cmf::MetaData const &md, std::vector< double > &bins)
int	KeyToOffset (long const &key, bool allSels=false)
int	EnergyToHistogramBin (double const &energy, cmf::MetaData const &md)
int	EnergyToBin (double const &energy, cmf::MetaData const &md)
double	BinToEnergy (int const &bin, bool allSels=false)
cmf::DetType_t	BinToDetector (int const &bin, bool allSels=false)
cmf::SelectionType_t	BinToSelection (int const &bin, bool allSels=false)
cmf::BeamType_t	BinToBeamType (int const &bin, bool allSels=false)
long	BinToKey (int const &bin, bool allSels=false)
size_t	TotalBins (bool allSels=false)
std::vector< int >	EmptyBins (TMatrixD mat)
TMatrixD	CompactMatrix (TMatrixD mat, std::vector< int > emptyBins)

Static Public Member Functions
static CovarianceBinUtility *	Instance ()

Private Member Functions
	CovarianceBinUtility ()
	~CovarianceBinUtility ()
std::map< double, int >	EnergyBinMap (std::vector< double > const &vec)
void	SetBinning (fhicl::ParameterSet const &fBinConfig, cmf::SelectionType_t const &selType)
void	FillOffsetAndBinMaps (std::map< long, int > &offsetMap, std::vector< long > &binMap, cmf::SelectionUtility::DetBeamSelSet const &detBeamSelSet)
int	NumSelectionBins (cmf::SelectionType_t const &selType, cmf::DetType_t const &detType)
int	KeyToOffset (long const &key, std::map< long, int > const &offsetMap)
std::map< int, double > const &	SelectionBinToEnergy (long key)

Static Private Member Functions
static void	FillBinToEnergyMap (std::map< double, int > const &selectionHighEdges, std::map< int, double > &selectionBinToEnergy, double lowestBound)

Private Attributes
double	fNuMuLowEdge
	low side of the first numu bin More...
std::map< double, int >	fNuMuHighEdges
	high side of the numu bins More...
double	fNuENDLowEdge
	low side of the first nue bin in ND More...
std::map< double, int >	fNuENDHighEdges
	high side of the nue bins in ND More...
double	fNuEFDLowEdge
	low side of the first nune bin in FD More...
std::map< double, int >	fNuEFDHighEdges
	high side of the nue bins in FD More...
double	fNuEPeriphLowEdge
	low side of the first nue bin in FD More...
std::map< double, int >	fNuEPeriphHighEdges
	high side of the nue bins in FD More...
double	fNCLowEdge
	low side of the first nc bin More...
std::map< double, int >	fNCHighEdges
	high side of the nc bins More...
double	fNCNDLowEdge
	low side of the first nc bin in ND More...
std::map< double, int >	fNCNDHighEdges
	high side of the nc bins in ND (2020ana) More...
double	fNCFDLowEdge
	low side of the first nc bin in FD More...
std::map< double, int >	fNCFDHighEdges
	high side of the nc bins in FD (2020ana) More...
std::map< int, double >	fNuMuBinToEnergy
	map energy bin to central value More...
std::map< int, double >	fNuENDBinToEnergy
	map energy bin to central value More...
std::map< int, double >	fNuEFDBinToEnergy
	map energy bin to central value More...
std::map< int, double >	fNuEPeriphBinToEnergy
	map energy bin to central value More...
std::map< int, double >	fNCNDBinToEnergy
	map energy bin to central value (2020ana) More...
std::map< int, double >	fNCFDBinToEnergy
	map energy bin to central value (2020ana) More...
std::map< long, int >	fOffsetMap
	map of keys to offset for different samples More...
std::map< long, int >	fAllSelsOffsetMap
	map of keys to offset for different samples More...
std::vector< long >	fBinToKey
	map bin numbers to keys for different samples More...
std::vector< long >	fAllSelsBinToKey
	map bin numbers to keys for different samples More...

Detailed Description

Definition at line 22 of file CovarianceBinUtility.h.

Constructor & Destructor Documentation

cmf::CovarianceBinUtility::CovarianceBinUtility ( )

private

Definition at line 26 of file CovarianceBinUtility.cxx.

Referenced by Instance(), and NuMuHighEdges().

  {

  }

cmf::CovarianceBinUtility::~CovarianceBinUtility ( )

private

Definition at line 223 of file CovarianceBinUtility.cxx.

References fAllSelsBinToKey, fAllSelsOffsetMap, fBinToKey, fNCFDBinToEnergy, fNCFDHighEdges, fNCHighEdges, fNCNDBinToEnergy, fNCNDHighEdges, fNuEFDBinToEnergy, fNuEFDHighEdges, fNuENDBinToEnergy, fNuENDHighEdges, fNuEPeriphBinToEnergy, fNuEPeriphHighEdges, fNuMuBinToEnergy, fNuMuHighEdges, and fOffsetMap.

Referenced by NuMuHighEdges().

  {
    fNuMuHighEdges   .clear();
    fNuENDHighEdges  .clear();
    fNuEFDHighEdges  .clear();
    fNCHighEdges     .clear();
    fNCNDHighEdges   .clear();
    fNCFDHighEdges   .clear();
    fAllSelsOffsetMap.clear();
    fOffsetMap       .clear();
    fBinToKey        .clear();
    fAllSelsBinToKey .clear();

    fNuEPeriphHighEdges  .clear();
    fNuMuBinToEnergy     .clear();
    fNuENDBinToEnergy    .clear();
    fNuEFDBinToEnergy    .clear();
    fNuEPeriphBinToEnergy.clear();
    fNCNDBinToEnergy     .clear();
    fNCFDBinToEnergy     .clear();
 
  }

Member Function Documentation

cmf::BeamType_t cmf::CovarianceBinUtility::BinToBeamType	(	int const &	bin,
		bool	allSels = false
	)

Definition at line 536 of file CovarianceBinUtility.cxx.

References BinToKey(), and cmf::KeyToBeamType().

Referenced by NuMuHighEdges().

  {
    return cmf::KeyToBeamType(this->BinToKey(bin, allSels));
  }

cmf::DetType_t cmf::CovarianceBinUtility::BinToDetector	(	int const &	bin,
		bool	allSels = false
	)

Definition at line 522 of file CovarianceBinUtility.cxx.

References BinToKey(), and cmf::KeyToDetectorType().

Referenced by NuMuHighEdges(), and cmf::RandomUniverses::writeResults().

  {
    return cmf::KeyToDetectorType(this->BinToKey(bin, allSels));
  }

double cmf::CovarianceBinUtility::BinToEnergy	(	int const &	bin,
		bool	allSels = false
	)

Definition at line 551 of file CovarianceBinUtility.cxx.

References BinToKey(), findDuplicateFiles::key, KeyToOffset(), cmf::KeyToString(), LOG_DEBUG, PandAna.reco_validation.add_data::offset, and SelectionBinToEnergy().

Referenced by cmf::CovarianceMatrixMaker::FillSpectrum(), cmf::PlotUtilities::MakeEnergySpectraFromBins(), and NuMuHighEdges().

  {
    auto key    = this->BinToKey   (bin, allSels);
    int  offset = this->KeyToOffset(key, allSels);

    auto const& binToEnergy = this->SelectionBinToEnergy(key);

    auto const btoEItr = binToEnergy.find(bin - offset);

    if(btoEItr == binToEnergy.end())
      throw cet::exception("CovarianceBinUtility")
          << "bin "
          << bin
          << " offset "
          << offset
          << " "
          << bin - offset
          << " is not in bin to energy map for key "
          << key
          << " "
          << cmf::KeyToString(key);

    LOG_DEBUG("CovarianceBinUtility")
        << "BinToEnergy: "
        << cmf::KeyToString(key)
        << " bin is "
        << bin
        << " offset "
        << offset
        << " energy "
        << btoEItr->second;

    return btoEItr->second;
  }

long cmf::CovarianceBinUtility::BinToKey	(	int const &	bin,
		bool	allSels = false
	)

Definition at line 543 of file CovarianceBinUtility.cxx.

References bin, fAllSelsBinToKey, and fBinToKey.

Referenced by BinToBeamType(), BinToDetector(), BinToEnergy(), BinToSelection(), cmf::CovarianceMatrixMaker::FillSpectrum(), cmf::ChiSqrCalculator::FindReducedMatrixBin(), and NuMuHighEdges().

  {
    if(allSels) return fAllSelsBinToKey[bin];
    return fBinToKey[bin];
  }

cmf::SelectionType_t cmf::CovarianceBinUtility::BinToSelection	(	int const &	bin,
		bool	allSels = false
	)

Definition at line 529 of file CovarianceBinUtility.cxx.

References BinToKey(), and cmf::KeyToSelectionType().

Referenced by NuMuHighEdges().

  {
    return cmf::KeyToSelectionType(this->BinToKey(bin, allSels));
  }

TMatrixD cmf::CovarianceBinUtility::CompactMatrix	(	TMatrixD	mat,
		std::vector< int >	emptyBins
	)

Definition at line 608 of file CovarianceBinUtility.cxx.

References MECModelEnuComparisons::i, and calib::j.

Referenced by NuMuHighEdges().

                                                                                    {

    int nRowsAndColumns = mat.GetNcols() - emptyBins.size();
    TMatrixD compactMatrix(nRowsAndColumns, nRowsAndColumns);
    TMatrixD semiCompactMatrix(nRowsAndColumns, mat.GetNcols());

    // correct in i
    int itest = 0;
    for (int i = 0; i < nRowsAndColumns; ++i){
      
      bool isXEmptyBin = false;
      for (size_t ieb = 0; ieb < emptyBins.size(); ++ieb){
        if (emptyBins.at(ieb) == itest)
          isXEmptyBin = true;
      }

      for (int j = 0; j < mat.GetNcols(); ++j){
        semiCompactMatrix(i,j) = mat(itest, j);
      }

      if (isXEmptyBin)
        --i;

      ++itest;

    }

    // correct in i
    int jtest = 0;
    for (int j = 0; j < nRowsAndColumns; ++j){
      
      bool isYEmptyBin = false;
      for (size_t ieb = 0; ieb < emptyBins.size(); ++ieb){
        if (emptyBins.at(ieb) == jtest)
          isYEmptyBin = true;
      }

      for (int i = 0; i < semiCompactMatrix.GetNrows(); ++i){
        compactMatrix(i,j) = semiCompactMatrix(i, jtest);
      }

      if (isYEmptyBin)
        --j;

      ++jtest;

    }

    return compactMatrix;

  }

std::vector< int > cmf::CovarianceBinUtility::EmptyBins

(

TMatrixD

mat

)

Definition at line 588 of file CovarianceBinUtility.cxx.

References std::abs(), MECModelEnuComparisons::i, and calib::j.

Referenced by NuMuHighEdges().

                                                            {

    std::vector<int> emptyBins;

    for (int i = 0; i < mat.GetNrows(); i++){

      // for every column, check every row is above some
      // tolerance. If it's not, count it as an empty bin
      double rowsum = 0;
      for (int j = 0; j < mat.GetNcols(); j++)
        rowsum += std::abs(mat(i,j)); 

      if (rowsum == 0)
        emptyBins.push_back(i);
    }

    return emptyBins;
  }

std::map< double, int > cmf::CovarianceBinUtility::EnergyBinMap ( std::vector< double > const &  vec )

private

Definition at line 107 of file CovarianceBinUtility.cxx.

References MECModelEnuComparisons::i.

Referenced by NuMuHighEdges(), and SetBinning().

  {
    std::map<double, int> ebMap;
    for (size_t i = 0; i < vec.size(); ++i)
      ebMap[vec.at(i)] = i;

    return ebMap;
  }

int cmf::CovarianceBinUtility::EnergyToBin	(	double const &	energy,
		cmf::MetaData const &	md
	)

Definition at line 480 of file CovarianceBinUtility.cxx.

References SelectionHighEdges().

Referenced by EnergyToHistogramBin(), cmf::ShifterAndWeighter::FileSystWeight(), cmf::FillSpectrumFromSingleList(), and NuMuHighEdges().

  {
    auto const& edgeMap = this->SelectionHighEdges(md);
    
    auto const& itr = edgeMap.upper_bound(energy);
    if(itr != edgeMap.end()) return itr->second;
    
    return 0;
  }

int cmf::CovarianceBinUtility::EnergyToHistogramBin	(	double const &	energy,
		cmf::MetaData const &	md
	)

Definition at line 473 of file CovarianceBinUtility.cxx.

References EnergyToBin().

Referenced by NuMuHighEdges().

  {
    return this->EnergyToBin(energy, md) + 1;
  }

void cmf::CovarianceBinUtility::FillBinToEnergyMap ( std::map< double, int > const &  selectionHighEdges, std::map< int, double > &  selectionBinToEnergy, double  lowestBound  )

staticprivate

Definition at line 247 of file CovarianceBinUtility.cxx.

References LOG_DEBUG.

Referenced by Initialize(), and NuMuHighEdges().

  {
    // fill the bin to energy maps for fast look up
    double lowerBound = lowestBound;
    for(auto const& itr : selectionHighEdges){
      selectionBinToEnergy[itr.second] = lowerBound + (itr.first - lowerBound) * 0.5;
      
      LOG_DEBUG("CovarianceBinUtility")
          << "FillBinToEnergyMap energy "
          << itr.first
          << " bin "
          << itr.second
          << " bin to energy "
          << selectionBinToEnergy.find(itr.second)->second;

      lowerBound = itr.first;
    }
  }

void cmf::CovarianceBinUtility::FillOffsetAndBinMaps ( std::map< long, int > &  offsetMap, std::vector< long > &  binMap, cmf::SelectionUtility::DetBeamSelSet const &  detBeamSelSet  )

private

Definition at line 55 of file CovarianceBinUtility.cxx.

References b, bins, cmf::DetectorBeamSelectionTypesToKey(), findDuplicateFiles::key, cmf::KeyToString(), cmf::kNEARDET, cmf::kNuESelectionPeripheral, LOG_DEBUG, NumSelectionBins(), and PandAna.reco_validation.add_data::offset.

Referenced by Initialize(), and NuMuHighEdges().

   {
     long key;
     int  offset = 0;
     int  bins;

     binMap   .clear();
     offsetMap.clear();
     
     // loop over detectors, beam type and selection to fill the maps
     for(auto const& dbsItr : detBeamSelSet){
       for(auto const& selItr : dbsItr.Selections() ){

         // there is no MidPID selection for nues any more and no nue peripheral
         // selection in the nd
         if(dbsItr.Detector() == cmf::kNEARDET &&
            selItr == cmf::kNuESelectionPeripheral)
           continue;

         // set the file type to 0 so that it does not contribute to the
         // key value in this case
         key = DetectorBeamSelectionTypesToKey(dbsItr.Detector(),
                                               dbsItr.BeamType(),
                                               selItr);

         offsetMap[key] = offset;

         bins = this->NumSelectionBins(selItr, dbsItr.Detector());
         for(int b = offset; b < bins + offset; ++b){
           binMap.emplace_back(key);
           LOG_DEBUG("CovarianceBinUtility")
               << cmf::KeyToString(key)
               << " key "
               << key
               << " offset "
               << offsetMap.find(key)->second
               << " "
               << b
               << " "
               << binMap[b];
         }

         offset += bins;
       } // end loop over selection types
     } // end loop over input DetBeamSels

   }

void cmf::CovarianceBinUtility::Initialize

(

fhicl::ParameterSet const &

pset

)

Definition at line 151 of file CovarianceBinUtility.cxx.

References fAllSelsBinToKey, fAllSelsOffsetMap, fBinToKey, FillBinToEnergyMap(), FillOffsetAndBinMaps(), fNCFDBinToEnergy, fNCFDHighEdges, fNCFDLowEdge, fNCNDBinToEnergy, fNCNDHighEdges, fNCNDLowEdge, fNuEFDBinToEnergy, fNuEFDHighEdges, fNuEFDLowEdge, fNuENDBinToEnergy, fNuENDHighEdges, fNuENDLowEdge, fNuEPeriphBinToEnergy, fNuEPeriphHighEdges, fNuEPeriphLowEdge, fNuMuBinToEnergy, fNuMuHighEdges, fNuMuLowEdge, fOffsetMap, fhicl::ParameterSet::get(), cmf::SelectionUtility::Instance(), cmf::kNCSelection, cmf::kNuESelection, cmf::kNuMuSelection, LOG_DEBUG, LOG_VERBATIM, and SetBinning().

Referenced by cmf::PredictionLibraryMaker::PredictionLibraryMaker(), cmf::CMFSpectraMaker::reconfigure(), cmf::CovarianceMatrixFitter::reconfigure(), cmf::CovarianceMatrixMaker::reconfigure(), cmf::RandomUniverses::reconfigure(), cmf::FitFeldmanCousinsPoint::reconfigure(), and cmf::ContourFromLibrary::reconfigure().

  {

    // grab configurations from parameter set
    auto fNuMuEdgeConfig = pset.get< fhicl::ParameterSet >("numu_bin_edges");
    auto fNuEEdgeConfig  = pset.get< fhicl::ParameterSet >("nue_bin_edges");
    auto fNCEdgeConfig   = pset.get< fhicl::ParameterSet >("nc_bin_edges");

    // using the pre-existing selection enum, doesn't matter
    // since we just need to differentiate between numu nue and nc
    this->SetBinning(fNuMuEdgeConfig, cmf::kNuMuSelection);
    this->SetBinning(fNuEEdgeConfig , cmf::kNuESelection);
    this->SetBinning(fNCEdgeConfig  , cmf::kNCSelection);

    LOG_DEBUG("CovarianceBinUtility")
      << "\nnum NuMuBins:     "  << fNuMuHighEdges.size()      
      << "\nnum NuENDBins:    "  << fNuENDHighEdges.size()     
      << "\nnum NuEFDBins:    "  << fNuEFDHighEdges.size()     
      << "\nnum NuEPeriphBins:"  << fNuEPeriphHighEdges.size() 
      << "\nnum NCNDBins:     "  << fNCNDHighEdges.size()      
      << "\nnum NCFDBins:     "  << fNCFDHighEdges.size();

    FillBinToEnergyMap(fNuMuHighEdges,
                       fNuMuBinToEnergy,
                       fNuMuLowEdge);

    FillBinToEnergyMap(fNuENDHighEdges,
                       fNuENDBinToEnergy,
                       fNuENDLowEdge);

    FillBinToEnergyMap(fNuEFDHighEdges,
                       fNuEFDBinToEnergy,
                       fNuEFDLowEdge);

    FillBinToEnergyMap(fNuEPeriphHighEdges,
                       fNuEPeriphBinToEnergy,
                       fNuEPeriphLowEdge);

    FillBinToEnergyMap(fNCNDHighEdges,
                       fNCNDBinToEnergy,
                       fNCNDLowEdge);

    FillBinToEnergyMap(fNCFDHighEdges,
                       fNCFDBinToEnergy,
                       fNCFDLowEdge);

    // call the function to check that the mapping from 
    // energy to bin and back again is consistent
    //this->CheckBinToEnergyMapping();

    LOG_DEBUG("CovarianceBinUtility")
    << "Fill All Selection Offset and Bin maps";

    this->FillOffsetAndBinMaps(fAllSelsOffsetMap,
                               fAllSelsBinToKey,
                               cmf::SelectionUtility::Instance()->AllAvailSelections());

    LOG_DEBUG("CovarianceBinUtility")
        << "Fill run time Selection Offset and Bin maps";

    this->FillOffsetAndBinMaps(fOffsetMap,
                               fBinToKey,
                               cmf::SelectionUtility::Instance()->SelectionsToUse());

    LOG_VERBATIM("CovarianceBinUtility")
      << "Using "
      << fBinToKey.size()
      << " bins for analysis out of a possible "
      << fAllSelsBinToKey.size();
  }

CovarianceBinUtility * cmf::CovarianceBinUtility::Instance ( )

static

Definition at line 18 of file CovarianceBinUtility.cxx.

References CovarianceBinUtility(), and cmf::gCovBinUtil.

Referenced by cmf::CovarianceMatrixMaker::beginJob(), cmf::RandomUniverses::beginJob(), cmf::ChiSqrCalculator::ChiSqrCalculator(), cmf::CMFDecorrelator::CMFDecorrelator(), cmf::CovarianceFitHelper::CovarianceFitHelper(), cmf::ContourFromLibrary::endJob(), cmf::ShifterAndWeighter::FileSystWeight(), cmf::CovarianceMatrixMaker::FillSpectrum(), cmf::FillSpectrum(), cmf::FillSpectrumFromSingleList(), cmf::ChiSqrCalculator::FindReducedMatrixBin(), cmf::ChiSqrCalculator::InitializeCovarianceMatrix(), cmf::CovarianceFitHelper::Make1DSpectra(), cmf::CovarianceFitHelper::MakeDataMCCanv(), cmf::PlotUtilities::MakeEnergySpectraFromBins(), cmf::SpectrumPredictionMaker::MakeHiddenParameterHistograms(), makeInitialDecorrelatedErrorPlot(), cmf::CovarianceFitHelper::MakeResultPlots(), cmf::RandomUniverseUtility::PoissonSpectrumBinByBin(), cmf::PredictionLibraryMaker::PredictionLibraryMaker(), cmf::CMFSpectraMaker::reconfigure(), cmf::CovarianceMatrixFitter::reconfigure(), cmf::CovarianceMatrixMaker::reconfigure(), cmf::RandomUniverses::reconfigure(), cmf::FitFeldmanCousinsPoint::reconfigure(), cmf::ContourFromLibrary::reconfigure(), cmf::CMFDecorrelator::SetComponentNameAndOffsets(), cmf::CMFDecorrelator::SetTemplateMatrix(), cmf::RandomUniverses::writeResults(), and cmf::FitFeldmanCousinsPoint::writeResults().

  {
    if(gCovBinUtil == nullptr) gCovBinUtil = new CovarianceBinUtility();

    return gCovBinUtil;
  }

int cmf::CovarianceBinUtility::KeyToOffset	(	long const &	key,
		bool	allSels = false
	)

Definition at line 446 of file CovarianceBinUtility.cxx.

References fAllSelsOffsetMap, and fOffsetMap.

Referenced by BinToEnergy(), cmf::FillSpectrumFromSingleList(), cmf::PlotUtilities::MakeEnergySpectraFromBins(), makeInitialDecorrelatedErrorPlot(), NuMuHighEdges(), and cmf::CMFDecorrelator::SetComponentNameAndOffsets().

  {
    if(allSels) return this->KeyToOffset(key, fAllSelsOffsetMap);
    return this->KeyToOffset(key, fOffsetMap);
  }

int cmf::CovarianceBinUtility::KeyToOffset ( long const &  key, std::map< long, int > const &  offsetMap  )

private

Definition at line 454 of file CovarianceBinUtility.cxx.

References cmf::cBeamType_Strings, cmf::cDetType_Strings, cmf::cSelectionType_Strings, cmf::KeyToBeamType(), cmf::KeyToDetectorType(), and cmf::KeyToSelectionType().

  {

    if(offsetMap.count(key) < 1)
      throw cet::exception("CovarianceBinUtility")
          << "could not find offset in map for key "
          << key
          << " " 
          << cmf::cBeamType_Strings[cmf::KeyToBeamType(key)]
          << " "
          << cmf::cDetType_Strings[cmf::KeyToDetectorType(key)]
          << " " 
          << cmf::cSelectionType_Strings[cmf::KeyToSelectionType(key)];

    return offsetMap.find(key)->second;
  }

std::map< double, int > const & cmf::CovarianceBinUtility::NCHighEdges

(

cmf::DetType_t const &

det

)

Definition at line 396 of file CovarianceBinUtility.cxx.

References fNCFDHighEdges, fNCNDHighEdges, and cmf::kNEARDET.

Referenced by NuMuHighEdges(), and SelectionHighEdges().

  {
    if(det == cmf::kNEARDET) return fNCNDHighEdges;
    
    return fNCFDHighEdges;
  }

double const & cmf::CovarianceBinUtility::NCLowEdge

(

cmf::DetType_t const &

det

)

Definition at line 380 of file CovarianceBinUtility.cxx.

References fNCFDLowEdge, fNCNDLowEdge, and cmf::kNEARDET.

Referenced by NuMuLowEdge(), and SelectionLowEdge().

  {
    if(det == cmf::kNEARDET) return fNCNDLowEdge;
    
    return fNCFDLowEdge;
  }

std::map< double, int > const & cmf::CovarianceBinUtility::NuEHighEdges

(

cmf::DetType_t const &

det

)

Definition at line 388 of file CovarianceBinUtility.cxx.

References fNuEFDHighEdges, fNuENDHighEdges, and cmf::kNEARDET.

Referenced by NuMuHighEdges(), and SelectionHighEdges().

  {
    if(det == cmf::kNEARDET) return fNuENDHighEdges;
    
    return fNuEFDHighEdges;
  }

double const & cmf::CovarianceBinUtility::NuELowEdge

(

cmf::MetaData const &

md

)

Definition at line 371 of file CovarianceBinUtility.cxx.

References cmf::MetaData::detector, fNuEFDLowEdge, fNuENDLowEdge, fNuEPeriphLowEdge, cmf::kNEARDET, cmf::kNuESelectionPeripheral, and cmf::MetaData::selectionType.

Referenced by NuMuLowEdge(), and SelectionLowEdge().

  {
    if(md.detector == cmf::kNEARDET) return fNuENDLowEdge;
    if(md.selectionType == cmf::kNuESelectionPeripheral) return fNuEPeriphLowEdge;

    return fNuEFDLowEdge;
  }

int cmf::CovarianceBinUtility::NumSelectionBins ( cmf::SelectionType_t const &  selType, cmf::DetType_t const &  detType  )

private

Definition at line 32 of file CovarianceBinUtility.cxx.

References fNCFDHighEdges, fNCNDHighEdges, fNuEFDHighEdges, fNuENDHighEdges, fNuEPeriphHighEdges, fNuMuHighEdges, cmf::IsNuESelected(), cmf::IsNuMuSelected(), cmf::kNCSelection, cmf::kNEARDET, and cmf::kNuESelectionPeripheral.

Referenced by FillOffsetAndBinMaps(), and NuMuHighEdges().

  {
    if(cmf::IsNuESelected(selType)){
      if(selType == cmf::kNuESelectionPeripheral) return fNuEPeriphHighEdges.size();
      return (detType == cmf::kNEARDET) ? fNuENDHighEdges.size() : fNuEFDHighEdges.size();
    }
    else if(cmf::IsNuMuSelected(selType)){
      return fNuMuHighEdges.size();
    }
    else if(selType == cmf::kNCSelection){
      return (detType == cmf::kNEARDET) ? fNCNDHighEdges.size() : fNCFDHighEdges.size();
    }

    throw cet::exception("CovarianceBinUtility")
        << "selection type "
        << selType
        << " is not recognized";

    return 0;
  }

std::map<double, int> const& cmf::CovarianceBinUtility::NuMuHighEdges ( )

inline

Definition at line 34 of file CovarianceBinUtility.h.

References bin, bins, BinToBeamType(), BinToDetector(), BinToEnergy(), BinToKey(), BinToSelection(), CompactMatrix(), CovarianceBinUtility(), EmptyBins(), energy, EnergyBinMap(), EnergyToBin(), EnergyToHistogramBin(), FillBinToEnergyMap(), FillOffsetAndBinMaps(), fNuMuHighEdges, findDuplicateFiles::key, KeyToOffset(), NCHighEdges(), NuEHighEdges(), NumSelectionBins(), SelectionBinToEnergy(), SelectionHighEdges(), SelectionHistBinning(), SetBinning(), TotalBins(), and ~CovarianceBinUtility().

Referenced by SelectionHighEdges().

{ return fNuMuHighEdges; }

double const& cmf::CovarianceBinUtility::NuMuLowEdge ( )

inline

Definition at line 30 of file CovarianceBinUtility.h.

References fillBadChanDBTables::det, fNuMuLowEdge, fetch_tb_beamline_files::md, NCLowEdge(), NuELowEdge(), and SelectionLowEdge().

Referenced by SelectionLowEdge().

{ return fNuMuLowEdge; }

std::map< int, double > const & cmf::CovarianceBinUtility::SelectionBinToEnergy ( long  key )

private

Definition at line 492 of file CovarianceBinUtility.cxx.

References cmf::cDetType_Strings, cmf::cSelectionType_Strings, fillBadChanDBTables::det, fNCFDBinToEnergy, fNCNDBinToEnergy, fNuEFDBinToEnergy, fNuENDBinToEnergy, fNuEPeriphBinToEnergy, fNuMuBinToEnergy, cmf::IsNuESelected(), cmf::IsNuMuSelected(), cmf::KeyToDetectorType(), cmf::KeyToSelectionType(), cmf::kNEARDET, cmf::kNuESelectionPeripheral, cmf::kUnknownDet, and cmf::kUnknownSelection.

Referenced by BinToEnergy(), and NuMuHighEdges().

  {
    auto sel = cmf::KeyToSelectionType(key);
    auto det = cmf::KeyToDetectorType(key);

    if(sel == cmf::kUnknownSelection ||
       det == cmf::kUnknownDet )
      throw cet::exception("Unknown Key")
      << "key "
      << key
      << "corresponds to "
      << cmf::cSelectionType_Strings[sel]
      << " and "
      << cmf::cDetType_Strings[det];

    if(cmf::IsNuESelected(sel)){
      if(sel == cmf::kNuESelectionPeripheral)  return fNuEPeriphBinToEnergy;
      if(det == cmf::kNEARDET) return fNuENDBinToEnergy;
      return fNuEFDBinToEnergy;
    }
    else if(cmf::IsNuMuSelected(sel)){
      return fNuMuBinToEnergy;
    }
    else{
      if (det == cmf::kNEARDET) return fNCNDBinToEnergy;
      else return fNCFDBinToEnergy;
    }
  }

std::map< double, int > const & cmf::CovarianceBinUtility::SelectionHighEdges

(

cmf::MetaData const &

md

)

Definition at line 404 of file CovarianceBinUtility.cxx.

References cmf::MetaData::detector, fNuEPeriphHighEdges, cmf::MetaData::IsNuESelected(), cmf::MetaData::IsNuMuSelected(), cmf::kNuESelectionPeripheral, NCHighEdges(), NuEHighEdges(), NuMuHighEdges(), and cmf::MetaData::selectionType.

Referenced by cmf::CovarianceMatrixMaker::beginJob(), EnergyToBin(), cmf::PlotUtilities::MakeEnergySpectraFromBins(), NuMuHighEdges(), SelectionHistBinning(), and cmf::CMFDecorrelator::SetTemplateMatrix().

  {
    if(md.IsNuMuSelected())
      return this->NuMuHighEdges();
    else if(md.IsNuESelected()){
      if(md.selectionType == cmf::kNuESelectionPeripheral) return fNuEPeriphHighEdges;
      return this->NuEHighEdges(md.detector);
    }
    
    return this->NCHighEdges(md.detector);
  }

std::map< double, int > const & cmf::CovarianceBinUtility::SelectionHighEdges

(

long const &

key

)

Definition at line 417 of file CovarianceBinUtility.cxx.

References fillBadChanDBTables::det, fNuEPeriphHighEdges, cmf::IsNuESelected(), cmf::IsNuMuSelected(), cmf::KeyToDetectorType(), cmf::KeyToSelectionType(), cmf::kNuESelectionPeripheral, NCHighEdges(), NuEHighEdges(), and NuMuHighEdges().

  {
    auto sel = cmf::KeyToSelectionType(key);
    auto det = cmf::KeyToDetectorType(key);

    if(cmf::IsNuMuSelected(sel))
      return this->NuMuHighEdges();
    else if(cmf::IsNuESelected(sel)){
      if(sel == cmf::kNuESelectionPeripheral) return fNuEPeriphHighEdges;
      return this->NuEHighEdges(det);
    }

    return this->NCHighEdges(det);
  }

void cmf::CovarianceBinUtility::SelectionHistBinning	(	cmf::MetaData const &	md,
		std::vector< double > &	bins
	)

Definition at line 433 of file CovarianceBinUtility.cxx.

References SelectionHighEdges(), SelectionLowEdge(), and std::swap().

Referenced by cmf::CovarianceFitHelper::Make1DSpectra(), cmf::PlotUtilities::MakeEnergySpectraFromBins(), and NuMuHighEdges().

  {
    std::vector<double> tmpbins;
    auto const& highEdgeMap = this->SelectionHighEdges(md);
    tmpbins.push_back(this->SelectionLowEdge(md));

    for(auto const& itr : highEdgeMap) tmpbins.push_back(itr.first);

    std::swap(tmpbins, bins);
  }

double const & cmf::CovarianceBinUtility::SelectionLowEdge

(

cmf::MetaData const &

md

)

Definition at line 359 of file CovarianceBinUtility.cxx.

References cmf::MetaData::detector, cmf::MetaData::IsNuESelected(), cmf::MetaData::IsNuMuSelected(), NCLowEdge(), NuELowEdge(), and NuMuLowEdge().

Referenced by cmf::CovarianceMatrixMaker::beginJob(), NuMuLowEdge(), and SelectionHistBinning().

  {
    if(md.IsNuMuSelected())
      return this->NuMuLowEdge();
    else if(md.IsNuESelected()){
      return this->NuELowEdge(md);
    }
    
    return this->NCLowEdge(md.detector);
  }

void cmf::CovarianceBinUtility::SetBinning ( fhicl::ParameterSet const &  fBinConfig, cmf::SelectionType_t const &  selType  )

private

Definition at line 118 of file CovarianceBinUtility.cxx.

References EnergyBinMap(), fNCFDHighEdges, fNCFDLowEdge, fNCNDHighEdges, fNCNDLowEdge, fNuEFDHighEdges, fNuEFDLowEdge, fNuENDHighEdges, fNuENDLowEdge, fNuEPeriphHighEdges, fNuEPeriphLowEdge, fNuMuHighEdges, fNuMuLowEdge, fhicl::ParameterSet::get(), cmf::IsNuESelected(), and cmf::IsNuMuSelected().

Referenced by Initialize(), and NuMuHighEdges().

  {
    auto ndConfig = fBinConfig.get<fhicl::ParameterSet>("nd");
    auto fdConfig = fBinConfig.get<fhicl::ParameterSet>("fd");

    // numu selection has the same binning in nd and fd
    // TODO if the numu selection is ever not the same in the detectors we
    // need to update this
    if(cmf::IsNuMuSelected(selType)){
      fNuMuLowEdge   = ndConfig.get<double>("low_edge");
      fNuMuHighEdges = this->EnergyBinMap(ndConfig.get<std::vector<double>>("high_edges"));
    }
    else if (cmf::IsNuESelected(selType)){
      fNuENDLowEdge   = ndConfig.get<double>("low_edge");
      fNuENDHighEdges = this->EnergyBinMap(ndConfig.get<std::vector<double>>("high_edges"));
      fNuEFDLowEdge   = fdConfig.get<double>("low_edge");
      fNuEFDHighEdges = this->EnergyBinMap(fdConfig.get<std::vector<double>>("high_edges"));

      // nue selection also has peripheral sample at the FD
      auto peConfig       = fBinConfig.get<fhicl::ParameterSet>("peripheral");
      fNuEPeriphLowEdge   = peConfig.get<double>("low_edge");
      fNuEPeriphHighEdges = this->EnergyBinMap(peConfig.get<std::vector<double>>("high_edges"));
    }
    else {
      fNCNDLowEdge   = ndConfig.get<double>("low_edge");
      fNCNDHighEdges = this->EnergyBinMap(ndConfig.get<std::vector<double>>("high_edges"));
      fNCFDLowEdge   = fdConfig.get<double>("low_edge");
      fNCFDHighEdges = this->EnergyBinMap(fdConfig.get<std::vector<double>>("high_edges"));
    }

  }

size_t cmf::CovarianceBinUtility::TotalBins

(

bool

allSels = false

)

Definition at line 330 of file CovarianceBinUtility.cxx.

References fAllSelsBinToKey, fBinToKey, fNCFDHighEdges, fNCNDHighEdges, fNuEFDHighEdges, fNuENDHighEdges, fNuEPeriphHighEdges, fNuMuHighEdges, and LOG_DEBUG.

Referenced by cmf::CovarianceMatrixMaker::beginJob(), cmf::RandomUniverses::beginJob(), cmf::ChiSqrCalculator::ChiSqrCalculator(), cmf::CMFDecorrelator::CMFDecorrelator(), cmf::ContourFromLibrary::endJob(), cmf::ChiSqrCalculator::InitializeCovarianceMatrix(), cmf::SpectrumPredictionMaker::MakeHiddenParameterHistograms(), cmf::CovarianceFitHelper::MakeResultPlots(), NuMuHighEdges(), and cmf::FitFeldmanCousinsPoint::writeResults().

  {

    LOG_DEBUG("CovarianceBinUtility")
      << " there are numu (each quartile): "
      << fNuMuHighEdges.size()
      << " nue low: "
      << fNuENDHighEdges.size()
      << " nue high: "
      << fNuEFDHighEdges.size()
      << " nue peripheral: "
      << fNuEPeriphHighEdges.size()
      << " nc ND: "
      << fNCNDHighEdges.size()
      << " nc FD: "
      << fNCFDHighEdges.size()
      << " total bins for all selections: "
      << fAllSelsBinToKey.size()
      << " total bins for run time selections: "
      << fBinToKey.size()
      << " "
      << allSels;
    
    if(allSels) return fAllSelsBinToKey.size();

    return fBinToKey.size();    
  }

Member Data Documentation

std::vector<long> cmf::CovarianceBinUtility::fAllSelsBinToKey

private

map bin numbers to keys for different samples

Definition at line 120 of file CovarianceBinUtility.h.

Referenced by BinToKey(), Initialize(), TotalBins(), and ~CovarianceBinUtility().

std::map<long, int> cmf::CovarianceBinUtility::fAllSelsOffsetMap

private

map of keys to offset for different samples

Definition at line 118 of file CovarianceBinUtility.h.

Referenced by Initialize(), KeyToOffset(), and ~CovarianceBinUtility().

std::vector<long> cmf::CovarianceBinUtility::fBinToKey

private

map bin numbers to keys for different samples

Definition at line 119 of file CovarianceBinUtility.h.

Referenced by BinToKey(), Initialize(), TotalBins(), and ~CovarianceBinUtility().

std::map<int, double> cmf::CovarianceBinUtility::fNCFDBinToEnergy

private

map energy bin to central value (2020ana)

Definition at line 116 of file CovarianceBinUtility.h.

Referenced by Initialize(), SelectionBinToEnergy(), and ~CovarianceBinUtility().

std::map<double, int> cmf::CovarianceBinUtility::fNCFDHighEdges

private

high side of the nc bins in FD (2020ana)

Definition at line 110 of file CovarianceBinUtility.h.

Referenced by Initialize(), NCHighEdges(), NumSelectionBins(), SetBinning(), TotalBins(), and ~CovarianceBinUtility().

double cmf::CovarianceBinUtility::fNCFDLowEdge

private

low side of the first nc bin in FD

Definition at line 109 of file CovarianceBinUtility.h.

Referenced by Initialize(), NCLowEdge(), and SetBinning().

std::map<double, int> cmf::CovarianceBinUtility::fNCHighEdges

private

high side of the nc bins

Definition at line 106 of file CovarianceBinUtility.h.

Referenced by ~CovarianceBinUtility().

double cmf::CovarianceBinUtility::fNCLowEdge

private

low side of the first nc bin

Definition at line 105 of file CovarianceBinUtility.h.

std::map<int, double> cmf::CovarianceBinUtility::fNCNDBinToEnergy

private

map energy bin to central value (2020ana)

Definition at line 115 of file CovarianceBinUtility.h.

Referenced by Initialize(), SelectionBinToEnergy(), and ~CovarianceBinUtility().

std::map<double, int> cmf::CovarianceBinUtility::fNCNDHighEdges

private

high side of the nc bins in ND (2020ana)

Definition at line 108 of file CovarianceBinUtility.h.

Referenced by Initialize(), NCHighEdges(), NumSelectionBins(), SetBinning(), TotalBins(), and ~CovarianceBinUtility().

double cmf::CovarianceBinUtility::fNCNDLowEdge

private

low side of the first nc bin in ND

Definition at line 107 of file CovarianceBinUtility.h.

Referenced by Initialize(), NCLowEdge(), and SetBinning().

std::map<int, double> cmf::CovarianceBinUtility::fNuEFDBinToEnergy

private

map energy bin to central value

Definition at line 113 of file CovarianceBinUtility.h.

Referenced by Initialize(), SelectionBinToEnergy(), and ~CovarianceBinUtility().

std::map<double, int> cmf::CovarianceBinUtility::fNuEFDHighEdges

private

high side of the nue bins in FD

Definition at line 102 of file CovarianceBinUtility.h.

Referenced by Initialize(), NuEHighEdges(), NumSelectionBins(), SetBinning(), TotalBins(), and ~CovarianceBinUtility().

double cmf::CovarianceBinUtility::fNuEFDLowEdge

private

low side of the first nune bin in FD

Definition at line 101 of file CovarianceBinUtility.h.

Referenced by Initialize(), NuELowEdge(), and SetBinning().

std::map<int, double> cmf::CovarianceBinUtility::fNuENDBinToEnergy

private

map energy bin to central value

Definition at line 112 of file CovarianceBinUtility.h.

Referenced by Initialize(), SelectionBinToEnergy(), and ~CovarianceBinUtility().

std::map<double, int> cmf::CovarianceBinUtility::fNuENDHighEdges

private

high side of the nue bins in ND

Definition at line 100 of file CovarianceBinUtility.h.

Referenced by Initialize(), NuEHighEdges(), NumSelectionBins(), SetBinning(), TotalBins(), and ~CovarianceBinUtility().

double cmf::CovarianceBinUtility::fNuENDLowEdge

private

low side of the first nue bin in ND

Definition at line 99 of file CovarianceBinUtility.h.

Referenced by Initialize(), NuELowEdge(), and SetBinning().

std::map<int, double> cmf::CovarianceBinUtility::fNuEPeriphBinToEnergy

private

map energy bin to central value

Definition at line 114 of file CovarianceBinUtility.h.

Referenced by Initialize(), SelectionBinToEnergy(), and ~CovarianceBinUtility().

std::map<double, int> cmf::CovarianceBinUtility::fNuEPeriphHighEdges

private

high side of the nue bins in FD

Definition at line 104 of file CovarianceBinUtility.h.

Referenced by Initialize(), NumSelectionBins(), SelectionHighEdges(), SetBinning(), TotalBins(), and ~CovarianceBinUtility().

double cmf::CovarianceBinUtility::fNuEPeriphLowEdge

private

low side of the first nue bin in FD

Definition at line 103 of file CovarianceBinUtility.h.

Referenced by Initialize(), NuELowEdge(), and SetBinning().

std::map<int, double> cmf::CovarianceBinUtility::fNuMuBinToEnergy

private

map energy bin to central value

Definition at line 111 of file CovarianceBinUtility.h.

Referenced by Initialize(), SelectionBinToEnergy(), and ~CovarianceBinUtility().

std::map<double, int> cmf::CovarianceBinUtility::fNuMuHighEdges

private

high side of the numu bins

Definition at line 98 of file CovarianceBinUtility.h.

Referenced by Initialize(), NumSelectionBins(), NuMuHighEdges(), SetBinning(), TotalBins(), and ~CovarianceBinUtility().

double cmf::CovarianceBinUtility::fNuMuLowEdge

private

low side of the first numu bin

Definition at line 97 of file CovarianceBinUtility.h.

Referenced by Initialize(), NuMuLowEdge(), and SetBinning().

std::map<long, int> cmf::CovarianceBinUtility::fOffsetMap

private

map of keys to offset for different samples

Definition at line 117 of file CovarianceBinUtility.h.

Referenced by Initialize(), KeyToOffset(), and ~CovarianceBinUtility().

---

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

• /cvmfs/nova-development.opensciencegrid.org/novasoft/releases/N21-02-25/CovarianceMatrixFit/utilities/CovarianceBinUtility.h
• /cvmfs/nova-development.opensciencegrid.org/novasoft/releases/N21-02-25/CovarianceMatrixFit/utilities/CovarianceBinUtility.cxx