#include <G4AdjointTrackingAction.hh>

Inheritance diagram for G4AdjointTrackingAction:

Public Member Functions
void	ClearEndOfAdjointTrackInfoVectors ()

	G4AdjointTrackingAction (G4AdjointSteppingAction *anAction)

G4double	GetCosthAtEndOfLastAdjointTrack (std::size_t i=0)

G4ThreeVector	GetDirectionAtEndOfLastAdjointTrack (std::size_t i=0)

G4double	GetEkinAtEndOfLastAdjointTrack (std::size_t i=0)

G4double	GetEkinNucAtEndOfLastAdjointTrack (std::size_t i=0)

const G4String &	GetFwdParticleNameAtEndOfLastAdjointTrack ()

G4int	GetFwdParticlePDGEncodingAtEndOfLastAdjointTrack (std::size_t i=0)

G4bool	GetIsAdjointTrackingMode ()

G4int	GetLastFwdParticleIndex (std::size_t i=0)

std::size_t	GetNbOfAdointTracksReachingTheExternalSurface ()

G4ThreeVector	GetPositionAtEndOfLastAdjointTrack (std::size_t i=0)

G4double	GetWeightAtEndOfLastAdjointTrack (std::size_t i=0)

virtual void	PostUserTrackingAction (const G4Track *)

virtual void	PreUserTrackingAction (const G4Track *)

void	RegisterAtEndOfAdjointTrack ()

void	SetListOfPrimaryFwdParticles (std::vector< G4ParticleDefinition * > *aListOfParticles)

virtual void	SetTrackingManagerPointer (G4TrackingManager *pValue)

void	SetUserForwardTrackingAction (G4UserTrackingAction *anAction)

virtual	~G4AdjointTrackingAction ()

Protected Attributes
G4TrackingManager *	fpTrackingManager = nullptr

Private Attributes
G4bool	is_adjoint_tracking_mode = false

G4double	last_cos_th = 0.0

std::vector< G4double >	last_cos_th_vec

G4ThreeVector	last_direction

std::vector< G4ThreeVector >	last_direction_vec

G4double	last_ekin = 0.0

G4double	last_ekin_nuc = 0.0

std::vector< G4double >	last_ekin_nuc_vec

std::vector< G4double >	last_ekin_vec

G4int	last_fwd_part_index = 0

std::vector< G4int >	last_fwd_part_index_vec

G4String	last_fwd_part_name

G4int	last_fwd_part_PDGEncoding = 0

std::vector< G4int >	last_fwd_part_PDGEncoding_vec

G4ThreeVector	last_pos

std::vector< G4ThreeVector >	last_pos_vec

G4double	last_weight = 0.0

std::vector< G4double >	last_weight_vec

std::vector< G4ParticleDefinition * > *	pListOfPrimaryFwdParticles = nullptr

G4AdjointSteppingAction *	theAdjointSteppingAction = nullptr

G4UserTrackingAction *	theUserFwdTrackingAction = nullptr

Detailed Description

Definition at line 49 of file G4AdjointTrackingAction.hh.

Constructor & Destructor Documentation

◆ G4AdjointTrackingAction()

G4AdjointTrackingAction::G4AdjointTrackingAction ( G4AdjointSteppingAction * anAction )

Definition at line 40 of file G4AdjointTrackingAction.cc.

  : theAdjointSteppingAction(anAction)
{
}

◆ ~G4AdjointTrackingAction()

G4AdjointTrackingAction::~G4AdjointTrackingAction ( )

virtual

Definition at line 48 of file G4AdjointTrackingAction.cc.

49{

50}

Member Function Documentation

◆ ClearEndOfAdjointTrackInfoVectors()

void G4AdjointTrackingAction::ClearEndOfAdjointTrackInfoVectors ( )

Definition at line 138 of file G4AdjointTrackingAction.cc.

{
  last_pos_vec.clear();
  last_direction_vec.clear();
  last_ekin_vec.clear();
  last_ekin_nuc_vec.clear();
  last_cos_th_vec.clear();
  last_weight_vec.clear();
  last_fwd_part_PDGEncoding_vec.clear();
  last_fwd_part_index_vec.clear();
}

References last_cos_th_vec, last_direction_vec, last_ekin_nuc_vec, last_ekin_vec, last_fwd_part_index_vec, last_fwd_part_PDGEncoding_vec, last_pos_vec, and last_weight_vec.

Referenced by G4AdjointSimManager::ClearEndOfAdjointTrackInfoVectors().

◆ GetCosthAtEndOfLastAdjointTrack()

G4double G4AdjointTrackingAction::GetCosthAtEndOfLastAdjointTrack ( std::size_t i = 0 )

inline

Definition at line 74 of file G4AdjointTrackingAction.hh.

75 { return last_cos_th_vec[i]; }

References last_cos_th_vec.

