TrackingVerboseAction Class Reference

#include <TrackingVerboseAction.h>

Inheritance diagram for TrackingVerboseAction:

Public Member Functions
	TrackingVerboseAction (edm::ParameterSet const &p)
void	update (const BeginOfEvent *) override
	This routine will be called when the appropriate signal arrives. More...
void	update (const BeginOfRun *) override
	This routine will be called when the appropriate signal arrives. More...
void	update (const BeginOfTrack *) override
	This routine will be called when the appropriate signal arrives. More...
void	update (const EndOfTrack *) override
	This routine will be called when the appropriate signal arrives. More...
void	update (const G4Step *) override
	This routine will be called when the appropriate signal arrives. More...
	~TrackingVerboseAction () override
Public Member Functions inherited from SimWatcher
	SimWatcher ()
virtual	~SimWatcher ()
Public Member Functions inherited from Observer< const BeginOfRun * >
	Observer ()
void	slotForUpdate (const BeginOfRun * iT)
virtual	~Observer ()
Public Member Functions inherited from Observer< const BeginOfEvent * >
	Observer ()
void	slotForUpdate (const BeginOfEvent * iT)
virtual	~Observer ()
Public Member Functions inherited from Observer< const BeginOfTrack * >
	Observer ()
void	slotForUpdate (const BeginOfTrack * iT)
virtual	~Observer ()
Public Member Functions inherited from Observer< const EndOfTrack * >
	Observer ()
void	slotForUpdate (const EndOfTrack * iT)
virtual	~Observer ()
Public Member Functions inherited from Observer< const G4Step * >
	Observer ()
void	slotForUpdate (const G4Step * iT)
virtual	~Observer ()

Private Member Functions
bool	checkTrackingVerbose (const G4Track *)
void	printTrackInfo (const G4Track *)
void	setTrackingVerbose (int verblev)

Private Attributes
bool	fDEBUG
bool	fG4Verbose
bool	fHighEtPhotons
int	fLarge
std::vector< int >	fPdgIds
bool	fTkVerbThisEventON
bool	fTrackingVerboseON
int	fTVEventMax
int	fTVEventMin
int	fTVEventStep
int	fTVTrackMax
int	fTVTrackMin
int	fTVTrackStep
G4VSteppingVerbose *	fVerbose
int	fVerboseLevel
G4TrackingManager *	theTrackingManager

Additional Inherited Members

Detailed Description

Definition at line 37 of file TrackingVerboseAction.h.

Constructor & Destructor Documentation

TrackingVerboseAction()

TrackingVerboseAction::TrackingVerboseAction

(

edm::ParameterSet const &

p

)

Definition at line 29 of file TrackingVerboseAction.cc.

    : theTrackingManager(nullptr), fVerbose(nullptr) {
  fLarge = int(1E10);
  fDEBUG = p.getUntrackedParameter<bool>("DEBUG", false);
  fHighEtPhotons = p.getUntrackedParameter<bool>("CheckForHighEtPhotons", false);
  fG4Verbose = p.getUntrackedParameter<bool>("G4Verbose", false);
 
  //----- Set which events are verbose
  fTVEventMin = p.getUntrackedParameter<int>("EventMin", 0);
  fTVEventMax = p.getUntrackedParameter<int>("EventMax", fLarge);
  fTVEventStep = p.getUntrackedParameter<int>("EventStep", 1);
 
  //----- Set which tracks of those events are verbose
  fTVTrackMin = p.getUntrackedParameter<int>("TrackMin", 0);
  fTVTrackMax = p.getUntrackedParameter<int>("TrackMax", fLarge);
  fTVTrackStep = p.getUntrackedParameter<int>("TrackStep", 1);
 
  //----- Set the verbosity level
  fVerboseLevel = p.getUntrackedParameter<int>("VerboseLevel", 1);
  fPdgIds = p.getUntrackedParameter<std::vector<int> >("PDGids");
  if (fDEBUG) {
    G4cout << "TV: fTVTrackMin " << fTVTrackMin << " fTVTrackMax " << fTVTrackMax << " fTVTrackStep " << fTVTrackStep
           << " fTVEventMin " << fTVEventMin << " fTVEventMax " << fTVEventMax << " fTVEventStep " << fTVEventStep
           << " fVerboseLevel " << fVerboseLevel << " fG4Verbose " << fG4Verbose << " PDGIds     " << fPdgIds.size()
           << G4endl;
    for (unsigned int ii = 0; ii < fPdgIds.size(); ++ii)
      G4cout << "TV: PDGId[" << ii << "] = " << fPdgIds[ii] << G4endl;
  }
 
  //----- Set verbosity off to start
  fTrackingVerboseON = false;
  fTkVerbThisEventON = false;
 
  G4cout << " TrackingVerbose constructed " << G4endl;
}

