| 1 |
pamelats |
1.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 |
|
|
} |
Parent Directory
| 
Revision Log