Changeset 1340 for trunk/source/processes/hadronic/models/incl/src/G4AblaEvaporation.cc

Timestamp:

Nov 5, 2010, 3:45:55 PM (14 years ago)

Author:

garnier

Message:

update ti head

File:

• trunk/source/processes/hadronic/models/incl/src/G4AblaEvaporation.cc (modified) (4 diffs)

trunk/source/processes/hadronic/models/incl/src/G4AblaEvaporation.cc

@@ -24 +24 @@
 // ********************************************************************
 //
-// $Id: G4AblaEvaporation.cc,v 1.4 2008/10/24 21:07:40 dennis Exp $
+// $Id: G4AblaEvaporation.cc,v 1.6 2010/10/26 02:47:59 kaitanie Exp $
 //
 #include <numeric>
@@ -82 +82 @@
   hazard->igraine[17] = 33759;
   hazard->igraine[18] = 13227;
+
+  G4VarNtp *evaporationResult = new G4VarNtp();
+  G4Volant *volant = new G4Volant();
+
+  // Initialize evaporation.
+  abla = new G4Abla(hazard, volant, evaporationResult);
+  abla->initEvapora();
 }
 
@@ -108 +115 @@
 }
 
-G4FragmentVector * G4AblaEvaporation::BreakItUp(const G4Fragment &theNucleus) {
- 
-
-  G4VarNtp *varntp = new G4VarNtp();
-  G4Volant *volant = new G4Volant();
-
-  G4Abla *abla = new G4Abla(hazard, volant, varntp);
-  G4cout <<"Initializing evaporation..." << G4endl;
-  abla->initEvapora();
-  G4cout <<"Initialization complete!" << G4endl;
-  
-  G4double nucleusA = theNucleus.GetA();
-  G4double nucleusZ = theNucleus.GetZ();
+G4FragmentVector * G4AblaEvaporation::BreakItUp(const G4Fragment &theNucleus) {  
+  G4int nucleusA = theNucleus.GetA_asInt();
+  G4int nucleusZ = theNucleus.GetZ_asInt();
   G4double nucleusMass = G4NucleiProperties::GetNuclearMass(nucleusA, nucleusZ);
   G4double excitationEnergy = theNucleus.GetExcitationEnergy();
@@ -137 +134 @@
 
   varntp->ntrack = -1;
-  varntp->massini = theNucleus.GetA();
-  varntp->mzini = theNucleus.GetZ();
+  varntp->massini = theNucleus.GetA_asInt();
+  varntp->mzini = theNucleus.GetZ_asInt();
 
   std::vector<G4DynamicParticle*> cascadeParticles;