Deprecated/UpgradeASCIITests/G4TankSimulatorASCII/G4TankTrackingAction.cc

Go to the documentation of this file.

#include "G4TankTrackingAction.h"
#include "G4TankTrackInformation.h"

#include <G4Track.hh>
#include <G4TrackingManager.hh>
#include <G4TrackStatus.hh>

#include <cstddef>
#include <iostream>

using namespace G4TankSimulatorASCII;

using std::cout;
using std::endl;


// Number of Cherenkov photons from injected particle
int G4TankTrackingAction::fNumC = 0;

// Number of Cherenkov photons NOT coming directly from
// the injected particle. (eg. delta rays)
//
int G4TankTrackingAction::fNumCDelta = 0;
int G4TankTrackingAction::fNumBounces = 0;
double G4TankTrackingAction::WCD_Track_p=0.;


double G4TankTrackingAction::ASCII_Track=0.;
double G4TankTrackingAction::ASCII_Track_p=0.;

double G4TankTrackingAction::ASCII_Edep=0.;
double G4TankTrackingAction::ASCII_Edep_p=0.;
double G4TankTrackingAction::ASCII_Edep_Michel=0.;

double G4TankTrackingAction::ASCII_Edep_i[27]={ 0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.};

int G4TankTrackingAction::flagDecay=0;
double G4TankTrackingAction::MichelElectronKE=0.,G4TankTrackingAction::MichelElectronTheta=0.,G4TankTrackingAction::MichelElectronZ=0.;
int G4TankTrackingAction::MichelElectronNPE[2]={0,0};
double G4TankTrackingAction::MichelElectronNPE_time[2]={0.,0.};

G4TankTrackingAction::G4TankTrackingAction()
{

}


G4TankTrackingAction::~G4TankTrackingAction()
{
}


// NOTE : no energy deposition is expected from the decay process along the step.
void
G4TankTrackingAction::PreUserTrackingAction(const G4Track* aTrack)
{
  BeforeTracking(aTrack);

  // Check if the track already has track information
  if(aTrack->GetUserInformation()!=0) return;
  if ( aTrack->GetParentID() !=1 ) return; // Only the decay of parent particle ( Particles considered mu+/mu-/e-/e+/gamma, it works for neutrons also )

  // Check if the track has the creator process (not the case for primaries)
  const G4VProcess* creatorProcess = aTrack->GetCreatorProcess();
  if(!creatorProcess) return;

  // Check if the track was generated by the interested process
  if(creatorProcess->GetProcessName()=="Decay" )
    {
      // If it is the case, create a track information and attach it to the track
      G4TankTrackInformation* anInfo = new G4TankTrackInformation(aTrack);
      G4Track* theTrack = (G4Track*)aTrack;
      theTrack->SetUserInformation(anInfo);
    }
}


void
G4TankTrackingAction::PostUserTrackingAction(const G4Track* aTrack)
{
  AfterTracking(aTrack);
  //G4cerr << "Number of Cerenkov photons : " << fNumC;

  // Check if the track already has track information
  G4TankTrackInformation* info
    = (G4TankTrackInformation*)(aTrack->GetUserInformation());
  if(info==0) return;

  // Copy parent's track information to its daughters
  G4TrackVector* secondaries = fpTrackingManager->GimmeSecondaries();
  if(secondaries)
    {
      size_t nSeco = secondaries->size();
      for(size_t i=0;i<nSeco;i++)
        { 
          G4TankTrackInformation* infoNew = new G4TankTrackInformation(info);
          (*secondaries)[i]->SetUserInformation(infoNew);
        }
    }
}


void
G4TankTrackingAction::BeforeTracking(const G4Track* aTrack)
{

  fNumBounces = 0;

  const int Code = aTrack->GetDefinition()->GetPDGEncoding();

  if (Code == 0) { // Optical photon

    const int Parent = aTrack->GetParentID();

    // Count number of Cerenkov photons
    if (Parent != 1)  // not directly produced by the injected particle (eg delta ray)
      ++fNumCDelta;
    else // Produced by injected particle
      ++fNumC;
  }
}


void
G4TankTrackingAction::AfterTracking(const G4Track* /*aTrack*/)
{
}