HiggsAnalysis/HiggsAnalysis-00-01-00/HiggsAnalysis/JetTruthTrackSelectorTool.h

#ifndef JETREC_JetTruthTrackSelectorTool_H
#define JETREC_JetTruthTrackSelectorTool_H

#include "Particle/TrackParticleContainer.h"
#include "JetRec/JetAlgToolBase.h"

#include <string>

class StoreGateSvc;

class GenAccessIO;

class JetTruthTrackSelectorTool : public JetAlgToolBase
{
 public:

  JetTruthTrackSelectorTool( const std::string& name, const std::string& type,
                     const IInterface* parent );

  virtual ~JetTruthTrackSelectorTool();

  virtual StatusCode initialize();
  virtual StatusCode execute(JetAlgToolBase::jetcollection_t* theJets); 
  static JetTruthTrackSelectorTool * getPointer(void);

 protected:

  // Properties //

  double m_ptMin, m_etaMax;
  bool   m_includeMuons;
  std::string m_outputCollectionName;
 
  // Store //

  StoreGateSvc* m_storeGate;
  GenAccessIO* m_truthAccess;

 private:
    
    bool m_ignore_multivertices; // option how to treat multivertex events
    bool m_remove_isol_leptons;  //flag to remove isolated e,mu from tracks list

    double m_electron_cone;      //cone size for electron isolation
    double m_electron_et_cut;    //et cut for electron isolation
    double m_muon_cone;          //cone size for muon isolation
    double m_muon_et_cut;        //et cut for muon isolation
    double m_inner_cone;         //size of the inner cone 
                                 //  -> et_in_cone: inner_cone < dR < m_lepton_cone
    

    double m_Et_min_vert; // minimum track Et in vertex selection
    double m_vertex_cut; // vertex cut
    
    

// private methods //
        Hep3Vector highest_Et_vertex(
                        std::vector<const HepMC::GenParticle*> theParticles);
                                                // get the primary vertex
                                                // with the highest transverse
                                                // momentum;
                                                // theParticles: truth particle
                                                // collection

    double get_energy_in_cone(const double & cone,
                              const HepLorentzVector & p, 
                              const Hep3Vector & vertex,
                              std::vector<const HepMC::GenParticle*> theParticles);
                                        // get the energy of particle emerging
                                        // from the vertex "vertex" in a cone
                                        // of "cone" in delta R of p

};
#endif

---