OptOSensor2D Class Reference

#include <OptOSensor2D.h>

Inheritance diagram for OptOSensor2D:

Public Member Functions
void	constructSolidShape () override
ALIdouble *	convertPointToLocalCoordinates (const CLHEP::Hep3Vector &point)
void	defaultBehaviour (LightRay &lightray, Measurement &meas) override
void	detailedTraversesLightRay (LightRay &lightray) override
void	fastTraversesLightRay (LightRay &lightray) override
void	fillExtraEntry (std::vector< ALIstring > &wordlist) override
void	makeMeasurement (LightRay &lightray, Measurement &meas) override
	OptOSensor2D ()
	OptOSensor2D (OpticalObject *parent, const ALIstring &type, const ALIstring &name, const ALIbool copy_data)
	~OptOSensor2D () override
Public Member Functions inherited from OpticalObject
void	addCoordinateEntryToList (Entry *entry)
void	addExtraEntryToList (Entry *entry)
void	addExtraEntryValueOriginalOriginalToList (ALIdouble entry_value)
void	addExtraEntryValueOriginalToList (ALIdouble entry_value)
void	addExtraEntryValueToList (ALIdouble entry_value)
double	addPii (double val)
double	approxTo0 (double val)
CLHEP::HepRotation	buildRmFromEntryValuesOriginalOriginal ()
const CLHEP::Hep3Vector &	centreGlob () const
const CLHEP::Hep3Vector &	centreGlobal () const
const CLHEP::Hep3Vector &	centreGlobOriginal () const
const CLHEP::Hep3Vector &	centreGlobOriginalOriginal () const
const CLHEP::Hep3Vector	centreLocal () const
int	checkMatrixEquations (double angleX, double angleY, double angleZ, CLHEP::HepRotation *rot=nullptr)
void	construct ()
void	constructFromOptAligInfo (const OpticalAlignInfo &oaInfo)
virtual void	constructMaterial ()
const std::vector< Entry * > &	CoordinateEntryList () const
void	createComponentOptOsFromOptAlignInfo ()
virtual void	detailedDeviatesLightRay (LightRay &lightray)
double	diff2pi (double ang1, double ang2)
void	displaceCentreGlob (const XYZcoor coor, const ALIdouble disp)
void	displaceCentreGlob (const CLHEP::Hep3Vector &dispVec)
void	displaceCentreGlobOriginal (const XYZcoor coor, const ALIdouble disp)
void	displaceCentreGlobOriginal (const CLHEP::Hep3Vector &dispVec)
void	displaceCentreGlobOriginalOriginal (const XYZcoor coor, const ALIdouble disp)
void	displaceCentreGlobOriginalOriginal (const CLHEP::Hep3Vector &dispVec)
void	displaceExtraEntry (const ALIuint entryNo, const ALIdouble disp)
void	displaceExtraEntryOriginal (const ALIuint entryNo, const ALIdouble disp)
void	displaceExtraEntryOriginalOriginal (const ALIuint entryNo, const ALIdouble disp)
void	displaceRmGlobAroundGlobal (OpticalObject *opto1stRotated, const XYZcoor coor, const ALIdouble disp)
void	displaceRmGlobAroundLocal (OpticalObject *opto1stRotated, const XYZcoor coor, const ALIdouble disp)
void	displaceRmGlobOriginal (const OpticalObject *opto1stRotated, const XYZcoor coor, const ALIdouble disp)
void	displaceRmGlobOriginalOriginal (const OpticalObject *opto1stRotated, const XYZcoor coor, const ALIdouble disp)
bool	eq2ang (double ang1, double ang2)
const std::vector< Entry * > &	ExtraEntryList () const
const ALIint	extraEntryNo (const ALIstring &entry_name) const
std::vector< ALIdouble > &	ExtraEntryValueList ()
const std::vector< ALIdouble > &	ExtraEntryValueOriginalList ()
const std::vector< ALIdouble > &	ExtraEntryValueOriginalOriginalList ()
virtual void	fastDeviatesLightRay (LightRay &lightray)
virtual void	fillIguana ()
virtual void	fillVRML ()
const ALIdouble	findExtraEntryValue (const ALIstring &eename) const
const ALIbool	findExtraEntryValueIfExists (const ALIstring &eename, ALIdouble &value) const
const ALIdouble	findExtraEntryValueMustExist (const ALIstring &eename) const
const ALIuint	getCmsswID () const
std::vector< ALIstring >	getCoordinateFromOptAlignParam (const OpticalAlignParam &oaParam)
CLHEP::Hep3Vector	getDisplacementInLocalCoordinates (const XYZcoor coor, const ALIdouble disp)
const double	getEntryCentre (const XYZcoor coor) const
const double	getEntryCentre (const ALIstring &coor) const
const double	getEntryRMangle (const XYZcoor coor) const
const double	getEntryRMangle (const ALIstring &coor) const
std::vector< double >	getLocalRotationAngles (const std::vector< Entry *> &entries) const
CocoaMaterialElementary *	getMaterial () const
ALIPlane	getPlate (const ALIbool forwardPlate, const ALIbool applyWedge)
std::vector< double >	getRotationAnglesFromMatrix (CLHEP::HepRotation &rmLocal, const std::vector< Entry *> &entries) const
std::vector< double >	getRotationAnglesInOptOFrame (const OpticalObject *optoAncestor, const std::vector< Entry *> &entries) const
CocoaSolidShape *	getSolidShape () const
CLHEP::Hep3Vector	getZAxis ()
const ALIuint	ID () const
const ALIstring	longName () const
Measurement *	meas ()
const ALIstring &	name () const
	OpticalObject ()
	OpticalObject (OpticalObject *parent, const ALIstring &type, const ALIstring &name, const ALIbool copy_data)
