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);
}