lem::Library Class Reference

Collection of events for matching. More...

#include "/cvmfs/nova-development.opensciencegrid.org/novasoft/releases/N21-05-06/LEM/func/Library.h"

Public Member Functions
	~Library ()
const LibrarySummary &	Summary () const
unsigned int	NEvents () const
const MatchableEvent &	Event (int i) const
const MatchableEvent &	DownsampledEvent (int i) const
const MatchableEvent *	Events () const
const MatchableEvent *	DownsampledEvents () const

Static Public Member Functions
static Library *	FromMMap (const std::string &libPath, bool touchAll)
static Library *	FromMMapOnDemand (const std::string &libPath)
static Library *	FromTrees (const std::string &libPath, int splitLibFactor, int splitLibIndex, int downsampleFactor)

Protected Member Functions
	Library ()
	Library (const Library &)
Library &	operator= (const Library &)

Static Protected Member Functions
static double	SimpleSurvivalProb (const EventSummary &evt)

Protected Attributes
LibrarySummary	fSummary
unsigned int	fLibSize
MatchableEvent *	fLib
MatchableEvent *	fLibDownsample
long	fMappingSize

Detailed Description

Collection of events for matching.

Definition at line 18 of file Library.h.

Constructor & Destructor Documentation

lem::Library::~Library

(

)

Definition at line 187 of file Library.cxx.

References fLib, fLibDownsample, and fMappingSize.

  {
    if(fMappingSize != 0){
      //      munmap(kLibAddr, fMappingSize);
    }
    else{
      // Was loaded from trees
      delete[] fLib;
      delete[] fLibDownsample;
    }
  }

lem::Library::Library ( )

inlineprotected

Definition at line 38 of file Library.h.

References evt, operator=(), and SimpleSurvivalProb().

Referenced by FromMMap(), FromMMapOnDemand(), and FromTrees().

{}

lem::Library::Library ( const Library &  )

protected

Member Function Documentation

const MatchableEvent& lem::Library::DownsampledEvent ( int  i ) const

inline

Definition at line 33 of file Library.h.

References fLibDownsample, and MECModelEnuComparisons::i.

Referenced by lem::FindMatchesTranspose::endRun().

{return fLibDownsample[i];}

const MatchableEvent* lem::Library::DownsampledEvents ( ) const

inline

Definition at line 35 of file Library.h.

References fLibDownsample.

Referenced by lem::LEM::produce().

{return fLibDownsample;}

const MatchableEvent& lem::Library::Event ( int  i ) const

inline

Definition at line 32 of file Library.h.

References fLib, and MECModelEnuComparisons::i.

Referenced by dump_event(), lem::FindMatchesTranspose::endRun(), lem::FindLEMMatches::produce(), lem::LEM::produce(), refresh(), and show_potential().

{return fLib[i];}

const MatchableEvent* lem::Library::Events ( ) const

inline

Definition at line 34 of file Library.h.

References fLib.

Referenced by lem::LEM::produce().

{return fLib;}

Library * lem::Library::FromMMap ( const std::string &  libPath, bool  touchAll  )

static

Definition at line 62 of file Library.cxx.

References fLib, fLibDownsample, fLibSize, fMappingSize, fSummary, cvnie::lib, Library(), lem::MMapFileAtAddress(), and lem::LibrarySummary::N.

Referenced by lem::LEM::beginJob(), dump_event(), lem::FindMatchesTranspose::endRun(), lem::FindMatchesAlg::FindMatchesAlg(), lem::LEMLibrary::LEMLibrary(), preload_lib(), show_event(), and show_potential().

  {
    Library* lib = new Library;

    lib->fMappingSize = MMapFileAtAddress(libPath+"/map.bin", kLibAddr,
                                          touchAll, "/tmp/lem_preload_lib.lock");

    lib->fSummary = *((LibrarySummary*)kLibAddr);

    lib->fLibSize = lib->fSummary.N;

    lib->fLib = (MatchableEvent*)((char*)kLibAddr+sizeof(LibrarySummary));

    lib->fLibDownsample = (MatchableEvent*)((char*)lib->fLib+lib->fLibSize*sizeof(MatchableEvent));

    return lib;
  }

Library * lem::Library::FromMMapOnDemand ( const std::string &  libPath )