const OpticalObject *	parent () const
virtual void	participateInMeasurement (LightRay &lightray, Measurement &meas, const ALIstring &behav)
void	propagateGlobalRMOriginalOriginalChangeToChildren (const CLHEP::HepRotation &rmorioriold, const CLHEP::HepRotation &rmoriorinew)
void	resetGlobalCoordinates ()
void	resetOriginalOriginalCoordinates ()
const CLHEP::HepRotation &	rmGlob () const
const CLHEP::HepRotation &	rmGlobOriginal () const
const CLHEP::HepRotation &	rmGlobOriginalOriginal () const
const CLHEP::HepRotation	rmLocal () const
void	setCmsswID (ALIuint id)
void	setExtraEntryValue (const ALIuint entryNo, const ALIdouble disp)
void	setGlobalCoordinates ()
void	setGlobalRMOriginalOriginal (const CLHEP::HepRotation &rmoriori)
void	setID (ALIuint id)
void	setMeas (Measurement *meas)
	set current measurement More...
void	setOriginalEntryValues ()
void	setRmGlobal (const CLHEP::HepRotation &rm)
void	setRmGlobalOriginal (const CLHEP::HepRotation &rm)
void	setType (const ALIstring &type)
const ALIstring	shortName () const
const ALIstring &	type () const
virtual void	userDefinedBehaviour (LightRay &lightray, Measurement &meas, const ALIstring &behav)
virtual	~OpticalObject ()

Private Attributes
DeviationsFromFileSensor2D *	deviFromFile
ALIbool	fdevi_from_file

Additional Inherited Members
Protected Attributes inherited from OpticalObject
CocoaMaterialElementary *	theMaterial
CocoaSolidShape *	theSolidShape
ALIint	verbose

Detailed Description

Definition at line 19 of file OptOSensor2D.h.

Constructor & Destructor Documentation

OptOSensor2D() [1/2]

OptOSensor2D::OptOSensor2D ( )

inline

Definition at line 22 of file OptOSensor2D.h.

{};

OptOSensor2D() [2/2]

OptOSensor2D::OptOSensor2D ( OpticalObject *  parent, const ALIstring &  type, const ALIstring &  name, const ALIbool  copy_data  )

inline

Definition at line 23 of file OptOSensor2D.h.

      : OpticalObject(parent, type, name, copy_data), fdevi_from_file(false){};

~OptOSensor2D()

OptOSensor2D::~OptOSensor2D ( )

inlineoverride

Definition at line 25 of file OptOSensor2D.h.

{};

Member Function Documentation

constructSolidShape()

void OptOSensor2D::constructSolidShape ( )

overridevirtual

Reimplemented from OpticalObject.

Definition at line 452 of file OptOSensor2D.cc.

