// Get event vertex
//
- TParticle * iparticle = (TParticle *) fParticles.At(0);
fVertex[0] = origin0[0];
fVertex[1] = origin0[1];
fVertex[2] = origin0[2];
//
// First select parent particles
//
-
+ TParticle * iparticle = 0;
for (i = 0; i < np; i++) {
iparticle = (TParticle *) fParticles.At(i);
Float_t* b = new Float_t[kMax];
Float_t* si1 = new Float_t[kMax];
Float_t* si2 = new Float_t[kMax];
-
+ for (i = 0; i < kMax; i++){
+ b[i] = 0.;
+ si1[i] = 0.;
+ si2[i] = 0.;
+ }
+
for (i = 0; i < kMax; i++)
{
Float_t xb = bMin+i*kdib;