Skip to content
Snippets Groups Projects
Commit 2870a72a authored by Kong Tu's avatar Kong Tu
Browse files

run dvcs-d analysis

parent 6ce24ea9
Branches
No related tags found
1 merge request!100Draft: Diffractive phi benchmark development
...@@ -154,7 +154,7 @@ int dvcs_d_analysis(const std::string& config_name) ...@@ -154,7 +154,7 @@ int dvcs_d_analysis(const std::string& config_name)
h_EtaPhi_proton_REC->Write(); h_EtaPhi_proton_REC->Write();
h_PxPy_proton_REC->Write(); h_PxPy_proton_REC->Write();
h_PxPy_neutron_REC->Write(); h_PxPy_neutron_REC->Write();
h_Momentum_neutron_REC->Wirte(); h_Momentum_neutron_REC->Write();
h_t_REC->Write(); h_t_REC->Write();
h_t_REC_A->Write(); h_t_REC_A->Write();
......
...@@ -189,7 +189,7 @@ auto findScatNeutronMC(const std::vector<dd4pod::Geant4ParticleData>& parts){ ...@@ -189,7 +189,7 @@ auto findScatNeutronMC(const std::vector<dd4pod::Geant4ParticleData>& parts){
std::vector<ROOT::Math::PxPyPzMVector> momenta; std::vector<ROOT::Math::PxPyPzMVector> momenta;
for(auto& i1 : parts){ for(auto& i1 : parts){
if(i1.genStatus==1&&i1.pdgID==2112){ if(i1.genStatus==1&&i1.pdgID==2112){
momenta.push_back(ROOT::Math::PxPyPzMVector{i1.p.x,i1.p.y,i1.p.z,i1.mass}); momenta.push_back(ROOT::Math::PxPyPzMVector{i1.ps.x,i1.ps.y,i1.ps.z,i1.mass});
} }
} }
return momenta; return momenta;
......
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Please register or to comment