/afs/cern.ch/work/a/aaltunda/public/www/CMSSW_6_2_5/src/Alignment/TrackerAlignment/interface/TPBNameSpace.h

Go to the documentation of this file.

#ifndef Alignment_TrackerAlignment_TPBNameSpace_H
#define Alignment_TrackerAlignment_TPBNameSpace_H

#include "CondFormats/Alignment/interface/Definitions.h"
#include "DataFormats/TrackerCommon/interface/TrackerTopology.h"
#include "Geometry/Records/interface/IdealGeometryRecord.h"

namespace align
{
  namespace tpb
  {
    const unsigned int lpqc[] = {5, 8, 11};

    inline unsigned int moduleNumber(align::ID, const TrackerTopology*);

    inline unsigned int ladderNumber(align::ID, const TrackerTopology*);

    inline unsigned int layerNumber(align::ID, const TrackerTopology*);

    inline unsigned int halfBarrelNumber(align::ID, const TrackerTopology*);
  }
}

unsigned int align::tpb::moduleNumber(align::ID id, const TrackerTopology* tTopo)
{
  return tTopo->pxbModule(id);
}

unsigned int align::tpb::ladderNumber(align::ID id, const TrackerTopology* tTopo)
{
  

  unsigned int l = tTopo->pxbLadder(id); // increases with phi
  unsigned int c = tTopo->pxbLayer(id) - 1;

// Ladder in 1st quadrant: number = lpqc + 1 - l     (1 to lpqc)
// Ladder in 2nd quadrant: number = l - lpqc         (1 to lpqc)
// Ladder in 3rd quadrant: number = l - lpqc         (lpqc + 1 to 2 * lpqc)
// Ladder in 4th quadrant: number = 5 * lpqc + 1 - l (lpqc + 1 to 2 * lpqc)

  return l > 3 * lpqc[c] ? // ladder in 4th quadrant
    5 * lpqc[c] + 1 - l :
    (l > lpqc[c] ? // ladder not in 1st quadrant
     l - lpqc[c] : lpqc[c] + 1 - l);
}

unsigned int align::tpb::layerNumber(align::ID id, const TrackerTopology* tTopo)
{
  return tTopo->pxbLayer(id);
}

unsigned int align::tpb::halfBarrelNumber(align::ID id, const TrackerTopology* tTopo)
{
  

  unsigned int l = tTopo->pxbLadder(id); // increases with phi
  unsigned int c = tTopo->pxbLayer(id) - 1;

  return l > lpqc[c] && l <= 3 * lpqc[c] ? 1 : 2;
}

#endif