---

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() != 0 )
            {
              G4cout << std::setw(10) << fTrack->GetVolume()->GetName();
            }
          else
            {
              G4cout << std::setw(10) << "OutOfWorld";
            }

          if ( fStep->GetPostStepPoint()->GetProcessDefinedStep() != 0 )
            {
              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() != 0)
        {
          G4cout << std::setw(10) << fTrack->GetVolume()->GetName();
        }
      else
        {
          G4cout << "OutOfWorld";
        }
      
      G4cout << "     initStep" << G4endl;
    }
  G4cout.precision(precision);
}

---