Changeset 1315 for trunk/source/particles/management/src/G4DecayProducts.cc
- Timestamp:
- Jun 18, 2010, 11:42:07 AM (14 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
trunk/source/particles/management/src/G4DecayProducts.cc
r1196 r1315 25 25 // 26 26 // 27 // $Id: G4DecayProducts.cc,v 1.1 7 2009/04/02 02:24:53kurasige Exp $28 // GEANT4 tag $Name: geant4-09-0 3-cand-01 $27 // $Id: G4DecayProducts.cc,v 1.18 2010/05/20 01:01:07 kurasige Exp $ 28 // GEANT4 tag $Name: geant4-09-04-beta-cand-01 $ 29 29 // 30 30 // … … 141 141 } else { 142 142 #ifdef G4VERBOSE 143 G4c out << "G4DecayProducts::PushProducts ";144 G4cout<< " exceeds MaxNumberOfProducts(="145 << G4int(MaxNumberOfProducts) << ")" ;146 G4cout<< G4endl;143 G4cerr << "G4DecayProducts::PushProducts " 144 << " exceeds MaxNumberOfProducts(=" 145 << G4int(MaxNumberOfProducts) << ")" 146 << G4endl; 147 147 #endif 148 148 } … … 242 242 if ( (parent_momentum.mag() >0.0) && (std::fabs(direction.mag()-1.0) >1.0e-6 ) ) { 243 243 #ifdef G4VERBOSE 244 G4cout << " Momentum Direction Vector of Parent is not normalized "; 245 G4cout << " (=" << direction.mag() << ")" << G4endl; 244 G4cerr << "G4DecayProducts::IsChecked():: " 245 << " Momentum Direction Vector of Parent is not normalized " 246 << " (=" << direction.mag() << ")" << G4endl; 246 247 #endif 247 248 returnValue = false; … … 263 264 if ( (momentum.mag()>0.0) && (std::fabs(direction.mag()-1.0) > 1.0e-6)) { 264 265 #ifdef G4VERBOSE 265 G4cout << " Momentum Direction Vector of Daughter [" << index; 266 G4cout << "] is not normalized (=" << direction.mag() << ")" << G4endl; 266 G4cerr << "G4DecayProducts::IsChecked():: " 267 << " Momentum Direction Vector of Daughter [" << index 268 << "] is not normalized (=" << direction.mag() << ")" << G4endl; 267 269 #endif 268 270 returnValue = false; … … 272 274 if (energy - mass < DBL_MIN ) { 273 275 #ifdef G4VERBOSE 274 G4cout << "Daughter [" << index << "] has no kinetic energy "<< G4endl; 276 G4cerr << "G4DecayProducts::IsChecked():: " 277 << " Daughter [" << index << "] has no kinetic energy "<< G4endl; 275 278 #endif 276 279 returnValue = false; … … 282 285 if ( (std::fabs(total_energy) >1.0e-9*MeV) || (total_momentum.mag() >1.0e-9*MeV ) ){ 283 286 #ifdef G4VERBOSE 284 G4cout << " Energy/Momentum is not conserved "<< G4endl; 285 G4cout << " difference between parent energy and sum of dughters' energy : " 287 G4cerr << "G4DecayProducts::IsChecked():: " 288 << " Energy/Momentum is not conserved "<< G4endl; 289 G4cerr << " difference between parent energy and sum of dughters' energy : " 286 290 << total_energy /MeV << "[MeV] " << G4endl; 287 G4c out<< " difference between parent momentum and sum of dughters' momentum : "291 G4cerr << " difference between parent momentum and sum of dughters' momentum : " 288 292 << " x:" << total_momentum.getX()/MeV 289 293 << " y:" << total_momentum.getY()/MeV
Note: See TracChangeset
for help on using the changeset viewer.