AlignmentPI::TkAlBarycenters Struct Reference

#include <AlignmentPayloadInspectorHelper.h>

Public Member Functions
void	computeBarycenters (const std::vector< AlignTransform > &input, const TrackerTopology &tTopo, const std::map< AlignmentPI::coordinate, float > &GPR)
const double	getNModules (AlignmentPI::PARTITION p)
GlobalPoint	getPartitionAvg (AlignmentPI::PARTITION p)
const std::array< double, 6 >	getX ()
const std::array< double, 6 >	getY ()
const std::array< double, 6 >	getZ ()
void	init ()
virtual	~TkAlBarycenters ()

Public Attributes
std::map< AlignmentPI::PARTITION, double >	nmodules
std::map< AlignmentPI::PARTITION, double >	Xbarycenters
std::map< AlignmentPI::PARTITION, double >	Ybarycenters
std::map< AlignmentPI::PARTITION, double >	Zbarycenters

Detailed Description

Definition at line 823 of file AlignmentPayloadInspectorHelper.h.

Constructor & Destructor Documentation

~TkAlBarycenters()

virtual AlignmentPI::TkAlBarycenters::~TkAlBarycenters ( )

inlinevirtual

Definition at line 876 of file AlignmentPayloadInspectorHelper.h.

{}

Member Function Documentation

computeBarycenters()

void AlignmentPI::TkAlBarycenters::computeBarycenters ( const std::vector< AlignTransform > &  input, const TrackerTopology &  tTopo, const std::map< AlignmentPI::coordinate, float > &  GPR  )

inline

Definition at line 887 of file AlignmentPayloadInspectorHelper.h.

