---

ExampleAnalysis.h

Go to the documentation of this file.

// ================================================
// ExampleAnalysis class descriptionn
// ================================================
//
// THIS TEXT TO BE REPLACED BY ATLAS STANDARD FORMAT
//
// Namespace Atlfast::
//
// class: ExampleAnalysis
//
// Description: 
//
//  This is an example Algorithm designed to show new users
//  of Athena/Atlfast how to access and use Atlfast objects from
//  the TES.
//
//
// ................................................................
//


#ifndef __Atlfast_example__
#define __Atlfast_example__

// STL
#include <string>
#include <vector>

// Athena/Gaudi 
#include "GaudiKernel/ISvcLocator.h"
#include "GaudiKernel/IAlgorithm.h"
#include "GaudiKernel/Algorithm.h"
#include "GaudiKernel/MsgStream.h"
#include "GaudiKernel/DataObject.h"
#include "GaudiKernel/ObjectVector.h"
#include "GaudiKernel/IHistogram1D.h"
#include "GaudiKernel/IHistogram2D.h"
#include "GaudiKernel/IHistogramSvc.h"

// Things needed for using MC truth information
#include "HepMC/GenParticle.h"
#include "AtlfastCode/HepMC_helper.h"

// Things needed to use objects through the IKinematic interface
#include "AtlfastCode/IKinematic.h"
#include "AtlfastCode/KinematicHelper.h"

// Atlfast output objects (and their collection types)
#include "AtlfastCode/ReconstructedParticle.h"
#include "AtlfastCode/ReconstructedParticleCollection.h"
#include "AtlfastCode/Cell.h"
#include "AtlfastCode/CellCollection.h"
#include "AtlfastCode/Cluster.h"
#include "AtlfastCode/ClusterCollection.h"
#include "AtlfastCode/TesIO.h"



//********************************************
//Default for parameters of this algorithm
//These can be overridden by job options.
//Descriptions are given in the class declaration
//********************************************

#define DEFAULT_rCellIsolation       0.2  

#define DEFAULT_particleLocation     "/Event/AtlfastIsolatedElectrons" 
#define DEFAULT_cellLocation         "/Event/AtlfastCells" 
#define DEFAULT_clusterLocation      "/Event/AtlfastClusters" 

#define DEFAULT_histStart            1000



//***************************************************************
//                 ExampleAnalysis class declaration
//
// An "Athena algorithm" is something which gets scheduled and controlled
// by the framework. In its simplest definition it is something which
//  - gets called for each event 
//  - can get anything it wants out of the Transient Event Store"
//  - can put anything it makes into the store.
//
// It must inherit from the  Algorithm base class
//****************************************************************

namespace Atlfast {
class ExampleAnalysis : public Algorithm
{

 public:

  //-------------------------
  // Constructors/Destructors
  //
  // Athena requires that the constructor takes certain arguments
  // (and passes them directly to the constructor of the base class)
  //-------------------------
  ExampleAnalysis( const std::string& name, ISvcLocator* pSvcLocator ) ; 
  virtual ~ExampleAnalysis();


  //------------------------------------------------------
  // Methods used by Athena to run the algorithm
  //------------------------------------------------------
  StatusCode initialize() ;
  StatusCode execute() ;
  StatusCode finalize() ;



 private:

  //---------------------------------------------------
  // Member variables
  //---------------------------------------------------

  IHistogram1D* m_hist_pT;

  IHistogram1D* m_hist_isolation;

  IHistogram2D* m_hist_cellFlow;
  IHistogram2D* m_hist_clusterFlow;

  TesIO* m_tesIO;

  //---------------------------------------------------
  // Parameters of this algorithm
  //---------------------------------------------------

  // Example of a parameter to be set via jobOptions.
  double        m_rCellIsolation;     

  std::string   m_particleLocation ;
  std::string   m_cellLocation ;
  std::string   m_clusterLocation ;

  int m_histStart ;


  //-------------------------------
  // Objects used by this class to do its job
  //-------------------------------

  // This contains some trivial but oft used procedures
  // See the WWW documentation for details.
  KinematicHelper m_kinehelp ; 



  //-------------------------------
  // Private methods
  //-------------------------------

  /* To iterate through a collection of ReconstructedParticles and dump out
   * to the MsgStream specified, and also plot pT() */
  void dumpParticles( MsgStream&,                       
                      std::vector<ReconstructedParticle*>&);


  void plotIsolation( MsgStream&, 
                      std::vector<ReconstructedParticle*>&,
                      std::vector<Cell*>& );


  void plotEnergyFlow(  MsgStream&, 
                        std::vector<Cluster*>&,
                        std::vector<Cell*>& );

};

} // end of namespace bracket


#endif

---