EntryAngleAffAngles Class Reference

#include <EntryAngleAffAngles.h>

Inheritance diagram for EntryAngleAffAngles:

Public Member Functions
ALIdouble	checkDiff (const CLHEP::Hep3Vector &axis, const CLHEP::Hep3Vector &axisOrig, const std::vector< double > &localrot, const std::vector< double > &localrotorig) const
void	displace (ALIdouble disp) override
void	displaceOriginal (ALIdouble disp) override
void	displaceOriginalOriginal (ALIdouble disp) override
	EntryAngleAffAngles (const ALIstring &type)
virtual void	FillName (const ALIstring &name)
ALIdouble	valueDisplaced () const override
	~EntryAngleAffAngles () override
Public Member Functions inherited from EntryAngle
	EntryAngle (const ALIstring &type)
ALIdouble	OutputSigmaDimensionFactor () const override
ALIdouble	OutputValueDimensionFactor () const override
ALIdouble	SigmaDimensionFactor () const override
ALIdouble	startingDisplacement () override
ALIdouble	ValueDimensionFactor () const override
	~EntryAngle () override
Public Member Functions inherited from Entry
void	addFittedDisplacementToValue (const ALIdouble val)
	Entry (const ALIstring &type)
void	fill (const std::vector< ALIstring > &wordlist)
virtual void	fillName (const ALIstring &name)
void	fillNull ()
ALIint	fitPos () const
EntryDim	getDimType () const
ALIdouble	lastAdditionToValueDisplacementByFitting () const
const ALIstring	longName () const
const ALIstring &	name () const
OpticalObject *	OptOCurrent () const
ALIint	quality () const
void	resetValueDisplacementByFitting ()
void	setFitPos (const ALIint fitpos)
void	setLastAdditionToValueDisplacementByFitting (const ALIdouble val)
void	setName (const ALIstring &name)
void	setOptOCurrent (OpticalObject *opto)
void	setQuality (ALIuint qual)
void	setSigma (ALIdouble sig)
void	setType (ALIstring type)
void	setValue (ALIdouble val)
ALIdouble	sigma () const
ALIdouble	sigmaOriginalOriginal () const
void	substractToHalfFittedDisplacementToValue ()
const ALIstring &	type () const
ALIdouble	value () const
ALIdouble	valueDisplacementByFitting () const
virtual ALIdouble	valueInGlobalReferenceFrame () const
ALIdouble	valueOriginalOriginal () const
virtual	~Entry ()

Additional Inherited Members
Protected Attributes inherited from Entry
ALIint	fitPos_
ALIdouble	lastAdditionToValueDisplacementByFitting_
ALIstring	name_
OpticalObject *	OptOCurrent_
ALIuint	quality_
ALIdouble	sigma_
ALIdouble	sigmaOriginalOriginal_
EntryDim	theDimType
ALIstring	type_
ALIdouble	value_
ALIdouble	valueDisplacementByFitting_
ALIdouble	valueOriginalOriginal_

Detailed Description

Definition at line 15 of file EntryAngleAffAngles.h.

Constructor & Destructor Documentation

EntryAngleAffAngles()

EntryAngleAffAngles::EntryAngleAffAngles

(

const ALIstring &

type

)

Definition at line 13 of file EntryAngleAffAngles.cc.

: EntryAngle(type) {}

~EntryAngleAffAngles()

EntryAngleAffAngles::~EntryAngleAffAngles ( )

inlineoverride

Definition at line 18 of file EntryAngleAffAngles.h.

{}

Member Function Documentation

checkDiff()

ALIdouble EntryAngleAffAngles::checkDiff	(	const CLHEP::Hep3Vector &	axis,
		const CLHEP::Hep3Vector &	axisOrig,
		const std::vector< double > &	localrot,
		const std::vector< double > &	localrotorig
	)		const

Definition at line 139 of file EntryAngleAffAngles.cc.

