}
if (fQuench) {
- fPythia->InitQuenching(0., 0.1, 1., 8., 0);
+ fPythia->InitQuenching(0., 0.05, 1.45, 8., 1);
}
}
// Run quenching routine
//
fPythia->Quench();
+// fPythia->Pylist(1);
//
// Switch fragmentation on
fPythia->SetMSTJ(1, 1);
//
// .. and perform fragmentation
fPythia->Pyexec();
+// fPythia->Pylist(1);
}
p[1] = iparticle->Py();
p[2] = iparticle->Pz();
p[3] = iparticle->Energy();
+
origin[0] = fVertex[0]+iparticle->Vx()/10; // [cm]
origin[1] = fVertex[1]+iparticle->Vy()/10; // [cm]