#include <HitPatternHelper.h>

Public Member Functions
std::vector< int >	binary ()

std::vector< float >	bonusFeatures ()

int	etaSector ()

int	findLayer (int layerId)

	HitPatternHelper (const Setup *setup, int hitpattern, double cot, double z0)

int	num2S ()

int	numExpLayer ()

int	numMissing2S ()

int	numMissingInterior1 ()

int	numMissingInterior2 ()

int	numMissingPS ()

int	numPS ()

int	reducedId (int layerId)

	~HitPatternHelper ()

Private Attributes
std::vector< int >	binary_

std::vector< float >	bonusFeatures_

int	cotBin_

int	etaBin_

std::vector< double >	etaRegions_

int	etaSector_

int	hitpattern_

bool	hphDebug_

std::vector< int >	layerEncoding_

std::map< int, const tt::SensorModule * >	layerEncodingMap_

std::map< int, std::map< int, std::vector< int > > >	layermap_

int	nKalmanLayers_

int	num2S_

int	numExpLayer_

int	numMissing2S_

int	numMissingInterior1_

int	numMissingInterior2_

int	numMissingLayer_

int	numMissingPS_

int	numPS_

const Setup *	setup_

bool	useNewKF_

int	zTBin_

Detailed Description

Definition at line 89 of file HitPatternHelper.h.

Constructor & Destructor Documentation

◆ HitPatternHelper()

hph::HitPatternHelper::HitPatternHelper	(	const Setup *	setup,
		int	hitpattern,
		double	cot,
		double	z0
	)

Definition at line 58 of file HitPatternHelper.cc.

