// small struct to hold <X/X0> and <rho> of the voxel

static constexpr int NParams = 2; // number of material parameters described

float meanRho; ///< mean density, g/cm^3

float meanX2X0; ///< fraction of radiaton lenght

GPUd() MatCell() : meanRho(0.f), meanX2X0(0.f) {}

GPUdDefault() MatCell(const MatCell& src) CON_DEFAULT;

GPUd() void scale(float scale)

{

meanRho *= scale;

meanX2X0 *= scale;

}

ClassDefNV(MatCell, 1);

// small struct to hold <X/X0>, <rho> and length traversed by track in the voxel

static constexpr int NParams = 3; // number of material parameters described

float length; ///< length in material

GPUd() MatBudget() : length(0.f) {}

GPUdDefault() MatBudget(const MatBudget& src) CON_DEFAULT;

GPUd() void scale(float scale)

{

MatCell::scale(scale);

length *= scale;

}

GPUd() float getXRho() const

{

return meanRho * length;

}

GPUd() float getXRho(int signCorr) const

{

return meanRho * (signCorr < 0 ? -length : length);

}

ClassDefNV(MatBudget, 1);