References GlobalOptionMgr::getGlobalOptionValue(), GlobalOptionMgr::getInstance(), and visualization-live-secondInstance_cfg::m.

                                       {
  ALIdouble go;
  GlobalOptionMgr* gomgr = GlobalOptionMgr::getInstance();
  gomgr->getGlobalOptionValue("VisScale", go);

  theSolidShape = new CocoaSolidShapeBox(
      "Box", go * 4. * cm / m, go * 4. * cm / m, go * 1. * cm / m);  //COCOA internal units are meters
}

convertPointToLocalCoordinates()

ALIdouble * OptOSensor2D::convertPointToLocalCoordinates

(

const CLHEP::Hep3Vector &

point

)

Definition at line 399 of file OptOSensor2D.cc.

References gather_cfg::cout, ALIUtils::debug, ALIUtils::dump3v(), ALIUtils::dumprm(), and point.

                                                                                    {
  ALIdouble* interslc = new ALIdouble[2];

  //----- X value
  CLHEP::HepRotation rmt = rmGlob();
  CLHEP::Hep3Vector XAxism(1., 0., 0.);
  XAxism *= rmt;
  if (ALIUtils::debug >= 5)
    ALIUtils::dump3v((this)->centreGlob(), "centre glob sensor2D");
  if (ALIUtils::debug >= 5)
    ALIUtils::dumprm(rmt, "rotation matrix sensor2D");
  //t  ALIUtils::dump3v(point - (this)->centreGlob() , "inters - (this)->centreGlob()");
  if (ALIUtils::debug >= 5)
    ALIUtils::dump3v(XAxism, "XAxism");
  interslc[0] = (point - (this)->centreGlob()) * XAxism;

  //----- Y value
  CLHEP::Hep3Vector YAxism(0., 1., 0.);
  YAxism *= rmt;
  if (ALIUtils::debug >= 5)
    ALIUtils::dump3v(YAxism, "YAxism");
  interslc[1] = (point - (this)->centreGlob()) * YAxism;

  if (ALIUtils::debug >= 5) {
    std::cout << " intersection in local coordinates: X= " << interslc[0] << "  Y= " << interslc[1] << std::endl;
    ALIUtils::dump3v(point - (this)->centreGlob(), " inters - centre ");
  }
  return interslc;
}

defaultBehaviour()

void OptOSensor2D::defaultBehaviour ( LightRay &  lightray, Measurement &  meas  )

overridevirtual

Reimplemented from OpticalObject.

Definition at line 35 of file OptOSensor2D.cc.

{ makeMeasurement(lightray, meas); }

detailedTraversesLightRay()

void OptOSensor2D::detailedTraversesLightRay ( LightRay &  lightray )

overridevirtual

Reimplemented from OpticalObject.

Definition at line 247 of file OptOSensor2D.cc.

References DeviationsFromFileSensor2D::apply(), EcnaPython_AdcPeg12_S1_10_R170298_1_0_150_Dee0::cerr, gather_cfg::cout, ALIUtils::debug, LightRay::dumpData(), LightRay::intersect(), Skims_PA_cff::name, LightRay::point(), LightRay::refract(), DeviationsFromFileSensor2D::setApply(), LightRay::setPoint(), and ApeEstimator_cff::width.

                                                               {
  if (ALIUtils::debug >= 4)
    std::cout << "%%% LR: DETAILED TRAVERSES SENSOR2D  " << name() << std::endl;
  if (DeviationsFromFileSensor2D::apply() && fdevi_from_file) {
    DeviationsFromFileSensor2D::setApply(false);
    //- std::cout << "fdeviFromFile" << fdevi_from_file << std::endl;
    if (ALIUtils::debug >= 0)
      std::cerr << "!!WARNING: sensor " << name()
                << " has read deviation from file and it will not be taken into account. Please use FAST TRAVERSES"
                << deviFromFile << std::endl;
  }

  //---------- If width is 0, just keep the intersection point
  ALIdouble width = findExtraEntryValue("width");
  if (width == 0) {
    //---------- Get intersection
    CLHEP::Hep3Vector ZAxis(0., 0, 1.);
    CLHEP::HepRotation rmt = rmGlob();
    ZAxis = rmt * ZAxis;
    lightray.intersect(ALIPlane(centreGlob(), ZAxis));
    CLHEP::Hep3Vector inters = lightray.point();
    lightray.setPoint(inters);
    if (ALIUtils::debug >= 2) {
      lightray.dumpData("LightRay Sensor2D traversed: ");
    }
    return;
  }

  if (ALIUtils::debug >= 4)
    std::cout << std::endl << "$$$ LR: REFRACTION IN FORWARD PLATE " << std::endl;
  //---------- Get forward plate
  ALIPlane plate = getPlate(true, true);
  //---------- Refract while entering object
  ALIdouble refra_ind1 = 1.;
  ALIdouble refra_ind2 = findExtraEntryValueMustExist("refra_ind");
  lightray.refract(plate, refra_ind1, refra_ind2);

  if (ALIUtils::debug >= 4)
    std::cout << std::endl << "$$$ LR: REFRACTION IN BACKWARD PLATE " << std::endl;
  //---------- Get backward plate
  plate = getPlate(false, true);
  //---------- Refract while exiting splitter
  lightray.refract(plate, refra_ind2, refra_ind1);

  CLHEP::Hep3Vector inters = lightray.point();
  lightray.setPoint(inters);

  if (ALIUtils::debug >= 4) {
    lightray.dumpData("LightRay Sensor2D traversed: ");
  }
}

