Changeset 1347 for trunk/source/processes/hadronic/models/abrasion/src/G4NuclearAbrasionGeometry.cc

Timestamp:

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

Author:

garnier

Message:

geant4 tag 9.4

File:

• trunk/source/processes/hadronic/models/abrasion/src/G4NuclearAbrasionGeometry.cc (modified) (3 diffs)

trunk/source/processes/hadronic/models/abrasion/src/G4NuclearAbrasionGeometry.cc

@@ -194 +194 @@
   return Es;
 }
-//////////////////////////////////////////////////////////////////////////////////////
-//
+
+
 G4double G4NuclearAbrasionGeometry::GetExcitationEnergyOfTarget ()
 {
-//
-//
-// This member function declares a new G4NuclearAbrasionGeometry object but with the
-// projectile and target exchanged to determine the values for F and P.  Determination
-// of the excess surface area and excitation energy is as above.
-//
-  G4NuclearAbrasionGeometry *revAbrasionGeometry =
+  // This member function declares a new G4NuclearAbrasionGeometry object 
+  // but with the projectile and target exchanged to determine the values
+  // for F and P.  Determination of the excess surface area and excitation
+  // energy is as above.
+
+  G4NuclearAbrasionGeometry* revAbrasionGeometry =
     new G4NuclearAbrasionGeometry(AT, AP, r);
   G4double F1 = revAbrasionGeometry->F();
@@ -212 +211 @@
   Es = 0.95 * MeV * 4.0 * pi * rT*rT/fermi/fermi *
        (1.0+P1-std::pow(1.0-F1,2.0/3.0));
+
 //  if (rP < rT && r < rT-rP)
-  if ((r-rT)/rP < rth)
-  {
+  if ((r-rT)/rP < rth) {
     G4double omega = 0.0;
     if      (AT < 12.0)  omega = 1500.0;
@@ -225 +224 @@
   else if (Es > B * AT)
     Es = B * AT;
+
+  delete revAbrasionGeometry;
+
   return Es;
 }
-////////////////////////////////////////////////////////////////////////////////
-//