---

CellHistogramMaker2.h

Go to the documentation of this file.

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

#ifndef ATLFAST_CELLHISTOGRAMMAKER2_H
#define ATLFAST_CELLHISTOGRAMMAKER2_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 "AtlfastCode/CellMaker.h"
#include "AtlfastCode/Cell.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/AtlfastCells" 
#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 "Cells "



 
namespace Atlfast {
class CellHistogramMaker2 : public Algorithm{

 public:

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

  //------------------------------------------------------
  // 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, double xmax);
  void sbook(const std::string title, const int nbins, const double xmin, double xmax);
  void fill(std::vector<IHistogram1D*> & start, const double rec, const double tru);
  void sfill(const double rec);
  
  // 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!
  
  IHistogram1D* m_h_phipart ;
  IHistogram1D* m_h_phicell ;
  IHistogram1D* m_h_etapart ;
  IHistogram1D* m_h_etacell ;
  
  IHistogram1D* m_h_phidiff_negEC ;
  IHistogram1D* m_h_etadiff_negEC ;
  
  IHistogram1D* m_h_phidiff_Bar ;
  IHistogram1D* m_h_etadiff_Bar ;
  
  IHistogram1D* m_h_phidiff_posEC ;
  IHistogram1D* m_h_etadiff_posEC ;
  
  IHistogram1D* m_h_etadiff ;
  
 
  int m_histStart;
  int m_nHist;

  int m_particleType ;

  std::string m_histTitle ;

};


} // end of namespace bracket

#endif

---