simcore::g4user::PrimaryGeneratorAction Class Reference

Implementation of Geant4 primary generator action. More...

#include <PrimaryGeneratorAction.h>

Public Member Functions
	PrimaryGeneratorAction (const framework::config::Parameters &parameters)
virtual	~PrimaryGeneratorAction ()=default
	Class destructor.
void	GeneratePrimaries (G4Event *event) override
	Generate primaries for the event.

Private Attributes
bool	time_shift_primaries_ {true}
	Should we time-shift so that the primary vertices arrive (or originate) at t=0ns at z_=0mm?

Detailed Description

Implementation of Geant4 primary generator action.

Definition at line 35 of file PrimaryGeneratorAction.h.

Constructor & Destructor Documentation

PrimaryGeneratorAction()

simcore::g4user::PrimaryGeneratorAction::PrimaryGeneratorAction

(

const framework::config::Parameters &

parameters

)

Definition at line 29 of file PrimaryGeneratorAction.cxx.

    : G4VUserPrimaryGeneratorAction() {
  time_shift_primaries_ = parameters.get<bool>("time_shift_primaries");
 
  auto generators{parameters.get<std::vector<framework::config::Parameters> >(
      "generators", {})};
  if (generators.empty()) {
    EXCEPTION_RAISE("MissingGenerator",
                    "Need to define some generator of primaries.");
  }
 
  for (auto& generator : generators) {
    if (not PrimaryGenerator::Factory::get().make(
            generator.get<std::string>("class_name"),
            generator.get<std::string>("instance_name"), generator)) {
      EXCEPTION_RAISE("UnableToCreate",
                      "Unable to create a PrimaryGenerator of type " +
                          generator.get<std::string>("class_name"));
    }
  }
}

Member Function Documentation

GeneratePrimaries()

void simcore::g4user::PrimaryGeneratorAction::GeneratePrimaries ( G4Event * event )

override

Generate primaries for the event.

This is called by the RunManager before being given to the EventManager to process. This means we must create the UserEventInformation here so that it is accessible for including the weights imported by the primary generators.

See also

G4RunManager::GenerateEvent for where this method is called

G4RunManager::ProcessOneEvent for where G4RunManager::GenerateEvent is called

Set UserInformation for primary vertices if they haven't been set before.

Some features downstream of the primaries require certain user info to function properly. This ensures that it happens.

Makes sure that each particle on each primary vertex has

1. A defined UserPrimaryParticleInformation member
2. The HepEvtStatus for this primary info is non-zero

If we are passed configuration to smear the beam spot, we smear the beam spot around the spot generated by the primary generator.

If we are configured to time-shift the primaries, we shift them so that t=0 coincides with primaries arriving at (or coming from) the target.

Parameters

event

The Geant4 event.

Definition at line 52 of file PrimaryGeneratorAction.cxx.

                                                             {
  /*
   * Create our Event information first so that it
   * can be accessed by everyone from now on.
   */
  // Make sure we aren't overwriting a different information container
  if (event->GetUserInformation()) {
    EXCEPTION_RAISE(
        "Misconfig",
        "There was a UserEventInformation attached before beginning event."
        "\nI don't know how this happend!!");
  }
 
  // Make our information container and give it to geant4
  // G4Event owns the event information and will delete it
  auto event_info = new UserEventInformation;
  event->SetUserInformation(event_info);
 
  PrimaryGenerator::Factory::get().apply([event](const auto& generator) {
    generator->GeneratePrimaryVertex(event);
  });
 
  // All beam spot smearing is handled by individual generators
  int n_pv = event->GetNumberOfPrimaryVertex();
  if (n_pv > 0) {
    // loop over all vertices generated
    for (int i_pv = 0; i_pv < n_pv; ++i_pv) {
      G4PrimaryVertex* primary_vertex = event->GetPrimaryVertex(i_pv);
 
      if (not primary_vertex) {
        EXCEPTION_RAISE(
            "BadGen",
            "One of the primary generators created a NULL primary vertex.");
      }
 
      // Loop over all particle associated with the primary vertex and
      // set the generator status to 1.
      for (int iparticle = 0; iparticle < primary_vertex->GetNumberOfParticle();
           ++iparticle) {
        G4PrimaryParticle* primary = primary_vertex->GetPrimary(iparticle);
 
        if (not primary) {
          EXCEPTION_RAISE(
              "BadGen",
              "One of the primary generators created a NULL primary particle.");
        }
 
        auto primary_info{dynamic_cast<UserPrimaryParticleInformation*>(
            primary->GetUserInformation())};
        if (not primary_info) {
          // no user info defined
          //  ==> make a new one
          primary_info = new UserPrimaryParticleInformation;
          primary->SetUserInformation(primary_info);
        }  // check if primaryinfo is defined
 
        int hep_status = primary_info->getHepEvtStatus();
        if (hep_status <= 0) {
          // undefined hepStatus ==> set to 1
          primary_info->setHepEvtStatus(1);
        }  // check if hepStatus defined
 
      }  // iparticle - loop over primary particles from this vertex
 
      // include the weight of this primary vertex in the event weight
      event_info->incWeight(primary_vertex->GetWeight());
 
      // shift so that t=0 coincides with primaries arriving at (or coming from)
      // the target
      if (time_shift_primaries_) {
        primary_vertex->SetT0(primary_vertex->GetT0() +
                              primary_vertex->GetZ0() / 299.702547);
      }
 
    }  // iPV - loop over primary vertices
  } else {
    EXCEPTION_RAISE(
        "NoPrimaries",
        "No primary vertices were produced by any of the generators.");
  }
}

Member Data Documentation

time_shift_primaries_

bool simcore::g4user::PrimaryGeneratorAction::time_shift_primaries_ {true}

private

Should we time-shift so that the primary vertices arrive (or originate) at t=0ns at z_=0mm?

Note

This should remain true unless the user knows what they are doing!

Definition at line 90 of file PrimaryGeneratorAction.h.

{true};

---

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

• SimCore/include/SimCore/G4User/PrimaryGeneratorAction.h
• SimCore/src/SimCore/G4User/PrimaryGeneratorAction.cxx