GenericMuonSensitiveDetector.cxx

Go to the documentation of this file.

/*
  Copyright (C) 2002-2019 CERN for the benefit of the ATLAS collaboration
*/
 
#include "GenericMuonSensitiveDetector.h"
#include "MCTruth/AtlasG4EventUserInfo.h"
#include "MCTruth/TrackHelper.h"
 
 
#include "G4Exception.hh"
#include "G4Track.hh"
 
#include <string>
 
#include "GeoPrimitives/CLHEPtoEigenConverter.h"
 
// construction/destruction
GenericMuonSensitiveDetector::GenericMuonSensitiveDetector(const std::string& name, const std::string& hitCollectionName)
  : G4VSensitiveDetector( name )
  , m_hitCollectionName( hitCollectionName )
{
  G4cout << " creating a GenericMuonSensitiveDetector: "<<name << G4endl;
}
 
// Implemenation of memebr functions
void GenericMuonSensitiveDetector::Initialize(G4HCofThisEvent*) 
{
  m_GenericMuonHitCollection = nullptr;
  if (auto* eventInfo = AtlasG4EventUserInfo::GetEventUserInfo()) {
    m_GenericMuonHitCollection = eventInfo->GetHitCollectionMap()->Find<GenericMuonSimHitCollection>(m_hitCollectionName);
    m_g4UserEventInfo = eventInfo;
  }
}
 
G4bool GenericMuonSensitiveDetector::ProcessHits(G4Step* aStep,G4TouchableHistory* /*ROHist*/) 
{
  if (!m_GenericMuonHitCollection) {
    G4Exception("GenericMuonSensitiveDetector::ProcessHits", "GenericMuonHitCollectionMissing", FatalException,
                "Hit collection not initialized; did SetupEvent run?");
    return false;
  }
  G4cout << "Hit in a sensitive layer!!!!! " << G4endl;
  G4Track* currentTrack = aStep->GetTrack();
  const G4AffineTransform trans = currentTrack->GetTouchable()->GetHistory()->GetTopTransform(); // from global to local
  G4StepPoint* postStep=aStep->GetPostStepPoint();
  G4StepPoint* preStep=aStep->GetPreStepPoint();
  const G4Step* post_Step=aStep->GetTrack()->GetStep();
 
  Amg::Vector3D position = Amg::Hep3VectorToEigen( postStep->GetPosition() );
  Amg::Vector3D preposition = Amg::Hep3VectorToEigen( preStep->GetPosition() );
  
  Amg::Vector3D local_position = Amg::Hep3VectorToEigen( trans.TransformPoint( postStep->GetPosition() ) ) ;
  Amg::Vector3D local_preposition = Amg::Hep3VectorToEigen( trans.TransformPoint( preStep->GetPosition() ) );
  
  int pdgCode=currentTrack->GetDefinition()->GetPDGEncoding();
 
  float globalTime=postStep->GetGlobalTime();
  float globalpreTime=preStep->GetGlobalTime();
  float eKin=postStep->GetKineticEnergy();
 
  Amg::Vector3D direction = Amg::Hep3VectorToEigen( postStep->GetMomentumDirection() );
  float depositEnergy=post_Step->GetTotalEnergyDeposit();
  float StepLength=post_Step->GetStepLength();
 
  TrackHelper trHelp(aStep->GetTrack());
 
  //G4cout << aHit->print() << G4endl;
  m_GenericMuonHitCollection->Emplace( 0 /* HitID id generic*/,globalTime,globalpreTime,position,local_position,preposition,local_preposition,pdgCode,eKin,direction,depositEnergy,StepLength,
                                      trHelp.GenerateParticleLink(m_g4UserEventInfo ? m_g4UserEventInfo->GetEventStore() : nullptr));
 
  return true;
}