fastTraversesLightRay()

void OptOSensor2D::fastTraversesLightRay ( LightRay &  lightray )

overridevirtual

Reimplemented from OpticalObject.

Definition at line 119 of file OptOSensor2D.cc.

References DeviationsFromFileSensor2D::apply(), gather_cfg::cout, ALIUtils::debug, hgcalPerformanceValidation::df, LightRay::dumpData(), Exception, cuy::ii, LightRay::intersect(), ALIUtils::LengthValueDimensionFactor(), Model::MeasurementList(), Skims_PA_cff::name, LightRay::point(), LightRay::setPoint(), LightRay::shiftAndDeviateWhileTraversing(), Measurement::value(), and verbose.

                                                           {
  verbose = ALIUtils::debug;
  if (ALIUtils::debug >= 2)
    std::cout << "LR: FAST TRAVERSES SENSOR2D  " << name() << std::endl;

  //---------- Shift and Deviate

  //---------- Get intersection
  CLHEP::Hep3Vector ZAxis(0., 0, 1.);
  CLHEP::HepRotation rmt = rmGlob();
  ZAxis = rmt * ZAxis;
  lightray.intersect(ALIPlane(centreGlob(), ZAxis));
  CLHEP::Hep3Vector inters = lightray.point();

  ALIdouble deviX, deviY, devi;
  // if deviationsFromFile are applied and this sensors has them get the deviations that corresponds to the intersection point.
  // Set this deviation as original one, as it will not be changed for derivatives
  if (DeviationsFromFileSensor2D::apply() && fdevi_from_file) {
    //-    std::cout << " DeviationsFromFileSensor2D::apply() " << DeviationsFromFileSensor2D::apply() << std::endl;
    if (ALIUtils::debug >= 4) {
      std::cout << "fdeviFromFile" << fdevi_from_file << std::endl;
      //-      std::cout << "deviFromFile" << deviFromFile << std::endl;
    }
    //--------- get measurement value of the current sensor
    std::vector<Measurement*>& measv = Model::MeasurementList();
    unsigned int ii;
    Measurement* omeas = nullptr;
    for (ii = 0; ii < measv.size(); ii++) {
      //-   std::cout << " sensor2d finding meas " <<  measv[ii]->sensorName() << " " << name() << std::endl;
      if (measv[ii]->sensorName() == name()) {
        omeas = measv[ii];
        break;
      }
    }
    if (omeas == nullptr) {
      throw cms::Exception("LogicError") << "@SUB=OptOSensor2D::fastTraversesLightRay\n"
                                         << "meas " << name() << " not found";
    }

    ALIdouble interslcx = omeas->value(0);
    ALIdouble interslcy = omeas->value(1);
    if (ALIUtils::debug >= 5)
      std::cout << " interslcx " << interslcx << " interslcy " << interslcy << std::endl;
    //----- transform in milimeters and positive
    //mum    interslcx = interslcx*1.E6 + 10000.;
    //mum    interslcy = interslcy*1.E6 + 10000.;
    ALIdouble df = ALIUtils::LengthValueDimensionFactor();
    interslcx = interslcx / df + 0.010 / df;
    interslcy = interslcy / df + 0.010 / df;
    if (ALIUtils::debug >= 5)
      std::cout << " interslcx " << interslcx << " interslcy " << interslcy << std::endl;

    //---------- Get deviations from file (they are in microrads)
    std::pair<ALIdouble, ALIdouble> devis = deviFromFile->getDevis(interslcx, interslcy);
    deviX = devis.second;
    deviY = -devis.first;
    // deviX = devis.first;
    //  deviY = devis.second;

    //o  deviX = *((deviFromFile->deviationsInX()).begin() + pointY*deviFromFile->nPoints() + pointX) * 1.E-6;;
    //o  deviY = *((deviFromFile->deviationsInY().begin()) + pointY*deviFromFile->nPoints() + pointX) * 1.E-6;

    //---------- Set this deviation value as original one, as it will not be changed for derivatives (change the entry and also the ExtraEntryValueOriginalList())
    ALIuint entryNo = extraEntryNo("deviX");
    if (verbose >= 3)
      std::cout << "entrynox" << entryNo << name() << verbose << std::endl;
    Entry* entryDeviX = *(ExtraEntryList().begin() + entryNo);
    entryDeviX->setValue(deviX);
    //-    std::vector< ALIdouble >::const_iterator eevolite = static_cast<std::vector< ALIdouble >::iterator>( ExtraEntryValueOriginalList().begin() );
    std::vector<ALIdouble> eevil = ExtraEntryValueOriginalList();
    //-    std::vector< ALIdouble >::const_iterator eevolite = ( ExtraEntryValueOriginalList().begin() );
    std::vector<ALIdouble>::iterator eevolite = eevil.begin();

    *(eevolite + entryNo) = deviX;
    if (verbose >= 3)
      std::cout << " entryDeviX name " << entryDeviX->name() << entryDeviX->value() << std::endl;
    entryNo = extraEntryNo("deviY");
    Entry* entryDeviY = *(ExtraEntryList().begin() + entryNo);
    //- std::cout << "entrynoy" << entryNo << name() << std::endl;
    entryDeviY->setValue(deviY);
    *(eevolite + entryNo) = deviY;
    //-  std::cout<< entryDeviY << " entryDeviY name " << entryDeviY->name() << entryDeviY->value() << std::endl;

  } else {
    deviX = findExtraEntryValue("deviX");
    deviY = findExtraEntryValue("deviY");

    //??? why previous does not work??
    if (fdevi_from_file) {
      if (ALIUtils::debug >= 5)
        std::cout << "fdeviFromFile" << fdevi_from_file << std::endl;
      ALIuint entryNo = extraEntryNo("deviX");
      Entry* entryDeviX = *(ExtraEntryList().begin() + entryNo);
      if (verbose >= 3)
        std::cout << entryDeviX << " entryDeviX name " << entryDeviX->name() << entryDeviX->value() << std::endl;
      deviX = entryDeviX->value();
      entryNo = extraEntryNo("deviY");
      Entry* entryDeviY = *(ExtraEntryList().begin() + entryNo);
      if (verbose >= 3)
        std::cout << entryDeviY << " entryDeviY name " << entryDeviY->name() << entryDeviY->value() << std::endl;
      deviY = entryDeviY->value();

    } else {
      ALIbool bb = findExtraEntryValueIfExists("devi", devi);
      if (bb) {
        deviX = devi;
        deviY = devi;
      }
      if (ALIUtils::debug >= 4) {
        std::cout << "devi " << devi << " devi x  " << deviX << " devi y  " << deviY << std::endl;
      }
    }
  }
  if (ALIUtils::debug >= 4) {
    std::cout << " devi x  " << deviX << " devi y  " << deviY << std::endl;
  }

  lightray.setPoint(inters);

  lightray.shiftAndDeviateWhileTraversing(this, 'T');
  if (ALIUtils::debug >= 2) {
    lightray.dumpData("Shifted and Deviated");
  }
}

