Inheritance diagram for evg::AddGenieEventsToArt:

Public Types
using	Parameters = art::EDProducer::Table< evg::AddGenieEventsToArtParams >

using	ModuleType = EDProducer

template<typename UserConfig , typename KeysToIgnore = void>
using	Table = Modifier::Table< UserConfig, KeysToIgnore >

Public Member Functions
	AddGenieEventsToArt (const Parameters &p)

	~AddGenieEventsToArt ()

	AddGenieEventsToArt (AddGenieEventsToArt const &)=delete

	AddGenieEventsToArt (AddGenieEventsToArt &&)=delete

AddGenieEventsToArt &	operator= (AddGenieEventsToArt const &)=delete

AddGenieEventsToArt &	operator= (AddGenieEventsToArt &&)=delete

void	produce (art::Event &e) override

void	doBeginJob (SharedResources const &resources)

void	doEndJob ()

void	doRespondToOpenInputFile (FileBlock const &fb)

void	doRespondToCloseInputFile (FileBlock const &fb)

void	doRespondToOpenOutputFiles (FileBlock const &fb)

void	doRespondToCloseOutputFiles (FileBlock const &fb)

bool	doBeginRun (RunPrincipal &rp, ModuleContext const &mc)

bool	doEndRun (RunPrincipal &rp, ModuleContext const &mc)

bool	doBeginSubRun (SubRunPrincipal &srp, ModuleContext const &mc)

bool	doEndSubRun (SubRunPrincipal &srp, ModuleContext const &mc)

bool	doEvent (EventPrincipal &ep, ModuleContext const &mc, std::atomic< std::size_t > &counts_run, std::atomic< std::size_t > &counts_passed, std::atomic< std::size_t > &counts_failed)

void	fillProductDescriptions ()

void	registerProducts (ProductDescriptions &productsToRegister)

ModuleDescription const &	moduleDescription () const

void	setModuleDescription (ModuleDescription const &)

std::array< std::vector< ProductInfo >, NumBranchTypes > const &	getConsumables () const

void	sortConsumables (std::string const &current_process_name)

std::unique_ptr< Worker >	makeWorker (WorkerParams const &wp)

template<typename T , BranchType BT>
ViewToken< T >	consumesView (InputTag const &tag)

template<typename T , BranchType BT>
ViewToken< T >	mayConsumeView (InputTag const &tag)

Protected Member Functions
ConsumesCollector &	consumesCollector ()

template<typename T , BranchType = InEvent>
ProductToken< T >	consumes (InputTag const &)

template<typename Element , BranchType = InEvent>
ViewToken< Element >	consumesView (InputTag const &)

template<typename T , BranchType = InEvent>
void	consumesMany ()

template<typename T , BranchType = InEvent>
ProductToken< T >	mayConsume (InputTag const &)

template<typename Element , BranchType = InEvent>
ViewToken< Element >	mayConsumeView (InputTag const &)

template<typename T , BranchType = InEvent>
void	mayConsumeMany ()

Private Types
enum	EDistrib { kUnknownDist, kFixed, kFlat, kPoisson, kPoissonMinus1, kGaussian, kRootino }

typedef enum evg::AddGenieEventsToArt::EDistrib	RndDist_t

Private Member Functions
void	ParseCountConfig ()

size_t	GetNumToAdd () const

void	ParseTimeConfig ()

void	ParseVtxOffsetConfig ()

Private Attributes
Parameters	fParams

std::vector< std::string >	fFileList

double	fGlobalTimeOffset

evgb::EvtTimeShiftI *	fTimeShifter

double	fXlo

double	fYlo

double	fZlo

double	fXhi

double	fYhi

double	fZhi

bool	fAddMCFlux

bool	fRandomEntries

std::string	fMyModuleType

std::string	fMyModuleLabel

int	fOutputPrintLevel

std::string	fOutputDumpFileName

std::ostream *	fOutputStream

std::string	fDistName

RndDist_t	fRndDist

double	fRndP1

double	fRndP2

TChain *	fGTreeChain

genie::NtpMCEventRecord *	fMCRec

size_t	fNumMCRec

size_t	fLastUsedMCRec

genie::flux::GNuMIFluxPassThroughInfo *	fGNuMIFluxPassThroughInfo

genie::flux::GSimpleNtpEntry *	fGSimpleNtpEntry

genie::flux::GSimpleNtpNuMI *	fGSimpleNtpNuMI

genie::flux::GSimpleNtpAux *	fGSimpleNtpAux

bsim::Dk2Nu *	fDk2Nu

bsim::NuChoice *	fNuChoice

unsigned int const	fSeed

CLHEP::HepRandomEngine &	fEngine

Detailed Description

Definition at line 193 of file AddGenieEventsToArt_module.cc.

Member Typedef Documentation

using art::EDProducer::ModuleType = EDProducer

inherited

Definition at line 17 of file EDProducer.h.

using evg::AddGenieEventsToArt::Parameters = art::EDProducer::Table<evg::AddGenieEventsToArtParams>

Definition at line 197 of file AddGenieEventsToArt_module.cc.

typedef enum evg::AddGenieEventsToArt::EDistrib evg::AddGenieEventsToArt::RndDist_t

private

template<typename UserConfig , typename KeysToIgnore = void>

using art::detail::Producer::Table = Modifier::Table<UserConfig, KeysToIgnore>

inherited

Definition at line 26 of file Producer.h.

Member Enumeration Documentation

enum evg::AddGenieEventsToArt::EDistrib

private

Enumerator
kUnknownDist
kFixed
kFlat
kPoisson
kPoissonMinus1
kGaussian
kRootino

Definition at line 233 of file AddGenieEventsToArt_module.cc.

                         { kUnknownDist,
                           kFixed,
                           kFlat,
                           kPoisson,
                           kPoissonMinus1,
                           kGaussian,
                           kRootino } RndDist_t;

Constructor & Destructor Documentation

evg::AddGenieEventsToArt::AddGenieEventsToArt ( const Parameters & p )

explicit

gtree->GetBranch("flux")->GetClassName() (const char* 0x31f41c0)"genie::flux::GNuMIFluxPassThroughInfo"

