LCTQualityControl Class Reference

#include <LCTQualityControl.h>

Inheritance diagram for LCTQualityControl:

Public Member Functions
void	checkMultiplicityBX (const std::vector< CSCALCTDigi > &alcts) const
void	checkMultiplicityBX (const std::vector< CSCCLCTDigi > &clcts) const
void	checkMultiplicityBX (const std::vector< CSCCorrelatedLCTDigi > &lcts) const
template<class T >
void	checkMultiplicityBX (const std::vector< T > &lcts, unsigned nLCT) const
void	checkRange (int parameter, int min_value, int max_value, const std::string &comment, unsigned &errors) const
void	checkValid (const CSCALCTDigi &alct, unsigned max_stubs=CSCConstants::MAX_ALCTS_PER_PROCESSOR) const
void	checkValid (const CSCCLCTDigi &lct, unsigned max_stubs=CSCConstants::MAX_CLCTS_PER_PROCESSOR) const
void	checkValid (const CSCCorrelatedLCTDigi &lct) const
void	checkValid (const CSCCorrelatedLCTDigi &lct, const unsigned station, const unsigned ring) const
void	checkValidReadout (const CSCALCTDigi &alct) const
std::pair< unsigned, unsigned >	get_csc_alct_min_max_quality () const
std::pair< unsigned, unsigned >	get_csc_clct_min_max_quality () const
std::pair< int, int >	get_csc_clct_min_max_slope () const
std::pair< unsigned, unsigned >	get_csc_lct_min_max_pattern () const
std::pair< unsigned, unsigned >	get_csc_lct_min_max_quality (unsigned station, unsigned ring) const
unsigned	get_csc_max_eighthstrip (unsigned station, unsigned ring) const
unsigned	get_csc_max_halfstrip (unsigned station, unsigned ring) const
unsigned	get_csc_max_quartstrip (unsigned station, unsigned ring) const
unsigned	get_csc_max_wiregroup (unsigned station, unsigned ring) const
std::pair< unsigned, unsigned >	get_csc_min_max_cfeb () const
std::pair< unsigned, unsigned >	get_csc_min_max_pattern () const
std::pair< unsigned, unsigned >	get_csc_min_max_pattern_run3 () const
int	getSlopePhase1 (unsigned pattern) const
	LCTQualityControl (unsigned endcap, unsigned station, unsigned sector, unsigned subsector, unsigned chamber, const edm::ParameterSet &conf)
template<class T >
void	reportErrors (const T &lct, const unsigned errors) const
	~LCTQualityControl () override=default
Public Member Functions inherited from CSCBaseboard
	CSCBaseboard (unsigned endcap, unsigned station, unsigned sector, unsigned subsector, unsigned chamber, const edm::ParameterSet &conf)
	CSCBaseboard ()
std::string	getCSCName () const
CSCDetId	id () const
void	setCSCGeometry (const CSCGeometry *g)
virtual	~CSCBaseboard ()=default

Private Attributes
unsigned	nplanes_clct_hit_pattern

Additional Inherited Members
Protected Member Functions inherited from CSCBaseboard
void	checkConfigParameters (unsigned int &var, const unsigned int var_max, const unsigned int var_def, const std::string &var_str)
Protected Attributes inherited from CSCBaseboard
edm::ParameterSet	alctParams_
edm::ParameterSet	clctParams_
edm::ParameterSet	commonParams_
const CSCChamber *	cscChamber_
const CSCGeometry *	cscGeometry_
CSCDetId	cscId_
bool	disableME1a_
bool	enableAlctPhase2_
bool	gangedME1a_
int	infoV
bool	isME11_
bool	isME12_
bool	isME13_
bool	isME21_
bool	isME22_
bool	isME31_
bool	isME32_
bool	isME41_
bool	isME42_
bool	runCCLUT_
bool	runCCLUT_OTMB_
bool	runCCLUT_TMB_
bool	runME11ILT_
bool	runME11Up_
bool	runME21ILT_
bool	runME21Up_
bool	runME31Up_
bool	runME41Up_
bool	runPhase2_
edm::ParameterSet	showerParams_
unsigned	theChamber
std::string	theCSCName_
const unsigned	theEndcap
int	theRegion
unsigned	theRing
const unsigned	theSector
const unsigned	theStation
const unsigned	theSubsector
const unsigned	theTrigChamber
edm::ParameterSet	tmbParams_

Detailed Description

Definition at line 22 of file LCTQualityControl.h.

Constructor & Destructor Documentation

LCTQualityControl::LCTQualityControl	(	unsigned	endcap,
		unsigned	station,
		unsigned	sector,
		unsigned	subsector,
		unsigned	chamber,
		const edm::ParameterSet &	conf
	)

Definition at line 6 of file LCTQualityControl.cc.

References CSCBaseboard::clctParams_, edm::ParameterSet::getParameter(), and nplanes_clct_hit_pattern.

    : CSCBaseboard(endcap, station, sector, subsector, chamber, conf) {
  nplanes_clct_hit_pattern = clctParams_.getParameter<unsigned int>("clctNplanesHitPattern");
}

LCTQualityControl::~LCTQualityControl ( )

overridedefault

Default destructor.

Member Function Documentation

void LCTQualityControl::checkMultiplicityBX

