CSCTriggerPrimitivesBuilder Class Reference

#include <CSCTriggerPrimitivesBuilder.h>

Public Types
enum	trig_cscs {   MAX_ENDCAPS = 2, MAX_STATIONS = 4, MAX_SECTORS = 6, MAX_SUBSECTORS = 2,   MAX_CHAMBERS = 9 }

Public Member Functions
void	build (const CSCBadChambers *badChambers, const CSCWireDigiCollection *wiredc, const CSCComparatorDigiCollection *compdc, const GEMPadDigiCollection *gemPads, const RPCDigiCollection *rpcDigis, CSCALCTDigiCollection &oc_alct, CSCCLCTDigiCollection &oc_clct, CSCCLCTPreTriggerCollection &oc_pretrig, CSCCorrelatedLCTDigiCollection &oc_lct, CSCCorrelatedLCTDigiCollection &oc_sorted_lct, GEMCoPadDigiCollection &oc_gemcopad, GEMCSCLCTDigiCollection &oc_gemcsclct)
	CSCTriggerPrimitivesBuilder (const edm::ParameterSet &)
void	runOnData (bool runOnData)
void	setConfigParameters (const CSCDBL1TPParameters *conf)
void	setCSCGeometry (const CSCGeometry *g)
	set CSC and GEM geometries for the matching needs More...
void	setGEMGeometry (const GEMGeometry *g)
void	setRPCGeometry (const RPCGeometry *g)
	~CSCTriggerPrimitivesBuilder ()

Private Attributes
bool	checkBadChambers_
	a flag whether to skip chambers from the bad chambers map More...
const CSCGeometry *	csc_g
bool	disableME1a
bool	disableME42
const GEMGeometry *	gem_g
int	m_maxBX
int	m_minBX
std::unique_ptr< CSCMuonPortCard >	m_muonportcard
const RPCGeometry *	rpc_g
bool	runME11ILT_
bool	runME21ILT_
bool	runME3141ILT_
bool	runOnData_
	temporary flag to run on data More...
bool	smartME1aME1b
std::unique_ptr< CSCMotherboard >	tmb_ [MAX_ENDCAPS][MAX_STATIONS][MAX_SECTORS][MAX_SUBSECTORS][MAX_CHAMBERS]

Static Private Attributes
static const int	max_chamber = CSCTriggerNumbering::maxTriggerCscId()
static const int	max_endcap = CSCDetId::maxEndcapId()
static const int	max_sector = CSCTriggerNumbering::maxTriggerSectorId()
static const int	max_station = CSCDetId::maxStationId()
static const int	max_subsector = CSCTriggerNumbering::maxTriggerSubSectorId()
static const int	min_chamber = CSCTriggerNumbering::minTriggerCscId()
static const int	min_endcap = CSCDetId::minEndcapId()
static const int	min_sector = CSCTriggerNumbering::minTriggerSectorId()
static const int	min_station = CSCDetId::minStationId()
static const int	min_subsector = CSCTriggerNumbering::minTriggerSubSectorId()

Detailed Description

Algorithm to build anode, cathode, and correlated LCTs from wire and comparator digis in endcap muon CSCs by implementing a 'build' function required by CSCTriggerPrimitivesProducer.

Configured via the Producer's ParameterSet.

Author

Slava Valuev, UCLA.

Definition at line 37 of file CSCTriggerPrimitivesBuilder.h.

Member Enumeration Documentation

enum CSCTriggerPrimitivesBuilder::trig_cscs

Max values of trigger labels for all CSCs; used to construct TMB processors.

Enumerator
MAX_ENDCAPS
MAX_STATIONS
MAX_SECTORS
MAX_SUBSECTORS
MAX_CHAMBERS

Definition at line 77 of file CSCTriggerPrimitivesBuilder.h.

                 {MAX_ENDCAPS = 2, MAX_STATIONS = 4, MAX_SECTORS = 6,
          MAX_SUBSECTORS = 2, MAX_CHAMBERS = 9};

Constructor & Destructor Documentation

CSCTriggerPrimitivesBuilder::CSCTriggerPrimitivesBuilder ( const edm::ParameterSet &  conf )

explicit

Configure the algorithm via constructor. Receives ParameterSet percolated down from EDProducer which owns this Builder.

Definition at line 49 of file CSCTriggerPrimitivesBuilder.cc.

References relativeConstraints::cham, checkBadChambers_, disableME1a, disableME42, edm::ParameterSet::existsAs(), edm::ParameterSet::getParameter(), m_maxBX, m_minBX, m_muonportcard, max_chamber, MAX_CHAMBERS, max_endcap, MAX_ENDCAPS, max_sector, MAX_SECTORS, max_station, MAX_STATIONS, max_subsector, MAX_SUBSECTORS, min_chamber, min_endcap, min_sector, min_station, min_subsector, relativeConstraints::ring, CSCTriggerNumbering::ringFromTriggerLabels(), runME11ILT_, runME21ILT_, runME3141ILT_, smartME1aME1b, and tmb_.

