---

TrackMaker.h

Go to the documentation of this file.

// ================================================
// TrackMaker class description
// ================================================
//
// THIS TEXT TO BE REPLACED BY ATLAS STANDARD FORMAT
//
//
// This version....
//
//
// Namespace Atlfast::
//
// class: TrackMaker
//
// Description: 
//
//  Default Algorithm which knows how to simulate Tracks from MC truth 
//
//
//
// ................................................................
//

#ifndef __ATLFAST_Trackmaker__
#define __ATLFAST_Trackmaker__

// 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 "CLHEP/Vector/ThreeVector.h"
#include "HepMC/GenEvent.h"
#include "HepMC/GenParticle.h"
#include "HepMC/ParticleDataTable.h"
#include "HepMC/IO_PDG_ParticleDataTable.h"
#include "AtlfastCode/HepMC_helper.h"


// Atlfast
#include "AtlfastCode/IKinematic.h"
#include "AtlfastCode/TrackTrajectory.h"
#include "AtlfastCode/Track.h"
#include "AtlfastCode/TrackCollection.h"
#include "AtlfastCode/TrackSmearer.h"
#include "AtlfastCode/TesIO.h"
#include "AtlfastCode/CommonData.h"
#include "AtlfastCode/HeaderPrinter.h"
#include "AtlfastCode/MCparticleCollection.h"

//*************************************
//Default parameters private to this algorithm
//*************************************

#define DEFAULT_mcPtMin     0.5   // Minimum pT  (applied to MC truth)
#define DEFAULT_mcEtaMax    2.5   // Maximum eta (applied to MC truth) 
#define DEFAULT_vtMax       3.0   // Maximum Vertex(trans)
#define DEFAULT_vlMax       40.0  // Maximum Vertex(long) 
#define DEFAULT_doSmearing  true  // Should the particle four vectors be smeared or not?
#define DEFAULT_bField      2.0   // magnetic field in Tesla
 
//htp random number generator should eventually be removed
#define DEFAULT_muConfig    "000"


//*****************************************
// Default paths in the Transient event store to get/put entities
//
//[Note:  This will change ]
//******************************************
#define DEFAULT_commonDataLocation "/Event/AtlfastCommonData" 
#define DEFAULT_MC_eventLocation    "/Event/McEventCollection" 
#define DEFAULT_outputLocation      "/Event/AtlfastTracks" 



//***************************************************************
// TrackMaker class declaration
//
//****************************************************************

namespace Atlfast {

class TrackMaker : public Algorithm
{

 public:

  //-------------------------
  // Constructors/Destructors
  //-------------------------

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


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

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


 private:

  // -------------------------------
  // Parameters of this algorithm
  // (descriptions above)
  //--------------------------------
  bool          m_notInit;
  double        m_vlMax;
  double        m_vtMax;
  double        m_mcPtMin;     
  double        m_mcEtaMax;       
  bool          m_doSmearing;
  double        m_bField ;
  string        m_muConfig;



  //--------------------------------------------
  // Selector for HepMC generator particles
  //---------------------------------------------


  TesIO* m_tesIO;
  HepMC_helper::IMCselector*  m_ncutter ;

  //--------------------------------------------
  // Smearer helper class which knows all about how to do smearing
  // for Tracks
  // There may need to be more of these added in order
  // to treat different species differently.
  // Alternately the single smearer object could be made to
  // dispatch on particle type. Not thought through yet.
  //---------------------------------------------


  TrackSmearer* m_smearer;
  HepMC::ParticleDataTable m_particleTable;


  //--------------------------------------------------
  // Paths in the Transient Event store to get/put things
  // These are currently set in member variables so that
  // they can be overwritten by the job options service
  //--------------------------------------------------

  std::string m_commonDataLocation ;
  std::string m_MC_eventLocation ;
  std::string m_outputLocation ;


  //-------------------------------
  // Private methods
  //-------------------------------

  // Create a Track 
  virtual Track* create( const HepMC::GenParticle* );


  // Test if Track passes kinematic requirements
  //  virtual bool isAcceptable( Track* );
  virtual bool isAcceptable(Track*);
  TrackTrajectory extractTrajectory( const HepMC::GenParticle* ) ;

};


}  //end of namespace bracket

#endif

---