References Entry::name().

                                                                                      {
  CLHEP::Hep3Vector axis = _axis;
  CLHEP::Hep3Vector axisOrig = _axisOrig;
  int inam = 0;
  if (name() == "angles_X") {
    inam = 1;
  } else if (name() == "angles_Y") {
    inam = 2;
  } else if (name() == "angles_Z") {
    inam = 3;
  }
  switch (inam) {
    case 1:
      axis.rotateX(localrot[0]);
      axisOrig.rotateX(localrotorig[0]);
      [[fallthrough]];
    case 2:
      axis.rotateY(localrot[1]);
      axisOrig.rotateY(localrotorig[1]);
      [[fallthrough]];
    case 3:
      axis.rotateZ(localrot[2]);
      axisOrig.rotateZ(localrotorig[2]);
      break;
  }

  ALIdouble ang = axis.angle(axisOrig);

  /*  }else if( name() == "angles_Y" ) {
    return localrot[1] - localrotorig[1];
  }else if( name() == "angles_Z" ) {
    return localrot[2] - localrotorig[2];
    }*/

  return ang;
}

displace()

void EntryAngleAffAngles::displace ( ALIdouble  disp )

overridevirtual

Reimplemented from Entry.

Definition at line 25 of file EntryAngleAffAngles.cc.

References OpticalObject::displaceRmGlobAroundGlobal(), OpticalObject::displaceRmGlobAroundLocal(), GlobalOptionMgr::getInstance(), GlobalOptionMgr::GlobalOptions(), Entry::name(), Entry::name_, Entry::OptOCurrent(), XCoor, YCoor, and ZCoor.

                                                 {
  XYZcoor coor = XCoor;
  ALIint namelength = name().length() - 1;
  //-  std::cout << this << "ENtryAnglesAffAngle" << name_ << namelength <<std::endl;
  if (name_[namelength] == 'X') {
    coor = XCoor;
  } else if (name_[namelength] == 'Y') {
    coor = YCoor;
  } else if (name_[namelength] == 'Z') {
    coor = ZCoor;
  }

  GlobalOptionMgr* gomgr = GlobalOptionMgr::getInstance();
  if (gomgr->GlobalOptions()["rotateAroundLocal"] == 0) {
    OptOCurrent()->displaceRmGlobAroundGlobal(OptOCurrent(), coor, disp);
  } else {
    OptOCurrent()->displaceRmGlobAroundLocal(OptOCurrent(), coor, disp);
  }
}

displaceOriginal()

void EntryAngleAffAngles::displaceOriginal ( ALIdouble  disp )

overridevirtual

Reimplemented from Entry.

Definition at line 46 of file EntryAngleAffAngles.cc.

References OpticalObject::displaceRmGlobOriginal(), Entry::name(), Entry::name_, Entry::OptOCurrent(), XCoor, YCoor, and ZCoor.

                                                         {
  ALIint namelength = name().length() - 1;
  //-  std::cout << this << "ENtryAnglesAffAngle" << name_ << namelength <<std::endl;
  if (name_[namelength] == 'X') {
    //-    std::cout << "displaX";
    OptOCurrent()->displaceRmGlobOriginal(OptOCurrent(), XCoor, disp);
  } else if (name_[namelength] == 'Y') {
    //-    std::cout << "displaY";
    OptOCurrent()->displaceRmGlobOriginal(OptOCurrent(), YCoor, disp);
  } else if (name_[namelength] == 'Z') {
    OptOCurrent()->displaceRmGlobOriginal(OptOCurrent(), ZCoor, disp);
  }
}

displaceOriginalOriginal()

void EntryAngleAffAngles::displaceOriginalOriginal ( ALIdouble  disp )

overridevirtual

Reimplemented from Entry.

Definition at line 61 of file EntryAngleAffAngles.cc.

References gather_cfg::cout, ALIUtils::debug, OpticalObject::displaceRmGlobOriginalOriginal(), Entry::name(), Entry::name_, Entry::OptOCurrent(), XCoor, YCoor, and ZCoor.

                                                                 {
  ALIint namelength = name().length() - 1;
  if (ALIUtils::debug >= 5)
    std::cout << this << "ENtryAnglesAffAngle displaceOriginalOriginal" << name_ << std::endl;
  if (name_[namelength] == 'X') {
    //-    std::cout << "displaX";
    OptOCurrent()->displaceRmGlobOriginalOriginal(OptOCurrent(), XCoor, disp);
  } else if (name_[namelength] == 'Y') {
    //-    std::cout << "displaY";
    OptOCurrent()->displaceRmGlobOriginalOriginal(OptOCurrent(), YCoor, disp);
  } else if (name_[namelength] == 'Z') {
    OptOCurrent()->displaceRmGlobOriginalOriginal(OptOCurrent(), ZCoor, disp);
  }
}

FillName()

void EntryAngleAffAngles::FillName ( const ALIstring &  name )

virtual

