#include <MuonCSCSeedFromRecHits.h>

Inheritance diagram for MuonCSCSeedFromRecHits:

Public Member Functions
ConstMuonRecHitPointer	bestEndcapHit (const MuonRecHitContainer &endcapHits) const

	MuonCSCSeedFromRecHits ()

virtual TrajectorySeed	seed () const

virtual	~MuonCSCSeedFromRecHits ()

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

Definition at line 16 of file MuonCSCSeedFromRecHits.cc.

 : MuonSeedFromRecHits()
 {
 }

virtual MuonCSCSeedFromRecHits::~MuonCSCSeedFromRecHits ( )

inlinevirtual

Definition at line 11 of file MuonCSCSeedFromRecHits.h.

References analyze(), bestEndcapHit(), computeDefaultPt(), createDefaultEndcapSeed(), cosmicBasicClusters_cfi::endcapHits, plotBeamSpotDB::last, makeDefaultSeed(), makeSeed(), makeSeed2(), seed(), and segmentQuality().

11 {}

Member Function Documentation

void MuonCSCSeedFromRecHits::analyze ( ) const

private

Definition at line 266 of file MuonCSCSeedFromRecHits.cc.

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

Referenced by ~MuonCSCSeedFromRecHits().

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

MuonCSCSeedFromRecHits::ConstMuonRecHitPointer MuonCSCSeedFromRecHits::bestEndcapHit ( const MuonRecHitContainer & endcapHits ) const

Definition at line 185 of file MuonCSCSeedFromRecHits.cc.

References MuonCkfTrajectoryBuilder_cfi::deltaPhi, M_PI, btvTracks_cfi::quality, and segmentQuality().

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

 {
   MuonRecHitPointer me1=0, meit=0;
   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;
 }

float MuonCSCSeedFromRecHits::computeDefaultPt ( ConstMuonRecHitPointer muon ) const

private

Referenced by ~MuonCSCSeedFromRecHits().

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

private

Definition at line 256 of file MuonCSCSeedFromRecHits.cc.

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

Referenced by makeDefaultSeed(), and ~MuonCSCSeedFromRecHits().

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

void MuonCSCSeedFromRecHits::makeDefaultSeed ( TrajectorySeed & seed ) const

private

Definition at line 238 of file MuonCSCSeedFromRecHits.cc.

References bestEndcapHit(), createDefaultEndcapSeed(), and MuonSeedFromRecHits::theRhits.

Referenced by seed(), and ~MuonCSCSeedFromRecHits().

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

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

private

Definition at line 103 of file MuonCSCSeedFromRecHits.cc.

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

Referenced by seed(), and ~MuonCSCSeedFromRecHits().

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

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

private

Referenced by ~MuonCSCSeedFromRecHits().

TrajectorySeed MuonCSCSeedFromRecHits::seed ( ) const

virtual

Definition at line 22 of file MuonCSCSeedFromRecHits.cc.

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

Referenced by MuonSeedFinder::seeds(), and ~MuonCSCSeedFromRecHits().

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

int MuonCSCSeedFromRecHits::segmentQuality ( ConstMuonRecHitPointer muon ) const

private

Definition at line 157 of file MuonCSCSeedFromRecHits.cc.

References relativeConstraints::chamber, MuonCkfTrajectoryBuilder_cfi::deltaPhi, nhits, btvTracks_cfi::quality, and relativeConstraints::ring.

Referenced by bestEndcapHit(), seed(), and ~MuonCSCSeedFromRecHits().

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

Public Member Functions

Private Member Functions

Additional Inherited Members

Detailed Description

Constructor & Destructor Documentation

Member Function Documentation