#include <EventShape.h>

Public Member Functions
	EventShape (const CandList &l, TLorentzVector cms)

	~EventShape ()

int	NParticles () const

int	NCharged () const

int	NNeutral () const

double	PmaxLab () const

double	PmaxCms () const

double	PminLab () const

double	PminCms () const

double	Ptmax () const

double	Ptmin () const

double	PRapmax () 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)

int	MultPtminLab (double ptmin)

int	MultPtmaxLab (double ptmax)

int	MultPtminCms (double ptmin)

int	MultPtmaxCms (double ptmax)

int	MultNeutEminLab (double emin)

int	MultNeutEmaxLab (double emax)

int	MultNeutEminCms (double emin)

int	MultNeutEmaxCms (double emax)

int	MultChrgPminLab (double pmin)

int	MultChrgPmaxLab (double pmax)

int	MultChrgPminCms (double pmin)

int	MultChrgPmaxCms (double pmax)

double	SumPminLab (double pmin)

double	SumPmaxLab (double pmax)

double	SumPminCms (double pmin)

double	SumPmaxCms (double pmax)

double	SumPtminLab (double ptmin)

double	SumPtmaxLab (double ptmax)

double	SumPtminCms (double ptmin)

double	SumPtmaxCms (double ptmax)

double	SumNeutEminLab (double emin)

double	SumNeutEmaxLab (double emax)

double	SumNeutEminCms (double emin)

double	SumNeutEmaxCms (double emax)

double	SumChrgPminLab (double pmin)

double	SumChrgPmaxLab (double pmax)

double	SumChrgPminCms (double pmin)

double	SumChrgPmaxCms (double pmax)

double	Sphericity ()

double	Aplanarity ()

double	Planarity ()

double	Circularity ()

double	FoxWolfMomH (int order)

double	FoxWolfMomR (int order)

double	Thrust ()

TVector3	ThrustVector ()

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

std::vector< TLorentzVector >	fCmsList

std::vector< int >	fCharge

int	fnChrg

int	fnNeut

int	fN

double	fpmaxlab

double	fpmaxcms

double	fpminlab

double	fpmincms

double	fptmax

double	fptmin

double	fprapmax

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	fcir

double	fFWmom [FWMAX]

bool	fFWready

double	fthr

TVector3	fThrVect

TVector3	fBoost

Detailed Description

Definition at line 9 of file EventShape.h.

Constructor & Destructor Documentation

EventShape::EventShape	(	const CandList &	l,
		TLorentzVector	cms
	)

Definition at line 152 of file EventShape.h.

References fBoost, fCharge, fchrgpsumcms, fchrgpsumlab, fchrgptsumcms, fchrgptsumlab, fCmsList, fFWmom, fLabList, fN, fnChrg, fneutesumcms, fneutesumlab, fneutetsumcms, fneutetsumlab, fnNeut, fpmaxcms, fpmaxlab, fpmincms, fpminlab, fprapmax, fptmax, fptmin, fptsumcms, fptsumlab, fThrVect, FWMAX, and i.

                                                             : 
   fnChrg(0), fnNeut(0), fN(0), fpmaxlab(0.), fpmaxcms(0.),fpminlab(0.), fpmincms(0.), fptmax(0.), fptmin(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.), fcir(-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.;
   double pmin=1000., ptmin=1000., pmincms=1000.;
   double prapmax=0.;
   
   fLabList.clear();
   fCmsList.clear();
   
   for (i=0;i<l.size();++i)
   {
         TLorentzVector lv(l[i].P4());
         int chrg(l[i].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 minimum momenta in lab
         if (lv.P()<pmin) pmin=lv.P();
         if (lv.Pt()<ptmin) ptmin=lv.Pt();
         
         // cache pseudorapidity
         if (lv.PseudoRapidity()>prapmax) prapmax = lv.PseudoRapidity();
         
         // 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();
         // cache minimum momenta in cms
         if (lv.P()<pmincms) pmincms=lv.P();
   }
   
   fpmaxlab = pmax;
   fptmax   = ptmax;
   fpmaxcms = pmaxcms;
   fpminlab = pmin;
   fptmin   = ptmin;
   fpmincms = pmincms;
   fprapmax = prapmax;
 }

EventShape::~EventShape ( )

inline

Definition at line 12 of file EventShape.h.

12 {};

Member Function Documentation

double EventShape::Aplanarity ( )

Definition at line 299 of file EventShape.h.

References ComputeSphericity(), and fapl.

Referenced by softtrigger_kin5(), and writeTuple().

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

double EventShape::ChrgPSumCms ( ) const

inline

Definition at line 40 of file EventShape.h.

References fchrgpsumcms.

Referenced by softtrigger_kin5(), softtrigger_toy12(), toy_core(), and writeTuple().

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

EventShape::fchrgpsumcms

double fchrgpsumcms

Definition: EventShape.h:134

double EventShape::ChrgPSumLab ( ) const

inline

Definition at line 33 of file EventShape.h.

References fchrgpsumlab.

Referenced by softtrigger_kin5().

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

EventShape::fchrgpsumlab

double fchrgpsumlab

Definition: EventShape.h:128

double EventShape::ChrgPtSumCms ( ) const

inline

Definition at line 39 of file EventShape.h.

References fchrgptsumcms.

Referenced by softtrigger_kin5().

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

EventShape::fchrgptsumcms

double fchrgptsumcms

Definition: EventShape.h:133

double EventShape::ChrgPtSumLab ( ) const

inline

Definition at line 32 of file EventShape.h.

References fchrgptsumlab.

Referenced by softtrigger_kin5(), softtrigger_toy12(), toy_core(), and writeTuple().

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

EventShape::fchrgptsumlab

double fchrgptsumlab

Definition: EventShape.h:127

double EventShape::Circularity ( )

Definition at line 308 of file EventShape.h.

References ComputeSphericity(), and fcir.

Referenced by softtrigger_kin5(), and writeTuple().

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

void EventShape::ComputeSphericity ( )

private

Definition at line 247 of file EventShape.h.

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

Referenced by Aplanarity(), Circularity(), 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));           // 3/2 (lam_2+lam_3)
   fapl = 1.5 * eiv(2,2);                                        // 3/2 lam_3
   fpla = eiv(1,1) - eiv(2,2);                           // lam_2-lam_3
   fcir = 2.*eiv(1,1)/(eiv(1,1)+eiv(0,0));       // 2* min(lam_1, lam_2)/(lam_1+lam_2)
 }

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

inlineprivate

Definition at line 105 of file EventShape.h.

Referenced by Thrust().

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

v1

TVector3 v1

Definition: bump_analys.C:40

v2

TVector3 v2

Definition: bump_analys.C:40

double EventShape::FoxWolfMomH ( int order )

Definition at line 404 of file EventShape.h.

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 EventShape::FoxWolfMomR ( int order )

Definition at line 462 of file EventShape.h.

References FoxWolfMomH().

Referenced by softtrigger_kin5(), softtrigger_toy12(), toy_core(), and writeTuple().

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

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

private

Definition at line 381 of file EventShape.h.

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 EventShape::MultChrgPmaxCms ( double pmax )

Definition at line 641 of file EventShape.h.

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

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

int EventShape::MultChrgPmaxLab ( double pmax )

Definition at line 619 of file EventShape.h.

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

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

int EventShape::MultChrgPminCms ( double pmin )

Definition at line 630 of file EventShape.h.

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

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

int EventShape::MultChrgPminLab ( double pmin )

Definition at line 608 of file EventShape.h.

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

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

int EventShape::MultNeutEmaxCms ( double emax )

Definition at line 595 of file EventShape.h.

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

 {
   int cnt=0;
   for (int i=0;i<fN;++i) 
         if (fCharge[i]==0 && fCmsList[i].E()<emax) cnt++;
  
   return cnt;
 }

int EventShape::MultNeutEmaxLab ( double emax )

Definition at line 573 of file EventShape.h.

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

 {
   int cnt=0;
   for (int i=0;i<fN;++i) 
         if (fCharge[i]==0 && fLabList[i].E()<emax) cnt++;
  
   return cnt;
 }

int EventShape::MultNeutEminCms ( double emin )

Definition at line 584 of file EventShape.h.

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

 {
   int cnt=0;
   for (int i=0;i<fN;++i) 
         if (fCharge[i]==0 && fCmsList[i].E()>emin) cnt++;
  
   return cnt;
 }

int EventShape::MultNeutEminLab ( double emin )

Definition at line 562 of file EventShape.h.

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

 {
   int cnt=0;
   for (int i=0;i<fN;++i) 
         if (fCharge[i]==0 && fLabList[i].E()>emin) cnt++;
  
   return cnt;
 }

