Skip to content
Snippets Groups Projects

Draft: Statuses

Closed Andrii Verbytskyi requested to merge averbyts/athena:statuses into main
Files
5
@@ -377,10 +377,12 @@ int ALFA_BeamTransport::TransportSelectedParticle(HepMC::GenEvent& evt, int evt_
ATH_MSG_INFO("Add transproted particle into HepMC event record Beam 1");
//Add Data for HepMC Collection
for(int i=0;i<(int)PosAtRP1.size();i++){//Beam1
//The factor 1000 comes from the fact that HepMC saves length in mm
HepMC::FourVector PositionVectorRP1 = HepMC::FourVector(PosAtRP1.at(i).x()*1000.,PosAtRP1.at(i).y()*1000.,PosAtRP1.at(i).z()*1000.,0.*1000.);
HepMC::FourVector MomentumVectorRP1 = HepMC::FourVector(MomAtPR1.at(i).x(),MomAtPR1.at(i).y(),MomAtPR1.at(i).z(),EnergyRP1.at(i));
for(size_t i = 0; i < PosAtRP1.size(); i++){//Beam1
//The factor 1000 comes from the fact that HepMC saves length in mm
HepMC::FourVector PositionVectorRP1 = HepMC::FourVector(PosAtRP1.at(i).x()*1000.,PosAtRP1.at(i).y()*1000.,PosAtRP1.at(i).z()*1000.,0.*1000.);
HepMC::FourVector MomentumVectorRP1 = HepMC::FourVector(MomAtPR1.at(i).x(),MomAtPR1.at(i).y(),MomAtPR1.at(i).z(),EnergyRP1.at(i));
HepMC::GenVertexPtr VertexRP1 = HepMC::newGenVertexPtr(PositionVectorRP1);
HepMC::GenParticlePtr ParticleRP1 = HepMC::newGenParticlePtr(MomentumVectorRP1,2212,1); //save the transported particle with status code 1 (added 120124) preview was 201
Loading