- Timestamp:
- Dec 22, 2010, 3:52:27 PM (14 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
trunk/source/processes/hadronic/models/util/src/G4KineticTrack.cc
r1196 r1347 493 493 G4cout << "DECAY Actual Mass " << theActualMass << G4endl; 494 494 */ 495 G4ParticleDefinition* thisDefinition = this->GetDefinition(); 496 if(!thisDefinition) 497 { 498 G4cerr << "Error condition encountered in G4KineticTrack::Decay()"<<G4endl; 499 G4cerr << " track has no particle definition associated."<<G4endl; 500 return 0; 501 } 502 G4DecayTable* theDecayTable = thisDefinition->GetDecayTable(); 503 if(!theDecayTable) 504 { 505 G4cerr << "Error condition encountered in G4KineticTrack::Decay()"<<G4endl; 506 G4cerr << " particle definiton has no decay table associated."<<G4endl; 507 G4cerr << " particle was "<<thisDefinition->GetParticleName()<<G4endl; 508 return 0; 509 } 495 510 496 511 G4int chargeBalance = G4lrint(theDefinition->GetPDGCharge() ); … … 512 527 // cout << "DECAY Total Width " << theTotalActualWidth << G4endl; 513 528 G4double r = theTotalActualWidth * G4UniformRand(); 514 G4ParticleDefinition* thisDefinition = this->GetDefinition(); 515 if(!thisDefinition) 516 { 517 G4cerr << "Error condition encountered in G4KineticTrack::Decay()"<<G4endl; 518 G4cerr << " track has no particle definition associated."<<G4endl; 519 return 0; 520 } 521 G4DecayTable* theDecayTable = thisDefinition->GetDecayTable(); 522 if(!theDecayTable) 523 { 524 G4cerr << "Error condition encountered in G4KineticTrack::Decay()"<<G4endl; 525 G4cerr << " particle definiton has no decay table associated."<<G4endl; 526 G4cerr << " particle was "<<thisDefinition->GetParticleName()<<G4endl; 527 return 0; 528 } 529 G4VDecayChannel* theDecayChannel=NULL; 529 G4VDecayChannel* theDecayChannel(0); 530 530 for (index = nChannels - 1; index >= 0; index--) 531 531 { … … 538 538 } 539 539 } 540 540 541 delete [] theCumActualWidth; 542 541 543 if(!theDecayChannel) 542 544 { … … 669 671 } 670 672 delete theDecayProducts; 671 delete [] theCumActualWidth;672 673 if(getenv("DecayEnergyBalanceCheck")) 673 674 std::cout << "DEBUGGING energy balance in cms and lab, charge baryon balance : "
Note: See TracChangeset
for help on using the changeset viewer.