HLLHCEvtVtxGenerator.h

Go to the documentation of this file.

#ifndef IOMC_HLLHCEvtVtxGenerator_H
#define IOMC_HLLHCEvtVtxGenerator_H
 
#include "IOMC/EventVertexGenerators/interface/BaseEvtVtxGenerator.h"
 
#include <string>
 
namespace CLHEP {
  class RandFlat;
}
 
namespace edm {
  class ConfigurationDescriptions;
}
 
class HLLHCEvtVtxGenerator : public BaseEvtVtxGenerator {
public:
  HLLHCEvtVtxGenerator(const edm::ParameterSet& p);
 
  ~HLLHCEvtVtxGenerator() override;
 
  static void fillDescriptions(edm::ConfigurationDescriptions& descriptions);
 
  HepMC::FourVector newVertex(CLHEP::HepRandomEngine*) const override;
 
  TMatrixD const* GetInvLorentzBoost() const override { return nullptr; };
 
private:
  HLLHCEvtVtxGenerator(const HLLHCEvtVtxGenerator& p) = delete;
 
  HLLHCEvtVtxGenerator& operator=(const HLLHCEvtVtxGenerator& rhs) = delete;
 
  //spatial and time offset for mean collision
  const double fMeanX, fMeanY, fMeanZ, fTimeOffset;
 
  //proton beam energy
  const double momeV;
  const double gamma;
  const double beta;
  const double betagamma;
 
  //crossing angle
  const double phi;
 
  //crab cavity frequency
  const double wcc;
 
  // 800 MHz RF?
  const bool RF800;
 
  //beta crossing plane (m)
  const double betx;
 
  //beta separation plane (m)
  const double bets;
 
  //horizontal emittance
  const double epsxn;
 
  //vertical emittance
  const double epssn;
 
  //bunch length
  const double sigs;
 
  //crabbing angle crossing
  const double alphax;
 
  //crabbing angle separation
  const double alphay;
 
  // ratio of crabbing angle to crossing angle
  const double oncc;
 
  //normalized crossing emittance
  const double epsx;
 
  //normlaized separation emittance
  const double epss;
 
  //size in x
  const double sigx;
 
  // crossing angle * crab frequency
  const double phiCR;
 
  //width for y plane
  double sigma(double z, double epsilon, double beta, double betagamma) const;
 
  //density with crabbing
  double integrandCC(double x, double z, double t) const;
 
  // 4D intensity
  double intensity(double x, double y, double z, double t) const;
};
 
#endif