int EventShape::MultPmaxCms ( double pmax )

Definition at line 504 of file EventShape.h.

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 EventShape::MultPmaxLab ( double pmax )

Definition at line 482 of file EventShape.h.

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 EventShape::MultPminCms ( double pmin )

Definition at line 493 of file EventShape.h.

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

Referenced by softtrigger_kin5(), softtrigger_toy12(), and toy_core().

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

int EventShape::MultPminLab ( double pmin )

Definition at line 471 of file EventShape.h.

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

Referenced by softtrigger_kin5().

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

int EventShape::MultPtmaxCms ( double ptmax )

Definition at line 549 of file EventShape.h.

References cnt, fCmsList, fN, and i.

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

int EventShape::MultPtmaxLab ( double ptmax )

Definition at line 527 of file EventShape.h.

References cnt, fLabList, fN, and i.

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

int EventShape::MultPtminCms ( double ptmin )

Definition at line 538 of file EventShape.h.

References cnt, fCmsList, fN, and i.

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

int EventShape::MultPtminLab ( double ptmin )

Definition at line 516 of file EventShape.h.

References cnt, fLabList, fN, and i.

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

int EventShape::NCharged ( ) const

inline

Definition at line 16 of file EventShape.h.

References fnChrg.

Referenced by softtrigger_kin5(), softtrigger_toy12(), toy_core(), and writeTuple().

16 {return fnChrg;} // number of charged candidates

EventShape::fnChrg

int fnChrg

Definition: EventShape.h:112

double EventShape::NeutESumCms ( ) const

inline

Definition at line 38 of file EventShape.h.

References fneutesumcms.

Referenced by softtrigger_kin5(), softtrigger_toy12(), and toy_core().

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

EventShape::fneutesumcms

double fneutesumcms

Definition: EventShape.h:132

double EventShape::NeutESumLab ( ) const

inline

Definition at line 31 of file EventShape.h.

References fneutesumlab.

Referenced by softtrigger_kin5().

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

EventShape::fneutesumlab

double fneutesumlab

Definition: EventShape.h:126

double EventShape::NeutEtSumCms ( ) const

inline

Definition at line 37 of file EventShape.h.

References fneutetsumcms.

Referenced by softtrigger_kin5().

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

EventShape::fneutetsumcms

double fneutetsumcms

Definition: EventShape.h:131

double EventShape::NeutEtSumLab ( ) const

inline

Definition at line 30 of file EventShape.h.

References fneutetsumlab.

Referenced by softtrigger_kin5().

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

EventShape::fneutetsumlab

double fneutetsumlab

Definition: EventShape.h:125

int EventShape::NNeutral ( ) const

inline

Definition at line 17 of file EventShape.h.

References fnNeut.

Referenced by softtrigger_kin5(), softtrigger_toy12(), toy_core(), and writeTuple().

17 {return fnNeut;} // number of neutral candidates

EventShape::fnNeut

int fnNeut

Definition: EventShape.h:113

int EventShape::NParticles ( ) const

inline

Definition at line 15 of file EventShape.h.

References fN.

Referenced by softtrigger_kin5(), softtrigger_toy12(), toy_core(), and writeTuple().

15 {return fN;} // number of particle candidates

EventShape::fN

int fN

Definition: EventShape.h:114

double EventShape::Planarity ( )

Definition at line 290 of file EventShape.h.

References ComputeSphericity(), and fpla.

Referenced by softtrigger_kin5(), and writeTuple().

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

double EventShape::PmaxCms ( ) const

inline

Definition at line 21 of file EventShape.h.

References fpmaxcms.

Referenced by softtrigger_kin5(), softtrigger_toy12(), toy_core(), and writeTuple().

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

EventShape::fpmaxcms

double fpmaxcms

Definition: EventShape.h:117

double EventShape::PmaxLab ( ) const

inline

Definition at line 20 of file EventShape.h.

References fpmaxlab.

Referenced by softtrigger_kin5(), softtrigger_toy12(), toy_core(), and writeTuple().

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

EventShape::fpmaxlab

double fpmaxlab

Definition: EventShape.h:116

double EventShape::PminCms ( ) const

inline

Definition at line 23 of file EventShape.h.

References fpmincms.

Referenced by softtrigger_kin5(), softtrigger_toy12(), toy_core(), and writeTuple().

23 {return fpmincms;} // min momentum im cms system

EventShape::fpmincms

double fpmincms

