/data/refman/pasoursint/CMSSW_5_3_10/src/DataFormats/SiStripDetId/src/SiStripSubStructure.cc

Go to the documentation of this file.

// -*- C++ -*-
//
// Package:     SiStripDetId
// Class  :     SiStripSubStructure
// 
// Implementation:
//     <Notes on implementation>
//
// Original Author:  dkcira
//         Created:  Wed Jan 25 07:19:38 CET 2006
// $Id: SiStripSubStructure.cc,v 1.5 2007/03/16 09:57:09 dkcira Exp $
//

#include "DataFormats/SiStripDetId/interface/StripSubdetector.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"
#include "DataFormats/SiStripDetId/interface/SiStripSubStructure.h"

using namespace std;

SiStripSubStructure::SiStripSubStructure(){
}

SiStripSubStructure::~SiStripSubStructure(){
}


void SiStripSubStructure::getTIBDetectors(const std::vector<uint32_t> & inputDetRawIds,
                                          std::vector<uint32_t> & tibDetRawIds,
                                          uint32_t requested_layer,
                                          uint32_t requested_bkw_frw,
                                          uint32_t requested_int_ext,
                                          uint32_t requested_string) const{
 // loop over all input detectors
  for(vector<uint32_t>::const_iterator it = inputDetRawIds.begin(); it!=inputDetRawIds.end();it++){
    uint32_t therawid = (*it);                  // raw id of single detector
    TIBDetId potentialTIB = TIBDetId(therawid); // build TIBDetId, at this point is just DetId, but do not want to cast twice
    if( potentialTIB.subdetId() ==  int (StripSubdetector::TIB) ){ // check if subdetector field is a TIB, both tested numbers are int
      if( // check if TIB is from the ones requested
         (    (potentialTIB.layer()==requested_layer) || requested_layer==0 )  // take everything if default value is 0
         && ( ((potentialTIB.string()).at(0)==(requested_bkw_frw)) || requested_bkw_frw==0 )
         && ( ((potentialTIB.string()).at(1)==(requested_int_ext)) || requested_int_ext==0 )
         && ( ((potentialTIB.string()).at(2)==requested_string) || requested_string==0 )
         ){
        tibDetRawIds.push_back(therawid);       // add detector to list of selected TIBdets
      }
    }
  }
}


void SiStripSubStructure::getTIDDetectors(const std::vector<uint32_t> & inputDetRawIds,
                                          std::vector<uint32_t> & tidDetRawIds,
                                          uint32_t requested_side,
                                          uint32_t requested_wheel,
                                          uint32_t requested_ring,
                                          uint32_t requested_ster) const{
 // loop over all input detectors
  for(vector<uint32_t>::const_iterator it = inputDetRawIds.begin(); it!=inputDetRawIds.end();it++){
    uint32_t therawid = (*it);                  // raw id of single detector
    TIDDetId potentialTID = TIDDetId(therawid); // build TIDDetId, at this point is just DetId, but do not want to cast twice
    if( potentialTID.subdetId() ==  int (StripSubdetector::TID) ){ // check if subdetector field is a TID, both tested numbers are int
      if( // check if TID is from the ones requested
         (    (potentialTID.side()==requested_side) || requested_side==0 )  // take everything if default value is 0
         && ( (potentialTID.wheel()==requested_wheel) || requested_wheel==0 )
         && ( (potentialTID.ring()==requested_ring) || requested_ring==0 )
         && ( (potentialTID.stereo()==requested_ster) || requested_ster==0 )
         ){
        tidDetRawIds.push_back(therawid);       // add detector to list of selected TIDdets
      }
    }
  }
}


void SiStripSubStructure::getTOBDetectors(const std::vector<uint32_t> & inputDetRawIds,
                                          std::vector<uint32_t> & tobDetRawIds,
                                          uint32_t requested_layer,
                                          uint32_t requested_bkw_frw,
                                          uint32_t requested_rod) const{
 // loop over all input detectors
  for(vector<uint32_t>::const_iterator it = inputDetRawIds.begin(); it!=inputDetRawIds.end();it++){
    uint32_t therawid = (*it);                  // raw id of single detector
    TOBDetId potentialTOB = TOBDetId(therawid); // build TOBDetId, at this point is just DetId, but do not want to cast twice
    if( potentialTOB.subdetId() ==  int (StripSubdetector::TOB) ){ // check if subdetector field is a TOB, both tested numbers are int
      if( // check if TOB is from the ones requested
         (    (potentialTOB.layer()==requested_layer) || requested_layer==0 )  // take everything if default value is 0
         && ( ((potentialTOB.rod()).at(0)==(requested_bkw_frw)) || requested_bkw_frw==0 )
         && ( ((potentialTOB.rod()).at(1)==requested_rod) || requested_rod==0 )
         ){
        tobDetRawIds.push_back(therawid);       // add detector to list of selected TOBdets
      }
    }
  }
}


void SiStripSubStructure::getTECDetectors(const std::vector<uint32_t> & inputDetRawIds,
                                          std::vector<uint32_t> & tecDetRawIds,
                                          uint32_t requested_side,
                                          uint32_t requested_wheel,
                                          uint32_t requested_petal_bkw_frw,
                                          uint32_t requested_petal,
                                          uint32_t requested_ring,
                                          uint32_t requested_ster) const{ 
 // loop over all input detectors
  for(vector<uint32_t>::const_iterator it = inputDetRawIds.begin(); it!=inputDetRawIds.end();it++){
    uint32_t therawid = (*it);                  // raw id of single detector
    TECDetId potentialTEC = TECDetId(therawid); // build TECDetId, at this point is just DetId, but do not want to cast twice
    if( potentialTEC.subdetId() ==  int (StripSubdetector::TEC) ){ // check if subdetector field is a TEC, both tested numbers are int
      if( // check if TEC is from the ones requested
         (    (potentialTEC.side()==requested_side) || requested_side==0 )  // take everything if default value is 0
         && ( (potentialTEC.wheel()==requested_wheel) || requested_wheel==0 )
         && ( ((potentialTEC.petal()).at(0)==(requested_petal_bkw_frw-1)) || requested_petal_bkw_frw==0 )
         && ( ((potentialTEC.petal()).at(1)==requested_petal) || requested_petal==0 )
         && ( (potentialTEC.ring()==requested_ring) || requested_ring==0 )
         && ( (potentialTEC.stereo()==requested_ster) || requested_ster==0 )
         ){
        tecDetRawIds.push_back(therawid);       // add detector to list of selected TECdets
      }
    }
  }
}