#include <HitPatternHelper.h>

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

int	etaSector ()

int	findLayer (int layerId)

	HitPatternHelper ()

	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 ()

Static Public Member Functions
static auto	equalID (tt::SensorModule lhs, tt::SensorModule rhs)

static auto	smallerID (tt::SensorModule lhs, tt::SensorModule rhs)

Private Attributes
std::vector< int >	binary_

float	chosenRofZ_

double	cot_

float	deltaTanL_

std::vector< double >	etaRegions_

int	etaSector_

int	hitpattern_

bool	hphDebug_

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

std::vector< tt::SensorModule >	layers_

int	nKalmanLayers_

int	num2S_

int	numExpLayer_

int	numMissing2S_

int	numMissingInterior1_

int	numMissingInterior2_

int	numMissingLayer_

int	numMissingPS_

int	numPS_

const Setup *	setup_

bool	useNewKF_

double	z0_

Detailed Description

Definition at line 72 of file HitPatternHelper.h.

Constructor & Destructor Documentation

◆ HitPatternHelper() [1/2]

hph::HitPatternHelper::HitPatternHelper ( )

inline

Definition at line 74 of file HitPatternHelper.h.

74 {}

◆ HitPatternHelper() [2/2]

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

Definition at line 44 of file HitPatternHelper.cc.

