char name[21];
for (j=0; j< 2; j++) {
- printf("\n A Z dens radL absL");
+ printf("\n A Z ZPos DZ dens radL absL");
printf("\n==================================================================================");
f0=f1=f2=0;
for (i=0; i< nl[j]; i++) {
f0 += dz/radl;
f1 += dz/(2.*radl)*(2.*zmin[j][i]+dz);
f2 += dz/(3.*radl)*(3.*zmin[j][i]*zmin[j][i]+3.*dz*zmin[j][i]+dz*dz);
- printf("\n %3d %14s %12.3f %12.3f %12.3f %12.3f %12.3f",
- i+1, name, a, z, dens, radl, absl);
+ printf("\n %3d %14s %12.3f %12.3f %12.3f %12.3f %12.3f %12.3f %12.3f",
+ i+1, name, a, z, (zmax[j][i]+zmin[j][i])/2., dz, dens, radl, absl);
}
//
Float_t deltaThetaB = 0.0136 * TMath::Sqrt(f0-f1*f1/f2);