References fDEBUG, fG4Verbose, fHighEtPhotons, fLarge, fPdgIds, fTkVerbThisEventON, fTrackingVerboseON, fTVEventMax, fTVEventMin, fTVEventStep, fTVTrackMax, fTVTrackMin, fTVTrackStep, fVerboseLevel, ecalTB2006H4_GenSimDigiReco_cfg::G4cout, cuy::ii, createfilelist::int, and AlCaHLTBitMon_ParallelJobs::p.

~TrackingVerboseAction()

TrackingVerboseAction::~TrackingVerboseAction ( )

override

Definition at line 65 of file TrackingVerboseAction.cc.

{}

Member Function Documentation

checkTrackingVerbose()

bool TrackingVerboseAction::checkTrackingVerbose ( const G4Track *  aTrack )

private

Definition at line 279 of file TrackingVerboseAction.cc.

                                                                      {
  int trackNo = aTrack->GetTrackID();
  bool trackingVerboseThisTrack = false;
  //----- Check if track is in the selected range
  if (trackNo >= fTVTrackMin && trackNo <= fTVTrackMax) {
    if ((trackNo - fTVTrackMin) % fTVTrackStep == 0)
      trackingVerboseThisTrack = true;
  }
  if (trackingVerboseThisTrack && (!fPdgIds.empty())) {
    int pdgId = aTrack->GetDefinition()->GetPDGEncoding();
    if (std::count(fPdgIds.begin(), fPdgIds.end(), pdgId) == 0)
      trackingVerboseThisTrack = false;
  }
  return trackingVerboseThisTrack;
}

References submitPVResolutionJobs::count, fPdgIds, fTVTrackMax, fTVTrackMin, fTVTrackStep, and EgammaValidation_cff::pdgId.

Referenced by update().

printTrackInfo()

void TrackingVerboseAction::printTrackInfo ( const G4Track *  aTrack )

private

Definition at line 295 of file TrackingVerboseAction.cc.

                                                                {
  G4cout << G4endl << "*******************************************************"
         << "**************************************************" << G4endl << "* G4Track Information: "
         << "  Particle = " << aTrack->GetDefinition()->GetParticleName() << ","
         << "   Track ID = " << aTrack->GetTrackID() << ","
         << "   Parent ID = " << aTrack->GetParentID() << G4endl
         << "*******************************************************"
         << "**************************************************" << G4endl << G4endl;
  if (fVerbose)
    fVerbose->TrackingStarted();
}

References fVerbose, and ecalTB2006H4_GenSimDigiReco_cfg::G4cout.

Referenced by update().

setTrackingVerbose()

void TrackingVerboseAction::setTrackingVerbose ( int  verblev )

private

Definition at line 272 of file TrackingVerboseAction.cc.

                                                          {
  if (fDEBUG)
    G4cout << " setting verbose level " << verblev << G4endl;
  if (theTrackingManager != nullptr)
    theTrackingManager->SetVerboseLevel(verblev);
}

References fDEBUG, ecalTB2006H4_GenSimDigiReco_cfg::G4cout, and theTrackingManager.

Referenced by update().

update() [1/5]

void TrackingVerboseAction::update ( const BeginOfEvent *  )

overridevirtual

This routine will be called when the appropriate signal arrives.

Implements Observer< const BeginOfEvent * >.