References binary_, relativeConstraints::empty, etaRegions_, etaSector_, findLayer(), hitpattern_, hphDebug_, mps_fire::i, createfilelist::int, dqmiolumiharvest::j, isotrackApplyRegressor::k, layerEncoding_, layerEncodingMap_, layermap_, nKalmanLayers_, num2S_, numExpLayer_, numMissing2S_, numMissingInterior1_, numMissingInterior2_, numMissingPS_, numPS_, reducedId(), cmsswSequenceInfo::seq, and useNewKF_.

       : setup_(setup),
         hphDebug_(setup_->hphDebug()),
         useNewKF_(setup_->useNewKF()),
         etaRegions_(setup_->etaRegions()),
         layermap_(setup_->layermap()),
         nKalmanLayers_(setup_->nKalmanLayers()),
         etaBin_(setup_->etaRegion(z0, cot, true)),
         cotBin_(setup_->digiCot(cot, etaBin_)),
         zTBin_(setup_->digiZT(z0, cot, etaBin_)),
         layerEncoding_(setup->layerEncoding(etaBin_, zTBin_, cotBin_)),
         layerEncodingMap_(setup->layerEncodingMap(etaBin_, zTBin_, cotBin_)),
         numExpLayer_(layerEncoding_.size()),
         hitpattern_(hitpattern),
         etaSector_(setup_->etaRegion(z0, cot, useNewKF_)),
         numMissingLayer_(0),
         numMissingPS_(0),
         numMissing2S_(0),
         numPS_(0),
         num2S_(0),
         numMissingInterior1_(0),
         numMissingInterior2_(0),
         binary_(11, 0),  //there are 11 unique layer IDs, as defined in variable "layerIds"
         bonusFeatures_() {
     int kf_eta_reg = etaSector_;
     if (kf_eta_reg < ((int)etaRegions_.size() - 1) / 2) {
       kf_eta_reg = ((int)etaRegions_.size() - 1) / 2 - 1 - kf_eta_reg;
     } else {
       kf_eta_reg = kf_eta_reg - (int)(etaRegions_.size() - 1) / 2;
     }
 
     int nbits = floor(log2(hitpattern_)) + 1;
     int lay_i = 0;
     bool seq = false;
     for (int i = 0; i < nbits; i++) {
       lay_i = ((1 << i) & hitpattern_) >> i;  //0 or 1 in ith bit (right to left)
 
       if (lay_i && !seq)
         seq = true;  //sequence starts when first 1 found
       if (!lay_i && seq) {
         numMissingInterior1_++;  //This is the same as the "tmp_trk_nlaymiss_interior" calculated in Trackquality.cc
       }
       if (!lay_i) {
         bool realhit = false;
         if (layermap_[kf_eta_reg][i].empty())
           continue;
         for (int j : layermap_[kf_eta_reg][i]) {
           int k = findLayer(j);
           if (k > 0)
             realhit = true;
         }
         if (realhit)
           numMissingInterior2_++;  //This variable doesn't make sense for new KF because it uses the layermap from Old KF
       }
     }
 
     if (hphDebug_) {
       if (useNewKF_) {
         edm::LogVerbatim("TrackTriggerHPH") << "Running with New KF";
       } else {
         edm::LogVerbatim("TrackTriggerHPH") << "Running with Old KF";
       }
       edm::LogVerbatim("TrackTriggerHPH") << "======================================================";
       edm::LogVerbatim("TrackTriggerHPH")
           << "Looking at hitpattern " << std::bitset<7>(hitpattern_) << "; Looping over KF layers:";
     }
 
     if (useNewKF_) {
       //New KF uses sensor modules to determine the hitmask already
       for (int i = 0; i < numExpLayer_; i++) {
         if (hphDebug_) {
           edm::LogVerbatim("TrackTriggerHPH") << "--------------------------";
           edm::LogVerbatim("TrackTriggerHPH") << "Looking at KF layer " << i;
           if (layerEncoding_[i] < 10) {
             edm::LogVerbatim("TrackTriggerHPH") << "KF expects L" << layerEncoding_[i];
           } else {
             edm::LogVerbatim("TrackTriggerHPH") << "KF expects D" << layerEncoding_[i] - 10;
           }
         }
 
         if (((1 << i) & hitpattern_) >> i) {
           if (hphDebug_) {
             edm::LogVerbatim("TrackTriggerHPH") << "Layer found in hitpattern";
           }
 
           binary_[reducedId(layerEncoding_[i])] = 1;
           if (layerEncodingMap_[layerEncoding_[i]]->psModule()) {
             numPS_++;
           } else {
             num2S_++;
           }
         } else {
           if (hphDebug_) {
             edm::LogVerbatim("TrackTriggerHPH") << "Layer missing in hitpattern";
           }
 
           if (layerEncodingMap_[layerEncoding_[i]]->psModule()) {
             numMissingPS_++;
           } else {
             numMissing2S_++;
           }
         }
       }
 
     } else {
       //Old KF uses the hard coded layermap to determien hitmask
       for (int i = 0; i < nKalmanLayers_; i++) {  //Loop over each digit of hitpattern
 
         if (hphDebug_) {
           edm::LogVerbatim("TrackTriggerHPH") << "--------------------------";
           edm::LogVerbatim("TrackTriggerHPH") << "Looking at KF layer " << i;
         }
 
         if (layermap_[kf_eta_reg][i].empty()) {
           if (hphDebug_) {
             edm::LogVerbatim("TrackTriggerHPH") << "KF does not expect this layer";
           }
 
           continue;
         }
 
         for (int j :
              layermap_[kf_eta_reg][i]) {  //Find out which layer the Old KF is dealing with when hitpattern is encoded
 
           if (hphDebug_) {
             if (j < 10) {
               edm::LogVerbatim("TrackTriggerHPH") << "KF expects L" << j;
             } else {
               edm::LogVerbatim("TrackTriggerHPH") << "KF expects D" << j - 10;
             }
           }
 
           int k = findLayer(j);
           if (k < 0) {
             //k<0 means even though layer j is predicted by Old KF, this prediction is rejected because it contradicts
             if (hphDebug_) {  //a more accurate prediction made with the help of information from sensor modules
               edm::LogVerbatim("TrackTriggerHPH") << "Rejected by sensor modules";
             }
 
             continue;
           }
 
           if (hphDebug_) {
             edm::LogVerbatim("TrackTriggerHPH") << "Confirmed by sensor modules";
           }
           //prediction is accepted
           if (((1 << i) & hitpattern_) >> i) {
             if (hphDebug_) {
               edm::LogVerbatim("TrackTriggerHPH") << "Layer found in hitpattern";
             }
 
             binary_[reducedId(j)] = 1;
             if (layerEncodingMap_[layerEncoding_[k]]->psModule()) {
               numPS_++;
             } else {
               num2S_++;
             }
           } else {
             if (hphDebug_) {
               edm::LogVerbatim("TrackTriggerHPH") << "Layer missing in hitpattern";
             }
 
             if (layerEncodingMap_[layerEncoding_[k]]->psModule()) {
               numMissingPS_++;
             } else {
               numMissing2S_++;
             }
           }
         }
       }
     }
 
     if (hphDebug_) {
       edm::LogVerbatim("TrackTriggerHPH") << "------------------------------";
       edm::LogVerbatim("TrackTriggerHPH") << "numPS = " << numPS_ << ", num2S = " << num2S_
                                           << ", missingPS = " << numMissingPS_ << ", missing2S = " << numMissing2S_;
       edm::LogVerbatim("TrackTriggerHPH") << "======================================================";
     }
   }

◆ ~HitPatternHelper()

hph::HitPatternHelper::~HitPatternHelper ( )

inline

Definition at line 92 of file HitPatternHelper.h.

92 {}

Member Function Documentation

◆ binary()

std::vector<int> hph::HitPatternHelper::binary ( )

inline

Definition at line 112 of file HitPatternHelper.h.