Definition at line 16 of file EntryAngleAffAngles.cc.

References Entry::name(), groupFilesInBlocks::nn, and Entry::setName().

                                                        {
  ALIstring nn = "Angles ";
  nn += name;
  setName(nn);
}

valueDisplaced()

ALIdouble EntryAngleAffAngles::valueDisplaced ( ) const

overridevirtual

Reimplemented from Entry.

Definition at line 77 of file EntryAngleAffAngles.cc.

References gather_cfg::cout, ALIUtils::debug, ALIUtils::dumprm(), OpticalObject::getRotationAnglesFromMatrix(), Entry::name(), OpticalObject::name(), Entry::OptOCurrent(), OpticalObject::parent(), class-composition::parent, OpticalObject::rmGlob(), and OpticalObject::rmGlobOriginal().

                                                    {
  //  std::cout << "  EntryAngleAffAngles::valueDisplaced() parent " << OptOCurrent()->parent()->name() << std::endl;
  //  ALIUtils::dumprm( OptOCurrent()->parent()->rmGlob() , " parent RmGlob ");
  if (ALIUtils::debug >= 5) {
    std::cout << "  EntryAngleAffAngles::valueDisplaced() parent opto  " << OptOCurrent()->parent()->name()
              << std::endl;
    ALIUtils::dumprm(OptOCurrent()->rmGlob(), " RmGlob ");
    ALIUtils::dumprm(OptOCurrent()->rmGlobOriginal(), " RmGlobOriginal ");
    ALIUtils::dumprm(OptOCurrent()->parent()->rmGlobOriginal(), " parent RmGlobOriginal ");
  }

  CLHEP::HepRotation diffRm = OptOCurrent()->rmGlob() * OptOCurrent()->rmGlobOriginal().inverse();
  CLHEP::HepRotation rmLocal = diffRm * OptOCurrent()->parent()->rmGlobOriginal().inverse();
  std::vector<double> localrot =
      OptOCurrent()->getRotationAnglesFromMatrix(rmLocal, OptOCurrent()->CoordinateEntryList());
  if (ALIUtils::debug >= 5) {
    ALIUtils::dumprm(diffRm, " diffRm ");
    ALIUtils::dumprm(rmLocal, " rmLocal ");
    std::cout << " localrot " << localrot[0] << " " << localrot[1] << " " << localrot[2] << std::endl;
  }

  if (name() == "angles_X") {
    return localrot[0];
  } else if (name() == "angles_Y") {
    return localrot[1];
  } else if (name() == "angles_Z") {
    return localrot[2];
  }

  /*
  CLHEP::HepRotation rmLocalOrig = OptOCurrent()->parent()->rmGlobOriginal().inverse() *  OptOCurrent()->rmGlobOriginal();

  CLHEP::HepRotation rmLocal = OptOCurrent()->parent()->rmGlob().inverse() * OptOCurrent()->rmGlob();
  std::vector<double> localrot = OptOCurrent()->getRotationAnglesFromMatrix(  rmLocal, OptOCurrent()->CoordinateEntryList() );

  std::cout << " localrot " << localrot[0] << " " << localrot[1] << " " << localrot[2] << std::endl;
  std::cout << " localrotorig " << localrotorig[0] << " " << localrotorig[1] << " " << localrotorig[2] << std::endl;
  ALIdouble diff;
  CLHEP::Hep3Vector Xaxis(0.,0.,1.);
  Xaxis = OptOCurrent()->parent()->rmGlob() * Xaxis;
  CLHEP::Hep3Vector XaxisOrig(0.,0.,1.);
  XaxisOrig = OptOCurrent()->parent()->rmGlobOriginal() * XaxisOrig;

  diff = fabs( checkDiff( Xaxis, XaxisOrig, localrot, localrotorig ) );

  //maybe X is not a good axis because the rotation is done precisely around X
  if( diff <= 1.E-9 ){
    CLHEP::Hep3Vector Yaxis(0.,1.,0.);
    Yaxis = OptOCurrent()->parent()->rmGlob() * Yaxis;
    CLHEP::Hep3Vector YaxisOrig(0.,1.,0.);
    YaxisOrig = OptOCurrent()->parent()->rmGlobOriginal() * YaxisOrig;

    diff = fabs( checkDiff( Yaxis, YaxisOrig, localrot, localrotorig ) );
  }

  return diff;
  */

  return 0.;  // to avoid warning
}