UpdateTProfile.cc

Go to the documentation of this file.

#include "DQM/SiStripCommon/interface/UpdateTProfile.h"
#include "TProfile.h"
#include <iostream>
#include <cmath>

// -----------------------------------------------------------------------------
UpdateTProfile::UpdateTProfile() { ; }

// -----------------------------------------------------------------------------
UpdateTProfile::~UpdateTProfile() { ; }

// -----------------------------------------------------------------------------
void UpdateTProfile::setBinContents(TProfile* const prof,
                                    const uint32_t& bin,
                                    const double& num_of_entries,
                                    const double& sum_of_contents,
                                    const double& sum_of_squares) {
  //   std::cout << "[UpdateTProfile::setBinContents]" << std::endl;

  double mean = 0.;
  double spread = 0.;
  if (num_of_entries) {
    mean = sum_of_contents / num_of_entries;
    spread = sqrt(sum_of_squares / num_of_entries - mean * mean);
  }

  //   LogTrace("TEST")
  //     << "[UpdateTProfile::setBinContents]"
  //     << " bin: " << bin
  //     << " entries: " << num_of_entries
  //     << " contents: " << sum_of_contents
  //     << " squared: " << sum_of_squares
  //     << " mean: " << mean
  //     << " spread: " << spread;

  UpdateTProfile::setBinContent(prof, bin, num_of_entries, mean, spread);
}

// -----------------------------------------------------------------------------
void UpdateTProfile::setBinContent(
    TProfile* const prof, const uint32_t& bin, const double& num_of_entries, const double& mean, const double& spread) {
  //   std::cout << "[UpdateTProfile::setBinContents]" << std::endl;

  // Check histo exists
  if (!prof) {
    std::cerr << "[UpdateTProfile::setBinContents]"
              << " NULL pointer to TProfile object!" << std::endl;
    return;
  }

  // Check bin number is valid
  if (bin == 0 || bin > static_cast<uint32_t>(prof->GetNbinsX())) {
    std::cerr << "[UpdateTProfile::setBinContents]"
              << " Unexpected bin number!" << std::endl;
    return;
  }

  // Check entries are present
  if (num_of_entries <= 0.) {
    return;
  }  //@@ what about negative weights???
  double entries = num_of_entries;

  // Check error option
  const char* spread_option = "s";
  const char* default_option = "";
  //   const char* option = prof->GetErrorOption();
  //   if ( option[0] != spread_option[0] ) {
  //     std::cout << "[UpdateTProfile::setBinContents]"
  //     << " Setting error option for TProfile to 's'!" << std::endl;
  //     prof->SetErrorOption( "s" );
  //   }
  const char* error_option = prof->GetErrorOption();

  // Calculate "weight" used for SetBinError() method
  double weight;
  if (error_option[0] == spread_option[0]) {
    weight = sqrt(mean * mean * entries + spread * spread * entries);
  } else if (error_option[0] == default_option[0]) {
    weight = sqrt(mean * mean * entries + spread * spread * entries * entries);
  } else {
    std::cerr << "[UpdateTProfile::setBinContents]"
              << " Unexpected error option for TProfile!" << std::endl;
    weight = 0.;
  }

  // Set bin entries, contents and error
  prof->SetBinEntries(bin, entries);
  prof->SetBinContent(bin, mean * entries);
  prof->SetBinError(bin, weight);

  //   LogTrace("TEST")
  //     << "[UpdateTProfile::" << __func__ << "]"
  //     << " bin/entries/mean/content/error: "
  //     << bin << "/"
  //     << entries << "/"
  //     << mean << "/"
  //     << mean*entries << "/"
  //     << weight;

  // Set total number of entries
  if (bin == 1) {
    prof->SetEntries(entries);
  } else {
    prof->SetEntries(prof->GetEntries() + entries);
  }
}