Ignore:
Timestamp:
Apr 6, 2009, 12:30:29 PM (15 years ago)
Author:
garnier
Message:

update processes

File:
1 edited

Legend:

Unmodified
Added
Removed
  • trunk/source/processes/hadronic/cross_sections/src/G4PiNuclearCrossSection.cc

    r819 r962  
    2727 #include "G4HadronicException.hh"
    2828 #include "G4HadTmpUtil.hh"
    29  #include "G4ping.hh"
     29// #include "G4ping.hh"
    3030
    3131 // by J.P Wellisch, Sun Sep 15 2002.
     
    350350{
    351351  // precondition
    352   G4ping debug("debug_PiNuclearCrossSection");
     352  //  G4ping debug("debug_PiNuclearCrossSection");
    353353  using namespace std;
    354354  G4bool ok = false;
     
    368368  G4double result = 0;
    369369  G4int Z=G4lrint(ZZ);
    370   debug.push_back(Z);
     370  //  debug.push_back(Z);
    371371  size_t it = 0;
    372372
    373373  while(it<theZ.size() && Z>theZ[it]) it++;
    374374
    375   debug.push_back(theZ[it]);
    376   debug.push_back(kineticEnergy);
     375  //  debug.push_back(theZ[it]);
     376  //  debug.push_back(kineticEnergy);
    377377
    378378  if(Z > theZ[it])
     
    390390      fTotalXsc = thePimData[it]->TotalXSection(kineticEnergy);
    391391
    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);
    395395    }
    396396    else
     
    406406      fTotalXsc = Interpolate(Z1, Z2, Z, xt1, xt2);
    407407
    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);
    417417    }
    418418  }
     
    430430      fTotalXsc = theData->operator[](it)->TotalXSection(kineticEnergy);
    431431
    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);
    435435    }
    436436    else
     
    456456      fTotalXsc = Interpolate(Z1, Z2, Z, xt1, xt2);
    457457
    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);
    467467    }
    468468  }
    469   debug.dump();
     469  //  debug.dump();
    470470
    471471  fElasticXsc = fTotalXsc - result;
Note: See TracChangeset for help on using the changeset viewer.