Skip to content
Snippets Groups Projects
Commit d2ba76e7 authored by Elidiano Tronchin's avatar Elidiano Tronchin
Browse files

* FIxing bug in Calorimeter Scorer

	- The Hit Level vector was not cleared at the begining of event
parent d2618c84
No related branches found
No related tags found
No related merge requests found
...@@ -61,11 +61,10 @@ G4bool PS_Calorimeter::ProcessHits(G4Step* aStep, G4TouchableHistory*){ ...@@ -61,11 +61,10 @@ G4bool PS_Calorimeter::ProcessHits(G4Step* aStep, G4TouchableHistory*){
static unsigned int mysize = m_NestingLevel.size(); static unsigned int mysize = m_NestingLevel.size();
t_Energy = aStep->GetTotalEnergyDeposit(); t_Energy = aStep->GetTotalEnergyDeposit();
t_Time = aStep->GetPreStepPoint()->GetGlobalTime(); t_Time = aStep->GetPreStepPoint()->GetGlobalTime();
t_Level.clear();
for(unsigned int i = 0 ; i < mysize ; i++){ for(unsigned int i = 0 ; i < mysize ; i++){
t_Level.push_back(aStep->GetPreStepPoint()->GetTouchableHandle()->GetCopyNumber(m_NestingLevel[i])); t_Level.push_back(aStep->GetPreStepPoint()->GetTouchableHandle()->GetCopyNumber(m_NestingLevel[i]));
} }
// Check if the particle has interact before, if yes, add up the energies. // Check if the particle has interact before, if yes, add up the energies.
vector<CalorimeterData>::iterator it; vector<CalorimeterData>::iterator it;
it = m_Data.find(CalorimeterData::CalculateIndex(t_Level)); it = m_Data.find(CalorimeterData::CalculateIndex(t_Level));
......
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment