Changeset 962 for trunk/source/processes/hadronic/models/cascade/cascade/src/G4PreCompoundCascadeInterface.cc

Timestamp:

Apr 6, 2009, 12:30:29 PM (15 years ago)

Author:

garnier

Message:

update processes

File:

• trunk/source/processes/hadronic/models/cascade/cascade/src/G4PreCompoundCascadeInterface.cc (modified) (6 diffs)

trunk/source/processes/hadronic/models/cascade/cascade/src/G4PreCompoundCascadeInterface.cc

@@ -130 +130 @@
   G4LorentzRotation toLabFrame = toZ.inverse();
 
-  std::vector<G4double> momentumBullet(4);
+  G4CascadeMomentum momentumBullet;
   momentumBullet[0] =0.;
   momentumBullet[1] =0;
@@ -153 +153 @@
   G4DynamicParticle* cascadeParticle = 0;
 
-  std::vector<G4double> targetMomentum(4, 0.0);
+  G4CascadeMomentum targetMomentum;
 
   G4double theNucleusA = theNucleus.GetN();
@@ -163 +163 @@
     target->setEnergy();
 
-    std::vector<G4double>  bmom = bullet->getMomentum();
+    const G4CascadeMomentum& bmom = bullet->getMomentum();
     eInit = std::sqrt(bmom[0] * bmom[0]);
-    std::vector<G4double> tmom = target->getMomentum();
+    const G4CascadeMomentum& tmom = target->getMomentum();
     eInit += std::sqrt(tmom[0] * tmom[0]);
 
@@ -240 +240 @@
     sumBaryon += 1;
 
-    std::vector<G4double> bmom = bullet->getMomentum();
+    const G4CascadeMomentum& bmom = bullet->getMomentum();
     eInit = std::sqrt(bmom[0] * bmom[0]);
-    std::vector<G4double> tmom = targetH->getMomentum();
+    const G4CascadeMomentum& tmom = targetH->getMomentum();
     eInit += std::sqrt(tmom[0] * tmom[0]);
 
@@ -290 +290 @@
     for (ipart = particles.begin(); ipart != particles.end(); ipart++) {
       outgoingParticle = ipart->type();
-      std::vector<G4double> mom = ipart->getMomentum();
+      const G4CascadeMomentum& mom = ipart->getMomentum();
       eTot   += std::sqrt(mom[0] * mom[0]);
 
@@ -421 +421 @@
       {
         G4double eKin = ifrag->getKineticEnergy() * GeV;
-        std::vector<G4double> mom = ifrag->getMomentum();
+        const G4CascadeMomentum& mom = ifrag->getMomentum();
         eTot   += std::sqrt(mom[0] * mom[0]);