DeDxHit.cc

Go to the documentation of this file.

#include "DataFormats/TrackReco/interface/DeDxHit.h"
#include "DataFormats/SiPixelDetId/interface/PixelSubdetector.h"
#include "DataFormats/SiStripDetId/interface/StripSubdetector.h"
#include "DataFormats/SiPixelDetId/interface/PXBDetId.h"
#include "DataFormats/SiPixelDetId/interface/PXFDetId.h"
#include "DataFormats/SiStripDetId/interface/TIBDetId.h"
#include "DataFormats/SiStripDetId/interface/TIDDetId.h"
#include "DataFormats/SiStripDetId/interface/TOBDetId.h"
#include "DataFormats/SiStripDetId/interface/TECDetId.h"
using namespace reco;

DeDxHit::DeDxHit(float ch, float mom, float len, DetId id):
    m_charge(ch),
    m_momentum(mom),
    m_pathLength(len)

{
    uint32_t subdet = id.subdetId();
    m_subDetId = (subdet & 0x7) << 5;

    uint32_t layer = 0;
    if (subdet == PixelSubdetector::PixelBarrel) {
        layer = PXBDetId(id).layer();
    } else if (subdet == PixelSubdetector::PixelEndcap) {
        layer = PXFDetId(id).disk();
    } else if (subdet == StripSubdetector::TIB) {
        layer = TIBDetId(id).layer();
    } else if (subdet == StripSubdetector::TID) {
        layer = (TIDDetId(id).wheel() & 0xF) + ((TIDDetId(id).side() - 1) << 4);
    } else if (subdet == StripSubdetector::TOB) {
        layer = TOBDetId(id).layer();
    } else if (subdet == StripSubdetector::TEC) {
        layer = (TECDetId(id).wheel() & 0xF) + ((TECDetId(id).side() - 1) << 4);
    }

    m_subDetId += (layer & 0x1F) ;
}