d5/d88/eventchannel_8h_source.html

//

//    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:: 255                         $: revision of last commit

// $Author:: jnystrand                $: author of last commit

// $Date:: 2016-04-06 14:01:46 +0100 #$: date of last commit

//

// Description:

//

//

//


#ifndef EVENTCHANNEL_H

#define EVENTCHANNEL_H


#include <vector>


#include "starlightconstants.h"

#include "readinluminosity.h"

#include "beambeamsystem.h"

#include "randomgenerator.h"

#include "eXevent.h"

#include "inputParameters.h"


class eventChannel : public readLuminosity

{


public:


  eventChannel(const inputParameters& input, beamBeamSystem&        bbsystem);

  virtual ~eventChannel();


  unsigned long nmbAttempts() const {return _nmbAttempts;}

  unsigned long nmbAccepted() const {return _nmbAccepted;}


  virtual eXEvent e_produceEvent() = 0;


  static void transform(const double betax,

                        const double betay,

                        const double betaz,

                        double&      E,

                        double&      px,

                        double&      py,

                        double&      pz,

                        int&         iFbadevent);


  beamBeamSystem  _bbs;

  randomGenerator _randy;

  static double pseudoRapidity(const double px,

                               const double py,

                               const double pz);


  double getTotalChannelCrossSection () const {return _totalChannelCrossSection;}

  void   setTotalChannelCrossSection (double sigma) {_totalChannelCrossSection = sigma;}


protected:


  unsigned long _nmbAttempts;

  unsigned long _nmbAccepted;


  bool   _ptCutEnabled;

  double _ptCutMin;

  double _ptCutMax;

  bool   _etaCutEnabled;

  double _etaCutMin;

  double _etaCutMax;


  double _totalChannelCrossSection;


};


#endif  // EVENTCHANNEL_H