RhoEventShapes Class Reference

#include <RhoEventShapes.h>

Public Member Functions
	RhoEventShapes (RhoCandList &l, TLorentzVector cms)
virtual	~RhoEventShapes ()
int	NParticles () const
int	NCharged () const
int	NNeutral () const
double	PmaxLab () const
double	PmaxCms () const
double	Ptmax () const
double	PtSumLab () const
double	NeutEtSumLab () const
double	NeutESumLab () const
double	ChrgPtSumLab () const
double	ChrgPSumLab () const
double	PtSumCms () const
double	NeutEtSumCms () const
double	NeutESumCms () const
double	ChrgPtSumCms () const
double	ChrgPSumCms () const
int	MultPminLab (double pmin)
int	MultPmaxLab (double pmax)
int	MultPminCms (double pmin)
int	MultPmaxCms (double pmax)
double	Sphericity ()
double	Aplanarity ()
double	Planarity ()
double	FoxWolfMomH (int order)
double	FoxWolfMomR (int order)
double	Thrust ()
TVector3	ThrustVector ()
	ClassDef (RhoEventShapes, 1)

Private Member Functions
void	ComputeSphericity ()
double	Eps (const TVector3 v1, const TVector3 v2)
double	Legendre (int l, double x)

Private Attributes
std::vector< TLorentzVector >	fLabList
	! List of 4-vectors in lab frame More...
std::vector< TLorentzVector >	fCmsList
	! List of 4-vectors in cms frame More...
std::vector< int >	fCharge
	! List of charges of particles More...
int	fnChrg
int	fnNeut
int	fN
double	fpmaxlab
double	fpmaxcms
double	fptmax
double	fptsumlab
double	fneutetsumlab
double	fneutesumlab
double	fchrgptsumlab
double	fchrgpsumlab
double	fptsumcms
double	fneutetsumcms
double	fneutesumcms
double	fchrgptsumcms
double	fchrgpsumcms
double	fsph
double	fapl
double	fpla
double	fFWmom [FWMAX]
bool	fFWready
double	fthr
TVector3	fThrVect
TVector3	fBoost
bool	fNeutralOnly
bool	fChargedOnly

Detailed Description

Definition at line 24 of file RhoEventShapes.h.

Constructor & Destructor Documentation

RhoEventShapes::RhoEventShapes	(	RhoCandList &	l,
		TLorentzVector	cms
	)

Definition at line 7 of file RhoEventShapes.cxx.

References RhoCandidate::Charge(), fBoost, fCharge, fchrgpsumcms, fchrgpsumlab, fchrgptsumcms, fchrgptsumlab, fCmsList, fFWmom, fLabList, fN, fnChrg, fneutesumcms, fneutesumlab, fneutetsumcms, fneutetsumlab, fnNeut, fpmaxcms, fpmaxlab, fptmax, fptsumcms, fptsumlab, fThrVect, FWMAX, RhoCandList::GetLength(), i, and RhoCandidate::P4().

                                                                 : 
  fnChrg(0), fnNeut(0), fN(0), fpmaxlab(0.), fpmaxcms(0.), fptmax(0.),
  fptsumlab(0.),fneutetsumlab(0.),fneutesumlab(0.),fchrgptsumlab(0.),fchrgpsumlab(0.),
  fptsumcms(0.),fneutetsumcms(0.),fneutesumcms(0.),fchrgptsumcms(0.),fchrgpsumcms(0.),
  fsph(-1.), fapl(-1.), fpla(-1.), fFWready(false), fthr(-1.)
{
  int i;

  // initialize more complex things
  fThrVect.SetXYZ(0.,0.,0.);
  fBoost=cms.BoostVector();  
  for (i=0;i<=FWMAX;++i) fFWmom[i]=0.;
  
  
  double pmax=0., ptmax=0., pmaxcms=0.;
  
  fLabList.clear();
  fCmsList.clear();
  
  for (i=0;i<l.GetLength();++i)
  {
    RhoCandidate* cand = l[i];
    TLorentzVector lv(cand->P4());
    int chrg(cand->Charge());

    fN++;
    
    // cache multiplicities
    if (chrg==0) fnNeut++;
    else fnChrg++;
    
    // cache unboosted 4-vectors
    fLabList.push_back(lv);
    // cache charges
    fCharge.push_back(chrg);
    
    // sum momentum variables (lab)
    fptsumlab += lv.Pt();
    if (chrg==0)
    {
      fneutetsumlab += lv.Pt();  
      fneutesumlab  += lv.E();
    }
    else
    {
      fchrgptsumlab += lv.Pt();
      fchrgpsumlab  += lv.P();
    }
    
    // cache maximum momenta in lab
    if (lv.P()>pmax) pmax=lv.P();
    if (lv.Pt()>ptmax) ptmax=lv.Pt();
    
    // cache boosted vectors
    lv.Boost(-fBoost);
    fCmsList.push_back(lv);

    // sum momentum variables (cms)
    fptsumcms += lv.Pt();
    if (chrg==0)
    {
      fneutetsumcms += lv.Pt();  
      fneutesumcms  += lv.E();
    }
    else
    {
      fchrgptsumcms += lv.Pt();
      fchrgpsumcms  += lv.P();
    }
    
    // cache maximum momenta in cms
    if (lv.P()>pmaxcms) pmaxcms=lv.P();
  }
  
  fpmaxlab = pmax;
  fptmax   = ptmax;
  fpmaxcms = pmaxcms;
}

