include/EvtGenBase/EvtSpinDensity.hh

0001
0002 /***********************************************************************
0003 * Copyright 1998-2020 CERN for the benefit of the EvtGen authors       *
0004 *                                                                      *
0005 * This file is part of EvtGen.                                         *
0006 *                                                                      *
0007 * EvtGen is free software: you can redistribute it and/or modify       *
0008 * it under the terms of the GNU General Public License as published by *
0009 * the Free Software Foundation, either version 3 of the License, or    *
0010 * (at your option) any later version.                                  *
0011 *                                                                      *
0012 * EvtGen is distributed in the hope that it will be useful,            *
0013 * but WITHOUT ANY WARRANTY; without even the implied warranty of       *
0014 * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the        *
0015 * GNU General Public License for more details.                         *
0016 *                                                                      *
0017 * You should have received a copy of the GNU General Public License    *
0018 * along with EvtGen.  If not, see <https://www.gnu.org/licenses/>.     *
0019 ***********************************************************************/
0020
0021 #ifndef EVTSPINDENSITY_HH
0022 #define EVTSPINDENSITY_HH
0023 #include "EvtGenBase/EvtComplex.hh"
0024
0025 // Description: This class holds a spin density matrix, it is
0026 //              a complex nxn matrix.
0027
0028 class EvtSpinDensity {
0029   public:
0030     EvtSpinDensity( const EvtSpinDensity& density );
0031     EvtSpinDensity& operator=( const EvtSpinDensity& density );
0032     virtual ~EvtSpinDensity();
0033
0034     EvtSpinDensity();
0035     void setDim( int n );
0036     int getDim() const;
0037     void set( int i, int j, const EvtComplex& rhoij );
0038     const EvtComplex& get( int i, int j ) const;
0039     double normalizedProb( const EvtSpinDensity& d );
0040     friend std::ostream& operator<<( std::ostream& s, const EvtSpinDensity& d );
0041     void setDiag( int n );
0042
0043     int check();
0044
0045   private:
0046     EvtComplexPtrPtr rho;
0047     int dim;
0048 };
0049
0050 #endif