MuonCSCSeedFromRecHits Class Reference

#include <MuonCSCSeedFromRecHits.h>

Inheritance diagram for MuonCSCSeedFromRecHits:

Public Member Functions
ConstMuonRecHitPointer	bestEndcapHit (const MuonRecHitContainer &endcapHits) const
	MuonCSCSeedFromRecHits ()
virtual TrajectorySeed	seed () const
	~MuonCSCSeedFromRecHits () override
Public Member Functions inherited from MuonSeedFromRecHits
void	add (MuonTransientTrackingRecHit::MuonRecHitPointer hit)
void	clear ()
TrajectorySeed	createSeed (float ptmean, float sptmean, MuonTransientTrackingRecHit::ConstMuonRecHitPointer last) const
MuonTransientTrackingRecHit::ConstMuonRecHitPointer	firstRecHit () const
	MuonSeedFromRecHits ()
unsigned int	nrhit () const
void	setBField (const MagneticField *field)
void	setPtExtractor (const MuonSeedPtExtractor *extractor)
virtual	~MuonSeedFromRecHits ()

Private Member Functions
void	analyze () const
float	computeDefaultPt (ConstMuonRecHitPointer muon) const
bool	createDefaultEndcapSeed (ConstMuonRecHitPointer last, TrajectorySeed &seed) const
void	makeDefaultSeed (TrajectorySeed &seed) const
bool	makeSeed (const MuonRecHitContainer &hits1, const MuonRecHitContainer &hits2, TrajectorySeed &seed) const
bool	makeSeed2 (const MuonRecHitContainer &hits1, const MuonRecHitContainer &hits2, TrajectorySeed &seed) const
int	segmentQuality (ConstMuonRecHitPointer muon) const

Additional Inherited Members
Protected Types inherited from MuonSeedFromRecHits
typedef MuonTransientTrackingRecHit::ConstMuonRecHitPointer	ConstMuonRecHitPointer
typedef MuonTransientTrackingRecHit::MuonRecHitContainer	MuonRecHitContainer
typedef MuonTransientTrackingRecHit::MuonRecHitPointer	MuonRecHitPointer
Protected Attributes inherited from MuonSeedFromRecHits
const MagneticField *	theField
const MuonSeedPtExtractor *	thePtExtractor
MuonTransientTrackingRecHit::MuonRecHitContainer	theRhits

Detailed Description

Definition at line 6 of file MuonCSCSeedFromRecHits.h.

Constructor & Destructor Documentation

MuonCSCSeedFromRecHits()

MuonCSCSeedFromRecHits::MuonCSCSeedFromRecHits

(

)

Definition at line 15 of file MuonCSCSeedFromRecHits.cc.

: MuonSeedFromRecHits() {}

~MuonCSCSeedFromRecHits()

MuonCSCSeedFromRecHits::~MuonCSCSeedFromRecHits ( )

inlineoverride

Definition at line 9 of file MuonCSCSeedFromRecHits.h.

{}

Member Function Documentation

analyze()

void MuonCSCSeedFromRecHits::analyze ( ) const

private

Definition at line 224 of file MuonCSCSeedFromRecHits.cc.

References gather_cfg::cout, SiPixelRawToDigiRegional_cfi::deltaPhi, and MuonSeedFromRecHits::theRhits.

                                           {
  for (MuonRecHitContainer::const_iterator iter = theRhits.begin(); iter != theRhits.end(); iter++) {
    if (!(*iter)->isCSC())
      continue;
    for (MuonRecHitContainer::const_iterator iter2 = iter + 1; iter2 != theRhits.end(); ++iter2) {
      if (!(*iter2)->isCSC())
        continue;

      CSCDetId cscId1((*iter)->geographicalId().rawId());
      CSCDetId cscId2((*iter2)->geographicalId().rawId());
      double dphi = deltaPhi((**iter).globalPosition().barePhi(), (**iter2).globalPosition().barePhi());

      int type1 = cscId1.iChamberType();
      int type2 = cscId2.iChamberType();

      // say the lower station first
      if (type1 < type2) {
        std::cout << "HITPAIRA," << type1 << type2 << "," << dphi << "," << (**iter).globalPosition().eta()
                  << std::endl;
      }
      if (type2 < type1) {
        std::cout << "HITPAIRB," << type2 << type1 << "," << -dphi << "," << (**iter2).globalPosition().eta()
                  << std::endl;
      }
      if (type1 == type2) {
        std::cout << "HITPAIRSAMESTATION," << type1 << cscId1.ring() << cscId2.ring() << "," << dphi << ","
                  << (**iter).globalPosition().eta() << std::endl;
      }
    }
  }
}

bestEndcapHit()

MuonCSCSeedFromRecHits::ConstMuonRecHitPointer MuonCSCSeedFromRecHits::bestEndcapHit

(

const MuonRecHitContainer &

endcapHits

)

const

Definition at line 154 of file MuonCSCSeedFromRecHits.cc.

References SiPixelRawToDigiRegional_cfi::deltaPhi, HLT_2023v12_cff::endcapHits, M_PI, hlt_dqm_clientPB-live_cfg::me1, quality, and segmentQuality().

