CbmFlesHistosTools.cxx

Go to the documentation of this file.

#include "CbmFlesHistosTools.h"
 
#include <cmath>
//#include <vector>
 
std::vector<double> GenerateLogBinArray(uint32_t uNbDecadesLog,
                                        uint32_t uNbStepsDecade,
                                        uint32_t uNbSubStepsInStep,
                                        uint32_t& uNbBinsLog) {
  uNbBinsLog =
    uNbStepsDecade + uNbStepsDecade * uNbSubStepsInStep * uNbDecadesLog;
 
  uint32_t uArrayLength = uNbBinsLog + 1;
  double dBinsLog[uArrayLength];
  for (uint32_t uSubU = 0; uSubU < uNbStepsDecade; uSubU++) {
    dBinsLog[uSubU] = 0.1 * (1 + uSubU);
  }
 
  double dSubstepSize = 1.0 / uNbSubStepsInStep;
  for (uint32_t uDecade = 0; uDecade < uNbDecadesLog; uDecade++) {
    double dBase = std::pow(10, uDecade);
    uint32_t uDecadeIdx =
      uNbStepsDecade + uDecade * uNbStepsDecade * uNbSubStepsInStep;
    for (uint32_t uStep = 0; uStep < uNbStepsDecade; uStep++) {
      uint32_t uStepIdx = uDecadeIdx + uStep * uNbSubStepsInStep;
      for (uint32_t uSubStep = 0; uSubStep < uNbSubStepsInStep; uSubStep++) {
        dBinsLog[uStepIdx + uSubStep] =
          dBase * (1 + uStep) + dBase * dSubstepSize * uSubStep;
      }  // for( uint32_t uSubStep = 0; uSubStep < uNbSubStepsInStep; uSubStep++ )
    }    // for( uint32_t uStep = 0; uStep < uNbStepsDecade; uStep++ )
  }      // for( uint32_t uDecade = 0; uDecade < uNbDecadesLog; uDecade ++)
  dBinsLog[uNbBinsLog] = std::pow(10, uNbDecadesLog);
 
  std::vector<double> dBinsLogVect;
 
  for (uint32_t i = 0; i < uArrayLength; ++i) {
    dBinsLogVect.push_back(dBinsLog[i]);
  }  // for( uint32_t i = 0; i < uArrayLength; ++i )
 
  return dBinsLogVect;
}