Changeset 1347 for trunk/source/processes/electromagnetic/xrays/test/testG4HeatedKleinNishina.cc

Timestamp:

Dec 22, 2010, 3:52:27 PM (14 years ago)

Author:

garnier

Message:

geant4 tag 9.4

File:

• trunk/source/processes/electromagnetic/xrays/test/testG4HeatedKleinNishina.cc (modified) (4 diffs)

trunk/source/processes/electromagnetic/xrays/test/testG4HeatedKleinNishina.cc

@@ -61 +61 @@
 #include "G4HeatedKleinNishinaCompton.hh"
 #include "G4PEEffectModel.hh"
+#include "G4PairProductionRelModel.hh"
+#include "G4BetheHeitlerModel.hh"
 
 
@@ -102 +104 @@
   G4int choice;
   // G4cin >> choice;
-  choice = 54;
+  choice = 6;
 
 
@@ -298 +300 @@
   kinEnergy = 25.0*GeV;
 
-  // G4VEmModel* comp = new G4KleinNishinaCompton();
-  G4VEmModel* comp = new G4HeatedKleinNishinaCompton();
+  G4VEmModel* comp = new G4KleinNishinaCompton();
+  // G4VEmModel* comp = new G4HeatedKleinNishinaCompton();
   G4VEmModel* photo = new G4PEEffectModel();
-
+  G4VEmModel* pair = new G4PairProductionRelModel(theParticleDefinition,"pp");
+  G4VEmModel* bhpair = new G4BetheHeitlerModel(theParticleDefinition,"bhpp");
 
   G4DynamicParticle*  theDynamicParticle = new G4DynamicParticle(theParticleDefinition,
@@ -326 +329 @@
 
 
+  //  energyMscXR = 2*MeV;
+  //xsc = pair->ComputeCrossSectionPerAtom(theParticleDefinition, energyMscXR, G4double(Z ) );
+
 
   iMax = 130;
 
   writef<<iMax<<G4endl;
+
+
+
 
   for( i = 0; i < iMax; i++ )
   {
     energyMscXR = std::exp(i*0.2)*0.001*keV;
-    //xsc = comp->ComputeCrossSectionPerAtom(theParticleDefinition,energyMscXR,G4double(Z ) );
-    xsc = photo->ComputeCrossSectionPerAtom(theParticleDefinition,energyMscXR,G4double(Z ) );
+    // xsc = comp->ComputeCrossSectionPerAtom(theParticleDefinition,energyMscXR,G4double(Z ) );
+    // xsc = photo->ComputeCrossSectionPerAtom(theParticleDefinition, energyMscXR, G4double(Z ) );
+    // xsc = pair->ComputeCrossSectionPerAtom(theParticleDefinition, energyMscXR, G4double(Z ) );
+    xsc = bhpair->ComputeCrossSectionPerAtom(theParticleDefinition, energyMscXR, G4double(Z ) );
     G4cout<<energyMscXR/MeV<<"\t\t"<<xsc/barn<<G4endl;
     writef<<energyMscXR/MeV<<"\t\t"<<xsc/barn<<G4endl;