(

const std::vector< CSCALCTDigi > &

alcts

)

const

Definition at line 284 of file LCTQualityControl.cc.

References CSCConstants::MAX_ALCTS_READOUT.

Referenced by checkMultiplicityBX().

                                                                                          {
  checkMultiplicityBX(collection, CSCConstants::MAX_ALCTS_READOUT);
}

void LCTQualityControl::checkMultiplicityBX

(

const std::vector< CSCCLCTDigi > &

clcts

)

const

Definition at line 288 of file LCTQualityControl.cc.

References checkMultiplicityBX(), and CSCConstants::MAX_CLCTS_READOUT.

                                                                                          {
  checkMultiplicityBX(collection, CSCConstants::MAX_CLCTS_READOUT);
}

void LCTQualityControl::checkMultiplicityBX

(

const std::vector< CSCCorrelatedLCTDigi > &

lcts

)

const

Definition at line 292 of file LCTQualityControl.cc.

References checkMultiplicityBX(), and CSCConstants::MAX_LCTS_PER_CSC.

                                                                                                   {
  checkMultiplicityBX(collection, CSCConstants::MAX_LCTS_PER_CSC);
}

template<class T >

void LCTQualityControl::checkMultiplicityBX	(	const std::vector< T > &	lcts,
		unsigned	nLCT
	)		const

Definition at line 100 of file LCTQualityControl.h.

References AlCaHLTBitMon_ParallelJobs::p.

                                                                                               {
  std::unordered_map<int, unsigned> freq;
  // check each BX
  for (const auto& p : collection) {
    if (p.isValid()) {
      freq[p.getBX()]++;

      // too many ALCTs, CLCTs or LCTs in this BX
      if (freq[p.getBX()] > nLCT) {
        edm::LogError("LCTQualityControl") << "Collection with more than " << nLCT << " in BX " << p.getBX();
        break;
      }
    }
  }
}

void LCTQualityControl::checkRange	(	int	parameter,
		int	min_value,
		int	max_value,
		const std::string &	comment,
		unsigned &	errors
	)		const

Definition at line 21 of file LCTQualityControl.cc.

Referenced by checkValid().

                                                                                               {
  if (value < min_value or value > max_value) {
    edm::LogError("LCTQualityControl") << comment << value << "; allowed [" << min_value << ", " << max_value << "]";
    errors++;
  }
}

void LCTQualityControl::checkValid	(	const CSCALCTDigi &	alct,
		unsigned	max_stubs = CSCConstants::MAX_ALCTS_PER_PROCESSOR
	)		const

Definition at line 37 of file LCTQualityControl.cc.

References checkRange(), get_csc_alct_min_max_quality(), get_csc_max_wiregroup(), CSCALCTDigi::getBX(), CSCALCTDigi::getCollisionB(), CSCALCTDigi::getKeyWG(), CSCALCTDigi::getQuality(), CSCALCTDigi::getTrknmb(), CSCALCTDigi::isValid(), CSCConstants::MAX_ALCT_TBINS, reportErrors(), CSCBaseboard::theRing, and CSCBaseboard::theStation.

Referenced by checkValid(), and checkValidReadout().

                                                                                    {
  const unsigned max_wiregroup = get_csc_max_wiregroup(theStation, theRing);
  const auto& [min_quality, max_quality] = get_csc_alct_min_max_quality();

  unsigned errors = 0;

  // stub must be valid
  checkRange(alct.isValid(), 1, 1, "CSCALCTDigi with invalid bit set: ", errors);

  // ALCT number is 1 or 2
  checkRange(alct.getTrknmb(), 1, max_stubs, "CSCALCTDigi with invalid track number: ", errors);

  // ALCT quality must be valid
  // number of layers - 3
  checkRange(alct.getQuality(), min_quality, max_quality, "CSCALCTDigi with invalid quality: ", errors);

  // ALCT key wire-group must be within bounds
  checkRange(alct.getKeyWG(), 0, max_wiregroup - 1, "CSCALCTDigi with invalid wire-group: ", errors);

  // ALCT with out-of-time BX
  checkRange(alct.getBX(), 0, CSCConstants::MAX_ALCT_TBINS - 1, "CSCALCTDigi with invalid BX: ", errors);

  // ALCT is neither accelerator or collision
  checkRange(alct.getCollisionB(), 0, 1, "CSCALCTDigi with invalid accel/coll biit: ", errors);

  reportErrors(alct, errors);
}

void LCTQualityControl::checkValid	(	const CSCCLCTDigi &	lct,
		unsigned	max_stubs = CSCConstants::MAX_CLCTS_PER_PROCESSOR
	)		const

Definition at line 66 of file LCTQualityControl.cc.

References checkRange(), get_csc_clct_min_max_quality(), get_csc_clct_min_max_slope(), get_csc_max_eighthstrip(), get_csc_max_halfstrip(), get_csc_max_quartstrip(), get_csc_min_max_cfeb(), get_csc_min_max_pattern(), get_csc_min_max_pattern_run3(), CSCCLCTDigi::getBend(), CSCCLCTDigi::getBX(), CSCCLCTDigi::getCFEB(), CSCCLCTDigi::getCompCode(), CSCCLCTDigi::getKeyStrip(), CSCCLCTDigi::getPattern(), CSCCLCTDigi::getQuality(), CSCCLCTDigi::getRun3Pattern(), CSCCLCTDigi::getSlope(), CSCCLCTDigi::getStrip(), CSCCLCTDigi::getTrknmb(), CSCCLCTDigi::isValid(), CSCConstants::MAX_CLCT_TBINS, CSCConstants::NUM_HALF_STRIPS_PER_CFEB, funct::pow(), reportErrors(), CSCBaseboard::runCCLUT_, CSCBaseboard::theRing, and CSCBaseboard::theStation.

                                                                                    {
  const unsigned max_strip = get_csc_max_halfstrip(theStation, theRing);
  const auto& [min_pattern_run2, max_pattern_run2] = get_csc_min_max_pattern();
  const auto& [min_pattern_run3, max_pattern_run3] = get_csc_min_max_pattern_run3();
  const auto& [min_slope, max_slope] = get_csc_clct_min_max_slope();
  const auto& [min_cfeb, max_cfeb] = get_csc_min_max_cfeb();
  const auto& [min_quality, max_quality] = get_csc_clct_min_max_quality();

  unsigned errors = 0;

  // CLCT must be valid
  checkRange(clct.isValid(), 1, 1, "CSCCLCTDigi with invalid bit set: ", errors);

  // CLCT number is 1 or max
  checkRange(clct.getTrknmb(), 1, max_stubs, "CSCCLCTDigi with invalid track number: ", errors);

  // CLCT quality must be valid
  // CLCTs require at least 4 layers hit
  // Run-3: ME1/1 CLCTs require only 3 layers
  // Run-4: ME2/1 CLCTs require only 3 layers
  checkRange(clct.getQuality(), min_quality, max_quality, "CSCCLCTDigi with invalid quality: ", errors);

  // CLCT half-strip must be within bounds
  checkRange(
      clct.getStrip(), 0, CSCConstants::NUM_HALF_STRIPS_PER_CFEB - 1, "CSCCLCTDigi with invalid half-strip: ", errors);

  // CLCT key half-strip must be within bounds
  checkRange(clct.getKeyStrip(), 0, max_strip - 1, "CSCCLCTDigi with invalid key half-strip: ", errors);

  // CLCT with out-of-time BX
  checkRange(clct.getBX(), 0, CSCConstants::MAX_CLCT_TBINS - 1, "CSCCLCTDigi with invalid BX: ", errors);

  // CLCT with neither left nor right bending
  checkRange(clct.getBend(), 0, 1, "CSCCLCTDigi with invalid bending: ", errors);

  // CLCT with an invalid pattern ID
  checkRange(
      clct.getPattern(), min_pattern_run2, max_pattern_run2, "CSCCLCTDigi with invalid Run-2 pattern ID: ", errors);

  // CLCT with an invalid pattern ID
  checkRange(
      clct.getRun3Pattern(), min_pattern_run3, max_pattern_run3, "CSCCLCTDigi with invalid Run-3 pattern ID: ", errors);

  // CLCT with an invalid slope
  checkRange(clct.getSlope(), min_slope, max_slope, "CSCCLCTDigi with invalid slope: ", errors);

  // CLCT with an invalid CFEB ID
  checkRange(clct.getCFEB(), min_cfeb, max_cfeb, "CSCCLCTDigi with invalid CFEB ID: ", errors);

  if (runCCLUT_) {
    // CLCT comparator code is invalid
    checkRange(clct.getCompCode(), 0, std::pow(2, 12) - 1, "CSCCLCTDigi with invalid comparator code: ", errors);

    const unsigned max_quartstrip = get_csc_max_quartstrip(theStation, theRing);
    const unsigned max_eighthstrip = get_csc_max_eighthstrip(theStation, theRing);

    // CLCT key quart-strip must be within bounds
    checkRange(clct.getKeyStrip(4), 0, max_quartstrip - 1, "CSCCLCTDigi with invalid key quart-strip: ", errors);

    // CLCT key eighth-strip must be within bounds
    checkRange(clct.getKeyStrip(8), 0, max_eighthstrip - 1, "CSCCLCTDigi with invalid key quart-strip: ", errors);
  }

  reportErrors(clct, errors);
}

void LCTQualityControl::checkValid

(

const CSCCorrelatedLCTDigi &

lct

)

const

Definition at line 132 of file LCTQualityControl.cc.

References checkValid(), CSCBaseboard::theRing, and CSCBaseboard::theStation.

{ checkValid(lct, theStation, theRing); }

void LCTQualityControl::checkValid	(	const CSCCorrelatedLCTDigi &	lct,
		const unsigned	station,
		const unsigned	ring
	)		const

Definition at line 134 of file LCTQualityControl.cc.

