genie::GEVGDriver Class Reference

GENIE Event Generation Driver. A minimalist user interface object for generating neutrino interactions. Each such object is configured for a given initial state and it drives all relevant GENIE neutrino interaction physics simulation code for that state. To set-up MC jobs involving a multitude of possible initial states, including arbitrarily complex neutrino flux and detector geometry descriptions, see the GMCJDriver object. More...

#include <GEVGDriver.h>

Collaboration diagram for genie::GEVGDriver:

Public Member Functions
	GEVGDriver ()
	~GEVGDriver ()
void	UseSplines (void)
void	SetEventGeneratorList (string listname)
void	SetUnphysEventMask (const TBits &mask)
void	Configure (int nu_pdgc, int Z, int A)
void	Configure (const InitialState &init_state)
EventRecord *	GenerateEvent (const TLorentzVector &nu4p)
const InteractionList *	Interactions (void) const
const EventGeneratorList *	EventGenerators (void) const
const EventGeneratorI *	FindGenerator (const Interaction *interaction) const
const Spline *	XSecSumSpline (void) const
const Spline *	XSecSpline (const Interaction *interaction) const
void	CreateSplines (int nknots=-1, double emax=-1, bool inLogE=true)
double	XSecSum (const TLorentzVector &nup4)
void	CreateXSecSumSpline (int nk, double Emin, double Emax, bool inlogE=true)
Range1D_t	ValidEnergyRange (void) const
void	Reset (void)
void	Print (ostream &stream) const

Private Member Functions
void	Init (void)
void	CleanUp (void)
void	BuildInitialState (const InitialState &init_state)
void	BuildGeneratorList (void)
void	BuildInteractionGeneratorMap (void)
void	BuildInteractionSelector (void)
void	AssertIsValidInitState (void) const

Private Attributes
InitialState *	fInitState
	initial state information for driver instance
EventRecord *	fCurrentRecord
	ptr to the event record being processed
EventGeneratorList *	fEvGenList
	all Event Generators available at this job
InteractionSelectorI *	fIntSelector
	interaction selector
InteractionGeneratorMap *	fIntGenMap
	interaction -> generator assosiative container
TBits *	fUnphysEventMask
	controls whether unphysical events are returned
bool	fUseSplines
	controls whether xsecs are computed or interpolated
Spline *	fXSecSumSpl
	sum{xsec(all interactions | this init state)}
unsigned int	fNRecLevel
	recursive mode depth counter
string	fEventGenList
	list of event generators loaded by this driver (what used to be the $GEVGL setting)

Friends
ostream &	operator<< (ostream &stream, const GEVGDriver &driver)

Detailed Description

GENIE Event Generation Driver. A minimalist user interface object for generating neutrino interactions. Each such object is configured for a given initial state and it drives all relevant GENIE neutrino interaction physics simulation code for that state. To set-up MC jobs involving a multitude of possible initial states, including arbitrarily complex neutrino flux and detector geometry descriptions, see the GMCJDriver object.

Author

Costas Andreopoulos <c.andreopoulos \at cern.ch> University of Liverpool

Created:\n August 06, 2004

License:\n Copyright (c) 2003-2025, The GENIE Collaboration

For the full text of the license visit http://copyright.genie-mc.org

Definition at line 54 of file GEVGDriver.h.

Constructor & Destructor Documentation

GEVGDriver()

GEVGDriver::GEVGDriver

(

)

Definition at line 57 of file GEVGDriver.cxx.

{
  this->Init();
}

References Init().

Referenced by operator<<.

~GEVGDriver()

GEVGDriver::~GEVGDriver

(

)

Definition at line 62 of file GEVGDriver.cxx.

{
  this->CleanUp();
}

References CleanUp().

Member Function Documentation

AssertIsValidInitState()

void GEVGDriver::AssertIsValidInitState ( void ) const

private

Definition at line 704 of file GEVGDriver.cxx.

{
  assert(fInitState);
  int ppdgc = fInitState->ProbePdg();
  bool isv = pdg::IsLepton(ppdgc) || pdg::IsDarkMatter(ppdgc) || pdg::IsAntiDarkMatter(ppdgc);
  assert(isv);
}

References fInitState, genie::pdg::IsAntiDarkMatter(), genie::pdg::IsDarkMatter(), and genie::pdg::IsLepton().

Referenced by BuildInitialState().

BuildGeneratorList()

void GEVGDriver::BuildGeneratorList ( void )

private

Load event generators. The list of event generators is named by fEventGenList.

Definition at line 176 of file GEVGDriver.cxx.

{
//! Load event generators.
//! The list of event generators is named by fEventGenList.
 
  LOG("GEVGDriver", pINFO)
    << "Building the event generator list (specified list name: "
    << fEventGenList << ")";
 
  EventGeneratorListAssembler evglist_assembler(fEventGenList.c_str());
  fEvGenList = evglist_assembler.AssembleGeneratorList();
}

References genie::EventGeneratorListAssembler::AssembleGeneratorList(), fEventGenList, fEvGenList, LOG, and pINFO.