References funct::abs(), binary_, chosenRofZ_, cot_, ztail::d, deltaTanL_, relativeConstraints::empty, equalID(), ALCARECOTkAlBeamHalo_cff::etaMax, etaRegions_, etaSector_, JetChargeProducer_cfi::exp, findLayer(), hitpattern_, hphDebug_, mps_fire::i, createfilelist::int, dqmiolumiharvest::j, dqmdumpme::k, layermap_, layers_, nKalmanLayers_, num2S_, numExpLayer_, numMissing2S_, numMissingInterior1_, numMissingInterior2_, numMissingPS_, numPS_, ReducedId(), hph::Setup::sensorModules(), cmsswSequenceInfo::seq, setup_, smallerID(), jetUpdater_cfi::sort, funct::tan(), tier0::unique(), useNewKF_, and z0_.

       : hitpattern_(hitpattern),
         numExpLayer_(0),
         numMissingLayer_(0),
         numMissingPS_(0),
         numMissing2S_(0),
         numPS_(0),
         num2S_(0),
         numMissingInterior1_(0),
         numMissingInterior2_(0),
         cot_(cot),
         z0_(z0),
         setup_(setup),
         binary_(11, 0),
         hphDebug_(setup_->hphDebug()),
         useNewKF_(setup_->useNewKF()),
         chosenRofZ_(setup_->chosenRofZ()),
         deltaTanL_(setup_->deltaTanL()),
         etaRegions_(setup_->etaRegions()),
         nKalmanLayers_(setup_->nKalmanLayers()),
         layermap_(setup_->layermap()) {
     //Calculating eta sector based on cot and z0
     float kfzRef = z0_ + chosenRofZ_ * cot_;
     int kf_eta_reg = 0;
     for (int iEtaSec = 1; iEtaSec < ((int)etaRegions_.size() - 1); iEtaSec++) {  // Doesn't apply eta < 2.4 cut.
       float etaMax = etaRegions_[iEtaSec];
       float zRefMax = chosenRofZ_ / tan(2. * atan(exp(-etaMax)));
       if (kfzRef > zRefMax) {
         kf_eta_reg = iEtaSec;
       }
     }
     etaSector_ = kf_eta_reg;
     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;
     }
     //Looping over sensor modules to make predictions on which layers particles are expected to hit
     for (const tt::SensorModule& sm : setup_->sensorModules()) {
       double d = (z0_ - sm.z() + sm.r() * cot_) / (sm.cosTilt() - sm.sinTilt() * cot_);
       double d_p =
           (z0_ - sm.z() + sm.r() * (cot_ + deltaTanL_ / 2)) / (sm.cosTilt() - sm.sinTilt() * (cot_ + deltaTanL_ / 2));
       double d_m =
           (z0_ - sm.z() + sm.r() * (cot_ - deltaTanL_ / 2)) / (sm.cosTilt() - sm.sinTilt() * (cot_ - deltaTanL_ / 2));
       if (useNewKF_ &&
           (abs(d_p) < sm.numColumns() * sm.pitchCol() / 2. || abs(d_m) < sm.numColumns() * sm.pitchCol() / 2.)) {
         layers_.push_back(sm);
       }
       if (!useNewKF_ && abs(d) < sm.numColumns() * sm.pitchCol() / 2.) {
         layers_.push_back(sm);
       }
     }
     //layers_ constains all the sensor modules that particles are expected to hit
     sort(layers_.begin(), layers_.end(), smallerID);
     layers_.erase(unique(layers_.begin(), layers_.end(), equalID), layers_.end());  //Keep only one sensor per layer
 
     numExpLayer_ = layers_.size();
 
     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 (layers_[i].layerId() < 10) {
             edm::LogVerbatim("TrackTriggerHPH") << "KF expects L" << layers_[i].layerId();
           } else {
             edm::LogVerbatim("TrackTriggerHPH") << "KF expects D" << layers_[i].layerId() - 10;
           }
         }
 
         if (((1 << i) & hitpattern_) >> i) {
           if (hphDebug_) {
             edm::LogVerbatim("TrackTriggerHPH") << "Layer found in hitpattern";
           }
 
           binary_[ReducedId(layers_[i].layerId())] = 1;
           if (layers_[i].psModule()) {
             numPS_++;
           } else {
             num2S_++;
           }
         } else {
           if (hphDebug_) {
             edm::LogVerbatim("TrackTriggerHPH") << "Layer missing in hitpattern";
           }
 
           if (layers_[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 (layers_[k].psModule()) {
               numPS_++;
             } else {
               num2S_++;
             }
           } else {
             if (hphDebug_) {
               edm::LogVerbatim("TrackTriggerHPH") << "Layer missing in hitpattern";
             }
 
             if (layers_[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 76 of file HitPatternHelper.h.

76 {}

Member Function Documentation

◆ binary()

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

inline

Definition at line 96 of file HitPatternHelper.h.

References binary_.

96 { 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:119

◆ equalID()

static auto hph::HitPatternHelper::equalID	(	tt::SensorModule	lhs,
		tt::SensorModule	rhs
	)

inlinestatic

Definition at line 98 of file HitPatternHelper.h.

References tt::SensorModule::layerId().

Referenced by HitPatternHelper().

98 { return lhs.layerId() == rhs.layerId(); }

tt::SensorModule::layerId

int layerId() const

Definition: SensorModule.h:43

◆ etaSector()

int hph::HitPatternHelper::etaSector ( )

inline

Definition at line 78 of file HitPatternHelper.h.

References etaSector_.

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

hph::HitPatternHelper::etaSector_

int etaSector_

Definition: HitPatternHelper.h:105

◆ findLayer()

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

Definition at line 264 of file HitPatternHelper.cc.

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

Referenced by HitPatternHelper().

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

◆ num2S()

int hph::HitPatternHelper::num2S ( )

inline

Definition at line 89 of file HitPatternHelper.h.

References num2S_.

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

hph::HitPatternHelper::num2S_

int num2S_

Definition: HitPatternHelper.h:112

◆ numExpLayer()

int hph::HitPatternHelper::numExpLayer ( )

inline

Definition at line 79 of file HitPatternHelper.h.

References numExpLayer_.

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

hph::HitPatternHelper::numExpLayer_

int numExpLayer_

Definition: HitPatternHelper.h:107

◆ numMissing2S()

int hph::HitPatternHelper::numMissing2S ( )

inline

Definition at line 85 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 90 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 93 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 80 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 88 of file HitPatternHelper.h.

References numPS_.

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

hph::HitPatternHelper::numPS_

int numPS_

Definition: HitPatternHelper.h:111

◆ ReducedId()

int hph::HitPatternHelper::ReducedId ( int layerId )

Definition at line 251 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;
     }
   };