Updated g4Config.C in vmctest:
authorihrivnac <ihrivnac@f7af4fe6-9843-0410-8265-dc069ae4e863>
Mon, 17 Jun 2013 16:19:21 +0000 (16:19 +0000)
committerihrivnac <ihrivnac@f7af4fe6-9843-0410-8265-dc069ae4e863>
Mon, 17 Jun 2013 16:19:21 +0000 (16:19 +0000)
- Commented out setting of MaxStepInLowDensityMaterials 1 cm
  (to speed up simulation)
- Inactivated setTrackSecondariesFirst for Cerenkov process
- Activated /mcEvent/printMemory
- Commented out setting /mcVerbose/composedPhysicsList 2
  (so that the default = 1 is used)

test/vmctest/gun/g4Config.C
test/vmctest/lego/g4Config.C
test/vmctest/ppbench/g4Config.C

index b09f009..9838af0 100644 (file)
@@ -47,13 +47,13 @@ void Config()
 //                                "FTFP_BERT_EMV+optical", 
 //                                "specialCuts+specialControls+stackPopper+stepLimiter",
 //                                 true);
-      
+
     geant4 = new TGeant4("TGeant4", 
                          "The Geant4 Monte Carlo : FTFP_BERT_EMV+optical", 
                          runConfiguration);
              // Repeat physics selection in the title; to be removed
              // with new geant4_vmc tag (1.13)            