gtree->GetBranch("simple")->GetClassName() (const char* 0x2a2a8e0)"genie::flux::GSimpleNtpEntry" gtree->GetBranch("numi")->GetClassName() (const char* 0x2a33d40)"genie::flux::GSimpleNtpNuMI gtree->GetBranch("aux")->GetClassName() (const char* 0x2a35320)"genie::flux::GSimpleNtpAux"

gtree->GetBranch("dk2nu")->GetClassName() (const char* 0x3457d59)"bsim::Dk2Nu" gtree->GetBranch("nuchoice")->GetClassName() (const char* 0x3479329)"bsim::NuChoice"

Definition at line 299 of file AddGenieEventsToArt_module.cc.

References art::detail::EngineCreator::createEngine(), fAddMCFlux, fDk2Nu, fEngine, fFileList, fGlobalTimeOffset, fGNuMIFluxPassThroughInfo, fGSimpleNtpAux, fGSimpleNtpEntry, fGSimpleNtpNuMI, fGTreeChain, fLastUsedMCRec, fMCRec, fMyModuleLabel, fMyModuleType, fNuChoice, fNumMCRec, fOutputDumpFileName, fOutputPrintLevel, fOutputStream, fParams, fRandomEntries, fSeed, art::ProducerTable< UserConfig, ImplicitConfig, UserKeysToIgnore >::get_PSet(), evgb::GetRandomNumberSeed(), ParseCountConfig(), ParseTimeConfig(), and ParseVtxOffsetConfig().

   : EDProducer(params)
   , fParams(params)
   , fGlobalTimeOffset(0)
   , fTimeShifter(0)
   , fXlo(0), fYlo(0), fZlo(0)
   , fXhi(0), fYhi(0), fZhi(0)
   , fAddMCFlux(false)
   , fRandomEntries(false)
   , fOutputPrintLevel(-1)
   , fOutputStream(0)
     //
   , fDistName("")
   , fRndDist(kUnknownDist)
   , fRndP1(-1)
   , fRndP2(-1)
   , fGTreeChain(new TChain("gtree"))
   , fMCRec(new genie::NtpMCEventRecord)
   , fNumMCRec(0)
   , fLastUsedMCRec(0)
   , fGNuMIFluxPassThroughInfo(0)
   , fGSimpleNtpEntry(0)
   , fGSimpleNtpNuMI(0)
   , fGSimpleNtpAux(0)
   , fDk2Nu(0)
   , fNuChoice(0)
   // get the random number seed, use a random default if not specified
   // in the configuration file.
   , fSeed{fParams().seed() == 0 ? evgb::GetRandomNumberSeed() : fParams().seed()}
   // only need sub-label if using more than one engine for each
   // instance of this module (already tagged by equiv of fMyModuleLabel)
   , fEngine{createEngine(fSeed/*, "HepJamesRandom", sub-label*/)}
 {
 
 #ifdef GENIE_PRE_R3
   // trigger early initialization of PDG database & GENIE message service
   // just to get it out of the way and not intermixed with other output
   genie::PDGLibrary::Instance();
 #else
   // get the GENIE banner out of the way
   // no longer can use genie::PDGLibrary::Instance() to do this
   // because that must happen, in some cases in v3_02_xx, after the tune
   // is determined
   // banner is triggered by first use of GENIE Messenger
   // avoid using GENIE macros (possible conflict with mf macros)
   // LOG("GENIE",pInfo) << "Trigger GENIE banner";
   (*genie::Messenger::Instance())("GENIE") << log4cpp::Priority::INFO
                                            << "Trigger GENIE banner";
 #endif
 
   fMyModuleType  = fParams.get_PSet().get<std::string>("module_type");
   fMyModuleLabel = fParams.get_PSet().get<std::string>("module_label");
 
   mf::LogDebug("AddGenieEventsToArt")
     << " ctor start " << fMyModuleLabel
     << " (" << fMyModuleType << ") " << std::endl << std::flush;
 
   fFileList         = fParams().fileList();
   fGlobalTimeOffset = fParams().globalTimeOffset();
   fAddMCFlux        = fParams().addMCFlux();
   fRandomEntries    = fParams().randomEntries();
 
   ParseCountConfig();
   ParseVtxOffsetConfig();
   ParseTimeConfig();
 
   produces< std::vector<simb::MCTruth> >();
   produces< std::vector<simb::GTruth>  >();
   produces< art::Assns<simb::MCTruth, simb::GTruth> >();
   if ( fAddMCFlux ) {
     produces< std::vector<simb::MCFlux>  >();
     // Associate every truth with the flux it came from
     produces< art::Assns<simb::MCTruth, simb::MCFlux> >();
   }
 
   //produces< sumdata::SpillData >();
   //produces< sumdata::POTSum, art::InSubRun  >();
   //produces< sumdata::RunData, art::InRun    >();
 
   std::string outFileList = "adding file pattern: ";
   for (size_t i=0; i < fFileList.size(); ++i) {
     outFileList += "\n\t";
     outFileList += fFileList[i];
     fGTreeChain->Add(fFileList[i].c_str());
   }
   mf::LogDebug("AddGenieEventsToArt") << outFileList;
 
   fNumMCRec = fGTreeChain->GetEntries();
   // fLastUsedMCRec is size_t so, unsigned we can't set it to -1
   // which we'd like to so that the first pre-increment gives us "0"
   size_t skip       = fParams().numberToSkip();
   fLastUsedMCRec    = (skip==0) ? fNumMCRec : skip - 1;
 
   // attach flux branches ...
   TObjArray* blist = fGTreeChain->GetListOfBranches();
   TIter    next(blist);
   TObject* obj;
   while ( ( obj = next() ) ) {
     std::string bname = obj->GetName();
     //  should be a list of TBranchElement or TBranchObject items
     //  TBranchObject are ancient ... should have been replaced by Elements
     const TBranchElement* belement = dynamic_cast<const TBranchElement*>(obj);
     const TBranchObject*  bobject  = dynamic_cast<const TBranchObject*>(obj);
     if ( ! belement && ! bobject ) {
       std::string reallyIsA = obj->ClassName();
       mf::LogError("AddGenieEventsToArt")
         << "### supposed branch element '" << bname
         << "' wasn't a TBranchElement/TBranchObject but instead a "
         << reallyIsA << std::endl;
       if ( bname == "gmcrec" ) {
         mf::LogError("AddGenieEventsToArt")
           << "### since this is '" << bname
           << "' this is likely to end very badly badly" << std::endl;
       }
       continue;
     }
     std::string bclass = (belement) ? belement->GetClassName()
                                     : bobject->GetClassName();
     if ( bclass == "genie::NtpMCEventRecord" ) {
       fGTreeChain->SetBranchAddress(bname.c_str(),&fMCRec);
     } else if ( bclass == "genie::flux::GNuMIFluxPassThroughInfo" ) {
       fGNuMIFluxPassThroughInfo = new genie::flux::GNuMIFluxPassThroughInfo;
       fGTreeChain->SetBranchAddress(bname.c_str(),&fGNuMIFluxPassThroughInfo);
     } else if ( bclass == "genie::flux::GSimpleNtpEntry" ) {
       fGSimpleNtpEntry = new genie::flux::GSimpleNtpEntry;
       fGTreeChain->SetBranchAddress(bname.c_str(),&fGSimpleNtpEntry);
     } else if ( bclass == "genie::flux::GSimpleNtpNuMI" ) {
       fGSimpleNtpNuMI = new genie::flux::GSimpleNtpNuMI;
       fGTreeChain->SetBranchAddress(bname.c_str(),&fGSimpleNtpNuMI);
     } else if ( bclass == "genie::flux::GSimpleNtpAux" ) {
       fGSimpleNtpAux = new genie::flux::GSimpleNtpAux;
       fGTreeChain->SetBranchAddress(bname.c_str(),&fGSimpleNtpAux);
     } else if ( bclass == "bsim::Dk2Nu" ) {
       fDk2Nu = new bsim::Dk2Nu;
       fGTreeChain->SetBranchAddress(bname.c_str(),&fDk2Nu);
     } else if ( bclass == "bsim::NuChoice" ) {
       fNuChoice = new bsim::NuChoice;
       fGTreeChain->SetBranchAddress(bname.c_str(),&fNuChoice);
     } else {
       mf::LogError("AddGenieEventsToArt")
         << "### branch element '" << bname
         << "' was unhandled '" << bclass << "' class" << std::endl;
     }
   } // while ( next )
 
   mf::LogInfo("AddGenieEventsToArt")
     << fMyModuleLabel
     << " (" << fMyModuleType << ") "
     << "chain has " << fNumMCRec << " entries"
     << std::endl;
   if ( fNumMCRec == 0 ) {
     throw cet::exception("badInput")
       << "input files have zero entries "
       << __FILE__ << ":" << __LINE__;
   }
 
   // setup to write out file, if requested
   if ( fOutputPrintLevel > 0 ) {
     if ( fOutputDumpFileName == ""          ||
          fOutputDumpFileName == "--"        ||
          fOutputDumpFileName == "cout"      ||
          fOutputDumpFileName == "std::cout"     ) {
       // standardize so we don't check all these again
       fOutputDumpFileName = "std::cout";
       fOutputStream       = &(std::cout);
     } else {
       size_t posl = fOutputDumpFileName.find("%l");
       if ( posl != std::string::npos ) {
         fOutputDumpFileName.replace(posl,2,fMyModuleLabel);
       }
       mf::LogDebug("AddGenieEventToArt")
         << "#### AddGenieEventsToArt::ctor open "
         << fOutputDumpFileName
         << std::endl << std::flush;
       fOutputStream =
         new std::ofstream(fOutputDumpFileName.c_str(),
                           std::ios_base::trunc|std::ios_base::out);
     }
   } // if ( fOutputPrintLevel > 0 )
 
 }