fillExtraEntry()

void OptOSensor2D::fillExtraEntry ( std::vector< ALIstring > &  wordlist )

overridevirtual

Reimplemented from OpticalObject.

Definition at line 302 of file OptOSensor2D.cc.

References gather_cfg::cout, ALIUtils::debug, OpticalObject::fillExtraEntry(), ALIUtils::getFloat(), ALIFileIn::getInstance(), ALIFileIn::getWordsInLine(), and LaserClient_cfi::wl.

                                                                {
  if (ALIUtils::debug >= 5)
    std::cout << "OptOSensor2D fillExtraEntry wordlist[1] " << wordlist[1] << std::endl;
  //---------- check if it is deviation read from file
  fdevi_from_file = false;
  //-  std::cout << "WL " << wordlist[1]<< "WL " << wordlist[2]<< "WL " << wordlist[3] << std::endl;
  if (wordlist[1] == ALIstring("devi") && wordlist[2] == ALIstring("from_file")) {
    //---------- Open file
    ALIstring fnam;
    if (wordlist.size() >= 4) {
      fnam = wordlist[3];
    } else {
      //----- Build up file name if it does not exists
      fnam = "dat/devi-";
      fnam += shortName();
      fnam += ".dat";
    }

    ALIFileIn& ifdevi = ALIFileIn::getInstance(fnam);

    //----- Check that file exists
    if (ALIUtils::debug >= 4)
      std::cout << "Opening deviation file: " << fnam << std::endl;
    /*-    if( !ifdevi ) {
      std::cerr << " !!! Sensor2D Deviation file not found: " << fnam << " of object " << name() << std::endl;
      exit(1);
      }*/

    deviFromFile = new DeviationsFromFileSensor2D();
    fdevi_from_file = true;
    if (ALIUtils::debug >= 5)
      std::cout << "deviFromFile " << deviFromFile << std::endl;
    //----- Read header
    ALIstring sensor1_name, sensor2_name;
    ALIdouble sensor_dist;
    ALIdouble prec_deviX, prec_deviY;

    std::vector<ALIstring> wl;
    ifdevi.getWordsInLine(wl);
    sensor1_name = wl[0];
    sensor2_name = wl[1];
    sensor_dist = atof(wl[2].c_str());
    // 'c' means that light is traversing the glass
    if (sensor1_name[sensor1_name.size() - 2] == 'c') {
      sensor1_name = sensor1_name.substr(0, sensor1_name.size() - 1);
    }
    if (sensor2_name[sensor2_name.size() - 2] == 'c') {
      sensor2_name = sensor2_name.substr(0, sensor2_name.size() - 1);
    }
    if (ALIUtils::debug >= 5)
      std::cout << "sensor1_name " << sensor1_name << " sensor2_name " << sensor2_name << " sensor_dist " << sensor_dist
                << " unknown " << wl[3] << std::endl;

    ifdevi.getWordsInLine(wl);
    prec_deviX = atof(wl[0].c_str());
    prec_deviY = atof(wl[1].c_str());

    if (ALIUtils::debug >= 5)
      std::cout << "prec_deviX " << prec_deviX << " prec_deviY " << prec_deviY << std::endl;

    deviFromFile = new DeviationsFromFileSensor2D();
    ALIdouble offsetX, offsetY;
    if (wl.size() == 5) {
      offsetX = ALIUtils::getFloat(wl[3]);
      offsetY = ALIUtils::getFloat(wl[4]);
      deviFromFile->setOffset(offsetX, offsetY);
    }
    deviFromFile->readFile(ifdevi);
    fdevi_from_file = true;
    if (ALIUtils::debug >= 5)
      std::cout << "deviFromFile " << deviFromFile << std::endl;

    //--- Fill extra entries 'deviX' & 'deviY' to compute derivatives
    std::vector<ALIstring> wlo;
    char chartmp[20];
    wlo.push_back(wordlist[0]);
    wlo.push_back("deviX");
    wlo.push_back("0");  // value is set to 0 as it is on the file and the point of intersection is not computed yet
    gcvt(prec_deviX, 10, chartmp);
    wlo.push_back(ALIstring(chartmp));
    wlo.push_back("cal");
    std::vector<ALIstring> wl2(wlo);
    OpticalObject::fillExtraEntry(wlo);

    wl2[1] = "deviY";
    gcvt(prec_deviY, 10, chartmp);
    wl2[3] = ALIstring(chartmp);
    OpticalObject::fillExtraEntry(wl2);

  } else {
    OpticalObject::fillExtraEntry(wordlist);
  }
}

