L1TMuonBarrelKalmanSectorProcessor Class Reference

#include <L1TMuonBarrelKalmanSectorProcessor.h>

Classes
struct	bmtf_out
class	TrackSorter

Public Member Functions
	L1TMuonBarrelKalmanSectorProcessor (const edm::ParameterSet &, int sector)
L1MuKBMTrackCollection	process (L1TMuonBarrelKalmanAlgo *, const L1MuKBMTCombinedStubRefVector &stubs, int bx)
void	verbose (L1TMuonBarrelKalmanAlgo *, const L1MuKBMTrackCollection &)
	~L1TMuonBarrelKalmanSectorProcessor ()

Private Member Functions
L1MuKBMTrackCollection	cleanNeighbor (const L1MuKBMTrackCollection &, const L1MuKBMTrackCollection &)
L1MuKBMTrackCollection	cleanNeighbors (const L1MuKBMTrackCollection &, const L1MuKBMTrackCollection &, const L1MuKBMTrackCollection &)
bmtf_out	makeWord (L1TMuonBarrelKalmanAlgo *, const L1MuKBMTrackCollection &)
L1MuKBMTrackCollection	wedgeSort (const L1MuKBMTrackCollection &, const L1MuKBMTrackCollection &, const L1MuKBMTrackCollection &, const L1MuKBMTrackCollection &, const L1MuKBMTrackCollection &)

Private Attributes
std::vector< L1TMuonBarrelKalmanRegionModule >	regions_
int	sector_
int	verbose_

Detailed Description

Definition at line 10 of file L1TMuonBarrelKalmanSectorProcessor.h.

Constructor & Destructor Documentation

L1TMuonBarrelKalmanSectorProcessor::L1TMuonBarrelKalmanSectorProcessor	(	const edm::ParameterSet &	iConfig,
		int	sector
	)

Definition at line 5 of file L1TMuonBarrelKalmanSectorProcessor.cc.

References edm::ParameterSet::getParameter(), regions_, and makeMuonMisalignmentScenario::wheel.

                                                                                                               : 
  verbose_(iConfig.getParameter<int>("verbose")),
  sector_(sector)
{
  std::vector<int> wheels = iConfig.getParameter<std::vector<int> >("wheelsToProcess");
    for (const auto wheel : wheels)
      regions_.push_back(L1TMuonBarrelKalmanRegionModule(iConfig.getParameter<edm::ParameterSet>("regionSettings"),wheel,sector));
}

L1TMuonBarrelKalmanSectorProcessor::~L1TMuonBarrelKalmanSectorProcessor

(

)

Definition at line 16 of file L1TMuonBarrelKalmanSectorProcessor.cc.

{}

Member Function Documentation

L1MuKBMTrackCollection L1TMuonBarrelKalmanSectorProcessor::cleanNeighbor ( const L1MuKBMTrackCollection &  coll1, const L1MuKBMTrackCollection &  coll2  )

private

Definition at line 234 of file L1TMuonBarrelKalmanSectorProcessor.cc.

References createfilelist::int, keep, MillePedeFileConverter_cfg::out, and verbose_.

Referenced by wedgeSort().

                                                                                                                                                {
  L1MuKBMTrackCollection out;

 for(const auto& track1 : coll1) {
    if (verbose_)
      printf("Pre Track charge=%d pt=%f eta=%f phi=%f curvature=%d curvature STA =%d stubs=%d bitmask=%d rank=%d chi=%d pts=%f %f\n",track1.charge(),track1.pt(),track1.eta(),track1.phi(),track1.curvatureAtVertex(),track1.curvatureAtMuon(),int(track1.stubs().size()),track1.hitPattern(),track1.rank(),track1.approxChi2(),track1.pt(),track1.ptUnconstrained()); 

    bool keep=true;
    for(const auto& track2 : coll2) {
      if (!track1.overlapTrack(track2))
    continue;

      if (track1.rank()<track2.rank())
    keep=false;

      if ((track1.rank()==track2.rank()) && (fabs(track1.stubs()[0]->whNum())<fabs(track2.stubs()[0]->whNum())))
    keep=false;
    }
    if (keep) 
      out.push_back(track1);
 }

 return out;
}

L1MuKBMTrackCollection L1TMuonBarrelKalmanSectorProcessor::cleanNeighbors ( const L1MuKBMTrackCollection &  coll1, const L1MuKBMTrackCollection &  coll2, const L1MuKBMTrackCollection &  coll3  )

