ChargeExchangeMC/src/CexmcSteppingAction.cc

0001 //
0002 // ********************************************************************
0003 // * License and Disclaimer                                           *
0004 // *                                                                  *
0005 // * The  Geant4 software  is  copyright of the Copyright Holders  of *
0006 // * the Geant4 Collaboration.  It is provided  under  the terms  and *
0007 // * conditions of the Geant4 Software License,  included in the file *
0008 // * LICENSE and available at  http://cern.ch/geant4/license .  These *
0009 // * include a list of copyright holders.                             *
0010 // *                                                                  *
0011 // * Neither the authors of this software system, nor their employing *
0012 // * institutes,nor the agencies providing financial support for this *
0013 // * work  make  any representation or  warranty, express or implied, *
0014 // * regarding  this  software system or assume any liability for its *
0015 // * use.  Please see the license in the file  LICENSE  and URL above *
0016 // * for the full disclaimer and the limitation of liability.         *
0017 // *                                                                  *
0018 // * This  code  implementation is the result of  the  scientific and *
0019 // * technical work of the GEANT4 collaboration.                      *
0020 // * By using,  copying,  modifying or  distributing the software (or *
0021 // * any work based  on the software)  you  agree  to acknowledge its *
0022 // * use  in  resulting  scientific  publications,  and indicate your *
0023 // * acceptance of all terms of the Geant4 Software license.          *
0024 // ********************************************************************
0025 //
0026 /*
0027  * ============================================================================
0028  *
0029  *       Filename:  CexmcSteppingAction.cc
0030  *
0031  *    Description:  stepping action
0032  *
0033  *        Version:  1.0
0034  *        Created:  27.10.2009 16:03:07
0035  *       Revision:  none
0036  *       Compiler:  gcc
0037  *
0038  *         Author:  Alexey Radkov (),
0039  *        Company:  PNPI
0040  *
0041  * ============================================================================
0042  */
0043
0044 #include <G4Step.hh>
0045 #include <G4Track.hh>
0046 #include <G4ParticleDefinition.hh>
0047 #include <G4VTouchable.hh>
0048 #include <G4NavigationHistory.hh>
0049 #include <G4AffineTransform.hh>
0050 #include <G4UnitsTable.hh>
0051 #include <G4RunManager.hh>
0052 #include "CexmcSteppingAction.hh"
0053 #include "CexmcPhysicsManager.hh"
0054 #include "CexmcSetup.hh"
0055 #include "CexmcIncidentParticleTrackInfo.hh"
0056 #include "CexmcCommon.hh"
0057
0058
0059 CexmcSteppingAction::CexmcSteppingAction(
0060                                     CexmcPhysicsManager *  physicsManager_ ) :
0061     physicsManager( physicsManager_ ), targetVolume( NULL )
0062 {
0063     G4RunManager *      runManager( G4RunManager::GetRunManager() );
0064     const CexmcSetup *  setup( static_cast< const CexmcSetup * >(
0065                                 runManager->GetUserDetectorConstruction() ) );
0066     targetVolume = setup->GetVolume( CexmcSetup::Target );
0067 }
0068
0069
0070 void  CexmcSteppingAction::UserSteppingAction( const G4Step *  step )
0071 {
0072     G4Track *         track( step->GetTrack() );
0073     CexmcTrackInfo *  trackInfo( static_cast< CexmcTrackInfo * >(
0074                                                 track->GetUserInformation() ) );
0075
0076     if ( ! trackInfo ||
0077          trackInfo->GetTypeInfo() != CexmcIncidentParticleTrackType )
0078         return;
0079
0080     CexmcIncidentParticleTrackInfo *  theTrackInfo(
0081                 static_cast< CexmcIncidentParticleTrackInfo * >( trackInfo ) );
0082
0083     G4StepPoint *         postStepPoint( step->GetPostStepPoint() );
0084     G4StepStatus          stepStatus( postStepPoint->GetStepStatus() );
0085     const G4VTouchable *  touchable( postStepPoint->GetTouchable() );
0086     G4VPhysicalVolume *   volume( touchable->GetVolume() );
0087
0088     if ( volume && volume->GetLogicalVolume() == targetVolume )
0089     {
0090         if ( ! theTrackInfo->IsStudiedProcessActivated() )
0091         {
0092             physicsManager->ResampleTrackLengthInTarget( track, postStepPoint );
0093             theTrackInfo->ActivateStudiedProcess();
0094         }
0095
0096         if ( stepStatus != fGeomBoundary )
0097         {
0098             if ( theTrackInfo->NeedsTrackLengthResampling() )
0099                 physicsManager->ResampleTrackLengthInTarget(
0100                                                         track, postStepPoint );
0101             else
0102                 theTrackInfo->AddTrackLengthInTarget( step->GetStepLength() );
0103         }
0104     }
0105
0106     G4StepPoint *  preStepPoint( step->GetPreStepPoint() );
0107     touchable = preStepPoint->GetTouchable();
0108     volume = touchable->GetVolume();
0109
0110     if ( volume && volume->GetLogicalVolume() == targetVolume )
0111     {
0112         if ( stepStatus == fGeomBoundary )
0113             theTrackInfo->ActivateStudiedProcess( false );
0114     }
0115 }
0116