Changeset 1347 for trunk/source/processes/electromagnetic/xrays/test
- Timestamp:
- Dec 22, 2010, 3:52:27 PM (14 years ago)
- Location:
- trunk/source/processes/electromagnetic/xrays/test
- Files:
-
- 2 edited
Legend:
- Unmodified
- Added
- Removed
-
trunk/source/processes/electromagnetic/xrays/test/testG4ForwardXrayTR.cc
r1315 r1347 26 26 // 27 27 // $Id: testG4ForwardXrayTR.cc,v 1.7 2006/06/29 19:56:33 gunter Exp $ 28 // GEANT4 tag $Name: geant4-09-04- beta-cand-01$28 // GEANT4 tag $Name: geant4-09-04-ref-00 $ 29 29 // 30 30 -
trunk/source/processes/electromagnetic/xrays/test/testG4HeatedKleinNishina.cc
r1199 r1347 61 61 #include "G4HeatedKleinNishinaCompton.hh" 62 62 #include "G4PEEffectModel.hh" 63 #include "G4PairProductionRelModel.hh" 64 #include "G4BetheHeitlerModel.hh" 63 65 64 66 … … 102 104 G4int choice; 103 105 // G4cin >> choice; 104 choice = 54;106 choice = 6; 105 107 106 108 … … 298 300 kinEnergy = 25.0*GeV; 299 301 300 //G4VEmModel* comp = new G4KleinNishinaCompton();301 G4VEmModel* comp = new G4HeatedKleinNishinaCompton();302 G4VEmModel* comp = new G4KleinNishinaCompton(); 303 // G4VEmModel* comp = new G4HeatedKleinNishinaCompton(); 302 304 G4VEmModel* photo = new G4PEEffectModel(); 303 305 G4VEmModel* pair = new G4PairProductionRelModel(theParticleDefinition,"pp"); 306 G4VEmModel* bhpair = new G4BetheHeitlerModel(theParticleDefinition,"bhpp"); 304 307 305 308 G4DynamicParticle* theDynamicParticle = new G4DynamicParticle(theParticleDefinition, … … 326 329 327 330 331 // energyMscXR = 2*MeV; 332 //xsc = pair->ComputeCrossSectionPerAtom(theParticleDefinition, energyMscXR, G4double(Z ) ); 333 328 334 329 335 iMax = 130; 330 336 331 337 writef<<iMax<<G4endl; 338 339 340 332 341 333 342 for( i = 0; i < iMax; i++ ) 334 343 { 335 344 energyMscXR = std::exp(i*0.2)*0.001*keV; 336 //xsc = comp->ComputeCrossSectionPerAtom(theParticleDefinition,energyMscXR,G4double(Z ) ); 337 xsc = photo->ComputeCrossSectionPerAtom(theParticleDefinition,energyMscXR,G4double(Z ) ); 345 // xsc = comp->ComputeCrossSectionPerAtom(theParticleDefinition,energyMscXR,G4double(Z ) ); 346 // xsc = photo->ComputeCrossSectionPerAtom(theParticleDefinition, energyMscXR, G4double(Z ) ); 347 // xsc = pair->ComputeCrossSectionPerAtom(theParticleDefinition, energyMscXR, G4double(Z ) ); 348 xsc = bhpair->ComputeCrossSectionPerAtom(theParticleDefinition, energyMscXR, G4double(Z ) ); 338 349 G4cout<<energyMscXR/MeV<<"\t\t"<<xsc/barn<<G4endl; 339 350 writef<<energyMscXR/MeV<<"\t\t"<<xsc/barn<<G4endl;
Note: See TracChangeset
for help on using the changeset viewer.