Referenced by Configure().

BuildInitialState()

void GEVGDriver::BuildInitialState ( const InitialState & init_state )

private

Definition at line 166 of file GEVGDriver.cxx.

{
  LOG("GEVGDriver", pINFO) << "Setting the initial state";
 
  if(fInitState) delete fInitState;
  fInitState = new InitialState(init_state);
 
  this->AssertIsValidInitState();
}

References AssertIsValidInitState(), fInitState, LOG, and pINFO.

Referenced by Configure().

BuildInteractionGeneratorMap()

void GEVGDriver::BuildInteractionGeneratorMap ( void )

private

Map each possible interaction, for the given initial state, to one of the generators loaded up

Definition at line 189 of file GEVGDriver.cxx.

{
//! Map each possible interaction, for the given initial state, to one
//! of the generators loaded up
 
  LOG("GEVGDriver", pINFO)
         << "Building the interaction -> generator associations...";
 
  fIntGenMap = new InteractionGeneratorMap;
  fIntGenMap->UseGeneratorList(fEvGenList);
  fIntGenMap->BuildMap(*fInitState);
 
  string mesgh = "Interaction -> Generator assignments for Initial State: ";
 
  LOG("GEVGDriver", pDEBUG)
    << utils::print::PrintFramedMesg(mesgh + fInitState->AsString(), 0, '-')
    << *fIntGenMap;
}

References fEvGenList, fInitState, fIntGenMap, LOG, pDEBUG, pINFO, and genie::utils::print::PrintFramedMesg().

Referenced by Configure().

BuildInteractionSelector()

void GEVGDriver::BuildInteractionSelector ( void )

private

Definition at line 208 of file GEVGDriver.cxx.

{
  LOG("GEVGDriver", pINFO) << "Building the interaction selector...";
 
  AlgFactory * algf = AlgFactory::Instance();
 
  if(fIntSelector) delete fIntSelector;
  fIntSelector = dynamic_cast<InteractionSelectorI *> (
        algf->AdoptAlgorithm("genie::PhysInteractionSelector","Default"));
}

References genie::AlgFactory::AdoptAlgorithm(), fIntSelector, genie::AlgFactory::Instance(), LOG, and pINFO.

Referenced by Configure().

CleanUp()

void GEVGDriver::CleanUp ( void )

private

Definition at line 121 of file GEVGDriver.cxx.

{
  if (fUnphysEventMask)  delete fUnphysEventMask;
  if (fInitState)        delete fInitState;
  if (fEvGenList)        delete fEvGenList;
  if (fIntSelector)      delete fIntSelector;
  if (fIntGenMap)        delete fIntGenMap;
  if (fXSecSumSpl)       delete fXSecSumSpl;
}

References fEvGenList, fInitState, fIntGenMap, fIntSelector, fUnphysEventMask, and fXSecSumSpl.

Referenced by Reset(), and ~GEVGDriver().

Configure() [1/2]

void GEVGDriver::Configure

(

const InitialState &

init_state

)

Definition at line 145 of file GEVGDriver.cxx.

{
  InitialState init_state(is.TgtPdg(), is.ProbePdg()); // filter any other init state info
 
  ostringstream mesg;
  mesg << "Configuring event generation driver for initial state: `"
       << init_state.AsString()
       << "' using event generator list: `"
       << fEventGenList << "'.";
 
  LOG("GEVGDriver", pNOTICE)
        << utils::print::PrintFramedMesg(mesg.str(), 0, '*');
 
  this -> BuildInitialState            (init_state);
  this -> BuildGeneratorList           ();
  this -> BuildInteractionGeneratorMap ();
  this -> BuildInteractionSelector     ();
 
  LOG("GEVGDriver", pINFO) << "Done configuring. \n";
}

References genie::InitialState::AsString(), BuildGeneratorList(), BuildInitialState(), BuildInteractionGeneratorMap(), BuildInteractionSelector(), fEventGenList, LOG, pINFO, pNOTICE, genie::utils::print::PrintFramedMesg(), genie::InitialState::ProbePdg(), and genie::InitialState::TgtPdg().

Configure() [2/2]

void GEVGDriver::Configure	(	int	nu_pdgc,
		int	Z,
		int	A )

Definition at line 137 of file GEVGDriver.cxx.

{
  Target target(Z, A);
  InitialState init_state(target, nu_pdgc);
 
  this->Configure(init_state);
}

References Configure().

Referenced by Configure(), Configure(), GenerateEventsAtFixedInitState(), GetEventGenDriver(), main(), and genie::GMCJDriver::PopulateEventGenDriverPool().

CreateSplines()

void GEVGDriver::CreateSplines	(	int	nknots = -1,
		double	emax = -1,
		bool	inLogE = true )

Definition at line 577 of file GEVGDriver.cxx.

