BeamProfileVtxGenerator Class Reference

#include <BeamProfileVtxGenerator.h>

Inheritance diagram for BeamProfileVtxGenerator:

Public Member Functions
void	beamPos (double m=0)
	set mean in Z in cm More...
	BeamProfileVtxGenerator (const edm::ParameterSet &p)
	BeamProfileVtxGenerator (const BeamProfileVtxGenerator &p)=delete
void	eta (double m=0)
	set eta More...
TMatrixD const *	GetInvLorentzBoost () const override
void	meanX (double m=0)
	set mean in X in cm More...
void	meanY (double m=0)
	set mean in Y in cm More...
HepMC::FourVector	newVertex (CLHEP::HepRandomEngine *) const override
	return a new event vertex More...
BeamProfileVtxGenerator &	operator= (const BeamProfileVtxGenerator &rhs)=delete
void	phi (double m=0)
	set phi in radian More...
void	psi (double m=999)
	set psi in radian More...
void	setType (bool m=true)
	set type More...
void	sigmaX (double s=1.0)
	set resolution in X in cm More...
void	sigmaY (double s=1.0)
	set resolution in Y in cm More...
	~BeamProfileVtxGenerator () override
Public Member Functions inherited from BaseEvtVtxGenerator
	BaseEvtVtxGenerator (const edm::ParameterSet &)
void	produce (edm::Event &, const edm::EventSetup &) override
	~BaseEvtVtxGenerator () override
Public Member Functions inherited from edm::stream::EDProducer<>
	EDProducer ()=default
	EDProducer (const EDProducer &)=delete
bool	hasAbilityToProduceInBeginLumis () const final
bool	hasAbilityToProduceInBeginProcessBlocks () const final
bool	hasAbilityToProduceInBeginRuns () const final
bool	hasAbilityToProduceInEndLumis () const final
bool	hasAbilityToProduceInEndProcessBlocks () const final
bool	hasAbilityToProduceInEndRuns () const final
const EDProducer &	operator= (const EDProducer &)=delete

Private Attributes
std::vector< double >	fdistn
double	fEta
bool	ffile
double	fMeanX
double	fMeanY
double	fMeanZ
double	fPhi
double	fPsi
double	fSigmaX
double	fSigmaY
double	fTheta
double	fTimeOffset
bool	fType
int	nBinx
int	nBiny

Additional Inherited Members
Public Types inherited from edm::stream::EDProducer<>
using	CacheTypes = CacheContexts< T... >
using	GlobalCache = typename CacheTypes::GlobalCache
using	HasAbility = AbilityChecker< T... >
using	InputProcessBlockCache = typename CacheTypes::InputProcessBlockCache
using	LuminosityBlockCache = typename CacheTypes::LuminosityBlockCache
using	LuminosityBlockContext = LuminosityBlockContextT< LuminosityBlockCache, RunCache, GlobalCache >
using	LuminosityBlockSummaryCache = typename CacheTypes::LuminosityBlockSummaryCache
using	RunCache = typename CacheTypes::RunCache
using	RunContext = RunContextT< RunCache, GlobalCache >
using	RunSummaryCache = typename CacheTypes::RunSummaryCache

Detailed Description

Definition at line 19 of file BeamProfileVtxGenerator.h.

Constructor & Destructor Documentation

BeamProfileVtxGenerator() [1/2]

BeamProfileVtxGenerator::BeamProfileVtxGenerator

(

const edm::ParameterSet &

p

)

Definition at line 19 of file BeamProfileVtxGenerator.cc.