-                         
+
     cout << "Geant4 has been created." << endl;
   } 
   else {
@@ -66,10 +66,19 @@ void Config()
     geant4->ProcessGeantCommand("/mcVerbose/all 1");  
     geant4->ProcessGeantCommand("/mcVerbose/geometryManager 1");  
     geant4->ProcessGeantCommand("/mcVerbose/opGeometryManager 1");  
-    geant4->ProcessGeantCommand("/mcTracking/loopVerbose 0");     
+    geant4->ProcessGeantCommand("/mcTracking/loopVerbose 1");     
     geant4->ProcessGeantCommand("/mcPhysics/rangeCuts 0.01 mm"); 
     geant4->ProcessGeantCommand("/mcTracking/skipNeutrino true");
-    geant4->ProcessGeantCommand("/mcDet/setMaxStepInLowDensityMaterials 1 cm");
+    
+    // Activate step limit defined in low density materials
+    // (the default value is 10 cm)
+    //geant4->ProcessGeantCommand("/mcDet/setIsMaxStepInLowDensityMaterials true");
+    //geant4->ProcessGeantCommand("/mcDet/setMaxStepInLowDensityMaterials 1 cm");
+
+    // Activate step limit defined in tracking media
+    // (Note: this slows down simulation significantly)
+    //geant4->ProcessGeantCommand("/mcDet/setIsUserMaxStep true");
+   
     // for Geant4 <= 9.4.p03
     //geant4->ProcessGeantCommand("/mcPhysics/selectOpProcess Scintillation");
     //geant4->ProcessGeantCommand("/mcPhysics/setOpProcessActivation false");
@@ -80,7 +89,9 @@ void Config()
     geant4->ProcessGeantCommand("/optics_engine/setOpProcessUse false");
     geant4->ProcessGeantCommand("/optics_engine/selectOpProcess OpMieHG");
     geant4->ProcessGeantCommand("/optics_engine/setOpProcessUse false");
-    
+    geant4->ProcessGeantCommand("/optics_engine/selectOpProcess Cerenkov");
+    geant4->ProcessGeantCommand("/optics_engine/setTrackSecondariesFirst false");   
+
     // Activate saving random engine status
     // (the file per event will be re-written with each new event)
     //gAlice->GetMCApp()->SetSaveRndmStatus(kTRUE);
@@ -91,6 +102,9 @@ void Config()
     //gAlice->GetMCApp()->SetSaveRndmStatusPerEvent(kTRUE);
     //geant4->ProcessGeantCommand("/mcRun/saveRandom true");
     //geant4->ProcessGeantCommand("/mcEvent/saveRandom true");
+    
+    // Activate printing size of used memory per event
+    geant4->ProcessGeantCommand("/mcEvent/printMemory true");
 
   // Uncomment this line to get a detail info from each step 
   //geant4->ProcessGeantCommand("/tracking/verbose 1");  
index 4d85f61..eb0ca2e 100644 (file)
@@ -49,7 +49,7 @@ void Config(const TString& det)
 //                                 true);
       
     geant4 = new TGeant4("TGeant4", 
-                         "The Geant4 Monte Carlo : QGSP_BERT_EMV+optical", 
+                         "The Geant4 Monte Carlo : FTFP_BERT_EMV+optical", 
                          runConfiguration);
              // Repeat physics selection in the title; to be removed
              // with new geant4_vmc tag (1.13)            
@@ -66,10 +66,19 @@ void Config(const TString& det)
     geant4->ProcessGeantCommand("/mcVerbose/all 1");  
     geant4->ProcessGeantCommand("/mcVerbose/geometryManager 1");  
     geant4->ProcessGeantCommand("/mcVerbose/opGeometryManager 1");  
-    geant4->ProcessGeantCommand("/mcTracking/loopVerbose 0");     
+    geant4->ProcessGeantCommand("/mcTracking/loopVerbose 1");     
     geant4->ProcessGeantCommand("/mcPhysics/rangeCuts 0.01 mm"); 
     geant4->ProcessGeantCommand("/mcTracking/skipNeutrino true");
-    geant4->ProcessGeantCommand("/mcDet/setMaxStepInLowDensityMaterials 1 cm");
+    
+    // Activate step limit defined in low density materials
+    // (the default value is 10 cm)
+    //geant4->ProcessGeantCommand("/mcDet/setIsMaxStepInLowDensityMaterials true");
+    //geant4->ProcessGeantCommand("/mcDet/setMaxStepInLowDensityMaterials 1 cm");
+
+    // Activate step limit defined in tracking media
+    // (Note: this slows down simulation significantly)
+    //geant4->ProcessGeantCommand("/mcDet/setIsUserMaxStep true");
+   
     // for Geant4 <= 9.4.p03
     //geant4->ProcessGeantCommand("/mcPhysics/selectOpProcess Scintillation");
     //geant4->ProcessGeantCommand("/mcPhysics/setOpProcessActivation false");
@@ -80,7 +89,9 @@ void Config(const TString& det)
     geant4->ProcessGeantCommand("/optics_engine/setOpProcessUse false");
     geant4->ProcessGeantCommand("/optics_engine/selectOpProcess OpMieHG");
     geant4->ProcessGeantCommand("/optics_engine/setOpProcessUse false");
-    
+    geant4->ProcessGeantCommand("/optics_engine/selectOpProcess Cerenkov");
+    geant4->ProcessGeantCommand("/optics_engine/setTrackSecondariesFirst false");   
+
     // Activate saving random engine status
     // (the file per event will be re-written with each new event)
     //gAlice->GetMCApp()->SetSaveRndmStatus(kTRUE);
@@ -91,7 +102,9 @@ void Config(const TString& det)
     //gAlice->GetMCApp()->SetSaveRndmStatusPerEvent(kTRUE);
     //geant4->ProcessGeantCommand("/mcRun/saveRandom true");
     //geant4->ProcessGeantCommand("/mcEvent/saveRandom true");
-
+    
+    // Activate printing size of used memory per event
+    geant4->ProcessGeantCommand("/mcEvent/printMemory true");
 
   // Uncomment this line to get a detail info from each step 
   //geant4->ProcessGeantCommand("/tracking/verbose 1");  
index 7e26935..63a6c41 100644 (file)
@@ -69,7 +69,16 @@ void Config()
     geant4->ProcessGeantCommand("/mcTracking/loopVerbose 1");     
     geant4->ProcessGeantCommand("/mcPhysics/rangeCuts 0.01 mm"); 
     geant4->ProcessGeantCommand("/mcTracking/skipNeutrino true");
-    geant4->ProcessGeantCommand("/mcDet/setMaxStepInLowDensityMaterials 1 cm");
+    
+    // Activate step limit defined in low density materials
+    // (the default value is 10 cm)
+    //geant4->ProcessGeantCommand("/mcDet/setIsMaxStepInLowDensityMaterials true");
+    //geant4->ProcessGeantCommand("/mcDet/setMaxStepInLowDensityMaterials 1 cm");
+
+    // Activate step limit defined in tracking media
+    // (Note: this slows down simulation significantly)
+    //geant4->ProcessGeantCommand("/mcDet/setIsUserMaxStep true");
+   
     // for Geant4 <= 9.4.p03
     //geant4->ProcessGeantCommand("/mcPhysics/selectOpProcess Scintillation");
     //geant4->ProcessGeantCommand("/mcPhysics/setOpProcessActivation false");
@@ -80,7 +89,9 @@ void Config()
     geant4->ProcessGeantCommand("/optics_engine/setOpProcessUse false");
     geant4->ProcessGeantCommand("/optics_engine/selectOpProcess OpMieHG");
     geant4->ProcessGeantCommand("/optics_engine/setOpProcessUse false");
-    
+    geant4->ProcessGeantCommand("/optics_engine/selectOpProcess Cerenkov");
+    geant4->ProcessGeantCommand("/optics_engine/setTrackSecondariesFirst false");   
+
     // Activate saving random engine status
     // (the file per event will be re-written with each new event)
     //gAlice->GetMCApp()->SetSaveRndmStatus(kTRUE);
@@ -91,6 +102,9 @@ void Config()
     //gAlice->GetMCApp()->SetSaveRndmStatusPerEvent(kTRUE);
     //geant4->ProcessGeantCommand("/mcRun/saveRandom true");
     //geant4->ProcessGeantCommand("/mcEvent/saveRandom true");
+    
+    // Activate printing size of used memory per event
+    geant4->ProcessGeantCommand("/mcEvent/printMemory true");
 
   // Uncomment this line to get a detail info from each step 
   //geant4->ProcessGeantCommand("/tracking/verbose 1");  
@@ -98,7 +112,7 @@ void Config()
   // More info from the physics list
   // the verbosity level is passed to all contained physics lists and their
   // physics builders
-  geant4->ProcessGeantCommand("/mcVerbose/composedPhysicsList 2");  
+  //geant4->ProcessGeantCommand("/mcVerbose/composedPhysicsList 2");  
   
   // More info from optical processes
   //geant4->ProcessGeantCommand("/mcVerbose/opticalPhysicsList 3");