---

JetMaker.h

Go to the documentation of this file.

// ================================================
// JetMaker class description
// ================================================
//
// THIS TEXT TO BE REPLACED BY ATLAS STANDARD FORMAT
//
//
// This version....
//
//
// Namespace Atlfast::
//
// class: JetMaker
//
// Description: 
//
//  Algorithm which knows how to make Jets
//
//  The resulting Jets are added to the Event
//
// ................................................................
//

#ifndef __ATLFAST_JetMaker__
#define __ATLFAST_JetMaker__

// 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"
#include "AtlfastCode/HepMC_helper.h"


// Atlfast
#include "AtlfastCode/KinematicHelper.h"
#include "AtlfastCode/Cluster.h"
#include "AtlfastCode/ISmearer.h"
#include "AtlfastCode/Jet.h"
#include "AtlfastCode/JetCollection.h"
#include "AtlfastCode/ReconstructedParticle.h"
#include "AtlfastCode/ReconstructedParticleCollection.h"
#include "AtlfastCode/MCparticleCollection.h"
#include "AtlfastCode/JetSmearer.h"
#include "AtlfastCode/FunctionObjects.h"
//#include "AtlfastCode/Muon.h"
#include "AtlfastCode/TesIO.h"

//*************************************
//Default parameters private to this algorithm
//*************************************
#define DEFAULT_minPT       10    // Minimum jet transverse momentum
#define DEFAULT_maxEta      5     // eta range to consider 
#define DEFAULT_doSmearing  true  // Smear the jet four vector?
 
#define DEFAULT_rconeb      0.400
#define DEFAULT_rconef      0.400
//*****************************************
//cuts on MC particles used to tag jets
#define DEFAULT_bPtMin          5.0
#define DEFAULT_cPtMin          5.0
#define DEFAULT_tauPtMin       10.0
#define DEFAULT_bMaxDeltaR      0.2
#define DEFAULT_cMaxDeltaR      0.2
#define DEFAULT_tauMaxDeltaR    0.3
#define DEFAULT_etaTagMax       2.5
#define DEFAULT_tauJetPtRatio   0.9
//*****************************************
// Default paths in the Transient event store to get/put entities
//
//******************************************

#define DEFAULT_inputLocation           "/Event/AtlfastClusters" 
#define DEFAULT_outputLocation          "/Event/AtlfastJets" 
#define DEFAULT_muonLocation            "/Event/AtlfastNonIsolatedMuons"
#define DEFAULT_commonDataLocation      "/Event/AtlfastCommonData"
#define DEFAULT_mcTruthLocation         "/Event/McEventCollection"
#define DEFAULT_unusedCellLocation      "/Event/AtlfastUnusedCells"
#define DEFAULT_missingMomentumLocation "/Event/AtlfastMissingMomentum"
//***************************************************************
// JetMaker class declaration
//
// A "Gaudi algorithm" is something which gets scheduled and controlled
// by the framework. In its simplest definition it is something which
//  - gets clled for each event 
//  - can get anything it wants out of the Transient Event Store"
//  - can put anything it makes into the store.

//
//****************************************************************

namespace Atlfast {
  class JetMaker : public Algorithm
{

 private:
  //------------------------------------
  // Types used internally by this class
  // for local collections 
  //------------------------------------

  typedef std::vector<Cluster*>                         localClusterCollection;
  typedef std::vector<Cluster*>::iterator               localClusterIterator;
  typedef std::vector<ReconstructedParticle*>           MuonCollection;
  typedef std::vector<ReconstructedParticle*>::iterator MuonIterator;

  //--------------------------------
  // Parameters of this algorithm
  // (descriptions above)
  //--------------------------------

  bool          m_notInit;
  double        m_minPT;     
  double        m_maxEta;       
  bool          m_doSmearing;

  double        m_barrelForwardEta;
  double        m_rconeb;
  double        m_rconef;

  double m_bPtMin;
  double m_bMaxDeltaR;
  double m_cPtMin;
  double m_cMaxDeltaR;
  double m_tauPtMin;
  double m_tauMaxDeltaR;
  double m_etaTagMax;
  double m_tauJetPtRatio;

  //--------------------------------------------
  // Smearer helper class
  // which knows all about how to do smearing
  // for my particles
  //---------------------------------------------


  ISmearer* m_smearer;
  ISmearer* m_cellSmearer;


  TesIO* m_tesIO;

  //--------------------------------------------------
  // 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_inputLocation ;
  std::string m_outputLocation ;
  std::string m_muonLocation;
  std::string m_mcTruthLocation;
  std::string m_unusedCellLocation;
  std::string m_commonDataLocation;
  std::string m_missingMomentumLocation;



  //-------------------------------
  // helper
  //-------------------------------

  KinematicHelper m_kinehelp ;

  int getJetTags(MCparticleCollection&  mcParticles,
                  const PartitionCondition::BelowThresholdDeltaR& tagSelector
                  )const;
  int getTauTags(MCparticleCollection&  mcParticles,
                 const PartitionCondition::BelowThresholdDeltaR& tagSelector,
                 double jetPt )const;
  void tag(Jet* jet, 
           MCparticleCollection&  mcParticles_b, 
           MCparticleCollection&  mcParticles_c, 
           MCparticleCollection&  mcParticles_tau) const;
 
  HepLorentzVector addUnusedCells(MsgStream&) const;

  //-------------------------------
  // Protected methods
  //-------------------------------
 public:

  //-------------------------
  // Constructors/Destructors
  //
  // Gaudi requires that the constructor takes certain arguments
  // (and passes them directly to the constructor of the base class)
  //-------------------------

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


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

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

  //------------------------------------------------------
  // Methods of IReconstructedParticleMaker which we are
  // implementing
  //------------------------------------------------------

   Jet* create( Cluster* );
   bool isAcceptable( Jet* );
   void addMuons(Jet*, MuonCollection& );
   MsgStream& dumpParameters(MsgStream& ) const; 

};


}  //end of namespace bracket

#endif

---