Line Code
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60 61 62 63 64 65 66 67 68 69 70 71 72 73 74 75 76 77 78 79 80 81 82 83 84 85 86 87 88 89 90 91 92 93 94 95 96 97
//   COCOA class implementation file
//Id:  EntryLengthAffCentre.C
//CAT: Model
//
//   History: v1.0
//   Pedro Arce

#include "Alignment/CocoaModel/interface/EntryLengthAffCentre.h"
#include "Alignment/CocoaModel/interface/OpticalObject.h"
#include "Alignment/CocoaUtilities/interface/ALIUtils.h"

//@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@
EntryLengthAffCentre::EntryLengthAffCentre(const ALIstring& type) : EntryLength(type) {}

//@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@
void EntryLengthAffCentre::FillName(const ALIstring& name) {
  ALIstring nn = "Centre ";
  nn += name;
  setName(nn);
}

//@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@
//@@
//@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@
void EntryLengthAffCentre::displace(ALIdouble disp) {
  if (ALIUtils::debug >= 9)
    std::cout << "EntryLengthAffCentre::Displace" << disp << std::endl;
  ALIint namelength = name().length() - 1;
  XYZcoor axisNo = XCoor;
  if (name_[namelength] == 'X') {
    axisNo = XCoor;
  } else if (name_[namelength] == 'Y') {
    axisNo = YCoor;
  } else if (name_[namelength] == 'Z') {
    axisNo = ZCoor;
  }
  OptOCurrent()->displaceCentreGlob(axisNo, disp);
}

//@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@
void EntryLengthAffCentre::displaceOriginal(ALIdouble disp) {
  if (ALIUtils::debug >= 9)
    std::cout << "EntryLengthAffCentre::DisplaceOriginal" << disp << std::endl;
  ALIint namelength = name().length() - 1;
  if (name_[namelength] == 'X') {
    OptOCurrent()->displaceCentreGlobOriginal(XCoor, disp);
  } else if (name_[namelength] == 'Y') {
    OptOCurrent()->displaceCentreGlobOriginal(YCoor, disp);
  } else if (name_[namelength] == 'Z') {
    OptOCurrent()->displaceCentreGlobOriginal(ZCoor, disp);
  }
}

//@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@
void EntryLengthAffCentre::displaceOriginalOriginal(ALIdouble disp) {
  if (ALIUtils::debug >= 9)
    std::cout << "EntryLengthAffCentre::DisplaceOriginalOriginal" << disp << std::endl;
  ALIint namelength = name().length() - 1;
  if (name_[namelength] == 'X') {
    OptOCurrent()->displaceCentreGlobOriginalOriginal(XCoor, disp);
  } else if (name_[namelength] == 'Y') {
    OptOCurrent()->displaceCentreGlobOriginalOriginal(YCoor, disp);
  } else if (name_[namelength] == 'Z') {
    OptOCurrent()->displaceCentreGlobOriginalOriginal(ZCoor, disp);
  }
}

//@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@
ALIdouble EntryLengthAffCentre::valueInGlobalReferenceFrame() const {
  ALIdouble a = 2.;

  return a;
}

//@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@
ALIdouble EntryLengthAffCentre::valueDisplaced() const {
  ALIdouble vdisp = 0.;

  CLHEP::Hep3Vector cdisp = OptOCurrent()->centreGlob() - OptOCurrent()->centreGlobOriginal();
  CLHEP::HepRotation rmParentInv = inverseOf(OptOCurrent()->parent()->rmGlob());
  cdisp = rmParentInv * cdisp;

  if (name() == "centre_X") {
    return cdisp.x();
  } else if (name() == "centre_Y") {
    return cdisp.y();
    //-   return OptOCurrent()->centreLocal().y() - value();
  } else if (name() == "centre_Z") {
    return cdisp.z();
  }

  if (ALIUtils::debug >= 5)
    std::cout << name() << " in OptO " << OptOCurrent()->name() << " valueDisplaced: " << vdisp << std::endl;

  return 0.;  // to avoid warning
}