function printDistributionStatistics(pbeDefinition)

    disp('--- DISTRIBUTION STATISTICS ---')
    % precalculations
    positionMatrix = pbeDefinition.seeding.disD.positionMatrix;
    numberVector = 1/pbeDefinition.seeding.nSample * ...
        pbeDefinition.seeding.fhandle(positionMatrix) ./ ...
        pbeDefinition.seeding.disD.probabilityVector;
        
    disp('- integral properties')
    totalNumber = sum(numberVector)
    
    totalVolume = (computeMeasureDec(pbeDefinition.dec, 'hC', ...
          positionMatrix, 'volume', [], [], -1)) * numberVector'
    totalSurface = computeMeasureDec(pbeDefinition.dec, 'hC', ...
        positionMatrix, 'surface', [], [], -1) * numberVector'
    totalMeanWidth = computeMeasureDec(pbeDefinition.dec, 'hC', ...
        positionMatrix, 'meanWidth', [], [], -1) * numberVector'

    disp('- plotting 1D distributions')
%     % particle number over mean width
%     binLower = 0; binUpper = 1e-17; binSpacing = 1e-18;
% %     binLower = 0; binUpper = 4e-6; binSpacing = 1e-7;
%     % bounds in [lower bounds; upper bounds]
%     boundDefinitionMatrix = [binLower:binSpacing:(binUpper - binSpacing);...
%         (binLower + binSpacing):binSpacing:binUpper]
%     % pivot position
%     pivotPositionVector = mean(boundDefinitionMatrix, 1)
%     nPivot = length(pivotPositionVector);
%     
%     % handles for integration
%     disp('..measure y')
% 	yVector = numberVector;
%     size(yVector)
%     disp('..measure x')
%     xVector = computeMeasureDec(pbeDefinition.dec, 'hC', ...
%         positionMatrix, 'volume', [], [], -1);
%     size(xVector)
%     disp('..measure finished')
%     
%     %   integrand looks like:
%     %     yProperty*(xProperty > lowerBound)*(xProperty <= upperBound)
%     %   Now consider that yProperty is the vector as given by
%     %   yAxisHandle. While yProperty is a row vector, the boundary
%     %   filters can be written as matrices that are multiplied from the
%     %   right:
%     %     yProperty*boundaryFilterMatrix
%     %   so that the result of the integral will become a column vector.
%     % Hence, we construct the boundary filter, first.
%     boundaryFilter = ...
%             bsxfun(@gt, xVector, boundDefinitionMatrix(1,:)') .* ...
%             bsxfun(@le, xVector, boundDefinitionMatrix(2,:)');
%     size(boundaryFilter)
%     % integrand
%     integrandVector = boundaryFilter .* ...
%                      (ones(size(boundaryFilter, 1), 1)*yVector) * ...
%                      numberVector';
%     size(integrandVector)
%     figure, plot(pivotPositionVector, integrandVector, 'bx-')    
% %     title('particle number distr. density over mean width')
%     ylabel('number distr. density in #/m^3')
%     xlabel('mean width in m^3')
%     % % set(gca,'xscale', 'log')

end