---

Isolator.h

Go to the documentation of this file.

// ================================================
// Isolator class description
// ================================================
//
// THIS TEXT TO BE REPLACED BY ATLAS STANDARD FORMAT
//
//
// This version....
//
//
// Namespace Atlfast::
//
// class: Isolator
//
// Description: 
//
//  Algorithm which knows how to test particles for isolation
//
//
// ................................................................
//


#ifndef __Atlfast_isolator__
#define __Atlfast_isolator__

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

// Gaudi 
#include "GaudiKernel/ISvcLocator.h"
#include "GaudiKernel/IAlgorithm.h"
#include "GaudiKernel/Algorithm.h"
#include "GaudiKernel/MsgStream.h"
#include "GaudiKernel/DataObject.h"
#include "GaudiKernel/ObjectVector.h"
// Other
#include "HepMC/GenParticle.h"

// Atlfast
#include "AtlfastCode/ReconstructedParticle.h"
#include "AtlfastCode/ReconstructedParticleCollection.h"
#include "AtlfastCode/Cell.h"
#include "AtlfastCode/CellCollection.h"
#include "AtlfastCode/Cluster.h"
#include "AtlfastCode/ClusterCollection.h"

#include "AtlfastCode/IKinematic.h"
#include "AtlfastCode/KinematicHelper.h"
#include "AtlfastCode/TesIO.h"
#include "AtlfastCode/CommonData.h"


//*********************************************
//Default parameters private to this algorithm
//Description in class declaration
//*********************************************
#define DEFAULT_rClusterMatch      0.150  
#define DEFAULT_rClusterIsolation  0.400  
#define DEFAULT_eClusterIsolation  0.0   
#define DEFAULT_rCellIsolation     0.200   
#define DEFAULT_eCellIsolation    10.0  
#define DEFAULT_inputLocation               "Unknown" 
#define DEFAULT_isolatedOutputLocation      "Unknown" 
#define DEFAULT_nonIsolatedOutputLocation   "Unknown" 
#define DEFAULT_cellLocation     "/Event/AtlfastCell" 
#define DEFAULT_clusterLocation  "/Event/AtlfastCluster" 
#define DEFAULT_commonDataLocation  "/Event/AtlfastCommonData" 




//***************************************************************
//                 Isolator class declaration
//
//****************************************************************


namespace Atlfast {
class Isolator : public Algorithm
{

 public:

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

  Isolator( const std::string& name, ISvcLocator* pSvcLocator ) ; 

  virtual ~Isolator();


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

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



 private:

  //---------------------------------------------------
  // Parameters of this algorithm 
  //---------------------------------------------------

  // R-cone for associating a Cluster to a particle
  double        m_rClusterMatch;   

  // R-cone in which to sum unassociated Cluster eT
  double        m_rClusterIsolation;     

  // eT threshold to apply to summed Cluster eT
  double        m_eClusterIsolation;     

  // R-cone in which to sum Cell eT
  double        m_rCellIsolation;     

  // eT threshold to apply to excess summed Cluster eT
  double        m_eCellIsolation;     

  // Location in TES to obtain input particles to be isolated
  std::string   m_inputLocation ;

  // Locations in TES to write out output lists.
  std::string   m_isolatedOutputLocation ;
  std::string   m_nonIsolatedOutputLocation ;

  // Locarions in TES to get Cells and Clusters from
  std::string   m_cellLocation ;
  std::string   m_clusterLocation ;

  std::string m_commonDataLocation;

  TesIO* m_tesIO;
  bool   m_notInit;
  //-------------------------------
  // Objects used by this class to do its job
  //-------------------------------

  // Encapsulates common operations on IKinematic types
  KinematicHelper m_kinehelp ; 
  HepMC_helper::IMCselector* m_visibleToCal;

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

  // This is the heart of this algorithm: the bit that does the isolation

  bool isIsolated( 
     MsgStream&, 
     ReconstructedParticleCollection::iterator,
     std::vector<Cell*>& , 
     std::vector<Cluster*>&  

     );
 
};

} // end of namespace bracket


#endif

---