TPENameSpace.h

Go to the documentation of this file.

#ifndef Alignment_CommonAlignment_TPENameSpace_H
#define Alignment_CommonAlignment_TPENameSpace_H

#include "CondFormats/Alignment/interface/Definitions.h"
#include "DataFormats/SiPixelDetId/interface/PXFDetId.h"

namespace align
{
  namespace tpe
  {
    const unsigned int bpqd = 6; // no. of blades per quarter disk

    inline unsigned int moduleNumber( align::ID );

    inline unsigned int panelNumber( align::ID );

    inline unsigned int bladeNumber( align::ID );

    inline unsigned int halfDiskNumber( align::ID );

    inline unsigned int halfCylinderNumber( align::ID );

    inline unsigned int endcapNumber( align::ID );
  }
}

unsigned int align::tpe::moduleNumber(align::ID id)
{
  return PXFDetId(id).module();
}

unsigned int align::tpe::panelNumber(align::ID id)
{
  return PXFDetId(id).panel();
}

unsigned int align::tpe::bladeNumber(align::ID id)
{
  unsigned int b = PXFDetId(id).blade(); // 1 to 24 in increasing phi

// Blade in 1st quadrant: number = bpqd + 1 - b     (1 to bpqd)
// Blade in 2nd quadrant: number = b - bpqd         (1 to bpqd)
// Blade in 3rd quadrant: number = b - bpqd         (bpqd + 1 to 2 * bpqd)
// Blade in 4th quadrant: number = 5 * bpqd + 1 - b (bpqd + 1 to 2 * bpqd)

  return b > 3 * bpqd ? // blade in 4th quadrant
    5 * bpqd + 1 - b :
    (b > bpqd ? // blade not in 1st quadrant
     b - bpqd : bpqd + 1 - b);
}

unsigned int align::tpe::halfDiskNumber(align::ID id)
{
  return PXFDetId(id).disk();
}

unsigned int align::tpe::halfCylinderNumber(align::ID id)
{
  unsigned int b = PXFDetId(id).blade(); // 1 to 24 in increasing phi

  return b > bpqd && b <= 3 * bpqd ? 1 : 2;
}

unsigned int align::tpe::endcapNumber(align::ID id)
{
  return PXFDetId(id).side();
}

#endif