8 |
|
|
9 |
gSystem->Load("libG4phys_builders"); |
gSystem->Load("libG4phys_builders"); |
10 |
gSystem->Load("libG4phys_lists"); |
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"); |
|
11 |
|
|
12 |
//set CUTS |
TString PAM_VMC=gSystem->Getenv("PAM_VMC"); |
13 |
// TG4PhysicsListOptions options; |
TString PLAT=gSystem->Getenv("PLATFORM"); |
|
// options.SetSpecialCutsPhysics(true); |
|
|
// options.SetEMPhysics(true); |
|
|
// options.SetMuonPhysics(true); |
|
|
// runConfiguration->SetPhysicsListOptions(options); |
|
14 |
|
|
15 |
|
gSystem->Load(PAM_VMC+"/lib/tgt_"+PLAT+"/libPamG4RunConfiguration"); |
16 |
|
cout<<PAM_VMC+"/lib/tgt_"+PLAT+"/libPamG4RunConfiguration"; |
17 |
|
PamG4RunConfiguration* runConfiguration |
18 |
|
= new PamG4RunConfiguration("geomRoot","QGSP_BIC_HP"); |
19 |
|
|
|
// Switch on optical physics (not activated by default) |
|
|
/// TG4PhysicsListOptions options; |
|
|
////runConfiguration->SetPhysicsListOptions(options); |
|
|
// TGeant4 |
|
20 |
TGeant4* geant4 |
TGeant4* geant4 |
21 |
= new TGeant4("TGeant4_PAMELA", "The Geant4 PAMELA Monte Carlo", runConfiguration); |
= new TGeant4("TGeant4_PAMELA", "The Geant4 PAMELA Monte Carlo", runConfiguration); |
22 |
// gMC->SetCut("CUTGAM", 1.0e-6); |
|
23 |
// gMC->SetCut("BCUTE", 1.0e-6); |
geant4->ProcessGeantMacro(PAM_VMC+"/config/g4config.in"); |
|
// 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, ..) |
|
|
|
|
24 |
|
|
25 |
} |
} |