evg::AddGenieEventsToArt::~AddGenieEventsToArt ( )

Definition at line 497 of file AddGenieEventsToArt_module.cc.

References fGTreeChain, fOutputDumpFileName, and fOutputStream.

 {
   // release resources
   if ( fGTreeChain ) delete fGTreeChain;
   if ( fOutputStream && ( fOutputDumpFileName != "std::cout" ) ) {
     std::ofstream* ofs = dynamic_cast<std::ofstream*>(fOutputStream);
     if ( ofs ) {
       mf::LogDebug("AddGenieEventToArt")
         << "#### AddGenieEventsToArt::dtor close "
         << fOutputDumpFileName
         << std::endl << std::flush;
       ofs->flush();
       ofs->close();
     }
     delete fOutputStream;
   }
 }

evg::AddGenieEventsToArt::AddGenieEventsToArt ( AddGenieEventsToArt const & )

delete

evg::AddGenieEventsToArt::AddGenieEventsToArt ( AddGenieEventsToArt && )

delete

Member Function Documentation

template<typename T , BranchType BT>

ProductToken< T > art::ModuleBase::consumes ( InputTag const & tag )

protectedinherited

Definition at line 61 of file ModuleBase.h.

References art::ModuleBase::collector_, and art::ConsumesCollector::consumes().

   {
     return collector_.consumes<T, BT>(tag);
   }

ConsumesCollector & art::ModuleBase::consumesCollector ( )

protectedinherited

Definition at line 57 of file ModuleBase.cc.

References art::ModuleBase::collector_.

   {
     return collector_;
   }

template<typename T , BranchType BT>

void art::ModuleBase::consumesMany ( )

protectedinherited

Definition at line 75 of file ModuleBase.h.

References art::ModuleBase::collector_, and art::ConsumesCollector::consumesMany().

   {
     collector_.consumesMany<T, BT>();
   }

template<typename Element , BranchType = InEvent>

ViewToken<Element> art::ModuleBase::consumesView ( InputTag const & )

protectedinherited

template<typename T , BranchType BT>

ViewToken<T> art::ModuleBase::consumesView ( InputTag const & tag )

inherited

Definition at line 68 of file ModuleBase.h.

References art::ModuleBase::collector_, and art::ConsumesCollector::consumesView().

   {
     return collector_.consumesView<T, BT>(tag);
   }

void art::detail::Producer::doBeginJob ( SharedResources const & resources )