private

Definition at line 259 of file L1TMuonBarrelKalmanSectorProcessor.cc.

References keep, and MillePedeFileConverter_cfg::out.

Referenced by wedgeSort().

                                                                                                                                                                                     {
  L1MuKBMTrackCollection out;

 for(const auto& track1 : coll1) {
    // if (verbose_)
    //   printf("Pre Track charge=%d pt=%f eta=%f phi=%f curvature=%d curvature STA =%d stubs=%d bitmask=%d rank=%d chi=%d pts=%f %f\n",track1.charge(),track1.pt(),track1.eta(),track1.phi(),track1.curvatureAtVertex(),track1.curvatureAtMuon(),int(track1.stubs().size()),track1.hitPattern(),track1.rank(),track1.approxChi2(),track1.pt(),track1.ptUnconstrained()); 

    bool keep=true;
    for(const auto& track2 : coll2) {
      if (!track1.overlapTrack(track2))
    continue;

      if (track1.rank()<track2.rank())
    keep=false;

      if ((track1.rank()==track2.rank()) && (fabs(track1.stubs()[0]->whNum())<fabs(track2.stubs()[0]->whNum())))
    keep=false;
    }

    for(const auto& track2 : coll3) {
      if (!track1.overlapTrack(track2))
    continue;

      if (track1.rank()<track2.rank())
    keep=false;

      if ((track1.rank()==track2.rank()) && (fabs(track1.stubs()[0]->whNum())<fabs(track2.stubs()[0]->whNum())))
    keep=false;
    }

    if (keep) 
      out.push_back(track1);
 }

 return out;
}

L1TMuonBarrelKalmanSectorProcessor::bmtf_out L1TMuonBarrelKalmanSectorProcessor::makeWord ( L1TMuonBarrelKalmanAlgo *  trackMaker, const L1MuKBMTrackCollection &  tracks  )

private

Definition at line 53 of file L1TMuonBarrelKalmanSectorProcessor.cc.

