RandomtXiGunProducer.h

Go to the documentation of this file.

#ifndef RandomtXiGunProducer_H
#define RandomtXiGunProducer_H

#include "IOMC/ParticleGuns/interface/BaseRandomtXiGunProducer.h"
#include "SimDataFormats/GeneratorProducts/interface/HepMCProduct.h"

namespace edm {

  class RandomtXiGunProducer : public BaseRandomtXiGunProducer {
  public:
    RandomtXiGunProducer(const ParameterSet&);
    ~RandomtXiGunProducer() override;

  private:
    void produce(Event& e, const EventSetup& es) override;

    HepMC::FourVector make_particle(double t, double Xi, double phi, int PartID, int direction);
    double Minimum_t(double xi) {
      double partE = fpEnergy * (1. - xi);
      double massSQ = pow(PData->mass().value(), 2);
      double partP = sqrt(partE * partE - massSQ);
      return -2. * (sqrt(fpEnergy * fpEnergy - massSQ) * partP - fpEnergy * partE + massSQ);
    };

  protected:
    // data members

    double fMint;
    double fMaxt;
    double fMinXi;
    double fMaxXi;
  };
}  // namespace edm
#endif