Definition: EventShape.h:119

double EventShape::PminLab ( ) const

inline

Definition at line 22 of file EventShape.h.

References fpminlab.

Referenced by softtrigger_kin5().

22 {return fpminlab;} // min momentum in lab system

EventShape::fpminlab

double fpminlab

Definition: EventShape.h:118

double EventShape::PRapmax ( ) const

inline

Definition at line 26 of file EventShape.h.

References fprapmax.

Referenced by softtrigger_kin5().

26 {return fprapmax;} // max pseudorapidity (lab)

EventShape::fprapmax

double fprapmax

Definition: EventShape.h:122

double EventShape::Ptmax ( ) const

inline

Definition at line 24 of file EventShape.h.

References fptmax.

Referenced by softtrigger_kin5(), softtrigger_toy12(), toy_core(), and writeTuple().

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

EventShape::fptmax

double fptmax

Definition: EventShape.h:120

double EventShape::Ptmin ( ) const

inline

Definition at line 25 of file EventShape.h.

References fptmin.

Referenced by softtrigger_kin5().

25 {return fptmin;} // min pt (same for lab and cms)

EventShape::fptmin

double fptmin

Definition: EventShape.h:121

double EventShape::PtSumCms ( ) const

inline

Definition at line 36 of file EventShape.h.

References fptsumcms.

Referenced by softtrigger_kin5().

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

EventShape::fptsumcms

double fptsumcms

Definition: EventShape.h:130

double EventShape::PtSumLab ( ) const

inline

Definition at line 29 of file EventShape.h.

References fptsumlab.

Referenced by softtrigger_kin5(), softtrigger_toy12(), toy_core(), and writeTuple().

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

EventShape::fptsumlab

double fptsumlab

Definition: EventShape.h:124

double EventShape::Sphericity ( )

Definition at line 281 of file EventShape.h.

References ComputeSphericity(), and fsph.

Referenced by softtrigger_kin5(), and writeTuple().

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

double EventShape::SumChrgPmaxCms ( double pmax )

Definition at line 825 of file EventShape.h.

References fCharge, fCmsList, fLabList, fN, i, and P.

 {
   double sum=0;
   for (int i=0;i<fN;++i) 
         if (fCharge[i]!=0 && fCmsList[i].P()<pmax) sum+=fLabList[i].P();
  
   return sum;
 }

double EventShape::SumChrgPmaxLab ( double pmax )

Definition at line 803 of file EventShape.h.

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

 {
   double sum=0;
   for (int i=0;i<fN;++i) 
         if (fCharge[i]!=0 && fLabList[i].P()<pmax) sum+=fLabList[i].P();
  
   return sum;
 }

double EventShape::SumChrgPminCms ( double pmin )

Definition at line 814 of file EventShape.h.

References fCharge, fCmsList, fLabList, fN, i, and P.

Referenced by softtrigger_kin5(), softtrigger_toy12(), and toy_core().

 {
   double sum=0;
   for (int i=0;i<fN;++i) 
         if (fCharge[i]!=0 && fCmsList[i].P()>pmin) sum+=fLabList[i].P();
  
   return sum;
 }

double EventShape::SumChrgPminLab ( double pmin )

Definition at line 792 of file EventShape.h.

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

Referenced by softtrigger_kin5().

 {
   double sum=0;
   for (int i=0;i<fN;++i) 
         if (fCharge[i]!=0 && fLabList[i].P()>pmin) sum+=fLabList[i].P();
  
   return sum;
 }

double EventShape::SumNeutEmaxCms ( double emax )

Definition at line 779 of file EventShape.h.

References fCharge, fCmsList, fLabList, fN, and i.

 {
   double sum=0;
   for (int i=0;i<fN;++i) 
         if (fCharge[i]==0 && fCmsList[i].E()<emax) sum+=fLabList[i].E();
  
   return sum;
 }

double EventShape::SumNeutEmaxLab ( double emax )

Definition at line 757 of file EventShape.h.

References fCharge, fLabList, fN, and i.

 {
   double sum=0;
   for (int i=0;i<fN;++i) 
         if (fCharge[i]==0 && fLabList[i].E()<emax) sum+=fLabList[i].E();
  
   return sum;
 }

double EventShape::SumNeutEminCms ( double emin )

Definition at line 768 of file EventShape.h.

References fCharge, fCmsList, fLabList, fN, and i.

