EvtSpinDensity Class Reference

#include <EvtSpinDensity.hh>

Public Member Functions
	EvtSpinDensity (const EvtSpinDensity &density)
EvtSpinDensity &	operator= (const EvtSpinDensity &density)
virtual	~EvtSpinDensity ()
	EvtSpinDensity ()
void	SetDim (int n)
int	GetDim () const
void	Set (int i, int j, const EvtComplex &rhoij)
const EvtComplex &	Get (int i, int j) const
double	NormalizedProb (const EvtSpinDensity &d)
void	SetDiag (int n)
int	Check ()

Friends
std::ostream &	operator<< (std::ostream &s, const EvtSpinDensity &d)

Detailed Description

Definition at line 27 of file EvtSpinDensity.hh.

Constructor & Destructor Documentation

EvtSpinDensity() [1/2]

EvtSpinDensity::EvtSpinDensity

(

const EvtSpinDensity &

density

)

Definition at line 33 of file EvtSpinDensity.cc.

                                                           {
  dim=0;
  rho=0;
 
  int i,j;
  SetDim(density.dim);
 
  for(i=0;i<dim;i++){
    for(j=0;j<dim;j++){
      rho[i][j]=density.rho[i][j];
    }
  }
}

~EvtSpinDensity()

EvtSpinDensity::~EvtSpinDensity ( )

virtual

Definition at line 61 of file EvtSpinDensity.cc.

                               {
  if (dim!=0){
    int i;
    for(i=0;i<dim;i++) delete [] rho[i];
  }
  
  delete [] rho;
}

EvtSpinDensity() [2/2]

EvtSpinDensity::EvtSpinDensity

(

)

Definition at line 70 of file EvtSpinDensity.cc.

                              {
  dim=0;
  rho=0;
}

Member Function Documentation

Check()

int EvtSpinDensity::Check

(

)

Definition at line 150 of file EvtSpinDensity.cc.

                         {
 
  if (dim<1) {
    report(ERROR,"EvtGen")<<"dim="<<dim<<"in SpinDensity::Check"<<endl;
  }
 
  int i,j;
 
  for(i=0;i<dim;i++){
 
    if (real(rho[i][i])<0.0) return 0;
    if (imag(rho[i][i])*1000000.0>abs(rho[i][i])) {
      report(INFO,"EvtGen") << "Failing 1"<<endl;
      return 0;
    }
  }
 
  for(i=0;i<dim;i++){
    for(j=i+1;j<dim;j++){
      if (fabs(real(rho[i][j]-rho[j][i]))>
      0.00000001*(abs(rho[i][i])+abs(rho[j][j]))) {
    report(INFO,"EvtGen") << "Failing 2"<<endl;
    return 0;
      }
      if (fabs(imag(rho[i][j]+rho[j][i]))>
      0.00000001*(abs(rho[i][i])+abs(rho[j][j]))) {
    report(INFO,"EvtGen") << "Failing 3"<<endl;
    return 0;
      }
    }
  }
 
  return 1;
}

Referenced by EvtDecayAmp::makeDecay().

Get()

const EvtComplex & EvtSpinDensity::Get	(	int	i,
		int	j
	)		const

Definition at line 104 of file EvtSpinDensity.cc.

                                                     {
  assert(i<dim&&j<dim);
  return rho[i][j];
}

Referenced by EvtAmp::contract(), and EvtParticle::setSpinDensityForwardHelicityBasis().

GetDim()

int EvtSpinDensity::GetDim

(

)

const

Definition at line 95 of file EvtSpinDensity.cc.

                                 {
  return dim;
}

Referenced by EvtParticle::setSpinDensityForwardHelicityBasis().

NormalizedProb()

double EvtSpinDensity::NormalizedProb

(

const EvtSpinDensity &

d

)