{
// Creates all the cross section splines that are needed by this driver.
// It will check for pre-loaded splines and it will skip the creation of the
// splines it already finds loaded.
 
  LOG("GEVGDriver", pINFO)
       << "Creating (missing) splines with [UseLogE: "
                                             << ((useLogE) ? "ON]" : "OFF]");
  // Get the list of spline objects
  XSecSplineList * xsl = XSecSplineList::Instance();
  xsl->SetLogE(useLogE);
 
  EventGeneratorList::const_iterator evgliter; // event generator list iter
  InteractionList::iterator          intliter; // interaction list iter
 
  // loop over all EventGenerator objects used in the current job
  for(evgliter = fEvGenList->begin();
                               evgliter != fEvGenList->end(); ++evgliter) {
     // current event generator
     const EventGeneratorI * evgen = *evgliter;
     LOG("GEVGDriver", pINFO)
            << "Querying [" << evgen->Id().Key()
                                            << "] for its InteractionList";
 
     // ask the event generator to produce a list of all interaction it can
     // generate for the input initial state
     const InteractionListGeneratorI * ilstgen = evgen->IntListGenerator();
     InteractionList * ilst = ilstgen->CreateInteractionList(*fInitState);
     if(!ilst) continue;
 
     // total cross section algorithm used by the current EventGenerator
     const XSecAlgorithmI * alg = evgen->CrossSectionAlg();
 
     // get the energy range of the spline from the EventGenerator
     // validity context
     double Emin = evgen->ValidityContext().Emin() ;
     xsl -> SetMinE( Emin ) ;
 
     double valid_Emax = evgen->ValidityContext().Emax();
 
     // if the user set a maximum energy, create the spline up to this
     // energy - otherwise use the upper limit of the validity range of
     // the current generator
     if ( emax > 0 ) {
       if ( emax > valid_Emax ) {
         LOG("GEVGDriver", pWARN)
           << "Refusing to exceed validity range: Emax = " << valid_Emax;
       }
       emax = TMath::Min(emax,valid_Emax); // don't exceed validity range
     } else {
       emax = valid_Emax;
     }
 
     assert( emax > Emin );
 
     xsl -> SetMaxE( emax ) ;
 
     // number of knots: use specified number. If not set, use 15 knots
     // per decade. Don't use less than 30 knots.
     if ( nknots < 0 ) {
       nknots = (int) (15 * TMath::Log10(emax-Emin));
     }
     nknots = TMath::Max(nknots,30);
 
     xsl -> SetNKnots( nknots ) ;
 
     // loop over all interactions that can be generated and ask the
     // appropriate cross section algorithm to compute its cross section
     for(intliter = ilst->begin(); intliter != ilst->end(); ++intliter) {
 
         // current interaction
         Interaction * interaction = *intliter;
         string code = interaction->AsString();
 
         SLOG("GEVGDriver", pINFO) << "Need xsec spline for " << code;
 
         // only create the spline if it does not already exists
         bool spl_exists = xsl->SplineExists(alg, interaction);
         if(!spl_exists) {
             SLOG("GEVGDriver", pDEBUG)
               << "The spline wasn't loaded at initialization. "
               << "I can build it now but it might take a while...";
             xsl->CreateSpline(alg, interaction) ;
         } else {
             SLOG("GEVGDriver", pDEBUG) << "Spline was found";
         }
     } // loop over interaction that can be generated by this generator
     delete ilst;
     ilst = 0;
  } // loop over event generators
 
  LOG("GEVGDriver", pINFO) << *xsl; // print list of splines
 
  fUseSplines = true;
}

References genie::Interaction::AsString(), genie::InteractionListGeneratorI::CreateInteractionList(), genie::XSecSplineList::CreateSpline(), genie::EventGeneratorI::CrossSectionAlg(), genie::GVldContext::Emax(), genie::GVldContext::Emin(), fEvGenList, fInitState, fUseSplines, genie::Algorithm::Id(), genie::XSecSplineList::Instance(), genie::EventGeneratorI::IntListGenerator(), genie::AlgId::Key(), LOG, pDEBUG, pINFO, pWARN, genie::XSecSplineList::SetLogE(), SLOG, genie::XSecSplineList::SplineExists(), and genie::EventGeneratorI::ValidityContext().

Referenced by genie::GMCJDriver::BootstrapXSecSplines(), GetEventGenDriver(), and main().

CreateXSecSumSpline()

void GEVGDriver::CreateXSecSumSpline	(	int	nk,
		double	Emin,
		double	Emax,
		bool	inlogE = true )

Definition at line 440 of file GEVGDriver.cxx.

