MuonSpectrometer/MuonPhaseII/MuonG4/MuonSensitiveDetectorsR4/MuonSensitiveDetectorsR4/Utils.h

Go to the documentation of this file.

/*
  Copyright (C) 2002-2023 CERN for the benefit of the ATLAS collaboration
*/
#ifndef MUONSENSITIVEDETECTORSR4_UTILS_H
#define MUONSENSITIVEDETECTORSR4_UTILS_H
#include <MuonReadoutGeometryR4/MuonDetectorDefs.h>
 
 
#include <G4VTouchable.hh>
#include <G4Track.hh>
 
 
namespace MuonG4R4 {
 
 
    inline Amg::Transform3D getTransform(const G4VTouchable* history, unsigned int level) {
            return Amg::Translation3D{Amg::Hep3VectorToEigen(history->GetTranslation(level))}*
                                      Amg::CLHEPRotationToEigen(*history->GetRotation(level)).inverse();
    }
}
    inline std::ostream& operator<<(std::ostream& ostr, const G4Track& step) {
        ostr<<"position: "<<Amg::toString(Amg::Hep3VectorToEigen(step.GetPosition()),2)<<", ";
        ostr<<"momentum: "<<Amg::toString(Amg::Hep3VectorToEigen(step.GetMomentum()),2)<<", ";
        ostr<<"velocity: "<<step.GetVelocity()<<", ";
        ostr<<"time: "<<step.GetGlobalTime()<<", ";
        ostr<<"mass: "<<step.GetDefinition()->GetPDGMass()<<", ";
        ostr<<"kinetic energy: "<<step.GetKineticEnergy()<<", ";
        ostr<<"charge: "<<step.GetDefinition()->GetPDGCharge();
        return ostr;
    }
    inline std::ostream& operator<<(std::ostream& ostr, const G4StepPoint& step) {
        ostr<<"position: "<<Amg::toString(Amg::Hep3VectorToEigen(step.GetPosition()),2)<<", ";
        ostr<<"momentum: "<<Amg::toString(Amg::Hep3VectorToEigen(step.GetMomentum()),2)<<", ";
        ostr<<"time: "<<step.GetGlobalTime()<<", ";
        ostr<<"mass: "<<step.GetMass()<<", ";
        ostr<<"kinetic energy: "<<step.GetKineticEnergy()<<", ";
        ostr<<"charge: "<<step.GetCharge();
        return ostr;
    }
 
 
#endif