{
  // Receives ParameterSet percolated down from EDProducer.

  // special configuration parameters for ME11 treatment
  edm::ParameterSet commonParams = conf.getParameter<edm::ParameterSet>("commonParam");
  smartME1aME1b = commonParams.getParameter<bool>("smartME1aME1b");
  disableME1a = commonParams.getParameter<bool>("disableME1a");
  disableME42 = commonParams.getParameter<bool>("disableME42");

  checkBadChambers_ = conf.getParameter<bool>("checkBadChambers");

  runME11ILT_ = commonParams.existsAs<bool>("runME11ILT")?commonParams.getParameter<bool>("runME11ILT"):false;
  runME21ILT_ = commonParams.existsAs<bool>("runME21ILT")?commonParams.getParameter<bool>("runME21ILT"):false;
  runME3141ILT_ = commonParams.existsAs<bool>("runME3141ILT")?commonParams.getParameter<bool>("runME3141ILT"):false;

  // ORCA way of initializing boards.
  for (int endc = min_endcap; endc <= max_endcap; endc++)
  {
    for (int stat = min_station; stat <= max_station; stat++)
    {
      int numsubs = ((stat == 1) ? max_subsector : 1);
      for (int sect = min_sector; sect <= max_sector; sect++)
      {
        for (int subs = min_subsector; subs <= numsubs; subs++)
        {
          for (int cham = min_chamber; cham <= max_chamber; cham++)
          {
            if ((endc <= 0 || endc > MAX_ENDCAPS)    ||
                (stat <= 0 || stat > MAX_STATIONS)   ||
                (sect <= 0 || sect > MAX_SECTORS)    ||
                (subs <= 0 || subs > MAX_SUBSECTORS) ||
                (cham <= 0 || stat > MAX_CHAMBERS))
            {
              edm::LogError("L1CSCTPEmulatorSetupError")
                << "+++ trying to instantiate TMB of illegal CSC id ["
                << " endcap = "  << endc << " station = "   << stat
                << " sector = "  << sect << " subsector = " << subs
                << " chamber = " << cham << "]; skipping it... +++\n";
              continue;
            }
            int ring = CSCTriggerNumbering::ringFromTriggerLabels(stat, cham);
            // When the motherboard is instantiated, it instantiates ALCT
            // and CLCT processors.
            if (stat==1 && ring==1 && smartME1aME1b && !runME11ILT_)
              tmb_[endc-1][stat-1][sect-1][subs-1][cham-1].reset( new CSCMotherboardME11(endc, stat, sect, subs, cham, conf) );
            else if (stat==1 && ring==1 && smartME1aME1b && runME11ILT_)
              tmb_[endc-1][stat-1][sect-1][subs-1][cham-1].reset( new CSCMotherboardME11GEM(endc, stat, sect, subs, cham, conf) );
            else if (stat==2 && ring==1 && runME21ILT_)
          tmb_[endc-1][stat-1][sect-1][subs-1][cham-1].reset( new CSCMotherboardME21GEM(endc, stat, sect, subs, cham, conf) );
            else if ((stat==3 || stat==4) && ring==1 && runME3141ILT_)
              tmb_[endc-1][stat-1][sect-1][subs-1][cham-1].reset( new CSCMotherboardME3141RPC(endc, stat, sect, subs, cham, conf) );
            else
              tmb_[endc-1][stat-1][sect-1][subs-1][cham-1].reset( new CSCMotherboard(endc, stat, sect, subs, cham, conf) );
          }
        }
      }
    }
  }

  // Get min and max BX to sort LCTs in MPC.
  m_minBX = conf.getParameter<int>("MinBX");
  m_maxBX = conf.getParameter<int>("MaxBX");

  // Init MPC
  m_muonportcard.reset( new CSCMuonPortCard(conf) );
}

CSCTriggerPrimitivesBuilder::~CSCTriggerPrimitivesBuilder

(

)

Definition at line 120 of file CSCTriggerPrimitivesBuilder.cc.

{
}

Member Function Documentation

void CSCTriggerPrimitivesBuilder::build	(	const CSCBadChambers *	badChambers,
		const CSCWireDigiCollection *	wiredc,
		const CSCComparatorDigiCollection *	compdc,
		const GEMPadDigiCollection *	gemPads,
		const RPCDigiCollection *	rpcDigis,
		CSCALCTDigiCollection &	oc_alct,
		CSCCLCTDigiCollection &	oc_clct,
		CSCCLCTPreTriggerCollection &	oc_pretrig,
		CSCCorrelatedLCTDigiCollection &	oc_lct,
		CSCCorrelatedLCTDigiCollection &	oc_sorted_lct,
		GEMCoPadDigiCollection &	oc_gemcopad,
		GEMCSCLCTDigiCollection &	oc_gemcsclct
	)

Build anode, cathode, and correlated LCTs in each chamber and fill them into output collections. Select up to three best correlated LCTs in each (sub)sector and put them into an output collection as well.

Definition at line 160 of file CSCTriggerPrimitivesBuilder.cc.

References CSCMotherboard::alct, relativeConstraints::cham, CSCTriggerGeomManager::chamber(), CSCTriggerNumbering::chamberFromTriggerLabels(), checkBadChambers_, CSCMotherboard::clct, CSCMotherboardME11::clct1a, CSCMotherboardME11GEM::clct1a, csc_g, CSCDetId, disableME1a, disableME42, alignCSCRings::e, gem_g, CSCTriggerGeometry::get(), CSCBadChambers::isInBadChamber(), LogDebug, LogTrace, m_maxBX, m_minBX, m_muonportcard, max_chamber, max_endcap, max_sector, max_station, max_subsector, min_chamber, min_endcap, min_sector, min_station, min_subsector, or, CSCMotherboardME11GEM::readoutCoPads(), CSCMotherboardME21GEM::readoutCoPads(), CSCMotherboard::readoutLCTs(), CSCMotherboardME3141RPC::readoutLCTs(), CSCMotherboardME21GEM::readoutLCTs(), CSCMotherboardME11::readoutLCTs1a(), CSCMotherboardME11GEM::readoutLCTs1a(), CSCMotherboardME11::readoutLCTs1b(), CSCMotherboardME11GEM::readoutLCTs1b(), mps_fire::result, relativeConstraints::ring, CSCTriggerNumbering::ringFromTriggerLabels(), rpc_g, CSCMotherboardME11::run(), CSCMotherboardME3141RPC::run(), CSCMotherboardME21GEM::run(), CSCMotherboardME11GEM::run(), CSCMotherboard::run(), runME11ILT_, runME21ILT_, runME3141ILT_, runOnData_, CSCMotherboardME21GEM::setCSCGeometry(), CSCMotherboardME3141RPC::setCSCGeometry(), CSCMotherboardME11GEM::setCSCGeometry(), CSCMotherboardME21GEM::setGEMGeometry(), CSCMotherboardME11GEM::setGEMGeometry(), CSCMotherboardME3141RPC::setRPCGeometry(), smartME1aME1b, tmb_, and CSCDetId::zendcap().

