Analysis/Triggers/MySISCone.H

0001 #ifndef ANALYSIS_Triggers_MySISCone_H
0002 #define ANALYSIS_Triggers_MySISCone_H
0003
0004 #include "AddOns/Analysis/Triggers/Kt_Algorithm.H"
0005 #include "ATOOLS/Phys/Particle_List.H"
0006 #include "ATOOLS/Math/Vector.H"
0007 #include "ATOOLS/Math/Histogram.H"
0008 #include <list>
0009 #include <vector>
0010
0011 namespace siscone {class Csiscone;}
0012 namespace ANALYSIS {
0013   class Primitive_Analysis;
0014
0015   class SISCone : public Jet_Algorithm_Base {
0016
0017
0018     //   temporary variable to store the jet momenta
0019     ATOOLS::Particle_List       *p_jets;
0020     const ATOOLS::Particle_List *p_orig;
0021     std::vector<double> * p_kts;
0022
0023     double m_f;
0024
0025     siscone::Csiscone* p_ConeCluster;
0026
0027     ATOOLS::Flavour GetBFlavour(const ATOOLS::Vec4D & mom);
0028     void AddToKtlist(double );
0029     void AddToJetlist(const ATOOLS::Vec4D &);
0030
0031   public:
0032     SISCone(ATOOLS::Particle_Qualifier_Base * const qualifier,double f=.5);
0033     ~SISCone();
0034
0035     bool ConstructJets(const ATOOLS::Particle_List * pl = 0,ATOOLS::Particle_List * jets = 0,
0036                std::vector<double> * kt2 = 0,double R=.7);
0037     int    NumberOfJets();
0038
0039   };
0040
0041
0042
0043   inline int    SISCone::NumberOfJets() { return p_jets->size(); }
0044
0045   //
0046 }
0047
0048 #endif