References L1TMuonBarrelKalmanSectorProcessor::bmtf_out::addr1_1, L1TMuonBarrelKalmanSectorProcessor::bmtf_out::addr1_2, L1TMuonBarrelKalmanSectorProcessor::bmtf_out::addr1_3, L1TMuonBarrelKalmanSectorProcessor::bmtf_out::addr2_1, L1TMuonBarrelKalmanSectorProcessor::bmtf_out::addr2_2, L1TMuonBarrelKalmanSectorProcessor::bmtf_out::addr2_3, L1TMuonBarrelKalmanSectorProcessor::bmtf_out::addr3_1, L1TMuonBarrelKalmanSectorProcessor::bmtf_out::addr3_2, L1TMuonBarrelKalmanSectorProcessor::bmtf_out::addr3_3, L1TMuonBarrelKalmanSectorProcessor::bmtf_out::addr4_1, L1TMuonBarrelKalmanSectorProcessor::bmtf_out::addr4_2, L1TMuonBarrelKalmanSectorProcessor::bmtf_out::addr4_3, L1TMuonBarrelKalmanSectorProcessor::bmtf_out::bx0_1, L1TMuonBarrelKalmanSectorProcessor::bmtf_out::bx0_2, L1TMuonBarrelKalmanSectorProcessor::bmtf_out::bx0_3, L1TMuonBarrelKalmanSectorProcessor::bmtf_out::charge_1, L1TMuonBarrelKalmanSectorProcessor::bmtf_out::charge_2, L1TMuonBarrelKalmanSectorProcessor::bmtf_out::charge_3, L1TMuonBarrelKalmanSectorProcessor::bmtf_out::chargeValid_1, L1TMuonBarrelKalmanSectorProcessor::bmtf_out::chargeValid_2, L1TMuonBarrelKalmanSectorProcessor::bmtf_out::chargeValid_3, L1TMuonBarrelKalmanAlgo::convertToBMTF(), L1TMuonBarrelKalmanSectorProcessor::bmtf_out::dxy_1, L1TMuonBarrelKalmanSectorProcessor::bmtf_out::dxy_2, L1TMuonBarrelKalmanSectorProcessor::bmtf_out::dxy_3, L1TMuonBarrelKalmanSectorProcessor::bmtf_out::eta_1, L1TMuonBarrelKalmanSectorProcessor::bmtf_out::eta_2, L1TMuonBarrelKalmanSectorProcessor::bmtf_out::eta_3, L1TMuonBarrelKalmanSectorProcessor::bmtf_out::HF_1, L1TMuonBarrelKalmanSectorProcessor::bmtf_out::HF_2, L1TMuonBarrelKalmanSectorProcessor::bmtf_out::HF_3, l1t::RegionalMuonCand::hwDXY(), l1t::RegionalMuonCand::hwEta(), l1t::RegionalMuonCand::hwHF(), l1t::RegionalMuonCand::hwPhi(), l1t::RegionalMuonCand::hwPt(), l1t::RegionalMuonCand::hwPt2(), l1t::RegionalMuonCand::hwQual(), l1t::RegionalMuonCand::hwSign(), l1t::RegionalMuonCand::hwSignValid(), l1t::RegionalMuonCand::kStat1, l1t::RegionalMuonCand::kStat2, l1t::RegionalMuonCand::kStat3, l1t::RegionalMuonCand::kStat4, l1t::RegionalMuonCand::kWheelNum, l1t::RegionalMuonCand::kWheelSide, RPCpg::mu, MillePedeFileConverter_cfg::out, L1TMuonBarrelKalmanSectorProcessor::bmtf_out::phi_1, L1TMuonBarrelKalmanSectorProcessor::bmtf_out::phi_2, L1TMuonBarrelKalmanSectorProcessor::bmtf_out::phi_3, L1TMuonBarrelKalmanSectorProcessor::bmtf_out::pt_1, L1TMuonBarrelKalmanSectorProcessor::bmtf_out::pt_2, L1TMuonBarrelKalmanSectorProcessor::bmtf_out::pt_3, L1TMuonBarrelKalmanSectorProcessor::bmtf_out::ptSTA_1, L1TMuonBarrelKalmanSectorProcessor::bmtf_out::ptSTA_2, L1TMuonBarrelKalmanSectorProcessor::bmtf_out::ptSTA_3, L1TMuonBarrelKalmanSectorProcessor::bmtf_out::qual_1, L1TMuonBarrelKalmanSectorProcessor::bmtf_out::qual_2, L1TMuonBarrelKalmanSectorProcessor::bmtf_out::qual_3, L1TMuonBarrelKalmanSectorProcessor::bmtf_out::reserved_1, L1TMuonBarrelKalmanSectorProcessor::bmtf_out::reserved_2, L1TMuonBarrelKalmanSectorProcessor::bmtf_out::reserved_3, L1TMuonBarrelKalmanSectorProcessor::bmtf_out::SE_1, L1TMuonBarrelKalmanSectorProcessor::bmtf_out::SE_2, L1TMuonBarrelKalmanSectorProcessor::bmtf_out::SE_3, l1t::RegionalMuonCand::trackSubAddress(), L1TMuonBarrelKalmanSectorProcessor::bmtf_out::wheel_1, L1TMuonBarrelKalmanSectorProcessor::bmtf_out::wheel_2, and L1TMuonBarrelKalmanSectorProcessor::bmtf_out::wheel_3.