References beamPos(), eta(), fdistn, fEta, ffile, geometryDiff::file, fMeanX, fMeanY, fMeanZ, fPhi, fPsi, fSigmaX, fSigmaY, fTheta, fTimeOffset, fType, mps_fire::i, recoMuon::in, dqmdumpme::last, meanX(), meanY(), nBinx, nBiny, AlCaHLTBitMon_ParallelJobs::p, phi(), psi(), setType(), sigmaX(), sigmaY(), and AlCaHLTBitMon_QueryRunRegistry::string.

                                                                         : BaseEvtVtxGenerator(p) {
  meanX(p.getParameter<double>("BeamMeanX") * cm);
  meanY(p.getParameter<double>("BeamMeanY") * cm);
  beamPos(p.getParameter<double>("BeamPosition") * cm);
  sigmaX(p.getParameter<double>("BeamSigmaX") * cm);
  sigmaY(p.getParameter<double>("BeamSigmaY") * cm);
  double fMinEta = p.getParameter<double>("MinEta");
  double fMaxEta = p.getParameter<double>("MaxEta");
  double fMinPhi = p.getParameter<double>("MinPhi");
  double fMaxPhi = p.getParameter<double>("MaxPhi");
  eta(0.5 * (fMinEta + fMaxEta));
  phi(0.5 * (fMinPhi + fMaxPhi));
  psi(p.getParameter<double>("Psi"));
  nBinx = p.getParameter<int>("BinX");
  nBiny = p.getParameter<int>("BinY");
  ffile = p.getParameter<bool>("UseFile");
  fTimeOffset = p.getParameter<double>("TimeOffset") * ns * c_light;

  if (ffile) {
    std::string file = p.getParameter<std::string>("File");
    std::ifstream is(file.c_str(), std::ios::in);
    if (is) {
      double elem, sum = 0;
      while (!is.eof()) {
        is >> elem;
        fdistn.push_back(elem);
        sum += elem;
      }
      if (((int)(fdistn.size())) >= nBinx * nBiny) {
        double last = 0;
        for (unsigned int i = 0; i < fdistn.size(); i++) {
          fdistn[i] /= sum;
          fdistn[i] += last;
          last = fdistn[i];
        }
        setType(false);
      } else {
        ffile = false;
      }
    } else {
      ffile = false;
    }
  }
  if (!ffile) {
    setType(p.getParameter<bool>("GaussianProfile"));
  }

  edm::LogInfo("VertexGenerator") << "BeamProfileVtxGenerator: with "
                                  << "beam along eta = " << fEta << " (Theta = " << fTheta / deg
                                  << ") phi = " << fPhi / deg << ") psi = " << fPsi / deg << " centred at (" << fMeanX
                                  << ", " << fMeanY << ", " << fMeanZ << ") "
                                  << "and spread (" << fSigmaX << ", " << fSigmaY << ") of type Gaussian = " << fType
                                  << " use file " << ffile;
  if (ffile)
    edm::LogInfo("VertexGenerator") << "With " << nBinx << " bins "
                                    << " along X and " << nBiny << " bins along Y";
}

BeamProfileVtxGenerator() [2/2]

BeamProfileVtxGenerator::BeamProfileVtxGenerator ( const BeamProfileVtxGenerator &  p )

delete

Copy constructor

~BeamProfileVtxGenerator()

BeamProfileVtxGenerator::~BeamProfileVtxGenerator ( )

override

Definition at line 77 of file BeamProfileVtxGenerator.cc.

{}

Member Function Documentation

beamPos()

void BeamProfileVtxGenerator::beamPos ( double  m = 0 )

inline

set mean in Z in cm

Definition at line 44 of file BeamProfileVtxGenerator.h.

References fMeanZ, and visualization-live-secondInstance_cfg::m.

Referenced by BeamProfileVtxGenerator().

{ fMeanZ = m; }

eta()

void BeamProfileVtxGenerator::eta

(

double

m = 0

)

set eta

Definition at line 179 of file BeamProfileVtxGenerator.cc.

References JetChargeProducer_cfi::exp, fEta, fTheta, and alignCSCRings::s.

Referenced by Particle.Particle::__str__(), BeamProfileVtxGenerator(), datamodel.Object::DeltaR(), Jet.Jet::jetID(), datamodel.Object::p4(), and Jet.Jet::puJetId().

                                          {
  fEta = s;
  fTheta = 2.0 * atan(exp(-fEta));
}

GetInvLorentzBoost()

TMatrixD const* BeamProfileVtxGenerator::GetInvLorentzBoost ( ) const

inlineoverridevirtual