Definition at line 79 of file TrackingVerboseAction.cc.

                                                          {
  if (evt == nullptr)
    return;
  const G4Event* anEvent = (*evt)();
  if (anEvent == nullptr)
    return;
 
  //----------- Set /tracking/verbose for this event
  int eventNo = anEvent->GetEventID();
  if (fDEBUG)
    G4cout << "TV: trackID: NEW EVENT " << eventNo << G4endl;
 
  fTkVerbThisEventON = false;
  //----- Check if event is in the selected range
  if (eventNo >= fTVEventMin && eventNo <= fTVEventMax) {
    if ((eventNo - fTVEventMin) % fTVEventStep == 0)
      fTkVerbThisEventON = true;
  }
 
  if (fDEBUG)
    G4cout << " TV: fTkVerbThisEventON " << fTkVerbThisEventON << " fTrackingVerboseON " << fTrackingVerboseON
           << " fTVEventMin " << fTVEventMin << " fTVEventMax " << fTVEventMax << G4endl;
  //----- check if verbosity has to be changed
  if ((fTkVerbThisEventON) && (!fTrackingVerboseON)) {
    if (fTVTrackMin == 0 && fTVTrackMax == fLarge && fTVTrackStep != 1) {
      setTrackingVerbose(fVerboseLevel);
      fTrackingVerboseON = true;
      if (fDEBUG)
        G4cout << "TV: VERBOSEet1 " << eventNo << G4endl;
    }
  } else if ((!fTkVerbThisEventON) && (fTrackingVerboseON)) {
    setTrackingVerbose(0);
    fTrackingVerboseON = false;
    if (fDEBUG)
      G4cout << "TV: VERBOSEet0 " << eventNo << G4endl;
  }
}

References fDEBUG, fLarge, fTkVerbThisEventON, fTrackingVerboseON, fTVEventMax, fTVEventMin, fTVEventStep, fTVTrackMax, fTVTrackMin, fTVTrackStep, fVerboseLevel, ecalTB2006H4_GenSimDigiReco_cfg::G4cout, and setTrackingVerbose().

Referenced by progressbar.ProgressBar::__next__(), MatrixUtil.Matrix::__setitem__(), MatrixUtil.Steps::__setitem__(), progressbar.ProgressBar::finish(), and MatrixUtil.Steps::overwrite().

update() [2/5]

void TrackingVerboseAction::update ( const BeginOfRun *  )

overridevirtual

This routine will be called when the appropriate signal arrives.

Implements Observer< const BeginOfRun * >.

Definition at line 67 of file TrackingVerboseAction.cc.

                                                        {
  /*
  TrackingAction * ta = 
    dynamic_cast<TrackingAction*>(G4EventManager::GetEventManager()->GetUserTrackingAction());
  theTrackingManager = ta->getTrackManager();
  */
  fVerbose = G4VSteppingVerbose::GetInstance();
  if (fDEBUG)
    G4cout << " TV: Get the Tracking Manager: " << theTrackingManager << " and the SteppingVerbose: " << fVerbose
           << G4endl;
}

References fDEBUG, fVerbose, ecalTB2006H4_GenSimDigiReco_cfg::G4cout, and theTrackingManager.

Referenced by progressbar.ProgressBar::__next__(), MatrixUtil.Matrix::__setitem__(), MatrixUtil.Steps::__setitem__(), progressbar.ProgressBar::finish(), and MatrixUtil.Steps::overwrite().

update() [3/5]

void TrackingVerboseAction::update ( const BeginOfTrack *  )

overridevirtual

This routine will be called when the appropriate signal arrives.

Implements Observer< const BeginOfTrack * >.