References CSCCorrelatedLCTDigi::ALCT2GEM, LCTQualityAssignment::ALCT_2GEM, LCTQualityAssignment::ALCT_CLCT, LCTQualityAssignment::ALCT_CLCT_1GEM_CSCBend, LCTQualityAssignment::ALCT_CLCT_1GEM_GEMCSCBend, LCTQualityAssignment::ALCT_CLCT_2GEM_CSCBend, LCTQualityAssignment::ALCT_CLCT_2GEM_GEMCSCBend, CSCCorrelatedLCTDigi::ALCTCLCT, CSCCorrelatedLCTDigi::ALCTCLCT2GEM, CSCCorrelatedLCTDigi::ALCTCLCTGEM, CSCCorrelatedLCTDigi::ALCTONLY, checkRange(), CSCCorrelatedLCTDigi::CLCT2GEM, LCTQualityAssignment::CLCT_2GEM, CSCCorrelatedLCTDigi::CLCTALCT, CSCCorrelatedLCTDigi::CLCTONLY, get_csc_lct_min_max_quality(), get_csc_max_eighthstrip(), get_csc_max_halfstrip(), get_csc_max_quartstrip(), get_csc_max_wiregroup(), get_csc_min_max_pattern(), get_csc_min_max_pattern_run3(), CSCCorrelatedLCTDigi::getBend(), CSCCorrelatedLCTDigi::getBX(), CSCCorrelatedLCTDigi::getCSCID(), CSCCorrelatedLCTDigi::getGEM1(), CSCCorrelatedLCTDigi::getGEM2(), CSCCorrelatedLCTDigi::getKeyWG(), CSCCorrelatedLCTDigi::getMPCLink(), CSCCorrelatedLCTDigi::getPattern(), CSCCorrelatedLCTDigi::getQuality(), CSCCorrelatedLCTDigi::getRun3Pattern(), CSCCorrelatedLCTDigi::getStrip(), CSCCorrelatedLCTDigi::getTrknmb(), CSCCorrelatedLCTDigi::getType(), LCTQualityAssignment::HighQ, CSCBaseboard::isME11_, CSCBaseboard::isME21_, GEMPadDigi::isValid(), CSCCorrelatedLCTDigi::isValid(), LCTQualityAssignment::LowQ, CSCConstants::MAX_LCT_TBINS, CSCConstants::MAX_LCTS_PER_MPC, CSCTriggerNumbering::maxTriggerCscId(), LCTQualityAssignment::MedQ, CSCTriggerNumbering::minTriggerCscId(), or, reportErrors(), CSCBaseboard::runCCLUT_, CSCBaseboard::runME11ILT_, and CSCBaseboard::runME21ILT_.

                                                                                                         {
  const unsigned max_strip = get_csc_max_halfstrip(station, ring);
  const unsigned max_quartstrip = get_csc_max_quartstrip(station, ring);
  const unsigned max_eighthstrip = get_csc_max_eighthstrip(station, ring);
  const unsigned max_wiregroup = get_csc_max_wiregroup(station, ring);
  const auto& [min_pattern_run2, max_pattern_run2] = get_csc_min_max_pattern();
  const auto& [min_pattern_run3, max_pattern_run3] = get_csc_min_max_pattern_run3();
  const auto& [min_quality, max_quality] = get_csc_lct_min_max_quality(station, ring);

  unsigned errors = 0;

  // LCT must be valid
  checkRange(lct.isValid(), 1, 1, "CSCCorrelatedLCTDigi with invalid bit set: ", errors);

  // LCT number is 1 or 2
  checkRange(lct.getTrknmb(), 1, 2, "CSCCorrelatedLCTDigi with invalid track number: ", errors);

  // LCT quality must be valid
  checkRange(lct.getQuality(), min_quality, max_quality, "CSCCorrelatedLCTDigi with invalid quality: ", errors);

  // LCT key half-strip must be within bounds
  checkRange(lct.getStrip(), 0, max_strip - 1, "CSCCorrelatedLCTDigi with invalid key half-strip: ", errors);

  // LCT key quart-strip must be within bounds
  checkRange(lct.getStrip(4), 0, max_quartstrip - 1, "CSCCorrelatedLCTDigi with invalid key quart-strip: ", errors);

  // LCT key eighth-strip must be within bounds
  checkRange(lct.getStrip(8), 0, max_eighthstrip - 1, "CSCCorrelatedLCTDigi with invalid key eighth-strip: ", errors);

  // LCT key wire-group must be within bounds
  checkRange(lct.getKeyWG(), 0, max_wiregroup - 1, "CSCCorrelatedLCTDigi with invalid wire-group: ", errors);

  // LCT with out-of-time BX
  checkRange(lct.getBX(), 0, CSCConstants::MAX_LCT_TBINS - 1, "CSCCorrelatedLCTDigi with invalid BX: ", errors);

  // LCT with neither left nor right bending
  checkRange(lct.getBend(), 0, 1, "CSCCorrelatedLCTDigi with invalid bending: ", errors);

  // LCT with invalid MPC link
  checkRange(lct.getMPCLink(), 0, CSCConstants::MAX_LCTS_PER_MPC, "CSCCorrelatedLCTDigi with MPC link: ", errors);

  // LCT with invalid CSCID
  checkRange(lct.getCSCID(),
             CSCTriggerNumbering::minTriggerCscId(),
             CSCTriggerNumbering::maxTriggerCscId(),
             "CSCCorrelatedLCTDigi with invalid CSCID: ",
             errors);

  // LCT with an invalid pattern ID
  checkRange(lct.getPattern(),
             min_pattern_run2,
             max_pattern_run2,
             "CSCCorrelatedLCTDigi with invalid Run-2 pattern ID: ",
             errors);

  checkRange(lct.getRun3Pattern(),
             min_pattern_run3,
             max_pattern_run3,
             "CSCCorrelatedLCTDigi with invalid Run-3 pattern ID: ",
             errors);

  // simulated LCT type must be valid
  if (lct.getType() == CSCCorrelatedLCTDigi::CLCTALCT or lct.getType() == CSCCorrelatedLCTDigi::CLCTONLY or
      lct.getType() == CSCCorrelatedLCTDigi::ALCTONLY) {
    edm::LogError("LCTQualityControl") << "CSCCorrelatedLCTDigi with invalid type (SIM): " << lct.getType()
                                       << "; allowed [" << CSCCorrelatedLCTDigi::ALCTCLCT << ", "
                                       << CSCCorrelatedLCTDigi::CLCT2GEM << "]";
    errors++;
  }

  // non-GEM-CSC stations ALWAYS send out ALCTCLCT type LCTs
  if (!isME11_ and !isME21_) {
    if (lct.getType() != CSCCorrelatedLCTDigi::ALCTCLCT) {
      edm::LogError("LCTQualityControl") << "CSCCorrelatedLCTDigi with invalid type (SIM) in this station: "
                                         << lct.getType() << "; allowed [" << CSCCorrelatedLCTDigi::ALCTCLCT << "]";
      errors++;
    }
  }

  // GEM-CSC stations can send out GEM-type LCTs ONLY when the ILT is turned on!
  if ((isME11_ and !runME11ILT_) or (isME21_ and !runME21ILT_)) {
    if (lct.getType() != CSCCorrelatedLCTDigi::ALCTCLCT) {
      edm::LogError("LCTQualityControl") << "CSCCorrelatedLCTDigi with invalid type (SIM) with GEM-CSC trigger not on: "
                                         << lct.getType() << "; allowed [" << CSCCorrelatedLCTDigi::ALCTCLCT << "]";
      errors++;
    }
  }

  // GEM-CSC types must have at least one valid GEM hit
  if ((lct.getType() == CSCCorrelatedLCTDigi::ALCTCLCTGEM or lct.getType() == CSCCorrelatedLCTDigi::ALCTCLCT2GEM or
       lct.getType() == CSCCorrelatedLCTDigi::ALCT2GEM or lct.getType() == CSCCorrelatedLCTDigi::CLCT2GEM) and
      !lct.getGEM1().isValid() and !lct.getGEM2().isValid()) {
    edm::LogError("LCTQualityControl") << "CSCCorrelatedLCTDigi with valid GEM-CSC type (SIM) has no valid GEM hits: "
                                       << lct.getType();
    errors++;
  }

  // LCT type does not agree with the LCT quality when CCLUT is on
  if (runCCLUT_) {
    const bool ME11ILT(isME11_ and runME11ILT_);
    const bool ME21ILT(isME21_ and runME21ILT_);

    // GEM-CSC cases
    if (ME11ILT or ME21ILT) {
      const bool case1(lct.getType() == CSCCorrelatedLCTDigi::ALCT2GEM and
                       lct.getQuality() == static_cast<unsigned>(LCTQualityAssignment::LCT_QualityRun3GEM::ALCT_2GEM));
      const bool case2(lct.getType() == CSCCorrelatedLCTDigi::CLCT2GEM and
                       lct.getQuality() == static_cast<unsigned>(LCTQualityAssignment::LCT_QualityRun3GEM::CLCT_2GEM));
      const bool case3(lct.getType() == CSCCorrelatedLCTDigi::ALCTCLCTGEM and
                       lct.getQuality() ==
                           static_cast<unsigned>(LCTQualityAssignment::LCT_QualityRun3GEM::ALCT_CLCT_1GEM_CSCBend));
      const bool case4(lct.getType() == CSCCorrelatedLCTDigi::ALCTCLCTGEM and
                       lct.getQuality() ==
                           static_cast<unsigned>(LCTQualityAssignment::LCT_QualityRun3GEM::ALCT_CLCT_1GEM_GEMCSCBend));
      const bool case5(lct.getType() == CSCCorrelatedLCTDigi::ALCTCLCT2GEM and
                       lct.getQuality() ==
                           static_cast<unsigned>(LCTQualityAssignment::LCT_QualityRun3GEM::ALCT_CLCT_2GEM_CSCBend));
      const bool case6(lct.getType() == CSCCorrelatedLCTDigi::ALCTCLCT2GEM and
                       lct.getQuality() ==
                           static_cast<unsigned>(LCTQualityAssignment::LCT_QualityRun3GEM::ALCT_CLCT_2GEM_GEMCSCBend));
      const bool case7(lct.getType() == CSCCorrelatedLCTDigi::ALCTCLCT and
                       lct.getQuality() == static_cast<unsigned>(LCTQualityAssignment::LCT_QualityRun3GEM::ALCT_CLCT));

      if (!(case1 or case2 or case3 or case4 or case5 or case6 or case7)) {
        edm::LogError("LCTQualityControl")
            << "CSCCorrelatedLCTDigi with valid GEM-CSC type (SIM) has no matching Run-3 quality: " << lct.getType()
            << " " << lct.getQuality();
        errors++;
      }
    }

    // regular cases
    else {
      const bool case1(lct.getType() == CSCCorrelatedLCTDigi::ALCTCLCT and
                       lct.getQuality() == static_cast<unsigned>(LCTQualityAssignment::LCT_QualityRun3::LowQ));
      const bool case2(lct.getType() == CSCCorrelatedLCTDigi::ALCTCLCT and
                       lct.getQuality() == static_cast<unsigned>(LCTQualityAssignment::LCT_QualityRun3::MedQ));
      const bool case3(lct.getType() == CSCCorrelatedLCTDigi::ALCTCLCT and
                       lct.getQuality() == static_cast<unsigned>(LCTQualityAssignment::LCT_QualityRun3::HighQ));
      if (!(case1 or case2 or case3)) {
        edm::LogError("LCTQualityControl")
            << "CSCCorrelatedLCTDigi with invalid CSC type (SIM) has no matching Run-3 quality: " << lct.getType()
            << " " << lct.getQuality();
        errors++;
      }
    }
  }
  reportErrors(lct, errors);
}

