EvtAmp.hh

Go to the documentation of this file.

//--------------------------------------------------------------------------
//
// Environment:
//      This software is part of the EvtGen package developed jointly
//      for the BaBar and CLEO collaborations.  If you use all or part
//      of it, please give an appropriate acknowledgement.
//
// Copyright Information: See EvtGen/COPYRIGHT
//      Copyright (C) 1998      Caltech, UCSB
//
// Module: EvtGen/EvtAmp.hh
//
// Description:This is the (new) amplitude class.
//
// Modification history:
//
//    DJL/RYD     May 29, 1997         Module created
//
//------------------------------------------------------------------------
 
#ifndef EVTAMP_HH
#define EVTAMP_HH
 
#include "EvtComplex.hh"
#include "EvtPatches.hh"
 
class EvtSpinDensity;
class EvtId;
 
class EvtAmp {
 
public:
  EvtAmp();
  virtual ~EvtAmp() {}
 
  EvtAmp( const EvtAmp& amp );
 
  void init( EvtId p, int ndaug, EvtId* daug );
 
  void setAmp( int* ind, const EvtComplex& amp );
 
  const EvtComplex& getAmp( int* ind ) const;
 
  EvtSpinDensity getSpinDensity();
 
  EvtSpinDensity contract( int i, const EvtAmp& a );
  EvtAmp         contract( int i, const EvtSpinDensity& rho );
 
  // sum over the i:th daugther of a1 and contract with parent of a2
  EvtAmp contract( int i, const EvtAmp& a1, const EvtAmp& a2 );
 
  EvtSpinDensity getForwardSpinDensity( EvtSpinDensity* rho_list, int k );
  EvtSpinDensity getBackwardSpinDensity( EvtSpinDensity* rho_list );
 
  EvtAmp& operator=( const EvtAmp& amp );
 
  /**
   * sets the amplitudes calculated in the decay objects
   */
  void vertex( const EvtComplex& amp );
 
  /**
   * sets the amplitudes calculated in the decay objects
   */
  void vertex( int i1, const EvtComplex& amp );
 
  /**
   * sets the amplitudes calculated in the decay objects
   */
  void vertex( int i1, int i2, const EvtComplex& amp );
 
  /**
   * sets the amplitudes calculated in the decay objects
   */
  void vertex( int i1, int i2, int i3, const EvtComplex& amp );
 
  /**
   * sets the amplitudes calculated in the decay objects
   */
  void vertex( int* i1, const EvtComplex& amp );
 
  void dump();
 
private:
  friend class EvtDecayAmp;
  friend class EvtKstarstargamma;
 
  void setNDaug( int n );
  void setNState( int parent_states, int* daug_states );
 
  // the amplitudes
  EvtComplex _amp[450]; // for run RaritaShwinger particle, the lowest dimension need 144, or
                        // crashed.  Deng z r and Ping rg
 
  // the number of daughters
  int _ndaug;
 
  // the number of states of the parent
  int _pstates;
 
  // number of states of the daughter
  int dstates[10];
 
  // the nontrivial index of the daughter
  int _dnontrivial[10];
 
  // number of nontrivial daugts+parent
  int _nontrivial;
 
  // compact nstates
  int _nstate[5];
};
 
#endif