virtual RhoEventShapes::~RhoEventShapes ( )

inlinevirtual

Definition at line 27 of file RhoEventShapes.h.

{};

Member Function Documentation

double RhoEventShapes::Aplanarity

(

)

Definition at line 139 of file RhoEventShapes.cxx.

References ComputeSphericity(), and fapl.

{
  if (fapl<0.) ComputeSphericity();
  
  return fapl;
}

double RhoEventShapes::ChrgPSumCms ( ) const

inline

Definition at line 51 of file RhoEventShapes.h.

References fchrgpsumcms.

{return fchrgpsumcms;}          // sum of momenta of charged (lab)

double RhoEventShapes::ChrgPSumLab ( ) const

inline

Definition at line 44 of file RhoEventShapes.h.

References fchrgpsumlab.

{return fchrgpsumlab;}          // sum of momenta of charged (lab)

double RhoEventShapes::ChrgPtSumCms ( ) const

inline

Definition at line 50 of file RhoEventShapes.h.

References fchrgptsumcms.

{return fchrgptsumcms;}         // sum of pt of charged (lab)

double RhoEventShapes::ChrgPtSumLab ( ) const

inline

Definition at line 43 of file RhoEventShapes.h.

References fchrgptsumlab.

{return fchrgptsumlab;}         // sum of pt of charged (lab)

RhoEventShapes::ClassDef	(	RhoEventShapes	,
		1
	)

void RhoEventShapes::ComputeSphericity ( )

private

Definition at line 88 of file RhoEventShapes.cxx.

References fapl, fCmsList, fN, fpla, fsph, i, and v.

Referenced by Aplanarity(), Planarity(), and Sphericity().

{
  if( fN==0 ) return;
  
  double stot=0, sxx=0, sxy=0, sxz=0, syy=0, syz=0, szz=0;  
  int i;
  
  for (i=0;i<fN;++i)
  {
    TVector3 v(fCmsList[i].Vect());
    
    sxx += v.X()*v.X(); sxy += v.X()*v.Y(); sxz += v.X()*v.Z();
    syy += v.Y()*v.Y(); syz += v.Y()*v.Z();
    szz += v.Z()*v.Z();
    
    stot += v.Mag2();
  }
  
  TMatrixD sm(3,3);
  sm(0,0) = sxx/stot; sm(0,1) = sxy/stot; sm(0,2) = sxz/stot; 
  sm(1,0) = sxy/stot; sm(1,1) = syy/stot; sm(1,2) = syz/stot; 
  sm(2,0) = sxz/stot; sm(2,1) = syz/stot; sm(2,2) = szz/stot; 
  
  TMatrixDEigen ei(sm);
  TMatrixD eiv=ei.GetEigenValues();
  
  fsph = 1.5 * (eiv(1,1) + eiv(2,2));
  fapl = 1.5 * eiv(2,2);
  fpla = eiv(1,1) - eiv(2,2);  
}

double RhoEventShapes::Eps ( const TVector3  v1, const TVector3  v2  )

inlineprivate

Definition at line 73 of file RhoEventShapes.h.

Referenced by Thrust().

{return (v1*v2)>0. ? 1. : -1.;}  // aux for Thrust

double RhoEventShapes::FoxWolfMomH

(

int

order

)

Definition at line 235 of file RhoEventShapes.cxx.

References cos(), fCmsList, fFWmom, fFWready, fN, FWMAX, i, Legendre(), p1, p2, and s.

Referenced by FoxWolfMomR().