{
// This method creates a spline with the *total* cross section vs E (or logE)
// for the initial state that this driver was configured with.
// This spline is used, for example, by the GMCJDriver to select a target
// material out of all the materials in a detector geometry (summing the
// cross sections again and again proved to be expensive...)
 
  LOG("GEVGDriver", pINFO)
     << "Creating spline (sum-xsec = f(" << ((inlogE) ? "logE" : "E")
     << ") in E = [" << Emin << ", " << Emax << "] using " << nk << " knots";
 
  if(!fUseSplines) {
     LOG("GEVGDriver", pFATAL) << "You haven't loaded any splines!! ";
  }
  assert(fUseSplines);
  assert(Emin<Emax && Emin>0 && nk>2);
 
  double logEmin=0, logEmax=0, dE=0;
 
  double * E    = new double[nk];
  double * xsec = new double[nk];
 
  if(inlogE) {
    logEmin = TMath::Log(Emin);
    logEmax = TMath::Log(Emax);
    dE = (logEmax-logEmin)/(nk-1);
  } else {
    dE = (Emax-Emin)/(nk-1);
  }
 
  TLorentzVector p4(0,0,0,0);
 
  for(int i=0; i<nk; i++) {
    double e = (inlogE) ? TMath::Exp(logEmin + i*dE) : Emin + i*dE;
    p4.SetPxPyPzE(0.,0.,e,e);
    double xs = this->XSecSum(p4);
 
    E[i]    = e;
    xsec[i] = xs;
  }
  if (fXSecSumSpl) delete fXSecSumSpl;
  fXSecSumSpl = new Spline(nk, E, xsec);
  delete [] E;
  delete [] xsec;
}

References e, fUseSplines, fXSecSumSpl, LOG, pFATAL, pINFO, and XSecSum().

Referenced by genie::GMCJDriver::BootstrapXSecSplineSummation(), CalcTotalXSec(), and GetEventGenDriver().

EventGenerators()

const EventGeneratorList * genie::GEVGDriver::EventGenerators ( void ) const

inline

Definition at line 80 of file GEVGDriver.h.

{ return fEvGenList; }

References fEvGenList.

FindGenerator()

const EventGeneratorI * GEVGDriver::FindGenerator

(

const Interaction *

interaction

)

const

Definition at line 357 of file GEVGDriver.cxx.

{
  if(!interaction) {
    LOG("GEVGDriver", pWARN) << "Null interaction!!";
    return 0;
  }
  if(!fIntGenMap) {
    LOG("GEVGDriver", pWARN)
       << "Interaction->Generator Map has not being built yet!";
    return 0;
  }
 
  const EventGeneratorI * evgen = fIntGenMap->FindGenerator(interaction);
  return evgen;
}

References fIntGenMap, LOG, and pWARN.

Referenced by main().

GenerateEvent()

EventRecord * GEVGDriver::GenerateEvent

(

const TLorentzVector &

nu4p

)

Definition at line 228 of file GEVGDriver.cxx.

{
  //-- Build initial state information from inputs
  LOG("GEVGDriver", pINFO) << "Creating the initial state";
  InitialState init_state(*fInitState);
  init_state.SetProbeP4(nu4p);
 
  //-- Select the interaction to be generated (amongst the entries of the
  //   InteractionList assembled by the EventGenerators) and bootstrap the
  //   event record
  LOG("GEVGDriver", pINFO)
     << "Selecting an Interaction & Bootstraping the EventRecord";
  fCurrentRecord = fIntSelector->SelectInteraction(fIntGenMap, nu4p);
 
  if(!fCurrentRecord) {
     LOG("GEVGDriver", pWARN)
         << "No interaction could be selected for: "
         << init_state.AsString() << " at E = " << nu4p.E() << " GeV";
     return 0;
  }
 
  //-- Get a ptr to the interaction summary
  LOG("GEVGDriver", pDEBUG) << "Getting the selected interaction";
  Interaction * interaction = fCurrentRecord->Summary();
 
  //-- Find the appropriate concrete EventGeneratorI implementation
  //   for generating this event.
  //
  //   The right EventGeneratorI will be selecting by iterating over the
  //   entries of the EventGeneratorList and compare the interaction
  //   against the ValidityContext declared by each EventGeneratorI
  //
  //   (note: use of the 'Chain of Responsibility' Design Pattern)
 
  LOG("GEVGDriver", pINFO) << "Finding an appropriate EventGenerator";
 
  const EventGeneratorI * evgen = fIntGenMap->FindGenerator(interaction);
  assert(evgen);
 
  RunningThreadInfo * rtinfo = RunningThreadInfo::Instance();
  rtinfo->UpdateRunningThread(evgen);
 
  //-- Generate the selected event
  //
  //   The selected EventGeneratorI subclass will start processing the
  //   event record (by sequentially asking each entry in its list of
  //   EventRecordVisitorI subclasses to visit and process the record).
  //   Most of the actual event generation takes place in this step.
  //
  //   (note: use of the 'Visitor' Design Pattern)
 
  string mesg = "Requesting from event generation thread: " +
         evgen->Id().Key() + " to generate the selected interaction";
 
  LOG("GEVGDriver", pNOTICE)
         << utils::print::PrintFramedMesg(mesg,1,'=');
 
  fCurrentRecord->SetUnphysEventMask(*fUnphysEventMask);
  evgen->ProcessEventRecord(fCurrentRecord);
 
  //-- Check the generated event flags. The default behaviour is
  //   to reject an unphysical event and enter in recursive mode
  //   and try to regenerate it. If an unphysical event mask has
  //   been set, error conditions may be ignored so that the
  //   requested classes of unphysical events can be passed-through.
 
  bool unphys = fCurrentRecord->IsUnphysical();
  if(!unphys) {
     LOG("GEVGDriver", pINFO) << "Returning the current event!";
     fNRecLevel = 0;
     return fCurrentRecord; // The client 'adopts' the event record
  } else {
     LOG("GEVGDriver", pWARN) << "An unphysical event was generated...";
     // Check whether the user wants to ignore the err
     bool accept = fCurrentRecord->Accept();
     if(accept) {
       LOG("GEVGDriver", pWARN)
          << "The generated unphysical event is accepted by the user";
       fNRecLevel = 0;
       return fCurrentRecord; // The client 'adopts' the event record
 
     } else {
       LOG("GEVGDriver", pWARN)
          << "The generated unphysical event is rejected";
       delete fCurrentRecord;
       fCurrentRecord = 0;
       fNRecLevel++; // increase the nested level counter
 
       if(fNRecLevel<=kRecursiveModeMaxDepth) {
          LOG("GEVGDriver", pWARN)
            << "Attempting to regenerate the event...";
          return this->GenerateEvent(nu4p);
       } else {
          LOG("GEVGDriver", pERROR)
               << "Could not produce a physical event after "
                      << kRecursiveModeMaxDepth << " attempts!";
          delete fCurrentRecord;
          fCurrentRecord = 0;
          fNRecLevel = 0;
          return 0;
//        exit(1);
       }
     }
  }
}

