MaterialAccountingTrack.cc

Go to the documentation of this file.

#include <iostream>

#include "DataFormats/GeometryVector/interface/GlobalPoint.h"
#include "G4VPhysicalVolume.hh"
#include "SimDataFormats/ValidationFormats/interface/MaterialAccountingTrack.h"

void MaterialAccountingTrack::enterDetector( const G4VPhysicalVolume* volume, const GlobalPoint& position, double cosTheta ) {
  if (m_current_volume != nullptr) {
    // error: entering a volume while inside an other (or the same) one !
    if (m_current_volume != volume)
      std::cerr << "MaterialAccountingTrack::leaveDetector(...): ERROR: entering volume (" << volume << ") while inside volume (" << m_current_volume << ")" << std::endl;
    else
      std::cerr << "MaterialAccountingTrack::leaveDetector(...): ERROR: entering volume (" << volume << ") twice" << std::endl;
    m_detector.clear();
    return;
  } 

  m_current_volume            = volume;
  m_detector.m_position       = position;
  m_detector.m_curvilinearIn  = m_total.length();
  m_detector.m_cosThetaIn     = cosTheta;
}

void MaterialAccountingTrack::leaveDetector( const G4VPhysicalVolume* volume, double cosTheta ) {
  if (m_current_volume != volume) {
    // error: leaving the wrong (or no) volume !
    if (m_current_volume)
      std::cerr << "MaterialAccountingTrack::leaveDetector(...): ERROR: leaving volume (" << volume << ") while inside volume (" << m_current_volume << ")" << std::endl;
    else
      std::cerr << "MaterialAccountingTrack::leaveDetector(...): ERROR: leaving volume (" << volume << ") while not inside any volume" << std::endl;
    m_detector.clear();
    return;
  } 

  m_current_volume = nullptr;
  m_detector.m_curvilinearOut = m_total.length();
  m_detector.m_cosThetaOut    = cosTheta;
  m_detectors.push_back( m_detector );
  m_detector.clear();
}