This method - and the comment - is a left-over from COBRA-OSCAR time : return the last generated event vertex. If no vertex has been generated yet, a NULL pointer is returned.

Implements BaseEvtVtxGenerator.

Definition at line 32 of file BeamProfileVtxGenerator.h.

{ return nullptr; }

meanX()

void BeamProfileVtxGenerator::meanX ( double  m = 0 )

inline

set mean in X in cm

Definition at line 40 of file BeamProfileVtxGenerator.h.

References fMeanX, and visualization-live-secondInstance_cfg::m.

Referenced by BeamProfileVtxGenerator().

{ fMeanX = m; }

meanY()

void BeamProfileVtxGenerator::meanY ( double  m = 0 )

inline

set mean in Y in cm

Definition at line 42 of file BeamProfileVtxGenerator.h.

References fMeanY, and visualization-live-secondInstance_cfg::m.

Referenced by BeamProfileVtxGenerator().

{ fMeanY = m; }

newVertex()

HepMC::FourVector BeamProfileVtxGenerator::newVertex ( CLHEP::HepRandomEngine *  engine ) const

overridevirtual

return a new event vertex

Implements BaseEvtVtxGenerator.

Definition at line 80 of file BeamProfileVtxGenerator.cc.

References funct::cos(), fdistn, ffile, fMeanX, fMeanY, fMeanZ, fPhi, fPsi, fSigmaX, fSigmaY, fTheta, fTimeOffset, fType, mps_fire::i, ALPAKA_ACCELERATOR_NAMESPACE::ecal::reconstruction::internal::endcap::ix(), ALPAKA_ACCELERATOR_NAMESPACE::ecal::reconstruction::internal::endcap::iy(), LogDebug, M_PI, nBinx, nBiny, MetAnalyzer::pv(), diffTwoXMLs::r1, and funct::sin().

                                                                                     {
  double aX, aY;
  if (ffile) {
    double r1 = engine->flat();
    int ixy = 0, ix, iy;
    for (unsigned int i = 0; i < fdistn.size(); i++) {
      if (r1 > fdistn[i])
        ixy = i + 1;
    }
    if (ixy >= nBinx * nBiny) {
      ix = nBinx - 1;
      iy = nBiny - 1;
    } else {
      ix = ixy % nBinx;
      iy = (ixy - ix) / nBinx;
    }
    aX = 0.5 * (2 * ix - nBinx + 2 * engine->flat()) * fSigmaX + fMeanX;
    aY = 0.5 * (2 * iy - nBiny + 2 * engine->flat()) * fSigmaY + fMeanY;
  } else {
    if (fType) {
      aX = fSigmaX * CLHEP::RandGaussQ::shoot(engine) + fMeanX;
      aY = fSigmaY * CLHEP::RandGaussQ::shoot(engine) + fMeanY;
    } else {
      aX = CLHEP::RandFlat::shoot(engine, -0.5 * fSigmaX, 0.5 * fSigmaX) + fMeanX;
      aY = CLHEP::RandFlat::shoot(engine, -0.5 * fSigmaY, 0.5 * fSigmaY) + fMeanY;
    }
  }

  double xp, yp, zp;
  if (2. * M_PI < fabs(fPsi)) {
    xp = -aX * cos(fTheta) * cos(fPhi) + aY * sin(fPhi) + fMeanZ * sin(fTheta) * cos(fPhi);
    yp = -aX * cos(fTheta) * sin(fPhi) - aY * cos(fPhi) + fMeanZ * sin(fTheta) * sin(fPhi);
    zp = aX * sin(fTheta) + fMeanZ * cos(fTheta);
  } else  // for endcap testbeam, use 3rd angle psi
  {
    const HepGeom::Vector3D<double> av(aX, aY, fMeanZ);

    /*
     static const double kRadToDeg ( 180./M_PI ) ;
     std::cout<<"theta = "<<std::setw(7) << std::setiosflags( std::ios::fixed ) << std::setprecision(5)
          <<fTheta*kRadToDeg<<", phi="<<std::setw(7) << std::setiosflags( std::ios::fixed ) << std::setprecision(5)
          << fPhi*kRadToDeg <<", PSI="<<std::setw(7) << std::setiosflags( std::ios::fixed ) << std::setprecision(5)
          <<fPsi*kRadToDeg<<std::endl ;
*/
    const HepGeom::RotateZ3D R1(fPhi - M_PI);
    const HepGeom::Point3D<double> xUnit(0, 1, 0);
    const HepGeom::Point3D<double> zUnit(0, 0, 1);
    const HepGeom::Transform3D RXRZ(HepGeom::Rotate3D(-fTheta, R1 * xUnit) * R1);
    const HepGeom::Transform3D TRF(HepGeom::Rotate3D(fPsi, RXRZ * zUnit) * RXRZ);
    /*
     std::cout<<"\n\n$$$$$$$$$$$Transform="
          <<" thetaZ = "<<std::setw(7) << std::setiosflags( std::ios::fixed ) << std::setprecision(5)
          <<TRF.getRotation().thetaZ()*kRadToDeg
          <<", phiZ = "<<std::setw(7) << std::setiosflags( std::ios::fixed ) << std::setprecision(5)
          <<TRF.getRotation().phiZ()*kRadToDeg
          <<", thetaY = "<<std::setw(7) << std::setiosflags( std::ios::fixed ) << std::setprecision(5)
          <<TRF.getRotation().thetaY()*kRadToDeg
          <<", phiY = "<<std::setw(7) << std::setiosflags( std::ios::fixed ) << std::setprecision(5)
          <<TRF.getRotation().phiY()*kRadToDeg
          <<", thetaX = "<<std::setw(7) << std::setiosflags( std::ios::fixed ) << std::setprecision(5)
          <<TRF.getRotation().thetaX()*kRadToDeg
          <<", phiX = "<<std::setw(7) << std::setiosflags( std::ios::fixed ) << std::setprecision(5)
          <<TRF.getRotation().phiX()*kRadToDeg
          <<std::setw(7) << std::setiosflags( std::ios::fixed ) << std::setprecision(5)
          <<TRF.getTranslation()<<std::endl ;
*/
    const HepGeom::Vector3D<double> pv(TRF * av);

    xp = pv.x();
    yp = pv.y();
    zp = pv.z();
  }

  LogDebug("VertexGenerator") << "BeamProfileVtxGenerator: Vertex created "
                              << "at (" << xp << ", " << yp << ", " << zp << ", " << fTimeOffset << ")";
  return HepMC::FourVector(xp, yp, zp, fTimeOffset);
  ;
}