Referenced by G4AdjointSimManager::GetCosthAtEndOfLastAdjointTrack().

◆ GetDirectionAtEndOfLastAdjointTrack()

G4ThreeVector G4AdjointTrackingAction::GetDirectionAtEndOfLastAdjointTrack ( std::size_t i = 0 )

inline

Definition at line 66 of file G4AdjointTrackingAction.hh.

67 { return last_direction_vec[i]; }

References last_direction_vec.

Referenced by G4AdjointSimManager::GetDirectionAtEndOfLastAdjointTrack().

◆ GetEkinAtEndOfLastAdjointTrack()

G4double G4AdjointTrackingAction::GetEkinAtEndOfLastAdjointTrack ( std::size_t i = 0 )

inline

Definition at line 68 of file G4AdjointTrackingAction.hh.

69 { return last_ekin_vec[i]; }

References last_ekin_vec.

Referenced by G4AdjointSimManager::GetEkinAtEndOfLastAdjointTrack().

◆ GetEkinNucAtEndOfLastAdjointTrack()

G4double G4AdjointTrackingAction::GetEkinNucAtEndOfLastAdjointTrack ( std::size_t i = 0 )

inline

Definition at line 70 of file G4AdjointTrackingAction.hh.

71 { return last_ekin_nuc_vec[i]; }

References last_ekin_nuc_vec.

Referenced by G4AdjointSimManager::GetEkinNucAtEndOfLastAdjointTrack().

◆ GetFwdParticleNameAtEndOfLastAdjointTrack()

const G4String & G4AdjointTrackingAction::GetFwdParticleNameAtEndOfLastAdjointTrack ( )

inline

Definition at line 76 of file G4AdjointTrackingAction.hh.

77 { return last_fwd_part_name; }

G4AdjointTrackingAction::last_fwd_part_name

G4String last_fwd_part_name

Definition: G4AdjointTrackingAction.hh:102

References last_fwd_part_name.

Referenced by G4AdjointSimManager::GetFwdParticleNameAtEndOfLastAdjointTrack().

◆ GetFwdParticlePDGEncodingAtEndOfLastAdjointTrack()

G4int G4AdjointTrackingAction::GetFwdParticlePDGEncodingAtEndOfLastAdjointTrack ( std::size_t i = 0 )

inline

Definition at line 78 of file G4AdjointTrackingAction.hh.

79 { return last_fwd_part_PDGEncoding_vec[i]; }

References last_fwd_part_PDGEncoding_vec.

Referenced by G4AdjointSimManager::GetFwdParticlePDGEncodingAtEndOfLastAdjointTrack().

◆ GetIsAdjointTrackingMode()

G4bool G4AdjointTrackingAction::GetIsAdjointTrackingMode ( )

inline

Definition at line 80 of file G4AdjointTrackingAction.hh.

81 { return is_adjoint_tracking_mode; }

G4AdjointTrackingAction::is_adjoint_tracking_mode

G4bool is_adjoint_tracking_mode

Definition: G4AdjointTrackingAction.hh:93

References is_adjoint_tracking_mode.

Referenced by G4AdjointSimManager::GetAdjointTrackingMode().

◆ GetLastFwdParticleIndex()

G4int G4AdjointTrackingAction::GetLastFwdParticleIndex ( std::size_t i = 0 )

inline

Definition at line 82 of file G4AdjointTrackingAction.hh.

83 { return last_fwd_part_index_vec[i]; }

References last_fwd_part_index_vec.

Referenced by G4AdjointSimManager::GetFwdParticleIndexAtEndOfLastAdjointTrack().

◆ GetNbOfAdointTracksReachingTheExternalSurface()

std::size_t G4AdjointTrackingAction::GetNbOfAdointTracksReachingTheExternalSurface ( )

inline

Definition at line 84 of file G4AdjointTrackingAction.hh.

85 { return last_pos_vec.size(); }

References last_pos_vec.

Referenced by G4AdjointStackingAction::ClassifyNewTrack(), and G4AdjointSimManager::GetNbOfAdointTracksReachingTheExternalSurface().

◆ GetPositionAtEndOfLastAdjointTrack()

G4ThreeVector G4AdjointTrackingAction::GetPositionAtEndOfLastAdjointTrack ( std::size_t i = 0 )

inline

Definition at line 64 of file G4AdjointTrackingAction.hh.

65 { return last_pos_vec[i]; }

References last_pos_vec.

Referenced by G4AdjointSimManager::GetPositionAtEndOfLastAdjointTrack().

◆ GetWeightAtEndOfLastAdjointTrack()

G4double G4AdjointTrackingAction::GetWeightAtEndOfLastAdjointTrack ( std::size_t i = 0 )

inline

Definition at line 72 of file G4AdjointTrackingAction.hh.

73 { return last_weight_vec[i]; }

References last_weight_vec.

Referenced by G4AdjointSimManager::GetWeightAtEndOfLastAdjointTrack().

