TotemTransport.h

Go to the documentation of this file.

#ifndef TOTEM_TRANSPORT
#define TOTEM_TRANSPORT
#include "SimTransport/PPSProtonTransport/interface/BaseProtonTransport.h"
#include "SimTransport/TotemRPProtonTransportParametrization/interface/LHCOpticsApproximator.h"
#include "FWCore/Framework/interface/Event.h"
#include "FWCore/Framework/interface/EventSetup.h"
#include "FWCore/Framework/interface/ESHandle.h"
#include "FWCore/MessageLogger/interface/MessageLogger.h"
 
#include "FWCore/ParameterSet/interface/ParameterSet.h"
#include "FWCore/ParameterSet/interface/ConfigurationDescriptions.h"
#include "FWCore/ParameterSet/interface/ParameterSetDescription.h"
#include "SimDataFormats/GeneratorProducts/interface/HepMCProduct.h"
#include "SimGeneral/HepPDTRecord/interface/ParticleDataTable.h"
 
#include <unordered_map>
#include <array>
 
namespace CLHEP {
  class HepRandomEngine;
}
 
class TotemTransport : public BaseProtonTransport {
public:
  TotemTransport(const edm::ParameterSet& ps);
  ~TotemTransport() override{};
  // look for scattered protons, propagates them, add them to the event
  void process(const HepMC::GenEvent* ev, const edm::EventSetup& es, CLHEP::HepRandomEngine* engine) override;
 
private:
  bool transportProton(const HepMC::GenParticle*);
  LHCOpticsApproximator* ReadParameterization(const std::string&, const std::string&);
 
  LHCOpticsApproximator* m_aprox_ip_150_r = nullptr;
  LHCOpticsApproximator* m_aprox_ip_150_l = nullptr;
 
  std::string m_model_ip_150_r_name;
  std::string m_model_ip_150_l_name;
 
  double m_beampipe_aperture_radius;
  double m_fEtacut;
  double m_fMomentumMin;
};
#endif