{
  // CSC geometry.
  CSCTriggerGeomManager* theGeom = CSCTriggerGeometry::get();

  for (int endc = min_endcap; endc <= max_endcap; endc++)
  {
    for (int stat = min_station; stat <= max_station; stat++)
    {
      int numsubs = ((stat == 1) ? max_subsector : 1);
      for (int sect = min_sector; sect <= max_sector; sect++)
      {
        for (int subs = min_subsector; subs <= numsubs; subs++)
        {
          for (int cham = min_chamber; cham <= max_chamber; cham++)
          {
            
            int ring = CSCTriggerNumbering::ringFromTriggerLabels(stat, cham);
            
            if (disableME42 && stat==4 && ring==2) continue;

            CSCMotherboard* tmb = tmb_[endc-1][stat-1][sect-1][subs-1][cham-1].get();

            // Run processors only if chamber exists in geometry.
            if (tmb == 0 || theGeom->chamber(endc, stat, sect, subs, cham) == 0) continue;

            int chid = CSCTriggerNumbering::chamberFromTriggerLabels(sect, subs, stat, cham);

            // 0th layer means whole chamber.
            CSCDetId detid(endc, stat, ring, chid, 0);

            // Skip chambers marked as bad (usually includes most of ME4/2 chambers;
            // also, there's no ME1/a-1/b separation, it's whole ME1/1)
            if (checkBadChambers_ && badChambers->isInBadChamber(detid)) continue;


            // running upgraded ME1/1 TMBs
            if (stat==1 && ring==1 && smartME1aME1b && !runME11ILT_)
            {
              CSCMotherboardME11* tmb11 = static_cast<CSCMotherboardME11*>(tmb);
 
              //LogTrace("CSCTriggerPrimitivesBuilder")<<"CSCTriggerPrimitivesBuilder::build in E:"<<endc<<" S:"<<stat<<" R:"<<ring;
 
              tmb11->run(wiredc,compdc);
              std::vector<CSCCorrelatedLCTDigi> lctV = tmb11->readoutLCTs1b();
              std::vector<CSCCorrelatedLCTDigi> lctV1a = tmb11->readoutLCTs1a();
 
              std::vector<CSCALCTDigi> alctV1a, alctV = tmb11->alct->readoutALCTs();

              std::vector<CSCCLCTDigi> clctV = tmb11->clct->readoutCLCTs();
              std::vector<int> preTriggerBXs = tmb11->clct->preTriggerBXs();
              std::vector<CSCCLCTDigi> clctV1a = tmb11->clct1a->readoutCLCTs();
              std::vector<int> preTriggerBXs1a = tmb11->clct1a->preTriggerBXs();

              // perform simple separation of ALCTs into 1/a and 1/b
              // for 'smart' case. Some duplication occurs for WG [10,15]
              std::vector<CSCALCTDigi> tmpV(alctV);
              alctV.clear();
              for (unsigned int al=0; al < tmpV.size(); al++)
              {
                if (tmpV[al].getKeyWG()<=15) alctV1a.push_back(tmpV[al]);
                if (tmpV[al].getKeyWG()>=10) alctV.push_back(tmpV[al]);
              }
              //LogTrace("CSCTriggerPrimitivesBuilder")<<"CSCTriggerPrimitivesBuilder:: a="<<alctV.size()<<" c="<<clctV.size()<<" l="<<lctV.size()
              //  <<"   1a: a="<<alctV1a.size()<<" c="<<clctV1a.size()<<" l="<<lctV1a.size();
 
              // ME1/b

              if (!(lctV.empty()&&alctV.empty()&&clctV.empty())) {
                LogTrace("L1CSCTrigger")
                  << "CSCTriggerPrimitivesBuilder results in " <<detid; 
              }

              // Correlated LCTs.
              if (!lctV.empty()) {
                LogTrace("L1CSCTrigger")
                  << "Put " << lctV.size() << " ME1b LCT digi"
                  << ((lctV.size() > 1) ? "s " : " ") << "in collection\n";
                oc_lct.put(std::make_pair(lctV.begin(),lctV.end()), detid);
              }
 
              // Anode LCTs.
              if (!alctV.empty()) {
                LogTrace("L1CSCTrigger")
                  << "Put " << alctV.size() << " ME1b ALCT digi"
                  << ((alctV.size() > 1) ? "s " : " ") << "in collection\n";
                oc_alct.put(std::make_pair(alctV.begin(),alctV.end()), detid);
              }
 
              // Cathode LCTs.
              if (!clctV.empty()) {
                LogTrace("L1CSCTrigger")
                  << "Put " << clctV.size() << " ME1b CLCT digi"
                  << ((clctV.size() > 1) ? "s " : " ") << "in collection\n";
                oc_clct.put(std::make_pair(clctV.begin(),clctV.end()), detid);
              }

              // Cathode LCTs pretriggers
              if (!preTriggerBXs.empty()) {
                LogTrace("L1CSCTrigger")
                  << "Put " << preTriggerBXs.size() << " CLCT pretrigger"
                  << ((preTriggerBXs.size() > 1) ? "s " : " ") << "in collection\n";
                oc_pretrig.put(std::make_pair(preTriggerBXs.begin(),preTriggerBXs.end()), detid);
              }            

              // ME1/a

              if (disableME1a) continue;

              CSCDetId detid1a(endc, stat, 4, chid, 0);

              if (!(lctV1a.empty()&&alctV1a.empty()&&clctV1a.empty())){
                LogTrace("L1CSCTrigger") << "CSCTriggerPrimitivesBuilder results in " <<detid1a;
              }
 
              // Correlated LCTs.
              if (!lctV1a.empty()) {
                LogTrace("L1CSCTrigger")
                  << "Put " << lctV1a.size() << " ME1a LCT digi"
                  << ((lctV1a.size() > 1) ? "s " : " ") << "in collection\n";
                oc_lct.put(std::make_pair(lctV1a.begin(),lctV1a.end()), detid1a);
              }
 
              // Anode LCTs.
              if (!alctV1a.empty()) {
                LogTrace("L1CSCTrigger")
                  << "Put " << alctV1a.size() << " ME1a ALCT digi"
                  << ((alctV1a.size() > 1) ? "s " : " ") << "in collection\n";
                oc_alct.put(std::make_pair(alctV1a.begin(),alctV1a.end()), detid1a);
              }
 
              // Cathode LCTs.
              if (!clctV1a.empty()) {
                LogTrace("L1CSCTrigger")
                  << "Put " << clctV1a.size() << " ME1a CLCT digi"
                  << ((clctV1a.size() > 1) ? "s " : " ") << "in collection\n";
                oc_clct.put(std::make_pair(clctV1a.begin(),clctV1a.end()), detid1a);
              }
              
              // Cathode LCTs pretriggers
              if (!preTriggerBXs1a.empty()) {
                LogTrace("L1CSCTrigger")
                  << "Put " << preTriggerBXs1a.size() << " CLCT pretrigger"
                  << ((preTriggerBXs1a.size() > 1) ? "s " : " ") << "in collection\n";
                oc_pretrig.put(std::make_pair(preTriggerBXs1a.begin(),preTriggerBXs1a.end()), detid1a);
              }
            } // upgraded TMB

        // running upgraded ME1/1 TMBs with GEMs
            else if (stat==1 && ring==1 && smartME1aME1b && runME11ILT_)
        {
          CSCMotherboardME11GEM* tmb11GEM = static_cast<CSCMotherboardME11GEM*>(tmb);
          
          tmb11GEM->setCSCGeometry(csc_g);
          tmb11GEM->setGEMGeometry(gem_g);
          //LogTrace("CSCTriggerPrimitivesBuilder")<<"CSCTriggerPrimitivesBuilder::build in E:"<<endc<<" S:"<<stat<<" R:"<<ring;
          tmb11GEM->run(wiredc, compdc, gemPads);
              
          std::vector<CSCCorrelatedLCTDigi> lctV = tmb11GEM->readoutLCTs1b();
          std::vector<CSCCorrelatedLCTDigi> lctV1a = tmb11GEM->readoutLCTs1a();
          
          std::vector<CSCALCTDigi> alctV1a, alctV = tmb11GEM->alct->readoutALCTs();
          
          std::vector<CSCCLCTDigi> clctV = tmb11GEM->clct->readoutCLCTs();
          std::vector<int> preTriggerBXs = tmb11GEM->clct->preTriggerBXs();
          std::vector<CSCCLCTDigi> clctV1a = tmb11GEM->clct1a->readoutCLCTs();
          std::vector<int> preTriggerBXs1a = tmb11GEM->clct1a->preTriggerBXs();
          
          std::vector<GEMCoPadDigi> copads = tmb11GEM->readoutCoPads();
          
          // perform simple separation of ALCTs into 1/a and 1/b
          // for 'smart' case. Some duplication occurs for WG [10,15]
          std::vector<CSCALCTDigi> tmpV(alctV);
          alctV.clear();
          for (unsigned int al=0; al < tmpV.size(); al++)
          {
        if (tmpV[al].getKeyWG()<=15) alctV1a.push_back(tmpV[al]);
        if (tmpV[al].getKeyWG()>=10) alctV.push_back(tmpV[al]);
          }
          //LogTrace("CSCTriggerPrimitivesBuilder")<<"CSCTriggerPrimitivesBuilder:: a="<<alctV.size()<<" c="<<clctV.size()<<" l="<<lctV.size()
          //  <<"   1a: a="<<alctV1a.size()<<" c="<<clctV1a.size()<<" l="<<lctV1a.size();
          
          // ME1/b
          
          if (!(lctV.empty()&&alctV.empty()&&clctV.empty())) {
        LogTrace("L1CSCTrigger")
          << "CSCTriggerPrimitivesBuilder results in " <<detid; 
          }
          
          // Correlated LCTs.
          if (!lctV.empty()) {
        LogTrace("L1CSCTrigger")
          << "Put " << lctV.size() << " ME1b LCT digi"
          << ((lctV.size() > 1) ? "s " : " ") << "in collection\n";
        oc_lct.put(std::make_pair(lctV.begin(),lctV.end()), detid);
          }
          
          // Anode LCTs.
          if (!alctV.empty()) {
        LogTrace("L1CSCTrigger")
          << "Put " << alctV.size() << " ME1b ALCT digi"
          << ((alctV.size() > 1) ? "s " : " ") << "in collection\n";
        oc_alct.put(std::make_pair(alctV.begin(),alctV.end()), detid);
          }
          
          // Cathode LCTs.
          if (!clctV.empty()) {
        LogTrace("L1CSCTrigger")
          << "Put " << clctV.size() << " ME1b CLCT digi"
          << ((clctV.size() > 1) ? "s " : " ") << "in collection\n";
        oc_clct.put(std::make_pair(clctV.begin(),clctV.end()), detid);
          }
          
          // Cathode LCTs pretriggers
          if (!preTriggerBXs.empty()) {
        LogTrace("L1CSCTrigger")
          << "Put " << preTriggerBXs.size() << " CLCT pretrigger"
          << ((preTriggerBXs.size() > 1) ? "s " : " ") << "in collection\n";
        oc_pretrig.put(std::make_pair(preTriggerBXs.begin(),preTriggerBXs.end()), detid);
          }            
          // 0th layer means whole chamber.
          GEMDetId gemId(detid.zendcap(), 1, 1, 1, chid, 0);
              
          // GEM coincidence pads
          if (!copads.empty()) {
        LogTrace("L1CSCTrigger")
          << "Put " << copads.size() << " GEM coincidence pad"
          << ((copads.size() > 1) ? "s " : " ") << "in collection\n";
        oc_gemcopad.put(std::make_pair(copads.begin(),copads.end()), gemId);
          }
          
          // ME1/a
          
          if (disableME1a) continue;
          
          CSCDetId detid1a(endc, stat, 4, chid, 0);
          
          if (!(lctV1a.empty()&&alctV1a.empty()&&clctV1a.empty())){
        LogTrace("L1CSCTrigger") << "CSCTriggerPrimitivesBuilder results in " <<detid1a;
          }
          
          // Correlated LCTs.
          if (!lctV1a.empty()) {
        LogTrace("L1CSCTrigger")
          << "Put " << lctV1a.size() << " ME1a LCT digi"
          << ((lctV1a.size() > 1) ? "s " : " ") << "in collection\n";
        oc_lct.put(std::make_pair(lctV1a.begin(),lctV1a.end()), detid1a);
          }
          
          // Anode LCTs.
          if (!alctV1a.empty()) {
        LogTrace("L1CSCTrigger")
          << "Put " << alctV1a.size() << " ME1a ALCT digi"
          << ((alctV1a.size() > 1) ? "s " : " ") << "in collection\n";
        oc_alct.put(std::make_pair(alctV1a.begin(),alctV1a.end()), detid1a);
          }
          
          // Cathode LCTs.
          if (!clctV1a.empty()) {
        LogTrace("L1CSCTrigger")
          << "Put " << clctV1a.size() << " ME1a CLCT digi"
          << ((clctV1a.size() > 1) ? "s " : " ") << "in collection\n";
        oc_clct.put(std::make_pair(clctV1a.begin(),clctV1a.end()), detid1a);
          }
          
          // Cathode LCTs pretriggers
          if (!preTriggerBXs1a.empty()) {
        LogTrace("L1CSCTrigger")
          << "Put " << preTriggerBXs.size() << " CLCT pretrigger"
          << ((preTriggerBXs.size() > 1) ? "s " : " ") << "in collection\n";
        oc_pretrig.put(std::make_pair(preTriggerBXs.begin(),preTriggerBXs.end()), detid);
          }
        } 

            // running upgraded ME2/1 TMBs
            else if (stat==2 && ring==1 && runME21ILT_)
        {
          CSCMotherboardME21GEM* tmb21GEM = static_cast<CSCMotherboardME21GEM*>(tmb);
          tmb21GEM->setCSCGeometry(csc_g);
          tmb21GEM->setGEMGeometry(gem_g);
          tmb21GEM->run(wiredc, compdc, gemPads);
          std::vector<CSCCorrelatedLCTDigi> lctV = tmb21GEM->readoutLCTs();
          std::vector<CSCALCTDigi> alctV = tmb21GEM->alct->readoutALCTs();
          std::vector<CSCCLCTDigi> clctV = tmb21GEM->clct->readoutCLCTs();
          std::vector<int> preTriggerBXs = tmb21GEM->clct->preTriggerBXs();
          
          std::vector<GEMCoPadDigi> copads = tmb21GEM->readoutCoPads();
          
          if (!(alctV.empty() && clctV.empty() && lctV.empty())) {
        LogTrace("L1CSCTrigger")
          << "CSCTriggerPrimitivesBuilder got results in " <<detid;
          }
          
          // Correlated LCTs.
          if (!lctV.empty()) {
        LogTrace("L1CSCTrigger")
          << "Put " << lctV.size() << " LCT digi"
          << ((lctV.size() > 1) ? "s " : " ") << "in collection\n";
        oc_lct.put(std::make_pair(lctV.begin(),lctV.end()), detid);
          }
          
          // Anode LCTs.
          if (!alctV.empty()) {
        LogTrace("L1CSCTrigger")
          << "Put " << alctV.size() << " ALCT digi"
          << ((alctV.size() > 1) ? "s " : " ") << "in collection\n";
        oc_alct.put(std::make_pair(alctV.begin(),alctV.end()), detid);
          }
          
          // Cathode LCTs.
          if (!clctV.empty()) {
        LogTrace("L1CSCTrigger")
          << "Put " << clctV.size() << " CLCT digi"
          << ((clctV.size() > 1) ? "s " : " ") << "in collection\n";
        oc_clct.put(std::make_pair(clctV.begin(),clctV.end()), detid);
          }
          
          // Cathode LCTs pretriggers
          if (!preTriggerBXs.empty()) {
        LogTrace("L1CSCTrigger")
          << "Put " << preTriggerBXs.size() << " CLCT pretrigger"
          << ((preTriggerBXs.size() > 1) ? "s " : " ") << "in collection\n";
        oc_pretrig.put(std::make_pair(preTriggerBXs.begin(),preTriggerBXs.end()), detid);
          }
          
          // 0th layer means whole chamber.
          GEMDetId gemId(detid.zendcap(), 1, 2, 1, chid, 0);
          
          // GEM coincidence pads
          if (!copads.empty()) {
        LogTrace("L1CSCTrigger")
          << "Put " << copads.size() << " GEM coincidence pad"
          << ((copads.size() > 1) ? "s " : " ") << "in collection\n";
        oc_gemcopad.put(std::make_pair(copads.begin(),copads.end()), gemId);
          }
        }
        // running upgraded ME3/1-ME4/1 TMBs
            else if ((stat==3 or stat==4) && ring==1 && runME3141ILT_)
        {
          CSCMotherboardME3141RPC* tmb3141RPC = static_cast<CSCMotherboardME3141RPC*>(tmb);
          tmb3141RPC->setCSCGeometry(csc_g);
          tmb3141RPC->setRPCGeometry(rpc_g);
          tmb3141RPC->run(wiredc, compdc, rpcDigis);
          std::vector<CSCCorrelatedLCTDigi> lctV = tmb3141RPC->readoutLCTs();
          std::vector<CSCALCTDigi> alctV = tmb3141RPC->alct->readoutALCTs();
          std::vector<CSCCLCTDigi> clctV = tmb3141RPC->clct->readoutCLCTs();
          std::vector<int> preTriggerBXs = tmb3141RPC->clct->preTriggerBXs();
          
          if (!(alctV.empty() && clctV.empty() && lctV.empty())) {
        LogTrace("L1CSCTrigger")
          << "CSCTriggerPrimitivesBuilder got results in " <<detid;
          }
          
          // Correlated LCTs.
          if (!lctV.empty()) {
        LogTrace("L1CSCTrigger")
          << "Put " << lctV.size() << " LCT digi"
          << ((lctV.size() > 1) ? "s " : " ") << "in collection\n";
        oc_lct.put(std::make_pair(lctV.begin(),lctV.end()), detid);
          }
          // Anode LCTs.
          if (!alctV.empty()) {
        LogTrace("L1CSCTrigger")
          << "Put " << alctV.size() << " ALCT digi"
          << ((alctV.size() > 1) ? "s " : " ") << "in collection\n";
        oc_alct.put(std::make_pair(alctV.begin(),alctV.end()), detid);
          }
          
          // Cathode LCTs.
          if (!clctV.empty()) {
        LogTrace("L1CSCTrigger")
          << "Put " << clctV.size() << " CLCT digi"
          << ((clctV.size() > 1) ? "s " : " ") << "in collection\n";
        oc_clct.put(std::make_pair(clctV.begin(),clctV.end()), detid);
          }
          
          // Cathode LCTs pretriggers
          if (!preTriggerBXs.empty()) {
        LogTrace("L1CSCTrigger")
          << "Put " << preTriggerBXs.size() << " CLCT pretrigger"
          << ((preTriggerBXs.size() > 1) ? "s " : " ") << "in collection\n";
        oc_pretrig.put(std::make_pair(preTriggerBXs.begin(),preTriggerBXs.end()), detid);
          }
        }       
        
            // running non-upgraded TMB
            else
            {
              tmb->run(wiredc,compdc);

              std::vector<CSCCorrelatedLCTDigi> lctV = tmb->readoutLCTs();
              std::vector<CSCALCTDigi> alctV = tmb->alct->readoutALCTs();
              std::vector<CSCCLCTDigi> clctV = tmb->clct->readoutCLCTs();
              std::vector<int> preTriggerBXs = tmb->clct->preTriggerBXs();

              if (!(alctV.empty() && clctV.empty() && lctV.empty())) {
                LogTrace("L1CSCTrigger")
                  << "CSCTriggerPrimitivesBuilder got results in " <<detid;
              }

              /*
              // tmp kludge: tightening of ME1a LCTs
              if (stat==1 && ring==1) {
                std::vector<CSCCorrelatedLCTDigi> lctV11;
                for (unsigned t=0;t<lctV.size();t++){
                  if (lctV[t].getStrip() < 127) lctV11.push_back(lctV[t]);
                  else if (lctV[t].getQuality() >= 14) lctV11.push_back(lctV[t]);
                }
                lctV = lctV11;
              }
              */

              // Correlated LCTs.
              if (!lctV.empty()) {
                LogTrace("L1CSCTrigger")
                  << "Put " << lctV.size() << " LCT digi"
                  << ((lctV.size() > 1) ? "s " : " ") << "in collection\n";
                oc_lct.put(std::make_pair(lctV.begin(),lctV.end()), detid);
              }

              // Anode LCTs.
              if (!alctV.empty()) {
                LogTrace("L1CSCTrigger")
                  << "Put " << alctV.size() << " ALCT digi"
                  << ((alctV.size() > 1) ? "s " : " ") << "in collection\n";
                oc_alct.put(std::make_pair(alctV.begin(),alctV.end()), detid);
              }

              // Cathode LCTs.
              if (!clctV.empty()) {
                LogTrace("L1CSCTrigger")
                  << "Put " << clctV.size() << " CLCT digi"
                  << ((clctV.size() > 1) ? "s " : " ") << "in collection\n";
                oc_clct.put(std::make_pair(clctV.begin(),clctV.end()), detid);
              }

              // Cathode LCTs pretriggers
              if (!preTriggerBXs.empty()) {
                LogTrace("L1CSCTrigger")
                  << "Put " << preTriggerBXs.size() << " CLCT pretrigger"
                  << ((preTriggerBXs.size() > 1) ? "s " : " ") << "in collection\n";
                oc_pretrig.put(std::make_pair(preTriggerBXs.begin(),preTriggerBXs.end()), detid);
              }
            } // non-upgraded TMB
          }
        }
      }
    }
  }

  // run MPC simulation
  m_muonportcard->loadDigis(oc_lct);

  // temporary hack to ensure that all MPC LCTs are read out
  if (runOnData_) {
    m_minBX = 5;
    m_maxBX = 11;
  }

  std::vector<csctf::TrackStub> result;
  for(int bx = m_minBX; bx <= m_maxBX; ++bx)
    for(int e = min_endcap; e <= max_endcap; ++e)
      for(int st = min_station; st <= max_station; ++st)
        for(int se = min_sector; se <= max_sector; ++se)
        {
          if(st == 1)
          {
            std::vector<csctf::TrackStub> subs1, subs2;
            subs1 = m_muonportcard->sort(e, st, se, 1, bx);
            subs2 = m_muonportcard->sort(e, st, se, 2, bx);
            result.insert(result.end(), subs1.begin(), subs1.end());
            result.insert(result.end(), subs2.begin(), subs2.end());
          }
          else
          {
            std::vector<csctf::TrackStub> sector;
            sector = m_muonportcard->sort(e, st, se, 0, bx);
            result.insert(result.end(), sector.begin(), sector.end());
          }
        }

  std::vector<csctf::TrackStub>::const_iterator itr = result.begin();
  for (; itr != result.end(); itr++)
  {
    oc_sorted_lct.insertDigi(CSCDetId(itr->getDetId().rawId()), *(itr->getDigi()));
    LogDebug("L1CSCTrigger")
      << "MPC " << *(itr->getDigi()) << " found in ME"
      << ((itr->endcap() == 1) ? "+" : "-") << itr->station() << "/"
      << CSCDetId(itr->getDetId().rawId()).ring() << "/"
      << CSCDetId(itr->getDetId().rawId()).chamber()
      << " (sector " << itr->sector()
      << " trig id. " << itr->cscid() << ")" << "\n";
  }
}

