- Timestamp:
- Apr 6, 2009, 12:30:29 PM (15 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
trunk/source/processes/hadronic/cross_sections/src/G4PiNuclearCrossSection.cc
r819 r962 27 27 #include "G4HadronicException.hh" 28 28 #include "G4HadTmpUtil.hh" 29 #include "G4ping.hh"29 // #include "G4ping.hh" 30 30 31 31 // by J.P Wellisch, Sun Sep 15 2002. … … 350 350 { 351 351 // precondition 352 G4ping debug("debug_PiNuclearCrossSection");352 // G4ping debug("debug_PiNuclearCrossSection"); 353 353 using namespace std; 354 354 G4bool ok = false; … … 368 368 G4double result = 0; 369 369 G4int Z=G4lrint(ZZ); 370 debug.push_back(Z);370 // debug.push_back(Z); 371 371 size_t it = 0; 372 372 373 373 while(it<theZ.size() && Z>theZ[it]) it++; 374 374 375 debug.push_back(theZ[it]);376 debug.push_back(kineticEnergy);375 // debug.push_back(theZ[it]); 376 // debug.push_back(kineticEnergy); 377 377 378 378 if(Z > theZ[it]) … … 390 390 fTotalXsc = thePimData[it]->TotalXSection(kineticEnergy); 391 391 392 debug.push_back("D1 ");393 debug.push_back(result);394 debug.push_back(fTotalXsc);392 // debug.push_back("D1 "); 393 // debug.push_back(result); 394 // debug.push_back(fTotalXsc); 395 395 } 396 396 else … … 406 406 fTotalXsc = Interpolate(Z1, Z2, Z, xt1, xt2); 407 407 408 debug.push_back("D2 ");409 debug.push_back(x1);410 debug.push_back(x2);411 debug.push_back(xt1);412 debug.push_back(xt2);413 debug.push_back(Z1);414 debug.push_back(Z2);415 debug.push_back(result);416 debug.push_back(fTotalXsc);408 // debug.push_back("D2 "); 409 // debug.push_back(x1); 410 // debug.push_back(x2); 411 // debug.push_back(xt1); 412 // debug.push_back(xt2); 413 // debug.push_back(Z1); 414 // debug.push_back(Z2); 415 // debug.push_back(result); 416 // debug.push_back(fTotalXsc); 417 417 } 418 418 } … … 430 430 fTotalXsc = theData->operator[](it)->TotalXSection(kineticEnergy); 431 431 432 debug.push_back("D3 ");433 debug.push_back(result);434 debug.push_back(fTotalXsc);432 // debug.push_back("D3 "); 433 // debug.push_back(result); 434 // debug.push_back(fTotalXsc); 435 435 } 436 436 else … … 456 456 fTotalXsc = Interpolate(Z1, Z2, Z, xt1, xt2); 457 457 458 debug.push_back("D4 ");459 debug.push_back(x1);460 debug.push_back(xt1);461 debug.push_back(x2);462 debug.push_back(xt2);463 debug.push_back(Z1);464 debug.push_back(Z2);465 debug.push_back(result);466 debug.push_back(fTotalXsc);458 // debug.push_back("D4 "); 459 // debug.push_back(x1); 460 // debug.push_back(xt1); 461 // debug.push_back(x2); 462 // debug.push_back(xt2); 463 // debug.push_back(Z1); 464 // debug.push_back(Z2); 465 // debug.push_back(result); 466 // debug.push_back(fTotalXsc); 467 467 } 468 468 } 469 debug.dump();469 // debug.dump(); 470 470 471 471 fElasticXsc = fTotalXsc - result;
Note: See TracChangeset
for help on using the changeset viewer.