◆ PostUserTrackingAction()

void G4AdjointTrackingAction::PostUserTrackingAction ( const G4Track * aTrack )

virtual

Reimplemented from G4UserTrackingAction.

Definition at line 75 of file G4AdjointTrackingAction.cc.

{
  // important to have it here !
  //
  last_weight = theAdjointSteppingAction->GetLastWeight();
  last_ekin = theAdjointSteppingAction->GetLastEkin();
 
  if(!is_adjoint_tracking_mode)
  {
    if (theUserFwdTrackingAction != nullptr)
    {
      theUserFwdTrackingAction->PostUserTrackingAction(aTrack);
    }
  }
  else if (theAdjointSteppingAction->GetDidAdjParticleReachTheExtSource())
  {
    last_pos = theAdjointSteppingAction->GetLastPosition();
    last_direction = theAdjointSteppingAction->GetLastMomentum();
    last_direction /=last_direction.mag();
    last_cos_th = last_direction.z();
    G4ParticleDefinition* aPartDef= theAdjointSteppingAction->GetLastPartDef();
    last_fwd_part_name= aPartDef->GetParticleName();
    last_fwd_part_name.erase(0,4);
    last_fwd_part_PDGEncoding=G4ParticleTable::GetParticleTable()
         ->FindParticle(last_fwd_part_name)->GetPDGEncoding();
    last_ekin = theAdjointSteppingAction->GetLastEkin();
    last_ekin_nuc = last_ekin;
    if (aPartDef->GetParticleType() == "adjoint_nucleus")
    {
      G4double  nb_nuc = G4double(aPartDef->GetBaryonNumber());
      last_ekin_nuc /= nb_nuc;
    }
 
    last_fwd_part_index = -1;
    std::size_t i = 0;
    while(i< pListOfPrimaryFwdParticles->size() && last_fwd_part_index<0)
    {
      if ((*pListOfPrimaryFwdParticles)[i]->GetParticleName()
         == last_fwd_part_name)
      {
        last_fwd_part_index=i;
      }
      ++i;
    }
 
    // Fill the vectors
    //
    last_pos_vec.push_back(last_pos);
    last_direction_vec.push_back(last_direction);
    last_ekin_vec.push_back(last_ekin);
    last_ekin_nuc_vec.push_back(last_ekin_nuc);
    last_cos_th_vec.push_back(last_cos_th);
    last_weight_vec.push_back(last_weight);
    last_fwd_part_PDGEncoding_vec.push_back(last_fwd_part_PDGEncoding);
    last_fwd_part_index_vec.push_back(last_fwd_part_index);
  }
  else
  {
  }
}

References G4ParticleTable::FindParticle(), G4ParticleDefinition::GetBaryonNumber(), G4AdjointSteppingAction::GetDidAdjParticleReachTheExtSource(), G4AdjointSteppingAction::GetLastEkin(), G4AdjointSteppingAction::GetLastMomentum(), G4AdjointSteppingAction::GetLastPartDef(), G4AdjointSteppingAction::GetLastPosition(), G4AdjointSteppingAction::GetLastWeight(), G4ParticleDefinition::GetParticleName(), G4ParticleTable::GetParticleTable(), G4ParticleDefinition::GetParticleType(), G4ParticleDefinition::GetPDGEncoding(), is_adjoint_tracking_mode, last_cos_th, last_cos_th_vec, last_direction, last_direction_vec, last_ekin, last_ekin_nuc, last_ekin_nuc_vec, last_ekin_vec, last_fwd_part_index, last_fwd_part_index_vec, last_fwd_part_name, last_fwd_part_PDGEncoding, last_fwd_part_PDGEncoding_vec, last_pos, last_pos_vec, last_weight, last_weight_vec, CLHEP::Hep3Vector::mag(), pListOfPrimaryFwdParticles, G4UserTrackingAction::PostUserTrackingAction(), theAdjointSteppingAction, theUserFwdTrackingAction, and CLHEP::Hep3Vector::z().

◆ PreUserTrackingAction()

void G4AdjointTrackingAction::PreUserTrackingAction ( const G4Track * aTrack )

virtual

Reimplemented from G4UserTrackingAction.

Definition at line 54 of file G4AdjointTrackingAction.cc.

{
  G4String partType = aTrack->GetParticleDefinition()->GetParticleType();
  if (G4StrUtil::contains(partType, "adjoint"))
  {
    is_adjoint_tracking_mode = true;
    theAdjointSteppingAction->SetPrimWeight(aTrack->GetWeight());
  }
  else
  {
    is_adjoint_tracking_mode = false;
    if (theUserFwdTrackingAction)
    {
      theUserFwdTrackingAction->PreUserTrackingAction(aTrack);
    }
  }
  theAdjointSteppingAction->SetAdjointTrackingMode(is_adjoint_tracking_mode);
}