Skip to content
Snippets Groups Projects

Implement detector benchmark for the barrel calorimeter for energy resolution

Compare and
2 files
+ 5
0
Compare changes
  • Side-by-side
  • Inline

Files

@@ -39,6 +39,7 @@ void emcal_barrel_electrons_analysis(const char* input_fname = "sim_output/sim_e
@@ -39,6 +39,7 @@ void emcal_barrel_electrons_analysis(const char* input_fname = "sim_output/sim_e
auto Ethr = [](std::vector<dd4pod::Geant4ParticleData> const& input) {
auto Ethr = [](std::vector<dd4pod::Geant4ParticleData> const& input) {
std::vector<double> result;
std::vector<double> result;
result.push_back(TMath::Sqrt(input[2].psx*input[2].psx + input[2].psy*input[2].psy + input[2].psz*input[2].psz + input[2].mass*input[2].mass));
result.push_back(TMath::Sqrt(input[2].psx*input[2].psx + input[2].psy*input[2].psy + input[2].psz*input[2].psz + input[2].mass*input[2].mass));
 
std::cout<<"Ethr: "<<result.size()<<endl;
return result;
return result;
};
};
@@ -52,6 +53,7 @@ void emcal_barrel_electrons_analysis(const char* input_fname = "sim_output/sim_e
@@ -52,6 +53,7 @@ void emcal_barrel_electrons_analysis(const char* input_fname = "sim_output/sim_e
for (const auto& i: evt)
for (const auto& i: evt)
total_edep += i.energyDeposit;
total_edep += i.energyDeposit;
result.push_back(total_edep);
result.push_back(total_edep);
 
std::cout<<"Esim: "<<result.size()<<endl;
return result;
return result;
};
};
@@ -63,6 +65,7 @@ void emcal_barrel_electrons_analysis(const char* input_fname = "sim_output/sim_e
@@ -63,6 +65,7 @@ void emcal_barrel_electrons_analysis(const char* input_fname = "sim_output/sim_e
for (const auto& E2 : sampled)
for (const auto& E2 : sampled)
result.push_back(E2 / E1);
result.push_back(E2 / E1);
}
}
 
std::cout<<"fsam: "<<result.size()<<endl;
return result;
return result;
};
};
Loading