#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 26 of file EvtSpinDensity.hh.

Constructor & Destructor Documentation

◆ EvtSpinDensity() [1/2]

EvtSpinDensity::EvtSpinDensity ( const EvtSpinDensity & density )

Definition at line 32 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]; }
  }
}

Referenced by EvtSpinDensity(), NormalizedProb(), operator<<, and operator=().

◆ ~EvtSpinDensity()

EvtSpinDensity::~EvtSpinDensity ( )

virtual

Definition at line 57 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 67 of file EvtSpinDensity.cc.

                               {
  dim = 0;
  rho = 0;
}

Member Function Documentation

◆ Check()

int EvtSpinDensity::Check ( )

Definition at line 138 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 96 of file EvtSpinDensity.cc.

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

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

◆ GetDim()

int EvtSpinDensity::GetDim ( ) const

Definition at line 89 of file EvtSpinDensity.cc.

89{ return dim; }

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

◆ NormalizedProb()

double EvtSpinDensity::NormalizedProb ( const EvtSpinDensity & d )

Definition at line 112 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 45 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 91 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(), and EvtParticle::setVectorSpinDensity().

◆ SetDiag()

void EvtSpinDensity::SetDiag ( int n )

Definition at line 101 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::setPolarizedSpinDensity(), and EvtParticle::setVectorSpinDensity().

◆ SetDim()

void EvtSpinDensity::SetDim ( int n )

Definition at line 72 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=(), and SetDiag().

◆ operator<<

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

friend

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

Public Member Functions

Friends