SimG4Core/Notification/src/TrackInformation.cc - master - githubmirror/cmssw

SimG4Core/Notification/src/TrackInformation.cc
46 строк3 KB
Fabio Cossutti
Switch off verbosity
17 дек 2024, 15:49
17 дек 2024, 15:4948ca33e
//#define EDM_ML_DEBUG

#include "SimG4Core/Notification/interface/TrackInformation.h"
#include "FWCore/MessageLogger/interface/MessageLogger.h"

#include "G4ThreeVector.hh"
#include <CLHEP/Units/SystemOfUnits.h>

#include <iostream>

G4ThreadLocal G4Allocator<TrackInformation>* fpTrackInformationAllocator = nullptr;

void TrackInformation::setCrossedBoundary(const G4Track* track) {
  const double invcm = 1.0 / CLHEP::cm;
  const double invgev = 1.0 / CLHEP::GeV;
  const double invsec = 1.0 / CLHEP::second;
  crossedBoundary_ = true;
  const G4ThreeVector& v = track->GetPosition();
  positionAtBoundary_ =
      math::XYZTLorentzVectorF(v.x() * invcm, v.y() * invcm, v.z() * invcm, track->GetGlobalTime() * invsec);
  const G4ThreeVector& p = track->GetMomentum();
  momentumAtBoundary_ =
      math::XYZTLorentzVectorF(p.x() * invgev, p.y() * invgev, p.z() * invgev, track->GetTotalEnergy() * invgev);
}

void TrackInformation::Print() const {
  LogDebug("TrackInformation") << " TrackInformation : storeTrack = " << storeTrack_ << "\n"
                               << "                    hasHits = " << hasHits_ << "\n"
                               << "                    isPrimary = " << isPrimary_ << "\n"
                               << "                    isGeneratedSecondary = " << isGeneratedSecondary_ << "\n"
                               << "                    mcTruthID = " << mcTruthID_ << "\n"
                               << "                    isInHistory = " << isInHistory_ << "\n"
                               << "                    idLastStoredAncestor = " << idLastStoredAncestor() << "\n"
                               << "                    idOnCaloSurface = " << getIDonCaloSurface() << "\n"
                               << "                    caloIDChecked = " << caloIDChecked() << "\n"
                               << "                    idCaloVolume = " << idCaloVolume_ << "\n"
                               << "                    idLastVolume = " << idLastVolume_ << "\n"
                               << "                    isFromTtoBTL = " << isFromTtoBTL() << "\n"
                               << "                    isFromBTLtoT = " << isFromBTLtoT() << "\n"
                               << "                    isBTLlooper = " << isBTLlooper() << "\n"
                               << "                    isInTrkFromBackscattering = " << isInTrkFromBackscattering()
                               << "\n"
                               << "                    isExtSecondary = " << isExtSecondary() << "\n"
                               << "                    isETLFromFront = " << isETLfromFront() << "\n"
                               << "                    isETLFromBack = " << isETLfromBack();
}