trieste/pamVMC/g4Config.C

// $Id: g4Config.C,v 1.0 2007/06/01  Exp $
//
// Configuration macro for Geant4 VirtualMC for PamG4VMC
// for geometry defined with Root and selected Geant4 native navigation

void Config()
{
  
  gSystem->Load("libG4phys_builders");
  gSystem->Load("libG4phys_lists");
  TString G4WD=gSystem->Getenv("G4WORKDIR");
  TString PLAT=gSystem->Getenv("PLATFORM");
  //gSystem->Load(G4WD+"/lib/tgt_"+PLAT+"/libPamG4RunConfiguration");
  //cout<<G4WD+"/lib/tgt_"+PLAT+"/libPamG4RunConfiguration";
    //PamG4RunConfiguration* runConfiguration 
    //= new PamG4RunConfiguration("geomRootToGeant4");
    TG4RunConfiguration* runConfiguration 
    = new TG4RunConfiguration("geomRootToGeant4","QGSP_BERT");
    //= new TG4RunConfiguration("geomRootToGeant4","QGSP_BIC_HP");
  
    //    TG4RunConfiguration* runConfiguration 
    //  = new TG4RunConfiguration("geomRootToGeant4");

   //set CUTS
//   TG4PhysicsListOptions options;
//   options.SetSpecialCutsPhysics(true);
//  options.SetEMPhysics(true);
//  options.SetMuonPhysics(true);
//    runConfiguration->SetPhysicsListOptions(options);


  // Switch on optical physics (not activated by default)
  /// TG4PhysicsListOptions options;
  ////runConfiguration->SetPhysicsListOptions(options);
  // TGeant4
  TGeant4* geant4
    = new TGeant4("TGeant4_PAMELA", "The Geant4 PAMELA Monte Carlo", runConfiguration);
//   gMC->SetCut("CUTGAM", 1.0e-6);
//   gMC->SetCut("BCUTE", 1.0e-6);
//   gMC->SetCut("CUTELE", 1.0e-6);
//   gMC->SetCut("DCUTE", 1.0e-6);
//   gMC->SetCut("CUTNEU", 1.0e-6);
//   gMC->SetCut("CUTHAD", 1.0e-6);
//   gMC->SetCut("CUTMUO", 1.0e-6);
//   gMC->SetCut("DCUTM", 1.0e-6);
//   gMC->SetCut("BCUTM", 1.0e-6);
//   gMC->SetCut("PPCUTM", 1.0e-6);
//   gMC->SetCut("TOFMAX", 1.0e-3);
  //  gMC->Gstpar(10,"LOSS",0.1);

  //  TVirtualMCGeometry *geo = TG4GeometryManager::Instance()->GetMCGeometry();  
  geant4->ProcessGeantMacro("g4config.in");
  // Customise Geant4 setting
  // (verbose level, global range cut, ..)


}
1	// $Id: g4Config.C,v 1.0 2007/06/01 Exp $
2	//
3	// Configuration macro for Geant4 VirtualMC for PamG4VMC
4	// for geometry defined with Root and selected Geant4 native navigation
5
6	void Config()
7	{
8
9	gSystem->Load("libG4phys_builders");
10	gSystem->Load("libG4phys_lists");
11	TString G4WD=gSystem->Getenv("G4WORKDIR");
12	TString PLAT=gSystem->Getenv("PLATFORM");
13	//gSystem->Load(G4WD+"/lib/tgt_"+PLAT+"/libPamG4RunConfiguration");
14	//cout<<G4WD+"/lib/tgt_"+PLAT+"/libPamG4RunConfiguration";
15	//PamG4RunConfiguration* runConfiguration
16	//= new PamG4RunConfiguration("geomRootToGeant4");
17	TG4RunConfiguration* runConfiguration
18	= new TG4RunConfiguration("geomRootToGeant4","QGSP_BERT");
19	//= new TG4RunConfiguration("geomRootToGeant4","QGSP_BIC_HP");
20
21	// TG4RunConfiguration* runConfiguration
22	// = new TG4RunConfiguration("geomRootToGeant4");
23
24	//set CUTS
25	// TG4PhysicsListOptions options;
26	// options.SetSpecialCutsPhysics(true);
27	// options.SetEMPhysics(true);
28	// options.SetMuonPhysics(true);
29	// runConfiguration->SetPhysicsListOptions(options);
30
31
32	// Switch on optical physics (not activated by default)
33	/// TG4PhysicsListOptions options;
34	////runConfiguration->SetPhysicsListOptions(options);
35	// TGeant4
36	TGeant4* geant4
37	= new TGeant4("TGeant4_PAMELA", "The Geant4 PAMELA Monte Carlo", runConfiguration);
38	// gMC->SetCut("CUTGAM", 1.0e-6);
39	// gMC->SetCut("BCUTE", 1.0e-6);
40	// gMC->SetCut("CUTELE", 1.0e-6);
41	// gMC->SetCut("DCUTE", 1.0e-6);
42	// gMC->SetCut("CUTNEU", 1.0e-6);
43	// gMC->SetCut("CUTHAD", 1.0e-6);
44	// gMC->SetCut("CUTMUO", 1.0e-6);
45	// gMC->SetCut("DCUTM", 1.0e-6);
46	// gMC->SetCut("BCUTM", 1.0e-6);
47	// gMC->SetCut("PPCUTM", 1.0e-6);
48	// gMC->SetCut("TOFMAX", 1.0e-3);
49	// gMC->Gstpar(10,"LOSS",0.1);
50
51	// TVirtualMCGeometry *geo = TG4GeometryManager::Instance()->GetMCGeometry();
52	geant4->ProcessGeantMacro("g4config.in");
53	// Customise Geant4 setting
54	// (verbose level, global range cut, ..)
55
56
57	}