#include <TPDGCode.h>
#include <TSystem.h>
+#include "AliLog.h"
#include "AliGenHaloProtvino.h"
#include "AliRun.h"
const Float_t kFlux = 1.e11 / 25.e-9; // protons/s
Float_t pFlux[5] = {0.2, 0.2, 0.3, 0.3, 1.0};
+ Int_t ncols = 0;
if (fRunPeriod < 5) {
//
// Ring 1
for (i = 0; i < fGPASize; i++)
{
- fscanf(file, "%f %f %f %f %f %f", &z, &p[0], &p[1], &p[2] , &p[3], &p[4]);
+ ncols = fscanf(file, "%f %f %f %f %f %f", &z, &p[0], &p[1], &p[2] , &p[3], &p[4]);
+ if (ncols<0) break;
+
fG1[i] = p[fRunPeriod];
if (i > 0) {
//
for (i = 0; i < fGPASize; i++)
{
- fscanf(file, "%f %f %f %f %f %f", &z, &p[0], &p[1], &p[2] , &p[3], &p[4]);
+ ncols = fscanf(file, "%f %f %f %f %f %f", &z, &p[0], &p[1], &p[2] , &p[3], &p[4]);
+ if (ncols<0) break;
+
fG2[i] = p[fRunPeriod];
if (i > 0) {
fZ2[i] = fZ2[i-1] + z;
} else {
for (i = 0; i < fGPASize; i++)
{
- fscanf(file, "%f %e %e %e %e %e", &z, &p[0], &p[1], &p[2], &p[3], &p[4]);
+ ncols = fscanf(file, "%f %e %e %e %e %e", &z, &p[0], &p[1], &p[2], &p[3], &p[4]);
+ if (ncols<0) break;
+
z /= 1000.;
fG1[i] = p[4] * kCrossSection * kFlux; // 1/m/s
// 1/3 of nominal intensity at startup
//
// Mean time between interactions
//
- Float_t dT = fTimePerEvent/nInt; // sec
+
+ Float_t dT = 0.; // sec
+ if (nInt > 0)
+ dT = fTimePerEvent/nInt;
Float_t t = 0; // sec
//