56 idCurrentLoop=nMaxLoops;
59 if (bCompareExp){nMaxLoops=-1;};
60 fileExperimentalData=FileExperimentalData;
63 if (bCompareExp && fileExperimentalData!=
"")
65 ML2ExpVoxels=
new CML2ExpVoxels(bCompareExp, saveEvents, seed, FileExperimentalData, FileExperimentalDataOut);
69 std::cout <<
"I don't have any convergence criteria set, I'll do " << nMaxLoops <<
" loop(s) for each rotation"<<
G4endl;
76 maxNumberOfEvents=maxNumEvents;
82 {
delete ML2ExpVoxels;}
91 {ML2ExpVoxels->
add(aStep);}
99 bStopRun=convergenceCriteria();
105 if (idCurrentLoop==0)
112 G4bool CML2Convergence::convergenceCriteria()
115 G4int nEventsAccumulated=0;
125 std::cout <<
"\n ++++++++++++++++++++ \n";
126 std::cout <<
"current geometry: " << nGeometry;
127 std::cout <<
"\nNumber of events accumulated in the current geometry:"<<
128 nEventsAccumulated<<
"\nNumber of events to be accumulated:" <<
129 maxNumberOfEvents<<
"\n -------------------------\n";
void add(const G4Step *aStep)
void resetNEventsInVoxels()
G4double GetTotalEnergyDeposit() const
G4int getMaxNumberOfEvents()
void setRecycling(int recycling)
void add(G4ThreeVector pos, G4double depEnergy, G4double density)