void CSCTriggerPrimitivesBuilder::runOnData ( bool  runOnData )

inline

Definition at line 58 of file CSCTriggerPrimitivesBuilder.h.

References runOnData(), and runOnData_.

Referenced by runOnData().

{runOnData_ = runOnData;}

void CSCTriggerPrimitivesBuilder::setConfigParameters

(

const CSCDBL1TPParameters *

conf

)

Sets configuration parameters obtained via EventSetup mechanism.

Definition at line 128 of file CSCTriggerPrimitivesBuilder.cc.

References relativeConstraints::cham, max_chamber, max_endcap, max_sector, max_station, max_subsector, min_chamber, min_endcap, min_sector, min_station, min_subsector, and tmb_.

{
  // Receives CSCDBL1TPParameters percolated down from ESProducer.

  for (int endc = min_endcap; endc <= max_endcap; endc++)
  {
    for (int stat = min_station; stat <= max_station; stat++)
    {
      int numsubs = ((stat == 1) ? max_subsector : 1);
      for (int sect = min_sector; sect <= max_sector; sect++)
      {
        for (int subs = min_subsector; subs <= numsubs; subs++)
        {
          for (int cham = min_chamber; cham <= max_chamber; cham++)
          {
            tmb_[endc-1][stat-1][sect-1][subs-1][cham-1]->setConfigParameters(conf);
          }
        }
      }
    }
  }
}

