photonNucleusCrossSection.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 PHOTONNUCLEUSCROSSSECTION_H
#define PHOTONNUCLEUSCROSSSECTION_H


#include "starlightconstants.h"
#include "beambeamsystem.h"


class photonNucleusCrossSection {

public:

        photonNucleusCrossSection(const inputParameters& input,
                                  const beamBeamSystem&  bbsystem);
        ~photonNucleusCrossSection();
  
        double         slopeParameter  () const { return _slopeParameter;   }  
        double         getChannelMass  () const { return _channelMass;      }  
        double         getBNORM        () const { return _BNORM;            }
        double         luminosity      () const { return _luminosity;       }  
        beamBeamSystem getbbs          () const { return _bbs;              }  
        double         vmPhotonCoupling() const { return _vmPhotonCoupling; }  
        double         getDefaultC     () const { return _defaultC;         }
        double         maxPhotonEnergy () const { return _maxPhotonEnergy;  }  

        void crossSectionCalculation(const double bwnormsave);
        // Will think about it...For VMs we just calculate it
        // So just use the wide or narrow constructor to calculate it
        // wide/narrow will inherit this.
        double getcsgA(const double Egamma,
                       const double W);
        double photonFlux(const double Egamma);
        double sigmagp(const double Wgp);
        double sigma_A(const double sig_N);
        double breitWigner(const double W,
                           const double C);

private:

        double nepoint(const double Egamma,
                       const double bmin);
  
        beamBeamSystem _bbs;
  
        // copied from inputParameters
        double                               _protonEnergy;
        double                               _beamLorentzGamma;    
        starlightConstants::particleTypeEnum _particleType;
        int                                  _beamBreakupMode;     
        bool                                 _coherentProduction;  
        double                               _incoherentFactor;    
        int                                  _sigmaNucleus;

        // locally defined parameters
        double _slopeParameter;    
        double _vmPhotonCoupling;  
        double _ANORM;
        double _BNORM;
        double _defaultC;
        double _luminosity;       
        double _maxPhotonEnergy;  
        double _width;            
        double _channelMass;      
};


#endif  // PHOTONNUCLEUSCROSSSECTION_H