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_H
#define ATLFAST_RECONSTRUCTEDPARTICLEHISTOGRAMMAKER_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"
 
namespace HepMC_helper{
  class IMCselector;
}
namespace Atlfast {
  using std::string;
  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;
    //holds the Stotegate MC location
    std::string m_mcLocation;

    //-------------------------------
    // 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 ;
    
    
  //-----------------------------------------------------------------
  // Selector to select only those truth particles which are required
  //-----------------------------------------------------------------
    HepMC_helper::IMCselector*  m_ncutter ;
  };
  
} // end of namespace bracket

#endif

---