Analysis/Triggers/ConeMCFM.H

0001 /* Feb 6th 2009
0002  *
0003  * Close copy of MCFM-5.3/src/User/genclust_cone.f cone algorithm implementation.
0004  * Implementation for up to 4 jets.
0005  * Uses sherpa infrastructure.
0006  */
0007
0008 #include "ATOOLS/Phys/Particle_List.H"
0009 #include "ATOOLS/Math/Poincare.H"
0010
0011 namespace ATOOLS {
0012   class ConeMCFM {
0013     double m_Rmin, m_Rsep;
0014
0015     void MergeSplit(std::vector<Vec4D> & qjet ,std::vector<Vec4D> & protoq, std::vector< std::vector<int> > & protoc, std::vector<Vec4D> & qfinal);
0016
0017     double deltarq(Vec4D p1,Vec4D p2);
0018     double etarap(Vec4D p1);
0019     double getet(Vec4D p1);
0020
0021   public:
0022     std::vector<Vec4D> m_pjets;                         /* jet momenta */
0023     void ConstructJets(std::vector<Vec4D> & p);                 /* constructs jet momenta stored in m_pjets.  */
0024     ConeMCFM(double Rmin, double Rsep);                         /* constructs on jet momenta stored.         */
0025     ~ConeMCFM() {}
0026
0027   };
0028 }