Definition at line 117 of file TrackingVerboseAction.cc.

                                                          {
  const G4Track* aTrack = (*trk)();
 
  //----- High ET photon printout
  //---------- Set /tracking/verbose
  //----- track is verbose only if event is verbose
  double tkP = aTrack->GetMomentum().mag();
  double tkPx = aTrack->GetMomentum().x();
  double tkPy = aTrack->GetMomentum().y();
  double tkPz = aTrack->GetMomentum().z();
 
  double tvtx = aTrack->GetVertexPosition().x();
  double tvty = aTrack->GetVertexPosition().y();
  double tvtz = aTrack->GetVertexPosition().z();
 
  double g4t_phi = atan2(tkPy, tkPx);
 
  double drpart = sqrt(tkPx * tkPx + tkPy * tkPy);
 
  double mythetapart = acos(tkPz / sqrt(drpart * drpart + tkPz * tkPz));
 
  double g4t_eta = -log(tan(mythetapart / 2.));
  G4int MytrackNo = aTrack->GetTrackID();
 
  if (fHighEtPhotons) {
    if (aTrack->GetDefinition()->GetParticleName() == "gamma" && aTrack->GetParentID() != 0) {
      if ((tkPx * tkPx + tkPy * tkPy + tkPz * tkPz) > 1000.0 * 1000.0 &&
          aTrack->GetCreatorProcess()->GetProcessName() == "LCapture") {
        G4cout << "MY NEW GAMMA " << G4endl;
        G4cout << "**********************************************************************" << G4endl;
        G4cout << "MY NEW TRACK ID = " << MytrackNo << "(" << aTrack->GetDefinition()->GetParticleName() << ")"
               << " PARENT =" << aTrack->GetParentID() << G4endl;
        G4cout << "Primary particle: " << aTrack->GetDynamicParticle()->GetPrimaryParticle() << G4endl;
        G4cout << "Process type: " << aTrack->GetCreatorProcess()->GetProcessType()
               << " Process name: " << aTrack->GetCreatorProcess()->GetProcessName() << G4endl;
        G4cout << "ToT E = " << aTrack->GetTotalEnergy() << " KineE = " << aTrack->GetKineticEnergy()
               << " Tot P = " << tkP << " Pt = " << sqrt(tkPx * tkPx + tkPy * tkPy) << " VTX=(" << tvtx << "," << tvty
               << "," << tvtz << ")" << G4endl;
        if (aTrack->GetKineticEnergy() > 1. * GeV && aTrack->GetCreatorProcess()->GetProcessName() != "LCapture")
          G4cout << " KineE > 1 GeV !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!" << G4endl;
        const G4VTouchable* touchable = aTrack->GetTouchable();
        if (touchable != nullptr && touchable->GetVolume() != nullptr &&
            touchable->GetVolume()->GetLogicalVolume() != nullptr) {
          G4Material* material = touchable->GetVolume()->GetLogicalVolume()->GetMaterial();
          G4cout << "G4LCapture Gamma E(GeV) " << aTrack->GetTotalEnergy() / GeV << "  " << material->GetName() << " "
                 << touchable->GetVolume()->GetName() << G4endl;
          G4cout << "G4LCapture Gamma position(m): " << aTrack->GetPosition() / m << G4endl;
          G4cout << "G4LCapture created Gamma direction " << aTrack->GetMomentumDirection() << G4endl;
          G4cout << "G4LCapture gamma (eta,phi) = "
                 << "(" << g4t_eta << "," << g4t_phi << ")" << G4endl;
        }
        aTrack->GetUserInformation()->Print();
        G4cout << "**********************************************************************" << G4endl;
      }
    }
 
    if (aTrack->GetDefinition()->GetParticleName() == "gamma") {
      const G4VProcess* proc = aTrack->GetCreatorProcess();
      double Tgamma = aTrack->GetKineticEnergy();
      std::string ProcName;
      const std::string nullstr("Null_prc");
      if (proc)
        ProcName = proc->GetProcessName();
      else
        ProcName = nullstr;
      if (Tgamma > 2.5 * GeV) {  //&& ProcName!="Decay" && ProcName!="eBrem")
        std::string volumeName("_Unknown_Vol_");
        std::string materialName("_Unknown_Mat_");
        G4Material* material = nullptr;
        G4VPhysicalVolume* pvolume = nullptr;
        G4LogicalVolume* lvolume = nullptr;
        const G4VTouchable* touchable = aTrack->GetTouchable();
        if (touchable)
          pvolume = touchable->GetVolume();
        if (pvolume) {
          volumeName = pvolume->GetName();
          lvolume = pvolume->GetLogicalVolume();
        }
        if (lvolume)
          material = lvolume->GetMaterial();
        if (material)
          materialName = material->GetName();
        G4cout << "**** ALL photons > 2.5 GeV ****" << G4endl;
        G4cout << ProcName << "**** ALL photons: gamma E(GeV) " << aTrack->GetTotalEnergy() / GeV << "  "
               << materialName << " " << volumeName << G4endl;
        G4cout << ProcName << "**** ALL photons: gamma position(m): " << aTrack->GetPosition() / m << G4endl;
        G4cout << ProcName << "**** ALL photons: gamma direction " << aTrack->GetMomentumDirection() << G4endl;
        G4cout << "**********************************************************************" << G4endl;
      }
    }
  }
 
  //---------- Set /tracking/verbose
  //----- track is verbose only if event is verbose
  if (fTkVerbThisEventON) {
    bool trackingVerboseThisTrack = checkTrackingVerbose(aTrack);
 
    //----- Set the /tracking/verbose for this track
    if ((trackingVerboseThisTrack) && (!fTrackingVerboseON)) {
      setTrackingVerbose(fVerboseLevel);
      fTrackingVerboseON = true;
      if (fDEBUG)
        G4cout << "TV: VERBOSEtt1 " << aTrack->GetTrackID() << G4endl;
      printTrackInfo(aTrack);
    } else if ((!trackingVerboseThisTrack) && (fTrackingVerboseON)) {
      setTrackingVerbose(0);
      fTrackingVerboseON = false;
      if (fDEBUG)
        G4cout << "TV: VERBOSEtt0 " << aTrack->GetTrackID() << G4endl;
    }
  }
}

