Changeset 1342 for trunk/examples/extended/electromagnetic/TestEm16/src

Timestamp:

Nov 5, 2010, 4:08:39 PM (14 years ago)

Author:

garnier

Message:

update ti head

File:

• trunk/examples/extended/electromagnetic/TestEm16/src/SteppingAction.cc (modified) (3 diffs)

trunk/examples/extended/electromagnetic/TestEm16/src/SteppingAction.cc

@@ -24 +24 @@
 // ********************************************************************
 //
-// $Id: SteppingAction.cc,v 1.5 2007/01/18 09:07:20 hbu Exp $
-// GEANT4 tag $Name: geant4-09-04-beta-01 $
+// $Id: SteppingAction.cc,v 1.6 2010/10/13 13:50:22 vnivanch Exp $
+// GEANT4 tag $Name: examples-V09-03-09 $
 //
 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
@@ -37 +37 @@
 #include "G4UnitsTable.hh"
 #include "HistoManager.hh"
+#include "G4EmProcessSubType.hh"
 
 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
@@ -53 +54 @@
 void SteppingAction::UserSteppingAction(const G4Step* aStep)
 {
- const G4VProcess* process = aStep->GetPostStepPoint()->GetProcessDefinedStep();
- if (process == 0) return;
- G4String processName = process->GetProcessName();
- static G4int iCalled=0;
- const  G4int nprint=0; // set to 10 to get debug print for first 10 calls
+  const G4VProcess* process = aStep->GetPostStepPoint()->GetProcessDefinedStep();
+  if (process == 0) { return; }
 
- if (processName == "SynRad")
- {
-   iCalled++;
-   G4StepPoint* PrePoint = aStep->GetPreStepPoint();
-   G4TrackVector* secondary = fpSteppingManager->GetSecondary();
+  static G4int iCalled=0;
+  const  G4int nprint=0; // set to 10 to get debug print for first 10 calls
 
-   //(*secondary)[lp-1]  points to the last photon generated
-   //   
-   size_t lp=(*secondary).size();
-   if (lp)
-   {
-     G4double  Egamma =  (*secondary)[lp-1]->GetTotalEnergy();
-     runAction->n_gam_sync++;
-     runAction->e_gam_sync += Egamma;
-     runAction->e_gam_sync2 += Egamma*Egamma;
-     if (Egamma > runAction->e_gam_sync_max) runAction->e_gam_sync_max = Egamma;
-     runAction->lam_gam_sync += aStep->GetStepLength();
-     if (iCalled<nprint)
-     {
-       G4double      Eelec  = PrePoint->GetTotalEnergy();
-       G4ThreeVector Pelec  = PrePoint->GetMomentum();
-       G4ThreeVector PGamma = (*secondary)[lp-1]->GetMomentum();
-       G4bool IsGamma = 
-            ((*secondary)[lp-1]->GetDefinition() == G4Gamma::GammaDefinition());
+  if (fSynchrotronRadiation == process->GetProcessSubType())
+    {
+      ++iCalled;
+      const G4StepPoint* PrePoint = aStep->GetPreStepPoint();
+      const G4TrackVector* secondary = fpSteppingManager->GetSecondary();
+
+      //(*secondary)[lp-1]  points to the last photon generated
+      //   
+      size_t lp=(*secondary).size();
+      if (lp)
+        {
+          G4double  Egamma =  (*secondary)[lp-1]->GetTotalEnergy();
+          runAction->n_gam_sync++;
+          runAction->e_gam_sync += Egamma;
+          runAction->e_gam_sync2 += Egamma*Egamma;
+          if (Egamma > runAction->e_gam_sync_max) 
+            { 
+              runAction->e_gam_sync_max = Egamma;
+            }
+          runAction->lam_gam_sync += aStep->GetStepLength();
+          if (iCalled<nprint)
+            {
+              G4double      Eelec  = PrePoint->GetTotalEnergy();
+              G4ThreeVector Pelec  = PrePoint->GetMomentum();
+              G4ThreeVector PGamma = (*secondary)[lp-1]->GetMomentum();
+              G4bool IsGamma = 
+                ((*secondary)[lp-1]->GetDefinition() == G4Gamma::Gamma());
             
-       G4cout << "UserSteppingAction processName=" << process->GetProcessName()
-         << " Step Length=" << std::setw(6) 
-                            << G4BestUnit(aStep->GetStepLength(),"Length")
-         << " Eelec=" << G4BestUnit(Eelec,"Energy")
-         << " Pelec=" << G4BestUnit(Pelec,"Energy")
-         << " IsGamma=" << IsGamma
-         << " Egamma=" << G4BestUnit(Egamma,"Energy")
-         << " PGamma=" << G4BestUnit(PGamma,"Energy")
-         << " #secondaries lp=" << lp
-         << '\n';
-         }
-     histoManager->FillHisto(1,Egamma);
-     histoManager->FillHisto(2,Egamma,Egamma/keV); // power spectrum : gamma weighted with its energy
-     histoManager->FillHisto(3,aStep->GetStepLength());
-   }
- }
+              G4cout << "UserSteppingAction processName=" << process->GetProcessName()
+                     << " Step Length=" << std::setw(6) 
+                     << G4BestUnit(aStep->GetStepLength(),"Length")
+                     << " Eelec=" << G4BestUnit(Eelec,"Energy")
+                     << " Pelec=" << G4BestUnit(Pelec,"Energy")
+                     << " IsGamma=" << IsGamma
+                     << " Egamma=" << G4BestUnit(Egamma,"Energy")
+                     << " PGamma=" << G4BestUnit(PGamma,"Energy")
+                     << " #secondaries lp=" << lp
+                     << '\n';
+            }
+          histoManager->FillHisto(1,Egamma);
 
+          // power spectrum : gamma weighted with its energy
+          histoManager->FillHisto(2,Egamma,Egamma/keV); 
+
+          histoManager->FillHisto(3,aStep->GetStepLength());
+        }
+    }
 }