Aida2Root.cpp

Go to the documentation of this file.

/***********************************************************************************\
* (c) Copyright 1998-2019 CERN for the benefit of the LHCb and ATLAS collaborations *
*                                                                                   *
* This software is distributed under the terms of the Apache version 2 licence,     *
* copied verbatim in the file "LICENSE".                                            *
*                                                                                   *
* In applying this licence, CERN does not waive the privileges and immunities       *
* granted to it by virtue of its status as an Intergovernmental Organization        *
* or submit itself to any jurisdiction.                                             *
\***********************************************************************************/
#ifdef __ICC
// disable icc remark #2259: non-pointer conversion from "X" to "Y" may lose significant bits
//   TODO: To be removed, since it comes from ROOT TMathBase.h
#  pragma warning( disable : 2259 )
#endif
 
#include "GaudiUtils/Aida2ROOT.h"
#include "AIDA/IHistogram1D.h"
#include "AIDA/IHistogram2D.h"
#include "AIDA/IHistogram3D.h"
#include "AIDA/IProfile1D.h"
#include "AIDA/IProfile2D.h"
#include "GaudiAlg/GaudiHistoAlg.h"
#include "GaudiUtils/HistoStats.h"
#include "TH1D.h"
#include "TH2D.h"
#include "TH3D.h"
#include "TProfile.h"
#include "TProfile2D.h"
#include <fmt/format.h>
 
// ============================================================================
class Aida2Root : public GaudiHistoAlg {
public:
  StatusCode execute() override { return StatusCode::SUCCESS; };
  StatusCode finalize() override;
 
public:
  // standard constructor
  Aida2Root( const std::string& name, ISvcLocator* pSvc ) : GaudiHistoAlg( name, pSvc ) {
    setProperty( "PropertiesPrint", true ).ignore();
  }
 
private:
  Gaudi::Property<std::vector<std::string>> m_1Ds{
      this,
      "Histos1D",
      { "SimpleHistos/Gaussian mean=0, sigma=1", "SimpleHistos/101", "SimpleHistos/102", "SimpleHistos/1111",
        "SimpleHistos/test1", "SimpleHistos/subdir2/bino", "SimpleHistos/subdir1/bino", "SimpleHistos/poisson" },
      "list of 1D-histograms" };
 
  Gaudi::Property<std::vector<std::string>> m_2Ds{ this,
                                                   "Histos2D",
                                                   { "SimpleHistos/Gaussian V Flat", "SimpleHistos/Exponential V Flat",
                                                     "SimpleHistos/binVpois", "SimpleHistos/expoVpois" },
                                                   "list of 2D-histograms" };
 
  Gaudi::Property<std::vector<std::string>> m_3Ds{
      this, "Histos3D", { "SimpleHistos/3D plot AutoID", "SimpleHistos/3d" }, "list of 3D-histograms" };
 
  Gaudi::Property<std::vector<std::string>> m_1Ps{
      this, "Profs1D", { "SimpleHistos/Expo V Gauss 1DProf" }, "list of 1D-profiles" };
 
  Gaudi::Property<std::vector<std::string>> m_2Ps{
      this, "Profs2D", { "SimpleHistos/321", "SimpleHistos/2dprof" }, "list of 2D-profiles" };
};
// ============================================================================
// ============================================================================
DECLARE_COMPONENT( Aida2Root )
 