References binary_.

112 { return binary_; } //11-bit hitmask needed by TrackQuality.cc (0~5->L1~L6;6~10->D1~D5)

hph::HitPatternHelper::binary_

std::vector< int > binary_

Definition: HitPatternHelper.h:141

◆ bonusFeatures()

std::vector<float> hph::HitPatternHelper::bonusFeatures ( )

inline

Definition at line 113 of file HitPatternHelper.h.

References bonusFeatures_.

113 { return bonusFeatures_; } //bonus features for track quality

hph::HitPatternHelper::bonusFeatures_

std::vector< float > bonusFeatures_

Definition: HitPatternHelper.h:142

◆ etaSector()

int hph::HitPatternHelper::etaSector ( )

inline

Definition at line 94 of file HitPatternHelper.h.

References etaSector_.

94 { return etaSector_; } //Eta sectors defined in KF

hph::HitPatternHelper::etaSector_

int etaSector_

Definition: HitPatternHelper.h:133

◆ findLayer()

int hph::HitPatternHelper::findLayer ( int layerId )

Definition at line 285 of file HitPatternHelper.cc.

References mps_fire::i, createfilelist::int, and layerEncoding_.

Referenced by HitPatternHelper().

                                              {
     for (int i = 0; i < (int)layerEncoding_.size(); i++) {
       if (layerId == (int)layerEncoding_[i]) {
         return i;
       }
     }
     return -1;
   }

◆ num2S()

int hph::HitPatternHelper::num2S ( )

inline

Definition at line 105 of file HitPatternHelper.h.

References num2S_.

105 { return num2S_; } //The number of 2S layers are found in hitpattern

hph::HitPatternHelper::num2S_

int num2S_

Definition: HitPatternHelper.h:138

◆ numExpLayer()

int hph::HitPatternHelper::numExpLayer ( )

inline

Definition at line 95 of file HitPatternHelper.h.

References numExpLayer_.

95 { return numExpLayer_; } //The number of layers KF expects

hph::HitPatternHelper::numExpLayer_

int numExpLayer_

Definition: HitPatternHelper.h:131

◆ numMissing2S()

int hph::HitPatternHelper::numMissing2S ( )

inline

Definition at line 101 of file HitPatternHelper.h.

References numMissing2S_.

                        {
       return numMissing2S_;
     }  //The number of 2S layers that are missing. It includes the two types of layers mentioned above.

◆ numMissingInterior1()

int hph::HitPatternHelper::numMissingInterior1 ( )

inline

Definition at line 106 of file HitPatternHelper.h.

References numMissingInterior1_.

                               {
       return numMissingInterior1_;
     }  //The number of missing interior layers (using only hitpattern)

◆ numMissingInterior2()

int hph::HitPatternHelper::numMissingInterior2 ( )

inline

Definition at line 109 of file HitPatternHelper.h.

References numMissingInterior2_.

                               {
       return numMissingInterior2_;
     }  //The number of missing interior layers (using hitpattern, layermap from Old KF and sensor modules)

◆ numMissingPS()

int hph::HitPatternHelper::numMissingPS ( )

inline

Definition at line 96 of file HitPatternHelper.h.

References numMissingPS_.

                        {
       return numMissingPS_;
     }  //The number of PS layers that are missing. It includes layers that are missing:

◆ numPS()

int hph::HitPatternHelper::numPS ( )

inline

Definition at line 104 of file HitPatternHelper.h.

References numPS_.

104 { return numPS_; } //The number of PS layers are found in hitpattern

hph::HitPatternHelper::numPS_

int numPS_

Definition: HitPatternHelper.h:137

◆ reducedId()

int hph::HitPatternHelper::reducedId ( int layerId )

Definition at line 272 of file HitPatternHelper.cc.

References hphDebug_.

Referenced by HitPatternHelper().

                                              {
     if (hphDebug_ && (layerId > 15 || layerId < 1)) {
       edm::LogVerbatim("TrackTriggerHPH") << "Warning: invalid layer id !";
     }
     if (layerId <= 6) {
       layerId = layerId - 1;
       return layerId;
     } else {
       layerId = layerId - 5;
       return layerId;
     }
   }

Member Data Documentation

◆ binary_

std::vector<int> hph::HitPatternHelper::binary_

private

Definition at line 141 of file HitPatternHelper.h.

Referenced by binary(), and HitPatternHelper().

◆ bonusFeatures_

std::vector<float> hph::HitPatternHelper::bonusFeatures_

private

Definition at line 142 of file HitPatternHelper.h.

Referenced by bonusFeatures().

◆ cotBin_

int hph::HitPatternHelper::cotBin_

private

Definition at line 127 of file HitPatternHelper.h.

◆ etaBin_

int hph::HitPatternHelper::etaBin_

private