References checkTrackingVerbose(), fDEBUG, fHighEtPhotons, fTkVerbThisEventON, fTrackingVerboseON, fVerboseLevel, ecalTB2006H4_GenSimDigiReco_cfg::G4cout, GeV, dqm-mbProfile::log, visualization-live-secondInstance_cfg::m, printTrackInfo(), ValidateTausOnZEEFastSim_cff::proc, setTrackingVerbose(), mathSSE::sqrt(), AlCaHLTBitMon_QueryRunRegistry::string, and funct::tan().

Referenced by progressbar.ProgressBar::__next__(), MatrixUtil.Matrix::__setitem__(), MatrixUtil.Steps::__setitem__(), progressbar.ProgressBar::finish(), and MatrixUtil.Steps::overwrite().

update() [4/5]

void TrackingVerboseAction::update ( const EndOfTrack *  )

overridevirtual

This routine will be called when the appropriate signal arrives.

Implements Observer< const EndOfTrack * >.

Definition at line 230 of file TrackingVerboseAction.cc.

                                                        {
  const G4Track* aTrack = (*trk)();
  if (fTkVerbThisEventON) {
    bool trackingVerboseThisTrack = checkTrackingVerbose(aTrack);
    if ((trackingVerboseThisTrack) && (fTrackingVerboseON) && (fTVTrackMax < fLarge || fTVTrackStep != 1)) {
      setTrackingVerbose(0);
      fTrackingVerboseON = false;
      if (fDEBUG)
        G4cout << "TV: VERBOSEtt0 " << aTrack->GetTrackID() << G4endl;
    }
  }
}

References checkTrackingVerbose(), fDEBUG, fLarge, fTkVerbThisEventON, fTrackingVerboseON, fTVTrackMax, fTVTrackStep, ecalTB2006H4_GenSimDigiReco_cfg::G4cout, and setTrackingVerbose().

Referenced by progressbar.ProgressBar::__next__(), MatrixUtil.Matrix::__setitem__(), MatrixUtil.Steps::__setitem__(), progressbar.ProgressBar::finish(), and MatrixUtil.Steps::overwrite().

update() [5/5]

void TrackingVerboseAction::update ( const G4Step *  )

overridevirtual

This routine will be called when the appropriate signal arrives.

Implements Observer< const G4Step * >.

Definition at line 243 of file TrackingVerboseAction.cc.

                                                      {
  if ((fG4Verbose) && (fTrackingVerboseON)) {
    G4Track* fTrack = fStep->GetTrack();
    G4cout << std::setw(5) << fTrack->GetCurrentStepNumber() << " " << std::setw(8)
           << G4BestUnit(fTrack->GetPosition().x(), "Length") << " " << std::setw(8)
           << G4BestUnit(fTrack->GetPosition().y(), "Length") << " " << std::setw(8)
           << G4BestUnit(fTrack->GetPosition().z(), "Length") << " " << std::setw(9)
           << G4BestUnit(fTrack->GetKineticEnergy(), "Energy") << " " << std::setw(8)
           << G4BestUnit(fStep->GetTotalEnergyDeposit(), "Energy") << " " << std::setw(8)
           << G4BestUnit(fStep->GetStepLength(), "Length") << " " << std::setw(9)
           << G4BestUnit(fTrack->GetTrackLength(), "Length") << " " << std::setw(9)
           << G4BestUnit(fTrack->GetGlobalTime(), "Time") << " ";
 
    // Put cut comment here
    if (fTrack->GetNextVolume() != nullptr) {
      G4cout << std::setw(11) << fTrack->GetNextVolume()->GetName() << " ";
    } else {
      G4cout << std::setw(11) << "OutOfWorld"
             << " ";
    }
    if (fStep->GetPostStepPoint()->GetProcessDefinedStep() != nullptr) {
      G4cout << fStep->GetPostStepPoint()->GetProcessDefinedStep()->GetProcessName();
    } else {
      G4cout << "User Limit";
    }
    G4cout << G4endl;
  }
}