Referenced by softtrigger_kin5().

 {
   double sum=0;
   for (int i=0;i<fN;++i) 
         if (fCharge[i]==0 && fCmsList[i].E()>emin) sum+=fLabList[i].E();
  
   return sum;
 }

double EventShape::SumNeutEminLab ( double emin )

Definition at line 746 of file EventShape.h.

References fCharge, fLabList, fN, and i.

Referenced by softtrigger_kin5().

 {
   double sum=0;
   for (int i=0;i<fN;++i) 
         if (fCharge[i]==0 && fLabList[i].E()>emin) sum+=fLabList[i].E();
  
   return sum;
 }

double EventShape::SumPmaxCms ( double pmax )

Definition at line 687 of file EventShape.h.

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

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

double EventShape::SumPmaxLab ( double pmax )

Definition at line 665 of file EventShape.h.

References fLabList, fN, i, and P.

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

double EventShape::SumPminCms ( double pmin )

Definition at line 676 of file EventShape.h.

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

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

double EventShape::SumPminLab ( double pmin )

Definition at line 654 of file EventShape.h.

References fLabList, fN, i, and P.

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

double EventShape::SumPtmaxCms ( double ptmax )

Definition at line 733 of file EventShape.h.

References fCmsList, fLabList, fN, and i.

 {
   double sum=0;
   for (int i=0;i<fN;++i) 
         if (fCmsList[i].Pt()<ptmax) sum+=fLabList[i].Pt();
  
   return sum;
 }

double EventShape::SumPtmaxLab ( double ptmax )

Definition at line 711 of file EventShape.h.

References fLabList, fN, and i.

 {
   double sum=0;
   for (int i=0;i<fN;++i) 
         if (fLabList[i].Pt()<ptmax) sum+=fLabList[i].Pt();
  
   return sum;
 }

double EventShape::SumPtminCms ( double ptmin )

Definition at line 722 of file EventShape.h.

References fCmsList, fLabList, fN, and i.

 {
   double sum=0;
   for (int i=0;i<fN;++i) 
         if (fCmsList[i].Pt()>ptmin) sum+=fLabList[i].Pt();
  
   return sum;
 }

double EventShape::SumPtminLab ( double ptmin )

Definition at line 700 of file EventShape.h.

References fLabList, fN, and i.

Referenced by softtrigger_kin5().

 {
   double sum=0;
   for (int i=0;i<fN;++i) 
         if (fLabList[i].Pt()>ptmin) sum+=fLabList[i].Pt();
  
   return sum;
 }

double EventShape::Thrust ( )

Definition at line 324 of file EventShape.h.

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

