BaseProtonTransport.h

Go to the documentation of this file.

#ifndef BASEPROTONTRANSPORT
#define BASEPROTONTRANSPORT
#include "FWCore/Framework/interface/EventSetup.h"
#include "FWCore/ParameterSet/interface/ParameterSet.h"
#include "HepMC/GenEvent.h"
#include "FWCore/MessageLogger/interface/MessageLogger.h"
#include "Utilities/PPS/interface/PPSUnitConversion.h"
#include "SimDataFormats/Forward/interface/LHCTransportLink.h"
#include <vector>
#include <map>
#include <string>
#include "TLorentzVector.h"

namespace CLHEP {
  class HepRandomEngine;
}
class BaseProtonTransport {
public:
  BaseProtonTransport(const edm::ParameterSet& iConfig);
  virtual ~BaseProtonTransport();

  std::vector<LHCTransportLink>& getCorrespondenceMap() { return m_CorrespondenceMap; }
  virtual void process(const HepMC::GenEvent* ev, const edm::EventSetup& es, CLHEP::HepRandomEngine* engine) = 0;

  void clear();

  void addPartToHepMC(const HepMC::GenEvent*, HepMC::GenEvent*);

  void ApplyBeamCorrection(HepMC::GenParticle* p);
  void ApplyBeamCorrection(TLorentzVector& p);

  double beamEnergy() { return beamEnergy_; };
  double beamMomentum() { return beamMomentum_; };

protected:
  void setBeamFileNames(const std::string& nam1, const std::string& nam2) {
    beam1Filename_ = nam1;
    beam2Filename_ = nam2;
  };

  void setBeamParameters(double stx, double sty, double sx, double sy, double se) {
    m_sigmaSTX = stx;
    m_sigmaSTY = sty;
    m_sigmaSX = sx;
    m_sigmaSY = sy;
    m_sig_E = se;
  };

  void setCrossingAngles(double cx45, double cx56, double cy45, double cy56) {
    fCrossingAngleX_45_ = cx45;
    fCrossingAngleX_56_ = cx56;
    fCrossingAngleY_45_ = cy45;
    fCrossingAngleY_56_ = cy56;
  };

  enum class TransportMode { HECTOR, TOTEM, OPTICALFUNCTIONS };
  TransportMode MODE;

  CLHEP::HepRandomEngine* engine_{nullptr};
  std::vector<LHCTransportLink> m_CorrespondenceMap;
  std::map<unsigned int, TLorentzVector> m_beamPart;
  std::map<unsigned int, double> m_xAtTrPoint;
  std::map<unsigned int, double> m_yAtTrPoint;

  bool verbosity_;
  bool bApplyZShift_;
  bool useBeamPositionFromLHCInfo_;
  bool produceHitsRelativeToBeam_;

  std::string beam1Filename_{""};
  std::string beam2Filename_{""};

  double fPPSRegionStart_45_;
  double fPPSRegionStart_56_;
  double fCrossingAngleX_45_{0.0};
  double fCrossingAngleX_56_{0.0};
  double fCrossingAngleY_45_{0.0};
  double fCrossingAngleY_56_{0.0};

  double etaCut_;
  double momentumCut_;

private:
  double beamMomentum_;
  double beamEnergy_;

  double m_sigmaSTX{0.0};
  double m_sigmaSTY{0.0};
  double m_sigmaSX{0.0};
  double m_sigmaSY{0.0};
  double m_sig_E{0.0};
};
#endif