evgb::EvtTimeFNALBeam Class Reference

configurable FNAL Beam time distribution More...

#include "EvtTimeFNALBeam.h"

Inheritance diagram for evgb::EvtTimeFNALBeam:

Public Member Functions
	EvtTimeFNALBeam (const std::string &config)
virtual	~EvtTimeFNALBeam ()
virtual void	Config (const std::string &config)
virtual double	TimeOffset ()
virtual double	TimeOffset (std::vector< double > bi)
virtual void	PrintConfig (bool verbose=true)
	provide a means of printing the configuration More...
void	SetTimeBetweenBuckets (double val)
	specific methods for this variant More...
double	GetTimeBetweenBuckets () const
void	SetBucketTimeSigma (double val)
double	GetBucketTimeSigma () const
void	SetNBucketsPerBatch (int ival)
int	GetNBucketsPerBatch () const
void	SetNFilledBucketsPerBatch (int ival)
int	GetNFilledBucketsPerBatch () const
void	SetBatchIntensities (std::vector< double > bi)
void	SetDisallowedBatchMask (std::vector< int > disallow)
void	SetGlobalOffset (double val)
double	GetGlobalOffset () const
TRandom *	GetRandomGenerator () const
bool	IsRandomGeneratorOwned () const
void	SetRandomGenerator (TRandom *gen, bool isOwned)

Protected Attributes
TRandom *	fRndmGen
bool	fIsOwned

Private Member Functions
void	CalculateCPDF (std::vector< double > batchi)

Private Attributes
double	fTimeBetweenBuckets
	time between buckets More...
double	fBucketTimeSigma
	how wide is distribution in bucket More...
int	fNBucketsPerBatch
int	fNFilledBucketsPerBatch
std::vector< double >	fCummulativeBatchPDF
	summed prob for batches More...
std::vector< int >	fDisallowedBatchMask
	disallow individual batches More...
double	fGlobalOffset
	always displaced by this (in ns) More...

Detailed Description

configurable FNAL Beam time distribution

Author

Robert Hatcher <rhatcher fnal.gov> Fermi National Accelerator Laboratory

2015-06-22

Version

Id

EvtTimeFNALBeam.h,v 1.1 2015/06/30 18:01:24 rhatcher Exp

Definition at line 22 of file EvtTimeFNALBeam.h.

Constructor & Destructor Documentation

evgb::EvtTimeFNALBeam::EvtTimeFNALBeam

(

const std::string &

config

)

Definition at line 20 of file EvtTimeFNALBeam.cxx.

References Config(), and SetBatchIntensities().

    : EvtTimeShiftI(config)
    , fTimeBetweenBuckets(1e9/53.103e6)
    , fBucketTimeSigma(0.750)
    , fNBucketsPerBatch(84)        // NOvA-era 81+3, MINOS-era 81+5
    , fNFilledBucketsPerBatch(81)  // 81 for both eras
    , fDisallowedBatchMask(6,0)    // don't disallow any
    , fGlobalOffset(0)
  {
    std::vector<double> bi(6,1.0); // 6 equal batches
    SetBatchIntensities(bi);
    Config(config);
  }

evgb::EvtTimeFNALBeam::~EvtTimeFNALBeam ( )

virtual

Definition at line 34 of file EvtTimeFNALBeam.cxx.

{ ; }

Member Function Documentation

void evgb::EvtTimeFNALBeam::CalculateCPDF ( std::vector< double >  batchi )

private

Definition at line 126 of file EvtTimeFNALBeam.cxx.

References fCummulativeBatchPDF, and fDisallowedBatchMask.

Referenced by GetGlobalOffset(), SetBatchIntensities(), and TimeOffset().

  {
    fCummulativeBatchPDF.clear();
    double sum = 0;
    size_t nbi = bi.size();
    for (size_t i=0; i < nbi; ++i) {
      sum += bi[i];
      fCummulativeBatchPDF.push_back(sum);
    }
    // normalize to unit probability
    for (size_t i=0; i < nbi; ++i) fCummulativeBatchPDF[i] /= sum;
    // make sure the mask vector keeps up (but never make it smaller)
    // allowing all new batches
    if ( nbi > fDisallowedBatchMask.size() )
      fDisallowedBatchMask.resize(nbi,0);

    /*
    for (size_t j=0; j<nbi; ++j) {
      std::cout << " CPDF[" << j << "] " << fCummulativeBatchPDF[j]
                << "  " << ((fDisallowedBatchMask[j])?"dis":"") << "allowed"
                << std::endl;
    }
    */

  }