void LCTQualityControl::checkValidReadout

(

const CSCALCTDigi &

alct

)

const

Definition at line 17 of file LCTQualityControl.cc.

References checkValid(), and CSCConstants::MAX_ALCTS_READOUT.

                                                                       {
  checkValid(alct, CSCConstants::MAX_ALCTS_READOUT);
}

std::pair< unsigned, unsigned > LCTQualityControl::get_csc_alct_min_max_quality

(

)

const

Definition at line 463 of file LCTQualityControl.cc.

References CSCBaseboard::isME21_, and CSCBaseboard::runME21ILT_.

Referenced by checkValid().

                                                                                  {
  // quality is number of layers - 3
  // at least 4 layers for CSC-only trigger
  unsigned min_quality = 1;
  if (isME21_ and runME21ILT_) {
    // at least 3 layers for GEM-CSC trigger in ME2/1
    min_quality = 0;
  }
  return std::make_pair(min_quality, 3);
}

std::pair< unsigned, unsigned > LCTQualityControl::get_csc_clct_min_max_quality

(

)

const

Definition at line 474 of file LCTQualityControl.cc.

References CSCBaseboard::isME11_, CSCBaseboard::isME21_, or, CSCBaseboard::runME11ILT_, and CSCBaseboard::runME21ILT_.

