PAMELA%20software/PamCAL/run_g4.cxx

#include "TROOT.h"
#include "TVirtualMC.h"
#include "TGeant4.h"
#include "TString.h"
#include "PamVMCApplication.h"

#include <iostream>
// $Id: run_g4.C,v 1.0 2007/06/01 
//
// Macro for running G4PAM_VMC with Geant4. 

using std::cout;
using std::endl;

//void run_g4(const TString& configMacro = "g4Config.C") 
  int main(){
    TROOT nik("Nikolay","Nikonov"); 

    TString configMacro = "g4Config.C"; 
    // Load basic libraries
    gROOT->LoadMacro("./macro/basiclibs.C");
    basiclibs();
    
    // Load PAMELA libs
    TString PAMLIB=gSystem->Getenv("PAM_LIB");
    gSystem->Load(PAMLIB+"/libyoda.so");
    gSystem->Load(PAMLIB+"/libDarthVader.so");
    gSystem->Load(PAMLIB+"/libPamLevel2.so");
    
    // Load Geant4 libraries
    
    gROOT->LoadMacro("./macro/PamG4Libs.C");
    PamG4Libs();
    
    // Load this example libraries
    TString G4WD=gSystem->Getenv("G4WORKDIR");
    TString PLAT=gSystem->Getenv("PLATFORM");
    gSystem->Load(G4WD+"/lib/tgt_"+PLAT+"/libPamVMC");
    
    // MC application
    PamVMCApplication* appl 
      = new PamVMCApplication("PAMG4_VMC", "PAMELA GEANT4 VMC application");
    
    appl->GetPrimaryGenerator()->SetNofPrimaries(1);  
    appl->SetVerboseLevel(1);  
    
    // Macro to run with primary = polarized optical photon
    //gROOT->LoadMacro("optPhoton.C");
    //optPhoton();  
    
    // Initialize MC
    appl->InitMC(configMacro);
    
    // Customise Geant4 setting after initialization:
    // Physics list
    ((TGeant4*)gMC)->ProcessGeantMacro("g4config2.in");
    // Visualization settings
    //   ((TGeant4*)gMC)->ProcessGeantMacro("g4vis.in");
    
    // Run MC
    appl->RunMC(5000);
    appl->FinishRun();
    delete appl;
    return 0;
  }  
1	#include "TROOT.h"
2	#include "TVirtualMC.h"
3	#include "TGeant4.h"
4	#include "TString.h"
5	#include "PamVMCApplication.h"
6
7	#include <iostream>
8	// $Id: run_g4.C,v 1.0 2007/06/01
9	//
10	// Macro for running G4PAM_VMC with Geant4.
11
12	using std::cout;
13	using std::endl;
14
15	//void run_g4(const TString& configMacro = "g4Config.C")
16	int main(){
17	TROOT nik("Nikolay","Nikonov");
18
19	TString configMacro = "g4Config.C";
20	// Load basic libraries
21	gROOT->LoadMacro("./macro/basiclibs.C");
22	basiclibs();
23
24	// Load PAMELA libs
25	TString PAMLIB=gSystem->Getenv("PAM_LIB");
26	gSystem->Load(PAMLIB+"/libyoda.so");
27	gSystem->Load(PAMLIB+"/libDarthVader.so");
28	gSystem->Load(PAMLIB+"/libPamLevel2.so");
29
30	// Load Geant4 libraries
31
32	gROOT->LoadMacro("./macro/PamG4Libs.C");
33	PamG4Libs();
34
35	// Load this example libraries
36	TString G4WD=gSystem->Getenv("G4WORKDIR");
37	TString PLAT=gSystem->Getenv("PLATFORM");
38	gSystem->Load(G4WD+"/lib/tgt_"+PLAT+"/libPamVMC");
39
40	// MC application
41	PamVMCApplication* appl
42	= new PamVMCApplication("PAMG4_VMC", "PAMELA GEANT4 VMC application");
43
44	appl->GetPrimaryGenerator()->SetNofPrimaries(1);
45	appl->SetVerboseLevel(1);
46
47	// Macro to run with primary = polarized optical photon
48	//gROOT->LoadMacro("optPhoton.C");
49	//optPhoton();
50
51	// Initialize MC
52	appl->InitMC(configMacro);
53
54	// Customise Geant4 setting after initialization:
55	// Physics list
56	((TGeant4*)gMC)->ProcessGeantMacro("g4config2.in");
57	// Visualization settings
58	// ((TGeant4*)gMC)->ProcessGeantMacro("g4vis.in");
59
60	// Run MC
61	appl->RunMC(5000);
62	appl->FinishRun();
63	delete appl;
64	return 0;
65	}