void evgb::EvtTimeFNALBeam::Config ( const std::string &  config )

virtual

each schema must take a string that configures it it is up to the individual model to parse said string and extract parameters

Implements evgb::EvtTimeShiftI.

Definition at line 36 of file EvtTimeFNALBeam.cxx.

References PrintConfig().

Referenced by EvtTimeFNALBeam().

  {
    // parse config string
    if ( config != "" ) {
      std::cerr
        << "!!!!! EvtTimeFNALBeam - not yet up to parsing Config string "
        << ", ignoring:"
        << std::endl
        << "\"" << config << "\""
        << std::endl
        << "Starting with: "
        << std::endl;
      PrintConfig();
    }
  }

double evgb::EvtTimeFNALBeam::GetBucketTimeSigma ( ) const

inline

Definition at line 51 of file EvtTimeFNALBeam.h.

References fBucketTimeSigma.

{ return fBucketTimeSigma; }

double evgb::EvtTimeFNALBeam::GetGlobalOffset ( ) const

inline

Definition at line 61 of file EvtTimeFNALBeam.h.

References CalculateCPDF(), and fGlobalOffset.

{ return fGlobalOffset; }

int evgb::EvtTimeFNALBeam::GetNBucketsPerBatch ( ) const

inline

Definition at line 53 of file EvtTimeFNALBeam.h.

References fNBucketsPerBatch.

{ return fNBucketsPerBatch; }

int evgb::EvtTimeFNALBeam::GetNFilledBucketsPerBatch ( ) const

inline

Definition at line 55 of file EvtTimeFNALBeam.h.

References fNFilledBucketsPerBatch, SetBatchIntensities(), and SetDisallowedBatchMask().

{ return fNFilledBucketsPerBatch; }

TRandom* evgb::EvtTimeShiftI::GetRandomGenerator ( ) const

inlineinherited

Allow users some control over random # sequences An "owned" object is expected to be deleted by the EvtTimeShift obj

Definition at line 61 of file EvtTimeShiftI.h.

References evgb::EvtTimeShiftI::fRndmGen.

{ return fRndmGen; }

double evgb::EvtTimeFNALBeam::GetTimeBetweenBuckets ( ) const

inline

Definition at line 49 of file EvtTimeFNALBeam.h.

References fTimeBetweenBuckets.

{ return fTimeBetweenBuckets; }

bool evgb::EvtTimeShiftI::IsRandomGeneratorOwned ( ) const

inlineinherited

Definition at line 62 of file EvtTimeShiftI.h.

References evgb::EvtTimeShiftI::fIsOwned, and evgb::EvtTimeShiftI::SetRandomGenerator().

{ return fIsOwned; }

void evgb::EvtTimeFNALBeam::PrintConfig ( bool  verbose = true )

virtual

provide a means of printing the configuration

Implements evgb::EvtTimeShiftI.

Definition at line 88 of file EvtTimeFNALBeam.cxx.

References fBucketTimeSigma, fCummulativeBatchPDF, fDisallowedBatchMask, fGlobalOffset, fNBucketsPerBatch, fNFilledBucketsPerBatch, and fTimeBetweenBuckets.

Referenced by Config().

  {
    std::cout << "EvtTimeFNALBeam config: " << std::endl
              << "  TimeBetweenBuckets:     " << fTimeBetweenBuckets << " ns" << std::endl
              << "  BucketTimeSigma:        " << fBucketTimeSigma << " ns" << std::endl
              << "  NBucketsPerBatch:       " << fNBucketsPerBatch << std::endl
              << "  NFilledBucketsPerBatch: " << fNFilledBucketsPerBatch << std::endl
              << "  Relative Fractions:    ";
    double prev=0;
    for (size_t i=0; i < fCummulativeBatchPDF.size(); ++i) {
      double frac = fCummulativeBatchPDF[i] - prev;
      bool skip = fDisallowedBatchMask[i];
      std::cout << " ";
      if (skip) std::cout << "{{";
      std::cout << frac;
      if (skip) std::cout << "}}";
      prev = fCummulativeBatchPDF[i];
    }
    std::cout << std::endl
              << "  GlobalOffset:           " << fGlobalOffset << " ns" << std::endl;
  }