Definition at line 126 of file HitPatternHelper.h.

◆ etaRegions_

std::vector<double> hph::HitPatternHelper::etaRegions_

private

Definition at line 123 of file HitPatternHelper.h.

Referenced by HitPatternHelper().

◆ etaSector_

int hph::HitPatternHelper::etaSector_

private

Definition at line 133 of file HitPatternHelper.h.

Referenced by etaSector(), and HitPatternHelper().

◆ hitpattern_

int hph::HitPatternHelper::hitpattern_

private

Definition at line 132 of file HitPatternHelper.h.

Referenced by HitPatternHelper().

◆ hphDebug_

bool hph::HitPatternHelper::hphDebug_

private

Definition at line 121 of file HitPatternHelper.h.

Referenced by HitPatternHelper(), and reducedId().

◆ layerEncoding_

std::vector<int> hph::HitPatternHelper::layerEncoding_

private

Definition at line 129 of file HitPatternHelper.h.

Referenced by findLayer(), and HitPatternHelper().

◆ layerEncodingMap_

std::map<int, const tt::SensorModule*> hph::HitPatternHelper::layerEncodingMap_

private

Definition at line 130 of file HitPatternHelper.h.

Referenced by HitPatternHelper().

◆ layermap_

std::map<int, std::map<int, std::vector<int> > > hph::HitPatternHelper::layermap_

private

Definition at line 124 of file HitPatternHelper.h.

Referenced by HitPatternHelper().

◆ nKalmanLayers_

int hph::HitPatternHelper::nKalmanLayers_

private

Definition at line 125 of file HitPatternHelper.h.

Referenced by HitPatternHelper().

◆ num2S_

int hph::HitPatternHelper::num2S_

private

Definition at line 138 of file HitPatternHelper.h.

Referenced by HitPatternHelper(), and num2S().

◆ numExpLayer_

int hph::HitPatternHelper::numExpLayer_

private

Definition at line 131 of file HitPatternHelper.h.

Referenced by HitPatternHelper(), and numExpLayer().

◆ numMissing2S_

int hph::HitPatternHelper::numMissing2S_

private

Definition at line 136 of file HitPatternHelper.h.

Referenced by HitPatternHelper(), and numMissing2S().

◆ numMissingInterior1_

int hph::HitPatternHelper::numMissingInterior1_

private

Definition at line 139 of file HitPatternHelper.h.

Referenced by HitPatternHelper(), and numMissingInterior1().

◆ numMissingInterior2_

int hph::HitPatternHelper::numMissingInterior2_

private

Definition at line 140 of file HitPatternHelper.h.

Referenced by HitPatternHelper(), and numMissingInterior2().

◆ numMissingLayer_

int hph::HitPatternHelper::numMissingLayer_

private

Definition at line 134 of file HitPatternHelper.h.

◆ numMissingPS_

int hph::HitPatternHelper::numMissingPS_

private

Definition at line 135 of file HitPatternHelper.h.

Referenced by HitPatternHelper(), and numMissingPS().

◆ numPS_

int hph::HitPatternHelper::numPS_

private

Definition at line 137 of file HitPatternHelper.h.

Referenced by HitPatternHelper(), and numPS().

◆ setup_

const Setup* hph::HitPatternHelper::setup_

private

Definition at line 120 of file HitPatternHelper.h.

Referenced by upgradeWorkflowComponents.UpgradeWorkflow::setup().

◆ useNewKF_

bool hph::HitPatternHelper::useNewKF_

private

Definition at line 122 of file HitPatternHelper.h.

Referenced by HitPatternHelper().

◆ zTBin_

int hph::HitPatternHelper::zTBin_

private

Definition at line 128 of file HitPatternHelper.h.

Public Member Functions

Private Attributes

Detailed Description

Constructor & Destructor Documentation

◆ HitPatternHelper()

◆ ~HitPatternHelper()

Member Function Documentation

◆ binary()

◆ bonusFeatures()

◆ etaSector()

◆ findLayer()

◆ num2S()

◆ numExpLayer()

◆ numMissing2S()

◆ numMissingInterior1()

◆ numMissingInterior2()

◆ numMissingPS()

◆ numPS()

◆ reducedId()

Member Data Documentation

◆ binary_

◆ bonusFeatures_

◆ cotBin_

◆ etaBin_

◆ etaRegions_

◆ etaSector_

◆ hitpattern_

◆ hphDebug_

◆ layerEncoding_

◆ layerEncodingMap_

◆ layermap_

◆ nKalmanLayers_

◆ num2S_

◆ numExpLayer_

◆ numMissing2S_

◆ numMissingInterior1_

◆ numMissingInterior2_

◆ numMissingLayer_

◆ numMissingPS_

◆ numPS_

◆ setup_

◆ useNewKF_

◆ zTBin_