// ============================================================================
// ============================================================================
StatusCode Aida2Root::finalize() {
 
  always() << "Get the native ROOT representation of histograms!" << endmsg;
 
  { // loop over all 1D-histograms
    for ( auto& path : m_1Ds ) {
      AIDA::IHistogram1D* aida = 0;
      StatusCode          sc   = histoSvc()->retrieveObject( path, aida );
      if ( sc.isFailure() || 0 == aida ) { return Error( "Unable to retrieve 1D-histogram '" + ( path ) + "'" ); }
      TH1D* root = Gaudi::Utils::Aida2ROOT::aida2root( aida );
      if ( 0 == root ) { return Error( "Unable to convert to ROOT the 1D-histogram '" + ( path ) + "'" ); }
      info() << "The native ROOT printout for 1D-histogram '" << ( path ) << "':" << endmsg;
      root->Print();
 
      info() << " |  Compare       | AIDA/HistoStats |     ROOT/TH1    |      Delta      | " << endmsg;
      auto print = []( double aida, double root, std::string_view name ) {
        return fmt::format( " | {:<14.14s} | {:>15.8g} | {:< 15.8g} | {:^ 15.8g} | ", name, aida, root,
                            ( aida - root ) );
      };
      info() << print( Gaudi::Utils::HistoStats::mean( aida ), root->GetMean(), "'mean'" ) << endmsg;
      info() << print( Gaudi::Utils::HistoStats::meanErr( aida ), root->GetMeanError(), "'meanErr'" ) << endmsg;
      info() << print( Gaudi::Utils::HistoStats::rms( aida ), root->GetRMS(), "'rms'" ) << endmsg;
      info() << print( Gaudi::Utils::HistoStats::rmsErr( aida ), root->GetRMSError(), "'rmsErr'" ) << endmsg;
      info() << print( Gaudi::Utils::HistoStats::skewness( aida ), root->GetSkewness(), "'skewness'" ) << endmsg;
      info() << print( Gaudi::Utils::HistoStats::skewnessErr( aida ), root->GetSkewness( 11 ), "'skewnessErr'" )
             << endmsg;
      info() << print( Gaudi::Utils::HistoStats::kurtosis( aida ), root->GetKurtosis(), "'kurtosis'" ) << endmsg;
      info() << print( Gaudi::Utils::HistoStats::kurtosisErr( aida ), root->GetKurtosis( 11 ), "'kurtosisErr'" )
             << endmsg;
    }
  }
 
  { // loop over all 2D-histograms
    for ( auto& path : m_2Ds ) {
      AIDA::IHistogram2D* aida = 0;
      StatusCode          sc   = histoSvc()->retrieveObject( path, aida );
      if ( sc.isFailure() || 0 == aida ) { return Error( "Unable to retrieve 2D-histogram '" + ( path ) + "'" ); }
      TH2D* root = Gaudi::Utils::Aida2ROOT::aida2root( aida );
      if ( 0 == root ) { return Error( "Unable to convert to ROOT the 2D-histogram '" + ( path ) + "'" ); }
      info() << "The native ROOT printout for 2D-histogram '" << ( path ) << "':" << endmsg;
      root->Print();
    }
  }
 
  { // loop over all 3D-histograms
    for ( auto& path : m_3Ds ) {
      AIDA::IHistogram3D* aida = 0;
      StatusCode          sc   = histoSvc()->retrieveObject( path, aida );
      if ( sc.isFailure() || 0 == aida ) { return Error( "Unable to retrieve 3D-histogram '" + ( path ) + "'" ); }
      TH3D* root = Gaudi::Utils::Aida2ROOT::aida2root( aida );
      if ( 0 == root ) { return Error( "Unable to convert to ROOT the 3D-histogram '" + ( path ) + "'" ); }
      info() << "The native ROOT printout for 3D-histogram '" << ( path ) << "':" << endmsg;
      root->Print();
    }
  }
 
  { // loop over all 1D-profiles
    for ( auto& path : m_1Ps ) {
      AIDA::IProfile1D* aida = 0;
      StatusCode        sc   = histoSvc()->retrieveObject( path, aida );
      if ( sc.isFailure() || 0 == aida ) { return Error( "Unable to retrieve 1D-profile '" + ( path ) + "'" ); }
      TProfile* root = Gaudi::Utils::Aida2ROOT::aida2root( aida );
      if ( 0 == root ) { return Error( "Unable to convert to ROOT the 1D-profile '" + ( path ) + "'" ); }
      info() << "The native ROOT printout for 1D-profile '" << ( path ) << "':" << endmsg;
      root->Print();
    }
  }
 
  { // loop over all 2D-profiles
    for ( auto& path : m_2Ps ) {
      AIDA::IProfile2D* aida = 0;
      StatusCode        sc   = histoSvc()->retrieveObject( path, aida );
      if ( sc.isFailure() || 0 == aida ) { Error( "Unable to retrieve 2D-profile '" + ( path ) + "'" ).ignore(); }
      TProfile2D* root = Gaudi::Utils::Aida2ROOT::aida2root( aida );
      if ( 0 == root ) { Error( "Unable to convert to ROOT the 2D-profile '" + ( path ) + "'" ).ignore(); }
      info() << "The native ROOT printout for 2D-profile '" << ( path ) << "':" << endmsg;
      root->Print();
    }
  }
 
  return GaudiHistoAlg::finalize();
}
 
// ============================================================================
// The END
// ============================================================================