Referenced by checkValid().

                                                                                  {
  // quality is number of layers
  unsigned min_quality = 4;
  if ((runME11ILT_ and isME11_) or (runME21ILT_ and isME21_)) {
    // at least 3 layers for GEM-CSC trigger in ME1/1 or ME2/1
    min_quality = 3;
  }
  return std::make_pair(min_quality, 6);
}

std::pair< int, int > LCTQualityControl::get_csc_clct_min_max_slope

(

)

const

Definition at line 303 of file LCTQualityControl.cc.

References CSCBaseboard::runCCLUT_.

Referenced by checkValid().

                                                                      {
  int min_slope, max_slope;
  // Run-3 case with CCLUT
  // 5-bit number (includes the L/R bending)
  if (runCCLUT_) {
    min_slope = 0;
    max_slope = 15;
  }
  // Run-1 or Run-2 case
  // Run-3 case without CCLUT
  else {
    min_slope = -10;
    max_slope = 10;
  }

  return std::make_pair(min_slope, max_slope);
}

std::pair< unsigned, unsigned > LCTQualityControl::get_csc_lct_min_max_pattern

(

)

const

Definition at line 448 of file LCTQualityControl.cc.

References CSCBaseboard::runCCLUT_.

                                                                                 {
  unsigned min_pattern, max_pattern;
  // Run-1 or Run-2 case
  if (!runCCLUT_) {
    min_pattern = 2;
    max_pattern = 10;
  }
  // Run-3 case: pattern Id is interpreted as the slope!
  else {
    min_pattern = 0;
    max_pattern = 15;
  }
  return std::make_pair(min_pattern, max_pattern);
}

std::pair< unsigned, unsigned > LCTQualityControl::get_csc_lct_min_max_quality	(	unsigned	station,
		unsigned	ring
	)		const

Definition at line 484 of file LCTQualityControl.cc.

References LCTQualityAssignment::ALCT_2GEM, LCTQualityAssignment::ALCT_CLCT_2GEM_GEMCSCBend, LCTQualityAssignment::HighQ, LCTQualityAssignment::HQ_ANODE_MARGINAL_CATHODE, LCTQualityAssignment::HQ_PATTERN_10, LCTQualityAssignment::HQ_PATTERN_2_3, CSCBaseboard::isME11_, CSCBaseboard::isME21_, LCTQualityAssignment::LowQ, or, CSCBaseboard::runCCLUT_, CSCBaseboard::runME11ILT_, and CSCBaseboard::runME21ILT_.

