Changeset 1196 for trunk/source/processes/hadronic/models/qmd
- Timestamp:
- Nov 25, 2009, 5:13:58 PM (15 years ago)
- Location:
- trunk/source/processes/hadronic/models/qmd
- Files:
-
- 3 edited
Legend:
- Unmodified
- Added
- Removed
-
trunk/source/processes/hadronic/models/qmd/History
r1055 r1196 14 14 --------------------------------------------------------------- 15 15 16 19-November-2009 Tatsumi Koi (hadr-qmd-V09-02-02) 17 - Fix bug for incidence of neutral particles; 18 Should be OK for neucleon-nucleus interaction 19 src/G4QMDReaction.cc 20 src/G4QMDGroundStateNucleus.cc 21 16 22 31-March-2009 Tatsumi Koi (hadr-qmd-V09-02-01) 17 23 - Fix bug in gamma(mass zero) partiticpants 18 24 src/G4QMDCollision.cc 19 - Chnage me ber object to pointer25 - Chnage member object to pointer 20 26 include/G4QMDReaction.hh 21 27 src/G4QMDReaction.cc -
trunk/source/processes/hadronic/models/qmd/src/G4QMDGroundStateNucleus.cc
r962 r1196 48 48 //std::cout << " G4QMDGroundStateNucleus( G4int z , G4int a ) Begin " << z << " " << a << std::endl; 49 49 50 // Hydrogen Case 51 if ( z == 1 && a == 1 ) 50 if ( z == 1 && a == 1 ) // Hydrogen Case or proton primary 52 51 { 53 52 SetParticipant( new G4QMDParticipant( G4Proton::Proton() , G4ThreeVector( 0.0 ) , G4ThreeVector( 0.0 ) ) ); 53 return; 54 } 55 else if ( z == 0 && a == 1 ) // Neutron primary 56 { 57 SetParticipant( new G4QMDParticipant( G4Neutron::Neutron() , G4ThreeVector( 0.0 ) , G4ThreeVector( 0.0 ) ) ); 54 58 return; 55 59 } -
trunk/source/processes/hadronic/models/qmd/src/G4QMDReaction.cc
r1055 r1196 33 33 // 090122 Be8 -> Alpha + Alpha 34 34 // 090331 Change member shenXS and genspaXS object to pointer 35 // 091119 Fix for incidence of neutral particles 35 36 // 36 37 #include "G4QMDReaction.hh" … … 636 637 G4double pccf = std::sqrt( pcca ); 637 638 638 G4double aas = 2.0 * eccm * b / double ( zp * zt ) / ccoul; 639 G4double bbs = 1.0 / std::sqrt ( 1.0 + aas*aas ); 640 G4double aas1 = ( 1.0 + aas * b / rmax ) * bbs; 639 //Fix for neutral particles 640 G4double aas1 = 0.0; 641 G4double bbs = 0.0; 642 643 if ( zp != 0 ) 644 { 645 G4double aas = 2.0 * eccm * b / double ( zp * zt ) / ccoul; 646 bbs = 1.0 / std::sqrt ( 1.0 + aas*aas ); 647 aas1 = ( 1.0 + aas * b / rmax ) * bbs; 648 } 641 649 642 650 G4double cost = 0.0;
Note: See TracChangeset
for help on using the changeset viewer.