Referenced by softtrigger_kin5(), ThrustVector(), and writeTuple().

 {
   // 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 EventShape::ThrustVector ( )

Definition at line 372 of file EventShape.h.

References fthr, fThrVect, and Thrust().

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

Member Data Documentation

double EventShape::fapl

private

Definition at line 137 of file EventShape.h.

Referenced by Aplanarity(), and ComputeSphericity().

TVector3 EventShape::fBoost

private

Definition at line 146 of file EventShape.h.

Referenced by EventShape().

std::vector<int> EventShape::fCharge

private

Definition at line 110 of file EventShape.h.

Referenced by EventShape(), MultChrgPmaxCms(), MultChrgPmaxLab(), MultChrgPminCms(), MultChrgPminLab(), MultNeutEmaxCms(), MultNeutEmaxLab(), MultNeutEminCms(), MultNeutEminLab(), SumChrgPmaxCms(), SumChrgPmaxLab(), SumChrgPminCms(), SumChrgPminLab(), SumNeutEmaxCms(), SumNeutEmaxLab(), SumNeutEminCms(), and SumNeutEminLab().

double EventShape::fchrgpsumcms

private

Definition at line 134 of file EventShape.h.

Referenced by ChrgPSumCms(), and EventShape().

double EventShape::fchrgpsumlab

private

Definition at line 128 of file EventShape.h.

Referenced by ChrgPSumLab(), and EventShape().

double EventShape::fchrgptsumcms

private

Definition at line 133 of file EventShape.h.

Referenced by ChrgPtSumCms(), and EventShape().

double EventShape::fchrgptsumlab

private

Definition at line 127 of file EventShape.h.

Referenced by ChrgPtSumLab(), and EventShape().

double EventShape::fcir

private

Definition at line 139 of file EventShape.h.

Referenced by Circularity(), and ComputeSphericity().

std::vector<TLorentzVector> EventShape::fCmsList

private

Definition at line 109 of file EventShape.h.

Referenced by ComputeSphericity(), EventShape(), FoxWolfMomH(), MultChrgPmaxCms(), MultChrgPminCms(), MultNeutEmaxCms(), MultNeutEminCms(), MultPmaxCms(), MultPminCms(), MultPtmaxCms(), MultPtminCms(), SumChrgPmaxCms(), SumChrgPminCms(), SumNeutEmaxCms(), SumNeutEminCms(), SumPmaxCms(), SumPminCms(), SumPtmaxCms(), SumPtminCms(), and Thrust().

double EventShape::fFWmom[FWMAX]

private

Definition at line 141 of file EventShape.h.

Referenced by EventShape(), and FoxWolfMomH().

bool EventShape::fFWready

private

Definition at line 142 of file EventShape.h.

Referenced by FoxWolfMomH().

std::vector<TLorentzVector> EventShape::fLabList

private

Definition at line 108 of file EventShape.h.

Referenced by EventShape(), MultChrgPmaxLab(), MultChrgPminLab(), MultNeutEmaxLab(), MultNeutEminLab(), MultPmaxLab(), MultPminLab(), MultPtmaxLab(), MultPtminLab(), SumChrgPmaxCms(), SumChrgPmaxLab(), SumChrgPminCms(), SumChrgPminLab(), SumNeutEmaxCms(), SumNeutEmaxLab(), SumNeutEminCms(), SumNeutEminLab(), SumPmaxCms(), SumPmaxLab(), SumPminCms(), SumPminLab(), SumPtmaxCms(), SumPtmaxLab(), SumPtminCms(), and SumPtminLab().

int EventShape::fN

private

int EventShape::fnChrg

private

Definition at line 112 of file EventShape.h.

Referenced by EventShape(), and NCharged().

double EventShape::fneutesumcms

private

Definition at line 132 of file EventShape.h.

Referenced by EventShape(), and NeutESumCms().

double EventShape::fneutesumlab

private

Definition at line 126 of file EventShape.h.

Referenced by EventShape(), and NeutESumLab().

double EventShape::fneutetsumcms

private

Definition at line 131 of file EventShape.h.

Referenced by EventShape(), and NeutEtSumCms().

double EventShape::fneutetsumlab

private

Definition at line 125 of file EventShape.h.

Referenced by EventShape(), and NeutEtSumLab().

int EventShape::fnNeut

private

Definition at line 113 of file EventShape.h.

Referenced by EventShape(), and NNeutral().

double EventShape::fpla

private

Definition at line 138 of file EventShape.h.

Referenced by ComputeSphericity(), and Planarity().

double EventShape::fpmaxcms

private

Definition at line 117 of file EventShape.h.

Referenced by EventShape(), and PmaxCms().

double EventShape::fpmaxlab

private

Definition at line 116 of file EventShape.h.

Referenced by EventShape(), and PmaxLab().

double EventShape::fpmincms

private

Definition at line 119 of file EventShape.h.

Referenced by EventShape(), and PminCms().

double EventShape::fpminlab

private

Definition at line 118 of file EventShape.h.

Referenced by EventShape(), and PminLab().

double EventShape::fprapmax

private

Definition at line 122 of file EventShape.h.

Referenced by EventShape(), and PRapmax().

double EventShape::fptmax

private

Definition at line 120 of file EventShape.h.

Referenced by EventShape(), and Ptmax().

double EventShape::fptmin

private

Definition at line 121 of file EventShape.h.

Referenced by EventShape(), and Ptmin().

double EventShape::fptsumcms

private

Definition at line 130 of file EventShape.h.

Referenced by EventShape(), and PtSumCms().

double EventShape::fptsumlab

private

Definition at line 124 of file EventShape.h.

Referenced by EventShape(), and PtSumLab().

double EventShape::fsph

private

Definition at line 136 of file EventShape.h.

Referenced by ComputeSphericity(), and Sphericity().

double EventShape::fthr

private

Definition at line 144 of file EventShape.h.

Referenced by Thrust(), and ThrustVector().

TVector3 EventShape::fThrVect

private

Definition at line 145 of file EventShape.h.

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

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

EventShape.h

Public Member Functions

Private Member Functions

Private Attributes

Detailed Description

Constructor & Destructor Documentation

Member Function Documentation

Member Data Documentation