beambeamsystem.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:31 $: date of last commit
//
// Description:
//     this class covers a coliding beam system SK
//
//
//


#ifndef BEAMBEAMSYSTEM_H
#define BEAMBEAMSYSTEM_H


#include "nucleus.h"
#include "beam.h"


class beamBeamSystem
{

public:

        //Better way to do this? Memory issues creating all of theses Beams?
        beamBeamSystem(const beam&            beam1,
                       const beam&            beam2,
                       const double           luminosity,
                       const inputParameters& parameters);
        beamBeamSystem(const beam&            beam1,
                       const beam&            beam2,
                       const inputParameters& parameters);
        beamBeamSystem(const inputParameters& parameters);
        ~beamBeamSystem();

        const beam& beam1() const { return _beam1; }  
        const beam& beam2() const { return _beam2; }  

        //      double getluminosity();
        double probabilityOfBreakup(const double D);

private:

        //      int _ibreakup;//temporary solution until read in parameters are done
        double probabilityOfHadronBreakup(const double impactparameter);
        double probabilityOfPhotonBreakup(const double impactparameter, const int mode);

        double _pHadronBreakup;
        double _pPhotonBreakup;
        //inputParameters inputbbs;
        //              double luminosity;

        double _beamLorentzGamma;  
        int    _beamBreakupMode;   
        beam   _beam1;             
        beam   _beam2;             

};


#endif  // BEAMBEAMSYSTEM_H