void evgb::EvtTimeFNALBeam::SetBatchIntensities

(

std::vector< double >

bi

)

Definition at line 111 of file EvtTimeFNALBeam.cxx.

References CalculateCPDF().

Referenced by EvtTimeFNALBeam(), and GetNFilledBucketsPerBatch().

  {
    CalculateCPDF(bi);
  }

void evgb::EvtTimeFNALBeam::SetBucketTimeSigma ( double  val )

inline

Definition at line 50 of file EvtTimeFNALBeam.h.

References fBucketTimeSigma.

{ fBucketTimeSigma=val; }

void evgb::EvtTimeFNALBeam::SetDisallowedBatchMask

(

std::vector< int >

disallow

)

Definition at line 116 of file EvtTimeFNALBeam.cxx.

References fCummulativeBatchPDF, and fDisallowedBatchMask.

Referenced by GetNFilledBucketsPerBatch().

  {
    size_t ndis = disallow.size();
    size_t nbi  = fCummulativeBatchPDF.size();
    fDisallowedBatchMask = disallow;
    // expand it so it's mirrors # of batch intensities
    // but allow all that haven't been set
    if ( nbi > ndis ) fDisallowedBatchMask.resize(nbi,0);
  }

void evgb::EvtTimeFNALBeam::SetGlobalOffset ( double  val )

inline

Definition at line 60 of file EvtTimeFNALBeam.h.

References fGlobalOffset.

{ fGlobalOffset=val; }

void evgb::EvtTimeFNALBeam::SetNBucketsPerBatch ( int  ival )

inline

Definition at line 52 of file EvtTimeFNALBeam.h.

References fNBucketsPerBatch.

{ fNBucketsPerBatch=ival; }

void evgb::EvtTimeFNALBeam::SetNFilledBucketsPerBatch ( int  ival )

inline

Definition at line 54 of file EvtTimeFNALBeam.h.

References fNFilledBucketsPerBatch.

{ fNFilledBucketsPerBatch=ival; }

void evgb::EvtTimeShiftI::SetRandomGenerator ( TRandom *  gen, bool  isOwned  )

inherited

Definition at line 28 of file EvtTimeShiftI.cxx.

References evgb::EvtTimeShiftI::fIsOwned, and evgb::EvtTimeShiftI::fRndmGen.

Referenced by evgb::EvtTimeShiftI::IsRandomGeneratorOwned().

  {
    // deal with what we might already have
    if ( fIsOwned ) { delete fRndmGen; fRndmGen = 0; fIsOwned = false; }

    fRndmGen = gen;
    fIsOwned = isOwned;
  }

void evgb::EvtTimeFNALBeam::SetTimeBetweenBuckets ( double  val )

inline

specific methods for this variant

Definition at line 48 of file EvtTimeFNALBeam.h.

References fTimeBetweenBuckets.

{ fTimeBetweenBuckets=val; }

double evgb::EvtTimeFNALBeam::TimeOffset ( )

virtual

return time within a 'record' in nanoseconds version taking array might be used for relative batch fractions that vary on a record-by-record basis

Implements evgb::EvtTimeShiftI.

Definition at line 52 of file EvtTimeFNALBeam.cxx.

References fBucketTimeSigma, fCummulativeBatchPDF, fDisallowedBatchMask, fGlobalOffset, fNBucketsPerBatch, fNFilledBucketsPerBatch, evgb::EvtTimeShiftI::fRndmGen, and fTimeBetweenBuckets.