inherited

Definition at line 22 of file Producer.cc.

References art::detail::Producer::beginJobWithFrame(), and art::detail::Producer::setupQueues().

   {
     setupQueues(resources);
     ProcessingFrame const frame{ScheduleID{}};
     beginJobWithFrame(frame);
   }

bool art::detail::Producer::doBeginRun	(	RunPrincipal &	rp,
		ModuleContext const &	mc
	)

inherited

Definition at line 65 of file Producer.cc.

References art::detail::Producer::beginRunWithFrame(), art::RangeSet::forRun(), art::RunPrincipal::makeRun(), r, art::RunPrincipal::runID(), and art::ModuleContext::scheduleID().

   {
     auto r = rp.makeRun(mc, RangeSet::forRun(rp.runID()));
     ProcessingFrame const frame{mc.scheduleID()};
     beginRunWithFrame(r, frame);
     r.commitProducts();
     return true;
   }

bool art::detail::Producer::doBeginSubRun	(	SubRunPrincipal &	srp,
		ModuleContext const &	mc
	)

inherited

Definition at line 85 of file Producer.cc.

References art::detail::Producer::beginSubRunWithFrame(), art::RangeSet::forSubRun(), art::SubRunPrincipal::makeSubRun(), art::ModuleContext::scheduleID(), and art::SubRunPrincipal::subRunID().

   {
     auto sr = srp.makeSubRun(mc, RangeSet::forSubRun(srp.subRunID()));
     ProcessingFrame const frame{mc.scheduleID()};
     beginSubRunWithFrame(sr, frame);
     sr.commitProducts();
     return true;
   }

void art::detail::Producer::doEndJob ( )

inherited

Definition at line 30 of file Producer.cc.

References art::detail::Producer::endJobWithFrame().

   {
     ProcessingFrame const frame{ScheduleID{}};
     endJobWithFrame(frame);
   }

bool art::detail::Producer::doEndRun	(	RunPrincipal &	rp,
		ModuleContext const &	mc
	)

inherited

Definition at line 75 of file Producer.cc.

References art::detail::Producer::endRunWithFrame(), art::RunPrincipal::makeRun(), r, art::ModuleContext::scheduleID(), and art::Principal::seenRanges().

   {
     auto r = rp.makeRun(mc, rp.seenRanges());
     ProcessingFrame const frame{mc.scheduleID()};
     endRunWithFrame(r, frame);
     r.commitProducts();
     return true;
   }

bool art::detail::Producer::doEndSubRun	(	SubRunPrincipal &	srp,
		ModuleContext const &	mc
	)

inherited

Definition at line 95 of file Producer.cc.

References art::detail::Producer::endSubRunWithFrame(), art::SubRunPrincipal::makeSubRun(), art::ModuleContext::scheduleID(), and art::Principal::seenRanges().

   {
     auto sr = srp.makeSubRun(mc, srp.seenRanges());
     ProcessingFrame const frame{mc.scheduleID()};
     endSubRunWithFrame(sr, frame);
     sr.commitProducts();
     return true;
   }

bool art::detail::Producer::doEvent	(	EventPrincipal &	ep,
		ModuleContext const &	mc,
		std::atomic< std::size_t > &	counts_run,
		std::atomic< std::size_t > &	counts_passed,
		std::atomic< std::size_t > &	counts_failed
	)

inherited

Definition at line 105 of file Producer.cc.

References art::detail::Producer::checkPutProducts_, e, art::EventPrincipal::makeEvent(), art::detail::Producer::produceWithFrame(), and art::ModuleContext::scheduleID().

   {
     auto e = ep.makeEvent(mc);
     ++counts_run;
     ProcessingFrame const frame{mc.scheduleID()};
     produceWithFrame(e, frame);
     e.commitProducts(checkPutProducts_, &expectedProducts<InEvent>());
     ++counts_passed;
     return true;
   }

void art::detail::Producer::doRespondToCloseInputFile ( FileBlock const & fb )

inherited

Definition at line 44 of file Producer.cc.

References art::detail::Producer::respondToCloseInputFileWithFrame().

   {
     ProcessingFrame const frame{ScheduleID{}};
     respondToCloseInputFileWithFrame(fb, frame);
   }

void art::detail::Producer::doRespondToCloseOutputFiles ( FileBlock const & fb )

inherited

Definition at line 58 of file Producer.cc.

References art::detail::Producer::respondToCloseOutputFilesWithFrame().

   {
     ProcessingFrame const frame{ScheduleID{}};
     respondToCloseOutputFilesWithFrame(fb, frame);
   }

void art::detail::Producer::doRespondToOpenInputFile ( FileBlock const & fb )

inherited

Definition at line 37 of file Producer.cc.

References art::detail::Producer::respondToOpenInputFileWithFrame().

   {
     ProcessingFrame const frame{ScheduleID{}};
     respondToOpenInputFileWithFrame(fb, frame);
   }

void art::detail::Producer::doRespondToOpenOutputFiles ( FileBlock const & fb )

inherited

Definition at line 51 of file Producer.cc.

References art::detail::Producer::respondToOpenOutputFilesWithFrame().

   {
     ProcessingFrame const frame{ScheduleID{}};
     respondToOpenOutputFilesWithFrame(fb, frame);
   }

void art::Modifier::fillProductDescriptions ( )

inherited

Definition at line 10 of file Modifier.cc.

References art::ProductRegistryHelper::fillDescriptions(), and art::ModuleBase::moduleDescription().

   {
     ProductRegistryHelper::fillDescriptions(moduleDescription());
   }

std::array< std::vector< ProductInfo >, NumBranchTypes > const & art::ModuleBase::getConsumables ( ) const

inherited

Definition at line 43 of file ModuleBase.cc.

References art::ModuleBase::collector_, and art::ConsumesCollector::getConsumables().

   {
     return collector_.getConsumables();
   }

size_t evg::AddGenieEventsToArt::GetNumToAdd ( ) const

private

Definition at line 842 of file AddGenieEventsToArt_module.cc.

References fEngine, fParams, fRndDist, fRndP1, fRndP2, kFixed, kFlat, kGaussian, kPoisson, kPoissonMinus1, kRootino, and tmp.

