---

ReconstructedParticleHistogramMaker.h

Go to the documentation of this file.

// ================================================
//ReconstructedParticleHistogramMaker class description
// ================================================
//
// THIS TEXT TO BE REPLACED BY ATLAS ReconstructedParticle FORMAT
//
//
// This version....
//
//
// Namespace Atlfast::
//
// class: ReconstructedParticleHistogramMaker
//
// Description: 
//
//  Algorithm which makes debug histograms for ReconstructedParticles
//
//
// ................................................................
//

#ifndef __ATLFAST_ReconstructedParticleHistogramMaker__
#define __ATLFAST_ReconstructedParticleHistogramMaker__

// 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 "AtlfastCode/ReconstructedParticle.h"
#include "AtlfastCode/TesIO.h"


//*****************************************
// Default paths in the Transient event store to get/put entities
//
//[Note:  This will change as the current system is not considered acceptable.]
//******************************************

//places to look in the TES
#define DEFAULT_inputLocation          "/Event/Electrons" 
#define DEFAULT_mcTruthLocation        "/Event/McEventCollection"

//which flavour of reconstructed particles are we looking at?
#define DEFAULT_particleType            11

//offset of histogram numbers
#define DEFAULT_histStart 0

//prefix to histogram titles
#define DEFAULT_histTitle "Electron "



 
namespace Atlfast {
class ReconstructedParticleHistogramMaker : public Algorithm
{

 public:

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

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


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

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



 private:
  TesIO* m_tesIO;

  //-------------------------------
  // Private methods
  //-------------------------------
   void book(
             std::vector<IHistogram1D*> & start, 
             const std::string title, 
             const int nbins, 
             const double xmin, 
             const double xmax,
             const double xminDiff,
             const double xmaxDiff
             );

   void fill(
             std::vector<IHistogram1D*> & start, 
             const double rec, 
             const double tru
             );

  // 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_inputLocation ;
  std::string m_mcTruthLocation;


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

  // Each quantity to be histogrammed has a vector[4] of histograms
  // these correspond to
  // [0] = Rec = reconstructed
  // [1] = Tru = true
  // [2] = Dif = difference, true-rec
  // [3] = Res = residual, (true-rec)/true
  
  // Multiplicity
  std::vector<IHistogram1D*> m_h_multiplicity ;
  // Energy
  std::vector<IHistogram1D*> m_h_energy ;
  // Pt
  std::vector<IHistogram1D*> m_h_pt ; 
  // eta
  std::vector<IHistogram1D*> m_h_eta ;
  // phi
  std::vector<IHistogram1D*> m_h_phi ;
  // theta
  std::vector<IHistogram1D*> m_h_theta ;

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

  // particle type to look at
  int m_particleType ;

  // prefix for histogram titles
  std::string m_histTitle ;

};


} // end of namespace bracket

#endif

---