| Geant4 Cross Reference |
1 // 1 //
2 // ******************************************* 2 // ********************************************************************
3 // * License and Disclaimer 3 // * License and Disclaimer *
4 // * 4 // * *
5 // * The Geant4 software is copyright of th 5 // * The Geant4 software is copyright of the Copyright Holders of *
6 // * the Geant4 Collaboration. It is provided 6 // * the Geant4 Collaboration. It is provided under the terms and *
7 // * conditions of the Geant4 Software License 7 // * conditions of the Geant4 Software License, included in the file *
8 // * LICENSE and available at http://cern.ch/ 8 // * LICENSE and available at http://cern.ch/geant4/license . These *
9 // * include a list of copyright holders. 9 // * include a list of copyright holders. *
10 // * 10 // * *
11 // * Neither the authors of this software syst 11 // * Neither the authors of this software system, nor their employing *
12 // * institutes,nor the agencies providing fin 12 // * institutes,nor the agencies providing financial support for this *
13 // * work make any representation or warran 13 // * work make any representation or warranty, express or implied, *
14 // * regarding this software system or assum 14 // * regarding this software system or assume any liability for its *
15 // * use. Please see the license in the file 15 // * use. Please see the license in the file LICENSE and URL above *
16 // * for the full disclaimer and the limitatio 16 // * for the full disclaimer and the limitation of liability. *
17 // * 17 // * *
18 // * This code implementation is the result 18 // * This code implementation is the result of the scientific and *
19 // * technical work of the GEANT4 collaboratio 19 // * technical work of the GEANT4 collaboration. *
20 // * By using, copying, modifying or distri 20 // * By using, copying, modifying or distributing the software (or *
21 // * any work based on the software) you ag 21 // * any work based on the software) you agree to acknowledge its *
22 // * use in resulting scientific publicati 22 // * use in resulting scientific publications, and indicate your *
23 // * acceptance of all terms of the Geant4 Sof 23 // * acceptance of all terms of the Geant4 Software license. *
24 // ******************************************* 24 // ********************************************************************
25 // 25 //
26 /* 26 /*
27 * =========================================== 27 * ============================================================================
28 * 28 *
29 * Filename: CexmcRunManager.cc 29 * Filename: CexmcRunManager.cc
30 * 30 *
31 * Description: run manager 31 * Description: run manager
32 * 32 *
33 * Version: 1.0 33 * Version: 1.0
34 * Created: 03.11.2009 20:27:46 34 * Created: 03.11.2009 20:27:46
35 * Revision: none 35 * Revision: none
36 * Compiler: gcc 36 * Compiler: gcc
37 * 37 *
38 * Author: Alexey Radkov (), 38 * Author: Alexey Radkov (),
39 * Company: PNPI 39 * Company: PNPI
40 * 40 *
41 * =========================================== 41 * ============================================================================
42 */ 42 */
43 43
44 #include <stdlib.h> 44 #include <stdlib.h>
45 #include <sys/stat.h> 45 #include <sys/stat.h>
46 #include <vector> 46 #include <vector>
47 #include <fstream> 47 #include <fstream>
48 #ifdef CEXMC_USE_PERSISTENCY 48 #ifdef CEXMC_USE_PERSISTENCY
49 #include <boost/archive/binary_iarchive.hpp> 49 #include <boost/archive/binary_iarchive.hpp>
50 #include <boost/archive/binary_oarchive.hpp> 50 #include <boost/archive/binary_oarchive.hpp>
51 #endif 51 #endif
52 #include <G4Eta.hh> 52 #include <G4Eta.hh>
53 #include <G4DigiManager.hh> 53 #include <G4DigiManager.hh>
54 #include <G4SDManager.hh> 54 #include <G4SDManager.hh>
55 #include <G4DecayTable.hh> 55 #include <G4DecayTable.hh>
56 #include <G4ParticleDefinition.hh> 56 #include <G4ParticleDefinition.hh>
57 #include <G4ParticleTable.hh> 57 #include <G4ParticleTable.hh>
58 #include <G4UImanager.hh> 58 #include <G4UImanager.hh>
59 #include <G4Timer.hh> 59 #include <G4Timer.hh>
60 #include <G4Region.hh> 60 #include <G4Region.hh>
61 #include <G4RegionStore.hh> 61 #include <G4RegionStore.hh>
62 #include <G4ProductionCuts.hh> 62 #include <G4ProductionCuts.hh>
63 #include <G4VisManager.hh> 63 #include <G4VisManager.hh>
64 #include <G4Scene.hh> 64 #include <G4Scene.hh>
65 #include <G4VModel.hh> 65 #include <G4VModel.hh>
66 #include <G4Version.hh> 66 #include <G4Version.hh>
67 #include "CexmcRunManager.hh" 67 #include "CexmcRunManager.hh"
68 #include "CexmcRunManagerMessenger.hh" 68 #include "CexmcRunManagerMessenger.hh"
69 #include "CexmcRunAction.hh" 69 #include "CexmcRunAction.hh"
70 #include "CexmcRun.hh" 70 #include "CexmcRun.hh"
71 #include "CexmcPhysicsManager.hh" 71 #include "CexmcPhysicsManager.hh"
72 #include "CexmcProductionModel.hh" 72 #include "CexmcProductionModel.hh"
73 #include "CexmcSimpleDecayTableStore.hh" 73 #include "CexmcSimpleDecayTableStore.hh"
74 #include "CexmcPrimaryGeneratorAction.hh" 74 #include "CexmcPrimaryGeneratorAction.hh"
75 #include "CexmcEnergyDepositDigitizer.hh" 75 #include "CexmcEnergyDepositDigitizer.hh"
76 #include "CexmcSimpleEnergyDeposit.hh" 76 #include "CexmcSimpleEnergyDeposit.hh"
77 #include "CexmcEnergyDepositInLeftRightSet.hh" 77 #include "CexmcEnergyDepositInLeftRightSet.hh"
78 #include "CexmcEnergyDepositInCalorimeter.hh" 78 #include "CexmcEnergyDepositInCalorimeter.hh"
79 #include "CexmcTrackPoints.hh" 79 #include "CexmcTrackPoints.hh"
80 #include "CexmcTrackPointsInLeftRightSet.hh" 80 #include "CexmcTrackPointsInLeftRightSet.hh"
81 #include "CexmcTrackPointsInCalorimeter.hh" 81 #include "CexmcTrackPointsInCalorimeter.hh"
82 #include "CexmcChargeExchangeReconstructor.hh" 82 #include "CexmcChargeExchangeReconstructor.hh"
83 #include "CexmcEventAction.hh" 83 #include "CexmcEventAction.hh"
84 #include "CexmcParticleGun.hh" 84 #include "CexmcParticleGun.hh"
85 #include "CexmcEnergyDepositStore.hh" 85 #include "CexmcEnergyDepositStore.hh"
86 #include "CexmcTrackPointsStore.hh" 86 #include "CexmcTrackPointsStore.hh"
87 #include "CexmcSetup.hh" 87 #include "CexmcSetup.hh"
88 #include "CexmcEventSObject.hh" 88 #include "CexmcEventSObject.hh"
89 #include "CexmcEventFastSObject.hh" 89 #include "CexmcEventFastSObject.hh"
90 #include "CexmcTrackPointInfo.hh" 90 #include "CexmcTrackPointInfo.hh"
91 #include "CexmcEventInfo.hh" 91 #include "CexmcEventInfo.hh"
92 #include "CexmcBasicPhysicsSettings.hh" 92 #include "CexmcBasicPhysicsSettings.hh"
93 #include "CexmcSensitiveDetectorsAttributes.hh 93 #include "CexmcSensitiveDetectorsAttributes.hh"
94 #include "CexmcCustomFilterEval.hh" 94 #include "CexmcCustomFilterEval.hh"
95 #include "CexmcScenePrimitives.hh" 95 #include "CexmcScenePrimitives.hh"
96 96
97 97
98 namespace 98 namespace
99 { 99 {
100 G4String gdmlFileExtension( ".gdml" ); 100 G4String gdmlFileExtension( ".gdml" );
101 G4String gdmlbz2FileExtension( ".gdml.bz2 101 G4String gdmlbz2FileExtension( ".gdml.bz2" );
102 } 102 }
103 103
104 104
105 CexmcRunManager::CexmcRunManager( const G4Stri << 105 CexmcRunManager::CexmcRunManager( const G4String & projectId,
106 const G4Stri << 106 const G4String & rProject,
107 G4bool over 107 G4bool overrideExistingProject ) :
108 basePhysicsUsed( CexmcPMFactoryInstance::G 108 basePhysicsUsed( CexmcPMFactoryInstance::GetBasePhysics() ),
109 productionModelType( CexmcUnknownProductio 109 productionModelType( CexmcUnknownProductionModel ),
110 gdmlFileName( "default.gdml" ), shouldGdml << 110 gdmlFileName( "default.gdml" ), shouldGdmlFileBeValidated( true ),
111 zipGdmlFile( false ), projectsDir( "." ), << 111 zipGdmlFile( false ), projectsDir( "." ), projectId( projectId ),
112 rProject( rProject_ ), guiMacroName( "" ), << 112 rProject( rProject ), guiMacroName( "" ), cfFileName( "" ),
113 eventCountPolicy( CexmcCountAllEvents ), << 113 eventCountPolicy( CexmcCountAllEvents ), areLiveHistogramsEnabled( false ),
114 skipInteractionsWithoutEDTonWrite( true ), 114 skipInteractionsWithoutEDTonWrite( true ),
115 evDataVerboseLevel( CexmcWriteEventDataOnE 115 evDataVerboseLevel( CexmcWriteEventDataOnEveryEDT ),
116 rEvDataVerboseLevel( CexmcWriteNoEventData 116 rEvDataVerboseLevel( CexmcWriteNoEventData ), numberOfEventsProcessed( 0 ),
117 numberOfEventsProcessedEffective( 0 ), cur 117 numberOfEventsProcessedEffective( 0 ), curEventRead( 0 ),
118 #ifdef CEXMC_USE_PERSISTENCY 118 #ifdef CEXMC_USE_PERSISTENCY
119 eventsArchive( NULL ), fastEventsArchive( 119 eventsArchive( NULL ), fastEventsArchive( NULL ),
120 #ifdef CEXMC_USE_CUSTOM_FILTER 120 #ifdef CEXMC_USE_CUSTOM_FILTER
121 customFilter( NULL ), 121 customFilter( NULL ),
122 #endif 122 #endif
123 #endif 123 #endif
124 physicsManager( NULL ), messenger( NULL ) 124 physicsManager( NULL ), messenger( NULL )
125 { 125 {
126 /* this exception must be caught before cr 126 /* this exception must be caught before creating the object! */
127 if ( rProject != "" && rProject == project 127 if ( rProject != "" && rProject == projectId )
128 throw CexmcException( CexmcWeirdExcept 128 throw CexmcException( CexmcWeirdException );
129 129
130 const char * projectsDirEnv( std::getenv( << 130 const char * projectsDirEnv( getenv( "CEXMC_PROJECTS_DIR" ) );
131 131
132 if ( projectsDirEnv ) 132 if ( projectsDirEnv )
133 projectsDir = projectsDirEnv; 133 projectsDir = projectsDirEnv;
134 134
135 struct stat tmp; 135 struct stat tmp;
136 if ( ProjectIsSaved() && 136 if ( ProjectIsSaved() &&
137 stat( ( projectsDir + "/" + projectId 137 stat( ( projectsDir + "/" + projectId + ".rdb" ).c_str(), &tmp ) == 0
138 && ! overrideExistingProject ) 138 && ! overrideExistingProject )
139 throw CexmcException( CexmcProjectExis 139 throw CexmcException( CexmcProjectExists );
140 140
141 messenger = new CexmcRunManagerMessenger( 141 messenger = new CexmcRunManagerMessenger( this );
142 142
143 #ifdef CEXMC_USE_PERSISTENCY 143 #ifdef CEXMC_USE_PERSISTENCY
144 if ( ProjectIsRead() ) 144 if ( ProjectIsRead() )
145 ReadPreinitProjectData(); 145 ReadPreinitProjectData();
146 #endif 146 #endif
147 } 147 }
148 148
149 149
150 CexmcRunManager::~CexmcRunManager() 150 CexmcRunManager::~CexmcRunManager()
151 { 151 {
152 #ifdef CEXMC_USE_CUSTOM_FILTER 152 #ifdef CEXMC_USE_CUSTOM_FILTER
153 delete customFilter; 153 delete customFilter;
154 #endif 154 #endif
155 delete messenger; 155 delete messenger;
156 } 156 }
157 157
158 158
159 #ifdef CEXMC_USE_PERSISTENCY 159 #ifdef CEXMC_USE_PERSISTENCY
160 160
161 void CexmcRunManager::ReadPreinitProjectData( 161 void CexmcRunManager::ReadPreinitProjectData( void )
162 { 162 {
163 if ( ! ProjectIsRead() ) 163 if ( ! ProjectIsRead() )
164 return; 164 return;
165 165
166 /* read run data */ 166 /* read run data */
167 std::ifstream runDataFile( ( projectsDir 167 std::ifstream runDataFile( ( projectsDir + "/" + rProject + ".rdb" ).
168 c_str() ); 168 c_str() );
169 if ( ! runDataFile ) 169 if ( ! runDataFile )
170 throw CexmcException( CexmcReadProject 170 throw CexmcException( CexmcReadProjectIncomplete );
171 171
172 { 172 {
173 boost::archive::binary_iarchive archi 173 boost::archive::binary_iarchive archive( runDataFile );
174 archive >> sObject; 174 archive >> sObject;
175 } 175 }
176 176
177 basePhysicsUsed = sObject.basePhysicsUsed; 177 basePhysicsUsed = sObject.basePhysicsUsed;
178 178
179 productionModelType = sObject.productionMo 179 productionModelType = sObject.productionModelType;
180 180
181 /* read gdml file */ 181 /* read gdml file */
182 G4String cmd; 182 G4String cmd;
183 if ( ProjectIsSaved() ) 183 if ( ProjectIsSaved() )
184 { 184 {
185 G4String fileExtension( zipGdmlFile ? 185 G4String fileExtension( zipGdmlFile ? gdmlbz2FileExtension :
186 186 gdmlFileExtension );
187 cmd = G4String( "cp " ) + projectsDir 187 cmd = G4String( "cp " ) + projectsDir + "/" + rProject + fileExtension +
188 " " + projectsDir 188 " " + projectsDir + "/" + projectId + fileExtension;
189 if ( system( cmd ) != 0 ) 189 if ( system( cmd ) != 0 )
190 throw CexmcException( CexmcReadPro 190 throw CexmcException( CexmcReadProjectIncomplete );
191 } 191 }
192 192
193 if ( zipGdmlFile ) 193 if ( zipGdmlFile )
194 { 194 {
195 cmd = G4String( "bunzip2 " ) + project 195 cmd = G4String( "bunzip2 " ) + projectsDir + "/" + rProject +
196 196 gdmlbz2FileExtension;
197 if ( system( cmd ) != 0 ) 197 if ( system( cmd ) != 0 )
198 throw CexmcException( CexmcFileCom 198 throw CexmcException( CexmcFileCompressException );
199 } 199 }
200 200
201 gdmlFileName = projectsDir + "/" + rProjec 201 gdmlFileName = projectsDir + "/" + rProject + gdmlFileExtension;
202 } 202 }
203 203
204 204
205 void CexmcRunManager::ReadProject( void ) 205 void CexmcRunManager::ReadProject( void )
206 { 206 {
207 if ( ! ProjectIsRead() ) 207 if ( ! ProjectIsRead() )
208 return; 208 return;
209 209
210 if ( ! physicsManager ) 210 if ( ! physicsManager )
211 throw CexmcException( CexmcWeirdExcept 211 throw CexmcException( CexmcWeirdException );
212 212
213 physicsManager->GetProductionModel()->SetA 213 physicsManager->GetProductionModel()->SetAngularRanges(
214 214 sObject.angularRanges );
215 G4DecayTable * etaDecayTable( G4Eta::Defi 215 G4DecayTable * etaDecayTable( G4Eta::Definition()->GetDecayTable() );
216 for ( CexmcDecayBranchesStore::const_itera 216 for ( CexmcDecayBranchesStore::const_iterator
217 k( sObject.etaDecayTable.GetDecayB 217 k( sObject.etaDecayTable.GetDecayBranches().begin() );
218 k != sObject.etaDecayTable.GetDeca 218 k != sObject.etaDecayTable.GetDecayBranches().end(); ++k )
219 { 219 {
220 etaDecayTable->GetDecayChannel( k->fir 220 etaDecayTable->GetDecayChannel( k->first )->SetBR( k->second );
221 } 221 }
222 222
223 physicsManager->GetProductionModel()->Appl 223 physicsManager->GetProductionModel()->ApplyFermiMotion(
224 sO 224 sObject.fermiMotionIsOn, false );
225 eventCountPolicy = sObject.eventCountPolic 225 eventCountPolicy = sObject.eventCountPolicy;
226 226
227 G4Region * region( G4RegionStore::GetInst 227 G4Region * region( G4RegionStore::GetInstance()->GetRegion(
228 228 CexmcCalorimeterRegionName ) );
229 if ( ! region ) 229 if ( ! region )
230 throw CexmcException( CexmcCalorimeter 230 throw CexmcException( CexmcCalorimeterRegionNotInitialized );
231 231
232 region->GetProductionCuts()->SetProduction 232 region->GetProductionCuts()->SetProductionCuts(
233 233 sObject.calorimeterRegCuts );
234 234
235 const CexmcPrimaryGeneratorAction * prima 235 const CexmcPrimaryGeneratorAction * primaryGeneratorAction(
236 static_cast< const 236 static_cast< const CexmcPrimaryGeneratorAction * >(
237 237 userPrimaryGeneratorAction ) );
238 CexmcPrimaryGeneratorAction * thePr 238 CexmcPrimaryGeneratorAction * thePrimaryGeneratorAction(
239 const_cast< CexmcP 239 const_cast< CexmcPrimaryGeneratorAction * >(
240 pr 240 primaryGeneratorAction ) );
241 CexmcParticleGun * particleGun( 241 CexmcParticleGun * particleGun(
242 thePrimaryGene 242 thePrimaryGeneratorAction->GetParticleGun() );
243 G4ParticleDefinition * particleDefinition 243 G4ParticleDefinition * particleDefinition(
244 G4ParticleTable::GetPartic 244 G4ParticleTable::GetParticleTable()->FindParticle(
245 245 sObject.beamParticle ) );
246 if ( ! particleDefinition ) 246 if ( ! particleDefinition )
247 throw CexmcException( CexmcWeirdExcept 247 throw CexmcException( CexmcWeirdException );
248 248
249 particleGun->SetParticleDefinition( partic 249 particleGun->SetParticleDefinition( particleDefinition );
250 particleGun->SetOrigPosition( sObject.beam 250 particleGun->SetOrigPosition( sObject.beamPos, false );
251 particleGun->SetOrigDirection( sObject.bea 251 particleGun->SetOrigDirection( sObject.beamDir, false );
252 particleGun->SetOrigMomentumAmp( sObject.b 252 particleGun->SetOrigMomentumAmp( sObject.beamMomentumAmp, false );
253 253
254 thePrimaryGeneratorAction->SetFwhmPosX( sO 254 thePrimaryGeneratorAction->SetFwhmPosX( sObject.beamFwhmPosX, false );
255 thePrimaryGeneratorAction->SetFwhmPosY( sO 255 thePrimaryGeneratorAction->SetFwhmPosY( sObject.beamFwhmPosY, false );
256 thePrimaryGeneratorAction->SetFwhmDirX( sO 256 thePrimaryGeneratorAction->SetFwhmDirX( sObject.beamFwhmDirX, false );
257 thePrimaryGeneratorAction->SetFwhmDirY( sO 257 thePrimaryGeneratorAction->SetFwhmDirY( sObject.beamFwhmDirY, false );
258 thePrimaryGeneratorAction->SetFwhmMomentum 258 thePrimaryGeneratorAction->SetFwhmMomentumAmp(
259 sObjec 259 sObject.beamFwhmMomentumAmp, false );
260 260
261 G4DigiManager * digiManager 261 G4DigiManager * digiManager( G4DigiManager::GetDMpointer() );
262 CexmcEnergyDepositDigitizer * edDigitizer 262 CexmcEnergyDepositDigitizer * edDigitizer(
263 static_cast< CexmcEnergyDepositDig 263 static_cast< CexmcEnergyDepositDigitizer * >(
264 digiManager->FindDigitizerModu 264 digiManager->FindDigitizerModule( CexmcEDDigitizerName ) ) );
265 if ( ! edDigitizer ) 265 if ( ! edDigitizer )
266 throw CexmcException( CexmcWeirdExcept 266 throw CexmcException( CexmcWeirdException );
267 267
268 edDigitizer->SetMonitorThreshold( sObject. 268 edDigitizer->SetMonitorThreshold( sObject.monitorEDThreshold, false );
269 edDigitizer->SetVetoCounterLeftThreshold( 269 edDigitizer->SetVetoCounterLeftThreshold(
270 sObject.vetoCo 270 sObject.vetoCounterEDLeftThreshold, false );
271 edDigitizer->SetVetoCounterRightThreshold( 271 edDigitizer->SetVetoCounterRightThreshold(
272 sObject.vetoCo 272 sObject.vetoCounterEDRightThreshold, false );
273 edDigitizer->SetCalorimeterLeftThreshold( 273 edDigitizer->SetCalorimeterLeftThreshold(
274 sObject.calori 274 sObject.calorimeterEDLeftThreshold, false );
275 edDigitizer->SetCalorimeterRightThreshold( 275 edDigitizer->SetCalorimeterRightThreshold(
276 sObject.calori 276 sObject.calorimeterEDRightThreshold, false );
277 edDigitizer->SetCalorimeterTriggerAlgorith 277 edDigitizer->SetCalorimeterTriggerAlgorithm(
278 sObject.calori 278 sObject.calorimeterTriggerAlgorithm, false );
279 edDigitizer->SetOuterCrystalsVetoAlgorithm 279 edDigitizer->SetOuterCrystalsVetoAlgorithm(
280 sObject.outerC 280 sObject.outerCrystalsVetoAlgorithm, false );
281 edDigitizer->SetOuterCrystalsVetoFraction( 281 edDigitizer->SetOuterCrystalsVetoFraction(
282 sObject.outerC 282 sObject.outerCrystalsVetoFraction, false );
283 edDigitizer->ApplyFiniteCrystalResolution( 283 edDigitizer->ApplyFiniteCrystalResolution(
284 sObject.applyF 284 sObject.applyFiniteCrystalResolution, false );
285 edDigitizer->SetCrystalResolutionData( sOb 285 edDigitizer->SetCrystalResolutionData( sObject.crystalResolutionData );
286 286
287 const CexmcEventAction * eventAction( 287 const CexmcEventAction * eventAction(
288 static_cast< const CexmcEventA 288 static_cast< const CexmcEventAction * >( userEventAction ) );
289 if ( ! eventAction ) 289 if ( ! eventAction )
290 throw CexmcException( CexmcWeirdExcept 290 throw CexmcException( CexmcWeirdException );
291 291
292 CexmcEventAction * theEventAction( 292 CexmcEventAction * theEventAction( const_cast< CexmcEventAction * >(
293 293 eventAction ) );
294 CexmcChargeExchangeReconstructor * recons 294 CexmcChargeExchangeReconstructor * reconstructor(
295 theEve 295 theEventAction->GetReconstructor() );
296 if ( ! reconstructor ) 296 if ( ! reconstructor )
297 throw CexmcException( CexmcWeirdExcept 297 throw CexmcException( CexmcWeirdException );
298 298
299 reconstructor->SetCalorimeterEntryPointDef 299 reconstructor->SetCalorimeterEntryPointDefinitionAlgorithm(
300 sObjec 300 sObject.epDefinitionAlgorithm );
301 reconstructor->SetCalorimeterEntryPointDep 301 reconstructor->SetCalorimeterEntryPointDepthDefinitionAlgorithm(
302 sObjec 302 sObject.epDepthDefinitionAlgorithm );
303 reconstructor->SetCrystalSelectionAlgorith 303 reconstructor->SetCrystalSelectionAlgorithm( sObject.csAlgorithm );
304 reconstructor->UseInnerRefCrystal( sObject 304 reconstructor->UseInnerRefCrystal( sObject.useInnerRefCrystal );
305 reconstructor->SetCalorimeterEntryPointDep 305 reconstructor->SetCalorimeterEntryPointDepth( sObject.epDepth );
306 reconstructor->UseTableMass( sObject.useTa 306 reconstructor->UseTableMass( sObject.useTableMass );
307 reconstructor->UseMassCut( sObject.useMass 307 reconstructor->UseMassCut( sObject.useMassCut );
308 reconstructor->SetMassCutOPCenter( sObject 308 reconstructor->SetMassCutOPCenter( sObject.mCutOPCenter );
309 reconstructor->SetMassCutNOPCenter( sObjec 309 reconstructor->SetMassCutNOPCenter( sObject.mCutNOPCenter );
310 reconstructor->SetMassCutOPWidth( sObject. 310 reconstructor->SetMassCutOPWidth( sObject.mCutOPWidth );
311 reconstructor->SetMassCutNOPWidth( sObject 311 reconstructor->SetMassCutNOPWidth( sObject.mCutNOPWidth );
312 reconstructor->SetMassCutEllipseAngle( sOb 312 reconstructor->SetMassCutEllipseAngle( sObject.mCutAngle );
313 reconstructor->UseAbsorbedEnergyCut( sObje 313 reconstructor->UseAbsorbedEnergyCut( sObject.useAbsorbedEnergyCut );
314 reconstructor->SetAbsorbedEnergyCutCLCente 314 reconstructor->SetAbsorbedEnergyCutCLCenter( sObject.aeCutCLCenter );
315 reconstructor->SetAbsorbedEnergyCutCRCente 315 reconstructor->SetAbsorbedEnergyCutCRCenter( sObject.aeCutCRCenter );
316 reconstructor->SetAbsorbedEnergyCutCLWidth 316 reconstructor->SetAbsorbedEnergyCutCLWidth( sObject.aeCutCLWidth );
317 reconstructor->SetAbsorbedEnergyCutCRWidth 317 reconstructor->SetAbsorbedEnergyCutCRWidth( sObject.aeCutCRWidth );
318 reconstructor->SetAbsorbedEnergyCutEllipse 318 reconstructor->SetAbsorbedEnergyCutEllipseAngle( sObject.aeCutAngle );
319 reconstructor->SetExpectedMomentumAmp( sOb 319 reconstructor->SetExpectedMomentumAmp( sObject.expectedMomentumAmp );
320 reconstructor->SetEDCollectionAlgorithm( s 320 reconstructor->SetEDCollectionAlgorithm( sObject.edCollectionAlgorithm );
321 321
322 physicsManager->SetProposedMaxIL( sObject. 322 physicsManager->SetProposedMaxIL( sObject.proposedMaxIL );
323 323
324 rEvDataVerboseLevel = sObject.evDataVerbos 324 rEvDataVerboseLevel = sObject.evDataVerboseLevel;
325 evDataVerboseLevel = rEvDataVerboseLevel; 325 evDataVerboseLevel = rEvDataVerboseLevel;
326 } 326 }
327 327
328 328
329 void CexmcRunManager::SaveProject( void ) 329 void CexmcRunManager::SaveProject( void )
330 { 330 {
331 if ( ! ProjectIsSaved() ) 331 if ( ! ProjectIsSaved() )
332 return; 332 return;
333 333
334 /* save run data */ 334 /* save run data */
335 if ( ! physicsManager ) 335 if ( ! physicsManager )
336 throw CexmcException( CexmcWeirdExcept 336 throw CexmcException( CexmcWeirdException );
337 337
338 CexmcSimpleDecayTableStore etaDecayTable( 338 CexmcSimpleDecayTableStore etaDecayTable(
339 G4Eta::Def 339 G4Eta::Definition()->GetDecayTable() );
340 const CexmcPrimaryGeneratorAction * prima 340 const CexmcPrimaryGeneratorAction * primaryGeneratorAction(
341 static_cast< const 341 static_cast< const CexmcPrimaryGeneratorAction * >(
342 342 userPrimaryGeneratorAction ) );
343 CexmcPrimaryGeneratorAction * thePrimaryG 343 CexmcPrimaryGeneratorAction * thePrimaryGeneratorAction(
344 const_cast< CexmcP 344 const_cast< CexmcPrimaryGeneratorAction * >(
345 pr 345 primaryGeneratorAction ) );
346 CexmcParticleGun * particleGun( 346 CexmcParticleGun * particleGun(
347 thePrimaryGenerato 347 thePrimaryGeneratorAction->GetParticleGun() );
348 348
349 G4DigiManager * digiManager 349 G4DigiManager * digiManager( G4DigiManager::GetDMpointer() );
350 CexmcEnergyDepositDigitizer * edDigitizer 350 CexmcEnergyDepositDigitizer * edDigitizer(
351 static_cast< CexmcEnergyDepositDig 351 static_cast< CexmcEnergyDepositDigitizer * >(
352 digiManager->FindDigitizerModu 352 digiManager->FindDigitizerModule( CexmcEDDigitizerName ) ) );
353 if ( ! edDigitizer ) 353 if ( ! edDigitizer )
354 throw CexmcException( CexmcWeirdExcept 354 throw CexmcException( CexmcWeirdException );
355 355
356 const CexmcEventAction * eventAction( 356 const CexmcEventAction * eventAction(
357 static_cast< const CexmcEventA 357 static_cast< const CexmcEventAction * >( userEventAction ) );
358 CexmcEventAction * theEventAction( 358 CexmcEventAction * theEventAction( const_cast< CexmcEventAction * >(
359 359 eventAction ) );
360 CexmcChargeExchangeReconstructor * recons 360 CexmcChargeExchangeReconstructor * reconstructor(
361 theEve 361 theEventAction->GetReconstructor() );
362 362
363 std::vector< G4double > calorimeterRegCut 363 std::vector< G4double > calorimeterRegCuts( 4 );
364 if ( ProjectIsRead() ) 364 if ( ProjectIsRead() )
365 { 365 {
366 calorimeterRegCuts = sObject.calorimet 366 calorimeterRegCuts = sObject.calorimeterRegCuts;
367 } 367 }
368 else 368 else
369 { 369 {
370 G4Region * region( G4RegionStore::Get 370 G4Region * region( G4RegionStore::GetInstance()->GetRegion(
371 371 CexmcCalorimeterRegionName ) );
372 if ( ! region ) 372 if ( ! region )
373 throw CexmcException( CexmcCalorim 373 throw CexmcException( CexmcCalorimeterRegionNotInitialized );
374 374
375 calorimeterRegCuts = region->GetProduc 375 calorimeterRegCuts = region->GetProductionCuts()->GetProductionCuts();
376 } 376 }
377 377
378 CexmcNmbOfHitsInRanges nmbOfHitsSampled; 378 CexmcNmbOfHitsInRanges nmbOfHitsSampled;
379 CexmcNmbOfHitsInRanges nmbOfHitsSampledFu 379 CexmcNmbOfHitsInRanges nmbOfHitsSampledFull;
380 CexmcNmbOfHitsInRanges nmbOfHitsTriggered 380 CexmcNmbOfHitsInRanges nmbOfHitsTriggeredRealRange;
381 CexmcNmbOfHitsInRanges nmbOfHitsTriggered 381 CexmcNmbOfHitsInRanges nmbOfHitsTriggeredRecRange;
382 CexmcNmbOfHitsInRanges nmbOfOrphanHits; 382 CexmcNmbOfHitsInRanges nmbOfOrphanHits;
383 G4int nmbOfFalseHitsTrig 383 G4int nmbOfFalseHitsTriggeredEDT( 0 );
384 G4int nmbOfFalseHitsTrig 384 G4int nmbOfFalseHitsTriggeredRec( 0 );
385 G4int nmbOfSavedEvents( 385 G4int nmbOfSavedEvents( 0 );
386 G4int nmbOfSavedFastEven 386 G4int nmbOfSavedFastEvents( 0 );
387 CexmcRun * run( static_cast< 387 CexmcRun * run( static_cast< CexmcRun * >( currentRun ) );
388 388
389 if ( run ) 389 if ( run )
390 { 390 {
391 nmbOfHitsSampled = run->GetNmbOfHitsSa 391 nmbOfHitsSampled = run->GetNmbOfHitsSampled();
392 nmbOfHitsSampledFull = run->GetNmbOfHi 392 nmbOfHitsSampledFull = run->GetNmbOfHitsSampledFull();
393 nmbOfHitsTriggeredRealRange = run->Get 393 nmbOfHitsTriggeredRealRange = run->GetNmbOfHitsTriggeredRealRange();
394 nmbOfHitsTriggeredRecRange = run->GetN 394 nmbOfHitsTriggeredRecRange = run->GetNmbOfHitsTriggeredRecRange();
395 nmbOfOrphanHits = run->GetNmbOfOrphanH 395 nmbOfOrphanHits = run->GetNmbOfOrphanHits();
396 nmbOfFalseHitsTriggeredEDT = run->GetN 396 nmbOfFalseHitsTriggeredEDT = run->GetNmbOfFalseHitsTriggeredEDT();
397 nmbOfFalseHitsTriggeredRec = run->GetN 397 nmbOfFalseHitsTriggeredRec = run->GetNmbOfFalseHitsTriggeredRec();
398 nmbOfSavedEvents = run->GetNmbOfSavedE 398 nmbOfSavedEvents = run->GetNmbOfSavedEvents();
399 nmbOfSavedFastEvents = run->GetNmbOfSa 399 nmbOfSavedFastEvents = run->GetNmbOfSavedFastEvents();
400 } 400 }
401 401
402 CexmcRunSObject sObjectToWrite = { 402 CexmcRunSObject sObjectToWrite = {
403 basePhysicsUsed, productionModelType, 403 basePhysicsUsed, productionModelType, gdmlFileName, etaDecayTable,
404 physicsManager->GetProductionModel()-> 404 physicsManager->GetProductionModel()->GetAngularRanges(),
405 physicsManager->GetProductionModel()-> 405 physicsManager->GetProductionModel()->IsFermiMotionOn(),
406 calorimeterRegCuts, eventCountPolicy, 406 calorimeterRegCuts, eventCountPolicy,
407 particleGun->GetParticleDefinition()-> 407 particleGun->GetParticleDefinition()->GetParticleName(),
408 particleGun->GetOrigPosition(), partic 408 particleGun->GetOrigPosition(), particleGun->GetOrigDirection(),
409 particleGun->GetOrigMomentumAmp(), 409 particleGun->GetOrigMomentumAmp(),
410 primaryGeneratorAction->GetFwhmPosX(), 410 primaryGeneratorAction->GetFwhmPosX(),
411 primaryGeneratorAction->GetFwhmPosY(), 411 primaryGeneratorAction->GetFwhmPosY(),
412 primaryGeneratorAction->GetFwhmDirX(), 412 primaryGeneratorAction->GetFwhmDirX(),
413 primaryGeneratorAction->GetFwhmDirY(), 413 primaryGeneratorAction->GetFwhmDirY(),
414 primaryGeneratorAction->GetFwhmMomentu 414 primaryGeneratorAction->GetFwhmMomentumAmp(),
415 edDigitizer->GetMonitorThreshold(), 415 edDigitizer->GetMonitorThreshold(),
416 edDigitizer->GetVetoCounterLeftThresho 416 edDigitizer->GetVetoCounterLeftThreshold(),
417 edDigitizer->GetVetoCounterRightThresh 417 edDigitizer->GetVetoCounterRightThreshold(),
418 edDigitizer->GetCalorimeterLeftThresho 418 edDigitizer->GetCalorimeterLeftThreshold(),
419 edDigitizer->GetCalorimeterRightThresh 419 edDigitizer->GetCalorimeterRightThreshold(),
420 edDigitizer->GetCalorimeterTriggerAlgo 420 edDigitizer->GetCalorimeterTriggerAlgorithm(),
421 edDigitizer->GetOuterCrystalsVetoAlgor 421 edDigitizer->GetOuterCrystalsVetoAlgorithm(),
422 edDigitizer->GetOuterCrystalsVetoFract 422 edDigitizer->GetOuterCrystalsVetoFraction(),
423 edDigitizer->IsFiniteCrystalResolution 423 edDigitizer->IsFiniteCrystalResolutionApplied(),
424 edDigitizer->GetCrystalResolutionData( 424 edDigitizer->GetCrystalResolutionData(),
425 reconstructor->GetCalorimeterEntryPoin 425 reconstructor->GetCalorimeterEntryPointDefinitionAlgorithm(),
426 reconstructor->GetCalorimeterEntryPoin 426 reconstructor->GetCalorimeterEntryPointDepthDefinitionAlgorithm(),
427 reconstructor->GetCrystalSelectionAlgo 427 reconstructor->GetCrystalSelectionAlgorithm(),
428 reconstructor->IsInnerRefCrystalUsed() 428 reconstructor->IsInnerRefCrystalUsed(),
429 reconstructor->GetCalorimeterEntryPoin 429 reconstructor->GetCalorimeterEntryPointDepth(),
430 reconstructor->IsTableMassUsed(), reco 430 reconstructor->IsTableMassUsed(), reconstructor->IsMassCutUsed(),
431 reconstructor->GetMassCutOPCenter(), 431 reconstructor->GetMassCutOPCenter(),
432 reconstructor->GetMassCutNOPCenter(), 432 reconstructor->GetMassCutNOPCenter(),
433 reconstructor->GetMassCutOPWidth(), re 433 reconstructor->GetMassCutOPWidth(), reconstructor->GetMassCutNOPWidth(),
434 reconstructor->GetMassCutEllipseAngle( 434 reconstructor->GetMassCutEllipseAngle(),
435 reconstructor->IsAbsorbedEnergyCutUsed 435 reconstructor->IsAbsorbedEnergyCutUsed(),
436 reconstructor->GetAbsorbedEnergyCutCLC 436 reconstructor->GetAbsorbedEnergyCutCLCenter(),
437 reconstructor->GetAbsorbedEnergyCutCRC 437 reconstructor->GetAbsorbedEnergyCutCRCenter(),
438 reconstructor->GetAbsorbedEnergyCutCLW 438 reconstructor->GetAbsorbedEnergyCutCLWidth(),
439 reconstructor->GetAbsorbedEnergyCutCRW 439 reconstructor->GetAbsorbedEnergyCutCRWidth(),
440 reconstructor->GetAbsorbedEnergyCutEll 440 reconstructor->GetAbsorbedEnergyCutEllipseAngle(),
441 nmbOfHitsSampled, nmbOfHitsSampledFull 441 nmbOfHitsSampled, nmbOfHitsSampledFull, nmbOfHitsTriggeredRealRange,
442 nmbOfHitsTriggeredRecRange, nmbOfOrpha 442 nmbOfHitsTriggeredRecRange, nmbOfOrphanHits, nmbOfFalseHitsTriggeredEDT,
443 nmbOfFalseHitsTriggeredRec, nmbOfSaved 443 nmbOfFalseHitsTriggeredRec, nmbOfSavedEvents, nmbOfSavedFastEvents,
444 numberOfEventsProcessed, numberOfEvent 444 numberOfEventsProcessed, numberOfEventsProcessedEffective,
445 numberOfEventToBeProcessed, rProject, 445 numberOfEventToBeProcessed, rProject, skipInteractionsWithoutEDTonWrite,
446 cfFileName, evDataVerboseLevel, physic 446 cfFileName, evDataVerboseLevel, physicsManager->GetProposedMaxIL(),
447 reconstructor->GetExpectedMomentumAmp( 447 reconstructor->GetExpectedMomentumAmp(),
448 reconstructor->GetEDCollectionAlgorith 448 reconstructor->GetEDCollectionAlgorithm(), 0 };
449 449
450 std::ofstream runDataFile( ( projectsDir 450 std::ofstream runDataFile( ( projectsDir + "/" + projectId + ".rdb" ).
451 c_str( 451 c_str() );
452 452
453 { 453 {
454 boost::archive::binary_oarchive archi 454 boost::archive::binary_oarchive archive( runDataFile );
455 archive << sObjectToWrite; 455 archive << sObjectToWrite;
456 } 456 }
457 } 457 }
458 458
459 #endif 459 #endif
460 460
461 461
462 void CexmcRunManager::DoCommonEventLoop( G4in 462 void CexmcRunManager::DoCommonEventLoop( G4int nEvent, const G4String & cmd,
463 G4in 463 G4int nSelect )
464 { 464 {
465 G4int iEvent( 0 ); 465 G4int iEvent( 0 );
466 G4int iEventEffective( 0 ); 466 G4int iEventEffective( 0 );
467 467
468 for ( iEvent = 0; iEventEffective < nEvent 468 for ( iEvent = 0; iEventEffective < nEvent; ++iEvent )
469 { 469 {
470 currentEvent = GenerateEvent( iEvent ) 470 currentEvent = GenerateEvent( iEvent );
471 eventManager->ProcessOneEvent( current 471 eventManager->ProcessOneEvent( currentEvent );
472 CexmcEventInfo * eventInfo( static_ca 472 CexmcEventInfo * eventInfo( static_cast< CexmcEventInfo * >(
473 curren 473 currentEvent->GetUserInformation() ) );
474 switch ( eventCountPolicy ) 474 switch ( eventCountPolicy )
475 { 475 {
476 case CexmcCountAllEvents : 476 case CexmcCountAllEvents :
477 ++iEventEffective; 477 ++iEventEffective;
478 break; 478 break;
479 case CexmcCountEventsWithInteraction : 479 case CexmcCountEventsWithInteraction :
480 if ( eventInfo->TpTriggerIsOk() ) 480 if ( eventInfo->TpTriggerIsOk() )
481 ++iEventEffective; 481 ++iEventEffective;
482 break; 482 break;
483 case CexmcCountEventsWithTrigger : 483 case CexmcCountEventsWithTrigger :
484 if ( eventInfo->EdTriggerIsOk() ) 484 if ( eventInfo->EdTriggerIsOk() )
485 ++iEventEffective; 485 ++iEventEffective;
486 break; 486 break;
487 default : 487 default :
488 ++iEventEffective; 488 ++iEventEffective;
489 break; 489 break;
490 } 490 }
491 AnalyzeEvent( currentEvent ); 491 AnalyzeEvent( currentEvent );
492 UpdateScoring(); 492 UpdateScoring();
493 if ( iEvent < nSelect ) 493 if ( iEvent < nSelect )
494 G4UImanager::GetUIpointer()->Apply 494 G4UImanager::GetUIpointer()->ApplyCommand( cmd );
495 StackPreviousEvent( currentEvent ); 495 StackPreviousEvent( currentEvent );
496 currentEvent = 0; 496 currentEvent = 0;
497 if ( runAborted ) 497 if ( runAborted )
498 break; 498 break;
499 } 499 }
500 500
501 numberOfEventsProcessed = iEvent; 501 numberOfEventsProcessed = iEvent;
502 numberOfEventsProcessedEffective = iEventE 502 numberOfEventsProcessedEffective = iEventEffective;
503 } 503 }
504 504
505 505
506 #ifdef CEXMC_USE_PERSISTENCY 506 #ifdef CEXMC_USE_PERSISTENCY
507 507
508 void CexmcRunManager::DoReadEventLoop( G4int 508 void CexmcRunManager::DoReadEventLoop( G4int nEvent )
509 { 509 {
510 G4int iEvent( 0 ); 510 G4int iEvent( 0 );
511 G4int iEventEffective( 0 ); 511 G4int iEventEffective( 0 );
512 G4int nEventCount( 0 ); 512 G4int nEventCount( 0 );
513 513
514 if ( ! ProjectIsRead() ) 514 if ( ! ProjectIsRead() )
515 return; 515 return;
516 516
517 if ( ! physicsManager ) 517 if ( ! physicsManager )
518 throw CexmcException( CexmcWeirdExcept 518 throw CexmcException( CexmcWeirdException );
519 519
520 CexmcProductionModel * productionModel( 520 CexmcProductionModel * productionModel(
521 physic 521 physicsManager->GetProductionModel() );
522 if ( ! productionModel ) 522 if ( ! productionModel )
523 throw CexmcException( CexmcWeirdExcept 523 throw CexmcException( CexmcWeirdException );
524 524
525 CexmcSetup * setup( static_cast< CexmcSet 525 CexmcSetup * setup( static_cast< CexmcSetup * >( userDetector ) );
526 if ( ! setup ) 526 if ( ! setup )
527 throw CexmcException( CexmcWeirdExcept 527 throw CexmcException( CexmcWeirdException );
528 528
529 CexmcEventSObject evSObject; 529 CexmcEventSObject evSObject;
530 CexmcEventFastSObject evFastSObject; 530 CexmcEventFastSObject evFastSObject;
531 531
532 /* read events data */ 532 /* read events data */
533 std::ifstream eventsDataFile( 533 std::ifstream eventsDataFile(
534 ( projectsDir + "/" + 534 ( projectsDir + "/" + rProject + ".edb" ).c_str() );
535 if ( ! eventsDataFile ) 535 if ( ! eventsDataFile )
536 throw CexmcException( CexmcReadProject 536 throw CexmcException( CexmcReadProjectIncomplete );
537 537
538 boost::archive::binary_iarchive evArchive 538 boost::archive::binary_iarchive evArchive( eventsDataFile );
539 539
540 std::ifstream eventsFastDataFile( 540 std::ifstream eventsFastDataFile(
541 ( projectsDir + "/" + 541 ( projectsDir + "/" + rProject + ".fdb" ).c_str() );
542 if ( ! eventsFastDataFile ) 542 if ( ! eventsFastDataFile )
543 throw CexmcException( CexmcReadProject 543 throw CexmcException( CexmcReadProjectIncomplete );
544 544
545 boost::archive::binary_iarchive evFastArc 545 boost::archive::binary_iarchive evFastArchive( eventsFastDataFile );
546 546
547 G4Event event; 547 G4Event event;
548 currentEvent = &event; 548 currentEvent = &event;
549 G4SDManager * sdManager( G4SDManager: 549 G4SDManager * sdManager( G4SDManager::GetSDMpointer() );
550 event.SetHCofThisEvent( sdManager->Prepare 550 event.SetHCofThisEvent( sdManager->PrepareNewEvent() );
551 G4HCofThisEvent * hcOfThisEvent( event.Ge 551 G4HCofThisEvent * hcOfThisEvent( event.GetHCofThisEvent() );
552 552
553 G4DigiManager * digiManager( G4DigiManage 553 G4DigiManager * digiManager( G4DigiManager::GetDMpointer() );
554 554
555 G4int hcId( digiManager->GetHitsCollectio 555 G4int hcId( digiManager->GetHitsCollectionID(
556 CexmcDetectorRoleName[ 556 CexmcDetectorRoleName[ CexmcMonitorDetectorRole ] +
557 "/" + CexmcDetectorTyp 557 "/" + CexmcDetectorTypeName[ CexmcEDDetector ] ) );
558 CexmcEnergyDepositCollection * monitorED( 558 CexmcEnergyDepositCollection * monitorED(
559 ne 559 new CexmcEnergyDepositCollection );
560 hcOfThisEvent->AddHitsCollection( hcId, mo 560 hcOfThisEvent->AddHitsCollection( hcId, monitorED );
561 hcId = digiManager->GetHitsCollectionID( 561 hcId = digiManager->GetHitsCollectionID(
562 CexmcDetectorRoleName[ 562 CexmcDetectorRoleName[ CexmcVetoCounterDetectorRole ] +
563 "/" + CexmcDetectorTyp 563 "/" + CexmcDetectorTypeName[ CexmcEDDetector ] );
564 CexmcEnergyDepositCollection * vetoCounte 564 CexmcEnergyDepositCollection * vetoCounterED(
565 ne 565 new CexmcEnergyDepositCollection );
566 hcOfThisEvent->AddHitsCollection( hcId, ve 566 hcOfThisEvent->AddHitsCollection( hcId, vetoCounterED );
567 hcId = digiManager->GetHitsCollectionID( 567 hcId = digiManager->GetHitsCollectionID(
568 CexmcDetectorRoleName[ 568 CexmcDetectorRoleName[ CexmcCalorimeterDetectorRole ] +
569 "/" + CexmcDetectorTyp 569 "/" + CexmcDetectorTypeName[ CexmcEDDetector ] );
570 CexmcEnergyDepositCollection * calorimete 570 CexmcEnergyDepositCollection * calorimeterED(
571 ne 571 new CexmcEnergyDepositCollection );
572 hcOfThisEvent->AddHitsCollection( hcId, ca 572 hcOfThisEvent->AddHitsCollection( hcId, calorimeterED );
573 hcId = digiManager->GetHitsCollectionID( 573 hcId = digiManager->GetHitsCollectionID(
574 CexmcDetectorRoleName[ 574 CexmcDetectorRoleName[ CexmcMonitorDetectorRole ] +
575 "/" + CexmcDetectorTyp 575 "/" + CexmcDetectorTypeName[ CexmcTPDetector ] );
576 CexmcTrackPointsCollection * monitorTP( n 576 CexmcTrackPointsCollection * monitorTP( new CexmcTrackPointsCollection );
577 hcOfThisEvent->AddHitsCollection( hcId, mo 577 hcOfThisEvent->AddHitsCollection( hcId, monitorTP );
578 hcId = digiManager->GetHitsCollectionID( 578 hcId = digiManager->GetHitsCollectionID(
579 CexmcDetectorRoleName[ 579 CexmcDetectorRoleName[ CexmcVetoCounterDetectorRole ] +
580 "/" + CexmcDetectorTyp 580 "/" + CexmcDetectorTypeName[ CexmcTPDetector ] );
581 CexmcTrackPointsCollection * vetoCounterT 581 CexmcTrackPointsCollection * vetoCounterTP(
582 ne 582 new CexmcTrackPointsCollection );
583 hcOfThisEvent->AddHitsCollection( hcId, ve 583 hcOfThisEvent->AddHitsCollection( hcId, vetoCounterTP );
584 hcId = digiManager->GetHitsCollectionID( 584 hcId = digiManager->GetHitsCollectionID(
585 CexmcDetectorRoleName[ 585 CexmcDetectorRoleName[ CexmcCalorimeterDetectorRole ] +
586 "/" + CexmcDetectorTyp 586 "/" + CexmcDetectorTypeName[ CexmcTPDetector ] );
587 CexmcTrackPointsCollection * calorimeterT 587 CexmcTrackPointsCollection * calorimeterTP(
588 ne 588 new CexmcTrackPointsCollection );
589 hcOfThisEvent->AddHitsCollection( hcId, ca 589 hcOfThisEvent->AddHitsCollection( hcId, calorimeterTP );
590 hcId = digiManager->GetHitsCollectionID( 590 hcId = digiManager->GetHitsCollectionID(
591 CexmcDetectorRoleName[ 591 CexmcDetectorRoleName[ CexmcTargetDetectorRole ] +
592 "/" + CexmcDetectorTyp 592 "/" + CexmcDetectorTypeName[ CexmcTPDetector ] );
593 CexmcTrackPointsCollection * targetTP( ne 593 CexmcTrackPointsCollection * targetTP( new CexmcTrackPointsCollection );
594 hcOfThisEvent->AddHitsCollection( hcId, ta 594 hcOfThisEvent->AddHitsCollection( hcId, targetTP );
595 595
596 #ifdef CEXMC_USE_CUSTOM_FILTER 596 #ifdef CEXMC_USE_CUSTOM_FILTER
597 if ( customFilter ) 597 if ( customFilter )
598 customFilter->SetAddressedData( &evFas 598 customFilter->SetAddressedData( &evFastSObject, &evSObject );
599 #endif 599 #endif
600 600
601 G4int nmbOfSavedEvents( rEvDataVerboseLe 601 G4int nmbOfSavedEvents( rEvDataVerboseLevel == CexmcWriteNoEventData ? 0 :
602 sObject.nmbOfSave 602 sObject.nmbOfSavedFastEvents );
603 G4bool eventDataWrittenOnEveryTPT( rEvDat 603 G4bool eventDataWrittenOnEveryTPT( rEvDataVerboseLevel ==
604 CexmcW 604 CexmcWriteEventDataOnEveryTPT );
605 605
606 for ( int i( 0 ); i < nmbOfSavedEvents; + 606 for ( int i( 0 ); i < nmbOfSavedEvents; ++i )
607 { 607 {
608 evFastArchive >> evFastSObject; 608 evFastArchive >> evFastSObject;
609 609
610 if ( nEventCount < curEventRead ) 610 if ( nEventCount < curEventRead )
611 { 611 {
612 if ( eventDataWrittenOnEveryTPT || 612 if ( eventDataWrittenOnEveryTPT ||
613 evFastSObject.edDigitizerHasT 613 evFastSObject.edDigitizerHasTriggered )
614 { 614 {
615 evArchive >> evSObject; 615 evArchive >> evSObject;
616 if ( evFastSObject.edDigitizer 616 if ( evFastSObject.edDigitizerHasTriggered )
617 ++nEventCount; 617 ++nEventCount;
618 } 618 }
619 continue; 619 continue;
620 } 620 }
621 621
622 ++iEvent; 622 ++iEvent;
623 623
624 productionModel->SetTriggeredAngularRa 624 productionModel->SetTriggeredAngularRanges(
625 625 evFastSObject.opCosThetaSCM );
626 const CexmcAngularRangeList & trigger 626 const CexmcAngularRangeList & triggeredAngularRanges(
627 productionMode 627 productionModel->GetTriggeredAngularRanges() );
628 628
629 if ( ! eventDataWrittenOnEveryTPT && 629 if ( ! eventDataWrittenOnEveryTPT &&
630 ! evFastSObject.edDigitizerHasTri 630 ! evFastSObject.edDigitizerHasTriggered )
631 { 631 {
632 #ifdef CEXMC_USE_CUSTOM_FILTER 632 #ifdef CEXMC_USE_CUSTOM_FILTER
633 /* user must be aware that using t 633 /* user must be aware that using tpt commands in custom filter
634 * scripts for poor event data set 634 * scripts for poor event data sets can easily lead to logical
635 * errors! This is because most of 635 * errors! This is because most of tpt data is only available for
636 * events with EDT trigger. There 636 * events with EDT trigger. There is no such problem if event data
637 * was written on every TPT event. 637 * was written on every TPT event. */
638 if ( customFilter && ! customFilte 638 if ( customFilter && ! customFilter->EvalTPT() )
639 continue; 639 continue;
640 #endif 640 #endif
641 SaveCurrentTPTEvent( evFastSObject 641 SaveCurrentTPTEvent( evFastSObject, triggeredAngularRanges,
642 ProjectIsSaved() && ! ski 642 ProjectIsSaved() && ! skipInteractionsWithoutEDTonWrite );
643 continue; 643 continue;
644 } 644 }
645 645
646 evArchive >> evSObject; 646 evArchive >> evSObject;
647 647
648 G4bool skipEDTOnThisEvent( false ); 648 G4bool skipEDTOnThisEvent( false );
649 649
650 #ifdef CEXMC_USE_CUSTOM_FILTER 650 #ifdef CEXMC_USE_CUSTOM_FILTER
651 if ( customFilter && ! customFilter->E 651 if ( customFilter && ! customFilter->EvalTPT() )
652 continue; 652 continue;
653 if ( customFilter && ! customFilter->E 653 if ( customFilter && ! customFilter->EvalEDT() )
654 { 654 {
655 if ( ! eventDataWrittenOnEveryTPT 655 if ( ! eventDataWrittenOnEveryTPT )
656 { 656 {
657 SaveCurrentTPTEvent( evFastSOb 657 SaveCurrentTPTEvent( evFastSObject, triggeredAngularRanges,
658 ProjectIs 658 ProjectIsSaved() );
659 continue; 659 continue;
660 } 660 }
661 skipEDTOnThisEvent = true; 661 skipEDTOnThisEvent = true;
662 } 662 }
663 #endif 663 #endif
664 664
665 event.SetEventID( evSObject.eventId ); 665 event.SetEventID( evSObject.eventId );
666 666
667 /* AAA: beware! If something throws an 667 /* AAA: beware! If something throws an exception between AAA and CCC
668 * then there would be segmentation vi 668 * then there would be segmentation violation in ~THitsMap() as it
669 * would try to delete local variable 669 * would try to delete local variable evSObject's fields.
670 * See also BBB */ 670 * See also BBB */
671 671
672 monitorED->GetMap()->operator[]( 0 ) = 672 monitorED->GetMap()->operator[]( 0 ) = &evSObject.monitorED;
673 vetoCounterED->GetMap()->operator[]( 0 673 vetoCounterED->GetMap()->operator[]( 0 ) = &evSObject.vetoCounterEDLeft;
674 vetoCounterED->GetMap()->operator[]( 1 674 vetoCounterED->GetMap()->operator[]( 1 <<
675 CexmcEnergyDepositInLeftRightS 675 CexmcEnergyDepositInLeftRightSet::GetLeftRightBitsOffset() ) =
676 676 &evSObject.vetoCounterEDRight;
677 G4int row( 0 ); 677 G4int row( 0 );
678 G4int column( 0 ); 678 G4int column( 0 );
679 for ( CexmcEnergyDepositCalorimeterCol 679 for ( CexmcEnergyDepositCalorimeterCollection::iterator
680 k( evSObject.calorimeterEDLeft 680 k( evSObject.calorimeterEDLeftCollection.begin() );
681 k != evSObject.calorimeter 681 k != evSObject.calorimeterEDLeftCollection.end(); ++k )
682 { 682 {
683 G4int index( row << 683 G4int index( row <<
684 CexmcEnergyDepositInCalorimete 684 CexmcEnergyDepositInCalorimeter::GetCopyDepth1BitsOffset() );
685 column = 0; 685 column = 0;
686 for ( CexmcEnergyDepositCrystalRow 686 for ( CexmcEnergyDepositCrystalRowCollection::iterator
687 l( k->begin() ); l != k->e 687 l( k->begin() ); l != k->end(); ++l )
688 { 688 {
689 calorimeterED->GetMap()->opera 689 calorimeterED->GetMap()->operator[]( index | column ) = &*l;
690 ++column; 690 ++column;
691 } 691 }
692 ++row; 692 ++row;
693 } 693 }
694 row = 0; 694 row = 0;
695 for ( CexmcEnergyDepositCalorimeterCol 695 for ( CexmcEnergyDepositCalorimeterCollection::iterator
696 k( evSObject.calorimeterEDRigh 696 k( evSObject.calorimeterEDRightCollection.begin() );
697 k != evSObject.calorimeter 697 k != evSObject.calorimeterEDRightCollection.end(); ++k )
698 { 698 {
699 G4int index( 699 G4int index(
700 1 << CexmcEnergyDepositInLeftR 700 1 << CexmcEnergyDepositInLeftRightSet::GetLeftRightBitsOffset()
701 | row << 701 | row <<
702 CexmcEnergyDepositInCalorimete 702 CexmcEnergyDepositInCalorimeter::GetCopyDepth1BitsOffset() );
703 column = 0; 703 column = 0;
704 for ( CexmcEnergyDepositCrystalRow 704 for ( CexmcEnergyDepositCrystalRowCollection::iterator
705 l( k->begin() ); l != k->e 705 l( k->begin() ); l != k->end(); ++l )
706 { 706 {
707 calorimeterED->GetMap()->opera 707 calorimeterED->GetMap()->operator[]( index | column ) = &*l;
708 ++column; 708 ++column;
709 } 709 }
710 ++row; 710 ++row;
711 } 711 }
712 712
713 CexmcTrackPointInfo monitorTPInfo( ev 713 CexmcTrackPointInfo monitorTPInfo( evSObject.monitorTP );
714 CexmcTrackPointInfo targetTPBeamParti 714 CexmcTrackPointInfo targetTPBeamParticleInfo(
715 evSObj 715 evSObject.targetTPBeamParticle );
716 CexmcTrackPointInfo targetTPOutputPar 716 CexmcTrackPointInfo targetTPOutputParticleInfo(
717 evSObj 717 evSObject.targetTPOutputParticle );
718 CexmcTrackPointInfo targetTPNucleusPa 718 CexmcTrackPointInfo targetTPNucleusParticleInfo(
719 evSObj 719 evSObject.targetTPNucleusParticle );
720 CexmcTrackPointInfo targetTPOutputPar 720 CexmcTrackPointInfo targetTPOutputParticleDecayProductParticle1Info(
721 evSObject.targetTPOutp 721 evSObject.targetTPOutputParticleDecayProductParticle1 );
722 CexmcTrackPointInfo targetTPOutputPar 722 CexmcTrackPointInfo targetTPOutputParticleDecayProductParticle2Info(
723 evSObject.targetTPOutp 723 evSObject.targetTPOutputParticleDecayProductParticle2 );
724 CexmcTrackPointInfo vetoCounterTPLeft 724 CexmcTrackPointInfo vetoCounterTPLeftInfo(
725 725 evSObject.vetoCounterTPLeft );
726 CexmcTrackPointInfo vetoCounterTPRigh 726 CexmcTrackPointInfo vetoCounterTPRightInfo(
727 727 evSObject.vetoCounterTPRight );
728 CexmcTrackPointInfo calorimeterTPLeft 728 CexmcTrackPointInfo calorimeterTPLeftInfo(
729 729 evSObject.calorimeterTPLeft );
730 CexmcTrackPointInfo calorimeterTPRigh 730 CexmcTrackPointInfo calorimeterTPRightInfo(
731 731 evSObject.calorimeterTPRight );
732 732
733 if ( monitorTPInfo.IsValid() ) 733 if ( monitorTPInfo.IsValid() )
734 monitorTP->GetMap()->operator[]( m 734 monitorTP->GetMap()->operator[]( monitorTPInfo.trackId ) =
735 735 &monitorTPInfo;
736 if ( targetTPBeamParticleInfo.IsValid( 736 if ( targetTPBeamParticleInfo.IsValid() )
737 targetTP->GetMap()->operator[]( 737 targetTP->GetMap()->operator[](
738 targetTPBeamParticleInfo.t 738 targetTPBeamParticleInfo.trackId ) =
739 739 &targetTPBeamParticleInfo;
740 if ( targetTPOutputParticleInfo.IsVali 740 if ( targetTPOutputParticleInfo.IsValid() )
741 targetTP->GetMap()->operator[]( 741 targetTP->GetMap()->operator[](
742 targetTPOutputParticleInfo 742 targetTPOutputParticleInfo.trackId ) =
743 743 &targetTPOutputParticleInfo;
744 if ( targetTPNucleusParticleInfo.IsVal 744 if ( targetTPNucleusParticleInfo.IsValid() )
745 targetTP->GetMap()->operator[]( 745 targetTP->GetMap()->operator[](
746 targetTPNucleusParticleInf 746 targetTPNucleusParticleInfo.trackId ) =
747 747 &targetTPNucleusParticleInfo;
748 if ( targetTPOutputParticleDecayProduc 748 if ( targetTPOutputParticleDecayProductParticle1Info.IsValid() )
749 targetTP->GetMap()->operator[]( 749 targetTP->GetMap()->operator[](
750 targetTPOutputParticleDeca 750 targetTPOutputParticleDecayProductParticle1Info.trackId ) =
751 &targetTPOutputPar 751 &targetTPOutputParticleDecayProductParticle1Info;
752 if ( targetTPOutputParticleDecayProduc 752 if ( targetTPOutputParticleDecayProductParticle2Info.IsValid() )
753 targetTP->GetMap()->operator[]( 753 targetTP->GetMap()->operator[](
754 targetTPOutputParticleDeca 754 targetTPOutputParticleDecayProductParticle2Info.trackId ) =
755 &targetTPOutputPar 755 &targetTPOutputParticleDecayProductParticle2Info;
756 if ( vetoCounterTPLeftInfo.IsValid() ) 756 if ( vetoCounterTPLeftInfo.IsValid() )
757 vetoCounterTP->GetMap()->operator[ 757 vetoCounterTP->GetMap()->operator[](
758 vetoCounterTPLeftInfo.trac 758 vetoCounterTPLeftInfo.trackId ) = &vetoCounterTPLeftInfo;
759 if ( vetoCounterTPRightInfo.IsValid() 759 if ( vetoCounterTPRightInfo.IsValid() )
760 vetoCounterTP->GetMap()->operator[ 760 vetoCounterTP->GetMap()->operator[](
761 1 << CexmcTrackPointsInLeftRig 761 1 << CexmcTrackPointsInLeftRightSet::GetLeftRightBitsOffset() |
762 vetoCounterTPRightInfo.tra 762 vetoCounterTPRightInfo.trackId ) = &vetoCounterTPRightInfo;
763 763
764 G4ThreeVector pos; 764 G4ThreeVector pos;
765 if ( calorimeterTPLeftInfo.IsValid() ) 765 if ( calorimeterTPLeftInfo.IsValid() )
766 { 766 {
767 pos = calorimeterTPLeftInfo.positi 767 pos = calorimeterTPLeftInfo.positionLocal;
768 setup->ConvertToCrystalGeometry( 768 setup->ConvertToCrystalGeometry(
769 calorimeterTPLeftInfo.posi 769 calorimeterTPLeftInfo.positionLocal, row, column, pos );
770 calorimeterTPLeftInfo.positionLoca 770 calorimeterTPLeftInfo.positionLocal = pos;
771 calorimeterTP->GetMap()->operator[ 771 calorimeterTP->GetMap()->operator[](
772 row << CexmcTrackPointsInCalor 772 row << CexmcTrackPointsInCalorimeter::
773 773 GetCopyDepth1BitsOffset() |
774 column << CexmcTrackPointsInCa 774 column << CexmcTrackPointsInCalorimeter::
775 775 GetCopyDepth0BitsOffset() |
776 calorimeterTPLeftInfo.trackId 776 calorimeterTPLeftInfo.trackId ) = &calorimeterTPLeftInfo;
777 } 777 }
778 if ( calorimeterTPRightInfo.IsValid() 778 if ( calorimeterTPRightInfo.IsValid() )
779 { 779 {
780 pos = calorimeterTPRightInfo.posit 780 pos = calorimeterTPRightInfo.positionLocal;
781 setup->ConvertToCrystalGeometry( 781 setup->ConvertToCrystalGeometry(
782 calorimeterTPRightInfo.pos 782 calorimeterTPRightInfo.positionLocal, row, column, pos );
783 calorimeterTPRightInfo.positionLoc 783 calorimeterTPRightInfo.positionLocal = pos;
784 calorimeterTP->GetMap()->operator[ 784 calorimeterTP->GetMap()->operator[](
785 1 << CexmcTrackPointsInLeftRig 785 1 << CexmcTrackPointsInLeftRightSet::GetLeftRightBitsOffset() |
786 row << CexmcTrackPointsInCalor 786 row << CexmcTrackPointsInCalorimeter::
787 787 GetCopyDepth1BitsOffset() |
788 column << CexmcTrackPointsInCa 788 column << CexmcTrackPointsInCalorimeter::
789 789 GetCopyDepth0BitsOffset() |
790 calorimeterTPRightInfo.trackId 790 calorimeterTPRightInfo.trackId ) = &calorimeterTPRightInfo;
791 } 791 }
792 792
793 productionModel->SetProductionModelDat 793 productionModel->SetProductionModelData(
794 ev 794 evSObject.productionModelData );
795 795
796 const CexmcEventAction * eventAction( 796 const CexmcEventAction * eventAction(
797 static_cast< const CexmcEventA 797 static_cast< const CexmcEventAction * >( userEventAction ) );
798 if ( ! eventAction ) 798 if ( ! eventAction )
799 { 799 {
800 /* BBB: all hits collections must 800 /* BBB: all hits collections must be cleared before throwing
801 * anything from here, otherwise ~ 801 * anything from here, otherwise ~THitsMap() will try to delete
802 * local variable evSObject's fiel 802 * local variable evSObject's fields like monitorED etc. */
803 monitorED->GetMap()->clear(); 803 monitorED->GetMap()->clear();
804 vetoCounterED->GetMap()->clear(); 804 vetoCounterED->GetMap()->clear();
805 calorimeterED->GetMap()->clear(); 805 calorimeterED->GetMap()->clear();
806 monitorTP->GetMap()->clear(); 806 monitorTP->GetMap()->clear();
807 targetTP->GetMap()->clear(); 807 targetTP->GetMap()->clear();
808 vetoCounterTP->GetMap()->clear(); 808 vetoCounterTP->GetMap()->clear();
809 calorimeterTP->GetMap()->clear(); 809 calorimeterTP->GetMap()->clear();
810 throw CexmcException( CexmcEventAc 810 throw CexmcException( CexmcEventActionIsNotInitialized );
811 } 811 }
812 812
813 if ( skipEDTOnThisEvent ) 813 if ( skipEDTOnThisEvent )
814 event.SetUserInformation( new Cexm 814 event.SetUserInformation( new CexmcEventInfo( false, false,
815 815 false ) );
816 816
817 CexmcEventAction * theEventAction( co 817 CexmcEventAction * theEventAction( const_cast< CexmcEventAction * >(
818 818 eventAction ) );
819 theEventAction->EndOfEventAction( &eve 819 theEventAction->EndOfEventAction( &event );
820 820
821 CexmcEventInfo * eventInfo( static_ca 821 CexmcEventInfo * eventInfo( static_cast< CexmcEventInfo * >(
822 822 event.GetUserInformation() ) );
823 823
824 if ( eventInfo->EdTriggerIsOk() ) 824 if ( eventInfo->EdTriggerIsOk() )
825 ++iEventEffective; 825 ++iEventEffective;
826 826
827 delete eventInfo; 827 delete eventInfo;
828 event.SetUserInformation( NULL ); 828 event.SetUserInformation( NULL );
829 829
830 monitorED->GetMap()->clear(); 830 monitorED->GetMap()->clear();
831 vetoCounterED->GetMap()->clear(); 831 vetoCounterED->GetMap()->clear();
832 calorimeterED->GetMap()->clear(); 832 calorimeterED->GetMap()->clear();
833 monitorTP->GetMap()->clear(); 833 monitorTP->GetMap()->clear();
834 targetTP->GetMap()->clear(); 834 targetTP->GetMap()->clear();
835 vetoCounterTP->GetMap()->clear(); 835 vetoCounterTP->GetMap()->clear();
836 calorimeterTP->GetMap()->clear(); 836 calorimeterTP->GetMap()->clear();
837 837
838 /* CCC: see AAA */ 838 /* CCC: see AAA */
839 839
840 if ( nEvent > 0 && iEventEffective == 840 if ( nEvent > 0 && iEventEffective == nEvent )
841 break; 841 break;
842 } 842 }
843 843
844 curEventRead = nEventCount + iEventEffecti 844 curEventRead = nEventCount + iEventEffective;
845 845
846 numberOfEventsProcessed = iEvent; 846 numberOfEventsProcessed = iEvent;
847 numberOfEventsProcessedEffective = iEventE 847 numberOfEventsProcessedEffective = iEventEffective;
848 848
849 #ifdef CEXMC_USE_CUSTOM_FILTER 849 #ifdef CEXMC_USE_CUSTOM_FILTER
850 if ( customFilter ) 850 if ( customFilter )
851 customFilter->SetAddressedData( NULL, 851 customFilter->SetAddressedData( NULL, NULL );
852 #endif 852 #endif
853 } 853 }
854 854
855 855
856 void CexmcRunManager::SaveCurrentTPTEvent( 856 void CexmcRunManager::SaveCurrentTPTEvent(
857 const CexmcEve 857 const CexmcEventFastSObject & evFastSObject,
858 const CexmcAng 858 const CexmcAngularRangeList & angularRanges,
859 G4bool writeT 859 G4bool writeToDatabase )
860 { 860 {
861 CexmcRun * run( static_cast< CexmcRun * > 861 CexmcRun * run( static_cast< CexmcRun * >( currentRun ) );
862 862
863 if ( ! run ) 863 if ( ! run )
864 return; 864 return;
865 865
866 for ( CexmcAngularRangeList::const_iterato 866 for ( CexmcAngularRangeList::const_iterator k( angularRanges.begin() );
867 867 k != angularRanges.end(); ++k )
868 { 868 {
869 run->IncrementNmbOfHitsSampledFull( k- 869 run->IncrementNmbOfHitsSampledFull( k->index );
870 if ( evFastSObject.edDigitizerMonitorH 870 if ( evFastSObject.edDigitizerMonitorHasTriggered )
871 run->IncrementNmbOfHitsSampled( k- 871 run->IncrementNmbOfHitsSampled( k->index );
872 } 872 }
873 873
874 if ( writeToDatabase ) 874 if ( writeToDatabase )
875 { 875 {
876 fastEventsArchive->operator<<( evFastS 876 fastEventsArchive->operator<<( evFastSObject );
877 run->IncrementNmbOfSavedFastEvents(); 877 run->IncrementNmbOfSavedFastEvents();
878 } 878 }
879 } 879 }
880 880
881 #endif 881 #endif
882 882
883 883
884 /* mostly adopted from G4RunManager::DoEventLo 884 /* mostly adopted from G4RunManager::DoEventLoop() */
885 void CexmcRunManager::DoEventLoop( G4int nEv 885 void CexmcRunManager::DoEventLoop( G4int nEvent, const char * macroFile,
886 G4int nSe 886 G4int nSelect )
887 { 887 {
888 if ( verboseLevel > 0 ) 888 if ( verboseLevel > 0 )
889 timer->Start(); 889 timer->Start();
890 890
891 G4String cmd; 891 G4String cmd;
892 if ( macroFile != 0 ) 892 if ( macroFile != 0 )
893 { 893 {
894 if ( nSelect < 0 ) 894 if ( nSelect < 0 )
895 nSelect = nEvent; 895 nSelect = nEvent;
896 cmd = "/control/execute "; 896 cmd = "/control/execute ";
897 cmd += macroFile; 897 cmd += macroFile;
898 } 898 }
899 else 899 else
900 { 900 {
901 nSelect = -1; 901 nSelect = -1;
902 } 902 }
903 903
904 numberOfEventsProcessed = 0; 904 numberOfEventsProcessed = 0;
905 numberOfEventsProcessedEffective = 0; 905 numberOfEventsProcessedEffective = 0;
906 906
907 #ifdef CEXMC_USE_PERSISTENCY 907 #ifdef CEXMC_USE_PERSISTENCY
908 eventsArchive = NULL; 908 eventsArchive = NULL;
909 fastEventsArchive = NULL; 909 fastEventsArchive = NULL;
910 if ( ProjectIsRead() ) 910 if ( ProjectIsRead() )
911 { 911 {
912 if ( ProjectIsSaved() ) 912 if ( ProjectIsSaved() )
913 { 913 {
914 std::ofstream eventsDataFile( 914 std::ofstream eventsDataFile(
915 ( projectsDir + "/" + 915 ( projectsDir + "/" + projectId + ".edb" ).c_str() );
916 boost::archive::binary_oarchive e 916 boost::archive::binary_oarchive eventsArchive_( eventsDataFile );
917 std::ofstream fastEventsDataFile 917 std::ofstream fastEventsDataFile(
918 ( projectsDir + "/" + 918 ( projectsDir + "/" + projectId + ".fdb" ).c_str() );
919 boost::archive::binary_oarchive f 919 boost::archive::binary_oarchive fastEventsArchive_(
920 920 fastEventsDataFile );
921 eventsArchive = &eventsArchive_; 921 eventsArchive = &eventsArchive_;
922 fastEventsArchive = &fastEventsArc 922 fastEventsArchive = &fastEventsArchive_;
923 DoReadEventLoop( nEvent ); 923 DoReadEventLoop( nEvent );
924 } 924 }
925 else 925 else
926 { 926 {
927 DoReadEventLoop( nEvent ); 927 DoReadEventLoop( nEvent );
928 } 928 }
929 } 929 }
930 else 930 else
931 { 931 {
932 if ( ProjectIsSaved() ) 932 if ( ProjectIsSaved() )
933 { 933 {
934 std::ofstream eventsDataFile( 934 std::ofstream eventsDataFile(
935 ( projectsDir + "/" + 935 ( projectsDir + "/" + projectId + ".edb" ).c_str() );
936 boost::archive::binary_oarchive e 936 boost::archive::binary_oarchive eventsArchive_( eventsDataFile );
937 std::ofstream fastEventsDataFile 937 std::ofstream fastEventsDataFile(
938 ( projectsDir + "/" + 938 ( projectsDir + "/" + projectId + ".fdb" ).c_str() );
939 boost::archive::binary_oarchive f 939 boost::archive::binary_oarchive fastEventsArchive_(
940 940 fastEventsDataFile );
941 eventsArchive = &eventsArchive_; 941 eventsArchive = &eventsArchive_;
942 fastEventsArchive = &fastEventsArc 942 fastEventsArchive = &fastEventsArchive_;
943 DoCommonEventLoop( nEvent, cmd, nS 943 DoCommonEventLoop( nEvent, cmd, nSelect );
944 } 944 }
945 else 945 else
946 { 946 {
947 DoCommonEventLoop( nEvent, cmd, nS 947 DoCommonEventLoop( nEvent, cmd, nSelect );
948 } 948 }
949 } 949 }
950 eventsArchive = NULL; 950 eventsArchive = NULL;
951 fastEventsArchive = NULL; 951 fastEventsArchive = NULL;
952 #else 952 #else
953 DoCommonEventLoop( nEvent, cmd, nSelect ); 953 DoCommonEventLoop( nEvent, cmd, nSelect );
954 #endif 954 #endif
955 955
956 if ( verboseLevel > 0 ) 956 if ( verboseLevel > 0 )
957 { 957 {
958 timer->Stop(); 958 timer->Stop();
959 G4cout << "Run terminated." << G4endl; 959 G4cout << "Run terminated." << G4endl;
960 G4cout << "Run Summary" << G4endl; 960 G4cout << "Run Summary" << G4endl;
961 if ( runAborted ) 961 if ( runAborted )
962 { 962 {
963 G4cout << " Run Aborted after " < 963 G4cout << " Run Aborted after " << numberOfEventsProcessed <<
964 " events processed." << 964 " events processed." << G4endl;
965 } 965 }
966 else 966 else
967 { 967 {
968 G4cout << " Number of events proc 968 G4cout << " Number of events processed : " <<
969 numberOfEventsProcessed 969 numberOfEventsProcessed << ", effectively: " <<
970 numberOfEventsProcessedE 970 numberOfEventsProcessedEffective << G4endl;
971 } 971 }
972 G4cout << " " << *timer << G4endl; 972 G4cout << " " << *timer << G4endl;
973 } 973 }
974 } 974 }
975 975
976 976
977 #ifdef CEXMC_USE_PERSISTENCY 977 #ifdef CEXMC_USE_PERSISTENCY
978 978
979 void CexmcRunManager::PrintReadRunData( void 979 void CexmcRunManager::PrintReadRunData( void ) const
980 { 980 {
981 if ( ! ProjectIsRead() ) 981 if ( ! ProjectIsRead() )
982 return; 982 return;
983 983
984 G4bool refCrystalInfoPrinted( false ); 984 G4bool refCrystalInfoPrinted( false );
985 985
986 G4cout << CEXMC_LINE_START << "Run data re 986 G4cout << CEXMC_LINE_START << "Run data read from project '" << rProject <<
987 "'" << G4endl; 987 "'" << G4endl;
988 G4cout << " (archive class v 988 G4cout << " (archive class version " <<
989 sObject.actualVersion << ")" << 989 sObject.actualVersion << ")" << G4endl;
990 if ( ! sObject.rProject.empty() ) 990 if ( ! sObject.rProject.empty() )
991 { 991 {
992 G4cout << " -- Based on project '" << 992 G4cout << " -- Based on project '" << sObject.rProject << "'" <<
993 G4endl; 993 G4endl;
994 if ( ! sObject.cfFileName.empty() ) 994 if ( ! sObject.cfFileName.empty() )
995 G4cout << " -- Custom filter scri 995 G4cout << " -- Custom filter script '" << sObject.cfFileName <<
996 "' was used" << G4endl; 996 "' was used" << G4endl;
997 } 997 }
998 G4cout << " -- Event data verbose level ( 998 G4cout << " -- Event data verbose level (0 - not saved, 1 - triggers, "
999 "2 - interactions): " << sObject 999 "2 - interactions): " << sObject.evDataVerboseLevel << G4endl;
1000 if ( ! sObject.rProject.empty() ) 1000 if ( ! sObject.rProject.empty() )
1001 { 1001 {
1002 if ( sObject.evDataVerboseLevel == Ce 1002 if ( sObject.evDataVerboseLevel == CexmcWriteEventDataOnEveryEDT )
1003 { 1003 {
1004 G4cout << " -- (fdb file contain 1004 G4cout << " -- (fdb file contains " <<
1005 ( sObject.interactionsW 1005 ( sObject.interactionsWithoutEDTWereSkipped ?
1006 "only interactions wh 1006 "only interactions when an event was triggered" :
1007 "all interactions" ) 1007 "all interactions" ) << ")" << std::endl;
1008 } 1008 }
1009 } 1009 }
1010 G4cout << " -- Base physics used" << 1010 G4cout << " -- Base physics used (1 - QGSP_BERT, 2 - QGSP_BIC_EMY): " <<
1011 "(1 - QGSP_BERT, 2 - QGSP_BIC_E <<
1012 sObject.basePhysicsUsed << G4en 1011 sObject.basePhysicsUsed << G4endl;
1013 G4cout << " -- Production model (1 - pi0 1012 G4cout << " -- Production model (1 - pi0, 2 - eta): " <<
1014 sObject.productionModelType << 1013 sObject.productionModelType << G4endl;
1015 G4cout << " -- Geometry definition file: 1014 G4cout << " -- Geometry definition file: " << sObject.gdmlFileName <<
1016 G4endl; 1015 G4endl;
1017 G4cout << " -- Angular ranges: " << sObj 1016 G4cout << " -- Angular ranges: " << sObject.angularRanges << G4endl;
1018 G4cout << " -- Eta decay modes: " << G4e 1017 G4cout << " -- Eta decay modes: " << G4endl;
1019 G4Eta::Definition()->GetDecayTable()->Dum 1018 G4Eta::Definition()->GetDecayTable()->DumpInfo();
1020 G4cout << " -- Fermi motion status (0 - 1019 G4cout << " -- Fermi motion status (0 - disabled, 1 - enabled): " <<
1021 sObject.fermiMotionIsOn << G4en 1020 sObject.fermiMotionIsOn << G4endl;
1022 if ( sObject.calorimeterRegCuts.size() < 1021 if ( sObject.calorimeterRegCuts.size() < 4 )
1023 throw CexmcException( CexmcWeirdExcep 1022 throw CexmcException( CexmcWeirdException );
1024 G4cout << " -- Production cuts in calori 1023 G4cout << " -- Production cuts in calorimeter (gamma, e-, e+, p): " <<
1025 G4BestUnit( sObject.calorimeter 1024 G4BestUnit( sObject.calorimeterRegCuts[ 0 ], "Length" ) << ", " <<
1026 G4BestUnit( sObject.calorimeter 1025 G4BestUnit( sObject.calorimeterRegCuts[ 1 ], "Length" ) << ", " <<
1027 G4BestUnit( sObject.calorimeter 1026 G4BestUnit( sObject.calorimeterRegCuts[ 2 ], "Length" ) << ", " <<
1028 G4BestUnit( sObject.calorimeter 1027 G4BestUnit( sObject.calorimeterRegCuts[ 3 ], "Length" ) << G4endl;
1029 G4cout << " -- Proposed max interaction 1028 G4cout << " -- Proposed max interaction length in the target: " <<
1030 G4BestUnit( sObject.proposedMax 1029 G4BestUnit( sObject.proposedMaxIL, "Length" ) << G4endl;
1031 G4cout << " -- Event count policy (0 - a 1030 G4cout << " -- Event count policy (0 - all, 1 - interaction, 2 - trigger)"
1032 ": " << sObject.eventCountPolic 1031 ": " << sObject.eventCountPolicy << G4endl;
1033 G4cout << " -- Number of events (process 1032 G4cout << " -- Number of events (processed / effective / ordered): " <<
1034 sObject.numberOfEventsProcessed 1033 sObject.numberOfEventsProcessed << " / " <<
1035 sObject.numberOfEventsProcessed 1034 sObject.numberOfEventsProcessedEffective << " / " <<
1036 sObject.numberOfEventsToBeProce 1035 sObject.numberOfEventsToBeProcessed << G4endl;
1037 G4cout << " -- Incident beam particle: " 1036 G4cout << " -- Incident beam particle: " << sObject.beamParticle << G4endl;
1038 G4cout << " position: " 1037 G4cout << " position: " <<
1039 G4BestUnit( sObject.beamPos, "L 1038 G4BestUnit( sObject.beamPos, "Length" ) << G4endl;
1040 G4cout << " direction: " 1039 G4cout << " direction: " <<
1041 G4ThreeVector( sObject.beamDir 1040 G4ThreeVector( sObject.beamDir ) << G4endl;
1042 G4cout << " momentum: " 1041 G4cout << " momentum: " <<
1043 G4BestUnit( sObject.beamMomentu 1042 G4BestUnit( sObject.beamMomentumAmp, "Energy" ) << G4endl;
1044 G4cout << " momentum fwhm: " 1043 G4cout << " momentum fwhm: " << sObject.beamFwhmMomentumAmp <<
1045 G4endl; 1044 G4endl;
1046 G4cout << " pos fwhm (x): " 1045 G4cout << " pos fwhm (x): " <<
1047 G4BestUnit( sObject.beamFwhmPos 1046 G4BestUnit( sObject.beamFwhmPosX, "Length" ) << G4endl;
1048 G4cout << " pos fwhm (y): " 1047 G4cout << " pos fwhm (y): " <<
1049 G4BestUnit( sObject.beamFwhmPos 1048 G4BestUnit( sObject.beamFwhmPosY, "Length" ) << G4endl;
1050 G4cout << " dir fwhm (x): " 1049 G4cout << " dir fwhm (x): " << sObject.beamFwhmDirX / deg <<
1051 " deg" << G4endl; 1050 " deg" << G4endl;
1052 G4cout << " dir fwhm (y): " 1051 G4cout << " dir fwhm (y): " << sObject.beamFwhmDirY / deg <<
1053 " deg" << G4endl; 1052 " deg" << G4endl;
1054 G4cout << " -- Monitor ED threshold: " < 1053 G4cout << " -- Monitor ED threshold: " <<
1055 G4BestUnit( sObject.monitorEDTh 1054 G4BestUnit( sObject.monitorEDThreshold, "Energy" ) << G4endl;
1056 G4cout << " -- Veto counter (l/r) ED thr 1055 G4cout << " -- Veto counter (l/r) ED threshold: " <<
1057 G4BestUnit( sObject.vetoCounter 1056 G4BestUnit( sObject.vetoCounterEDLeftThreshold, "Energy" ) <<
1058 " / " << 1057 " / " <<
1059 G4BestUnit( sObject.vetoCounter 1058 G4BestUnit( sObject.vetoCounterEDRightThreshold, "Energy" ) <<
1060 G4endl; 1059 G4endl;
1061 G4cout << " -- Calorimeter (l/r) ED thre 1060 G4cout << " -- Calorimeter (l/r) ED threshold: " <<
1062 G4BestUnit( sObject.calorimeter 1061 G4BestUnit( sObject.calorimeterEDLeftThreshold, "Energy" ) <<
1063 " / " << 1062 " / " <<
1064 G4BestUnit( sObject.calorimeter 1063 G4BestUnit( sObject.calorimeterEDRightThreshold, "Energy" ) <<
1065 G4endl; 1064 G4endl;
1066 G4cout << " -- Calorimeter trigger algor 1065 G4cout << " -- Calorimeter trigger algorithm (0 - all, 1 - inner): " <<
1067 sObject.calorimeterTriggerAlgor 1066 sObject.calorimeterTriggerAlgorithm << G4endl;
1068 G4cout << " -- Outer crystals veto algor 1067 G4cout << " -- Outer crystals veto algorithm "
1069 "(0 - none, 1 - max, 2 - fracti 1068 "(0 - none, 1 - max, 2 - fraction): " <<
1070 sObject.outerCrystalsVetoAlgori 1069 sObject.outerCrystalsVetoAlgorithm << G4endl;
1071 if ( sObject.outerCrystalsVetoAlgorithm = 1070 if ( sObject.outerCrystalsVetoAlgorithm ==
1072 CexmcFractionOfEDInOuterCrystalsVeto 1071 CexmcFractionOfEDInOuterCrystalsVeto )
1073 { 1072 {
1074 G4cout << " -- Outer crystals veto f 1073 G4cout << " -- Outer crystals veto fraction: " <<
1075 sObject.outerCrystalsVetoFr 1074 sObject.outerCrystalsVetoFraction << G4endl;
1076 } 1075 }
1077 G4cout << " -- Finite crystal resolution 1076 G4cout << " -- Finite crystal resolution applied (0 - no, 1 - yes): " <<
1078 sObject.applyFiniteCrystalResol 1077 sObject.applyFiniteCrystalResolution << G4endl;
1079 if ( sObject.applyFiniteCrystalResolution 1078 if ( sObject.applyFiniteCrystalResolution )
1080 { 1079 {
1081 G4cout << " -- Crystal resolution da 1080 G4cout << " -- Crystal resolution data: " <<
1082 sObject.crystalResolutionDa 1081 sObject.crystalResolutionData;
1083 } 1082 }
1084 G4cout << " -- Reconstructor settings: " 1083 G4cout << " -- Reconstructor settings: " << G4endl;
1085 if ( sObject.expectedMomentumAmp > 0 ) 1084 if ( sObject.expectedMomentumAmp > 0 )
1086 { 1085 {
1087 G4cout << " -- expected momentum 1086 G4cout << " -- expected momentum in the target: " <<
1088 G4BestUnit( sObject.expecte 1087 G4BestUnit( sObject.expectedMomentumAmp, "Energy" ) << G4endl;
1089 } 1088 }
1090 G4cout << " -- ed collection algorith 1089 G4cout << " -- ed collection algorithm (0 - all, 1 - adjacent): " <<
1091 sObject.edCollectionAlgorithm < 1090 sObject.edCollectionAlgorithm << G4endl;
1092 if ( sObject.edCollectionAlgorithm == Cex 1091 if ( sObject.edCollectionAlgorithm == CexmcCollectEDInAdjacentCrystals )
1093 { 1092 {
1094 G4cout << 1093 G4cout <<
1095 " -- inner crystal used as re 1094 " -- inner crystal used as reference (0 - no, 1 - yes): " <<
1096 sObject.useInnerRefCrystal << G4e 1095 sObject.useInnerRefCrystal << G4endl;
1097 refCrystalInfoPrinted = true; 1096 refCrystalInfoPrinted = true;
1098 } 1097 }
1099 G4cout << " -- entry point definition 1098 G4cout << " -- entry point definition algorithm " << G4endl;
1100 G4cout << " (0 - center of calorim 1099 G4cout << " (0 - center of calorimeter, 1 - center of crystal with "
1101 "max ED," << G4endl; 1100 "max ED," << G4endl;
1102 G4cout << " 2 - linear, 3 - squar 1101 G4cout << " 2 - linear, 3 - square): " <<
1103 sObject.epDefinitionAlgorithm < 1102 sObject.epDefinitionAlgorithm << G4endl;
1104 G4cout << " -- entry point depth defi 1103 G4cout << " -- entry point depth definition algorithm "
1105 "(0 - plain, 1 - sphere 1104 "(0 - plain, 1 - sphere): " <<
1106 sObject.epDepthDefi 1105 sObject.epDepthDefinitionAlgorithm << G4endl;
1107 G4cout << " -- entry point depth: " < 1106 G4cout << " -- entry point depth: " <<
1108 G4BestUnit( sObject.epDepth, "L 1107 G4BestUnit( sObject.epDepth, "Length" ) << G4endl;
1109 if ( sObject.epDefinitionAlgorithm == Cex 1108 if ( sObject.epDefinitionAlgorithm == CexmcEntryPointByLinearEDWeights ||
1110 sObject.epDefinitionAlgorithm == Cex 1109 sObject.epDefinitionAlgorithm == CexmcEntryPointBySqrtEDWeights )
1111 { 1110 {
1112 G4cout << 1111 G4cout <<
1113 " -- crystal selection algori 1112 " -- crystal selection algorithm (0 - all, 1 - adjacent): " <<
1114 sObject.csAlgorithm << G4endl; 1113 sObject.csAlgorithm << G4endl;
1115 } 1114 }
1116 if ( ! refCrystalInfoPrinted && 1115 if ( ! refCrystalInfoPrinted &&
1117 ( sObject.epDefinitionAlgorithm == 1116 ( sObject.epDefinitionAlgorithm ==
1118 CexmcEntryPoi 1117 CexmcEntryPointInTheCenterOfCrystalWithMaxED ||
1119 ( ( sObject.epDefinitionAlgorithm == 1118 ( ( sObject.epDefinitionAlgorithm == CexmcEntryPointBySqrtEDWeights ||
1120 sObject.epDefinitionAlgorithm == 1119 sObject.epDefinitionAlgorithm ==
1121 Cexmc 1120 CexmcEntryPointByLinearEDWeights ) &&
1122 sObject.csAlgorithm == CexmcSe 1121 sObject.csAlgorithm == CexmcSelectAdjacentCrystals ) ) )
1123 { 1122 {
1124 G4cout << 1123 G4cout <<
1125 " -- inner crystal used as re 1124 " -- inner crystal used as reference (0 - no, 1 - yes): " <<
1126 sObject.useInnerRefCrystal << G4e 1125 sObject.useInnerRefCrystal << G4endl;
1127 } 1126 }
1128 G4cout << " -- table mass of output p 1127 G4cout << " -- table mass of output particle used "
1129 "(0 - no, 1 - yes): " < 1128 "(0 - no, 1 - yes): " << sObject.useTableMass << G4endl;
1130 G4cout << " -- mass cut is enabled (0 1129 G4cout << " -- mass cut is enabled (0 - no, 1 - yes): " <<
1131 sObject.useMassCut << G4endl; 1130 sObject.useMassCut << G4endl;
1132 if ( sObject.useMassCut ) 1131 if ( sObject.useMassCut )
1133 { 1132 {
1134 G4cout << " -- mass cut output pa 1133 G4cout << " -- mass cut output particle center: " <<
1135 G4BestUnit( sObject.mCutOPC 1134 G4BestUnit( sObject.mCutOPCenter, "Energy" ) << G4endl;
1136 G4cout << " -- mass cut nucleus o 1135 G4cout << " -- mass cut nucleus output particle center: " <<
1137 G4BestUnit( sObject.mCutNOP 1136 G4BestUnit( sObject.mCutNOPCenter, "Energy" ) << G4endl;
1138 G4cout << " -- mass cut output pa 1137 G4cout << " -- mass cut output particle width of the ellipse: " <<
1139 G4BestUnit( sObject.mCutOPW 1138 G4BestUnit( sObject.mCutOPWidth, "Energy" ) << G4endl;
1140 G4cout << " -- mass cut nucleus o 1139 G4cout << " -- mass cut nucleus output particle width of the "
1141 "ellipse: " 1140 "ellipse: "
1142 << G4BestUnit( sObject.mCutNOP 1141 << G4BestUnit( sObject.mCutNOPWidth, "Energy" ) << G4endl;
1143 G4cout << " -- mass cut angle of 1142 G4cout << " -- mass cut angle of the ellipse: " <<
1144 sObject.mCutAngle / deg << 1143 sObject.mCutAngle / deg << " deg" << G4endl;
1145 } 1144 }
1146 G4cout << " -- absorbed energy cut is 1145 G4cout << " -- absorbed energy cut is enabled (0 - no, 1 - yes): " <<
1147 sObject.useAbsorbedEnergyCut << 1146 sObject.useAbsorbedEnergyCut << G4endl;
1148 if ( sObject.useAbsorbedEnergyCut ) 1147 if ( sObject.useAbsorbedEnergyCut )
1149 { 1148 {
1150 G4cout << " -- absorbed energy cu 1149 G4cout << " -- absorbed energy cut left calorimeter center: " <<
1151 G4BestUnit( sObject.aeCutCL 1150 G4BestUnit( sObject.aeCutCLCenter, "Energy" ) << G4endl;
1152 G4cout << " -- absorbed energy cu 1151 G4cout << " -- absorbed energy cut right calorimeter center: " <<
1153 G4BestUnit( sObject.aeCutCR 1152 G4BestUnit( sObject.aeCutCRCenter, "Energy" ) << G4endl;
1154 G4cout << " -- absorbed energy cu 1153 G4cout << " -- absorbed energy cut left calorimeter width of the "
1155 "ellipse: " << 1154 "ellipse: " <<
1156 G4BestUnit( sObject.aeCutCL 1155 G4BestUnit( sObject.aeCutCLWidth, "Energy" ) << G4endl;
1157 G4cout << " -- absorbed energy cu 1156 G4cout << " -- absorbed energy cut right calorimeter width of the "
1158 "ellipse: " 1157 "ellipse: "
1159 << G4BestUnit( sObject.aeCutCR 1158 << G4BestUnit( sObject.aeCutCRWidth, "Energy" ) << G4endl;
1160 G4cout << " -- absorbed energy cu 1159 G4cout << " -- absorbed energy cut angle of the ellipse: " <<
1161 sObject.aeCutAngle / deg << 1160 sObject.aeCutAngle / deg << " deg" << G4endl;
1162 } 1161 }
1163 G4cout << G4endl; 1162 G4cout << G4endl;
1164 CexmcRunAction::PrintResults( sObject.nmb 1163 CexmcRunAction::PrintResults( sObject.nmbOfHitsSampled,
1165 sObject.nmb 1164 sObject.nmbOfHitsSampledFull,
1166 sObject.nmb 1165 sObject.nmbOfHitsTriggeredRealRange,
1167 sObject.nmb 1166 sObject.nmbOfHitsTriggeredRecRange,
1168 sObject.nmb 1167 sObject.nmbOfOrphanHits,
1169 sObject.ang 1168 sObject.angularRanges,
1170 sObject.nmb 1169 sObject.nmbOfFalseHitsTriggeredEDT,
1171 sObject.nmb 1170 sObject.nmbOfFalseHitsTriggeredRec );
1172 G4cout << G4endl; 1171 G4cout << G4endl;
1173 } 1172 }
1174 1173
1175 1174
1176 void CexmcRunManager::ReadAndPrintEventsData 1175 void CexmcRunManager::ReadAndPrintEventsData( void ) const
1177 { 1176 {
1178 if ( ! ProjectIsRead() ) 1177 if ( ! ProjectIsRead() )
1179 return; 1178 return;
1180 1179
1181 CexmcEventSObject evSObject; 1180 CexmcEventSObject evSObject;
1182 1181
1183 /* read events data */ 1182 /* read events data */
1184 std::ifstream eventsDataFile( 1183 std::ifstream eventsDataFile(
1185 ( projectsDir + "/" + 1184 ( projectsDir + "/" + rProject + ".edb" ).c_str() );
1186 if ( ! eventsDataFile ) 1185 if ( ! eventsDataFile )
1187 throw CexmcException( CexmcReadProjec 1186 throw CexmcException( CexmcReadProjectIncomplete );
1188 1187
1189 boost::archive::binary_iarchive evArchiv 1188 boost::archive::binary_iarchive evArchive( eventsDataFile );
1190 1189
1191 for ( int i( 0 ); i < sObject.nmbOfSaved 1190 for ( int i( 0 ); i < sObject.nmbOfSavedEvents; ++i )
1192 { 1191 {
1193 evArchive >> evSObject; 1192 evArchive >> evSObject;
1194 1193
1195 if ( ! evSObject.edDigitizerMonitorHa 1194 if ( ! evSObject.edDigitizerMonitorHasTriggered )
1196 continue; 1195 continue;
1197 1196
1198 CexmcEnergyDepositStore edStore( evS 1197 CexmcEnergyDepositStore edStore( evSObject.monitorED,
1199 evSObject.vetoCounterEDLeft, evSO 1198 evSObject.vetoCounterEDLeft, evSObject.vetoCounterEDRight,
1200 evSObject.calorimeterEDLeft, evSO 1199 evSObject.calorimeterEDLeft, evSObject.calorimeterEDRight,
1201 0, 0, 0, 0, evSObject.calorimeter 1200 0, 0, 0, 0, evSObject.calorimeterEDLeftCollection,
1202 evSObject.calorimeterEDRightColle 1201 evSObject.calorimeterEDRightCollection );
1203 1202
1204 CexmcTrackPointsStore tpStore( evS 1203 CexmcTrackPointsStore tpStore( evSObject.monitorTP,
1205 evSObject.targetTPBeamParticle, e 1204 evSObject.targetTPBeamParticle, evSObject.targetTPOutputParticle,
1206 evSObject.targetTPNucleusParticle 1205 evSObject.targetTPNucleusParticle,
1207 evSObject.targetTPOutputParticleD 1206 evSObject.targetTPOutputParticleDecayProductParticle1,
1208 evSObject.targetTPOutputParticleD 1207 evSObject.targetTPOutputParticleDecayProductParticle2,
1209 evSObject.vetoCounterTPLeft, evSO 1208 evSObject.vetoCounterTPLeft, evSObject.vetoCounterTPRight,
1210 evSObject.calorimeterTPLeft, evSO 1209 evSObject.calorimeterTPLeft, evSObject.calorimeterTPRight );
1211 1210
1212 const CexmcProductionModelData & pmD 1211 const CexmcProductionModelData & pmData(
1213 e 1212 evSObject.productionModelData );
1214 1213
1215 G4cout << "Event " << evSObject.event 1214 G4cout << "Event " << evSObject.eventId << G4endl;
1216 CexmcEventAction::PrintTrackPoints( & 1215 CexmcEventAction::PrintTrackPoints( &tpStore );
1217 G4cout << " --- Production model data 1216 G4cout << " --- Production model data: " << pmData;
1218 CexmcEventAction::PrintEnergyDeposit( 1217 CexmcEventAction::PrintEnergyDeposit( &edStore );
1219 } 1218 }
1220 } 1219 }
1221 1220
1222 1221
1223 void CexmcRunManager::PrintReadData( 1222 void CexmcRunManager::PrintReadData(
1224 const CexmcOutput 1223 const CexmcOutputDataTypeSet & outputData ) const
1225 { 1224 {
1226 if ( ! ProjectIsRead() ) 1225 if ( ! ProjectIsRead() )
1227 return; 1226 return;
1228 1227
1229 G4bool addSpace( false ); 1228 G4bool addSpace( false );
1230 1229
1231 CexmcOutputDataTypeSet::const_iterator f 1230 CexmcOutputDataTypeSet::const_iterator found(
1232 outputDat 1231 outputData.find( CexmcOutputGeometry ) );
1233 if ( found != outputData.end() ) 1232 if ( found != outputData.end() )
1234 { 1233 {
1235 G4String cmd( G4String( "cat " ) + p 1234 G4String cmd( G4String( "cat " ) + projectsDir + "/" + rProject +
1236 gdmlFileExtension ); 1235 gdmlFileExtension );
1237 if ( system( cmd ) != 0 ) 1236 if ( system( cmd ) != 0 )
1238 throw CexmcException( CexmcReadPr 1237 throw CexmcException( CexmcReadProjectIncomplete );
1239 1238
1240 if ( zipGdmlFile ) 1239 if ( zipGdmlFile )
1241 { 1240 {
1242 cmd = G4String( "bzip2 " ) + proj 1241 cmd = G4String( "bzip2 " ) + projectsDir + "/" + rProject +
1243 1242 gdmlFileExtension;
1244 if ( system( cmd ) != 0 ) 1243 if ( system( cmd ) != 0 )
1245 throw CexmcException( CexmcFi 1244 throw CexmcException( CexmcFileCompressException );
1246 } 1245 }
1247 1246
1248 addSpace = true; 1247 addSpace = true;
1249 } 1248 }
1250 1249
1251 found = outputData.find( CexmcOutputEvent 1250 found = outputData.find( CexmcOutputEvents );
1252 if ( found != outputData.end() ) 1251 if ( found != outputData.end() )
1253 { 1252 {
1254 if ( addSpace ) 1253 if ( addSpace )
1255 G4cout << G4endl << G4endl; 1254 G4cout << G4endl << G4endl;
1256 1255
1257 ReadAndPrintEventsData(); 1256 ReadAndPrintEventsData();
1258 1257
1259 addSpace = true; 1258 addSpace = true;
1260 } 1259 }
1261 1260
1262 found = outputData.find( CexmcOutputRun ) 1261 found = outputData.find( CexmcOutputRun );
1263 if ( found != outputData.end() ) 1262 if ( found != outputData.end() )
1264 { 1263 {
1265 if ( addSpace ) 1264 if ( addSpace )
1266 G4cout << G4endl << G4endl; 1265 G4cout << G4endl << G4endl;
1267 1266
1268 G4DecayTable * etaDecayTable( G4Eta: 1267 G4DecayTable * etaDecayTable( G4Eta::Definition()->GetDecayTable() );
1269 for ( CexmcDecayBranchesStore::const_ 1268 for ( CexmcDecayBranchesStore::const_iterator
1270 k( sObject.etaDecayTable.GetD 1269 k( sObject.etaDecayTable.GetDecayBranches().begin() );
1271 k != sObject.etaDecayTable.Ge 1270 k != sObject.etaDecayTable.GetDecayBranches().end(); ++k )
1272 { 1271 {
1273 etaDecayTable->GetDecayChannel( k 1272 etaDecayTable->GetDecayChannel( k->first )->SetBR( k->second );
1274 } 1273 }
1275 1274
1276 PrintReadRunData(); 1275 PrintReadRunData();
1277 } 1276 }
1278 } 1277 }
1279 1278
1280 1279
1281 #ifdef CEXMC_USE_CUSTOM_FILTER 1280 #ifdef CEXMC_USE_CUSTOM_FILTER
1282 1281
1283 void CexmcRunManager::SetCustomFilter( const 1282 void CexmcRunManager::SetCustomFilter( const G4String & cfFileName_ )
1284 { 1283 {
1285 if ( customFilter ) 1284 if ( customFilter )
1286 { 1285 {
1287 delete customFilter; 1286 delete customFilter;
1288 customFilter = NULL; 1287 customFilter = NULL;
1289 } 1288 }
1290 1289
1291 if ( cfFileName_.empty() ) 1290 if ( cfFileName_.empty() )
1292 return; 1291 return;
1293 1292
1294 /* should not get here */ 1293 /* should not get here */
1295 if ( ! ProjectIsRead() ) 1294 if ( ! ProjectIsRead() )
1296 throw CexmcException( CexmcCmdIsNotAl 1295 throw CexmcException( CexmcCmdIsNotAllowed );
1297 1296
1298 cfFileName = cfFileName_; 1297 cfFileName = cfFileName_;
1299 1298
1300 customFilter = new CexmcCustomFilterEval( 1299 customFilter = new CexmcCustomFilterEval( cfFileName );
1301 } 1300 }
1302 1301
1303 #endif 1302 #endif
1304 1303
1305 #endif 1304 #endif
1306 1305
1307 1306
1308 void CexmcRunManager::RegisterScenePrimitive 1307 void CexmcRunManager::RegisterScenePrimitives( void )
1309 { 1308 {
1310 G4VisManager * visManager( static_cast< 1309 G4VisManager * visManager( static_cast< G4VisManager * >(
1311 G4VVisMan 1310 G4VVisManager::GetConcreteInstance() ) );
1312 if ( ! visManager ) 1311 if ( ! visManager )
1313 return; 1312 return;
1314 1313
1315 G4Scene * curScene( visManager->Get 1314 G4Scene * curScene( visManager->GetCurrentScene() );
1316 if ( ! curScene ) 1315 if ( ! curScene )
1317 return; 1316 return;
1318 1317
1319 /* G4Scene declarations lack this kind of 1318 /* G4Scene declarations lack this kind of typedef */
1320 #if G4VERSION_NUMBER < 960 1319 #if G4VERSION_NUMBER < 960
1321 typedef std::vector< G4VModel * > ML 1320 typedef std::vector< G4VModel * > MList;
1322 #else 1321 #else
1323 typedef std::vector< G4Scene::Model > ML 1322 typedef std::vector< G4Scene::Model > MList;
1324 #endif 1323 #endif
1325 const MList & mList( curScene->GetRunDur 1324 const MList & mList( curScene->GetRunDurationModelList() );
1326 1325
1327 for ( MList::const_iterator k( mList.beg 1326 for ( MList::const_iterator k( mList.begin() ); k != mList.end(); ++k )
1328 { 1327 {
1329 #if G4VERSION_NUMBER < 960 1328 #if G4VERSION_NUMBER < 960
1330 const G4String & modelDesc( ( *k )-> 1329 const G4String & modelDesc( ( *k )->GetGlobalDescription() );
1331 #else 1330 #else
1332 const G4String & modelDesc( k->fpMod 1331 const G4String & modelDesc( k->fpModel->GetGlobalDescription() );
1333 #endif 1332 #endif
1334 if ( modelDesc == CexmcScenePrimitive 1333 if ( modelDesc == CexmcScenePrimitivesDescription )
1335 return; 1334 return;
1336 } 1335 }
1337 1336
1338 CexmcSetup * setup( static_cast< CexmcSe 1337 CexmcSetup * setup( static_cast< CexmcSetup * >( userDetector ) );
1339 if ( ! setup ) 1338 if ( ! setup )
1340 throw CexmcException( CexmcWeirdExcep 1339 throw CexmcException( CexmcWeirdException );
1341 1340
1342 /* BEWARE: looks like G4Scene won't delet 1341 /* BEWARE: looks like G4Scene won't delete models from its lists upon
1343 * termination! Hence destructor of the n 1342 * termination! Hence destructor of the new model won't be called */
1344 curScene->AddRunDurationModel( new CexmcS 1343 curScene->AddRunDurationModel( new CexmcScenePrimitives( setup ) );
1345 } 1344 }
1346 1345
1347 1346
1348 void CexmcRunManager::BeamParticleChangeHook 1347 void CexmcRunManager::BeamParticleChangeHook( void )
1349 { 1348 {
1350 const CexmcEventAction * eventAction( 1349 const CexmcEventAction * eventAction(
1351 static_cast< const CexmcEvent 1350 static_cast< const CexmcEventAction * >( userEventAction ) );
1352 if ( ! eventAction ) 1351 if ( ! eventAction )
1353 throw CexmcException( CexmcWeirdExcep 1352 throw CexmcException( CexmcWeirdException );
1354 1353
1355 CexmcEventAction * theEventAction( 1354 CexmcEventAction * theEventAction( const_cast< CexmcEventAction * >(
1356 1355 eventAction ) );
1357 theEventAction->BeamParticleChangeHook(); 1356 theEventAction->BeamParticleChangeHook();
1358 } 1357 }
1359 1358
1360 1359
1361 void CexmcRunManager::SetupConstructionHook( 1360 void CexmcRunManager::SetupConstructionHook( void )
1362 { 1361 {
1363 #ifdef CEXMC_USE_PERSISTENCY 1362 #ifdef CEXMC_USE_PERSISTENCY
1364 /* save gdml file */ 1363 /* save gdml file */
1365 G4String cmd( "" ); 1364 G4String cmd( "" );
1366 CexmcExceptionType exceptionType( CexmcS 1365 CexmcExceptionType exceptionType( CexmcSystemException );
1367 1366
1368 if ( zipGdmlFile ) 1367 if ( zipGdmlFile )
1369 { 1368 {
1370 if ( ProjectIsRead() ) 1369 if ( ProjectIsRead() )
1371 { 1370 {
1372 cmd = G4String( "bzip2 " ) + proj 1371 cmd = G4String( "bzip2 " ) + projectsDir + "/" + rProject +
1373 1372 gdmlFileExtension;
1374 } 1373 }
1375 else 1374 else
1376 { 1375 {
1377 if ( ProjectIsSaved() ) 1376 if ( ProjectIsSaved() )
1378 cmd = G4String( "bzip2 -c " ) 1377 cmd = G4String( "bzip2 -c " ) + gdmlFileName + " > " +
1379 projectsDir + "/" + p 1378 projectsDir + "/" + projectId + gdmlbz2FileExtension;
1380 } 1379 }
1381 exceptionType = CexmcFileCompressExce 1380 exceptionType = CexmcFileCompressException;
1382 } 1381 }
1383 else 1382 else
1384 { 1383 {
1385 if ( ! ProjectIsRead() && ProjectIsSa 1384 if ( ! ProjectIsRead() && ProjectIsSaved() )
1386 cmd = G4String( "cp " ) + gdmlFil 1385 cmd = G4String( "cp " ) + gdmlFileName + " " + projectsDir + "/" +
1387 1386 projectId + gdmlFileExtension;
1388 /* else already saved in ReadPreinitP 1387 /* else already saved in ReadPreinitProjectData() */
1389 } 1388 }
1390 1389
1391 if ( ! cmd.empty() && system( cmd ) != 0 1390 if ( ! cmd.empty() && system( cmd ) != 0 )
1392 throw CexmcException( exceptionType ) 1391 throw CexmcException( exceptionType );
1393 #endif 1392 #endif
1394 } 1393 }
1395 1394
1396 1395