Referenced by produce().

 {
 
   size_t nchosen = 0;
 
   switch ( fRndDist ) {
   case kFixed:
     {
       nchosen = fRndP1; // nothing random about it
     }
     break;
   case kFlat:
     {
       CLHEP::RandFlat flat(fEngine);   // pass by ref, doesn't own
       // couldn't find good documentation on how this worked ... empirically
       //   fireInt(0) & fireInt(1) give 0 always
       //   fireInt(2) gives 0 or 1
       //   fireInt(3) gives 0, 1, 2
       //   fireInt(4) gives 0, 1, 2, 3
       // ...
       // so for p1=5, p2=7 we want 5 + [0:2] i.e 5+0=5, 5+1=6, 5+2=7
       //    and thus range should be 3
       int range = (int)(fRndP2-fRndP1) + 1;
       nchosen = fRndP1 + flat.fireInt(range);
     }
     break;
   case kPoisson:
   case kPoissonMinus1:
     {
       CLHEP::RandPoisson poisson(fEngine);
       nchosen = poisson.fire(fRndP1);
       if ( fRndDist == kPoissonMinus1 ) {
         if ( nchosen > 0 ) --nchosen;
         else {
           mf::LogError("AddGenieEventsToArt")
             << "fRndDist[type=" << (int)fRndDist
             << "] '" << fParams().timeConfig() << "' "
             << " nchosen " << nchosen
             << " can't subtract 1 for kPoissonMinus1"
             << std::endl;
         }
       }
     }
     break;
   case kGaussian:
     {
       CLHEP::RandGauss gauss(fEngine);
       double tmp = gauss.fire(fRndP1,fRndP2);
       if ( tmp > 0 ) nchosen = (size_t)(tmp);
       else {
         nchosen = 0;
           mf::LogError("AddGenieEventsToArt")
             << "fRndDist[type=" << (int)fRndDist
             << "] '" << fParams().timeConfig() << "' "
             << " tmp " << tmp
             << "; can't return < 0 for kGaussian, return 0"
             << std::endl;
       }
     }
     break;
   case kRootino:
     {
       nchosen = fRndP1; // pre-chosen maximum
     }
     break;
   default:
     {
       nchosen = 0;
       mf::LogError("AddGenieEventsToArt")
         << "fRndDist[type=" << (int)fRndDist
         << "] '" << fParams().timeConfig() << "' not handled"
         << std::endl;
     }
   }
 
   /*
   mf::LogInfo("AddGenieEventsToArt")
     << "fRndDist[type=" << (int)fRndDist
     << "] '" << fParams().timeConfig() << "' "
     << " nchosen " << nchosen << std::endl;
   */
 
   return nchosen;
 }

std::unique_ptr< Worker > art::ModuleBase::makeWorker ( WorkerParams const & wp )

inherited

Definition at line 37 of file ModuleBase.cc.

References art::ModuleBase::doMakeWorker(), and art::NumBranchTypes.

   {
     return doMakeWorker(wp);
   }

template<typename T , BranchType BT>

ProductToken< T > art::ModuleBase::mayConsume ( InputTag const & tag )

protectedinherited

Definition at line 82 of file ModuleBase.h.

References art::ModuleBase::collector_, and art::ConsumesCollector::mayConsume().

   {
     return collector_.mayConsume<T, BT>(tag);
   }

template<typename T , BranchType BT>

void art::ModuleBase::mayConsumeMany ( )

protectedinherited

Definition at line 96 of file ModuleBase.h.

References art::ModuleBase::collector_, and art::ConsumesCollector::mayConsumeMany().

   {
     collector_.mayConsumeMany<T, BT>();
   }

template<typename Element , BranchType = InEvent>

ViewToken<Element> art::ModuleBase::mayConsumeView ( InputTag const & )

protectedinherited

template<typename T , BranchType BT>

ViewToken<T> art::ModuleBase::mayConsumeView ( InputTag const & tag )

inherited

Definition at line 89 of file ModuleBase.h.

References art::ModuleBase::collector_, and art::ConsumesCollector::mayConsumeView().

   {
     return collector_.mayConsumeView<T, BT>(tag);
   }

ModuleDescription const & art::ModuleBase::moduleDescription ( ) const

inherited

Definition at line 13 of file ModuleBase.cc.

References art::errors::LogicError.

Referenced by art::OutputModule::doRespondToOpenInputFile(), art::OutputModule::doWriteEvent(), art::Modifier::fillProductDescriptions(), art::OutputModule::makePlugins_(), art::OutputWorker::OutputWorker(), reco::shower::LArPandoraModularShowerCreation::produce(), art::Modifier::registerProducts(), and art::OutputModule::registerProducts().

   {
     if (md_.has_value()) {
       return *md_;
     }
 
     throw Exception{errors::LogicError,
                     "There was an error while calling moduleDescription().\n"}
       << "The moduleDescription() base-class member function cannot be called\n"
          "during module construction.  To determine which module is "
          "responsible\n"
          "for calling it, find the '<module type>:<module "
          "label>@Construction'\n"
          "tag in the message prefix above.  Please contact artists@fnal.gov\n"
          "for guidance.\n";
   }

AddGenieEventsToArt& evg::AddGenieEventsToArt::operator= ( AddGenieEventsToArt const & )

delete

AddGenieEventsToArt& evg::AddGenieEventsToArt::operator= ( AddGenieEventsToArt && )

delete

void evg::AddGenieEventsToArt::ParseCountConfig ( )

private

Definition at line 723 of file AddGenieEventsToArt_module.cc.

References fDistName, fEngine, fMyModuleLabel, fParams, fRandomEntries, fRndDist, fRndP1, fRndP2, kFixed, kFlat, kGaussian, kPoisson, kPoissonMinus1, and kRootino.