void CSCTriggerPrimitivesBuilder::setCSCGeometry ( const CSCGeometry *  g )

inline

set CSC and GEM geometries for the matching needs

Definition at line 53 of file CSCTriggerPrimitivesBuilder.h.

References csc_g, and g.

{ csc_g = g; }

void CSCTriggerPrimitivesBuilder::setGEMGeometry ( const GEMGeometry *  g )

inline

Definition at line 54 of file CSCTriggerPrimitivesBuilder.h.

References g, and gem_g.

{ gem_g = g; }

void CSCTriggerPrimitivesBuilder::setRPCGeometry ( const RPCGeometry *  g )

inline

Definition at line 55 of file CSCTriggerPrimitivesBuilder.h.

References g, and rpc_g.

{ rpc_g = g; }

Member Data Documentation

bool CSCTriggerPrimitivesBuilder::checkBadChambers_

private

a flag whether to skip chambers from the bad chambers map

Definition at line 98 of file CSCTriggerPrimitivesBuilder.h.

Referenced by build(), and CSCTriggerPrimitivesBuilder().

const CSCGeometry* CSCTriggerPrimitivesBuilder::csc_g

private

Definition at line 121 of file CSCTriggerPrimitivesBuilder.h.

Referenced by build(), and setCSCGeometry().