Referenced by checkValid().

                                                                                                                {
  // Run-1 or Run-2
  unsigned min_quality = static_cast<unsigned>(LCTQualityAssignment::LCT_QualityRun2::HQ_PATTERN_2_3);
  unsigned max_quality = static_cast<unsigned>(LCTQualityAssignment::LCT_QualityRun2::HQ_PATTERN_10);

  const bool GEMCSC = (isME11_ and runME11ILT_) or (isME21_ and runME21ILT_);

  // Run-3 with CCLUT on
  if (runCCLUT_ and !GEMCSC) {
    min_quality = static_cast<unsigned>(LCTQualityAssignment::LCT_QualityRun3::LowQ);
    max_quality = static_cast<unsigned>(LCTQualityAssignment::LCT_QualityRun3::HighQ);
  }

  // Run-3 with GEM-CSC on (low-quality CLCTs are permitted, but use Run-2 data format)
  if (!runCCLUT_ and GEMCSC) {
    min_quality = static_cast<unsigned>(LCTQualityAssignment::LCT_QualityRun2::HQ_ANODE_MARGINAL_CATHODE);
    max_quality = static_cast<unsigned>(LCTQualityAssignment::LCT_QualityRun2::HQ_PATTERN_10);
  }

  // Run-3 CSC with GEM-CSC on and CCLUT on
  if (runCCLUT_ and GEMCSC) {
    min_quality = static_cast<unsigned>(LCTQualityAssignment::LCT_QualityRun3GEM::ALCT_2GEM);
    max_quality = static_cast<unsigned>(LCTQualityAssignment::LCT_QualityRun3GEM::ALCT_CLCT_2GEM_GEMCSCBend);
  }
  return std::make_pair(min_quality, max_quality);
}

unsigned LCTQualityControl::get_csc_max_eighthstrip	(	unsigned	station,
		unsigned	ring
	)		const

Definition at line 398 of file LCTQualityControl.cc.

References get_csc_max_halfstrip().

Referenced by checkValid().

                                                                                         {
  return get_csc_max_halfstrip(station, ring) * 4;
}

unsigned LCTQualityControl::get_csc_max_halfstrip	(	unsigned	station,
		unsigned	ring
	)		const

Definition at line 355 of file LCTQualityControl.cc.

References CSCBaseboard::gangedME1a_, CSCConstants::NUM_HALF_STRIPS_ME11_GANGED, CSCConstants::NUM_HALF_STRIPS_ME11_UNGANGED, CSCConstants::NUM_HALF_STRIPS_ME12, CSCConstants::NUM_HALF_STRIPS_ME13, CSCConstants::NUM_HALF_STRIPS_ME1A_GANGED, CSCConstants::NUM_HALF_STRIPS_ME1A_UNGANGED, CSCConstants::NUM_HALF_STRIPS_ME21, CSCConstants::NUM_HALF_STRIPS_ME22, and CSCConstants::NUM_HALF_STRIPS_ME31.

Referenced by checkValid(), get_csc_max_eighthstrip(), and get_csc_max_quartstrip().

                                                                                       {
  unsigned max_half_strip = 0;
  // ME1/1a
  if (station == 1 && ring == 4 and gangedME1a_) {
    max_half_strip = CSCConstants::NUM_HALF_STRIPS_ME1A_GANGED;
  } else if (station == 1 && ring == 4 and !gangedME1a_) {
    max_half_strip = CSCConstants::NUM_HALF_STRIPS_ME1A_UNGANGED;
  }
  // ME1/1b
  // In the CSC local trigger
  // ME1/a is taken together with ME1/b
  else if (station == 1 && ring == 1 and gangedME1a_) {
    max_half_strip = CSCConstants::NUM_HALF_STRIPS_ME11_GANGED;
  } else if (station == 1 && ring == 1 and !gangedME1a_) {
    max_half_strip = CSCConstants::NUM_HALF_STRIPS_ME11_UNGANGED;
  }
  // ME1/2
  else if (station == 1 && ring == 2) {
    max_half_strip = CSCConstants::NUM_HALF_STRIPS_ME12;
  }
  // ME1/3
  else if (station == 1 && ring == 3) {
    max_half_strip = CSCConstants::NUM_HALF_STRIPS_ME13;
  }
  // ME2/1
  else if (station == 2 && ring == 1) {
    max_half_strip = CSCConstants::NUM_HALF_STRIPS_ME21;
  }
  // ME3/1, ME4/1
  else if (station >= 3 && ring == 1) {
    max_half_strip = CSCConstants::NUM_HALF_STRIPS_ME31;
  }
  // ME2/2, ME3/2, ME4/2
  else if (station >= 2 && ring == 2) {
    max_half_strip = CSCConstants::NUM_HALF_STRIPS_ME22;
  }
  return max_half_strip;
}

unsigned LCTQualityControl::get_csc_max_quartstrip	(	unsigned	station,
		unsigned	ring
	)		const

Definition at line 394 of file LCTQualityControl.cc.

References get_csc_max_halfstrip().

Referenced by checkValid().

                                                                                        {
  return get_csc_max_halfstrip(station, ring) * 2;
}

unsigned LCTQualityControl::get_csc_max_wiregroup	(	unsigned	station,
		unsigned	ring
	)		const

Definition at line 335 of file LCTQualityControl.cc.

References CSCConstants::NUM_WIREGROUPS_ME11, CSCConstants::NUM_WIREGROUPS_ME12, CSCConstants::NUM_WIREGROUPS_ME13, CSCConstants::NUM_WIREGROUPS_ME21, CSCConstants::NUM_WIREGROUPS_ME22, and CSCConstants::NUM_WIREGROUPS_ME31.