Referenced by AddGenieEventsToArt().

 {
   // Parse countConfig to get parameters on how many to add
   // Should be one of these:
   //    "fixed:"     <N>
   //    "flat:       <Nmin> <Nmax>"
   //    "poisson:    <Nmean>"
   //    "poisson-1:  <Nmean>"
   //    "gauss:      <mean> <rms>"
   //    "rootino:    <maxentries>"  # pick an upper limit
 
   std::string str = fParams().countConfig();
 
   // let user use whatever case they like
   std::transform(str.begin(),str.end(),str.begin(),::tolower);
   // trim any leading whitespace
   if( str.find_first_not_of(" \t\n") != 0)
       str.erase( 0, str.find_first_not_of(" \t\n")  );
 
   size_t i = str.find_first_of(" \t\n");
   fDistName = str.substr(0,i);
   str.erase(0,i);
 
   // now 'str' should have 1 or 2 numerical values
   // special case rootino allowed to have no numerical values
 
   int nf = sscanf(str.c_str(),"%lf %lf",&fRndP1,&fRndP2);
 
   if ( nf == 0 && fDistName.find("rootino") != 0 ) {
     mf::LogError("AddGenieEventsToArt")
       << "ParseCountConfig " << str
       << " had " << nf << " args, expected something '"
       << str << "'"<< std::endl;
     throw cet::exception("badDist countConfig")
       << __FILE__ << ":" << __LINE__
       << " badDist '" << str << "'";
   }
 
   if ( fDistName.find("fix") == 0 || fDistName == "n" || fDistName == "n:" ) {
     fRndDist = kFixed;
     if ( nf != 1 ) {
       mf::LogError("AddGenieEventsToArt")
         << "ParseCountConfig " << fDistName
         << " had 2 args, expected 1: '"
         << str << "', ignoring 2nd" << std::endl;
     }
   } else if ( fDistName.find("flat") == 0 ) {
     fRndDist = kFlat;
     if ( nf == 1 ) fRndP2 = fRndP1;
     // make sure they're ordered
     if ( fRndP2 < fRndP1 ) std::swap(fRndP1,fRndP2);
 
   } else if ( fDistName.find("poiss") == 0 ) {
     fRndDist = kPoisson;
     if ( fDistName.find("-1") != std::string::npos ) fRndDist = kPoissonMinus1;
     if ( nf != 1 ) {
       mf::LogError("AddGenieEventsToArt")
         << "ParseCountConfig " << fDistName
         << " had 2 args, expected 1: '"
         << str << "', ignoring 2nd" << std::endl;
     }
   } else if ( fDistName.find("gaus") == 0 ) {
     fRndDist = kGaussian;
     if ( nf != 2 ) {
       mf::LogError("AddGenieEventsToArt")
         << "ParseCountConfig " << fDistName
         << " had " << nf << " args, expected 2: '"
         << str << "'"<< std::endl;
       throw cet::exception("badDist countConfig")
         << __FILE__ << ":" << __LINE__
         << " badDist '" << fDistName << "'";
     }
   } else if ( fDistName.find("rootino") == 0 ) {
     fRndDist = kRootino;
     if ( fRndP1 < 1 ) fRndP1 = 9999;
     if ( fRandomEntries ) {
       mf::LogError("AddGenieEventsToArt")
         << "ParseCountConfig " << fDistName
         << " 'distribution' is incompatible with randomEntries=true"
         << std::endl;
       throw cet::exception("badDist incompatible DistConfig")
         << __FILE__ << ":" << __LINE__
         << " badDist '" << fDistName << "' w/ randomEntries=true";
     }
   } else {
       mf::LogError("AddGenieEventsToArt")
         << "ParseCountConfig unknown fDistName " << fDistName
         << " had' " << nf << " args '"
         << str << "'"<< std::endl;
     throw cet::exception("unknownDist countConfig")
       << __FILE__ << ":" << __LINE__
       << " unknown '" << fDistName << "'";
   }
 
   mf::LogInfo("AddGenieEventsToArt")
     << "ParseCountConfig label='" << fMyModuleLabel << "'"
     << " fDistName='" << fDistName  << "' (int)"
     << (int)fRndDist << "; cfgstr '" << str << "' --> "
     << " p1 " << fRndP1 << " p2 " << fRndP2 << " nf " << nf
     << std::endl;
 
 #if 0
   // test how fireInt() works ...
   static bool first = true;
   if ( first ) {
     CLHEP::RandFlat flatTest(fEngine);   // pass by ref, doesn't own
     first = false;
     for (int rtest = 0; rtest < 5; ++rtest ) {
       std::cout << " ======= testing CLHEP::RandFlat::fireInt("
                 << rtest << ") =======" << std::endl;
       for (int i=0; i<100; ++i)
         std::cout << "  " << flatTest.fireInt(rtest);
       std::cout << std::endl;
     }
   }
 #endif
 
 }

void evg::AddGenieEventsToArt::ParseTimeConfig ( )

private

Definition at line 928 of file AddGenieEventsToArt_module.cc.

References fMyModuleLabel, fParams, fTimeShifter, evgb::EvtTimeShiftFactory::GetEvtTimeShift(), evgb::EvtTimeShiftFactory::Instance(), evgb::EvtTimeShiftFactory::Print(), and evgb::EvtTimeShiftI::PrintConfig().

Referenced by AddGenieEventsToArt().

 {
   std::string timeConfig = fParams().timeConfig();
   // trim any leading whitespace
   if( timeConfig.find_first_not_of(" \t\n") != 0)
       timeConfig.erase( 0, timeConfig.find_first_not_of(" \t\n")  );
 
   size_t i = timeConfig.find_first_of(": \t\n");
   std::string timeName   = timeConfig.substr(0,i);
   timeConfig.erase(0,i);
 
   mf::LogInfo("AddGenieEventsToArt")
     << "ParseTimeConfig label='" << fMyModuleLabel << "'"
     << " name='" << timeName << "'"
     << " cfg='" << timeConfig << "'" << std::endl;
   if ( timeName == "none" ) timeName = "evgb::EvtTimeNone";
   if ( timeName == "flat" ) timeName = "evgb::EvtTimeFlat";
   if ( timeName == "numi" || timeName == "NuMI" ||
        timeName == "fnal" || timeName == "FNAL"    )
     timeName = "evgb::EvtTimeFNALBeam";
   if ( timeName == "Booster" || timeName == "booster" )
     timeName = "evgb::EvtTimeFNALBeam Booster";
 
   evgb::EvtTimeShiftFactory& timeFactory =
     evgb::EvtTimeShiftFactory::Instance();
   fTimeShifter = timeFactory.GetEvtTimeShift(timeName,timeConfig);
 
   if ( fTimeShifter ) {
     fTimeShifter->PrintConfig();
   } else {
     timeFactory.Print();
     throw cet::exception("BAD TimeShifter")
       << __FILE__ << ":" << __LINE__
       << " unknown '" << timeName << "'";
   }
 }

