Float_t polar[3] = {0., 0., 0.};
Int_t np = fDecayer->ImportParticles(particles);
fNpartProd += (np-1);
- Int_t nt;
+ Int_t nt = 0;
for (Int_t i = 1; i < np; i++)
{
TParticle* iParticle = (TParticle *) particles->At(i);
// 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;
if ( key < 1 || key > 200) {
printf("ERROR in THijing::SetPARJ(key,parm):\n");
printf(" key=%i is out of range [1..200]\n",key);
+ } else {
+ LUDAT1_HIJING.parj[key-1] = parm;
}
-
- LUDAT1_HIJING.parj[key-1] = parm;
}
if ( key < 1 || key > 200) {
printf("ERROR in THijing::SetMSTJ(key,parm):\n");
printf(" key=%i is out of range [1..200]\n",key);
+ } else {
+ LUDAT1_HIJING.mstj[key-1] = parm;
}
-
- LUDAT1_HIJING.mstj[key-1] = parm;
}