BetafuncEvtVtxGenerator.h

Go to the documentation of this file.

#ifndef IOMC_BetafuncEvtVtxGenerator_H
#define IOMC_BetafuncEvtVtxGenerator_H

// $Id: BetafuncEvtVtxGenerator.h,v 1.9 2013/02/27 18:41:06 wmtan Exp $
/*
________________________________________________________________________

 BetafuncEvtVtxGenerator

 Smear vertex according to the Beta function on the transverse plane
 and a Gaussian on the z axis. It allows the beam to have a crossing
 angle (dx/dz and dy/dz).

 Based on GaussEvtVtxGenerator.h
 implemented by Francisco Yumiceva (yumiceva@fnal.gov)

 FERMILAB
 2006
________________________________________________________________________
*/

#include "IOMC/EventVertexGenerators/interface/BaseEvtVtxGenerator.h"
#include "FWCore/Framework/interface/ESWatcher.h"
#include "CondFormats/DataRecord/interface/SimBeamSpotObjectsRcd.h"

namespace CLHEP {
   class RandGaussQ;
}

class BetafuncEvtVtxGenerator : public BaseEvtVtxGenerator 
{
public:
  BetafuncEvtVtxGenerator(const edm::ParameterSet & p);
  virtual ~BetafuncEvtVtxGenerator();

  virtual void beginRun(const edm::Run & , const edm::EventSetup&) override;
  virtual void beginLuminosityBlock(edm::LuminosityBlock const&, edm::EventSetup const&) override;

  //virtual CLHEP::Hep3Vector * newVertex();
  virtual HepMC::FourVector* newVertex() ;

  virtual TMatrixD* GetInvLorentzBoost();

    
  void sigmaZ(double s=1.0);

  void X0(double m=0) { fX0=m; }
  void Y0(double m=0) { fY0=m; }
  void Z0(double m=0) { fZ0=m; }

  void Phi(double m=0) { phi_=m; }
  void Alpha(double m=0) { alpha_=m; }

  void betastar(double m=0) { fbetastar=m; }
  void emittance(double m=0) { femittance=m; }

  double BetaFunction(double z, double z0);
    
private:
  BetafuncEvtVtxGenerator(const BetafuncEvtVtxGenerator &p);
  BetafuncEvtVtxGenerator&  operator = (const BetafuncEvtVtxGenerator & rhs );
  
private:

  bool readDB_;

  double alpha_, phi_;
  //TMatrixD boost_;
  
  double fX0, fY0, fZ0;
  double fSigmaZ;
  //double fdxdz, fdydz;
  double fbetastar, femittance;
  //  double falpha;
  double fTimeOffset;
    
  CLHEP::RandGaussQ*  fRandom ;
  
  void update(const edm::EventSetup& iEventSetup);
  edm::ESWatcher<SimBeamSpotObjectsRcd> parameterWatcher_;
};

#endif