#include "tracking.hh"
#include "G4Track.hh"
#include "G4ParticleDefinition.hh"
#include "G4ParticleTypes.hh"
#include "G4SystemOfUnits.hh"
#include "G4AnalysisMAnager.hh"
#include "G4Step.hh"
#include "stepping.hh"



MyTrackingAction::MyTrackingAction()
{
    
}

MyTrackingAction::~MyTrackingAction()
{
    
}

void MyTrackingAction::PreUserTrackingAction(const G4Track*)
{

}

void MyTrackingAction::PostUserTrackingAction(const G4Track* track, G4double energyDeposit)
{
    // Check if the particle is a proton
    const G4ParticleDefinition* particle = track->GetDefinition();
    if (particle != G4Proton::ProtonDefinition()) return;

    // Retrieve the position of the particle and print it to the console
    G4ThreeVector position = track->GetPosition();
   
   // G4cout << "Proton position: (" << position.x() / cm << ", "
     //   << position.y() / cm << ", " << position.z() / cm << ") cm" <<G4endl;

    

   // G4double energy = GetEnergyDeposit();
    G4cout << position.x() / cm << ", "
        << position.y() / cm << ", " << position.z() / cm <<"," << energyDeposit << G4endl;
  

    
}