Referenced by checkValid().

                                                                                       {
  unsigned max_wiregroup = 0;
  if (station == 1 && ring == 4) {  // ME1/1a
    max_wiregroup = CSCConstants::NUM_WIREGROUPS_ME11;
  } else if (station == 1 && ring == 1) {  // ME1/1b
    max_wiregroup = CSCConstants::NUM_WIREGROUPS_ME11;
  } else if (station == 1 && ring == 2) {  // ME1/2
    max_wiregroup = CSCConstants::NUM_WIREGROUPS_ME12;
  } else if (station == 1 && ring == 3) {  // ME1/3
    max_wiregroup = CSCConstants::NUM_WIREGROUPS_ME13;
  } else if (station == 2 && ring == 1) {  // ME2/1
    max_wiregroup = CSCConstants::NUM_WIREGROUPS_ME21;
  } else if (station >= 3 && ring == 1) {  // ME3/1, ME4/1
    max_wiregroup = CSCConstants::NUM_WIREGROUPS_ME31;
  } else if (station >= 2 && ring == 2) {  // ME2/2, ME3/2, ME4/2
    max_wiregroup = CSCConstants::NUM_WIREGROUPS_ME22;
  }
  return max_wiregroup;
}

std::pair< unsigned, unsigned > LCTQualityControl::get_csc_min_max_cfeb

(

)

const

Definition at line 402 of file LCTQualityControl.cc.

References CSCBaseboard::gangedME1a_, CSCConstants::NUM_CFEBS_ME11_GANGED, CSCConstants::NUM_CFEBS_ME11_UNGANGED, CSCConstants::NUM_CFEBS_ME12, CSCConstants::NUM_CFEBS_ME13, CSCConstants::NUM_CFEBS_ME1A_GANGED, CSCConstants::NUM_CFEBS_ME1A_UNGANGED, CSCConstants::NUM_CFEBS_ME21, CSCConstants::NUM_CFEBS_ME22, CSCConstants::NUM_CFEBS_ME31, CSCBaseboard::theRing, and CSCBaseboard::theStation.

Referenced by checkValid().

                                                                          {
  // counts at 0!
  unsigned min_cfeb = 0;
  unsigned max_cfeb = 0;

  // ME1/1a
  if (theStation == 1 && theRing == 4 and gangedME1a_) {
    max_cfeb = CSCConstants::NUM_CFEBS_ME1A_GANGED;
  } else if (theStation == 1 && theRing == 4 and !gangedME1a_) {
    max_cfeb = CSCConstants::NUM_CFEBS_ME1A_UNGANGED;
  }
  // ME1/1b
  // In the CSC local trigger
  // ME1/a is taken together with ME1/b
  else if (theStation == 1 && theRing == 1 and gangedME1a_) {
    max_cfeb = CSCConstants::NUM_CFEBS_ME11_GANGED;
  } else if (theStation == 1 && theRing == 1 and !gangedME1a_) {
    max_cfeb = CSCConstants::NUM_CFEBS_ME11_UNGANGED;
  }
  // ME1/2
  else if (theStation == 1 && theRing == 2) {
    max_cfeb = CSCConstants::NUM_CFEBS_ME12;
  }
  // ME1/3
  else if (theStation == 1 && theRing == 3) {
    max_cfeb = CSCConstants::NUM_CFEBS_ME13;
  }
  // ME2/1
  else if (theStation == 2 && theRing == 1) {
    max_cfeb = CSCConstants::NUM_CFEBS_ME21;
  }
  // ME3/1, ME4/1
  else if (theStation >= 3 && theRing == 1) {
    max_cfeb = CSCConstants::NUM_CFEBS_ME31;
  }
  // ME2/2, ME3/2, ME4/2
  else if (theStation >= 2 && theRing == 2) {
    max_cfeb = CSCConstants::NUM_CFEBS_ME22;
  }
  return std::make_pair(min_cfeb, max_cfeb - 1);
}

std::pair< unsigned, unsigned > LCTQualityControl::get_csc_min_max_pattern

(

)

const

Definition at line 444 of file LCTQualityControl.cc.

Referenced by checkValid().

{ return std::make_pair(2, 10); }

std::pair< unsigned, unsigned > LCTQualityControl::get_csc_min_max_pattern_run3

(

)

const

Definition at line 446 of file LCTQualityControl.cc.

Referenced by checkValid().

{ return std::make_pair(0, 4); }

int LCTQualityControl::getSlopePhase1

(

unsigned

pattern

)

const

Definition at line 296 of file LCTQualityControl.cc.

References CSCConstants::NUM_CLCT_PATTERNS.

                                                            {
  // PID 2 is actually a left-bending pattern with a negative slope
  // PID 3 is actually a right-bending pattern with a positive slope
  int slopeList[CSCConstants::NUM_CLCT_PATTERNS] = {0, 0, -8, 8, -6, 6, -4, 4, -2, 2, 0};
  return slopeList[pattern];
}

template<class T >

void LCTQualityControl::reportErrors	(	const T &	lct,
		const unsigned	errors
	)		const

Definition at line 30 of file LCTQualityControl.cc.

References CSCBaseboard::cscId_.

Referenced by checkValid().

                                                                              {
  if (errors > 0) {
    edm::LogError("LCTQualityControl") << "Invalid stub in " << cscId_ << " (" << errors << " errors):\n" << lct;
  }
}

Member Data Documentation

unsigned LCTQualityControl::nplanes_clct_hit_pattern

private

Definition at line 96 of file LCTQualityControl.h.

Referenced by LCTQualityControl().