Referenced by MuonOverlapSeedFromRecHits::bestHit(), makeDefaultSeed(), and makeSeed().

                                                 {
  MuonRecHitPointer me1 = nullptr, meit = nullptr;
  float dPhiGloDir = .0;          //  +v
  float bestdPhiGloDir = M_PI;    //  +v
  int quality1 = 0, quality = 0;  //  +v  I= 5,6-p. / II= 4p.  / III= 3p.

  for (MuonRecHitContainer::const_iterator iter = endcapHits.begin(); iter != endcapHits.end(); iter++) {
    if (!(*iter)->isCSC() && !(*iter)->isME0())
      continue;

    // tmp compar. Glob-Dir for the same tr-segm:

    meit = *iter;

    quality = segmentQuality(meit);

    dPhiGloDir = fabs(deltaPhi(meit->globalPosition().barePhi(), meit->globalDirection().barePhi()));

    if (!me1) {
      me1 = meit;
      quality1 = quality;
      bestdPhiGloDir = dPhiGloDir;
    }

    if (me1) {
      if (!me1->isValid()) {
        me1 = meit;
        quality1 = quality;
        bestdPhiGloDir = dPhiGloDir;
      }

      if (me1->isValid() && quality < quality1) {
        me1 = meit;
        quality1 = quality;
        bestdPhiGloDir = dPhiGloDir;
      }

      if (me1->isValid() && bestdPhiGloDir > .03) {
        if (dPhiGloDir < bestdPhiGloDir - .01 && quality == quality1) {
          me1 = meit;
          quality1 = quality;
          bestdPhiGloDir = dPhiGloDir;
        }
      }
    }

  }  //  iter

  return me1;
}

computeDefaultPt()

float MuonCSCSeedFromRecHits::computeDefaultPt ( ConstMuonRecHitPointer  muon ) const

private

createDefaultEndcapSeed()

bool MuonCSCSeedFromRecHits::createDefaultEndcapSeed ( ConstMuonRecHitPointer  last, TrajectorySeed &  seed  ) const

private

Definition at line 217 of file MuonCSCSeedFromRecHits.cc.

References MuonSeedFromRecHits::createSeed(), dqmdumpme::last, MuonSeedPtExtractor::pT_extract(), seed(), and MuonSeedFromRecHits::thePtExtractor.

Referenced by makeDefaultSeed().

                                                                                                            {
  //float momentum = computeDefaultPt(last);
  std::vector<double> momentum = thePtExtractor->pT_extract(last, last);
  seed = createSeed(momentum[0], momentum[1], last);
  return true;
}

makeDefaultSeed()

void MuonCSCSeedFromRecHits::makeDefaultSeed ( TrajectorySeed &  seed ) const

private

Definition at line 206 of file MuonCSCSeedFromRecHits.cc.

References bestEndcapHit(), createDefaultEndcapSeed(), hlt_dqm_clientPB-live_cfg::me1, seed(), and MuonSeedFromRecHits::theRhits.

Referenced by seed().

                                                                       {
  //Search ME1  ...
  ConstMuonRecHitPointer me1 = bestEndcapHit(theRhits);
  //  bool good=false;

  if (me1 && me1->isValid()) {
    //revert if a LogTrace or smth is necessary    good =
    createDefaultEndcapSeed(me1, seed);
  }
}

makeSeed()

bool MuonCSCSeedFromRecHits::makeSeed ( const MuonRecHitContainer &  hits1, const MuonRecHitContainer &  hits2, TrajectorySeed &  seed  ) const

private

Definition at line 78 of file MuonCSCSeedFromRecHits.cc.

References bestEndcapHit(), MuonSeedFromRecHits::createSeed(), HiEvtPlane_cfi::maxpt, HiEvtPlane_cfi::minpt, DiDispStaMuonMonitor_cfi::pt, MuonSeedPtExtractor::pT_extract(), RPCpg::pts, seed(), MuonSeedFromRecHits::thePtExtractor, and MuonSeedFromRecHits::theRhits.

Referenced by seed().

                                                                  {
  for (MuonRecHitContainer::const_iterator itr1 = hits1.begin(), end1 = hits1.end(); itr1 != end1; ++itr1) {
    CSCDetId cscId1((*itr1)->geographicalId().rawId());
    //int type1 = CSCChamberSpecs::whatChamberType(cscId1.station(), cscId1.ring());

    for (MuonRecHitContainer::const_iterator itr2 = hits2.begin(), end2 = hits2.end(); itr2 != end2; ++itr2) {
      CSCDetId cscId2((*itr2)->geographicalId().rawId());
      //int type2 = CSCChamberSpecs::whatChamberType(cscId2.station(), cscId2.ring());

      // take the first pair that comes along.  Probably want to rank them later
      std::vector<double> pts = thePtExtractor->pT_extract(*itr1, *itr2);

      double pt = pts[0];
      double sigmapt = pts[1];
      double minpt = 3.;

      // if too small, probably an error.  Keep trying.
      if (fabs(pt) > minpt) {
        double maxpt = 2000.;
        if (pt > maxpt) {
          pt = maxpt;
          sigmapt = maxpt;
        }
        if (pt < -maxpt) {
          pt = -maxpt;
          sigmapt = maxpt;
        }

        // get the position and direction from the higher-quality segment
        ConstMuonRecHitPointer bestSeg = bestEndcapHit(theRhits);
        seed = createSeed(pt, sigmapt, bestSeg);

        //std::cout << "FITTED TIMESPT " << pt << " dphi " << dphi << " eta " << eta << std::endl;
        return true;
      }
    }
  }

  // guess it didn't find one
  //std::cout << "NOTHING FOUND ! " << std::endl;
  return false;
}

