StubWindowSuggest.cc

Go to the documentation of this file.

#include "L1Trigger/TrackFindingTMTT/interface/StubWindowSuggest.h"
#include "L1Trigger/TrackFindingTMTT/interface/Stub.h"
#include "L1Trigger/TrackFindingTMTT/interface/TrackerModule.h"
#include "L1Trigger/TrackFindingTMTT/interface/PrintL1trk.h"
 
#include "DataFormats/DetId/interface/DetId.h"
#include "FWCore/Utilities/interface/Exception.h"
 
#include <sstream>
#include <mutex>
 
using namespace std;
 
namespace tmtt {
 
  //=== Get FE window size arrays (via copy) used with stub producer, but set to zero.
 
  void StubWindowSuggest::setFEWindows(const StubFEWindows* sw) {
    static std::mutex myMutex;
    std::lock_guard<std::mutex> myGuard(myMutex);  // Allow only one thread.
    // Only need to create FE windows once.
    if (not sw_) {
      sw_ = std::make_unique<StubFEWindows>(*sw);  // Copy
      sw_->setZero();
    }
  }
 
  //=== Analyse stub window required for this stub.
 
  void StubWindowSuggest::process(const TrackerTopology* trackerTopo, const Stub* stub) {
    static std::mutex myMutex;
    std::lock_guard<std::mutex> myGuard(myMutex);  // Allow only one thread.
 
    // Half-size of FE chip bend window corresponding to Pt range in which tracks are to be found.
    const double invPtMax = 1 / ptMin_;
    double bendHalfWind = invPtMax / std::abs(stub->qOverPtOverBend());
    // Increase half-indow size to allow for resolution in bend.
    bendHalfWind += stub->bendCutInFrontend();
    // Stub bend is measured here in half-integer values.
    bendHalfWind = int(2 * bendHalfWind) / 2.;
 
    // Compare with half-size of FE bend window stored in arrays.
    this->updateStoredWindow(trackerTopo, stub, bendHalfWind);
  }
 
  //===  Update stored stub window size with this stub.
 
  void StubWindowSuggest::updateStoredWindow(const TrackerTopology* trackerTopo,
                                             const Stub* stub,
                                             double bendHalfWind) {
    // Code accessing geometry inspired by L1Trigger/TrackTrigger/src/TTStubAlgorithm_official.cc
 
    DetId stDetId(stub->trackerModule()->detId());
 
    double* storedHalfWindow = sw_->storedWindowSize(trackerTopo, stDetId);
 
    if (*storedHalfWindow < bendHalfWind)
      *storedHalfWindow = bendHalfWind;
  }
 
  //=== Print results (should be done in endJob();
 
  void StubWindowSuggest::printResults() const {
    PrintL1trk(1) << "==============================================================================";
    PrintL1trk(1) << " Stub window sizes that TMTT suggests putting inside ";
    PrintL1trk(1) << "   /L1Trigger/TrackTrigger/python/TTStubAlgorithmRegister_cfi.py";
    PrintL1trk(1) << " (These should give good efficiency, but tighter windows may be needed to";
    PrintL1trk(1) << "  limit the data rate from the FE tracker electronics).";
    PrintL1trk(1) << "==============================================================================";
 
    std::stringstream text;
    string div;
 
    text << "BarrelCut = cms.vdouble( ";
    div = "";
    for (const auto& cut : sw_->windowSizeBarrelLayers()) {
      text << div << cut;
      div = ", ";
    }
    text << "),";
    PrintL1trk(1) << text.str();
 
    PrintL1trk(1) << "TiltedBarrelCutSet = cms.VPSET( ";
    for (const auto& cutVec : sw_->windowSizeTiltedLayersRings()) {
      text.str("");
      text << "     cms.PSet( TiltedCut = cms.vdouble(";
      if (cutVec.empty())
        text << "0";
      div = "";
      for (const auto& cut : cutVec) {
        text << div << cut;
        div = ", ";
      }
      text << "), ),";
      PrintL1trk(1) << text.str();
    }
    PrintL1trk(1) << "),";
 
    PrintL1trk(1) << "EndcapCutSet = cms.VPSET( ";
    for (const auto& cutVec : sw_->windowSizeEndcapDisksRings()) {
      text.str("");
      text << "     cms.PSet( EndcapCut = cms.vdouble(";
      if (cutVec.empty())
        text << "0";
      div = "";
      for (const auto& cut : cutVec) {
        text << div << cut;
        div = ", ";
      }
      text << "), ),";
      PrintL1trk(1) << text.str();
    }
    PrintL1trk(1) << ")";
 
    PrintL1trk(1) << "==============================================================================";
  }
 
}  // namespace tmtt