//
// ********************************************************************
// * License and Disclaimer                                           *
// *                                                                  *
// * The  Geant4 software  is  copyright of the Copyright Holders  of *
// * the Geant4 Collaboration.  It is provided  under  the terms  and *
// * conditions of the Geant4 Software License,  included in the file *
// * LICENSE and available at  http://cern.ch/geant4/license .  These *
// * include a list of copyright holders.                             *
// *                                                                  *
// * Neither the authors of this software system, nor their employing *
// * institutes,nor the agencies providing financial support for this *
// * work  make  any representation or  warranty, express or implied, *
// * regarding  this  software system or assume any liability for its *
// * use.  Please see the license in the file  LICENSE  and URL above *
// * for the full disclaimer and the limitation of liability.         *
// *                                                                  *
// * This  code  implementation is the result of  the  scientific and *
// * technical work of the GEANT4 collaboration.                      *
// * By using,  copying,  modifying or  distributing the software (or *
// * any work based  on the software)  you  agree  to acknowledge its *
// * use  in  resulting  scientific  publications,  and indicate your *
// * acceptance of all terms of the Geant4 Software license.          *
// ********************************************************************
//

#include "PlhSteppingAction.hh"
#include "G4SteppingManager.hh"
#include "G4Track.hh"
#include "G4Step.hh"
#include "G4StepPoint.hh"
#include "G4TrackStatus.hh"
#include "G4TrackVector.hh"
#include "G4VPhysicalVolume.hh"
#include "G4ParticleDefinition.hh"
#include "G4ParticleTypes.hh"
#include "G4ios.hh"
#include <iomanip>

PlhSteppingAction::PlhSteppingAction()
{;}

PlhSteppingAction::~PlhSteppingAction()
{;}

void PlhSteppingAction::UserSteppingAction()
{
  G4SteppingManager * SM = GetOmnipotentSteppingManager();
  G4Step * theStep = SM->GetStep();
  G4Track * theTrack = theStep->GetTrack();

  // check if it is alive
  if(theTrack->GetTrackStatus()==fAlive) { return; }

      G4cout << std::setw( 5) << "#Step#" << " "
        << std::setw( 9) << "X(mm)" << " "
          << std::setw( 9) << "Y(mm)" << " "
            << std::setw( 9) << "Z(mm)" << " "
              << std::setw( 9) << "KineE(MeV)" << " "
                << std::setw( 9) << "dE(MeV)" << " "
                  << std::setw( 9) << "StepLeng" << " "
                    << std::setw( 9) << "TrackLeng" << " "
                        << std::setw(10) << "ProcName" << G4endl;
    G4cout.precision(3);
    G4cout << std::setw( 5) << theTrack->GetCurrentStepNumber() << " "
      << std::setw( 9) << theTrack->GetPosition().x() / mm << " "
        << std::setw( 9) << theTrack->GetPosition().y() / mm << " "
          << std::setw( 9) << theTrack->GetPosition().z() / mm << " "
             << std::setw( 9) << theTrack->GetKineticEnergy() / MeV << " "
              << std::setw( 9) << theStep->GetTotalEnergyDeposit() /MeV << " "
                << std::setw( 9) << theStep->GetStepLength() / mm << " "
                  << std::setw( 9) << theTrack->GetTrackLength() / mm << " ";
    if(theStep->GetPostStepPoint()->GetProcessDefinedStep() != NULL){
      G4cout << theStep->GetPostStepPoint()->GetProcessDefinedStep()
        ->GetProcessName();
    } else {
      G4cout << "User Limit";
    }
    G4cout << G4endl;

   G4TrackVector* fSecondary = SM->GetfSecondary();
       G4cout << "   -- List of secondaries generated : "
         << "(x,y,z,kE,t,PID) --" << G4endl;
       for( G4int lp1=0;lp1<(*fSecondary).entries(); lp1++){
         G4cout << "      "
           << std::setw( 9)
             << (*fSecondary)[lp1]->GetPosition().x() / mm << " "
               << std::setw( 9)
                 << (*fSecondary)[lp1]->GetPosition().y() / mm << " "
                   << std::setw( 9)
                     << (*fSecondary)[lp1]->GetPosition().z() / mm << " "
                       << std::setw( 9)
                         << (*fSecondary)[lp1]->GetKineticEnergy() / MeV << " "
                           << std::setw( 9)
                             << (*fSecondary)[lp1]->GetGlobalTime() / ns << " "
                               << std::setw(18)
                                 << (*fSecondary)[lp1]->GetDefinition()
                                   ->GetParticleName();
         G4cout << G4endl;
       }

}