MuonOverlapSeedFromRecHits Class Reference

#include <MuonOverlapSeedFromRecHits.h>

Inheritance diagram for MuonOverlapSeedFromRecHits:

Public Member Functions
bool	makeSeed (MuonTransientTrackingRecHit::ConstMuonRecHitPointer barrelHit, MuonTransientTrackingRecHit::ConstMuonRecHitPointer endcapHit, MuonTransientTrackingRecHit::ConstMuonRecHitPointer bestSegment, TrajectorySeed &result) const
	MuonOverlapSeedFromRecHits ()
std::vector< TrajectorySeed >	seeds () const
	~MuonOverlapSeedFromRecHits () 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
ConstMuonRecHitPointer	bestHit (const MuonRecHitContainer &barrelHits, const MuonRecHitContainer &endcapHits) 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 MuonOverlapSeedFromRecHits.h.

Constructor & Destructor Documentation

MuonOverlapSeedFromRecHits()

MuonOverlapSeedFromRecHits::MuonOverlapSeedFromRecHits

(

)

Definition at line 13 of file MuonOverlapSeedFromRecHits.cc.

: MuonSeedFromRecHits() {}

~MuonOverlapSeedFromRecHits()

MuonOverlapSeedFromRecHits::~MuonOverlapSeedFromRecHits ( )

inlineoverride

Definition at line 9 of file MuonOverlapSeedFromRecHits.h.

{}

Member Function Documentation

bestHit()

MuonTransientTrackingRecHit::ConstMuonRecHitPointer MuonOverlapSeedFromRecHits::bestHit ( const MuonRecHitContainer &  barrelHits, const MuonRecHitContainer &  endcapHits  ) const

private

Definition at line 83 of file MuonOverlapSeedFromRecHits.cc.

References HLT_2022v11_cff::barrelHits, MuonDTSeedFromRecHits::bestBarrelHit(), MuonCSCSeedFromRecHits::bestEndcapHit(), HLT_2022v11_cff::endcapHits, and mps_fire::result.

Referenced by seeds().

                                                                            {
  MuonDTSeedFromRecHits dtSeeder;
  MuonCSCSeedFromRecHits cscSeeder;

  ConstMuonRecHitPointer result;
  if (barrelHits.size() > endcapHits.size()) {
    result = dtSeeder.bestBarrelHit(barrelHits);
    if (result->dimension() == 2)
      result = cscSeeder.bestEndcapHit(endcapHits);
  } else {
    result = cscSeeder.bestEndcapHit(endcapHits);
  }
  return result;
}

makeSeed()

bool MuonOverlapSeedFromRecHits::makeSeed	(	MuonTransientTrackingRecHit::ConstMuonRecHitPointer	barrelHit,
		MuonTransientTrackingRecHit::ConstMuonRecHitPointer	endcapHit,
		MuonTransientTrackingRecHit::ConstMuonRecHitPointer	bestSegment,
		TrajectorySeed &	result
	)		const

Definition at line 54 of file MuonOverlapSeedFromRecHits.cc.

References MuonSeedFromRecHits::createSeed(), HiEvtPlane_cfi::maxpt, HiEvtPlane_cfi::minpt, DiDispStaMuonMonitor_cfi::pt, MuonSeedPtExtractor::pT_extract(), RPCpg::pts, mps_fire::result, and MuonSeedFromRecHits::thePtExtractor.

Referenced by seeds().

                                                                        {
  std::vector<double> pts = thePtExtractor->pT_extract(barrelHit, endcapHit);
  double minpt = 3.;
  double pt = pts[0];
  double sigmapt = pts[1];
  // if too small, probably an error.  Keep trying.
  if (pt != 0) {
    if (fabs(pt) > minpt) {
      double maxpt = 2000.;
      if (pt > maxpt) {
        pt = maxpt;
        sigmapt = maxpt;
      }
      if (pt < -maxpt) {
        pt = -maxpt;
        sigmapt = maxpt;
      }
    }

    result = createSeed(pt, sigmapt, bestSegment);
    //std::cout << "OVERLAPFITTED PT " << pt << " dphi " << dphi << " eta " << eta << std::endl;
    return true;
  }
  return false;
}

seeds()

std::vector< TrajectorySeed > MuonOverlapSeedFromRecHits::seeds

(

)

const

Definition at line 15 of file MuonOverlapSeedFromRecHits.cc.

References HLT_2022v11_cff::barrelHits, bestHit(), mps_fire::end, HLT_2022v11_cff::endcapHits, good, makeSeed(), mps_fire::result, fileCollector::seed, and MuonSeedFromRecHits::theRhits.

Referenced by MuonSeedFinder::seeds().

                                                                  {
  std::vector<TrajectorySeed> result;
  //@@ doesn't handle overlap between ME11 and ME12 correctly
  // sort by station
  MuonRecHitContainer barrelHits, endcapHits;
  for (MuonRecHitContainer::const_iterator iter = theRhits.begin(), end = theRhits.end(); iter != end; ++iter) {
    if ((*iter)->isDT()) {
      DTChamberId dtId((**iter).geographicalId().rawId());
      // try not doing seeds that start in DT station 2, if there'as a good single segment seed
      if (dtId.station() == 1 || (dtId.station() == 2 && (*iter)->dimension() == 2)) {
        barrelHits.push_back(*iter);
      }
    } else {
      endcapHits.push_back(*iter);
    }
  }

  ConstMuonRecHitPointer bestSegment = bestHit(barrelHits, endcapHits);
  for (MuonRecHitContainer::const_iterator barrelHitItr = barrelHits.begin(), lastBarrelHit = barrelHits.end();
       barrelHitItr != lastBarrelHit;
       ++barrelHitItr) {
    for (MuonRecHitContainer::const_iterator endcapHitItr = endcapHits.begin(), lastEndcapHit = endcapHits.end();
         endcapHitItr != lastEndcapHit;
         ++endcapHitItr) {
      TrajectorySeed seed;
      bool good = makeSeed(*barrelHitItr, *endcapHitItr, bestSegment, seed);
      if (good)
        result.push_back(seed);
      // try just one seed
      return result;
    }
  }

  //std::cout << "Overlap hits " << barrelHits.size() << " "
  //                             << endcapHits.size() << std::endl;

  return result;
}