Changeset 1315 for trunk/source/processes/cuts/src/G4ProductionCutsTable.cc
- Timestamp:
- Jun 18, 2010, 11:42:07 AM (14 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
trunk/source/processes/cuts/src/G4ProductionCutsTable.cc
r1228 r1315 25 25 // 26 26 // 27 // $Id: G4ProductionCutsTable.cc,v 1.2 5 2009/11/11 03:20:22kurasige Exp $28 // GEANT4 tag $Name: geant4-09-0 3$27 // $Id: G4ProductionCutsTable.cc,v 1.27 2010/01/22 11:10:26 kurasige Exp $ 28 // GEANT4 tag $Name: geant4-09-04-beta-cand-01 $ 29 29 // 30 30 // … … 252 252 if (verboseLevel>2) { 253 253 timer.Stop(); 254 std::cout << "G4ProductionCutsTable::UpdateCoupleTable "254 G4cout << "G4ProductionCutsTable::UpdateCoupleTable " 255 255 << " elapsed time for calculation of energy cuts " << G4endl; 256 std::cout << timer <<G4endl;256 G4cout << timer <<G4endl; 257 257 } 258 258 … … 291 291 292 292 // check range 293 if (range <=0.0) return -1.0; 293 if (range ==0.0) return 0.0; 294 if (range <0.0) return -1.0; 294 295 295 296 // check particle … … 922 923 G4bool fRatio = true; 923 924 for (size_t j=0; j< NumberOfG4CutIndex; j++) { 924 G4double ratio = cutValues[j]/aCut->GetProductionCut(j); 925 fRatio = fRatio && (0.999<ratio) && (ratio<1.001) ; 926 } 925 // check ratio only if values are not the same 926 if (cutValues[j] != aCut->GetProductionCut(j)) { 927 G4double ratio = cutValues[j]/aCut->GetProductionCut(j); 928 fRatio = fRatio && (0.999<ratio) && (ratio<1.001) ; 929 } 930 } 927 931 if (!fRatio) continue; 928 932 // MCC matched
Note: See TracChangeset
for help on using the changeset viewer.