bool CSCTriggerPrimitivesBuilder::disableME1a

private

Definition at line 101 of file CSCTriggerPrimitivesBuilder.h.

Referenced by build(), and CSCTriggerPrimitivesBuilder().

bool CSCTriggerPrimitivesBuilder::disableME42

private

SLHC: special switch for disabling ME42

Definition at line 104 of file CSCTriggerPrimitivesBuilder.h.

Referenced by build(), and CSCTriggerPrimitivesBuilder().

const GEMGeometry* CSCTriggerPrimitivesBuilder::gem_g

private

Definition at line 122 of file CSCTriggerPrimitivesBuilder.h.

Referenced by build(), and setGEMGeometry().

int CSCTriggerPrimitivesBuilder::m_maxBX

private

Definition at line 115 of file CSCTriggerPrimitivesBuilder.h.

Referenced by build(), and CSCTriggerPrimitivesBuilder().

int CSCTriggerPrimitivesBuilder::m_minBX

private

Definition at line 115 of file CSCTriggerPrimitivesBuilder.h.

Referenced by build(), and CSCTriggerPrimitivesBuilder().

std::unique_ptr<CSCMuonPortCard> CSCTriggerPrimitivesBuilder::m_muonportcard

private

Pointer to MPC processor.

Definition at line 126 of file CSCTriggerPrimitivesBuilder.h.

