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::L1TMuonBarrelKalmanSectorProcessor	(	const edm::ParameterSet &	iConfig,
		int	sector
	)

Definition at line 3 of file L1TMuonBarrelKalmanSectorProcessor.cc.

References edm::ParameterSet::getParameter(), regions_, nano_mu_digi_cff::sector, makeMuonMisalignmentScenario::wheel, and twinMuxStage2Digis_cfi::wheels.

    : 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::~L1TMuonBarrelKalmanSectorProcessor

(

)

Definition at line 11 of file L1TMuonBarrelKalmanSectorProcessor.cc.

{}

Member Function Documentation

cleanNeighbor()

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

private

Definition at line 221 of file L1TMuonBarrelKalmanSectorProcessor.cc.

References runTheMatrix::keep, MillePedeFileConverter_cfg::out, and groupFilesInBlocks::temp.

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);
    else {
      L1MuKBMTrack temp = track1;
      temp.setPtEtaPhi(0, 0, 0);
      temp.setRank(0);
      out.push_back(temp);
    }
  }

  return out;
}

cleanNeighbors()

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

private

Definition at line 268 of file L1TMuonBarrelKalmanSectorProcessor.cc.

References runTheMatrix::keep, MillePedeFileConverter_cfg::out, and groupFilesInBlocks::temp.

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);
    else {
      L1MuKBMTrack temp = track1;
      temp.setPtEtaPhi(0, 0, 0);
      temp.setRank(0);
      out.push_back(temp);
    }
  }

  return out;
}

makeWord()

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

private

Definition at line 43 of file L1TMuonBarrelKalmanSectorProcessor.cc.

References L1TMuonBarrelKalmanAlgo::convertToBMTF(), l1t::RegionalMuonCand::kStat1, l1t::RegionalMuonCand::kStat2, l1t::RegionalMuonCand::kStat3, l1t::RegionalMuonCand::kStat4, l1t::RegionalMuonCand::kWheelNum, l1t::RegionalMuonCand::kWheelSide, amptDefaultParameters_cff::mu, MillePedeFileConverter_cfg::out, and DiMuonV_cfg::tracks.

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.hwPtUnconstrained();
  }

  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.hwPtUnconstrained();
  }

  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.hwPtUnconstrained();
  }
  return out;
}

process()

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

Definition at line 13 of file L1TMuonBarrelKalmanSectorProcessor.cc.

References nano_mu_digi_cff::bx, MillePedeFileConverter_cfg::out, nano_mu_digi_cff::region, regions_, createJobs::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;
}

verbose()

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

Definition at line 163 of file L1TMuonBarrelKalmanSectorProcessor.cc.

References gather_cfg::cout, makeWord(), MillePedeFileConverter_cfg::out, sector_, and DiMuonV_cfg::tracks.

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

wedgeSort()

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

private

Definition at line 332 of file L1TMuonBarrelKalmanSectorProcessor.cc.

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

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

regions_

std::vector<L1TMuonBarrelKalmanRegionModule> L1TMuonBarrelKalmanSectorProcessor::regions_

private

Definition at line 33 of file L1TMuonBarrelKalmanSectorProcessor.h.

Referenced by L1TMuonBarrelKalmanSectorProcessor(), and process().

sector_

int L1TMuonBarrelKalmanSectorProcessor::sector_

private

Definition at line 31 of file L1TMuonBarrelKalmanSectorProcessor.h.

Referenced by verbose().

verbose_

int L1TMuonBarrelKalmanSectorProcessor::verbose_

private

Definition at line 30 of file L1TMuonBarrelKalmanSectorProcessor.h.

Referenced by process().