In vmctest:
authorivana <ivana@f7af4fe6-9843-0410-8265-dc069ae4e863>
Tue, 5 Jul 2011 12:56:44 +0000 (12:56 +0000)
committerivana <ivana@f7af4fe6-9843-0410-8265-dc069ae4e863>
Tue, 5 Jul 2011 12:56:44 +0000 (12:56 +0000)
Updating G3/G4 configs to be consistent with the configuration used in
G4 test productions

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

index 00e5eeb..c707bba 100644 (file)
@@ -12,7 +12,7 @@ void Config()
   // AliRoot setup
   //
   gROOT->LoadMacro("$ALICE_ROOT/test/vmctest/gun/commonConfig.C");
-  commonConfig(kFALSE);
+  commonConfig();
 
   // Load Geant3 + Geant3 VMC libraries
   //
index 7d780e7..579b2eb 100644 (file)
@@ -13,7 +13,16 @@ void Config()
   // AliRoot setup
   //
   gROOT->LoadMacro("$ALICE_ROOT/test/vmctest/gun/commonConfig.C");
-  commonConfig(kTRUE);
+  commonConfig();
+
+  // TPC primary ionization 
+  AliTPC* TPC = (AliTPC*)gAlice->GetDetector("TPC");
+  if ( ! TPC )
+    cerr << "Cannot get TPC detector" << endl;
+  else  { 
+    cerr << "Setting TPC primary ionization" << endl;
+    TPC->SetPrimaryIonisation(); // not used with Geant3
+  }  
 
   // Load Geant4 + Geant4 VMC libraries
   //
@@ -50,19 +59,18 @@ void Config()
   // Customization of Geant4 VMC
   //
 
-  //geant4->ProcessGeantCommand("/mcPhysics/rangeCuts 0.01 mm"); 
     geant4->ProcessGeantCommand("/mcVerbose/all 1");  
-    geant4->ProcessGeantCommand("/mcVerbose/geometryManager 2");  
+    geant4->ProcessGeantCommand("/mcVerbose/geometryManager 1");  
     geant4->ProcessGeantCommand("/mcVerbose/opGeometryManager 1");  
     geant4->ProcessGeantCommand("/mcTracking/loopVerbose 0");     
     geant4->ProcessGeantCommand("/mcPhysics/rangeCuts 0.01 mm"); 
-    //geant4->ProcessGeantCommand("/mcPhysics/rangeCuts 1 mm"); 
     geant4->ProcessGeantCommand("/mcPhysics/selectOpProcess Scintillation");
     geant4->ProcessGeantCommand("/mcPhysics/setOpProcessActivation false");
     geant4->ProcessGeantCommand("/mcTracking/skipNeutrino true");
+    geant4->ProcessGeantCommand("/mcDet/setMaxStepInLowDensityMaterials 1 cm");
 
   // Uncomment this line to get a detail info from each step 
-  // geant4->ProcessGeantCommand("/tracking/verbose 1");  
+  //geant4->ProcessGeantCommand("/tracking/verbose 1");  
   
   // More info from the physics list
   // the verbosity level is passed to all contained physics lists and their
@@ -82,6 +90,7 @@ void Config()
   // More info about regions construction 
   // and conversion of VMC cuts in cuts in range per regions 
   // geant4->ProcessGeantCommand("/mcVerbose/regionsManager 2");
+  // geant4->ProcessGeantCommand("/mcRegions/print true");
   
   // Suppress verbose info from tracks which reached maximum number of steps
   // (default value is 30000)  
index ed44b3c..eea0031 100644 (file)
@@ -15,6 +15,15 @@ void Config()
   gROOT->LoadMacro("$ALICE_ROOT/test/vmctest/ppbench/commonConfig.C");
   commonConfig();
 
+  // TPC primary ionization 
+  AliTPC* TPC = (AliTPC*)gAlice->GetDetector("TPC");
+  if ( ! TPC )
+    cerr << "Cannot get TPC detector" << endl;
+  else  { 
+    cerr << "Setting TPC primary ionization" << endl;
+    TPC->SetPrimaryIonisation(); // not used with Geant3
+  }  
+
   // Load Geant4 + Geant4 VMC libraries
   //
   if (gClassTable->GetID("TGeant4") == -1) {
@@ -50,18 +59,15 @@ void Config()
   // Customization of Geant4 VMC
   //
 
-  //geant4->ProcessGeantCommand("/mcPhysics/rangeCuts 0.01 mm"); 
     geant4->ProcessGeantCommand("/mcVerbose/all 1");  
-    geant4->ProcessGeantCommand("/mcVerbose/geometryManager 2");  
+    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("/mcPhysics/rangeCuts 1 mm"); 
     geant4->ProcessGeantCommand("/mcPhysics/selectOpProcess Scintillation");
     geant4->ProcessGeantCommand("/mcPhysics/setOpProcessActivation false");
     geant4->ProcessGeantCommand("/mcTracking/skipNeutrino true");
-    //geant4->ProcessGeantCommand("/mcMagField/stepperType HelixImplicitEuler");
-    //geant4->ProcessGeantCommand("/mcMagField/stepperType HelixSimpleRunge");
+    geant4->ProcessGeantCommand("/mcDet/setMaxStepInLowDensityMaterials 1 cm");
 
   // Uncomment this line to get a detail info from each step 
   // geant4->ProcessGeantCommand("/tracking/verbose 1");