Referenced by build(), and CSCTriggerPrimitivesBuilder().

const int CSCTriggerPrimitivesBuilder::max_chamber = CSCTriggerNumbering::maxTriggerCscId()

staticprivate

Definition at line 92 of file CSCTriggerPrimitivesBuilder.h.

Referenced by build(), CSCTriggerPrimitivesBuilder(), and setConfigParameters().

const int CSCTriggerPrimitivesBuilder::max_endcap = CSCDetId::maxEndcapId()

staticprivate

Definition at line 84 of file CSCTriggerPrimitivesBuilder.h.

Referenced by build(), CSCTriggerPrimitivesBuilder(), and setConfigParameters().

const int CSCTriggerPrimitivesBuilder::max_sector = CSCTriggerNumbering::maxTriggerSectorId()

staticprivate

Definition at line 88 of file CSCTriggerPrimitivesBuilder.h.

Referenced by build(), CSCTriggerPrimitivesBuilder(), and setConfigParameters().

const int CSCTriggerPrimitivesBuilder::max_station = CSCDetId::maxStationId()

staticprivate

Definition at line 86 of file CSCTriggerPrimitivesBuilder.h.

Referenced by build(), CSCTriggerPrimitivesBuilder(), and setConfigParameters().

const int CSCTriggerPrimitivesBuilder::max_subsector = CSCTriggerNumbering::maxTriggerSubSectorId()

