Geant4 Cross Reference |
1 // 1 // 2 // ******************************************* 2 // ******************************************************************** 3 // * License and Disclaimer 3 // * License and Disclaimer * 4 // * 4 // * * 5 // * The Geant4 software is copyright of th 5 // * The Geant4 software is copyright of the Copyright Holders of * 6 // * the Geant4 Collaboration. It is provided 6 // * the Geant4 Collaboration. It is provided under the terms and * 7 // * conditions of the Geant4 Software License 7 // * conditions of the Geant4 Software License, included in the file * 8 // * LICENSE and available at http://cern.ch/ 8 // * LICENSE and available at http://cern.ch/geant4/license . These * 9 // * include a list of copyright holders. 9 // * include a list of copyright holders. * 10 // * 10 // * * 11 // * Neither the authors of this software syst 11 // * Neither the authors of this software system, nor their employing * 12 // * institutes,nor the agencies providing fin 12 // * institutes,nor the agencies providing financial support for this * 13 // * work make any representation or warran 13 // * work make any representation or warranty, express or implied, * 14 // * regarding this software system or assum 14 // * regarding this software system or assume any liability for its * 15 // * use. Please see the license in the file 15 // * use. Please see the license in the file LICENSE and URL above * 16 // * for the full disclaimer and the limitatio 16 // * for the full disclaimer and the limitation of liability. * 17 // * 17 // * * 18 // * This code implementation is the result 18 // * This code implementation is the result of the scientific and * 19 // * technical work of the GEANT4 collaboratio 19 // * technical work of the GEANT4 collaboration. * 20 // * By using, copying, modifying or distri 20 // * By using, copying, modifying or distributing the software (or * 21 // * any work based on the software) you ag 21 // * any work based on the software) you agree to acknowledge its * 22 // * use in resulting scientific publicati 22 // * use in resulting scientific publications, and indicate your * 23 // * acceptance of all terms of the Geant4 Sof 23 // * acceptance of all terms of the Geant4 Software license. * 24 // ******************************************* 24 // ******************************************************************** 25 // 25 // 26 /* 26 /* 27 * =========================================== 27 * ============================================================================ 28 * 28 * 29 * Filename: CexmcProductionModel.cc 29 * Filename: CexmcProductionModel.cc 30 * 30 * 31 * Description: interface to production mo 31 * Description: interface to production model 32 * 32 * 33 * Version: 1.0 33 * Version: 1.0 34 * Created: 03.11.2009 16:56:51 34 * Created: 03.11.2009 16:56:51 35 * Revision: none 35 * Revision: none 36 * Compiler: gcc 36 * Compiler: gcc 37 * 37 * 38 * Author: Alexey Radkov (), 38 * Author: Alexey Radkov (), 39 * Company: PNPI 39 * Company: PNPI 40 * 40 * 41 * =========================================== 41 * ============================================================================ 42 */ 42 */ 43 43 44 #include "CexmcRunManager.hh" 44 #include "CexmcRunManager.hh" 45 #include "CexmcProductionModel.hh" 45 #include "CexmcProductionModel.hh" 46 #include "CexmcProductionModelMessenger.hh" 46 #include "CexmcProductionModelMessenger.hh" 47 47 48 48 49 CexmcProductionModel::CexmcProductionModel( co << 49 CexmcProductionModel::CexmcProductionModel( const G4String & name, 50 G4 << 50 G4bool fermiMotionIsOn ) : 51 name( name_ ), fermiMotionIsOn( fermiMotio << 51 name( name ), fermiMotionIsOn( fermiMotionIsOn ), incidentParticle( NULL ), 52 incidentParticle( NULL ), nucleusParticle( << 52 nucleusParticle( NULL ), outputParticle( NULL ), 53 nucleusOutputParticle( NULL ), messenger( 53 nucleusOutputParticle( NULL ), messenger( NULL ) 54 { 54 { 55 angularRanges.push_back( CexmcAngularRange 55 angularRanges.push_back( CexmcAngularRange( 1.0, -1.0, 0 ) ); 56 messenger = new CexmcProductionModelMessen 56 messenger = new CexmcProductionModelMessenger( this ); 57 } 57 } 58 58 59 59 60 CexmcProductionModel::~CexmcProductionModel() 60 CexmcProductionModel::~CexmcProductionModel() 61 { 61 { 62 delete messenger; 62 delete messenger; 63 } 63 } 64 64 65 65 66 void CexmcProductionModel::SetAngularRange( G 66 void CexmcProductionModel::SetAngularRange( G4double top, G4double bottom, 67 G 67 G4int nmbOfDivs ) 68 { 68 { 69 if ( ! IsValidCandidateForAngularRange( to 69 if ( ! IsValidCandidateForAngularRange( top, bottom, nmbOfDivs ) ) 70 throw CexmcException( CexmcInvalidAngu 70 throw CexmcException( CexmcInvalidAngularRange ); 71 71 72 if ( ! IsGoodCandidateForAngularRange( top 72 if ( ! IsGoodCandidateForAngularRange( top, bottom ) ) 73 throw CexmcException( CexmcBadAngularR 73 throw CexmcException( CexmcBadAngularRange ); 74 74 75 angularRanges.clear(); 75 angularRanges.clear(); 76 G4double curBottom( top ); 76 G4double curBottom( top ); 77 for ( int i( 0 ); i < nmbOfDivs; ++i ) 77 for ( int i( 0 ); i < nmbOfDivs; ++i ) 78 { 78 { 79 G4double binWidth( ( top - bottom ) / 79 G4double binWidth( ( top - bottom ) / nmbOfDivs ); 80 G4double curTop( curBottom ); 80 G4double curTop( curBottom ); 81 curBottom -= binWidth; 81 curBottom -= binWidth; 82 angularRanges.push_back( CexmcAngularR 82 angularRanges.push_back( CexmcAngularRange( curTop, curBottom, i ) ); 83 } 83 } 84 #ifdef CEXMC_USE_ROOT 84 #ifdef CEXMC_USE_ROOT 85 CexmcHistoManager::Instance()->SetupARHist 85 CexmcHistoManager::Instance()->SetupARHistos( angularRanges ); 86 #endif 86 #endif 87 } 87 } 88 88 89 89 90 void CexmcProductionModel::AddAngularRange( G 90 void CexmcProductionModel::AddAngularRange( G4double top, G4double bottom, 91 G 91 G4int nmbOfDivs ) 92 { 92 { 93 if ( ! IsValidCandidateForAngularRange( to 93 if ( ! IsValidCandidateForAngularRange( top, bottom, nmbOfDivs ) ) 94 throw CexmcException( CexmcInvalidAngu 94 throw CexmcException( CexmcInvalidAngularRange ); 95 95 96 if ( ! IsGoodCandidateForAngularRange( top 96 if ( ! IsGoodCandidateForAngularRange( top, bottom ) ) 97 throw CexmcException( CexmcBadAngularR 97 throw CexmcException( CexmcBadAngularRange ); 98 98 99 G4int curIndex( angularRanges.size() ); 99 G4int curIndex( angularRanges.size() ); 100 G4double curBottom( top ); 100 G4double curBottom( top ); 101 for ( int i( 0 ); i < nmbOfDivs; ++i ) 101 for ( int i( 0 ); i < nmbOfDivs; ++i ) 102 { 102 { 103 G4double binWidth( ( top - bottom ) / 103 G4double binWidth( ( top - bottom ) / nmbOfDivs ); 104 G4double curTop( curBottom ); 104 G4double curTop( curBottom ); 105 curBottom -= binWidth; 105 curBottom -= binWidth; 106 CexmcAngularRange aRange( curTop, cur 106 CexmcAngularRange aRange( curTop, curBottom, curIndex + i ); 107 angularRanges.push_back( aRange ); 107 angularRanges.push_back( aRange ); 108 #ifdef CEXMC_USE_ROOT 108 #ifdef CEXMC_USE_ROOT 109 CexmcHistoManager::Instance()->AddARHi 109 CexmcHistoManager::Instance()->AddARHistos( aRange ); 110 #endif 110 #endif 111 } 111 } 112 } 112 } 113 113 114 114 115 void CexmcProductionModel::SetTriggeredAngula 115 void CexmcProductionModel::SetTriggeredAngularRanges( G4double opCosThetaSCM ) 116 { 116 { 117 triggeredAngularRanges.clear(); 117 triggeredAngularRanges.clear(); 118 118 119 for ( CexmcAngularRangeList::iterator k( 119 for ( CexmcAngularRangeList::iterator k( angularRanges.begin() ); 120 k 120 k != angularRanges.end(); ++k ) 121 { 121 { 122 if ( opCosThetaSCM <= k->top && opCosT 122 if ( opCosThetaSCM <= k->top && opCosThetaSCM > k->bottom ) 123 triggeredAngularRanges.push_back( 123 triggeredAngularRanges.push_back( CexmcAngularRange( 124 k- 124 k->top, k->bottom, k->index ) ); 125 } 125 } 126 } 126 } 127 127 128 128 129 void CexmcProductionModel::FermiMotionStatusC 129 void CexmcProductionModel::FermiMotionStatusChangeHook( void ) 130 { 130 { 131 } 131 } 132 132 133 133 134 G4bool CexmcProductionModel::IsGoodCandidateF 134 G4bool CexmcProductionModel::IsGoodCandidateForAngularRange( G4double top, 135 135 G4double bottom ) const 136 { 136 { 137 CexmcRunManager * runManager( static_cast 137 CexmcRunManager * runManager( static_cast< CexmcRunManager * >( 138 G4R 138 G4RunManager::GetRunManager() ) ); 139 139 140 if ( ! runManager->ProjectIsRead() ) 140 if ( ! runManager->ProjectIsRead() ) 141 return true; 141 return true; 142 142 143 CexmcAngularRangeList normalizedARanges; 143 CexmcAngularRangeList normalizedARanges; 144 GetNormalizedAngularRange( angularRangesRe 144 GetNormalizedAngularRange( angularRangesRef, normalizedARanges ); 145 145 146 for ( CexmcAngularRangeList::iterator k( 146 for ( CexmcAngularRangeList::iterator k( normalizedARanges.begin() ); 147 k 147 k != normalizedARanges.end(); ++k ) 148 { 148 { 149 if ( top <= k->top && bottom >= k->bot 149 if ( top <= k->top && bottom >= k->bottom ) 150 return true; 150 return true; 151 } 151 } 152 152 153 return false; 153 return false; 154 } 154 } 155 155 156 156