1 |
|
---|
2 | // Geant4:
|
---|
3 | #include <G4ios.hh>
|
---|
4 | #include <G4RunManager.hh>
|
---|
5 | #include <G4UImanager.hh>
|
---|
6 | #include <G4SDManager.hh>
|
---|
7 |
|
---|
8 | //AIDA :
|
---|
9 | #include <AIDA/IAnalysisFactory.h>
|
---|
10 |
|
---|
11 | // std::
|
---|
12 | #include <iostream>
|
---|
13 |
|
---|
14 | // Slash :
|
---|
15 | #include <Slash/Core/IManager.h>
|
---|
16 | #include <Slash/Core/ISession.h>
|
---|
17 | #include <Slash/Data/IProcessor.h>
|
---|
18 |
|
---|
19 | // Lib :
|
---|
20 | #include <Lib/Manager.h>
|
---|
21 | #include <Lib/System.h>
|
---|
22 | #include <Lib/Out.h>
|
---|
23 | #include <Lib/Cast.h>
|
---|
24 |
|
---|
25 | // G4Lab :
|
---|
26 | #include <G4Lab/UIOnX.h>
|
---|
27 | #include <G4Lab/TrackingAction.h>
|
---|
28 | #include <G4Lab/Interfaces/IGeant4Manager.h>
|
---|
29 | //#include <G4Lab/DigitsCollectionAccessor.h>
|
---|
30 |
|
---|
31 | //ELYSE:
|
---|
32 | #include "ELYSE/Analysis.hh"
|
---|
33 | #include "ELYSE/SteppingVerbose.hh"
|
---|
34 | #include "ELYSE/DetectorConstruction.hh"
|
---|
35 | #include "ELYSE/PhysicsList.hh"
|
---|
36 | #include "ELYSE/PrimaryGeneratorAction.hh"
|
---|
37 | #include "ELYSE/EventAction.hh"
|
---|
38 | #include "ELYSE/RunAction.hh"
|
---|
39 | #include "ELYSE/StackingAction.hh"
|
---|
40 | #include "ELYSE/TrackingAction.hh"
|
---|
41 | #include "ELYSE/SteppingAction.hh"
|
---|
42 | #include <ELYSE/HitsCollectionAccessor.h>
|
---|
43 |
|
---|
44 | // Create a manager to access the ELYSEAnalysis from the callbacks :
|
---|
45 |
|
---|
46 | //JEC 10/1/06 introduce ELYSE namespace
|
---|
47 | namespace ELYSE {
|
---|
48 |
|
---|
49 | class AppManager
|
---|
50 | :public virtual Slash::Core::IManager
|
---|
51 | ,public Analysis {
|
---|
52 | public: //Slash::Core::IManager
|
---|
53 | virtual std::string name() const { return fName;}
|
---|
54 | virtual void* cast(const std::string& aClass) const {
|
---|
55 | if_Lib_SCast(ELYSE::AppManager)
|
---|
56 | else if_Lib_SCast(ELYSE::IAppManager)
|
---|
57 | else if_Lib_SCast(ELYSE::Analysis)
|
---|
58 | else if_Lib_SCast(Slash::Core::IManager)
|
---|
59 | else return 0;
|
---|
60 | }
|
---|
61 | public://IAppManager
|
---|
62 | virtual bool initialize() {
|
---|
63 | //GB :
|
---|
64 | // Used in the ELYSE_Initialize callback.
|
---|
65 | // Do all the below here because of platforms having
|
---|
66 | // not Geant4 shared libs.
|
---|
67 | IGeant4Manager* g4Manager =
|
---|
68 | Lib_findManager(fSession,"Geant4Manager",IGeant4Manager);
|
---|
69 | if(!g4Manager) {
|
---|
70 | Lib::Out out(fSession.printer());
|
---|
71 | out << "ELYSE_Initialize :"
|
---|
72 | << " Geant4Manager not found."
|
---|
73 | << Lib::endl;
|
---|
74 | return false;
|
---|
75 | }
|
---|
76 |
|
---|
77 | // To initialize G4Lab Types (PV, Trajectory, HitsCollection, etc.. types).
|
---|
78 | // To initialize G4Lab SoNodes.
|
---|
79 | g4Manager->initialize();
|
---|
80 |
|
---|
81 | // Overload the HitsCollectionAccessor of G4Lab :
|
---|
82 | Slash::Data::IProcessor* accessorManager =
|
---|
83 | Lib_findManager(fSession,"AccessorManager",Slash::Data::IProcessor);
|
---|
84 | if(accessorManager) {
|
---|
85 | G4SDManager* sdManager = g4Manager->sdManager();
|
---|
86 | if(sdManager) {
|
---|
87 | G4HCtable* hcTable = sdManager->GetHCtable();
|
---|
88 | if(hcTable) {
|
---|
89 | int number = hcTable->entries();
|
---|
90 | for(int index=0;index<number;index++) {
|
---|
91 | G4String hcName = hcTable->GetHCname(index);
|
---|
92 | accessorManager->removeAccessor(hcName);
|
---|
93 | accessorManager->addAccessor
|
---|
94 | (new ELYSE::HitsCollectionAccessor(fSession,hcName));
|
---|
95 | }
|
---|
96 | }
|
---|
97 | }
|
---|
98 | }
|
---|
99 | return true;
|
---|
100 | }
|
---|
101 | public:
|
---|
102 | AppManager(Slash::Core::ISession& aSession,AIDA::IAnalysisFactory* aAIDA,bool aBatch)
|
---|
103 | :Analysis(aAIDA,"",aBatch)
|
---|
104 | ,fSession(aSession)
|
---|
105 | ,fName("ELYSE::AppManager") {}
|
---|
106 | virtual ~AppManager(){}
|
---|
107 | private:
|
---|
108 | Slash::Core::ISession& fSession;
|
---|
109 | std::string fName;
|
---|
110 | };
|
---|
111 | }
|
---|
112 |
|
---|
113 | //////////////////////////////////////////////////////////////////////////////
|
---|
114 | int main(
|
---|
115 | int aArgc
|
---|
116 | ,char** aArgv
|
---|
117 | )
|
---|
118 | //////////////////////////////////////////////////////////////////////////////
|
---|
119 | //!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!//
|
---|
120 | {
|
---|
121 | //Lib::Debug::checkByClass(true);
|
---|
122 |
|
---|
123 |
|
---|
124 | // User Verbose output class
|
---|
125 | //
|
---|
126 | G4VSteppingVerbose* verbosity = new ELYSE::SteppingVerbose;
|
---|
127 | G4VSteppingVerbose::SetInstance(verbosity);
|
---|
128 |
|
---|
129 |
|
---|
130 | G4RunManager* runManager = new G4RunManager;
|
---|
131 |
|
---|
132 | // Need to pass the G4RunManager at UIOnX creation
|
---|
133 | std::string gui = "$ELYSEROOT/scripts/OnX/ELYSE_session.onx";
|
---|
134 | G4Lab::UIOnX* session = new G4Lab::UIOnX(*runManager,gui,aArgc,aArgv);
|
---|
135 | if(!session->isValid()) {
|
---|
136 | std::cout << "ELYSE_session : problem starting OnX." << std::endl;
|
---|
137 | return EXIT_FAILURE;
|
---|
138 | }
|
---|
139 |
|
---|
140 | ELYSE::AppManager* appManager = 0;
|
---|
141 |
|
---|
142 | AIDA::IAnalysisFactory* aida = session->findAIDA();
|
---|
143 | if(!aida) {
|
---|
144 | std::cout << "ELYSE_session : AIDA not found." << std::endl;
|
---|
145 | }
|
---|
146 |
|
---|
147 | appManager = new ELYSE::AppManager(*session->session(),aida,false);
|
---|
148 | session->addManager(appManager);
|
---|
149 |
|
---|
150 | ELYSE::DetectorConstruction* ELYSEdetector = new ELYSE::DetectorConstruction();
|
---|
151 |
|
---|
152 | runManager->SetUserInitialization(ELYSEdetector);
|
---|
153 |
|
---|
154 | runManager->SetUserInitialization(new ELYSE::PhysicsList);
|
---|
155 |
|
---|
156 | G4UImanager* UI = G4UImanager::GetUIpointer();
|
---|
157 | if (UI) {
|
---|
158 | std::string file =
|
---|
159 | Lib::System::getenv("ELYSEROOT")+"/scripts/Geant4/jobOptions.mac";
|
---|
160 | UI->ApplyCommand("/control/execute "+file);
|
---|
161 | }
|
---|
162 |
|
---|
163 |
|
---|
164 | // Set user action classes
|
---|
165 | ELYSE::PrimaryGeneratorAction* myGeneratorAction = new ELYSE::PrimaryGeneratorAction();
|
---|
166 | runManager->SetUserAction(myGeneratorAction);
|
---|
167 |
|
---|
168 | ELYSE::RunAction* myRunAction = new ELYSE::RunAction();
|
---|
169 |
|
---|
170 | runManager->SetUserAction(myRunAction);
|
---|
171 |
|
---|
172 |
|
---|
173 | runManager->SetUserAction(new ELYSE::EventAction(*appManager,
|
---|
174 | myRunAction,
|
---|
175 | ELYSEdetector,
|
---|
176 | myGeneratorAction));
|
---|
177 | runManager->SetUserAction(new ELYSE::TrackingAction);
|
---|
178 |
|
---|
179 | runManager->SetUserAction(new ELYSE::StackingAction);
|
---|
180 |
|
---|
181 | runManager->SetUserAction(new ELYSE::SteppingAction);
|
---|
182 |
|
---|
183 | // Initialize G4 kernel
|
---|
184 | runManager->Initialize();
|
---|
185 |
|
---|
186 | if (UI) {
|
---|
187 | std::string file =
|
---|
188 | Lib::System::getenv("ELYSEROOT")+"/scripts/Geant4/init.g4m";
|
---|
189 | UI->ApplyCommand("/control/execute "+file);
|
---|
190 | }
|
---|
191 |
|
---|
192 |
|
---|
193 | // Main UI loop :
|
---|
194 | // GUI is really constructed here.
|
---|
195 | // Create callbacks are executed here too.
|
---|
196 | session->SessionStart();
|
---|
197 |
|
---|
198 | delete session; //It will delete the appManager.
|
---|
199 |
|
---|
200 | delete runManager;
|
---|
201 |
|
---|
202 | return EXIT_SUCCESS;
|
---|
203 | }
|
---|
204 |
|
---|