SiStripTKNumbers.cc

Go to the documentation of this file.

//#include <map>

#include "DataFormats/DetId/interface/DetId.h"
#include "DataFormats/SiStripDetId/interface/StripSubdetector.h"
#include "DPGAnalysis/SiStripTools/interface/SiStripTKNumbers.h"

SiStripTKNumbers::SiStripTKNumbers() {

  DetId tk(DetId::Tracker,0);
  _nmodules[tk.rawId()] = (3540-816)+816+5208+6400;
  _nfibres[tk.rawId()] = (9192-2208)+2208+12906+15104;
  _nmodules[tk.subdetId()] = _nmodules[tk.rawId()];
  _nfibres[tk.subdetId()] = _nfibres[tk.rawId()];

  DetId tib(DetId::Tracker, StripSubdetector::TIB);
  _nmodules[tib.rawId()] = 3540-816;
  _nfibres[tib.rawId()] = 9192-2208;
  _nmodules[tib.subdetId()] = _nmodules[tib.rawId()];
  _nfibres[tib.subdetId()] = _nfibres[tib.rawId()];

  DetId tid(DetId::Tracker, StripSubdetector::TID);
  _nmodules[tid.rawId()] = 816;
  _nfibres[tid.rawId()] = 2208;
  _nmodules[tid.subdetId()] = _nmodules[tid.rawId()];
  _nfibres[tid.subdetId()] = _nfibres[tid.rawId()];

  DetId tob(DetId::Tracker, StripSubdetector::TOB);
  _nmodules[tob.rawId()] = 5208;
  _nfibres[tob.rawId()] = 12906;
  _nmodules[tob.subdetId()] = _nmodules[tob.rawId()];
  _nfibres[tob.subdetId()] = _nfibres[tob.rawId()];

  DetId tec(DetId::Tracker, StripSubdetector::TEC);
  _nmodules[tec.rawId()] = 6400;
  _nfibres[tec.rawId()] = 15104;
  _nmodules[tec.subdetId()] = _nmodules[tec.rawId()];
  _nfibres[tec.subdetId()] = _nfibres[tec.rawId()];

  DetId tecp(DetId(DetId::Tracker, StripSubdetector::TEC).rawId() | ( (1& 0x3) << 18)); 
  _nmodules[tecp.rawId()] = 3200;
  _nfibres[tecp.rawId()] = 7552;

  DetId tecm(DetId(DetId::Tracker, StripSubdetector::TEC).rawId() | ( (2& 0x3) << 18)); 
  _nmodules[tecm.rawId()] = 3200;
  _nfibres[tecm.rawId()] = 7552;

}
  
int SiStripTKNumbers::nmodules(const DetId& detid) const {

  int subd = detid.subdetId();
  if(_nmodules.find(subd)!=_nmodules.end()) return _nmodules.find(subd)->second;

  return 0;
}

int SiStripTKNumbers::nmodules(const int id) const {

  if(_nmodules.find(id)!=_nmodules.end()) return _nmodules.find(id)->second;

  return 0;
}

int SiStripTKNumbers::nfibres(const DetId& detid) const {

  int subd = detid.subdetId();
  if(_nfibres.find(subd)!=_nfibres.end()) return _nfibres.find(subd)->second;

  return 0;
}

int SiStripTKNumbers::nfibres(const int id) const {

  if(_nfibres.find(id)!=_nfibres.end()) return _nfibres.find(id)->second;

  return 0;
}

int SiStripTKNumbers::napvs(const DetId& detid) const {
  
  return nfibres(detid)*_apvsperfibre;
}

int SiStripTKNumbers::napvs(const int id) const {

  return nfibres(id)*_apvsperfibre;
}

int SiStripTKNumbers::nstrips(const DetId& detid) const {
  
  return nfibres(detid)*_apvsperfibre*_stripsperapv;
}

int SiStripTKNumbers::nstrips(const int id) const {

  return nfibres(id)*_apvsperfibre*_stripsperapv;
}