Commit 7a7b0151 authored by Frank Siegert's avatar Frank Siegert
Browse files

merge r19796 from trunk: partonic hadron decays vertex offset fixed.

parent 6b60da24
......@@ -241,15 +241,20 @@ bool Perturbative_Interface::FillBlobs(ATOOLS::Blob_List *blobs)
sblob->SetStatus(blob_status::needs_showers);
sblob->SetId();
if (p_shower->On()) {
if (!p_hd)
if (!p_hd) {
for (int i(0);i<p_hard->NInP();++i)
sblob->AddToOutParticles(p_hard->InParticle(i));
for (size_t j(0);j<blobs->size();++j) {
Blob *cb((*blobs)[j]);
if (cb->Has(blob_status::needs_showers))
for (int i(0);i<cb->NOutP();++i)
if (cb->OutParticle(i)->DecayBlob()==NULL)
sblob->AddToInParticles(cb->OutParticle(i));
for (size_t j(0);j<blobs->size();++j) {
Blob *cb((*blobs)[j]);
if (cb->Has(blob_status::needs_showers))
for (int i(0);i<cb->NOutP();++i)
if (cb->OutParticle(i)->DecayBlob()==NULL)
sblob->AddToInParticles(cb->OutParticle(i));
}
}
else {
for (int i(0);i<p_hard->NOutP();++i)
sblob->AddToInParticles(p_hard->OutParticle(i));
}
}
blobs->push_back(sblob);
......
......@@ -414,9 +414,14 @@ Cluster_Amplitude* Decay_Handler_Base::ClusterConfiguration(Blob *const bl)
p_ampl->CreateLeg(cl->Mom(),cl->Flav(),cl->Col(),cl->Id());
}
}
double mu2=p_ampl->Leg(0)->Mom().Abs2();
p_ampl->SetMuF2(mu2);
p_ampl->SetKT2(mu2);
msg_Debugging()<<*p_ampl<<"\n";
while (p_ampl->Prev()) {
p_ampl=p_ampl->Prev();
p_ampl->SetMuF2(mu2);
p_ampl->SetKT2(mu2);
}
msg_Debugging()<<"}\n";
return p_ampl;
......
Markdown is supported
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment