ME0TriggerPseudoBuilder Class Reference

#include <ME0TriggerPseudoBuilder.h>

Public Types
enum	trig_me0s { trig_me0s::MAX_ENDCAPS = 2, trig_me0s::MAX_CHAMBERS = 18 }

Public Member Functions
void	build (const ME0SegmentCollection *me0segments, ME0TriggerDigiCollection &oc_trig)
void	dumpAllME0Segments (const ME0SegmentCollection &segments) const
	ME0TriggerPseudoBuilder (const edm::ParameterSet &)
void	setME0Geometry (const ME0Geometry *g)
	~ME0TriggerPseudoBuilder ()

Private Member Functions
ME0TriggerDigi	segmentConversion (const ME0Segment segment)

Private Attributes
edm::ParameterSet	config_
double	dphiresolution_
int	info_
const ME0Geometry *	me0_g

Static Private Attributes
static const int	max_chamber
static const int	max_endcap
static const unsigned int	ME0KeyLayer = 3
static const int	ME0TriggerCentralBX = 8
static const int	min_chamber
static const int	min_endcap

Detailed Description

Builds ME0 trigger objects from ME0 segment

Author

Tao Huang (TAMU)

Definition at line 20 of file ME0TriggerPseudoBuilder.h.

Member Enumeration Documentation

trig_me0s

enum ME0TriggerPseudoBuilder::trig_me0s

strong

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

Enumerator
MAX_ENDCAPS
MAX_CHAMBERS

Definition at line 42 of file ME0TriggerPseudoBuilder.h.

{ MAX_ENDCAPS = 2, MAX_CHAMBERS = 18 };

Constructor & Destructor Documentation

ME0TriggerPseudoBuilder()

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

explicit

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

Definition at line 12 of file ME0TriggerPseudoBuilder.cc.

References config_, dphiresolution_, edm::ParameterSet::getUntrackedParameter(), and info_.

                                                                            {
  config_ = conf;
  info_ = config_.getUntrackedParameter<int>("info", 0);
  dphiresolution_ = config_.getUntrackedParameter<double>("DeltaPhiResolution", 0.25);
}

~ME0TriggerPseudoBuilder()

ME0TriggerPseudoBuilder::~ME0TriggerPseudoBuilder

(

)

Definition at line 18 of file ME0TriggerPseudoBuilder.cc.

{}

Member Function Documentation

build()

void ME0TriggerPseudoBuilder::build	(	const ME0SegmentCollection *	me0segments,
		ME0TriggerDigiCollection &	oc_trig
	)

Build Triggers from ME0 segment in each chamber and fill them into output collections.

Definition at line 20 of file ME0TriggerPseudoBuilder.cc.

References relativeConstraints::cham, dumpAllME0Segments(), info_, LogTrace, MAX_CHAMBERS, MAX_ENDCAPS, me0Segments_cfi::me0Segments, nano_mu_digi_cff::region, segmentConversion(), and ClusterTask_cfi::trig.

                                                                                                              {
  if (info_ > 2)
    dumpAllME0Segments(*me0Segments);

  for (int endc = 0; endc < static_cast<int>(trig_me0s::MAX_ENDCAPS); endc++) {
    for (int cham = 0; cham < static_cast<int>(trig_me0s::MAX_CHAMBERS); cham++) {
      // 0th layer means whole chamber.
      // chamber counts from 1 to 18 in ME0ID
      const int region(endc == 0 ? -1 : 1);
      ME0DetId detid(region, 0, cham + 1, 0);

      const auto& drange = me0Segments->get(detid);
      std::vector<ME0TriggerDigi> trigV;
      for (auto digiIt = drange.first; digiIt != drange.second; digiIt++) {
        if (info_ > 1)
          LogTrace("L1ME0Trigger") << "ME0TriggerPseudoBuilder id " << detid << " ME0 segment " << *digiIt
                                   << " to be converted into trigger digi\n";
        ME0TriggerDigi trig = segmentConversion(*digiIt);
        if (trig.isValid())
          trigV.push_back(trig);
        if (info_ > 1 and trig.isValid())
          LogTrace("L1ME0Trigger") << " ME0trigger " << trig << "\n";
        else if (info_ > 1)
          LogTrace("L1ME0Trigger") << " ME0trigger is not valid. Conversion failed \n";
      }

      if (!trigV.empty()) {
        LogTrace("L1ME0Trigger") << "ME0TriggerPseudoBuilder got results in " << detid << std::endl
                                 << "Put " << trigV.size() << " Trigger digi" << ((trigV.size() > 1) ? "s " : " ")
                                 << "in collection\n";
        oc_trig.put(std::make_pair(trigV.begin(), trigV.end()), detid);
      }
    }
  }
}

dumpAllME0Segments()

void ME0TriggerPseudoBuilder::dumpAllME0Segments

(

const ME0SegmentCollection &

segments

)

const

Definition at line 143 of file ME0TriggerPseudoBuilder.cc.

References runTauDisplay::gp, ME0Geometry::idToDet(), LogTrace, me0_g, ME0RecHit::me0Id(), FastTrackerRecHitMaskProducer_cfi::recHits, and GeomDet::toGlobal().

Referenced by build().

                                                                                           {
  LogTrace("L1ME0Trigger") << "dumpt all ME0 Segments" << std::endl;
  for (auto iC = segments.id_begin(); iC != segments.id_end(); ++iC) {
    auto ch_segs = segments.get(*iC);
    for (auto iS = ch_segs.first; iS != ch_segs.second; ++iS) {
      GlobalPoint gp = me0_g->idToDet(iS->me0DetId())->surface().toGlobal(iS->localPosition());
      LogTrace("L1ME0Trigger") << "ME0Detid " << iS->me0DetId() << " segment " << *iS << " eta " << gp.eta() << " phi "
                               << gp.phi() << std::endl;
      auto recHits(iS->recHits());
      LogTrace("L1ME0Trigger") << "\t has " << recHits.size() << " me0 rechits" << std::endl;
      for (auto& rh : recHits) {
        const ME0RecHit* me0rh(dynamic_cast<const ME0RecHit*>(rh));
        LogTrace("L1ME0Trigger") << "\t  detid " << me0rh->me0Id() << " rechit " << *me0rh << std::endl;
      }
    }
  }
}

segmentConversion()

ME0TriggerDigi ME0TriggerPseudoBuilder::segmentConversion ( const ME0Segment  segment )

private

Definition at line 56 of file ME0TriggerPseudoBuilder.cc.

References cms::cuda::assert(), L1TStage2uGTEmulatorClient_cff::BX, ME0EtaPartition::centreOfStrip(), ME0Geometry::chamber(), relativeConstraints::chamber, dphiresolution_, PV3DBase< T, PVType, FrameType >::eta(), ME0Layer::etaPartition(), spr::find(), runTauDisplay::gp, ME0Geometry::idToDet(), info_, ME0Chamber::layer(), ME0Segment::localDirection(), ME0EtaPartition::localPitch(), ME0Segment::localPosition(), LogTrace, me0_g, ME0Segment::me0DetId(), ME0KeyLayer, ME0TriggerCentralBX, ME0Layer::nEtaPartitions(), ME0Segment::nRecHits(), ME0EtaPartition::nstrips(), or, l1ctLayer1_patternWriters_cff::partition, PV3DBase< T, PVType, FrameType >::phi(), quality, mps_fire::result, ME0Segment::specificRecHits(), nano_mu_digi_cff::strip, ME0EtaPartition::strip(), GeomDet::surface(), hcalRecHitTable_cff::time, ME0Segment::time(), GeomDet::toGlobal(), and GloballyPositioned< T >::toLocal().

Referenced by build().

                                                                                  {
  auto detid = segment.me0DetId();
  const ME0Chamber* chamber = me0_g->chamber(detid);
  const ME0Layer* keylayer = me0_g->chamber(detid)->layer(ME0TriggerPseudoBuilder::ME0KeyLayer);
  int chamberid = detid.chamber() % 2;
  int totRolls = keylayer->nEtaPartitions();
  float dphi = chamber->computeDeltaPhi(segment.localPosition(), segment.localDirection());
  float time = segment.time();
  int sign_time = (time > 0) ? 1 : -1;
  int nrechits = segment.nRecHits();
  std::vector<int> rolls;
  for (const auto& rechit : segment.specificRecHits()) {
    if (std::find(rolls.begin(), rolls.end(), rechit.me0Id().roll()) == rolls.end())
      rolls.push_back(rechit.me0Id().roll());
  }
  if (rolls.size() > 2 or rolls.empty())
    LogTrace("L1ME0Trigger") << " ME0 segment is crossing " << rolls.size() << " roll !!! \n";
  //assert(rolls.size() <= 2);   // we did found very few ME0 segments crossing 3 rolls!!! this cut is applied offline
  if (rolls.empty())
    return ME0TriggerDigi();
  if (rolls[0] < 1)
    LogTrace("L1ME0Trigger") << " ME0 segment has wrong roll number " << rolls[0] << " which should be >= 1 \n !!!";
  assert(rolls[0] >= 1);
  int partition = (rolls[0] - 1) << 1;  //roll from detid counts from 1
  if (rolls.size() == 2 and rolls[0] > rolls[1])
    partition = partition - 1;
  else if (rolls.size() == 2 and rolls[0] < rolls[1])
    partition = partition + 1;

  if (partition < 0 or partition >= 2 * totRolls) {
    LogTrace("L1ME0Trigger") << " ME0 segment rolls size of all hits " << rolls.size() << " rolls[0] " << rolls[0]
                             << " rolls.back() " << rolls.back() << " and ME0 trigger roll is " << partition
                             << " max expected " << 2 * totRolls - 1 << "\n";
    return ME0TriggerDigi();
  }

  //globalpoint from ME0 segment
  GlobalPoint gp = me0_g->idToDet(segment.me0DetId())->surface().toGlobal(segment.localPosition());
  const ME0EtaPartition* etapart = keylayer->etaPartition(rolls[0]);
  LocalPoint segment_lp = etapart->surface().toLocal(gp);  // convert segment gp into lp in etapartition coordinate
  float strippitch = etapart->localPitch(segment_lp);
  float strip = etapart->strip(segment_lp);
  int totstrip = etapart->nstrips();
  int istrip = static_cast<int>(strip);
  int phiposition = istrip;
  if (phiposition > totstrip)
    LogTrace("L1ME0Trigger") << " ME0 segment strip number is " << phiposition << " larger than nstrip " << totstrip
                             << " !!! \n";
  float phi_resolution = 0.5;                                                         //halfstrip
  int phiposition2 = (static_cast<int>((strip - phiposition) / phi_resolution) & 1);  // half-strip resolution
  phiposition = (phiposition << 1) | phiposition2;

  //gloablpoint from ME0 trigger digi
  float centreOfStrip = istrip + 0.25 + phiposition2 * 0.5;
  GlobalPoint gp_digi = etapart->toGlobal(etapart->centreOfStrip(centreOfStrip));

  float strippitch_rad = strippitch / gp.perp();  //unit in rad

  int idphi = static_cast<int>(fabs(dphi) / (strippitch_rad * dphiresolution_));
  const int max_idphi = 512;
  if (idphi >= max_idphi) {
    LogTrace("L1ME0Trigger") << " ME0 segment dphi " << dphi << " and int type: " << idphi
                             << " larger than max allowed: " << max_idphi << " !!! \n";
    idphi = max_idphi - 1;
  }
  int quality = nrechits;  // attention: not the same as discussed in meeting
  int BX = (static_cast<int>(fabs(time) / 25.0)) * sign_time + ME0TriggerPseudoBuilder::ME0TriggerCentralBX;
  int bend = (dphi > 0.0) ? 0 : 1;
  if (info_ > 2)
    LogTrace("L1ME0Trigger") << " ME0trigger in conversion function:\n"
                             << "\t chamber(1-18) " << detid.chamber() << " chamber id " << chamberid << " \n"
                             << "\t rolls size of all hits " << rolls.size() << " rolls[0] " << rolls[0]
                             << " rolls.back() " << rolls.back() << " roll " << partition << " \n"
                             << "\t nRechits " << nrechits << " quality " << quality << " \n"
                             << "\t strip(float) " << strip << " (int) " << istrip << " phiposition " << phiposition
                             << " resolution (in term of strip) " << phi_resolution << " \n"
                             << "\t deltaphi(float) " << dphi << " (int) " << idphi << " resolution "
                             << strippitch_rad * dphiresolution_ << " bend " << bend << " \n"
                             << "\t global point eta " << gp.eta() << " phi " << gp.phi() << " trigger digi eta "
                             << gp_digi.eta() << " phi " << gp_digi.phi() << " \n"
                             << "\t time (ns, float) " << time << " BX " << BX << " \n";

  ME0TriggerDigi result = ME0TriggerDigi(chamberid, quality, phiposition, partition, idphi, bend, BX);
  result.setStrip(istrip);
  return result;
}

