EvtColAlg.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.                                             *
\***********************************************************************************/
// ============================================================================
// Include files
// ============================================================================
// GaudiKernel
// ============================================================================
#include "GaudiKernel/DataObject.h"
#include "GaudiKernel/INTupleSvc.h"
#include "GaudiKernel/IOpaqueAddress.h"
#include "GaudiKernel/IRegistry.h"
#include "GaudiKernel/IRndmGenSvc.h"
#include "GaudiKernel/RndmGenerators.h"
// ============================================================================
// GaudiAlg
// ============================================================================
#include "GaudiAlg/GaudiTupleAlg.h"
#include "GaudiAlg/Tuple.h"
#include "GaudiAlg/TupleID.h"
#include "GaudiAlg/Tuples.h"
// ============================================================================
// ============================================================================
namespace Gaudi {
  namespace Examples {
    // ========================================================================
    class EvtColAlg : public GaudiTupleAlg {
    public:
      StatusCode initialize() override {
        StatusCode sc = GaudiTupleAlg::initialize();
        if ( sc.isFailure() ) { return sc; }
        // check for random numbers service
        Assert( randSvc() != 0, "Random Service is not available!" );
        //
        return StatusCode::SUCCESS;
      };
      StatusCode execute() override;
      EvtColAlg( const std::string& name, ISvcLocator* pSvc ) : GaudiTupleAlg( name, pSvc ) {
        setProperty( "NTupleProduce", false ).ignore();
        setProperty( "NTuplePrint", false ).ignore();
        setProperty( "HistogramProduce", false ).ignore();
        setProperty( "HistogramPrint", false ).ignore();
        setProperty( "PropertiesPrint", true ).ignore();
        setProperty( "TypePrint", false ).ignore();
        setProperty( "EvtColsProduce", true ).ignore();
        setProperty( "EvtColsPrint", true ).ignore();
      }
 
    private:
      // default constructor is private
      EvtColAlg();
      // copy constructor is disabled
      EvtColAlg( const EvtColAlg& );
      // assignement operator is disabled
      EvtColAlg& operator=( const EvtColAlg& );
    };
    // ========================================================================
  } // namespace Examples
} // end of namespace Gaudi
// ============================================================================
// ============================================================================
StatusCode Gaudi::Examples::EvtColAlg::execute() {
  static int s_nEvt = 0;
  static int s_nRun = 0;
 
  if ( 1 == ++s_nEvt % 50 ) { ++s_nRun; }
 
  DataObject* event = get<DataObject>( "/Event" );
  if ( 0 == event ) { return StatusCode::FAILURE; }
  IRegistry* registry = event->registry();
  if ( 0 == registry ) { return Error( "IRegistry* point to NULL!" ); }
  IOpaqueAddress* address = registry->address();
  if ( 0 == address ) { return Error( "Address points to NULL!" ); }
 
  Rndm::Numbers gauss( randSvc(), Rndm::Gauss( 0.0, 1.0 ) );
  Rndm::Numbers flat( randSvc(), Rndm::Flat( -10.0, 10.0 ) );
  Rndm::Numbers expo( randSvc(), Rndm::Exponential( 1.0 ) );
  Rndm::Numbers breit( randSvc(), Rndm::BreitWigner( 0.0, 1.0 ) );
  Rndm::Numbers poisson( randSvc(), Rndm::Poisson( 2.0 ) );
  Rndm::Numbers binom( randSvc(), Rndm::Binomial( 8, 0.25 ) );
 
  // get the event tag collection
  Tuple tuple = evtCol( TupleID( "COL1" ), "The most trivial Event Tag Collection" );
 
  // event address (the most important information)
  tuple->column( "Address", address ).ignore();
  // put 'event' and 'run' number
  tuple->column( "evtNum", s_nEvt ).ignore();
  tuple->column( "runNum", s_nRun ).ignore();
  // put some 'data'
  tuple->column( "gauss", gauss() ).ignore();
  tuple->column( "flat", flat() ).ignore();
  tuple->column( "expo", expo() ).ignore();
  tuple->column( "breit", breit() ).ignore();
  tuple->column( "poisson", (int)poisson() ).ignore();
  tuple->column( "binom", (int)binom() ).ignore();
  tuple->column( "flag", 0 > gauss() ).ignore(); // boolean
 
  // The action of is restored again!
  return tuple->write();
}
// ============================================================================
using Gaudi::Examples::EvtColAlg;
DECLARE_COMPONENT( EvtColAlg )
// ============================================================================
// The END
// ============================================================================