static

Definition at line 82 of file Library.cxx.

References lem::CreateOnDemandMapping(), fLib, fLibDownsample, fLibSize, fSummary, cvnie::lib, Library(), and lem::LibrarySummary::N.

Referenced by lem::LEM::beginJob(), lem::FindMatchesTranspose::endRun(), and lem::LEMLibrary::LEMLibrary().

  {
    CreateOnDemandMapping((libPath+"/map.bin").c_str(), kLibAddr);

    Library* lib = new Library;

    lib->fSummary = *((LibrarySummary*)kLibAddr);

    lib->fLibSize = lib->fSummary.N;

    lib->fLib = (MatchableEvent*)((char*)kLibAddr+sizeof(LibrarySummary));

    lib->fLibDownsample = (MatchableEvent*)((char*)lib->fLib+lib->fLibSize*sizeof(MatchableEvent));

    return lib;
  }

Library * lem::Library::FromTrees ( const std::string &  libPath, int  splitLibFactor, int  splitLibIndex, int  downsampleFactor  )

static

Definition at line 100 of file Library.cxx.

References ana::assert(), lem::EventSummary::ccnc, om::cout, lem::MatchableEvent::Downsampled(), allTimeWatchdog::endl, evt, fLib, fLibDownsample, fLibSize, fMappingSize, lem::EventSummary::FromTree(), fSummary, lem::EventSummary::id, lem::EventSummary::InitFromTree(), lem::EventSummary::IsSig(), cvnie::lib, Library(), lem::LibrarySummary::nSwapNC, lem::LibrarySummary::nTrueNC, lem::EventSummary::origPdg, P, lem::EventSummary::pdg, SimpleSurvivalProb(), sum, lem::LibrarySummary::totBkg, lem::LibrarySummary::totEnrich, lem::LibrarySummary::totSig, and lem::MatchableEvent::weight.

Referenced by lem::LEM::beginJob(), and lem::LEMLibrary::LEMLibrary().

  {
    // TODO: this order won't match what heads is expecting

    TChain ch("makelib/tree");
    ch.Add((libPath+"/lem_hist_sig.root").c_str());
    ch.Add((libPath+"/lem_hist_bkg.root").c_str());
    const int Nunenrich = ch.GetEntries();
    ch.Add((libPath+"/lem_hist_enrich.root").c_str());

    EventSummary::InitFromTree(&ch);

    const int startIdx = (splitLibIndex*ch.GetEntries())/splitLibFactor;
    const int endIdx = ((splitLibIndex+1)*ch.GetEntries())/splitLibFactor;

    const int N = endIdx-startIdx;

    assert(N > 0);

    std::cout << "Loading libraries..." << std::endl;
    int nSig = 0, nBkg = 0;
    double totSig = 0, totBkg = 0, totEnrich = 0;
    int nTrueNC = 0, nSwapNC = 0;

    Library* lib = new Library;
    lib->fMappingSize = 0;

    lib->fLibSize = N;
    lib->fLib = new MatchableEvent[N];
    lib->fLibDownsample = new MatchableEvent[N];

    for(int n = startIdx; n < endIdx; ++n){
      if((n-startIdx)%10000 == 0)
        std::cout << (n-startIdx)/10000 << "/" << N/10000 << std::endl;

      const bool isEnrich = n >= Nunenrich;
      EventSummary sum = EventSummary::FromTree(&ch, n, isEnrich);
      sum.id = n;

      // This has to happen here, and not in MakePID where the other
      // reweighting happens, because it feeds into the totals for the
      // different event classes that we write out.
      const double P = isEnrich ? 1 : SimpleSurvivalProb(sum);
      // Screws up indexing, don't do it. FindMatches can cope
      //      if(P == 0) continue; // beam nue

      lem::MatchableEvent evt(sum);
      evt.weight = P;

      if(!isEnrich && sum.ccnc == 1){
        if(sum.pdg == sum.origPdg) ++nTrueNC; else ++nSwapNC;
      }

      lib->fLib[n-startIdx] = evt;

      lib->fLibDownsample[n-startIdx] = evt.Downsampled(downsampleFactor);

      if(isEnrich){
        ++totEnrich;
      }
      else{
        if(evt.IsSig()){
          ++nSig;
          totSig += P;
        }
        else{
          ++nBkg;
          totBkg += P;
        }
      }
    } // end for n

    lib->fSummary.nTrueNC = nTrueNC;
    lib->fSummary.nSwapNC = nSwapNC;
    lib->fSummary.totBkg = totBkg;
    lib->fSummary.totSig = totSig;
    lib->fSummary.totEnrich = totEnrich;

    std::cout << "Loaded " << lib->fSummary << std::endl;

    return lib;
  }

