---

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


// *********************************************************************
// Default Parameter of this algorithm
// (see main class for description)
#define DEFAULT_particleType 11   
#define DEFAULT_mcPtMin      0.0  
#define DEFAULT_mcEtaMax   100.0  
#define DEFAULT_PtMin        5.0  
#define DEFAULT_EtaMax       2.5  
#define DEFAULT_doSmearing  true  
#define DEFAULT_MC_eventLocation   "/Event/McEventCollection" 
#define DEFAULT_commonDataLocation "/Event/AtlfastCommonData" 
#define DEFAULT_outputLocation     "/Event/AtlfastReconstructedParticle" 
#define DEFAULT_muSmearKey    1   
//********************************************************************


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

namespace HepMC_helper{
  class IMCselector;
}

namespace Atlfast {

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
  //--------------------------------
  bool m_notInit;
  // 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_MC_eventLocation ;
  std::string m_commonDataLocation ;
  std::string m_outputLocation ;

  // Things which should not be in this class eventually, but have to be now
  int           m_muSmearKey;

  //-----------------------------------------------------------------
  // 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;

  ReconstructedParticle * lnkReconstructedParticle;
  //--------------------------
  // Private  methods 
  //--------------------------
  // virtual StatusCode getMC ( MsgStream&, t_MC_particleCollection& );

  TesIO* m_tesIO;

  void getSmearer(int lumi, int seed, MsgStream& log );
  virtual ReconstructedParticle* create( MsgStream&, const HepMC::GenParticle* );
  virtual bool isAcceptable( MsgStream&, const ReconstructedParticle* );

};


}  //end of namespace bracket

#endif

---