References genie::InitialState::AsString(), fCurrentRecord, fInitState, fIntGenMap, fIntSelector, fNRecLevel, fUnphysEventMask, GenerateEvent(), genie::Algorithm::Id(), genie::RunningThreadInfo::Instance(), genie::AlgId::Key(), genie::controls::kRecursiveModeMaxDepth, LOG, pDEBUG, pERROR, pINFO, pNOTICE, genie::utils::print::PrintFramedMesg(), genie::EventRecordVisitorI::ProcessEventRecord(), pWARN, genie::InitialState::SetProbeP4(), and genie::RunningThreadInfo::UpdateRunningThread().

Referenced by GenerateEvent(), GenerateEvent(), genie::GMCJDriver::GenerateEventKinematics(), and GenerateEventsAtFixedInitState().

Init()

void GEVGDriver::Init ( void )

private

Definition at line 67 of file GEVGDriver.cxx.

{
  // initial state for which this driver is configured
  fInitState = 0;
 
  // current event record (ownership is transfered at GenerateEvent())
  fCurrentRecord = 0;
 
  // list of Event Generator objects loaded into the driver
  fEvGenList = 0;
 
  // interaction selector
  fIntSelector = 0;
 
  // flag instructing the driver whether, for each interaction, to compute
  // cross section by running their corresponding XSecAlgorithm or by
  // evaluating their corresponding xsec spline
  fUseSplines = false;
 
  // set of event generators to be loaded by this driver
  fEventGenList = "Default";
 
  // 'depth' counter when entering a recursive mode to re-generate a failed/
  // unphysical event - the driver is not allowed to go into arbitrarily large
  // depths
  fNRecLevel = 0;
 
  // an "interaction" -> "generator" associative contained built for all
  // simulated interactions (from the loaded Event Generators and for the
  // input initial state)
  fIntGenMap = 0;
 
  // A spline describing the sum of all interaction cross sections given an
  // initial state (the init state with which this driver was configured).
  // Create it using the CreateXSecSumSpline() method
  // The sum of all interaction cross sections is used, for example, by
  // GMCJDriver for selecting an initial state.
  fXSecSumSpl = 0;
 
  // Default driver behaviour is to filter out unphysical events
  // If needed, set the fUnphysEventMask bitfield to get pre-selected types of
  // unphysical events (just set to 1 the bit you want ignored from the check).
  // Be warned that the event record for unphysical events might be incomplete
  // depending on the processing step that event generation was stopped.
  fUnphysEventMask = new TBits(GHepFlags::NFlags());
  //fUnphysEventMask->ResetAllBits(true);
  for(unsigned int i = 0; i < GHepFlags::NFlags(); i++) {
   fUnphysEventMask->SetBitNumber(i, true);
  }
  LOG("GEVGDriver", pNOTICE)
    << "Initializing unphysical event mask (bits: " << GHepFlags::NFlags()-1
    << " -> 0) : " << *fUnphysEventMask;
}

References fCurrentRecord, fEventGenList, fEvGenList, fInitState, fIntGenMap, fIntSelector, fNRecLevel, fUnphysEventMask, fUseSplines, fXSecSumSpl, LOG, genie::GHepFlags::NFlags(), and pNOTICE.

Referenced by GEVGDriver(), and Reset().

Interactions()

const InteractionList * GEVGDriver::Interactions

(

void

)

const

Definition at line 334 of file GEVGDriver.cxx.