Referenced by verbose().

                                                                                                                                                                 {
  L1TMuonBarrelKalmanSectorProcessor::bmtf_out out;
    out.pt_1=0;
    out.qual_1=0;
    out.eta_1=0;
    out.HF_1=0;
    out.phi_1=0;
    out.bx0_1=0;
    out.charge_1=0;
    out.chargeValid_1=0;
    out.dxy_1=0;
    out.addr1_1=3;
    out.addr2_1=15;
    out.addr3_1=15;
    out.addr4_1=15;
    out.reserved_1=0;
    out.wheel_1=0;
    out.ptSTA_1=0;
    out.SE_1=0;

    out.pt_2=0;
    out.qual_2=0;
    out.eta_2=0;
    out.HF_2=0;
    out.phi_2=0;
    out.bx0_2=0;
    out.charge_2=0;
    out.chargeValid_2=0;
    out.dxy_2=0;
    out.addr1_2=3;
    out.addr2_2=15;
    out.addr3_2=15;
    out.addr4_2=15;
    out.reserved_2=0;
    out.wheel_2=0;
    out.ptSTA_2=0;
    out.SE_2=0;

    out.pt_3=0;
    out.qual_3=0;
    out.eta_3=0;
    out.HF_3=0;
    out.phi_3=0;
    out.bx0_3=0;
    out.charge_3=0;
    out.chargeValid_3=0;
    out.dxy_3=0;
    out.addr1_3=3;
    out.addr2_3=15;
    out.addr3_3=15;
    out.addr4_3=15;
    out.reserved_3=0;
    out.wheel_3=0;
    out.ptSTA_3=0;
    out.SE_3=0;

    if (!tracks.empty()) {
      l1t::RegionalMuonCand mu = trackMaker->convertToBMTF(tracks[0]);
      out.pt_1=mu.hwPt();
      out.qual_1=mu.hwQual();
      out.eta_1=mu.hwEta();
      out.HF_1=mu.hwHF();
      out.phi_1=mu.hwPhi();
      out.bx0_1=0;
      out.charge_1=mu.hwSign();
      out.chargeValid_1=mu.hwSignValid();
      out.dxy_1=mu.hwDXY();
      out.addr1_1=mu.trackSubAddress(l1t::RegionalMuonCand::kStat1);
      out.addr2_1=mu.trackSubAddress(l1t::RegionalMuonCand::kStat2);
      out.addr3_1=mu.trackSubAddress(l1t::RegionalMuonCand::kStat3);
      out.addr4_1=mu.trackSubAddress(l1t::RegionalMuonCand::kStat4);
      out.wheel_1 = mu.trackSubAddress(l1t::RegionalMuonCand::kWheelSide)*(1<<2)+mu.trackSubAddress(l1t::RegionalMuonCand::kWheelNum);
      out.ptSTA_1=mu.hwPt2();
    }

    if (tracks.size()>1) {
      l1t::RegionalMuonCand mu = trackMaker->convertToBMTF(tracks[1]);
      out.pt_2=mu.hwPt();
      out.qual_2=mu.hwQual();
      out.eta_2=mu.hwEta();
      out.HF_2=mu.hwHF();
      out.phi_2=mu.hwPhi();
      out.bx0_2=0;
      out.charge_2=mu.hwSign();
      out.chargeValid_2=mu.hwSignValid();
      out.dxy_2=mu.hwDXY();
      out.addr1_2=mu.trackSubAddress(l1t::RegionalMuonCand::kStat1);
      out.addr2_2=mu.trackSubAddress(l1t::RegionalMuonCand::kStat2);
      out.addr3_2=mu.trackSubAddress(l1t::RegionalMuonCand::kStat3);
      out.addr4_2=mu.trackSubAddress(l1t::RegionalMuonCand::kStat4);
      out.wheel_2 = mu.trackSubAddress(l1t::RegionalMuonCand::kWheelSide)*(1<<2)+mu.trackSubAddress(l1t::RegionalMuonCand::kWheelNum);


      out.ptSTA_2=mu.hwPt2();
    }

    if (tracks.size()>2) {
      l1t::RegionalMuonCand mu = trackMaker->convertToBMTF(tracks[2]);
      out.pt_3=mu.hwPt();
      out.qual_3=mu.hwQual();
      out.eta_3=mu.hwEta();
      out.HF_3=mu.hwHF();
      out.phi_3=mu.hwPhi();
      out.bx0_3=0;
      out.charge_3=mu.hwSign();
      out.chargeValid_3=mu.hwSignValid();
      out.dxy_3=mu.hwDXY();
      out.addr1_3=mu.trackSubAddress(l1t::RegionalMuonCand::kStat1);
      out.addr2_3=mu.trackSubAddress(l1t::RegionalMuonCand::kStat2);
      out.addr3_3=mu.trackSubAddress(l1t::RegionalMuonCand::kStat3);
      out.addr4_3=mu.trackSubAddress(l1t::RegionalMuonCand::kStat4);
      out.wheel_3 = mu.trackSubAddress(l1t::RegionalMuonCand::kWheelSide)*(1<<2)+mu.trackSubAddress(l1t::RegionalMuonCand::kWheelNum);
      out.ptSTA_3=mu.hwPt2();
    }
    return out;
  
}

L1MuKBMTrackCollection L1TMuonBarrelKalmanSectorProcessor::process	(	L1TMuonBarrelKalmanAlgo *	trackMaker,
		const L1MuKBMTCombinedStubRefVector &	stubs,
		int	bx
	)

Definition at line 18 of file L1TMuonBarrelKalmanSectorProcessor.cc.