operator=()

BeamProfileVtxGenerator& BeamProfileVtxGenerator::operator= ( const BeamProfileVtxGenerator &  rhs )

delete

Copy assignment operator

phi()

void BeamProfileVtxGenerator::phi ( double  m = 0 )

inline

set phi in radian

Definition at line 49 of file BeamProfileVtxGenerator.h.

References fPhi, and visualization-live-secondInstance_cfg::m.

Referenced by Particle.Particle::__str__(), BeamProfileVtxGenerator(), datamodel.Object::DeltaR(), datamodel.Object::p4(), and ntupleDataFormat.Track::phiPull().

{ fPhi = m; }

psi()

void BeamProfileVtxGenerator::psi ( double  m = 999 )

inline

set psi in radian

Definition at line 51 of file BeamProfileVtxGenerator.h.

References fPsi, and visualization-live-secondInstance_cfg::m.

Referenced by BeamProfileVtxGenerator().

{ fPsi = m; }

setType()

void BeamProfileVtxGenerator::setType

(

bool

m = true

)

set type

Definition at line 184 of file BeamProfileVtxGenerator.cc.

References fType, and alignCSCRings::s.

Referenced by BeamProfileVtxGenerator().

{ fType = s; }

sigmaX()

void BeamProfileVtxGenerator::sigmaX

(

double

s = 1.0

)

set resolution in X in cm

Definition at line 159 of file BeamProfileVtxGenerator.cc.

References fSigmaX, and alignCSCRings::s.

