Changeset 961 for trunk/source/processes/electromagnetic/utils/src/G4EmProcessOptions.cc

Timestamp:

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

Author:

garnier

Message:

update processes

File:

• trunk/source/processes/electromagnetic/utils/src/G4EmProcessOptions.cc (modified) (5 diffs)

trunk/source/processes/electromagnetic/utils/src/G4EmProcessOptions.cc

@@ -24 +24 @@
 // ********************************************************************
 //
-// $Id: G4EmProcessOptions.cc,v 1.22 2007/11/07 18:38:49 vnivanch Exp $
-// GEANT4 tag $Name: geant4-09-01-patch-02 $
+// $Id: G4EmProcessOptions.cc,v 1.26 2009/02/18 14:43:27 vnivanch Exp $
+// GEANT4 tag $Name: geant4-09-02-ref-02 $
 //
 // -------------------------------------------------------------------
@@ -59 +59 @@
 #include "G4VMultipleScattering.hh"
 #include "G4Region.hh"
+#include "G4RegionStore.hh"
 
 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
@@ -392 +393 @@
 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
 
-void G4EmProcessOptions::ActivateDeexcitation(G4bool val, const G4Region* r)
-{
-  const std::vector<G4VEnergyLossProcess*>& v =
-        theManager->GetEnergyLossProcessVector();
-  std::vector<G4VEnergyLossProcess*>::const_iterator itr;
-  for(itr = v.begin(); itr != v.end(); itr++) {
-    G4VEnergyLossProcess* p = *itr;
-    if(p) p->ActivateDeexcitation(val,r);
-  }
-  const std::vector<G4VEmProcess*>& w =
-        theManager->GetEmProcessVector();
-  std::vector<G4VEmProcess*>::const_iterator itp;
-  for(itp = w.begin(); itp != w.end(); itp++) {
-    G4VEmProcess* q = *itp;
-    if(q) q->ActivateDeexcitation(val,r);
+void G4EmProcessOptions::ActivateDeexcitation(const G4String& pname, 
+                                              G4bool val, 
+                                              const G4String& reg)
+{
+  G4RegionStore* regionStore = G4RegionStore::GetInstance();
+  const G4Region* r = 0;
+  if(reg == "" || reg == "World") {
+    r = regionStore->GetRegion("DefaultRegionForTheWorld", false);
+  } else {
+    r = regionStore->GetRegion(reg, false);
+  }
+  if(!r) {
+    G4cout << "G4EmProcessOptions::ActivateDeexcitation ERROR: G4Region <"
+           << reg << "> not found, the command ignored" << G4endl;
+    return;
+  }
+
+  const std::vector<G4VEnergyLossProcess*>& v =
+        theManager->GetEnergyLossProcessVector();
+  std::vector<G4VEnergyLossProcess*>::const_iterator itr;
+  for(itr = v.begin(); itr != v.end(); itr++) {
+    G4VEnergyLossProcess* p = *itr;
+    if(p) {
+      if(pname == p->GetProcessName()) p->ActivateDeexcitation(val,r);
+    }
+  }
+  const std::vector<G4VEmProcess*>& w =
+        theManager->GetEmProcessVector();
+  std::vector<G4VEmProcess*>::const_iterator itp;
+  for(itp = w.begin(); itp != w.end(); itp++) {
+    G4VEmProcess* q = *itp;
+    if(q) {
+      if(pname == q->GetProcessName()) q->ActivateDeexcitation(val,r);
+    }
   }
 }
@@ -477 +497 @@
 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
 
+void G4EmProcessOptions::SetPolarAngleLimit(G4double val)
+{
+  const std::vector<G4VMultipleScattering*>& u =
+        theManager->GetMultipleScatteringVector();
+  std::vector<G4VMultipleScattering*>::const_iterator itm;
+  for(itm = u.begin(); itm != u.end(); itm++) {
+    if(*itm) (*itm)->SetPolarAngleLimit(val);
+  }
+  const std::vector<G4VEmProcess*>& w =
+        theManager->GetEmProcessVector();
+  std::vector<G4VEmProcess*>::const_iterator itp;
+  for(itp = w.begin(); itp != w.end(); itp++) {
+    G4VEmProcess* q = *itp;
+    if(q) q->SetPolarAngleLimit(val);
+  }
+}
+
+//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
+
 void G4EmProcessOptions::SetLPMFlag(G4bool val)
 {
@@ -484 +523 @@
 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
 
+void G4EmProcessOptions::SetSplineFlag(G4bool val)
+{
+  theManager->SetSplineFlag(val);
+}
+
+//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
+
 void G4EmProcessOptions::SetLinearLossLimit(G4double val)
 {