References fG4Verbose, fTrackingVerboseON, and ecalTB2006H4_GenSimDigiReco_cfg::G4cout.

Referenced by progressbar.ProgressBar::__next__(), MatrixUtil.Matrix::__setitem__(), MatrixUtil.Steps::__setitem__(), progressbar.ProgressBar::finish(), and MatrixUtil.Steps::overwrite().

Member Data Documentation

fDEBUG

bool TrackingVerboseAction::fDEBUG

private

Definition at line 59 of file TrackingVerboseAction.h.

Referenced by setTrackingVerbose(), TrackingVerboseAction(), and update().

fG4Verbose

bool TrackingVerboseAction::fG4Verbose

private

Definition at line 60 of file TrackingVerboseAction.h.

Referenced by TrackingVerboseAction(), and update().

fHighEtPhotons

bool TrackingVerboseAction::fHighEtPhotons

private

Definition at line 61 of file TrackingVerboseAction.h.

Referenced by TrackingVerboseAction(), and update().

fLarge

int TrackingVerboseAction::fLarge

private

Definition at line 58 of file TrackingVerboseAction.h.

Referenced by TrackingVerboseAction(), and update().

fPdgIds

std::vector<int> TrackingVerboseAction::fPdgIds

private

Definition at line 71 of file TrackingVerboseAction.h.

Referenced by checkTrackingVerbose(), and TrackingVerboseAction().

fTkVerbThisEventON

bool TrackingVerboseAction::fTkVerbThisEventON

private

Definition at line 70 of file TrackingVerboseAction.h.

Referenced by TrackingVerboseAction(), and update().

fTrackingVerboseON

bool TrackingVerboseAction::fTrackingVerboseON

private

Definition at line 69 of file TrackingVerboseAction.h.

Referenced by TrackingVerboseAction(), and update().

fTVEventMax

int TrackingVerboseAction::fTVEventMax

private

Definition at line 66 of file TrackingVerboseAction.h.

Referenced by TrackingVerboseAction(), and update().

fTVEventMin

int TrackingVerboseAction::fTVEventMin

private

Definition at line 65 of file TrackingVerboseAction.h.

Referenced by TrackingVerboseAction(), and update().

fTVEventStep

int TrackingVerboseAction::fTVEventStep

private

Definition at line 67 of file TrackingVerboseAction.h.

Referenced by TrackingVerboseAction(), and update().

fTVTrackMax

int TrackingVerboseAction::fTVTrackMax

private

Definition at line 63 of file TrackingVerboseAction.h.

Referenced by checkTrackingVerbose(), TrackingVerboseAction(), and update().

fTVTrackMin

int TrackingVerboseAction::fTVTrackMin

private

Definition at line 62 of file TrackingVerboseAction.h.

Referenced by checkTrackingVerbose(), TrackingVerboseAction(), and update().

fTVTrackStep

int TrackingVerboseAction::fTVTrackStep

private

Definition at line 64 of file TrackingVerboseAction.h.

Referenced by checkTrackingVerbose(), TrackingVerboseAction(), and update().

fVerbose

G4VSteppingVerbose* TrackingVerboseAction::fVerbose

private

Definition at line 73 of file TrackingVerboseAction.h.

Referenced by printTrackInfo(), and update().

fVerboseLevel

int TrackingVerboseAction::fVerboseLevel

private

Definition at line 68 of file TrackingVerboseAction.h.

Referenced by TrackingVerboseAction(), and update().

theTrackingManager

G4TrackingManager* TrackingVerboseAction::theTrackingManager

private

Definition at line 72 of file TrackingVerboseAction.h.

Referenced by setTrackingVerbose(), and update().