Changeset 1315 for trunk/source/processes/hadronic/models/cascade/cascade/include/G4InuclNuclei.hh

Timestamp:

Jun 18, 2010, 11:42:07 AM (14 years ago)

Author:

garnier

Message:

update geant4-09-04-beta-cand-01 interfaces-V09-03-09 vis-V09-03-08

File:

• trunk/source/processes/hadronic/models/cascade/cascade/include/G4InuclNuclei.hh (modified) (2 diffs)

trunk/source/processes/hadronic/models/cascade/cascade/include/G4InuclNuclei.hh

@@ +1 @@
+#ifndef G4INUCL_NUCLEI_HH
+#define G4INUCL_NUCLEI_HH
 //
 // ********************************************************************
@@ -23 +25 @@
 // * acceptance of all terms of the Geant4 Software license.          *
 // ********************************************************************
+// $Id: G4InuclNuclei.hh,v 1.16 2010/04/09 19:33:11 mkelsey Exp $
+// Geant4 tag: $Name: geant4-09-04-beta-cand-01 $
 //
-#ifndef G4INUCL_NUCLEI_HH
-#define G4INUCL_NUCLEI_HH
+// 20100112  Michael Kelsey -- Replace G4CascadeMomentum with G4LorentzVector
+// 20100301  M. Kelsey -- Add function to create unphysical nuclei for use
+//           as temporary final-state fragments.
+// 20100319  M. Kelsey -- Remove "using" directory and unnecessary #includes.
+// 20100409  M. Kelsey -- Drop unused string argument from ctors.
 
-#ifndef G4INUCL_PARTICLE_HH
 #include "G4InuclParticle.hh"
-#endif
 #include "G4ExitonConfiguration.hh"
-#include "G4InuclSpecialFunctions.hh"
 
-using namespace G4InuclSpecialFunctions;
+class G4ParticleDefinition;
 
 class G4InuclNuclei : public G4InuclParticle {
+public:
+  G4InuclNuclei() : G4InuclParticle() {}
 
-public:
+  G4InuclNuclei(G4double a, G4double z)
+    : G4InuclParticle(makeDefinition(a,z)),
+      exitationEnergy(0.0) {}
 
-  G4InuclNuclei() {};
+  G4InuclNuclei(const G4LorentzVector& mom, G4double a, G4double z)
+    : G4InuclParticle(makeDefinition(a,z), mom),
+      exitationEnergy(0.0) {}
 
-  G4InuclNuclei(G4double a, 
-                G4double z)
-    : A(a),
-      Z(z) { 
+  G4InuclNuclei(G4double ekin, G4double a, G4double z) 
+    : G4InuclParticle(makeDefinition(a,z), ekin),
+      exitationEnergy(0.0) {}
 
-    setNucleiMass();
-    exitationEnergy = 0.0; 
-  };
+  virtual ~G4InuclNuclei() {}
 
-  G4InuclNuclei(const G4CascadeMomentum& mom, 
-                G4double a, 
-                G4double z)
-    : G4InuclParticle(mom),
-      A(a),
-      Z(z) { 
+  // Copy and assignment constructors for use with std::vector<>
+  G4InuclNuclei(const G4InuclNuclei& right)
+    : G4InuclParticle(right), exitationEnergy(right.exitationEnergy),
+      theExitonConfiguration(right.theExitonConfiguration) {}
 
-    setNucleiMass();  
-    exitationEnergy = 0.0;
-  };
+  G4InuclNuclei& operator=(const G4InuclNuclei& right);
 
-  G4InuclNuclei(G4double ekin, 
-                G4double a, 
-                G4double z) 
-    : A(a), 
-    Z(z) { 
-
-    setNucleiMass();
-    G4CascadeMomentum mom;
-    mom[0] = ekin + nucleiMass;
-    mom[3] = std::sqrt(mom[0] * mom[0] - nucleiMass * nucleiMass);
-    G4InuclParticle::setMomentum(mom);
-    exitationEnergy = 0.0;
-  };
-        
-  void setA(G4double a) { 
-
-    A = a; 
-  };
-
-  void setZ(G4double z) { 
-
-    Z = z; 
-  };
-
-  void setExitationEnergy(G4double e) { 
-
-    exitationEnergy = e; 
-  };
-
-  void setEnergy() {
-
-    momentum[0] = std::sqrt(momentum[1] * momentum[1] + momentum[2] * momentum[2] +
-                       momentum[3] * momentum[3] + nucleiMass * nucleiMass);  
-  };
-
-  void setNucleiMass() {
-
-    nucleiMass = 0.93827 * Z + 0.93957 * (A - Z) - 0.001 * bindingEnergy(A, Z);
-  };
+  void setExitationEnergy(G4double e) { exitationEnergy = e; }
 
   void setExitonConfiguration(const G4ExitonConfiguration& config) { 
+    theExitonConfiguration = config;
+  }
 
-    theExitonConfiguration = config;
-  };
+  G4double getA() const { return getDefinition()->GetAtomicMass(); }
 
-  G4double getA() const { 
+  G4double getZ() const { return getDefinition()->GetAtomicNumber(); }
 
-    return A; 
-  };
-
-  G4double getZ() const { 
-
-    return Z; 
-  };
-
-  G4double getExitationEnergy() const { 
-
-    return exitationEnergy; 
-  };
+  G4double getExitationEnergy() const { return exitationEnergy; }
 
   G4double getExitationEnergyInGeV() const { 
+    return 0.001 * exitationEnergy; 
+  }
 
-    return 0.001 * exitationEnergy; 
-  };
+  const G4ExitonConfiguration& getExitonConfiguration() const {
+    return theExitonConfiguration;
+  }
 
-  G4ExitonConfiguration getExitonConfiguration() const {
-
-    return theExitonConfiguration;
-  };
-
-  G4double getMass() const { 
-
-    return nucleiMass; 
-  };
-
-  G4double getKineticEnergy() const { 
-
-    return momentum[0] - nucleiMass; 
-  };
-
-  G4double getNucleiMass(G4double a, 
-                         G4double z) const {
-
-    return 0.93827 * z + 0.93957 * (a - z) - 0.001 * bindingEnergy(a, z);
-  };
+  static G4double getNucleiMass(G4double a, G4double z);
 
   virtual void printParticle() const {
+    G4cout << " A " << getA() << " Z " << getZ() << " mass " 
+           << getMass() << " Eex (MeV) " << exitationEnergy << G4endl;
+    G4InuclParticle::printParticle();
+  }
 
-    G4cout << " A " << A << " Z " << Z << " mass " << nucleiMass << 
-      " Eex (MeV) " << exitationEnergy << G4endl;
+protected:
+  // Convert nuclear configuration to standard GEANT4 pointer
+  static G4ParticleDefinition*
+  makeDefinition(G4double a, G4double z, G4double exc=0.);
 
-    G4cout << " Px " << momentum[1] << " Py " << momentum[2] << " Pz " <<  
-        momentum[3] <<  " E " << momentum[0] << G4endl;
-  };
+  static G4ParticleDefinition* 
+  makeNuclearFragment(G4double a, G4double z, G4double exc=0.);
 
 private: 
-
-  G4double A;
-  G4double Z;
   G4double exitationEnergy;
-  G4double nucleiMass;
   G4ExitonConfiguration theExitonConfiguration;
-
 };