References MillePedeFileConverter_cfg::out, regions_, tmp, verbose(), verbose_, and wedgeSort().

                                                                                                                                                            {


  L1MuKBMTrackCollection tracksM2;
  L1MuKBMTrackCollection tracksM1;
  L1MuKBMTrackCollection tracks0;
  L1MuKBMTrackCollection tracksP1;
  L1MuKBMTrackCollection tracksP2;

  for (auto& region: regions_) {
    L1MuKBMTrackCollection tmp = region.process(trackMaker,stubsAll,bx);
    if (region.wheel()==-2)
      tracksM2.insert(tracksM2.end(),tmp.begin(),tmp.end());
    if (region.wheel()==-1)
      tracksM1.insert(tracksM1.end(),tmp.begin(),tmp.end());
    if (region.wheel()==0)
      tracks0.insert(tracks0.end(),tmp.begin(),tmp.end());
    if (region.wheel()==1)
      tracksP1.insert(tracksP1.end(),tmp.begin(),tmp.end());
    if (region.wheel()==2)
      tracksP2.insert(tracksP2.end(),tmp.begin(),tmp.end());
   
  } 

  L1MuKBMTrackCollection out =wedgeSort(tracksM2,tracksM1,tracks0,tracksP1,tracksP2);
  if (verbose_==1)
    verbose(trackMaker,out);

  return out;

}

void L1TMuonBarrelKalmanSectorProcessor::verbose	(	L1TMuonBarrelKalmanAlgo *	algo,
		const L1MuKBMTrackCollection &	tracks
	)

Definition at line 175 of file L1TMuonBarrelKalmanSectorProcessor.cc.

References L1TMuonBarrelKalmanSectorProcessor::bmtf_out::addr1_1, L1TMuonBarrelKalmanSectorProcessor::bmtf_out::addr1_2, L1TMuonBarrelKalmanSectorProcessor::bmtf_out::addr1_3, L1TMuonBarrelKalmanSectorProcessor::bmtf_out::addr2_1, L1TMuonBarrelKalmanSectorProcessor::bmtf_out::addr2_2, L1TMuonBarrelKalmanSectorProcessor::bmtf_out::addr2_3, L1TMuonBarrelKalmanSectorProcessor::bmtf_out::addr3_1, L1TMuonBarrelKalmanSectorProcessor::bmtf_out::addr3_2, L1TMuonBarrelKalmanSectorProcessor::bmtf_out::addr3_3, L1TMuonBarrelKalmanSectorProcessor::bmtf_out::addr4_1, L1TMuonBarrelKalmanSectorProcessor::bmtf_out::addr4_2, L1TMuonBarrelKalmanSectorProcessor::bmtf_out::addr4_3, L1TMuonBarrelKalmanSectorProcessor::bmtf_out::charge_1, L1TMuonBarrelKalmanSectorProcessor::bmtf_out::charge_2, L1TMuonBarrelKalmanSectorProcessor::bmtf_out::charge_3, L1TMuonBarrelKalmanSectorProcessor::bmtf_out::chargeValid_1, L1TMuonBarrelKalmanSectorProcessor::bmtf_out::chargeValid_2, L1TMuonBarrelKalmanSectorProcessor::bmtf_out::chargeValid_3, gather_cfg::cout, L1TMuonBarrelKalmanSectorProcessor::bmtf_out::dxy_1, L1TMuonBarrelKalmanSectorProcessor::bmtf_out::dxy_2, L1TMuonBarrelKalmanSectorProcessor::bmtf_out::dxy_3, L1TMuonBarrelKalmanSectorProcessor::bmtf_out::eta_1, L1TMuonBarrelKalmanSectorProcessor::bmtf_out::eta_2, L1TMuonBarrelKalmanSectorProcessor::bmtf_out::eta_3, L1TMuonBarrelKalmanSectorProcessor::bmtf_out::HF_1, L1TMuonBarrelKalmanSectorProcessor::bmtf_out::HF_2, L1TMuonBarrelKalmanSectorProcessor::bmtf_out::HF_3, makeWord(), MillePedeFileConverter_cfg::out, L1TMuonBarrelKalmanSectorProcessor::bmtf_out::phi_1, L1TMuonBarrelKalmanSectorProcessor::bmtf_out::phi_2, L1TMuonBarrelKalmanSectorProcessor::bmtf_out::phi_3, L1TMuonBarrelKalmanSectorProcessor::bmtf_out::pt_1, L1TMuonBarrelKalmanSectorProcessor::bmtf_out::pt_2, L1TMuonBarrelKalmanSectorProcessor::bmtf_out::pt_3, L1TMuonBarrelKalmanSectorProcessor::bmtf_out::ptSTA_1, L1TMuonBarrelKalmanSectorProcessor::bmtf_out::ptSTA_2, L1TMuonBarrelKalmanSectorProcessor::bmtf_out::ptSTA_3, L1TMuonBarrelKalmanSectorProcessor::bmtf_out::qual_1, L1TMuonBarrelKalmanSectorProcessor::bmtf_out::qual_2, L1TMuonBarrelKalmanSectorProcessor::bmtf_out::qual_3, sector_, L1TMuonBarrelKalmanSectorProcessor::bmtf_out::wheel_1, L1TMuonBarrelKalmanSectorProcessor::bmtf_out::wheel_2, and L1TMuonBarrelKalmanSectorProcessor::bmtf_out::wheel_3.