{
// Returns the list of all interactions that can be generated by this driver
 
  if(!fIntGenMap) {
    LOG("GEVGDriver", pWARN)
       << "Interaction->Generator Map has not being built yet!";
    return 0;
  }
 
  const InteractionList & ilst = fIntGenMap->GetInteractionList();
  return &ilst;
}

References fIntGenMap, LOG, and pWARN.

Referenced by main(), SaveGraphsToRootFile(), and SaveToPsFile().

Print()

void GEVGDriver::Print

(

ostream &

stream

)

const

Definition at line 712 of file GEVGDriver.cxx.

{
  stream
    << "\n\n *********************** GEVGDriver ***************************";
 
  int ppdg   = fInitState->ProbePdg();
  int tgtpdg = fInitState->Tgt().Pdg();
 
  stream << "\n  |---o Probe PDG-code ......: " << ppdg;
  stream << "\n  |---o Target PDG-code .....: " << tgtpdg;
 
  stream << "\n  |---o Using cross section splines is turned "
                                << utils::print::BoolAsIOString(fUseSplines);
  stream << "\n  |---o Unphysical event filter mask ("
         << GHepFlags::NFlags() << "->0) = " << *fUnphysEventMask;
 
  stream << "\n *********************************************************\n";
}

References genie::utils::print::BoolAsIOString(), fInitState, fUnphysEventMask, fUseSplines, and genie::GHepFlags::NFlags().

Referenced by operator<<.

Reset()

void GEVGDriver::Reset

(

void

)

Definition at line 131 of file GEVGDriver.cxx.

{
  this->CleanUp();
  this->Init();
}

References CleanUp(), and Init().

SetEventGeneratorList()

void GEVGDriver::SetEventGeneratorList

(

string

listname

)

Definition at line 348 of file GEVGDriver.cxx.

{
  LOG("GEVGDriver", pNOTICE)
       << "Setting event generator list: " << listname;
 
  fEventGenList = listname;
}

References fEventGenList, LOG, and pNOTICE.

Referenced by GenerateEventsAtFixedInitState(), GetEventGenDriver(), main(), and genie::GMCJDriver::PopulateEventGenDriverPool().

SetUnphysEventMask()

void GEVGDriver::SetUnphysEventMask

(

const TBits &

mask

)

Definition at line 219 of file GEVGDriver.cxx.

{
  *fUnphysEventMask = mask;
 
  LOG("GEVGDriver", pNOTICE)
    << "Setting unphysical event mask (bits: " << GHepFlags::NFlags() - 1
    << " -> 0) : " << *fUnphysEventMask;
}

References fUnphysEventMask, LOG, genie::GHepFlags::NFlags(), and pNOTICE.

Referenced by genie::GMCJDriver::GenerateEventKinematics(), and GenerateEventsAtFixedInitState().

UseSplines()

void GEVGDriver::UseSplines

(

void

)

Definition at line 508 of file GEVGDriver.cxx.

{
// Instructs the driver to use cross section splines rather than computing
// cross sections by integrating the differential cross section model which
// can be very time-consuming.
// **Note**
// -- If you called GEVGDriver::CreateSplines() already the driver would
//    a) assume that you want to use them and b) would be assured that it
//    has all the splines it needs, so you do not need to call this method.
// -- If you populated the XSecSplineList in another way without this driver
//    knowing about it, eg from an external XML file, do call this method
//    to let the driver know that you want to use the splines. However, note
//    that the driver would **explicitly check** that you have loaded all the
//    splines it needs. If not, then its fiery personality will take over and
//    it will refuse your request, reverting back to not using splines.
 
  fUseSplines = true;
 
  // Get the list of spline objects
  // Should have been constructed at the job initialization
  XSecSplineList * xsl = XSecSplineList::Instance();
 
  // If the user wants to use splines, make sure that all the splines needed
  // have been computed or loaded
  if(fUseSplines) {
 
    // Get the list of all interactions that can be generated by this driver
    const InteractionList & ilst = fIntGenMap->GetInteractionList();
 
    // Loop over all interactions & check that all splines have been loaded
    InteractionList::const_iterator intliter;
    for(intliter = ilst.begin(); intliter != ilst.end(); ++intliter) {
 
       // current interaction
       Interaction * interaction = *intliter;
 
       // corresponding cross section algorithm
       const XSecAlgorithmI * xsec_alg =
                 fIntGenMap->FindGenerator(interaction)->CrossSectionAlg();
       assert(xsec_alg);
 
       // spline exists in spline list?
       bool spl_exists = xsl->SplineExists(xsec_alg, interaction);
 
       // update the 'use splines' flag
       fUseSplines = fUseSplines && spl_exists;
 
       if(!spl_exists) {
          if(!xsec_alg) {
            LOG("GEVGDriver", pWARN)
              << "Null cross-section algorithm! Can not load cross-section spline.";
            return;
          }
          if(!interaction) {
            LOG("GEVGDriver", pWARN)
              << "Null interaction! Can not load cross-section spline.";
            return;
          }
          LOG("GEVGDriver", pWARN)
             << "*** At least a spline (algorithm: "
             << xsec_alg->Id().Key() << ", interaction: "
             << interaction->AsString() << " doesn't exist. "
             << "Reverting back to not using splines";
          return;
       }
     } // loop over interaction list
  }//use-splines?
}

