Geant4 Cross Reference

Cross-Referencing   Geant4
Geant4/examples/advanced/ChargeExchangeMC/src/CexmcProductionModel.cc

Version: [ ReleaseNotes ] [ 1.0 ] [ 1.1 ] [ 2.0 ] [ 3.0 ] [ 3.1 ] [ 3.2 ] [ 4.0 ] [ 4.0.p1 ] [ 4.0.p2 ] [ 4.1 ] [ 4.1.p1 ] [ 5.0 ] [ 5.0.p1 ] [ 5.1 ] [ 5.1.p1 ] [ 5.2 ] [ 5.2.p1 ] [ 5.2.p2 ] [ 6.0 ] [ 6.0.p1 ] [ 6.1 ] [ 6.2 ] [ 6.2.p1 ] [ 6.2.p2 ] [ 7.0 ] [ 7.0.p1 ] [ 7.1 ] [ 7.1.p1 ] [ 8.0 ] [ 8.0.p1 ] [ 8.1 ] [ 8.1.p1 ] [ 8.1.p2 ] [ 8.2 ] [ 8.2.p1 ] [ 8.3 ] [ 8.3.p1 ] [ 8.3.p2 ] [ 9.0 ] [ 9.0.p1 ] [ 9.0.p2 ] [ 9.1 ] [ 9.1.p1 ] [ 9.1.p2 ] [ 9.1.p3 ] [ 9.2 ] [ 9.2.p1 ] [ 9.2.p2 ] [ 9.2.p3 ] [ 9.2.p4 ] [ 9.3 ] [ 9.3.p1 ] [ 9.3.p2 ] [ 9.4 ] [ 9.4.p1 ] [ 9.4.p2 ] [ 9.4.p3 ] [ 9.4.p4 ] [ 9.5 ] [ 9.5.p1 ] [ 9.5.p2 ] [ 9.6 ] [ 9.6.p1 ] [ 9.6.p2 ] [ 9.6.p3 ] [ 9.6.p4 ] [ 10.0 ] [ 10.0.p1 ] [ 10.0.p2 ] [ 10.0.p3 ] [ 10.0.p4 ] [ 10.1 ] [ 10.1.p1 ] [ 10.1.p2 ] [ 10.1.p3 ] [ 10.2 ] [ 10.2.p1 ] [ 10.2.p2 ] [ 10.2.p3 ] [ 10.3 ] [ 10.3.p1 ] [ 10.3.p2 ] [ 10.3.p3 ] [ 10.4 ] [ 10.4.p1 ] [ 10.4.p2 ] [ 10.4.p3 ] [ 10.5 ] [ 10.5.p1 ] [ 10.6 ] [ 10.6.p1 ] [ 10.6.p2 ] [ 10.6.p3 ] [ 10.7 ] [ 10.7.p1 ] [ 10.7.p2 ] [ 10.7.p3 ] [ 10.7.p4 ] [ 11.0 ] [ 11.0.p1 ] [ 11.0.p2 ] [ 11.0.p3, ] [ 11.0.p4 ] [ 11.1 ] [ 11.1.1 ] [ 11.1.2 ] [ 11.1.3 ] [ 11.2 ] [ 11.2.1 ] [ 11.2.2 ] [ 11.3.0 ]

  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:  CexmcProductionModel.cc
 30  *
 31  *    Description:  interface to production model
 32  *
 33  *        Version:  1.0
 34  *        Created:  03.11.2009 16:56:51
 35  *       Revision:  none
 36  *       Compiler:  gcc
 37  *
 38  *         Author:  Alexey Radkov (), 
 39  *        Company:  PNPI
 40  *
 41  * ============================================================================
 42  */
 43 
 44 #include "CexmcRunManager.hh"
 45 #include "CexmcProductionModel.hh"
 46 #include "CexmcProductionModelMessenger.hh"
 47 
 48 
 49 CexmcProductionModel::CexmcProductionModel( const G4String &  name_,
 50                                             G4bool  fermiMotionIsOn_ ) :
 51     name( name_ ), fermiMotionIsOn( fermiMotionIsOn_ ),
 52     incidentParticle( NULL ), nucleusParticle( NULL ), outputParticle( NULL ),
 53     nucleusOutputParticle( NULL ), messenger( NULL )
 54 {
 55     angularRanges.push_back( CexmcAngularRange( 1.0, -1.0, 0 ) );
 56     messenger = new CexmcProductionModelMessenger( this );
 57 }
 58 
 59 
 60 CexmcProductionModel::~CexmcProductionModel()
 61 {
 62     delete messenger;
 63 }
 64 
 65 
 66 void  CexmcProductionModel::SetAngularRange( G4double  top, G4double  bottom,
 67                                              G4int  nmbOfDivs )
 68 {
 69     if ( ! IsValidCandidateForAngularRange( top, bottom, nmbOfDivs ) )
 70         throw CexmcException( CexmcInvalidAngularRange );
 71 
 72     if ( ! IsGoodCandidateForAngularRange( top, bottom ) )
 73         throw CexmcException( CexmcBadAngularRange );
 74 
 75     angularRanges.clear();
 76     G4double  curBottom( top );
 77     for ( int  i( 0 ); i < nmbOfDivs; ++i )
 78     {
 79         G4double  binWidth( ( top - bottom ) / nmbOfDivs );
 80         G4double  curTop( curBottom );
 81         curBottom -= binWidth;
 82         angularRanges.push_back( CexmcAngularRange( curTop, curBottom, i ) );
 83     }
 84 #ifdef CEXMC_USE_ROOT
 85     CexmcHistoManager::Instance()->SetupARHistos( angularRanges );
 86 #endif
 87 }
 88 
 89 
 90 void  CexmcProductionModel::AddAngularRange( G4double  top, G4double  bottom,
 91                                              G4int  nmbOfDivs )
 92 {
 93     if ( ! IsValidCandidateForAngularRange( top, bottom, nmbOfDivs ) )
 94         throw CexmcException( CexmcInvalidAngularRange );
 95 
 96     if ( ! IsGoodCandidateForAngularRange( top, bottom ) )
 97         throw CexmcException( CexmcBadAngularRange );
 98 
 99     G4int  curIndex( angularRanges.size() );
100     G4double  curBottom( top );
101     for ( int  i( 0 ); i < nmbOfDivs; ++i )
102     {
103         G4double  binWidth( ( top - bottom ) / nmbOfDivs );
104         G4double  curTop( curBottom );
105         curBottom -= binWidth;
106         CexmcAngularRange  aRange( curTop, curBottom, curIndex + i );
107         angularRanges.push_back( aRange );
108 #ifdef CEXMC_USE_ROOT
109         CexmcHistoManager::Instance()->AddARHistos( aRange );
110 #endif
111     }
112 }
113 
114 
115 void  CexmcProductionModel::SetTriggeredAngularRanges( G4double  opCosThetaSCM )
116 {
117     triggeredAngularRanges.clear();
118 
119     for ( CexmcAngularRangeList::iterator  k( angularRanges.begin() );
120                                             k != angularRanges.end(); ++k )
121     {
122         if ( opCosThetaSCM <= k->top && opCosThetaSCM > k->bottom )
123             triggeredAngularRanges.push_back( CexmcAngularRange(
124                                             k->top, k->bottom, k->index ) );
125     }
126 }
127 
128 
129 void  CexmcProductionModel::FermiMotionStatusChangeHook( void )
130 {
131 }
132 
133 
134 G4bool  CexmcProductionModel::IsGoodCandidateForAngularRange( G4double  top,
135                                                       G4double  bottom ) const
136 {
137     CexmcRunManager *  runManager( static_cast< CexmcRunManager * >(
138                                            G4RunManager::GetRunManager() ) );
139 
140     if ( ! runManager->ProjectIsRead() )
141         return true;
142 
143     CexmcAngularRangeList  normalizedARanges;
144     GetNormalizedAngularRange( angularRangesRef, normalizedARanges );
145 
146     for ( CexmcAngularRangeList::iterator  k( normalizedARanges.begin() );
147                                             k != normalizedARanges.end(); ++k )
148     {
149         if ( top <= k->top && bottom >= k->bottom )
150             return true;
151     }
152 
153     return false;
154 }
155 
156