unsigned int lem::Library::NEvents ( ) const

inline

Definition at line 31 of file Library.h.

References fLibSize.

Referenced by lem::FindMatchesTranspose::endRun(), and lem::LEM::produce().

{return fLibSize;}

Library& lem::Library::operator= ( const Library &  )

protected

Referenced by Library().

double lem::Library::SimpleSurvivalProb ( const EventSummary &  evt )

staticprotected

Definition at line 28 of file Library.cxx.

References abs(), lem::EventSummary::ccnc, e, lem::EventSummary::IsSig(), lem::EventSummary::origPdg, lem::EventSummary::pdg, std::sin(), util::sqr(), and lem::EventSummary::trueEVis.

Referenced by FromTrees(), and Library().

  {
    if(evt.IsSig()){
      // Scale up to keep as much signal as possible in sample, weighted
      // correctly internally. The missing prefactor is (ss2th13*ssth23)
      return util::sqr(sin(1.267*2.35e-3*810/evt.trueEVis));
    }
    if(evt.ccnc == 0 && abs(evt.pdg) == 14 && abs(evt.origPdg) == 14){
      // CC mu
      return 1-util::sqr(sin(1.267*2.35e-3*810/evt.trueEVis));
    }
    if(evt.ccnc == 0 && abs(evt.pdg) == 12 && abs(evt.origPdg) == 12){
      // Beam nue
      //      return 1-.1*util::sqr(sin(1.267*2.35e-3*810/evt.trueEVis));
      // Want library matching to ignore nues
      return 0;
    }
    if(evt.ccnc == 1 && evt.pdg == evt.origPdg){
      // NC
      return 1;
    }
    if(evt.ccnc == 1 && evt.pdg != evt.origPdg){
      // Swapped NC. We play clever reweighting games so this is OK
      return 1;
    }
    if(evt.ccnc == 0 && abs(evt.pdg) == 14 && abs(evt.origPdg) == 12){
      // CC mu appearance, take ss2th13=0.1, maximal th23
      return .1*.5*util::sqr(sin(1.267*2.35e-3*810/evt.trueEVis));
    }
    // Don't know what this is
    return 0;
  }

const LibrarySummary& lem::Library::Summary ( ) const

inline

Definition at line 30 of file Library.h.

References fSummary.

Referenced by lem::FindLEMMatches::beginRun(), lem::LEM::beginRun(), and lem::FindMatchesTranspose::endRun().

{return fSummary;}

Member Data Documentation

MatchableEvent* lem::Library::fLib

protected

Definition at line 46 of file Library.h.

Referenced by Event(), Events(), FromMMap(), FromMMapOnDemand(), FromTrees(), and ~Library().

MatchableEvent* lem::Library::fLibDownsample

protected

Definition at line 47 of file Library.h.

Referenced by DownsampledEvent(), DownsampledEvents(), FromMMap(), FromMMapOnDemand(), FromTrees(), and ~Library().

unsigned int lem::Library::fLibSize

protected

Definition at line 45 of file Library.h.

Referenced by FromMMap(), FromMMapOnDemand(), FromTrees(), and NEvents().

long lem::Library::fMappingSize

protected

Definition at line 49 of file Library.h.

Referenced by FromMMap(), FromTrees(), and ~Library().

LibrarySummary lem::Library::fSummary

protected

Definition at line 44 of file Library.h.

Referenced by FromMMap(), FromMMapOnDemand(), FromTrees(), and Summary().

---

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

• /cvmfs/nova-development.opensciencegrid.org/novasoft/releases/N21-05-06/LEM/func/Library.h
• /cvmfs/nova-development.opensciencegrid.org/novasoft/releases/N21-05-06/LEM/func/Library.cxx