---

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_H
#define ATLFAST_ISOLATOR_H

// 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 "AtlfastEvent/ReconstructedParticle.h"
//#include "AtlfastEvent/ReconstructedParticleCollection.h"
#include "AtlfastEvent/CollectionDefs.h"
#include "AtlfastEvent/Cell.h"
#include "AtlfastEvent/Cluster.h"

#ifndef ATLFAST_COLLECTIONDEFS_H
#include "AtlfastEvent/CollectionDefs.h"
#endif

#include "AtlfastEvent/IKinematic.h"
#include "AtlfastUtils/KinematicHelper.h"
#include "AtlfastUtils/TesIO.h"


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


namespace Atlfast {
  using std::string;
  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 ;
    
    
    
    TesIO* m_tesIO;
    //-------------------------------
    // 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

---