Referenced by process().

                                                                                                                   {
  L1TMuonBarrelKalmanSectorProcessor::bmtf_out out = makeWord(algo,tracks);
    std::cout << "O "<<sector_<<" "<< out.pt_1 << " " << out.qual_1 << " " << out.eta_1 << " " << out.HF_1 << " " << out.phi_1 << " " << out.charge_1 << " " << out.chargeValid_1 << " " << out.dxy_1 <<" " << out.addr1_1<< " " << out.addr2_1 <<" " << out.addr3_1 << " " << out.addr4_1 << " " << out.wheel_1 << " " << out.ptSTA_1 << " " << out.pt_2 << " " << out.qual_2 << " " << out.eta_2 << " " << out.HF_2 << " " << out.phi_2 << " " << out.charge_2 << " " << out.chargeValid_2 << " " << out.dxy_2 <<" " << out.addr1_2<< " " << out.addr2_2 <<" " << out.addr3_2 << " " << out.addr4_2 << " " << out.wheel_2 << " " << out.ptSTA_2 << " " << out.pt_3 << " " << out.qual_3 << " " << out.eta_3 << " " << out.HF_3 << " " << out.phi_3 << " " << out.charge_3 << " " << out.chargeValid_3 << " " << out.dxy_3 <<" " << out.addr1_3<< " " << out.addr2_3 <<" " << out.addr3_3 << " " << out.addr4_3 << " " << out.wheel_3 << " " << out.ptSTA_3 << std::endl;


}

L1MuKBMTrackCollection L1TMuonBarrelKalmanSectorProcessor::wedgeSort ( const L1MuKBMTrackCollection &  preminus2, const L1MuKBMTrackCollection &  preminus1, const L1MuKBMTrackCollection &  prezero, const L1MuKBMTrackCollection &  preone, const L1MuKBMTrackCollection &  pretwo  )

private

Definition at line 321 of file L1TMuonBarrelKalmanSectorProcessor.cc.

References cleanNeighbor(), cleanNeighbors(), mps_fire::i, MillePedeFileConverter_cfg::out, swap(), and parallelization::uint().