Referenced by TimeOffset().

  {
    // calculate in small to large

    // pick a time within a bucket
    double offset = fRndmGen->Gaus(0.0,fBucketTimeSigma);

    // pick a bucket within a batch
    // assume all ~ buckets constant in batch until we have another model
    offset +=  fTimeBetweenBuckets *
               (double)fRndmGen->Integer(fNFilledBucketsPerBatch);

    // pick a bucket
    bool   disallowed = true;
    size_t ibatch = 0;
    size_t nbatch = fCummulativeBatchPDF.size();
    double r = 2;
    while ( disallowed ) {
      r = fRndmGen->Uniform();
      for (ibatch=0; ibatch<nbatch; ++ibatch) {
        if ( r <= fCummulativeBatchPDF[ibatch] ) break;
      }
      disallowed = ( fDisallowedBatchMask[ibatch] != 0 );
    }
    offset += fTimeBetweenBuckets*(double)fNBucketsPerBatch*(double)ibatch;

    // finally the global offset
    return offset + fGlobalOffset;
  }

double evgb::EvtTimeFNALBeam::TimeOffset ( std::vector< double >  bi )

virtual

Implements evgb::EvtTimeShiftI.

Definition at line 82 of file EvtTimeFNALBeam.cxx.

References CalculateCPDF(), and TimeOffset().

  {
    CalculateCPDF(bi);
    return TimeOffset();
  }

Member Data Documentation

double evgb::EvtTimeFNALBeam::fBucketTimeSigma

private

how wide is distribution in bucket

Definition at line 68 of file EvtTimeFNALBeam.h.

Referenced by GetBucketTimeSigma(), PrintConfig(), SetBucketTimeSigma(), and TimeOffset().

std::vector<double> evgb::EvtTimeFNALBeam::fCummulativeBatchPDF

private

summed prob for batches

Definition at line 71 of file EvtTimeFNALBeam.h.

Referenced by CalculateCPDF(), PrintConfig(), SetDisallowedBatchMask(), and TimeOffset().

std::vector<int> evgb::EvtTimeFNALBeam::fDisallowedBatchMask

private

disallow individual batches

Definition at line 72 of file EvtTimeFNALBeam.h.

Referenced by CalculateCPDF(), PrintConfig(), SetDisallowedBatchMask(), and TimeOffset().

double evgb::EvtTimeFNALBeam::fGlobalOffset

private

always displaced by this (in ns)

Definition at line 73 of file EvtTimeFNALBeam.h.

Referenced by GetGlobalOffset(), PrintConfig(), SetGlobalOffset(), and TimeOffset().

bool evgb::EvtTimeShiftI::fIsOwned

protectedinherited

Definition at line 69 of file EvtTimeShiftI.h.

Referenced by evgb::EvtTimeShiftI::IsRandomGeneratorOwned(), evgb::EvtTimeShiftI::SetRandomGenerator(), and evgb::EvtTimeShiftI::~EvtTimeShiftI().

int evgb::EvtTimeFNALBeam::fNBucketsPerBatch

private

Definition at line 69 of file EvtTimeFNALBeam.h.

Referenced by GetNBucketsPerBatch(), PrintConfig(), SetNBucketsPerBatch(), and TimeOffset().

int evgb::EvtTimeFNALBeam::fNFilledBucketsPerBatch

private

Definition at line 70 of file EvtTimeFNALBeam.h.

Referenced by GetNFilledBucketsPerBatch(), PrintConfig(), SetNFilledBucketsPerBatch(), and TimeOffset().

TRandom* evgb::EvtTimeShiftI::fRndmGen

protectedinherited

Definition at line 68 of file EvtTimeShiftI.h.

Referenced by evgb::EvtTimeShiftI::GetRandomGenerator(), evgb::EvtTimeShiftI::SetRandomGenerator(), evgb::EvtTimeFlat::TimeOffset(), TimeOffset(), and evgb::EvtTimeShiftI::~EvtTimeShiftI().

double evgb::EvtTimeFNALBeam::fTimeBetweenBuckets

private

time between buckets

Definition at line 67 of file EvtTimeFNALBeam.h.

Referenced by GetTimeBetweenBuckets(), PrintConfig(), SetTimeBetweenBuckets(), and TimeOffset().

---

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

• nutools/EventGeneratorBase/GENIE/EvtTimeFNALBeam.h
• nutools/EventGeneratorBase/GENIE/EvtTimeFNALBeam.cxx