{
  if (order>FWMAX) return -1.;
 
  // already computed FW moments
  if (fFWready) return fFWmom[order];
  
  if( fN==0 ) return -1.;

  double s = 0.;
  int i,j,l;

  for (i=0; i<fN-1; ++i)
  {
     // this candidate's 3-momentum
      TVector3 p1(fCmsList[i].Vect());
      double pmag1 = p1.Mag();

      // loop over other candidates, starting at the next one in the list
      
      for (j=i+1; j<fN; ++j)
      {
        // this candidate's 3-momentum
        TVector3 p2(fCmsList[j].Vect());
        double pmag2 = p2.Mag();

        // the cosine of the angle between the two candidates
        double cosPhi =  cos ( p1.Angle(p2) );

        // the contribution of this pair of track 
        // (note the factor 2 : the pair enters the sum twice)
        for( l=0; l<=FWMAX; l++ )
          fFWmom[l] += 2 * pmag1 * pmag2 * Legendre( l, cosPhi );
          
      }
      // contribution of this track
      for( l=0; l<=FWMAX; l++ )
        fFWmom[l] += pmag1 * pmag1 * Legendre( l, 1. );
      
      // total energy
      s += fCmsList[i].Energy();
   }
  
  // well ...
  if( s<=0. ) return -1.;
  double s2=s*s;
  
  // normalize Fox Wolfram Moments
  for( i=0; i<=FWMAX; i++) fFWmom[i]/=s2 ;  
  
  fFWready = true;
  
  return fFWmom[order];
}

double RhoEventShapes::FoxWolfMomR

(

int

order

)

Definition at line 293 of file RhoEventShapes.cxx.

References FoxWolfMomH().

{
  return FoxWolfMomH(order)/FoxWolfMomH(0);
}

double RhoEventShapes::Legendre ( int  l, double  x  )

private

Definition at line 212 of file RhoEventShapes.cxx.

References fabs(), and x.

Referenced by FoxWolfMomH().

{
    if (fabs(x)>1.) return -999.;

    if (l==0) return 1.;
    
    double pmm = 1.;
    double pmmp1 = x;

    if (l>1)
    {
        for(int ll=2; ll<=l; ll++)
        {
          double pll = (x * (2 * ll - 1) * pmmp1 - (ll - 1) * pmm) / (ll);
          pmm = pmmp1;
          pmmp1 = pll;
        }
    }
   return pmmp1;
}

int RhoEventShapes::MultPmaxCms

(

double

pmax

)

Definition at line 333 of file RhoEventShapes.cxx.

References cnt, fCmsList, fN, i, and P.

{
  int cnt=0;
  for (int i=0;i<fN;++i) 
    if (fCmsList[i].P()<pmax) cnt++;
 
  return cnt;
}

int RhoEventShapes::MultPmaxLab

(

double

pmax

)

Definition at line 311 of file RhoEventShapes.cxx.

References cnt, fLabList, fN, i, and P.

{
  int cnt=0;
  for (int i=0;i<fN;++i) 
    if (fLabList[i].P()<pmax) cnt++;
 
  return cnt;
}

int RhoEventShapes::MultPminCms

(

double

pmin

)

Definition at line 322 of file RhoEventShapes.cxx.

References cnt, fCmsList, fN, i, and P.

{
  int cnt=0;
  for (int i=0;i<fN;++i) 
    if (fCmsList[i].P()>pmin) cnt++;
 
  return cnt;
}

int RhoEventShapes::MultPminLab

(

double

pmin

)

Definition at line 300 of file RhoEventShapes.cxx.

References cnt, fLabList, fN, i, and P.

{
  int cnt=0;
  for (int i=0;i<fN;++i) 
    if (fLabList[i].P()>pmin) cnt++;
 
  return cnt;
}

int RhoEventShapes::NCharged ( ) const

inline

Definition at line 31 of file RhoEventShapes.h.

References fnChrg.

{return fnChrg;}                // number of charged candidates

double RhoEventShapes::NeutESumCms ( ) const

inline

Definition at line 49 of file RhoEventShapes.h.

References fneutesumcms.

{return fneutesumcms;}          // sum of energys of neutrals (lab)

double RhoEventShapes::NeutESumLab ( ) const

inline

Definition at line 42 of file RhoEventShapes.h.

References fneutesumlab.

{return fneutesumlab;}          // sum of energys of neutrals (lab)

double RhoEventShapes::NeutEtSumCms ( ) const

inline

Definition at line 48 of file RhoEventShapes.h.