References cms::cuda::assert(), AlignmentPI::BPIX, COUT, AlignmentPI::FPIXm, AlignmentPI::FPIXp, input, nmodules, AlCaHLTBitMon_ParallelJobs::p, AlignmentPI::PARTITIONS, PixelSubdetector::PixelBarrel, PixelSubdetector::PixelEndcap, TrackerTopology::pxfSide(), AlignmentPI::t_x, AlignmentPI::t_y, AlignmentPI::t_z, StripSubdetector::TEC, AlignmentPI::TECm, AlignmentPI::TECp, TrackerTopology::tecSide(), StripSubdetector::TIB, AlignmentPI::TIB, StripSubdetector::TID, AlignmentPI::TIDm, AlignmentPI::TIDp, TrackerTopology::tidSide(), StripSubdetector::TOB, AlignmentPI::TOB, DetId::Tracker, Xbarycenters, Ybarycenters, and Zbarycenters.

  {
    // zero in the n. modules per partition...
    for (const auto& p : PARTITIONS) {
      nmodules[p] = 0.;
    }

    for (const auto& ali : input) {
      if (DetId(ali.rawId()).det() != DetId::Tracker) {
        edm::LogWarning("TkAlBarycenters::computeBarycenters")
            << "Encountered invalid Tracker DetId:" << ali.rawId() << " " << DetId(ali.rawId()).det()
            << " is different from " << DetId::Tracker << "  - terminating ";
        assert(DetId(ali.rawId()).det() != DetId::Tracker);
      }

      int subid = DetId(ali.rawId()).subdetId();
      switch (subid) {
        case PixelSubdetector::PixelBarrel:
          Xbarycenters[PARTITION::BPIX] += (ali.translation().x());
          Ybarycenters[PARTITION::BPIX] += (ali.translation().y());
          Zbarycenters[PARTITION::BPIX] += (ali.translation().z());
          nmodules[PARTITION::BPIX]++;
          break;
        case PixelSubdetector::PixelEndcap:

          // minus side
          if (tTopo.pxfSide(DetId(ali.rawId())) == 1) {
            Xbarycenters[PARTITION::FPIXm] += (ali.translation().x());
            Ybarycenters[PARTITION::FPIXm] += (ali.translation().y());
            Zbarycenters[PARTITION::FPIXm] += (ali.translation().z());
            nmodules[PARTITION::FPIXm]++;
          }  // plus side
          else {
            Xbarycenters[PARTITION::FPIXp] += (ali.translation().x());
            Ybarycenters[PARTITION::FPIXp] += (ali.translation().y());
            Zbarycenters[PARTITION::FPIXp] += (ali.translation().z());
            nmodules[PARTITION::FPIXp]++;
          }
          break;
        case StripSubdetector::TIB:
          Xbarycenters[PARTITION::TIB] += (ali.translation().x());
          Ybarycenters[PARTITION::TIB] += (ali.translation().y());
          Zbarycenters[PARTITION::TIB] += (ali.translation().z());
          nmodules[PARTITION::TIB]++;
          break;
        case StripSubdetector::TID:
          // minus side
          if (tTopo.tidSide(DetId(ali.rawId())) == 1) {
            Xbarycenters[PARTITION::TIDm] += (ali.translation().x());
            Ybarycenters[PARTITION::TIDm] += (ali.translation().y());
            Zbarycenters[PARTITION::TIDm] += (ali.translation().z());
            nmodules[PARTITION::TIDm]++;
          }  // plus side
          else {
            Xbarycenters[PARTITION::TIDp] += (ali.translation().x());
            Ybarycenters[PARTITION::TIDp] += (ali.translation().y());
            Zbarycenters[PARTITION::TIDp] += (ali.translation().z());
            nmodules[PARTITION::TIDp]++;
          }
          break;
        case StripSubdetector::TOB:
          Xbarycenters[PARTITION::TOB] += (ali.translation().x());
          Ybarycenters[PARTITION::TOB] += (ali.translation().y());
          Zbarycenters[PARTITION::TOB] += (ali.translation().z());
          nmodules[PARTITION::TOB]++;
          break;
        case StripSubdetector::TEC:
          // minus side
          if (tTopo.tecSide(DetId(ali.rawId())) == 1) {
            Xbarycenters[PARTITION::TECm] += (ali.translation().x());
            Ybarycenters[PARTITION::TECm] += (ali.translation().y());
            Zbarycenters[PARTITION::TECm] += (ali.translation().z());
            nmodules[PARTITION::TECm]++;
          }  // plus side
          else {
            Xbarycenters[PARTITION::TECp] += (ali.translation().x());
            Ybarycenters[PARTITION::TECp] += (ali.translation().y());
            Zbarycenters[PARTITION::TECp] += (ali.translation().z());
            nmodules[PARTITION::TECp]++;
          }
          break;
        default:
          edm::LogError("TrackerAlignment_PayloadInspector") << "Unrecognized partition " << subid << std::endl;
          break;
      }
    }

    for (const auto& p : PARTITIONS) {
      // take the arithmetic mean
      Xbarycenters[p] /= nmodules[p];
      Ybarycenters[p] /= nmodules[p];
      Zbarycenters[p] /= nmodules[p];

      // add the Tracker Global Position Record
      Xbarycenters[p] += GPR.at(AlignmentPI::t_x);
      Ybarycenters[p] += GPR.at(AlignmentPI::t_y);
      Zbarycenters[p] += GPR.at(AlignmentPI::t_z);

      COUT << "Partition: " << p << " n. modules: " << nmodules[p] << "|"
           << " X: " << std::right << std::setw(12) << Xbarycenters[p] << " Y: " << std::right << std::setw(12)
           << Ybarycenters[p] << " Z: " << std::right << std::setw(12) << Zbarycenters[p] << std::endl;
    }
  }

getNModules()

const double AlignmentPI::TkAlBarycenters::getNModules ( AlignmentPI::PARTITION  p )

inline

Definition at line 835 of file AlignmentPayloadInspectorHelper.h.

References nmodules, and AlCaHLTBitMon_ParallelJobs::p.

{ return nmodules[p]; };

getPartitionAvg()

GlobalPoint AlignmentPI::TkAlBarycenters::getPartitionAvg ( AlignmentPI::PARTITION  p )

inline

Definition at line 880 of file AlignmentPayloadInspectorHelper.h.

References AlCaHLTBitMon_ParallelJobs::p, Xbarycenters, Ybarycenters, and Zbarycenters.

  {
    return GlobalPoint(Xbarycenters[p], Ybarycenters[p], Zbarycenters[p]);
  }

getX()

const std::array<double, 6> AlignmentPI::TkAlBarycenters::getX ( )

inline