Referenced by BeamProfileVtxGenerator().

                                             {
  if (s >= 0) {
    fSigmaX = s;
  } else {
    edm::LogWarning("VertexGenerator") << "Warning BeamProfileVtxGenerator:"
                                       << " Illegal resolution in X " << s << "- set to default value 0 cm";
    fSigmaX = 0;
  }
}

sigmaY()

void BeamProfileVtxGenerator::sigmaY

(

double

s = 1.0

)

set resolution in Y in cm

Definition at line 169 of file BeamProfileVtxGenerator.cc.

References fSigmaY, and alignCSCRings::s.

Referenced by BeamProfileVtxGenerator().

                                             {
  if (s >= 0) {
    fSigmaY = s;
  } else {
    edm::LogWarning("VertexGenerator") << "Warning BeamProfileVtxGenerator:"
                                       << " Illegal resolution in Y " << s << "- set to default value 0 cm";
    fSigmaY = 0;
  }
}

Member Data Documentation

fdistn

std::vector<double> BeamProfileVtxGenerator::fdistn

private

Definition at line 64 of file BeamProfileVtxGenerator.h.

Referenced by BeamProfileVtxGenerator(), and newVertex().

fEta

double BeamProfileVtxGenerator::fEta

private

Definition at line 58 of file BeamProfileVtxGenerator.h.

Referenced by BeamProfileVtxGenerator(), and eta().

ffile

bool BeamProfileVtxGenerator::ffile

private

Definition at line 62 of file BeamProfileVtxGenerator.h.

Referenced by BeamProfileVtxGenerator(), and newVertex().

fMeanX

double BeamProfileVtxGenerator::fMeanX

private

Definition at line 57 of file BeamProfileVtxGenerator.h.

Referenced by BeamProfileVtxGenerator(), meanX(), and newVertex().

fMeanY

double BeamProfileVtxGenerator::fMeanY

private

Definition at line 57 of file BeamProfileVtxGenerator.h.

Referenced by BeamProfileVtxGenerator(), meanY(), and newVertex().

fMeanZ

double BeamProfileVtxGenerator::fMeanZ

private

Definition at line 57 of file BeamProfileVtxGenerator.h.

Referenced by beamPos(), BeamProfileVtxGenerator(), and newVertex().

fPhi

double BeamProfileVtxGenerator::fPhi

private

Definition at line 58 of file BeamProfileVtxGenerator.h.

Referenced by BeamProfileVtxGenerator(), newVertex(), and phi().

fPsi

double BeamProfileVtxGenerator::fPsi

private

Definition at line 60 of file BeamProfileVtxGenerator.h.

Referenced by BeamProfileVtxGenerator(), newVertex(), and psi().

fSigmaX

double BeamProfileVtxGenerator::fSigmaX

private

Definition at line 56 of file BeamProfileVtxGenerator.h.

Referenced by BeamProfileVtxGenerator(), newVertex(), and sigmaX().

fSigmaY

double BeamProfileVtxGenerator::fSigmaY

private

Definition at line 56 of file BeamProfileVtxGenerator.h.

Referenced by BeamProfileVtxGenerator(), newVertex(), and sigmaY().

fTheta

double BeamProfileVtxGenerator::fTheta

private

Definition at line 58 of file BeamProfileVtxGenerator.h.

Referenced by BeamProfileVtxGenerator(), eta(), and newVertex().

fTimeOffset

double BeamProfileVtxGenerator::fTimeOffset

private

Definition at line 65 of file BeamProfileVtxGenerator.h.

Referenced by BeamProfileVtxGenerator(), and newVertex().

fType

bool BeamProfileVtxGenerator::fType

private

Definition at line 62 of file BeamProfileVtxGenerator.h.

Referenced by BeamProfileVtxGenerator(), newVertex(), and setType().

nBinx

int BeamProfileVtxGenerator::nBinx

private

Definition at line 63 of file BeamProfileVtxGenerator.h.

Referenced by BeamProfileVtxGenerator(), and newVertex().

nBiny

int BeamProfileVtxGenerator::nBiny

private

Definition at line 63 of file BeamProfileVtxGenerator.h.

Referenced by BeamProfileVtxGenerator(), and newVertex().