#include <TGeoManager.h>
#include <TGeoMatrix.h>
+#include "AliLog.h"
#include "AliMagF.h"
#include "AliTracker.h"
#include "AliGeomManager.h"
for (Int_t i=0;i<6;i++) bparam[i]=0;
if (!gGeoManager) {
- printf("ERROR: no TGeo\n");
+ AliErrorClass("No TGeo\n");
return 0.;
}
//
// Initialize start point and direction
TGeoNode *currentnode = 0;
TGeoNode *startnode = gGeoManager->InitTrack(start, dir);
- //printf("%s length=%f\n",gGeoManager->GetPath(),length);
if (!startnode) {
AliErrorClass(Form("start point out of geometry: x %f, y %f, z %f",
start[0],start[1],start[2]));
if (nzero>3) {
// This means navigation has problems on one boundary
// Try to cross by making a small step
- printf("ERROR: cannot cross boundary\n");
+ AliErrorClass("Cannot cross boundary\n");
mparam[0] = bparam[0]/step;
mparam[1] = bparam[1];
mparam[2] = bparam[2]/step;
if (snext>=length) break;
if (!currentnode) break;
length -= snext;
- //printf("%s snext=%f length=%f\n", currentnode->GetName(),snext,length);
material = currentnode->GetVolume()->GetMedium()->GetMaterial();
lparam[0] = material->GetDensity();
lparam[1] = material->GetRadLen();
lparam[2] = material->GetA();
lparam[3] = material->GetZ();
- //printf(" %f %f %f %f\n",lparam[0],lparam[1],lparam[2],lparam[3]);
lparam[5] = lparam[3]/lparam[2];
if (material->IsMixture()) {
TGeoMixture * mixture = (TGeoMixture*)material;
}
gGeoManager->FindNextBoundaryAndStep(length, kFALSE);
snext = gGeoManager->GetStep();
- //printf("snext %f\n",snext);
}
mparam[0] = bparam[0]/step;
mparam[1] = bparam[1];