makeSeed2()

bool MuonCSCSeedFromRecHits::makeSeed2 ( const MuonRecHitContainer &  hits1, const MuonRecHitContainer &  hits2, TrajectorySeed &  seed  ) const

private

seed()

TrajectorySeed MuonCSCSeedFromRecHits::seed ( ) const

virtual

Definition at line 17 of file MuonCSCSeedFromRecHits.cc.

References mps_fire::end, makeDefaultSeed(), makeSeed(), mps_fire::result, segmentQuality(), relativeConstraints::station, and MuonSeedFromRecHits::theRhits.

Referenced by createDefaultEndcapSeed(), makeDefaultSeed(), makeSeed(), and MuonSeedFinder::seeds().

                                                  {
  TrajectorySeed result;
  if (theRhits.size() == 1) {
    //return createSeed(100., 100., theRhits[0]);
    makeDefaultSeed(result);
    return result;
  }
  //@@ doesn't handle overlap between ME11 and ME12 correctly
  // sort by station
  MuonRecHitContainer station1Hits, station2Hits, station3Hits, station4Hits;
  for (MuonRecHitContainer::const_iterator iter = theRhits.begin(), end = theRhits.end(); iter != end; ++iter) {
    int station = CSCDetId((*iter)->geographicalId().rawId()).station();
    if ((*iter)->isME0()) {
      station = 1;  //ME0DetId((*iter)->geographicalId().rawId()).station();
    }

    if (station == 1) {
      station1Hits.push_back(*iter);
    } else if (station == 2) {
      station2Hits.push_back(*iter);
    } else if (station == 3) {
      station3Hits.push_back(*iter);
    } else if (station == 4) {
      station4Hits.push_back(*iter);
    }
  }

  //std::cout << "Station hits " << station1Hits.size() << " "
  //                            << station2Hits.size() << " "
  //                             << station3Hits.size() << std::endl;

  // see whether station 2 or station 3 is better
  MuonRecHitContainer* betterSecondHits = &station2Hits;
  MuonRecHitContainer* notAsGoodSecondHits = &station3Hits;
  if (!station2Hits.empty() && !station3Hits.empty()) {
    // swap if station 3 has better quailty
    if (segmentQuality(station3Hits[0]) < segmentQuality(station2Hits[0])) {
      betterSecondHits = &station3Hits;
      notAsGoodSecondHits = &station2Hits;
    }
  }

  // now try to make pairs
  if (makeSeed(station1Hits, *betterSecondHits, result)) {
    return result;
  }
  if (makeSeed(station1Hits, *notAsGoodSecondHits, result)) {
    return result;
  }
  if (makeSeed(station2Hits, station3Hits, result)) {
    return result;
  }
  if (makeSeed(station1Hits, station4Hits, result)) {
    return result;
  }

  // no luck
  makeDefaultSeed(result);
  return result;
}

segmentQuality()

int MuonCSCSeedFromRecHits::segmentQuality ( ConstMuonRecHitPointer  muon ) const

private

Definition at line 124 of file MuonCSCSeedFromRecHits.cc.

References relativeConstraints::chamber, SiPixelRawToDigiRegional_cfi::deltaPhi, nhits, quality, and relativeConstraints::ring.

Referenced by bestEndcapHit(), and seed().

                                                                               {
  int Nchi2 = 0;
  int quality = 0;
  int nhits = segment->recHits().size();
  if (segment->chi2() / (nhits * 2. - 4.) > 3.)
    Nchi2 = 1;
  if (segment->chi2() / (nhits * 2. - 4.) > 9.)
    Nchi2 = 2;

  if (nhits > 4)
    quality = 1 + Nchi2;
  if (nhits == 4)
    quality = 3 + Nchi2;
  if (nhits == 3)
    quality = 5 + Nchi2;

  float dPhiGloDir = fabs(deltaPhi(segment->globalPosition().barePhi(), segment->globalDirection().barePhi()));

  if (dPhiGloDir > .2)
    ++quality;
  // add a penalty for being ME1A if the chamber is ganged
  if (segment->isCSC() and CSCDetId(segment->geographicalId()).ring() == 4) {
    const auto chamber = dynamic_cast<const CSCChamber*>(segment->det());
    if (chamber->specs()->gangedStrips())
      ++quality;
  }

  return quality;
}