setME0Geometry()

void ME0TriggerPseudoBuilder::setME0Geometry ( const ME0Geometry *  g )

inline

set geometry for the matching needs

Definition at line 34 of file ME0TriggerPseudoBuilder.h.

References g, and me0_g.

{ me0_g = g; }

Member Data Documentation

config_

edm::ParameterSet ME0TriggerPseudoBuilder::config_

private

Definition at line 60 of file ME0TriggerPseudoBuilder.h.

Referenced by ME0TriggerPseudoBuilder().

dphiresolution_

double ME0TriggerPseudoBuilder::dphiresolution_

private

Definition at line 56 of file ME0TriggerPseudoBuilder.h.

Referenced by ME0TriggerPseudoBuilder(), and segmentConversion().

info_

int ME0TriggerPseudoBuilder::info_

private

Definition at line 54 of file ME0TriggerPseudoBuilder.h.

Referenced by build(), ME0TriggerPseudoBuilder(), and segmentConversion().

max_chamber

const int ME0TriggerPseudoBuilder::max_chamber

staticprivate

Definition at line 48 of file ME0TriggerPseudoBuilder.h.

max_endcap

const int ME0TriggerPseudoBuilder::max_endcap

staticprivate

Definition at line 46 of file ME0TriggerPseudoBuilder.h.

me0_g

const ME0Geometry* ME0TriggerPseudoBuilder::me0_g

private

Definition at line 52 of file ME0TriggerPseudoBuilder.h.

Referenced by dumpAllME0Segments(), segmentConversion(), and setME0Geometry().

ME0KeyLayer

const unsigned int ME0TriggerPseudoBuilder::ME0KeyLayer = 3

staticprivate

Definition at line 49 of file ME0TriggerPseudoBuilder.h.

Referenced by segmentConversion().

ME0TriggerCentralBX

const int ME0TriggerPseudoBuilder::ME0TriggerCentralBX = 8

staticprivate

Definition at line 50 of file ME0TriggerPseudoBuilder.h.

Referenced by segmentConversion().

min_chamber

const int ME0TriggerPseudoBuilder::min_chamber

staticprivate

Definition at line 47 of file ME0TriggerPseudoBuilder.h.

min_endcap

const int ME0TriggerPseudoBuilder::min_endcap

staticprivate

Definition at line 45 of file ME0TriggerPseudoBuilder.h.