Main Page | Namespace List | Class Hierarchy | Alphabetical List | Class List | Directories | File List | Namespace Members | Class Members | File Members | Related Pages

GammaShape Class Reference

#include <GammaShape.h>

List of all members.

Public Member Functions

double de (const double eorigin, const double peak)
double de (const double eorigin, const double peak)
double fitm (double x, double pk, double width, double tail, double mean)
double fitm (double x, double pk, double width, double tail, double mean)
 GammaShape (const HepLorentzVector p4, const double dphi, const double dthe, const double de)
 GammaShape (const double eorigin, const double epeak)
 GammaShape (const GammaShape &)
 GammaShape ()
 GammaShape (const HepLorentzVector p4, const double dphi, const double dthe, const double de)
 GammaShape (const double eorigin, const double epeak)
 GammaShape (const GammaShape &)
 GammaShape ()
double getde () const
double getde () const
double getdphi () const
double getdphi () const
double getdthe () const
double getdthe () const
double gete () const
double gete () const
double getpeak () const
double getpeak () const
double gettail () const
double gettail () const
double getwidth () const
double getwidth () const
double mean (const double epeak)
double mean (const double epeak)
GammaShape & operator= (const GammaShape &)
GammaShape & operator= (const GammaShape &)
double peak (const double fitvalue)
double peak (const double fitvalue)
void setde (const double de)
void setde (const double de)
void setdphi (const double dphi)
void setdphi (const double dphi)
void setdtheta (const double dtheta)
void setdtheta (const double dtheta)
void sete (const double e)
void sete (const double e)
void setpeak (const double peak)
void setpeak (const double peak)
void settail (const double tail)
void settail (const double tail)
void setwidth (const double width)
void setwidth (const double width)
double tail (const double epeak)
double tail (const double epeak)
double width (const double epeak)
double width (const double epeak)
 ~GammaShape ()
 ~GammaShape ()

Private Attributes

double m_de
double m_dphi
double m_dtheta
double m_e
double m_peak
double m_tail
double m_width

Constructor & Destructor Documentation

GammaShape::GammaShape (   ) 
  

00005                        {
00006   m_tail = 0;
00007   m_width = 0;
00008   m_peak = 0;
00009   m_e = 0;
00010   m_de = 0;
00011   m_dphi = 0;
00012   m_dtheta = 0;
00013 }

GammaShape::~GammaShape (   )  [inline]
  

00032 {;}

GammaShape::GammaShape (  const GammaShape &   ) 
  

00015                                                    {
00016   m_tail = gammashape.m_tail;
00017   m_width = gammashape.m_width;
00018   m_peak = gammashape.m_peak;
00019   m_e = gammashape.m_e;
00020   m_de = gammashape.m_de;
00021   m_dphi = gammashape.m_dphi;
00022   m_dtheta = gammashape.m_dtheta;
00023 }

GammaShape::GammaShape (  const double  eorigin,
const double  epeak
) 
  

00036                                                               {
00037   m_tail = tail(peak);
00038   m_width = width(peak);
00039   m_peak = peak;
00040   m_e = eorigin;
00041   m_de = de(eorigin, peak);
00042   m_dphi = 0;
00043   m_dtheta = 0;
00044 }

GammaShape::GammaShape (  const HepLorentzVector  p4,
const double  dphi,
const double  dthe,
const double  de
) 
  

00047                                                                                                        {
00048   m_tail = 0;
00049   m_width = 0;
00050   m_peak = 0;
00051   m_e = 0;
00052   m_de = 0;
00053   m_dphi = dphi;
00054   m_dtheta = dthe;
00055 }

GammaShape::GammaShape (   ) 
 

GammaShape::~GammaShape (   )  [inline]
  

00032 {;}

GammaShape::GammaShape (  const GammaShape &   ) 
 

GammaShape::GammaShape (  const double  eorigin,
const double  epeak
) 
 

GammaShape::GammaShape (  const HepLorentzVector  p4,
const double  dphi,
const double  dthe,
const double  de
) 
 

Member Function Documentation

double GammaShape::de (  const double  eorigin,
const double  peak
) 
 

double GammaShape::de (  const double  eorigin,
const double  peak
) 
  

00092                                                  {
00093   double f1;
00094   double tail1 = tail(peak);
00095   double width1 = width(peak);
00096   double mean1 = mean(peak);
00097   double qa = tail1*sqrt(log(4));
00098   double qb = sinh(qa)/qa;
00099   double qx = (eorigin-peak)/width1*qb*tail1;
00100 //cout<<"qx ="<<qx<<endl;
00101 //  f1 = (eorigin-peak)/fitm(eorigin,peak, width1, tail1);
00102 
00103   if(fabs(qx)>0.001){
00104     f1 = (eorigin-peak)/fitm(eorigin,peak, width1, tail1, mean1);
00105 //    cout<<"dynamic"<<endl;
00106   }
00107   else { 
00108 //    f1=width1/qb;
00109     f1 = width1/(qb*tail1);
00110 //    cout<<"gaus"<<endl;
00111   }
00112 
00113  // else f1 = 0;
00114 //  return f1*eorigin;
00115   return f1;
00116 }