References fneutetsumcms.

{return fneutetsumcms;}         // sum of transvers energys of neutrals (lab)

double RhoEventShapes::NeutEtSumLab ( ) const

inline

Definition at line 41 of file RhoEventShapes.h.

References fneutetsumlab.

{return fneutetsumlab;}         // sum of transvers energys of neutrals (lab)

int RhoEventShapes::NNeutral ( ) const

inline

Definition at line 32 of file RhoEventShapes.h.

References fnNeut.

{return fnNeut;}                // number of neutral candidates

int RhoEventShapes::NParticles ( ) const

inline

Definition at line 30 of file RhoEventShapes.h.

References fN.

{return fN;}                    // number of particle candidates

double RhoEventShapes::Planarity

(

)

Definition at line 130 of file RhoEventShapes.cxx.

References ComputeSphericity(), and fpla.

{
  if (fpla<0.) ComputeSphericity();
  
  return fpla;
}

double RhoEventShapes::PmaxCms ( ) const

inline

Definition at line 36 of file RhoEventShapes.h.

References fpmaxcms.

{return fpmaxcms;}              // max momentum im cms system

double RhoEventShapes::PmaxLab ( ) const

inline

Definition at line 35 of file RhoEventShapes.h.

References fpmaxlab.

{return fpmaxlab;}              // max momentum in lab system

double RhoEventShapes::Ptmax ( ) const

inline

Definition at line 37 of file RhoEventShapes.h.

References fptmax.

{return fptmax;}                        // max pt (same for lab and cms)

double RhoEventShapes::PtSumCms ( ) const

inline

Definition at line 47 of file RhoEventShapes.h.

References fptsumcms.

{return fptsumcms;}                             // sum of pt in (lab)

double RhoEventShapes::PtSumLab ( ) const

inline

Definition at line 40 of file RhoEventShapes.h.

References fptsumlab.

{return fptsumlab;}                             // sum of pt in (lab)

double RhoEventShapes::Sphericity

(

)

Definition at line 121 of file RhoEventShapes.cxx.

References ComputeSphericity(), and fsph.

{
  if (fsph<0.) ComputeSphericity();
  
  return fsph;
}

double RhoEventShapes::Thrust

(

)

Definition at line 155 of file RhoEventShapes.cxx.

References Eps(), fabs(), fCmsList, fN, fthr, fThrVect, and i.

Referenced by ThrustVector().

{
  // did we already compute?
  if (fthr>-1.) return fthr; 
  
  TVector3 n0(0,0,0);
  if( fN==0 ) return -1.;
  
  int i,j;
  double pmax=0;
  
  //find starting vector as maximum momentum vector
  for (i=0;i<fN;++i)
  {
    if (fCmsList[i].Vect().Mag()>pmax)
    {
      n0=fCmsList[i].Vect();
      pmax=fCmsList[i].Vect().Mag();
    }
  }
  
  TVector3 nnew(0,0,0);

  // find thrust axis (10 iterations)
  for (i=0;i<10;++i)
  {
    for (j=0;j<fN;++j) 
      nnew += Eps(n0,fCmsList[j].Vect()) * fCmsList[j].Vect();
    
    n0=nnew.Unit();
  }
  
  fThrVect = n0;
  
  double thr=0, sum=0;
  for (i=0;i<fN;++i)
  {
    thr += fabs(fThrVect.Dot(fCmsList[i].Vect()));
    sum += fCmsList[i].Vect().Mag();
  }
  
  fthr = thr/sum;
  
  return fthr;
}

TVector3 RhoEventShapes::ThrustVector

(

)

Definition at line 203 of file RhoEventShapes.cxx.

References fthr, fThrVect, and Thrust().

{
  if (fthr<0.) Thrust();
  
  return fThrVect;
}

Member Data Documentation

double RhoEventShapes::fapl

private

Definition at line 101 of file RhoEventShapes.h.

Referenced by Aplanarity(), and ComputeSphericity().

TVector3 RhoEventShapes::fBoost

private

Definition at line 109 of file RhoEventShapes.h.

Referenced by RhoEventShapes().

std::vector<int> RhoEventShapes::fCharge

private

! List of charges of particles

Definition at line 78 of file RhoEventShapes.h.

Referenced by RhoEventShapes().

bool RhoEventShapes::fChargedOnly

private

Definition at line 112 of file RhoEventShapes.h.

double RhoEventShapes::fchrgpsumcms

private

Definition at line 98 of file RhoEventShapes.h.