Referenced by process().

                                                                                                                                                                                                                                                                    {

  //first clean

  L1MuKBMTrackCollection minus2 = cleanNeighbor(preminus2,preminus1);
  L1MuKBMTrackCollection minus1 = cleanNeighbors(preminus1,preminus2,prezero);
  L1MuKBMTrackCollection zero   = cleanNeighbors(prezero,preminus1,preone);
  L1MuKBMTrackCollection plus1   = cleanNeighbors(preone,prezero,pretwo);
  L1MuKBMTrackCollection plus2   = cleanNeighbor(pretwo,preone);

  std::map<uint,L1MuKBMTrack> trackInfo;
  std::map<uint,double> ptInfo;




  for (uint i=0;i<10;++i)
    ptInfo[i]=0.0;


  if (minus2.size()>1) {
    ptInfo[0] = minus2[0].pt();
    trackInfo[0] = minus2[0];
    ptInfo[1] = minus2[1].pt();
    trackInfo[1] = minus2[1];
  }
  else if (minus2.size()==1) {
    if (minus2[0].stubs()[0]->tag()) {
      ptInfo[1] = minus2[0].pt();
      trackInfo[1] = minus2[0];
    }
    else {
      ptInfo[0] = minus2[0].pt();
      trackInfo[0] = minus2[0];
    }
  }

  if (minus1.size()>1) {
    ptInfo[2] = minus1[0].pt();
    trackInfo[2] = minus1[0];
    ptInfo[3] = minus1[1].pt();
    trackInfo[3] = minus1[1];
  }
  else if (minus1.size()==1) {
    if (minus1[0].stubs()[0]->tag()) {
      ptInfo[3] = minus1[0].pt();
      trackInfo[3] = minus1[0];
    }
    else {
      ptInfo[2] = minus1[0].pt();
      trackInfo[2] = minus1[0];
    }
  }

  if (zero.size()>1) {
    ptInfo[4] = zero[0].pt();
    trackInfo[4] = zero[0];
    ptInfo[5] = zero[1].pt();
    trackInfo[5] = zero[1];
  }
  else if (zero.size()==1) {
    if (zero[0].stubs()[0]->tag()) {
      ptInfo[5] = zero[0].pt();
      trackInfo[5] = zero[0];
    }
    else {
      ptInfo[4] = zero[0].pt();
      trackInfo[4] = zero[0];
    }
  }


  if (plus1.size()>1) {
    ptInfo[6] = plus1[0].pt();
    trackInfo[6] = plus1[0];
    ptInfo[7] = plus1[1].pt();
    trackInfo[7] = plus1[1];
  }
  else if (plus1.size()==1) {
    if (plus1[0].stubs()[0]->tag()) {
      ptInfo[7] = plus1[0].pt();
      trackInfo[7] = plus1[0];
    }
    else {
      ptInfo[6] = plus1[0].pt();
      trackInfo[6] = plus1[0];
    }
  }

  if (plus2.size()>1) {
    ptInfo[8] = plus2[0].pt();
    trackInfo[8] = plus2[0];
    ptInfo[9] = plus2[1].pt();
    trackInfo[9] = plus2[1];
  }
  else if (plus2.size()==1) {
    if (plus2[0].stubs()[0]->tag()) {
      ptInfo[9] = plus2[0].pt();
      trackInfo[9] = plus2[0];
    }
    else {
      ptInfo[8] = plus2[0].pt();
      trackInfo[8] = plus2[0];
    }
  }




    //Now my glorious partial bitonic != power of two sorter 
    swap(ptInfo,trackInfo,0,5);
    swap(ptInfo,trackInfo,1,6);
    swap(ptInfo,trackInfo,2,7);
    swap(ptInfo,trackInfo,3,8);
    swap(ptInfo,trackInfo,4,9);

    swap(ptInfo,trackInfo,0,3);
    swap(ptInfo,trackInfo,1,4);
    swap(ptInfo,trackInfo,5,8);
    swap(ptInfo,trackInfo,6,9);

    swap(ptInfo,trackInfo,0,2);
    swap(ptInfo,trackInfo,3,6);
    swap(ptInfo,trackInfo,7,9);

    swap(ptInfo,trackInfo,0,1);
    swap(ptInfo,trackInfo,2,4);
    swap(ptInfo,trackInfo,5,7);
    swap(ptInfo,trackInfo,8,9);


    swap(ptInfo,trackInfo,1,2);
    swap(ptInfo,trackInfo,3,5);
    swap(ptInfo,trackInfo,4,6);
    swap(ptInfo,trackInfo,7,8);

    swap(ptInfo,trackInfo,1,3);
    swap(ptInfo,trackInfo,2,5);
    swap(ptInfo,trackInfo,4,7);
    swap(ptInfo,trackInfo,6,8);

    swap(ptInfo,trackInfo,2,3);


    
    L1MuKBMTrackCollection out;
    if (ptInfo[0]!=0.0)
      out.push_back(trackInfo[0]);
    if (ptInfo[1]!=0.0)
      out.push_back(trackInfo[1]);
    if (ptInfo[2]!=0.0)
      out.push_back(trackInfo[2]);

    return out;
}

Member Data Documentation

std::vector<L1TMuonBarrelKalmanRegionModule> L1TMuonBarrelKalmanSectorProcessor::regions_

private

Definition at line 38 of file L1TMuonBarrelKalmanSectorProcessor.h.

Referenced by L1TMuonBarrelKalmanSectorProcessor(), and process().

int L1TMuonBarrelKalmanSectorProcessor::sector_

private

Definition at line 36 of file L1TMuonBarrelKalmanSectorProcessor.h.

Referenced by verbose().

int L1TMuonBarrelKalmanSectorProcessor::verbose_

private

Definition at line 35 of file L1TMuonBarrelKalmanSectorProcessor.h.

Referenced by cleanNeighbor(), and process().