Skip to content
Snippets Groups Projects

Implement detector benchmark for the barrel calorimeter for energy resolution

Merged Maria Zurek requested to merge 10-calo-resolution into master
Compare and
2 files
+ 5
0
Compare changes
  • Side-by-side
  • Inline
Files
2
@@ -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) {
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));
std::cout<<"Ethr: "<<result.size()<<endl;
return result;
};
@@ -52,6 +53,7 @@ void emcal_barrel_electrons_analysis(const char* input_fname = "sim_output/sim_e
for (const auto& i: evt)
total_edep += i.energyDeposit;
result.push_back(total_edep);
std::cout<<"Esim: "<<result.size()<<endl;
return result;
};
@@ -63,6 +65,7 @@ void emcal_barrel_electrons_analysis(const char* input_fname = "sim_output/sim_e
for (const auto& E2 : sampled)
result.push_back(E2 / E1);
}
std::cout<<"fsam: "<<result.size()<<endl;
return result;
};
Loading