staticprivate

Definition at line 90 of file CSCTriggerPrimitivesBuilder.h.

Referenced by build(), CSCTriggerPrimitivesBuilder(), and setConfigParameters().

const int CSCTriggerPrimitivesBuilder::min_chamber = CSCTriggerNumbering::minTriggerCscId()

staticprivate

Definition at line 91 of file CSCTriggerPrimitivesBuilder.h.

Referenced by build(), CSCTriggerPrimitivesBuilder(), and setConfigParameters().

const int CSCTriggerPrimitivesBuilder::min_endcap = CSCDetId::minEndcapId()

staticprivate

Min and max allowed values for various CSC elements, defined in CSCDetId and CSCTriggerNumbering classes.

Definition at line 83 of file CSCTriggerPrimitivesBuilder.h.

Referenced by build(), CSCTriggerPrimitivesBuilder(), and setConfigParameters().

const int CSCTriggerPrimitivesBuilder::min_sector = CSCTriggerNumbering::minTriggerSectorId()

staticprivate

Definition at line 87 of file CSCTriggerPrimitivesBuilder.h.

Referenced by build(), CSCTriggerPrimitivesBuilder(), and setConfigParameters().

const int CSCTriggerPrimitivesBuilder::min_station = CSCDetId::minStationId()

staticprivate

Definition at line 85 of file CSCTriggerPrimitivesBuilder.h.

Referenced by build(), CSCTriggerPrimitivesBuilder(), and setConfigParameters().

const int CSCTriggerPrimitivesBuilder::min_subsector = CSCTriggerNumbering::minTriggerSubSectorId()

staticprivate

Definition at line 89 of file CSCTriggerPrimitivesBuilder.h.

Referenced by build(), CSCTriggerPrimitivesBuilder(), and setConfigParameters().

const RPCGeometry* CSCTriggerPrimitivesBuilder::rpc_g

private

Definition at line 123 of file CSCTriggerPrimitivesBuilder.h.

Referenced by build(), and setRPCGeometry().

bool CSCTriggerPrimitivesBuilder::runME11ILT_

private

SLHC: special switch for the upgrade ME1/1 TMB

Definition at line 107 of file CSCTriggerPrimitivesBuilder.h.

Referenced by build(), and CSCTriggerPrimitivesBuilder().

bool CSCTriggerPrimitivesBuilder::runME21ILT_

private

SLHC: special switch for the upgrade ME2/1 TMB

Definition at line 110 of file CSCTriggerPrimitivesBuilder.h.

Referenced by build(), and CSCTriggerPrimitivesBuilder().

bool CSCTriggerPrimitivesBuilder::runME3141ILT_

private

SLHC: special switch for the upgrade ME3/1 and ME4/1 TMB

Definition at line 113 of file CSCTriggerPrimitivesBuilder.h.

Referenced by build(), and CSCTriggerPrimitivesBuilder().

bool CSCTriggerPrimitivesBuilder::runOnData_

private

temporary flag to run on data

Definition at line 95 of file CSCTriggerPrimitivesBuilder.h.

Referenced by build(), and runOnData().

bool CSCTriggerPrimitivesBuilder::smartME1aME1b

private

SLHC: special configuration parameters for ME11 treatment.

Definition at line 101 of file CSCTriggerPrimitivesBuilder.h.

Referenced by build(), and CSCTriggerPrimitivesBuilder().

std::unique_ptr<CSCMotherboard> CSCTriggerPrimitivesBuilder::tmb_[MAX_ENDCAPS][MAX_STATIONS][MAX_SECTORS][MAX_SUBSECTORS][MAX_CHAMBERS]

private

Pointers to TMB processors for all possible chambers.

Definition at line 119 of file CSCTriggerPrimitivesBuilder.h.

Referenced by build(), CSCTriggerPrimitivesBuilder(), and setConfigParameters().