G4UA::CosmicPerigeeAction Class Reference

NEEDS DOCUMENTATION. More...

#include <CosmicPerigeeAction.h>

Inheritance diagram for G4UA::CosmicPerigeeAction:

Collaboration diagram for G4UA::CosmicPerigeeAction:

Classes
struct	Config

Public Member Functions
	CosmicPerigeeAction (const Config &config)
virtual void	UserSteppingAction (const G4Step *) override
virtual void	EndOfEventAction (const G4Event *) override
virtual void	BeginOfEventAction (const G4Event *) override
virtual void	PreUserTrackingAction (const G4Track *) override

Private Attributes
Config	m_config
	Configuration options. More...
SG::WriteHandle< TrackRecordCollection >	m_trackRecordCollection
	Output track track record. More...
bool	m_hasBeenSaved
	State members. More...
double	m_idZ
double	m_idR

Detailed Description

NEEDS DOCUMENTATION.

Definition at line 21 of file CosmicPerigeeAction.h.

Constructor & Destructor Documentation

CosmicPerigeeAction()

G4UA::CosmicPerigeeAction::CosmicPerigeeAction

(

const Config &

config

)

Definition at line 20 of file CosmicPerigeeAction.cxx.

    : m_config(config)
    , m_trackRecordCollection("CosmicPerigee")
    , m_hasBeenSaved(false)
    , m_idZ(3490.) // ID maximum Z coordiate by default.
    , m_idR(1150.) // ID outer radius by default.
  {
  }

Member Function Documentation

BeginOfEventAction()

void G4UA::CosmicPerigeeAction::BeginOfEventAction ( const G4Event *  )

overridevirtual

Definition at line 30 of file CosmicPerigeeAction.cxx.

  {
    if (!m_trackRecordCollection.isValid()) {
      m_trackRecordCollection = std::make_unique<TrackRecordCollection>(
          m_trackRecordCollection.name());
    }
 
    //FIXME need a nice way of getting the maximum size of the ID envelope in R and Z.
    //EnvelopeGeometryManager *gm=EnvelopeGeometryManager::GetGeometryManager();
    //m_idR = gm->IdetOuterRadius();
    //m_idZ = gm->IdetMaxZ();
  }

EndOfEventAction()

void G4UA::CosmicPerigeeAction::EndOfEventAction ( const G4Event *  )

overridevirtual

Definition at line 44 of file CosmicPerigeeAction.cxx.

  {
  }

PreUserTrackingAction()

void G4UA::CosmicPerigeeAction::PreUserTrackingAction ( const G4Track *  )

overridevirtual

Definition at line 49 of file CosmicPerigeeAction.cxx.

  {
    // reset the field
    m_hasBeenSaved = false;
  }

UserSteppingAction()

void G4UA::CosmicPerigeeAction::UserSteppingAction ( const G4Step *  aStep )

overridevirtual

Definition at line 56 of file CosmicPerigeeAction.cxx.

  {
    // See if this is a new track
    if (aStep->GetPreStepPoint()->GetStepStatus() == fUndefined)
      m_hasBeenSaved = false;
 
    // See if we've already saved it
    if (m_hasBeenSaved) return;
 
    // Only save muons or tracks in the ID
    if (aStep->GetTrack()->GetDefinition() != G4MuonPlus::Definition() &&
        aStep->GetTrack()->GetDefinition() != G4MuonMinus::Definition() &&
        ( aStep->GetPostStepPoint()->GetPosition().rho() > m_idR ||
          aStep->GetPostStepPoint()->GetPosition().z() > m_idZ   ||
          aStep->GetPostStepPoint()->GetPosition().z() < -m_idZ   ) ){
      // Both not a muon and not in the ID
      return;
    }
 
    // Check momentum
    if (aStep->GetTrack()->GetMomentum().mag() < m_config.pMinPrimary) return;
 
    // First order approximation of the answer to "is this the perigee"
    // if it is now moving away from the perigee and has not been saved, save it
    // That is true if mom dot pos > 0
    if ( aStep->GetTrack()->GetMomentum().x() * aStep->GetPostStepPoint()->GetPosition().x() +
         aStep->GetTrack()->GetMomentum().y() * aStep->GetPostStepPoint()->GetPosition().y() < 0) return;
 
    // Save the vertex...
    m_hasBeenSaved = true;
 
    // Decide whether to save the prestep or poststep point
    // - which one is closer to Perigee?
    G4StepPoint* preStep = aStep->GetPreStepPoint();
    G4StepPoint* postStep = aStep->GetPostStepPoint();
    G4StepPoint* theStep = preStep;
    if ( fabs( preStep->GetMomentumDirection().x() * preStep->GetPosition().x() +
               preStep->GetMomentumDirection().y() * preStep->GetPosition().y() ) >
         fabs( postStep->GetMomentumDirection().x() * postStep->GetPosition().x() +
               postStep->GetMomentumDirection().y() * postStep->GetPosition().y() ) )
    {
      // Using the post step point
      theStep = postStep;
    }
 
    // Could the following code be optimized?
    // There seems to be a bit of object copying.
    int pdgcode = aStep->GetTrack()->GetDefinition()->GetPDGEncoding();
    double ener = theStep->GetTotalEnergy();
    G4ThreeVector pos = theStep->GetPosition();
    G4ThreeVector mom = theStep->GetMomentum();
    double time       = theStep->GetGlobalTime();
    G4VPhysicalVolume *preVol = theStep->GetPhysicalVolume();
 
    if (aStep->GetTrack()->GetDefinition() == G4Geantino::Definition() ) pdgcode=999;
    if (aStep->GetTrack()->GetDefinition() == G4ChargedGeantino::Definition() ) pdgcode=998;
 
    // Create the TimedTrackRecord
    TrackHelper trHelp(aStep->GetTrack());
    const int barcode = trHelp.GetBarcode(); // FIXME barcode based
    const int id = trHelp.GetUniqueID();
    const int status = trHelp.GetStatus();
    m_trackRecordCollection->Emplace(
                                     pdgcode,
                                     status,
                                     ener,
                                     mom,
                                     pos,
                                     time,
                                     barcode, // FIXME barcode based
                                     id,
                                     preVol->GetName());
  }

Member Data Documentation

m_config

Config G4UA::CosmicPerigeeAction::m_config

private

Configuration options.

Definition at line 43 of file CosmicPerigeeAction.h.

m_hasBeenSaved

bool G4UA::CosmicPerigeeAction::m_hasBeenSaved

private

State members.

Definition at line 49 of file CosmicPerigeeAction.h.

m_idR

double G4UA::CosmicPerigeeAction::m_idR

private

Definition at line 51 of file CosmicPerigeeAction.h.

m_idZ

double G4UA::CosmicPerigeeAction::m_idZ

private

Definition at line 50 of file CosmicPerigeeAction.h.

m_trackRecordCollection

SG::WriteHandle<TrackRecordCollection> G4UA::CosmicPerigeeAction::m_trackRecordCollection

private

Output track track record.

Definition at line 46 of file CosmicPerigeeAction.h.

---

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

• CosmicPerigeeAction.h
• CosmicPerigeeAction.cxx