Geant4 Cross Reference |
1 // 2 // ******************************************************************** 3 // * License and Disclaimer * 4 // * * 5 // * The Geant4 software is copyright of the Copyright Holders of * 6 // * the Geant4 Collaboration. It is provided under the terms and * 7 // * conditions of the Geant4 Software License, included in the file * 8 // * LICENSE and available at http://cern.ch/geant4/license . These * 9 // * include a list of copyright holders. * 10 // * * 11 // * Neither the authors of this software system, nor their employing * 12 // * institutes,nor the agencies providing financial support for this * 13 // * work make any representation or warranty, express or implied, * 14 // * regarding this software system or assume any liability for its * 15 // * use. Please see the license in the file LICENSE and URL above * 16 // * for the full disclaimer and the limitation of liability. * 17 // * * 18 // * This code implementation is the result of the scientific and * 19 // * technical work of the GEANT4 collaboration. * 20 // * By using, copying, modifying or distributing the software (or * 21 // * any work based on the software) you agree to acknowledge its * 22 // * use in resulting scientific publications, and indicate your * 23 // * acceptance of all terms of the Geant4 Software license. * 24 // ******************************************************************** 25 // 26 /* 27 * ============================================================================= 28 * 29 * Filename: CexmcStudiedPhysics.hh 30 * 31 * Description: studied physics in the target 32 * 33 * Version: 1.0 34 * Created: 18.10.2009 16:10:52 35 * Revision: none 36 * Compiler: gcc 37 * 38 * Author: Alexey Radkov (), 39 * Company: PNPI 40 * 41 * ============================================================================= 42 */ 43 44 #ifndef CEXMC_STUDIED_PHYSICS_HH 45 #define CEXMC_STUDIED_PHYSICS_HH 46 47 #include <G4VPhysicsConstructor.hh> 48 #include <G4ProcessManager.hh> 49 #include <G4ParticleDefinition.hh> 50 #include "CexmcFakeCrossSectionData.hh" 51 #include "CexmcStudiedProcess.hh" 52 #include "CexmcPhysicsManager.hh" 53 #include "CexmcProductionModel.hh" 54 55 class G4VProcess; 56 57 58 template < typename Process > 59 class CexmcStudiedPhysics : public G4VPhysicsConstructor 60 { 61 public: 62 explicit CexmcStudiedPhysics( CexmcPhysicsManager * physicsManager ); 63 64 virtual ~CexmcStudiedPhysics(); 65 66 public: 67 void ConstructParticle( void ); 68 69 void ConstructProcess( void ); 70 71 public: 72 CexmcProductionModel * GetProductionModel( void ); 73 74 protected: 75 virtual void ApplyInteractionModel( G4VProcess * process ); 76 77 protected: 78 CexmcPhysicsManager * physicsManager; 79 80 CexmcProductionModel * productionModel; 81 82 private: 83 G4bool wasActivated; 84 }; 85 86 87 template < typename Process > 88 CexmcStudiedPhysics< Process >::CexmcStudiedPhysics( 89 CexmcPhysicsManager * physicsManager_ ) : 90 G4VPhysicsConstructor( "studiedPhysics" ), 91 physicsManager( physicsManager_ ), productionModel( NULL ), 92 wasActivated( false ) 93 { 94 } 95 96 97 template < typename Process > 98 CexmcStudiedPhysics< Process >::~CexmcStudiedPhysics() 99 { 100 } 101 102 103 template < typename Process > 104 void CexmcStudiedPhysics< Process >::ConstructParticle( void ) 105 { 106 if ( productionModel ) 107 productionModel->GetIncidentParticle(); 108 } 109 110 111 template < typename Process > 112 void CexmcStudiedPhysics< Process >::ConstructProcess( void ) 113 { 114 if ( wasActivated ) 115 return; 116 117 wasActivated = true; 118 119 Process * process( new Process ); 120 121 process->AddDataSet( new CexmcFakeCrossSectionData ); 122 123 CexmcStudiedProcess * studiedProcess( new CexmcStudiedProcess( 124 physicsManager, process->GetProcessType() ) ); 125 126 ApplyInteractionModel( process ); 127 128 studiedProcess->RegisterProcess( process ); 129 130 G4ParticleDefinition * particle( NULL ); 131 132 if ( productionModel ) 133 particle = productionModel->GetIncidentParticle(); 134 135 if ( particle ) 136 { 137 G4ProcessManager * processManager( particle->GetProcessManager() ); 138 processManager->AddDiscreteProcess( studiedProcess ); 139 } 140 } 141 142 143 template < typename Process > 144 CexmcProductionModel * CexmcStudiedPhysics< Process >::GetProductionModel( 145 void ) 146 { 147 return productionModel; 148 } 149 150 151 template < typename Process > 152 void CexmcStudiedPhysics< Process >::ApplyInteractionModel( G4VProcess * ) 153 { 154 } 155 156 157 #endif 158 159