#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

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 44 of file ME0TriggerPseudoBuilder.h.

44 {MAX_ENDCAPS = 2, MAX_CHAMBERS = 18};

Constructor & Destructor Documentation

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 ( )

Definition at line 20 of file ME0TriggerPseudoBuilder.cc.

21 {

22 }

Member Function Documentation

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 24 of file ME0TriggerPseudoBuilder.cc.

References relativeConstraints::cham, dumpAllME0Segments(), info_, ME0TriggerDigi::isValid(), LogTrace, MAX_CHAMBERS, MAX_ENDCAPS, and segmentConversion().

 {
   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.
       const int region(endc == 0 ? -1 : 1);
       ME0DetId detid(region, 0, cham, 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);
       }
     } 
   }
 }

void ME0TriggerPseudoBuilder::dumpAllME0Segments ( const ME0SegmentCollection & segments ) const

Definition at line 134 of file ME0TriggerPseudoBuilder.cc.

References PV3DBase< T, PVType, FrameType >::eta(), runTauDisplay::gp, ME0Geometry::idToDet(), LogTrace, me0_g, ME0RecHit::me0Id(), PV3DBase< T, PVType, FrameType >::phi(), and GeomDet::toGlobal().

Referenced by build(), and setME0Geometry().

 {
   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;
           }
       }
   }
 
 }

ME0TriggerDigi ME0TriggerPseudoBuilder::segmentConversion ( const ME0Segment segment )

private

Definition at line 60 of file ME0TriggerPseudoBuilder.cc.

References rpcdqm::BX, ME0EtaPartition::centreOfStrip(), ME0Geometry::chamber(), relativeConstraints::chamber, ME0Chamber::computeDeltaPhi(), 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, PV3DBase< T, PVType, FrameType >::phi(), jets_cff::quality, mps_fire::result, ME0TriggerDigi::setStrip(), ME0Segment::specificRecHits(), digitizers_cfi::strip, ME0EtaPartition::strip(), ME0Segment::time(), ntuplemaker::time, and GeomDet::toGlobal().

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 (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 );
   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();
   }
   
   const ME0EtaPartition* etapart = keylayer->etaPartition(rolls[0]);
   float strippitch = etapart->localPitch(segment.localPosition());
   float strip = etapart->strip(segment.localPosition());
   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;
   
   //globalpoint from ME0 segment
   GlobalPoint gp = me0_g->idToDet(segment.me0DetId())->surface().toGlobal(segment.localPosition());
   //gloablpoint from ME0 trigger digi
   float centreOfStrip = istrip + 0.25 +phiposition2*0.5;
   GlobalPoint gp_digi = etapart->toGlobal(etapart->centreOfStrip(centreOfStrip));
 
 
   int idphi = static_cast<int>(fabs(dphi)/(strippitch*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*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;
 }

void ME0TriggerPseudoBuilder::setME0Geometry ( const ME0Geometry * g )

inline

set geometry for the matching needs

Definition at line 36 of file ME0TriggerPseudoBuilder.h.

References dumpAllME0Segments(), g, and me0_g.

36 { me0_g = g; }

g

The Signals That Services Can Subscribe To This is based on ActivityRegistry and is current per Services can connect to the signals distributed by the ActivityRegistry in order to monitor the activity of the application Each possible callback has some defined which we here list in angle e g

Definition: Activities.doc:4

ME0TriggerPseudoBuilder::me0_g

const ME0Geometry * me0_g

Definition: ME0TriggerPseudoBuilder.h:55