Changeset 961 for trunk/source/processes/decay/src/G4DecayWithSpin.cc

Timestamp:

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

Author:

garnier

Message:

update processes

File:

• trunk/source/processes/decay/src/G4DecayWithSpin.cc (modified) (4 diffs)

trunk/source/processes/decay/src/G4DecayWithSpin.cc

@@ -29 +29 @@
 //      History:
 //      17 August 2004  P. Gumplinger, T. MacPhail
+//      11 April  2008  Kamil Sedlak (PSI), Toni Shiroka (PSI)
 // ------------------------------------------------------------
 //
@@ -47 +48 @@
 #include "G4Transform3D.hh"
 
-G4DecayWithSpin::G4DecayWithSpin(const G4String& processName):G4Decay(processName){}
+G4DecayWithSpin::G4DecayWithSpin(const G4String& processName):G4Decay(processName)
+{
+  // set Process Sub Type   
+  SetProcessSubType(static_cast<int>(DECAY_WithSpin));
+
+}
 
 G4DecayWithSpin::~G4DecayWithSpin(){}
@@ -97 +103 @@
     if(fieldMgr)field = fieldMgr->GetDetectorField();
 
-    if (field && !(fieldMgr->DoesFieldChangeEnergy())) {
+    if (field) {
 
        G4double point[4];
@@ -105 +111 @@
        point[3] = aTrack.GetGlobalTime();
 
-       G4double fieldValue[3];
+       G4double fieldValue[6];
        field -> GetFieldValue(point,fieldValue);
 
        G4ThreeVector B(fieldValue[0],fieldValue[1],fieldValue[2]);
 
-       parent_polarization = Spin_Precession(aStep,B,fRemainderLifeTime);
+       // Call the spin precession only for non-zero mag. field
+       if (B.mag2() > 0.) parent_polarization = 
+                                 Spin_Precession(aStep,B,fRemainderLifeTime);
 
     }