#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)

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...

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 Member Functions
	BeamProfileVtxGenerator (const BeamProfileVtxGenerator &p)=delete

BeamProfileVtxGenerator &	operator= (const BeamProfileVtxGenerator &rhs)=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.

                                                                          : 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";
 }

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.

◆ ~BeamProfileVtxGenerator()

BeamProfileVtxGenerator::~BeamProfileVtxGenerator ( )

override

Definition at line 77 of file BeamProfileVtxGenerator.cc.

77 {}

◆ BeamProfileVtxGenerator() [2/2]

BeamProfileVtxGenerator::BeamProfileVtxGenerator ( const BeamProfileVtxGenerator & p )

privatedelete

Copy constructor

Member Function Documentation

◆ beamPos()

void BeamProfileVtxGenerator::beamPos ( double m = 0 )

inline

set mean in Z in cm

Definition at line 40 of file BeamProfileVtxGenerator.h.

40 { fMeanZ = m; }

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

Referenced by BeamProfileVtxGenerator().

◆ eta()

void BeamProfileVtxGenerator::eta ( double m = 0 )

set eta

Definition at line 179 of file BeamProfileVtxGenerator.cc.

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

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

Referenced by Particle.Particle::__str__(), BeamProfileVtxGenerator(), Jet.Jet::jetID(), and Jet.Jet::puJetId().

◆ 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 28 of file BeamProfileVtxGenerator.h.

28 { return nullptr; }

◆ meanX()

void BeamProfileVtxGenerator::meanX ( double m = 0 )

inline

set mean in X in cm

Definition at line 36 of file BeamProfileVtxGenerator.h.

36 { fMeanX = m; }

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

Referenced by BeamProfileVtxGenerator().

◆ meanY()

void BeamProfileVtxGenerator::meanY ( double m = 0 )

inline

set mean in Y in cm

Definition at line 38 of file BeamProfileVtxGenerator.h.

38 { fMeanY = m; }

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

Referenced by BeamProfileVtxGenerator().

◆ 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.

                                                                                      {
   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);
   ;
 }

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

◆ operator=()

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

privatedelete

Copy assignment operator

◆ phi()

void BeamProfileVtxGenerator::phi ( double m = 0 )

inline

set phi in radian

Definition at line 45 of file BeamProfileVtxGenerator.h.

45 { fPhi = m; }

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

Referenced by Particle.Particle::__str__(), BeamProfileVtxGenerator(), and ntupleDataFormat.Track::phiPull().

◆ psi()

void BeamProfileVtxGenerator::psi ( double m = 999 )

inline

set psi in radian

Definition at line 47 of file BeamProfileVtxGenerator.h.

47 { fPsi = m; }

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

Referenced by BeamProfileVtxGenerator().

◆ setType()

void BeamProfileVtxGenerator::setType ( bool m = true )

set type

Definition at line 184 of file BeamProfileVtxGenerator.cc.

184 { fType = s; }

References fType, and alignCSCRings::s.

Referenced by BeamProfileVtxGenerator().

◆ sigmaX()

void BeamProfileVtxGenerator::sigmaX ( double s = 1.0 )

set resolution in X in cm

Definition at line 159 of file BeamProfileVtxGenerator.cc.

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

References fSigmaX, and alignCSCRings::s.

Referenced by BeamProfileVtxGenerator().

◆ sigmaY()

void BeamProfileVtxGenerator::sigmaY ( double s = 1.0 )

set resolution in Y in cm

Definition at line 169 of file BeamProfileVtxGenerator.cc.

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