void evg::AddGenieEventsToArt::ParseVtxOffsetConfig ( )

private

Definition at line 966 of file AddGenieEventsToArt_module.cc.

References DEFINE_ART_MODULE, fMyModuleLabel, fParams, fXhi, fXlo, fYhi, fYlo, fZhi, and fZlo.

Referenced by AddGenieEventsToArt().

 {
   fXlo = fParams().vtxOffsets().xlo();
   fYlo = fParams().vtxOffsets().ylo();
   fZlo = fParams().vtxOffsets().zlo();
   fXhi = fParams().vtxOffsets().xhi();
   fYhi = fParams().vtxOffsets().yhi();
   fZhi = fParams().vtxOffsets().zhi();
 
   if ( fXlo != 0 && fYlo != 0 && fZlo != 0 &&
        fXhi != 0 && fYhi != 0 && fZhi != 0    ) {
     mf::LogInfo("AddGenieEventsToArt")
       << "ParseVtxOffsetConfig label='" << fMyModuleLabel << "' \n"
       << " x [" << std::setw(11) << fXlo << " "
       <<           std::setw(11) << fXhi << " ]\n"
       << " y [" << std::setw(11) << fYlo << " "
       <<           std::setw(11) << fYhi << " ]\n"
       << " z [" << std::setw(11) << fZlo << " "
       <<           std::setw(11) << fZhi << " ]";
   }
 
 }

void evg::AddGenieEventsToArt::produce ( art::Event & e )

overridevirtual

Implements art::EDProducer.

Definition at line 515 of file AddGenieEventsToArt_module.cc.