makeMeasurement()

void OptOSensor2D::makeMeasurement ( LightRay &  lightray, Measurement &  meas  )

overridevirtual

Reimplemented from OpticalObject.

Definition at line 40 of file OptOSensor2D.cc.

References funct::abs(), gather_cfg::cout, ALIUtils::debug, LightRay::direction(), ALIUtils::dump3v(), MillePedeFileConverter_cfg::e, LightRay::intersect(), Skims_PA_cff::name, class-composition::parent, LightRay::point(), Measurement::setLightRayDirection(), Measurement::setLightRayPosition(), Measurement::setValueSimulated(), Measurement::value(), Measurement::valueSimulated(), and Measurement::valueType().

                                                                        {
  //---------- check is last and not the only one
  /*    if( vocite - OptOList().begin() == 0 ||
    OptOList().end() - vocite != 1) {
    std::cerr << " last and not only one Optical Object should be 'sensor2D' (unless you specify (:T)raverse) " <<
    OptOList().end() - vocite << " size " <<OptOList().size() << std::endl;
    DumpBadOrderOptOs(); 
    }*/

  //---------- Get simulated value
  //---------- Get intersection
  CLHEP::Hep3Vector ZAxis(0., 0, 1.);
  CLHEP::HepRotation rmt = rmGlob();
  ZAxis = rmt * ZAxis;
  if (ALIUtils::debug >= 4) {
    ALIUtils::dump3v(centreGlob(), " sensor2D centre Glob ");
    ALIUtils::dump3v(ZAxis, " snsor2D normal ");
  }
  //-  ALIUtils::dumprm( rmt, "rot global " );
  lightray.intersect(ALIPlane(centreGlob(), ZAxis));
  CLHEP::Hep3Vector inters = lightray.point();

  ALIdouble* interslc;
  interslc = convertPointToLocalCoordinates(inters);
  ALIdouble interslcx = interslc[0];
  ALIdouble interslcy = interslc[1];
  meas.setValueSimulated(0, interslcx);
  meas.setValueSimulated(1, interslcy);

  //----- Dump info
  if (ALIUtils::debug >= 2) {
    //--- Special for range studies
    ALIstring chrg = "";
    /*t    if(Model::Ncmslinkrange >= 1 && Model::Ncmslinkrange <= 8 ) {
      chrg = "RG";
    } else {
      chrg = "";
      } t*/
    CLHEP::Hep3Vector measvv(meas.value()[0], meas.value()[1], 0.);
    measvv = rmt * measvv;
    ALIUtils::dump3v(measvv, " $$$$$$MEAS IN LOCAL FRAME");
    ALIUtils::dump3v(measvv + centreGlob(), " $$$$$$MEAS IN GLOBAL FRAME");

    ALIdouble detH = 1000 * meas.valueSimulated(0);
    if (std::abs(detH) <= 1.e-9)
      detH = 0.;
    ALIdouble detV = 1000 * meas.valueSimulated(1);
    if (std::abs(detV) <= 1.e-9)
      detV = 0.;
    std::cout << "REAL value: " << chrg << meas.valueType(0) << ": " << 1000 * meas.value()[0] << chrg << " "
              << meas.valueType(1) << ": " << 1000 * meas.value()[1] << " (mm)  " << (this)->name()
              << "   DIFF= " << detH - 1000 * meas.value()[0] << " " << detV - 1000 * meas.value()[1] << std::endl;
    std::cout << "SIMU value: " << chrg << " " << meas.valueType(0)
              << ": "
              // << setprecision(3) << setw(4)
              << detH << chrg << " " << meas.valueType(1) << ": " << detV << " (mm)  " << (this)->name() << std::endl;
    /*-    std::cout << "SIMU value: " << chrg << " " << meas.valueType(0) << ": "
      // << setprecision(3) << setw(4)
              << detH / 0.3125
          << chrg << " " << meas.valueType(1) << ": " << detV / 0.3125
          << " STRIPS  " << (this)->name() << std::endl; */
    //   << detH
    //   << chrg << " V: " << detV
    //   << " (mm)  " << (this)->name() << std::endl;
    ALIUtils::dump3v(1000. * (inters - parent()->centreGlob()), " $$$$$$SIMU inters - parent centre");
    ALIUtils::dump3v(1000. * (inters - centreGlob()), " $$$$$$SIMU inters - centre");
  }
  //t    delete &lightray;

  // store the lightray position and direction
  meas.setLightRayPosition(lightray.point());
  meas.setLightRayDirection(lightray.direction());

  delete[] interslc;
}

Member Data Documentation

deviFromFile

DeviationsFromFileSensor2D* OptOSensor2D::deviFromFile

private

Definition at line 50 of file OptOSensor2D.h.

fdevi_from_file

ALIbool OptOSensor2D::fdevi_from_file

private

Definition at line 51 of file OptOSensor2D.h.