---

StandardHistogramMaker.h

Go to the documentation of this file.

// ================================================
//StandardHistogramMaker class description
// ================================================
//
// THIS TEXT TO BE REPLACED BY ATLAS STANDARD FORMAT
//
//
// This version....
//
//
// Namespace Atlfast::
//
// class: StandardHistogramMaker
//
// Description: 
//
//  Algorithm which makes some standard Atlfast debug histograms
//
//
// ................................................................
//

#ifndef ATLFAST_STANDARDHISTOGRAMMAKER_H
#define ATLFAST_STANDARDHISTOGRAMMAKER_H

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

// Gaudi 
#include "GaudiKernel/ISvcLocator.h"
#include "GaudiKernel/Algorithm.h"
#include "GaudiKernel/MsgStream.h"
#include "GaudiKernel/DataObject.h"

class IHistogram1D;  // forward declaration

// Other
#include "CLHEP/Vector/LorentzVector.h"
#include "HepMC/GenEvent.h"
#include "HepMC/GenParticle.h"

// Atlfast
#include "AtlfastEvent/ReconstructedParticle.h"
#include "AtlfastUtils/TesIO.h"

//***************************************************************
//                 StandardHistogramMaker class declaration
//
// A "Gaudi algorithm" is something which gets scheduled and controlled
// by the framework. In its simplest definition it is something which
//  - gets clled 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 {
  using std::string;
class StandardHistogramMaker : public Algorithm
{

 public:

  //-------------------------
  // Constructors/Destructors
  //
  // Gaudi requires that the constructor takes certain arguments
  // (and passes them directly to the constructor of the base class)
  //-------------------------

  StandardHistogramMaker( const std::string& name, ISvcLocator* pSvcLocator ) ; 
  ~StandardHistogramMaker();


  //------------------------------------------------------
  // Methods used by Gaudi to run the algorithm
  //------------------------------------------------------

  StatusCode initialize() ;
  StatusCode execute() ;
  StatusCode finalize() ;



 private:


  //-------------------------------
  // Private methods
  //-------------------------------
 
  void bookElectronHistograms();
  void bookPhotonHistograms();
  void bookIsolatedElectronHistograms();
  void bookIsolatedPhotonHistograms();
  void bookCellHistograms();
  void bookClusterHistograms();

  void fillElectronHistograms() ;
  void fillPhotonHistograms() ;
  void fillIsolatedElectronHistograms();
  void fillIsolatedPhotonHistograms();
  void fillCellHistograms();
  void fillClusterHistograms();

  TesIO* m_tesIO;


  // the following data members tell the
  // algorithm where to retrieve its input particles from
  // in the TES. These are written by the job Options service.

  std::string m_electronLocation ;
  std::string m_isolatedElectronLocation;
  std::string m_photonLocation ;
  std::string m_isolatedPhotonLocation;
  std::string m_cellLocation;
  std::string m_clusterLocation;
  std::string m_mcTruthLocation;


  // the following data members tell the
  // algorithm whether or not certain sets of histograms
  // are required

  bool m_doElectronHistograms ;
  bool m_doIsolatedElectronHistograms ;
  bool m_doPhotonHistograms ;
  bool m_doIsolatedPhotonHistograms ;
  bool m_doCellHistograms ;
  bool m_doClusterHistograms ;

  // The following are the actual histograms
  // htp has added his own naming convention here
  // m_h_ means a histogram!

  // electron histograms
  IHistogram1D* m_h_electronMultiplicity;
  IHistogram1D* m_h_electronPt;

  // isolated electron histograms
  IHistogram1D* m_h_isolatedElectronMultiplicity;
  IHistogram1D* m_h_isolatedElectronPt;

  // photon histograms
  IHistogram1D* m_h_photonMultiplicity;
  IHistogram1D* m_h_photonPt;

  // isolated photon histograms
  IHistogram1D* m_h_isolatedPhotonMultiplicity;
  IHistogram1D* m_h_isolatedPhotonPt;

  // cell histograms
  IHistogram1D* m_h_cellMultiplicity;
  IHistogram1D* m_h_cellPt;

  // cluster histograms
  IHistogram1D* m_h_clusterMultiplicity;
  IHistogram1D* m_h_clusterPt;


  // histogram numbers, start at...
  int m_histStart;
  // and current count
  int m_nHist;
};


} // end of namespace bracket

#endif

---