nucleus.h

//
//    Copyright 2010
//
//    This file is part of starlight.
//
//    starlight is free software: you can redistribute it and/or modify
//    it under the terms of the GNU General Public License as published by
//    the Free Software Foundation, either version 3 of the License, or
//    (at your option) any later version.
//
//    starlight is distributed in the hope that it will be useful,
//    but WITHOUT ANY WARRANTY; without even the implied warranty of
//    MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
//    GNU General Public License for more details.
//
//    You should have received a copy of the GNU General Public License
//    along with starlight. If not, see <http://www.gnu.org/licenses/>.
//
//
// File and Version Information:
// $Rev::                             $: revision of last commit
// $Author: jwebb $: author of last commit
// $Date: 2012/11/27 22:27:32 $: date of last commit
//
// Description:
//
//
//


#ifndef NUCLEUS_H
#define NUCLEUS_H


#include <cmath>


//This class holds the information for a target nucleus
class nucleus
{

public:

        nucleus(const int    Z,
                const int    A,
                const double deuteronSlopePar,
                const bool   dAuCoherentProduction);
        ~nucleus();
 
        int    Z              () const { return _Z;                     }  
        int    A              () const { return _A;                     }  
        double woodSaxonRadius() const { return 1.2 * pow(_A, 1. / 3.); }  
        double nuclearRadius  () const; 

        double formFactor(const double t) const;  
        double thickness (const double b) const;  

        double Q0  () const { return _Q0;   }
        double rho0() const { return _rho0; }

        
private:

        double woodSaxonSkinDepth() const { return 0.53;                   }  
        double fritiofR0         () const { return _r0 * pow(_A, (1./3.)); }  

        double rws(const double r) const
        { return 1.0 / (1. + exp((r - fritiofR0()) / woodSaxonSkinDepth())); } 

        int    _Z;                      
        int    _A;                      
        double _deuteronSlopePar;       
        bool   _dAuCoherentProduction;  

        double _r0;
        double _rho0;
        double _Q0;

};


#endif  // NUCLEUS_H