Definition at line 121 of file EvtSpinDensity.cc.

                                                            {
 
  int i,j;
  EvtComplex prob(0.0,0.0);
  double norm=0.0;
 
  if (dim!=d.dim) {
    report(ERROR,"EvtGen")<<"Not matching dimensions in NormalizedProb"<<endl;
    ::abort();
  }
 
  for(i=0;i<dim;i++){
    norm+=real(rho[i][i]);
    for(j=0;j<dim;j++){
      prob+=rho[i][j]*d.rho[i][j];
    }
  }
 
  if (imag(prob)>0.00000001*real(prob)) {
    report(ERROR,"EvtGen")<<"Imaginary probability:"<<prob<<" "<<norm<<endl;
  }
  if (real(prob)<0.0) {
    report(ERROR,"EvtGen")<<"Negative probability:"<<prob<<" "<<norm<<endl;
  }
 
  return real(prob)/norm;
 
}

Referenced by EvtbTosllAmp::CalcMaxProb(), EvtSemiLeptonicAmp::CalcMaxProb(), and EvtDecayAmp::makeDecay().

operator=()

EvtSpinDensity & EvtSpinDensity::operator=

(

const EvtSpinDensity &

density

)

Definition at line 47 of file EvtSpinDensity.cc.

                                                                      {
  int i,j;
  SetDim(density.dim);
 
  for(i=0;i<dim;i++){
    for(j=0;j<dim;j++){
      rho[i][j]=density.rho[i][j];
    }
  }
 
  return *this;
 
}

Set()

void EvtSpinDensity::Set	(	int	i,
		int	j,
		const EvtComplex &	rhoij
	)

Definition at line 99 of file EvtSpinDensity.cc.

                                                           {
  assert(i<dim&&j<dim);
  rho[i][j]=rhoij;
}

Referenced by EvtAmp::contract(), EvtVectorIsr::decay(), EvtAmp::getBackwardSpinDensity(), EvtAmp::getForwardSpinDensity(), EvtAmp::getSpinDensity(), EvtDecayAmp::makeDecay(), EvtParticle::setPolarizedSpinDensity(), EvtParticle::setSpinDensityForwardHelicityBasis(), and EvtParticle::setVectorSpinDensity().

SetDiag()

void EvtSpinDensity::SetDiag

(

int

n

)

Definition at line 109 of file EvtSpinDensity.cc.

                                 {
  SetDim(n);
  int i,j;
 
  for(i=0;i<n;i++){
    for(j=0;j<n;j++){
      rho[i][j]=EvtComplex(0.0);
    }
    rho[i][i]=EvtComplex(1.0);
  }
}

Referenced by EvtbTosllAmp::CalcMaxProb(), EvtSemiLeptonicAmp::CalcMaxProb(), EvtParticle::decay(), EvtDecayIncoherent::makeDecay(), EvtDecayProb::makeDecay(), EvtParticleFactory::particleFactory(), EvtParticle::setDiagonalSpinDensity(), EvtParticle::setPolarizedSpinDensity(), and EvtParticle::setVectorSpinDensity().

SetDim()

void EvtSpinDensity::SetDim

(

int

n

)

Definition at line 75 of file EvtSpinDensity.cc.

                                {
  if (dim==n) return;
  if (dim!=0){
    int i;
    for(i=0;i<dim;i++) delete [] rho[i];
    delete [] rho;
    rho=0;
    dim=0;
  }
  if (n==0) return;
  dim=n;
  rho=new EvtComplexPtr[n];
  int i;
  for(i=0;i<n;i++){
    rho[i]=new EvtComplex[n];
  }
 
 
}

Referenced by EvtAmp::contract(), EvtVectorIsr::decay(), EvtSpinDensity(), EvtAmp::getBackwardSpinDensity(), EvtAmp::getForwardSpinDensity(), EvtAmp::getSpinDensity(), EvtDecayAmp::makeDecay(), operator=(), SetDiag(), and EvtParticle::setSpinDensityForwardHelicityBasis().

Friends And Related Function Documentation

operator<<

std::ostream & operator<< ( std::ostream &  s, const EvtSpinDensity &  d  )

friend

---

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

• EvtSpinDensity.hh
• EvtSpinDensity.cc