LaserSteppingVerbose.cc

Go to the documentation of this file.

#include "Alignment/LaserAlignmentSimulation/interface/LaserSteppingVerbose.h"
#include "G4UnitsTable.hh"

LaserSteppingVerbose::LaserSteppingVerbose()
{
}

LaserSteppingVerbose::~LaserSteppingVerbose()
{
}

void LaserSteppingVerbose::StepInfo()
{
  //  G4cout << "<LaserSteppingVerbose::StepInfo()>: verboseLevel = " << verboseLevel << endl;

  CopyState();

  G4int precision = G4cout.precision(3);

  if ( verboseLevel >= 1 )
    {
      if ( verboseLevel >= 4 ) VerboseTrack();
      if ( verboseLevel >= 3 )
    {
      G4cout << G4endl;
      G4cout << std::setw( 5) << "#Step#"      << " "
         << std::setw( 6) << "X"           << "   "
         << std::setw( 6) << "Y"           << "   "
         << std::setw( 6) << "Z"           << "   "
         << std::setw( 9) << "KineE"       << "   "
         << std::setw( 9) << "dEStep"      << "   "
         << std::setw(10) << "StepLength"  << " "
         << std::setw(10) << "TrackLength" << " "
         << std::setw(10) << "Volume"      << " "
         << std::setw(10) << "Process"     << G4endl;

      G4cout << std::setw( 5) << fTrack->GetCurrentStepNumber() << " "
         << std::setw( 6) << G4BestUnit(fTrack->GetPosition().x(),"Length")
         << std::setw( 6) << G4BestUnit(fTrack->GetPosition().y(),"Length")
         << std::setw( 6) << G4BestUnit(fTrack->GetPosition().z(),"Length")
         << std::setw( 6) << G4BestUnit(fTrack->GetKineticEnergy(),"Energy")
         << std::setw( 6) << G4BestUnit(fStep->GetTotalEnergyDeposit(),"Energy")
         << std::setw( 6) << G4BestUnit(fStep->GetStepLength(),"Length")
         << std::setw( 6) << G4BestUnit(fTrack->GetTrackLength(),"Length")
         << " ";

      if ( fTrack->GetNextVolume() != nullptr )
        {
          G4cout << std::setw(10) << fTrack->GetVolume()->GetName();
        }
      else
        {
          G4cout << std::setw(10) << "OutOfWorld";
        }

      if ( fStep->GetPostStepPoint()->GetProcessDefinedStep() != nullptr )
        {
          G4cout << " " << std::setw(10) << fStep->GetPostStepPoint()->GetProcessDefinedStep()
        ->GetProcessName();
        }
      else
        {
          G4cout << "    UserLimit";
        }

      G4cout << G4endl;

      if ( verboseLevel == 2 )
        {
          // total number of secondaries
          G4int tN2ndariesTot = fN2ndariesAtRestDoIt +
        fN2ndariesAlongStepDoIt + fN2ndariesPostStepDoIt;

          if ( tN2ndariesTot > 0 ) 
        {
          G4cout << "   :---- List of Secondaries - "
             << "#SpawnInStep = " << std::setw(3) << tN2ndariesTot
             << "(Rest = " << std::setw(2) << fN2ndariesAtRestDoIt
             << ", Along = " << std::setw(2) << fN2ndariesAlongStepDoIt
             << ", Post = " << std::setw(2) << fN2ndariesPostStepDoIt
             << "), "
             << "#SpawnTotal = " << std::setw(3) << (*fSecondary).size()
             << " --------- " << G4endl;

          for ( size_t lp1 = (*fSecondary).size() - tN2ndariesTot; lp1 < (*fSecondary).size(); lp1++)
            {
              G4cout << "   : " 
                 << std::setw(6) 
                 << G4BestUnit((*fSecondary)[lp1]->GetPosition().x(),"Length")
                 << std::setw(6)
                 << G4BestUnit((*fSecondary)[lp1]->GetPosition().y(),"Length")
                 << std::setw(6)
                 << G4BestUnit((*fSecondary)[lp1]->GetPosition().z(),"Length")
                 << std::setw(6)
                 << G4BestUnit((*fSecondary)[lp1]->GetKineticEnergy(),"Energy")
                 << std::setw(10)
                 << (*fSecondary)[lp1]->GetDefinition()->GetParticleName();
              G4cout << G4endl;
            }

          G4cout << "   :----------------------"
             << "--------------------------"
             << "-- End of Secondaries Info --------------------------------- "
             << G4endl;
        }
        }
    }
    }
  G4cout.precision(precision);
}

void LaserSteppingVerbose::TrackingStarted()
{
  CopyState();
  
  G4int precision = G4cout.precision(3);
  if ( verboseLevel > 0 )
    {
      G4cout << std::setw( 5) << "Step#"       << " " 
         << std::setw( 6) << "X"           << "     "
         << std::setw( 6) << "Y"           << "     "
         << std::setw( 6) << "Z"           << "     "
         << std::setw( 9) << "KineE"       << " "
         << std::setw( 9) << "dEStep"      << " "
         << std::setw(10) << "StepLength"  << " "
         << std::setw(10) << "TrackLength" << " "
         << std::setw(10) << "Volume"      << "  "
         << std::setw(10) << "Process "    << G4endl;
      
      G4cout << std::setw( 5) << fTrack->GetCurrentStepNumber() << " "
         << std::setw( 6) << G4BestUnit(fTrack->GetPosition().x(),"Length")
         << std::setw( 6) << G4BestUnit(fTrack->GetPosition().y(),"Length")
         << std::setw( 6) << G4BestUnit(fTrack->GetPosition().z(),"Length")
         << std::setw( 6) << G4BestUnit(fTrack->GetKineticEnergy(),"Energy")
         << std::setw( 6) << G4BestUnit(fStep->GetTotalEnergyDeposit(),"Energy")
         << std::setw( 6) << G4BestUnit(fStep->GetStepLength(),"Length")
         << std::setw( 6) << G4BestUnit(fTrack->GetTrackLength(),"Length") << "   ";
      
      if (fTrack->GetNextVolume() != nullptr)
    {
      G4cout << std::setw(10) << fTrack->GetVolume()->GetName();
    }
      else
    {
      G4cout << "OutOfWorld";
    }
      
      G4cout << "     initStep" << G4endl;
    }
  G4cout.precision(precision);
}