References evgb::util::CreateAssn(), fAddMCFlux, fDistName, fDk2Nu, fEngine, fGlobalTimeOffset, fGNuMIFluxPassThroughInfo, fGSimpleNtpAux, fGSimpleNtpEntry, fGSimpleNtpNuMI, fGTreeChain, evgb::FillGTruth(), evgb::FillMCFlux(), evgb::FillMCTruth(), fLastUsedMCRec, fMCRec, fNuChoice, fNumMCRec, fOutputPrintLevel, fOutputStream, fParams, fRandomEntries, fRndDist, fTimeShifter, fXhi, fXlo, fYhi, fYlo, fZhi, fZlo, GetNumToAdd(), kRootino, n, art::Event::put(), and evgb::EvtTimeShiftI::TimeOffset().

 {
 
   //std::cerr << "AddGenieEventsToArt::produce start" << std::endl << std::flush;
 
   std::unique_ptr< std::vector<simb::MCTruth> >
     mctruthcol(new std::vector<simb::MCTruth>);
 
   std::unique_ptr< std::vector<simb::GTruth> >
      gtruthcol(new std::vector<simb::GTruth>);
 
   std::unique_ptr< std::vector<simb::MCFlux> >
      mcfluxcol(new std::vector<simb::MCFlux>);
 
   std::unique_ptr< art::Assns<simb::MCTruth, simb::MCFlux> >
         tfassn(new art::Assns<simb::MCTruth, simb::MCFlux>);
   std::unique_ptr< art::Assns<simb::MCTruth, simb::GTruth> >
         tgassn(new art::Assns<simb::MCTruth, simb::GTruth>);
 
   // number of interactions to add to _this_ record/"event"
   size_t n = GetNumToAdd();
 
   // make a list of entries in TChain to use for this overlay
   // same entry should never be in the list twice ...
   std::vector<size_t> entries;
 
   CLHEP::RandFlat flat(fEngine);
 
   mf::LogDebug("AddGeniEventsToArt") << "#### AddGenieEventsToArt::produce "
                                     << "attempt to get " << n << " entries "
                                     << "from " << fDistName << " distribution";
   //std::ostringstream msg;
   //msg << "add indx to list: \n";
   while ( entries.size() != n ) {
     if ( ! fRandomEntries ) {
       // going through file sequentially
       ++fLastUsedMCRec;
       // fLsatUseMCRec is size_t so never less than zero
       if ( fLastUsedMCRec >= fNumMCRec ) fLastUsedMCRec = 0;
       entries.push_back(fLastUsedMCRec);
       //msg << " [" << entries.size()-1 << "] = "
       //    << fLastUsedMCRec << '\n';
     } else {
       size_t indx = flat.fireInt(fNumMCRec);
       // ensure it isn't already there ..
       if ( std::find(entries.begin(),entries.end(),indx) != entries.end() ) {
         // mf::LogInfo("AddGeniEventsToArt") << "rejecting "
         //                                   << indx << " as already there";
       } else {
         entries.push_back(indx);
         //msg << " [" << entries.size()-1 << "] = "
         //    << indx << '\n';
       }
     }
   }
   //mf::LogInfo("AddGeniEventsToArt") << "entries.size " << entries.size()
   //                                  << " " << msg.str();
 
   for (size_t i=0; i<n; ++i) {
     simb::MCTruth mctruth;
     simb::GTruth  gtruth;
     simb::MCFlux  mcflux;
 
     size_t ientry = entries[i];
     mf::LogDebug("AddGenieEventsToArt")
       << "Event: " << i << " - Using entry " << ientry << std::endl;
     fMCRec->Clear(); // don't leak previously fetched info
     // fetch a single entry from GENIE input file
     fGTreeChain->GetEntry(ientry);
     genie::EventRecord* grec = fMCRec->event;
 
     if (fRndDist == kRootino) {
       // check for end-condition
       //  exactly 1 particle which is a rootino (pdg=0)
       if ( grec->GetEntries() == 1 ) {
         genie::GHepParticle* p = grec->Particle(0);
         if ( p && p->Pdg() == 0 ) {
           // update where we left off
           fLastUsedMCRec = ientry;
           // we're done with the loop, don't go on
           // and don't add this marker to the art event record
           break;
         }
       }
     }
 
     /*
     mf::LogInfo("AddGenieEventsToArt")
       << "#### AddGenieEventsToArt::produce " << i+1 << " of " << n
       << " using entry " << ientry
       << std::endl;
     */
 
     if ( fOutputStream ) {
       // alas not currently able to get current setting
 #if __GENIE_RELEASE_CODE__ >= GRELCODE(3,0,2)
       int plevel = genie::GHepRecord::GetPrintLevel(); //
 #endif
       genie::GHepRecord::SetPrintLevel(fOutputPrintLevel); //
       //std::cerr << "####  AddGenieEventsToArt::produce() writing to "
       //          << fOutputDumpFileName
       //          << std::endl << std::flush;
       *fOutputStream << *fMCRec;
       fOutputStream->flush();
 #if __GENIE_RELEASE_CODE__ >= GRELCODE(3,0,2)
       genie::GHepRecord::SetPrintLevel(plevel);
 #endif
     }
 
     // generate offset in time
     double evtTimeOffset = fGlobalTimeOffset + fTimeShifter->TimeOffset();
 
     // offset vertex position
     double xoff = flat.fire(fXlo,fXhi);
     double yoff = flat.fire(fYlo,fYhi);
     double zoff = flat.fire(fZlo,fZhi);
 
     TLorentzVector vtxOffset(xoff,yoff,zoff,evtTimeOffset);
 
     // convert to simb:: ART objects using GENIE2ART functions
     evgb::FillMCTruth(grec,vtxOffset,mctruth,
                       fParams().inputGenieVersion(),
                       fParams().inputGenieTune(),
                       fParams().addGenieVtxTime());
     evgb::FillGTruth(grec,gtruth);
 
     if ( fAddMCFlux ) {
       if ( fGNuMIFluxPassThroughInfo ) {
         double dk2gen = -99999.;
         evgb::FillMCFlux(fGNuMIFluxPassThroughInfo,dk2gen,mcflux);
       } else if ( fGSimpleNtpEntry ) {
         // cheat at meta data for now ...
         // should be pulling this from input file ... no sure that
         // it is being copied correctly
         // (TChain isn't made of GSimpleNtpMeta objs ... )
         // we need to know how to interpret the Aux variablees
         static genie::flux::GSimpleNtpMeta* meta = 0;
         if ( ! meta ) {
           meta = new genie::flux::GSimpleNtpMeta;
           // hopefully this is the layout
           // aux ints:     tgen
           // aux doubles:  fgXYWgt nimpwt muparpx muparpy muparpz mupare necm
           meta->auxintname.push_back("tgen");
           meta->auxdblname.push_back("fgXYWgt");
           meta->auxdblname.push_back("nimpwt");
           meta->auxdblname.push_back("muparpx");
           meta->auxdblname.push_back("muparpy");
           meta->auxdblname.push_back("muparpz");
           meta->auxdblname.push_back("mupare");
           meta->auxdblname.push_back("necm");
         }
         evgb::FillMCFlux(fGSimpleNtpEntry,fGSimpleNtpNuMI,
                          fGSimpleNtpAux,meta,mcflux);
       } else if ( fDk2Nu ) {
         evgb::FillMCFlux(fDk2Nu,fNuChoice,mcflux);
       }
     }
 
     /*
     genie::NtpMCRecHeader rec_header = fMCRec->hdr;
     //LOG("gevdump", pNOTICE)
     std::cout
       << " ** Event: " << rec_header.ievent // implicit newline in print
       << *grec;
     // add to our collections
     */
 
     mctruthcol->push_back(mctruth);
     gtruthcol->push_back(gtruth);
     if (fAddMCFlux) {
       // deterimine if there is something to fetch
       mcfluxcol->push_back(mcflux);
     }
 
     // LArSoft #include "lardata/Utilities/AssociationUtil.h"
     // NOVA    #include "Utilities/AssociationUtil.h"
     // these util::CreateAssn are taken from LArSoft's GENIEGen_module
     //   ~/Work/DUNE/code/larsim/larsim/EventGenerator/GENIE/GENIEGen_module.cc
     // this seems to be MARK CreateAssn_07 ??  but that's one-to-many
     //    implicit in this is a indx=UNIT_MAX arg so this is only acting
     //    on the last element of truthcol
 
     evgb::util::CreateAssn(*this, evt, *mctruthcol, *gtruthcol, *tgassn,
                            gtruthcol->size()-1, gtruthcol->size());
 
     if (fAddMCFlux) {
       evgb::util::CreateAssn(*this, evt, *mctruthcol, *mcfluxcol, *tfassn,
                              mcfluxcol->size()-1, mcfluxcol->size());
     }
     //
 
   } // done collecting input
 
   //std::cerr << "AddGenieEventsToArt::produce put into event"
   //          << std::endl << std::flush;
 
   // put the collections in the event
   evt.put(std::move(mctruthcol));
   evt.put(std::move(gtruthcol));
   evt.put(std::move(tgassn));
   if ( fAddMCFlux ) {
     evt.put(std::move(mcfluxcol));
     evt.put(std::move(tfassn));
   }
 
 }

void art::Modifier::registerProducts ( ProductDescriptions & productsToRegister )

inherited

Definition at line 16 of file Modifier.cc.

References art::ModuleBase::moduleDescription(), and art::ProductRegistryHelper::registerProducts().

   {
     ProductRegistryHelper::registerProducts(productsToRegister,
                                             moduleDescription());
   }

void art::ModuleBase::setModuleDescription ( ModuleDescription const & md )

inherited

Definition at line 31 of file ModuleBase.cc.

References art::ModuleBase::md_.

   {
     md_ = md;
   }

void art::ModuleBase::sortConsumables ( std::string const & current_process_name )

inherited

Definition at line 49 of file ModuleBase.cc.

References art::ModuleBase::collector_, and art::ConsumesCollector::sortConsumables().

   {
     // Now that we know we have seen all the consumes declarations,
     // sort the results for fast lookup later.
     collector_.sortConsumables(current_process_name);
   }