double GammaShape::fitm (  double  x,
double  pk,
double  width,
double  tail,
double  mean
) 
 

double GammaShape::fitm (  double  x,
double  pk,
double  width,
double  tail,
double  mean
) 
  

00122                                                                                    {
00123     
00124   double qa = tail*sqrt(log(4));
00125   double qb = sinh(qa)/qa;
00126   double qx = ((x-pk)/pk)/width*qb;
00127   double qy = 1 + tail * qx;
00128   double qc;
00129   if(qy < 1.e-7) qc = 15;
00130   else
00131     qc = log(qy)/tail;
00132   double val = qc;
00133   return val-mean;
00134  }

double GammaShape::getde (   )  const [inline]
  

00061 {return m_de;}

double GammaShape::getde (   )  const [inline]
  

00061 {return m_de;}

double GammaShape::getdphi (   )  const [inline]
  

00062 {return m_dphi;}

double GammaShape::getdphi (   )  const [inline]
  

00062 {return m_dphi;}

double GammaShape::getdthe (   )  const [inline]
  

00063 {return m_dtheta;}

double GammaShape::getdthe (   )  const [inline]
  

00063 {return m_dtheta;}

double GammaShape::gete (   )  const [inline]
  

00060 {return m_e;}

double GammaShape::gete (   )  const [inline]
  

00060 {return m_e;}

double GammaShape::getpeak (   )  const [inline]
  

00059 {return m_peak;}

double GammaShape::getpeak (   )  const [inline]
  

00059 {return m_peak;}

double GammaShape::gettail (   )  const [inline]
  

00057 {return m_tail;}

double GammaShape::gettail (   )  const [inline]
  

00057 {return m_tail;}

double GammaShape::getwidth (   )  const [inline]
  

00058 {return m_width;}

double GammaShape::getwidth (   )  const [inline]
  

00058 {return m_width;}

double GammaShape::mean (  const double  epeak  ) 
 

double GammaShape::mean (  const double  epeak  ) 
  

00088                                          {
00089      return -0.7984 + 0.7759*peak - 0.6356*peak*peak + 0.1476*peak*peak*peak; 
00090 }

GammaShape& GammaShape::operator= (  const GammaShape &   ) 
 

GammaShape & GammaShape::operator= (  const GammaShape &   ) 
  

00025                                                                    {
00026   m_tail = gammashape.m_tail;
00027   m_width = gammashape.m_width;
00028   m_peak = gammashape.m_peak;
00029   m_e = gammashape.m_e;
00030   m_de = gammashape.m_de;
00031   m_dphi = gammashape.m_dphi;
00032   m_dtheta = gammashape.m_dtheta;
00033   return (*this);
00034 }

double GammaShape::peak (  const double  fitvalue  ) 
 

double GammaShape::peak (  const double  fitvalue  ) 
 

void GammaShape::setde (  const double  de  )  [inline]
  

00052 {m_de = de;}

void GammaShape::setde (  const double  de  )  [inline]
  

00052 {m_de = de;}

void GammaShape::setdphi (  const double  dphi  )  [inline]
  

00053 {m_dphi = dphi;}

void GammaShape::setdphi (  const double  dphi  )  [inline]
  

00053 {m_dphi = dphi;}

void GammaShape::setdtheta (  const double  dtheta  )  [inline]
  

00054 {m_dtheta = dtheta;}

void GammaShape::setdtheta (  const double  dtheta  )  [inline]
  

00054 {m_dtheta = dtheta;}

void GammaShape::sete (  const double  e  )  [inline]
  

00051 {m_e = e;}

void GammaShape::sete (  const double  e  )  [inline]
  

00051 {m_e = e;}

void GammaShape::setpeak (  const double  peak  )  [inline]
  

00050 {m_peak = peak;}

void GammaShape::setpeak (  const double  peak  )  [inline]
  

00050 {m_peak = peak;}

void GammaShape::settail (  const double  tail  )  [inline]
  

00048 {m_tail = tail;}

void GammaShape::settail (  const double  tail  )  [inline]
  

00048 {m_tail = tail;}

void GammaShape::setwidth (  const double  width  )  [inline]
  

00049 {m_width = width;}

void GammaShape::setwidth (  const double  width  )  [inline]
  

00049 {m_width = width;}

double GammaShape::tail (  const double  epeak  ) 
 

double GammaShape::tail (  const double  epeak  ) 
  

00080                                         {
00081     return -0.7931 + 0.8026*peak - 0.6694*peak*peak + 0.1476*peak*peak*peak;
00082 }

double GammaShape::width (  const double  epeak  ) 
 

double GammaShape::width (  const double  epeak  ) 
  

00084                                          {
00085      return 0.002496/sqrt(peak) + 0.01861 + 0.0001788*peak;
00086 }

Member Data Documentation

double GammaShape::m_de [private]
 

double GammaShape::m_dphi [private]
 

double GammaShape::m_dtheta [private]
 

double GammaShape::m_e [private]
 

double GammaShape::m_peak [private]
 

double GammaShape::m_tail [private]
 

double GammaShape::m_width [private]
 

The documentation for this class was generated from the following files:
Generated on Wed Feb 2 16:16:27 2011 for BOSS6.5.5 by  doxygen 1.3.9.1