References genie::Interaction::AsString(), fIntGenMap, fUseSplines, genie::Algorithm::Id(), genie::XSecSplineList::Instance(), genie::AlgId::Key(), LOG, pWARN, and genie::XSecSplineList::SplineExists().

Referenced by Configure(), and genie::GMCJDriver::PopulateEventGenDriverPool().

ValidEnergyRange()

Range1D_t GEVGDriver::ValidEnergyRange

(

void

)

const

Definition at line 674 of file GEVGDriver.cxx.

{
// loops over all loaded event generation threads, queries for the energy
// range at their validity context and builds the valid energy range for
// this driver
 
  Range1D_t E;
  E.min =  9999;
  E.max = -9999;
 
  EventGeneratorList::const_iterator evgliter; // event generator list iter
 
  // loop over all EventGenerator objects used in the current job
  for(evgliter = fEvGenList->begin();
                               evgliter != fEvGenList->end(); ++evgliter) {
     // current event generator
     const EventGeneratorI * evgen = *evgliter;
 
     // Emin, Emax as declared in current generator's validity context
     double Emin = TMath::Max(0.001,evgen->ValidityContext().Emin());
     double Emax = evgen->ValidityContext().Emax();
 
     // combined Emin, Emax
     E.min = TMath::Min(E.min, Emin);
     E.max = TMath::Max(E.max, Emax);
  }
  assert(E.min<E.max && E.min>=0);
  return E;
}

References genie::GVldContext::Emax(), genie::GVldContext::Emin(), fEvGenList, genie::Range1D_t::max, genie::Range1D_t::min, and genie::EventGeneratorI::ValidityContext().

Referenced by genie::GMCJDriver::BootstrapXSecSplineSummation().

XSecSpline()

const Spline * GEVGDriver::XSecSpline

(

const Interaction *

interaction

)

const

Definition at line 488 of file GEVGDriver.cxx.

{
// Returns the cross section spline for the input interaction as was
// computed from the cross section model associated with that interaction.
 
  if (!fUseSplines) return 0;
 
  // Get the list of spline objects
  // Should have been constructed at the job initialization
  XSecSplineList * xssl = XSecSplineList::Instance();
 
  // get corresponding cross section algorithm for the input interaction
  const XSecAlgorithmI * xsec_alg =
               fIntGenMap->FindGenerator(interaction)->CrossSectionAlg();
  assert(xsec_alg);
 
  const Spline * spl = xssl->GetSpline(xsec_alg,interaction);
  return spl;
}

References fIntGenMap, fUseSplines, genie::XSecSplineList::GetSpline(), and genie::XSecSplineList::Instance().

Referenced by SaveGraphsToRootFile(), and SaveToPsFile().

XSecSum()

double GEVGDriver::XSecSum

(

const TLorentzVector &

nup4

)

Definition at line 373 of file GEVGDriver.cxx.

{
// Computes the sum of the cross sections for all the interactions that can
// be simulated for the given initial state and for the input neutrino energy
//
  LOG("GEVGDriver", pDEBUG) << "Computing the cross section sum";
 
  double xsec_sum = 0;
 
  // Get the list of spline objects
  // Should have been constructed at the job initialization
  XSecSplineList * xssl = XSecSplineList::Instance();
 
  // Get the list of all interactions that can be generated by this driver
  const InteractionList & ilst = fIntGenMap->GetInteractionList();
 
  // Loop over all interactions & compute cross sections
  InteractionList::const_iterator intliter;
  for(intliter = ilst.begin(); intliter != ilst.end(); ++intliter) {
 
     // get current interaction
     Interaction * interaction = new Interaction(**intliter);
     interaction->InitStatePtr()->SetProbeP4(nup4);
 
     string code = interaction->AsString();
     SLOG("GEVGDriver", pDEBUG)
             << "Compute cross section for interaction: \n" << code;
 
     // get corresponding cross section algorithm
     const XSecAlgorithmI * xsec_alg =
               fIntGenMap->FindGenerator(interaction)->CrossSectionAlg();
     assert(xsec_alg);
 
     // compute (or evaluate) the cross section
     double xsec = 0;
     bool spline_exists = xssl->SplineExists(xsec_alg, interaction);
     if (spline_exists && fUseSplines) {
        double E = nup4.Energy();
        xsec = xssl->GetSpline(xsec_alg,interaction)->Evaluate(E);
     } else
        xsec = xsec_alg->Integral(interaction);
 
     xsec = TMath::Max(0., xsec);
 
     // sum-up and report
     xsec_sum += xsec;
     LOG("GEVGDriver", pDEBUG)
            << "\nInteraction   = " << code
            << "\nCross Section "
            << (fUseSplines ? "*interpolated*" : "*computed*")
            << " = " << (xsec/units::cm2) << " cm2";
 
     delete interaction;
  } // loop over event generators
 
  PDGLibrary * pdglib = PDGLibrary::Instance();
  LOG("GEVGDriver", pINFO)
    << "SumXSec("
    << pdglib->Find(fInitState->ProbePdg())->GetName() << "+"
    << pdglib->Find(fInitState->Tgt().Pdg())->GetName() << "->X, "
    << "E = " << nup4.Energy() << " GeV)"
    << (fUseSplines ? "*interpolated*" : "*computed*")
    << " = " << (xsec_sum/units::cm2) << " cm2";
 
  return xsec_sum;
}

