beambeamsystem.h

Go to the documentation of this file. Or view the newest version in sPHENIX GitHub for file 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:: 213                         $: revision of last commit
// $Author:: butter                   $: author of last commit
// $Date:: 2015-08-15 22:08:02 +0100 #$: 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"
#include <vector>
#include "inputParameters.h"

class beamBeamSystem
{

public:

  beamBeamSystem(const inputParameters& input,
           const beam&            electronBeam,
                 const beam&            targetBeam);
  beamBeamSystem(const inputParameters& input);
  ~beamBeamSystem();

  const beam& electronBeam() const { return _electronBeam; }  
  const beam& targetBeam() const { return _targetBeam; }  

  double probabilityOfBreakup(const double D) const;
  
  double cmsBoost() const { return _cmsBoost; }
  
  double beamLorentzGamma() const { return _beamLorentzGamma; }
  
  void init();

private:
  void generateBreakupProbabilities();
  double probabilityOfHadronBreakup(const double impactparameter);
  double probabilityOfPhotonBreakup(const double impactparameter, const int mode);

  double _pHadronBreakup;
  double _pPhotonBreakup;

  double _beamLorentzGamma;  
        const double _electronBeamLorentzGamma;  
        const double _targetBeamLorentzGamma;  
  const int    _beamBreakupMode;   






                 
  beam   _electronBeam;             
  beam   _targetBeam;             

  double _cmsBoost;    
  
  std::vector<double> _breakupProbabilities; 
//  double _breakupImpactParameterStep; ///< Step size in the calculation of the breakup probs
//  double _breakupCutOff;  ///< Cut off for minimum impact parameter probability
};


#endif  // BEAMBEAMSYSTEM_H