Definition at line 842 of file AlignmentPayloadInspectorHelper.h.

References AlignmentPI::BPIX, AlignmentPI::FPIXm, AlignmentPI::FPIXp, AlignmentPI::TECm, AlignmentPI::TECp, AlignmentPI::TIB, AlignmentPI::TIDm, AlignmentPI::TIDp, AlignmentPI::TOB, and Xbarycenters.

    {
      return {{Xbarycenters[PARTITION::BPIX],
               (Xbarycenters[PARTITION::FPIXm] + Xbarycenters[PARTITION::FPIXp]) / 2.,
               Xbarycenters[PARTITION::TIB],
               (Xbarycenters[PARTITION::TIDm] + Xbarycenters[PARTITION::TIDp]) / 2,
               Xbarycenters[PARTITION::TOB],
               (Xbarycenters[PARTITION::TECm] + Xbarycenters[PARTITION::TECp]) / 2}};
    };

getY()

const std::array<double, 6> AlignmentPI::TkAlBarycenters::getY ( )

inline

Definition at line 854 of file AlignmentPayloadInspectorHelper.h.

References AlignmentPI::BPIX, AlignmentPI::FPIXm, AlignmentPI::FPIXp, AlignmentPI::TECm, AlignmentPI::TECp, AlignmentPI::TIB, AlignmentPI::TIDm, AlignmentPI::TIDp, AlignmentPI::TOB, and Ybarycenters.

    {
      return {{Ybarycenters[PARTITION::BPIX],
               (Ybarycenters[PARTITION::FPIXm] + Ybarycenters[PARTITION::FPIXp]) / 2.,
               Ybarycenters[PARTITION::TIB],
               (Ybarycenters[PARTITION::TIDm] + Ybarycenters[PARTITION::TIDp]) / 2,
               Ybarycenters[PARTITION::TOB],
               (Ybarycenters[PARTITION::TECm] + Ybarycenters[PARTITION::TECp]) / 2}};
    };

getZ()

const std::array<double, 6> AlignmentPI::TkAlBarycenters::getZ ( )

inline

Definition at line 866 of file AlignmentPayloadInspectorHelper.h.

References AlignmentPI::BPIX, AlignmentPI::FPIXm, AlignmentPI::FPIXp, AlignmentPI::TECm, AlignmentPI::TECp, AlignmentPI::TIB, AlignmentPI::TIDm, AlignmentPI::TIDp, AlignmentPI::TOB, and Zbarycenters.

    {
      return {{Zbarycenters[PARTITION::BPIX],
               (Zbarycenters[PARTITION::FPIXm] + Zbarycenters[PARTITION::FPIXp]) / 2.,
               Zbarycenters[PARTITION::TIB],
               (Zbarycenters[PARTITION::TIDm] + Zbarycenters[PARTITION::TIDp]) / 2,
               Zbarycenters[PARTITION::TOB],
               (Zbarycenters[PARTITION::TECm] + Zbarycenters[PARTITION::TECp]) / 2}};
    };

init()

void AlignmentPI::TkAlBarycenters::init

(

)

Member Data Documentation

nmodules

std::map<AlignmentPI::PARTITION, double> AlignmentPI::TkAlBarycenters::nmodules

Definition at line 827 of file AlignmentPayloadInspectorHelper.h.

Referenced by computeBarycenters(), and getNModules().

Xbarycenters

std::map<AlignmentPI::PARTITION, double> AlignmentPI::TkAlBarycenters::Xbarycenters

Definition at line 824 of file AlignmentPayloadInspectorHelper.h.

Referenced by computeBarycenters(), getPartitionAvg(), and getX().

Ybarycenters

std::map<AlignmentPI::PARTITION, double> AlignmentPI::TkAlBarycenters::Ybarycenters

Definition at line 825 of file AlignmentPayloadInspectorHelper.h.

Referenced by computeBarycenters(), getPartitionAvg(), and getY().

Zbarycenters

std::map<AlignmentPI::PARTITION, double> AlignmentPI::TkAlBarycenters::Zbarycenters

Definition at line 826 of file AlignmentPayloadInspectorHelper.h.

Referenced by computeBarycenters(), getPartitionAvg(), and getZ().