References genie::Interaction::AsString(), genie::units::cm2, genie::Spline::Evaluate(), genie::PDGLibrary::Find(), fInitState, fIntGenMap, fUseSplines, genie::XSecSplineList::GetSpline(), genie::Interaction::InitStatePtr(), genie::PDGLibrary::Instance(), genie::XSecSplineList::Instance(), genie::XSecAlgorithmI::Integral(), LOG, pDEBUG, pINFO, genie::InitialState::SetProbeP4(), SLOG, and genie::XSecSplineList::SplineExists().

Referenced by CreateXSecSumSpline().

XSecSumSpline()

const Spline * genie::GEVGDriver::XSecSumSpline ( void ) const

inline

Definition at line 87 of file GEVGDriver.h.

{ return fXSecSumSpl; }

References fXSecSumSpl.

Referenced by CalcTotalXSec(), genie::GMCJDriver::ComputeInteractionProbabilities(), genie::GMCJDriver::ComputeProbScales(), and SaveToPsFile().

Friends And Related Symbol Documentation

operator<<

ostream & operator<< ( ostream & stream, const GEVGDriver & driver )

friend

Definition at line 50 of file GEVGDriver.cxx.

 {
   driver.Print(stream);
   return stream;
 }

References GEVGDriver(), and Print().

Member Data Documentation

fCurrentRecord

EventRecord* genie::GEVGDriver::fCurrentRecord

private

ptr to the event record being processed

Definition at line 119 of file GEVGDriver.h.

Referenced by GenerateEvent(), and Init().

fEventGenList

string genie::GEVGDriver::fEventGenList

private

list of event generators loaded by this driver (what used to be the $GEVGL setting)

Definition at line 127 of file GEVGDriver.h.

Referenced by BuildGeneratorList(), Configure(), Init(), and SetEventGeneratorList().

fEvGenList

EventGeneratorList* genie::GEVGDriver::fEvGenList

private

all Event Generators available at this job

Definition at line 120 of file GEVGDriver.h.

Referenced by BuildGeneratorList(), BuildInteractionGeneratorMap(), CleanUp(), CreateSplines(), EventGenerators(), Init(), and ValidEnergyRange().

fInitState

InitialState* genie::GEVGDriver::fInitState

private

initial state information for driver instance

Definition at line 118 of file GEVGDriver.h.

Referenced by AssertIsValidInitState(), BuildInitialState(), BuildInteractionGeneratorMap(), CleanUp(), CreateSplines(), GenerateEvent(), Init(), Print(), and XSecSum().

fIntGenMap

InteractionGeneratorMap* genie::GEVGDriver::fIntGenMap

private

interaction -> generator assosiative container

Definition at line 122 of file GEVGDriver.h.

Referenced by BuildInteractionGeneratorMap(), CleanUp(), FindGenerator(), GenerateEvent(), Init(), Interactions(), UseSplines(), XSecSpline(), and XSecSum().

fIntSelector

InteractionSelectorI* genie::GEVGDriver::fIntSelector

private

interaction selector

Definition at line 121 of file GEVGDriver.h.

Referenced by BuildInteractionSelector(), CleanUp(), GenerateEvent(), and Init().

fNRecLevel

unsigned int genie::GEVGDriver::fNRecLevel

private

recursive mode depth counter

Definition at line 126 of file GEVGDriver.h.

Referenced by GenerateEvent(), and Init().

fUnphysEventMask

TBits* genie::GEVGDriver::fUnphysEventMask

private

controls whether unphysical events are returned

Definition at line 123 of file GEVGDriver.h.

Referenced by CleanUp(), GenerateEvent(), Init(), Print(), and SetUnphysEventMask().

fUseSplines

bool genie::GEVGDriver::fUseSplines

private

controls whether xsecs are computed or interpolated

Definition at line 124 of file GEVGDriver.h.

Referenced by CreateSplines(), CreateXSecSumSpline(), Init(), Print(), UseSplines(), XSecSpline(), and XSecSum().

fXSecSumSpl

Spline* genie::GEVGDriver::fXSecSumSpl

private

sum{xsec(all interactions | this init state)}

Definition at line 125 of file GEVGDriver.h.

Referenced by CleanUp(), CreateXSecSumSpline(), Init(), and XSecSumSpline().

---

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

• GEVGDriver.h
• GEVGDriver.cxx