PrimitiveConverterRPC Class Reference

#include <PrimitiveConverterRPC.h>

Public Member Functions
l1t::EMTFHit2016ExtraCollection	convert (std::vector< L1TMuon::TriggerPrimitive > TriggPrim, int SectIndex, edm::ESHandle< RPCGeometry > rpc_geom)
std::vector< ConvertedHit >	fillConvHits (l1t::EMTFHit2016ExtraCollection exHits)
	PrimitiveConverterRPC ()
bool	sameRpcChamber (l1t::EMTFHit2016Extra hitA, l1t::EMTFHit2016Extra hitB)

Detailed Description

Definition at line 20 of file PrimitiveConverterRPC.h.

Constructor & Destructor Documentation

PrimitiveConverterRPC::PrimitiveConverterRPC

(

)

Definition at line 9 of file PrimitiveConverterRPC.cc.

                                             {
}

Member Function Documentation

l1t::EMTFHit2016ExtraCollection PrimitiveConverterRPC::convert	(	std::vector< L1TMuon::TriggerPrimitive >	TriggPrim,
		int	SectIndex,
		edm::ESHandle< RPCGeometry >	rpc_geom
	)

Get all the input variables

Definition at line 13 of file PrimitiveConverterRPC.cc.

References funct::abs(), L1TMuon::TriggerPrimitive::RPCData::bx, gather_cfg::cout, L1TMuon::TriggerPrimitive::detId(), PV3DBase< T, PVType, FrameType >::eta(), L1TMuon::TriggerPrimitive::getRPCData(), l1t::EMTFHit2016::ImportRPCDetId(), l1t::EMTFHit2016::ImportRPCDigi(), L1TMuon::TriggerPrimitive::RPCData::layer, PV3DBase< T, PVType, FrameType >::phi(), Geom::pi(), RPCDetId::region(), RPCGeometry::roll(), l1t::EMTFHit2016::RPC_DetId(), sameRpcChamber(), RPCDetId::sector(), l1t::EMTFHit2016Extra::set_eta(), l1t::EMTFHit2016Extra::set_layer(), l1t::EMTFHit2016Extra::set_phi_glob_deg(), l1t::EMTFHit2016Extra::set_phi_glob_rad(), l1t::EMTFHit2016::set_sector_index(), l1t::EMTFHit2016::set_strip_hi(), l1t::EMTFHit2016::set_strip_low(), L1TMuon::TriggerPrimitive::RPCData::strip, l1t::EMTFHit2016::Strip_hi(), and l1t::EMTFHit2016::Strip_low().

                                                                   {

  // bool verbose = true;
  bool verbose = false;


  if (verbose) std::cout << "\n========== RPC Primitive Converter ==========" << std::endl;

  l1t::EMTFHit2016ExtraCollection tmpHits;
  for (std::vector<L1TMuon::TriggerPrimitive>::iterator iter = TrigPrim.begin(); iter != TrigPrim.end(); iter++) {

    L1TMuon::TriggerPrimitive prim = *iter; // Eventually would like to deprecate TriggerPrimitive entirely - AWB 03.06.16
    RPCDetId detID = prim.detId<RPCDetId>();
    RPCDigi digi = RPCDigi( prim.getRPCData().strip, prim.getRPCData().bx );

    // Only include RPC hits from correct sector in endcap
    if ( abs(detID.region()) != 1 ) continue;
    if ( SectIndex != (detID.sector() - 1) + (detID.region() == -1)*6 ) continue;

    l1t::EMTFHit2016Extra thisHit;
    thisHit.ImportRPCDetId( detID );
    thisHit.ImportRPCDigi( digi );
    thisHit.set_sector_index( SectIndex );
    thisHit.set_layer( prim.getRPCData().layer ); // In RE1/2 there are two layers of chambers: 1 is inner (front) and 2 is outer (rear)
    
    tmpHits.push_back( thisHit );
  }

  l1t::EMTFHit2016ExtraCollection clustHits;
  for (unsigned int iHit = 0; iHit < tmpHits.size(); iHit++) {
    l1t::EMTFHit2016Extra hit1 = tmpHits.at(iHit);

    // Skip hit if it is already in a cluster
    bool hit_in_cluster = false;
    for (unsigned int jHit = 0; jHit < clustHits.size(); jHit++) {
      l1t::EMTFHit2016Extra clustHit = clustHits.at(jHit);
      if ( sameRpcChamber(hit1, clustHit) && hit1.Strip_hi() <= clustHit.Strip_hi() && 
       hit1.Strip_low() >= clustHit.Strip_low() ) hit_in_cluster = true;
    }
    if (hit_in_cluster) continue;

    // Cluster adjascent strip hits.  Phi of cluster corresponds to central strip.
    // Does this clustering catch all the hits in large clusters? - AWB 03.06.16
    int prevHi = -999, prevLow = -999;
    while (hit1.Strip_hi() != prevHi || hit1.Strip_low() != prevLow) {
      prevHi = hit1.Strip_hi();
      prevLow = hit1.Strip_low();

      for (unsigned int jHit = 0; jHit < tmpHits.size(); jHit++) {
    if (iHit == jHit) continue;
    l1t::EMTFHit2016Extra hit2 = tmpHits.at(jHit);

    if (not sameRpcChamber(hit1, hit2)) continue;
    if (hit2.Strip_hi()  == hit1.Strip_hi()  + 1) hit1.set_strip_hi ( hit2.Strip_hi()  );
    if (hit2.Strip_low() == hit1.Strip_low() - 1) hit1.set_strip_low( hit2.Strip_low() );
      }
    }
    
    // Get phi and eta
    std::unique_ptr<const RPCRoll> roll(rpc_geom->roll( hit1.RPC_DetId() ));
    const LocalPoint lpHi = roll->centreOfStrip(hit1.Strip_hi());
    const GlobalPoint gpHi = roll->toGlobal(lpHi);
    const LocalPoint lpLow = roll->centreOfStrip(hit1.Strip_low());
    const GlobalPoint gpLow = roll->toGlobal(lpLow);
    roll.release();

    float glob_phi_hi  = gpHi.phi();  // Averaging global point phi's directly does weird things,
    float glob_phi_low = gpLow.phi(); // e.g. 2*pi/3 + 2*pi/3 = 4*pi/3 = -2*pi/3, so avg. = -pi/3, not 2*pi/3
    float glob_phi = (glob_phi_hi + glob_phi_low) / 2.0;
    float glob_eta = (gpHi.eta() + gpLow.eta()) / 2.0;

    if (verbose) std::cout << "RPC cluster phi = " << glob_phi << " (" << glob_phi_hi << ", " << glob_phi_low 
               << "), eta = " << glob_eta << " (" << gpHi.eta() << ", " << gpLow.eta() << ")" << std::endl;

    hit1.set_phi_glob_rad( glob_phi );
    hit1.set_phi_glob_deg( glob_phi*180/Geom::pi() );
    hit1.set_eta( glob_eta );
    clustHits.push_back( hit1 );

    // !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
    // Need to set phi_loc_int, theta_int, phi_hit, zone, csc_ID, quality,
    // pattern, wire, strip, phi_zone, and zone_contribution
    // !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
    
  } // End loop over iHit

  return clustHits;

} // End PrimitiveConverterRPC::convert

std::vector< ConvertedHit > PrimitiveConverterRPC::fillConvHits

(

l1t::EMTFHit2016ExtraCollection

exHits

)

Definition at line 106 of file PrimitiveConverterRPC.cc.

References l1t::EMTFHit2016::BX(), l1t::EMTFHit2016::CSC_ID(), l1t::EMTFHit2016Extra::Layer(), l1t::EMTFHit2016::Pattern(), l1t::EMTFHit2016Extra::Phi_hit(), l1t::EMTFHit2016Extra::Phi_loc_int(), l1t::EMTFHit2016Extra::Phi_zone(), l1t::EMTFHit2016::Quality(), l1t::EMTFHit2016::RPC_DetId(), l1t::EMTFHit2016::Sector_index(), ConvertedHit::SetNeighbor(), ConvertedHit::SetSectorIndex(), ConvertedHit::SetTP(), ConvertedHit::SetValues(), ConvertedHit::SetZhit(), l1t::EMTFHit2016::Station(), l1t::EMTFHit2016::Strip(), l1t::EMTFHit2016::Subsector(), l1t::EMTFHit2016Extra::Theta_int(), l1t::EMTFHit2016::Wire(), and l1t::EMTFHit2016Extra::Zone().

                                                                        {

  std::vector<ConvertedHit> ConvHits;
  for (unsigned int iHit = 0; iHit < exHits.size(); iHit++) {
    l1t::EMTFHit2016Extra exHit = exHits.at(iHit);

    // // Replace with SetZoneWord - AWB 04.09.16
    // std::vector<int> zone_contribution;

    ConvertedHit ConvHit;
    ConvHit.SetValues(exHit.Phi_loc_int(), exHit.Theta_int(), exHit.Phi_hit(), exHit.Zone(), 
              exHit.Station(), exHit.Subsector(), exHit.CSC_ID(), exHit.Quality(), 
              exHit.Pattern(), exHit.Wire(), exHit.Strip(), exHit.BX() + 6);
    ConvHit.SetTP( L1TMuon::TriggerPrimitive( exHit.RPC_DetId(), exHit.Strip(), exHit.Layer(), exHit.BX() + 6 ) );
    ConvHit.SetZhit( exHit.Phi_zone() );
    // // Replace with SetZoneWord - AWB 04.09.16
    // ConvHit.SetZoneContribution(zone_contribution);
    ConvHit.SetSectorIndex( exHit.Sector_index() );
    ConvHit.SetNeighbor(0);
    ConvHits.push_back(ConvHit);
  }
  return ConvHits;
}

bool PrimitiveConverterRPC::sameRpcChamber	(	l1t::EMTFHit2016Extra	hitA,
		l1t::EMTFHit2016Extra	hitB
	)

Definition at line 130 of file PrimitiveConverterRPC.cc.

References l1t::EMTFHit2016::Endcap(), l1t::EMTFHit2016Extra::Layer(), l1t::EMTFHit2016::Ring(), l1t::EMTFHit2016::Roll(), l1t::EMTFHit2016::Sector(), l1t::EMTFHit2016::Station(), and l1t::EMTFHit2016::Subsector().

Referenced by convert().

                                                                                               {

  if ( hitA.Endcap() == hitB.Endcap() && hitA.Station() == hitB.Station() && hitA.Ring() == hitB.Ring() &&
       hitA.Roll() == hitB.Roll() && hitA.Sector() == hitB.Sector() && hitA.Subsector() == hitB.Subsector() &&
       hitA.Layer() == hitB.Layer() ) return true;
  else return false;
}