Referenced by ChrgPSumCms(), and RhoEventShapes().

double RhoEventShapes::fchrgpsumlab

private

Definition at line 92 of file RhoEventShapes.h.

Referenced by ChrgPSumLab(), and RhoEventShapes().

double RhoEventShapes::fchrgptsumcms

private

Definition at line 97 of file RhoEventShapes.h.

Referenced by ChrgPtSumCms(), and RhoEventShapes().

double RhoEventShapes::fchrgptsumlab

private

Definition at line 91 of file RhoEventShapes.h.

Referenced by ChrgPtSumLab(), and RhoEventShapes().

std::vector<TLorentzVector> RhoEventShapes::fCmsList

private

! List of 4-vectors in cms frame

Definition at line 77 of file RhoEventShapes.h.

Referenced by ComputeSphericity(), FoxWolfMomH(), MultPmaxCms(), MultPminCms(), RhoEventShapes(), and Thrust().

double RhoEventShapes::fFWmom[FWMAX]

private

Definition at line 104 of file RhoEventShapes.h.

Referenced by FoxWolfMomH(), and RhoEventShapes().

bool RhoEventShapes::fFWready

private

Definition at line 105 of file RhoEventShapes.h.

Referenced by FoxWolfMomH().

std::vector<TLorentzVector> RhoEventShapes::fLabList

private

! List of 4-vectors in lab frame

Definition at line 76 of file RhoEventShapes.h.

Referenced by MultPmaxLab(), MultPminLab(), and RhoEventShapes().

int RhoEventShapes::fN

private

Definition at line 82 of file RhoEventShapes.h.

Referenced by ComputeSphericity(), FoxWolfMomH(), MultPmaxCms(), MultPmaxLab(), MultPminCms(), MultPminLab(), NParticles(), RhoEventShapes(), and Thrust().

int RhoEventShapes::fnChrg

private

Definition at line 80 of file RhoEventShapes.h.

Referenced by NCharged(), and RhoEventShapes().

double RhoEventShapes::fneutesumcms

private

Definition at line 96 of file RhoEventShapes.h.

Referenced by NeutESumCms(), and RhoEventShapes().

double RhoEventShapes::fneutesumlab

private

Definition at line 90 of file RhoEventShapes.h.

Referenced by NeutESumLab(), and RhoEventShapes().

double RhoEventShapes::fneutetsumcms

private

Definition at line 95 of file RhoEventShapes.h.

Referenced by NeutEtSumCms(), and RhoEventShapes().

double RhoEventShapes::fneutetsumlab

private

Definition at line 89 of file RhoEventShapes.h.

Referenced by NeutEtSumLab(), and RhoEventShapes().

bool RhoEventShapes::fNeutralOnly

private

Definition at line 111 of file RhoEventShapes.h.

int RhoEventShapes::fnNeut

private

Definition at line 81 of file RhoEventShapes.h.

Referenced by NNeutral(), and RhoEventShapes().

double RhoEventShapes::fpla

private

Definition at line 102 of file RhoEventShapes.h.

Referenced by ComputeSphericity(), and Planarity().

double RhoEventShapes::fpmaxcms

private

Definition at line 85 of file RhoEventShapes.h.

Referenced by PmaxCms(), and RhoEventShapes().

double RhoEventShapes::fpmaxlab

private

Definition at line 84 of file RhoEventShapes.h.

Referenced by PmaxLab(), and RhoEventShapes().

double RhoEventShapes::fptmax

private

Definition at line 86 of file RhoEventShapes.h.

Referenced by Ptmax(), and RhoEventShapes().

double RhoEventShapes::fptsumcms

private

Definition at line 94 of file RhoEventShapes.h.

Referenced by PtSumCms(), and RhoEventShapes().

double RhoEventShapes::fptsumlab

private

Definition at line 88 of file RhoEventShapes.h.

Referenced by PtSumLab(), and RhoEventShapes().

double RhoEventShapes::fsph

private

Definition at line 100 of file RhoEventShapes.h.

Referenced by ComputeSphericity(), and Sphericity().

double RhoEventShapes::fthr

private

Definition at line 107 of file RhoEventShapes.h.

Referenced by Thrust(), and ThrustVector().

TVector3 RhoEventShapes::fThrVect

private

Definition at line 108 of file RhoEventShapes.h.

Referenced by RhoEventShapes(), Thrust(), and ThrustVector().

---

The documentation for this class was generated from the following files:

• RhoEventShapes.h
• RhoEventShapes.cxx