DefaultReconstructedParticleMaker.h

Go to the documentation of this file.

// ================================================
// DefaultReconstructedParticleMaker class description
// ================================================
//
// THIS TEXT TO BE REPLACED BY ATLAS STANDARD FORMAT
//
//
// This version....
//
//
// Namespace Atlfast::
//
// class: DefaultReconstructedParticleMaker
//
// Description: 
//
//  Default Algorithm to make smeared 'Atlfast::ReconstructedParticles' from 
//  MC truth 'HepMC::GenParticles'
//
//  HepMC::GenParticles are accessed in the TES. Only those satisfying
//  specified pre-selection requirements are used by this algorithm. 
//  These requirements are:
//     - particle type (given by pdg_id) 
//     - minimum pT
//     - maximum eta
//  These are set via job options.
//
//  The different treatement of different species of particle is embodied
//  completely in a smearer object.  
//  Arbitrarily complex smearers can be supplied if desired.
//
//  For each HepMC::GenParticle, a candidate ReconstructedParticle is
//  created. This is only kept if
//  it satisfies the following requirements:
//      - minimum pT
//      - maximum eta
//   
//  The surviving ReconstructedParticles are stored in the TES
//
// ................................................................


#ifndef ATLFAST_DEFAULTRECONTRUCTEDPARTICLEMAKER_H
#define ATLFAST_DEFAULTRECONTRUCTEDPARTICLEMAKER_H

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

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

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

// Atlfast
#include "AtlfastEvent/ReconstructedParticle.h"
#include "AtlfastEvent/CollectionDefs.h"
#include "AtlfastAlgs/ISmearer.h"
#include "AtlfastAlgs/IAcceptor.h"
#include "AtlfastUtils/TesIO.h"


namespace HepMC_helper{
  class IMCselector;
}

namespace Atlfast {
  using std::string;

class DefaultReconstructedParticleMaker : 
    public virtual Algorithm 
{

 public:

  //-------------------------
  // Constructors/Destructors
  //-------------------------
  DefaultReconstructedParticleMaker(const std::string& name, 
                                    ISvcLocator* pSvcLocator);
  virtual ~DefaultReconstructedParticleMaker();
  
  
  //------------------------------------------------------
  // Mandatory Athena methods
  //------------------------------------------------------
  virtual StatusCode initialize() ;
  virtual StatusCode execute() ;
  virtual StatusCode finalize() ;
  

 private:

  //------------------------------------
  // Typedefs used internally by this class 
  //------------------------------------

  typedef std::vector<HepMC::GenParticle*>           t_MC_particleCollection ;
  typedef std::vector<HepMC::GenParticle*>::iterator t_MC_particleIterator ;


  //--------------------------------
  // Parameters of this algorithm
  //--------------------------------
  // General parameters
  // Pre-selection criteria to apply to the HepMC::GenParticles  
  int           m_particleType;
  double        m_mcPtMin ;     
  double        m_mcEtaMax ;    

  // Post creation criteria to apply to created ReconstructedParticles 
  double        m_PtMin; 
  double        m_EtaMax;      

  bool          m_doSmearing;    
  
  // Paths in the Transient event store to get/put entities
  std::string m_outputLocation;
  std::string m_muonResFile;

  int           m_muSmearKey;

  //cct: array to hold smearing parameters loaded from JobOptions
  std::vector<double>      m_smearParamArray;
  //cct: integer to hold smearSchema loaded from JobOptions
  int           m_smearParamSchema;

  // Apply efficiencies at particle creation stage?
  bool    m_applyEfficiencies;
  
  //-----------------------------------------------------------------
  // Selector to select only those truth particles which are required
  //-----------------------------------------------------------------
  HepMC_helper::IMCselector*  m_ncutter ;


  //--------------------------------------------
  // Smearer helper class which specialises this algorithm for a
  // specific particle species
  //---------------------------------------------
  ISmearer* m_smearer;

  // SES resolution factor ********************
  //CorrectionFactor* m_CorrFactor;

  // Class to accept or reject particles
  IAcceptor *m_acceptor;

  ReconstructedParticle * m_lnkReconstructedParticle;
  //--------------------------
  // Private  methods 
  //--------------------------

  TesIO* m_tesIO;

  MsgStream m_log;
  
  void getAcceptor();

  virtual ReconstructedParticle* create( const HepMC::GenParticle* );

  virtual bool isAcceptable( const ReconstructedParticle* );
  
  std::string m_mcLocation;

};


}  //end of namespace bracket

#endif

---