#include <PFMuonAlgo.h>

Classes
class	IndexPtComparator

class	TrackMETComparator

class	TrackPtErrorSorter

Public Member Functions
void	addMissingMuons (edm::Handle< reco::MuonCollection >, reco::PFCandidateCollection *cands)

void	changeTrack (reco::PFCandidate &, const MuonTrackTypePair &)

	PFMuonAlgo ()
	constructor More...

void	postClean (reco::PFCandidateCollection *)

bool	reconstructMuon (reco::PFCandidate &, const reco::MuonRef &, bool allowLoose=false)

void	setInputsForCleaning (const reco::VertexCollection *)

void	setParameters (const edm::ParameterSet &)

std::auto_ptr < reco::PFCandidateCollection > &	transferAddedMuonCandidates ()

std::auto_ptr < reco::PFCandidateCollection > &	transferCleanedCosmicCandidates ()

std::auto_ptr < reco::PFCandidateCollection > &	transferCleanedFakeCandidates ()

std::auto_ptr < reco::PFCandidateCollection > &	transferCleanedTrackerAndGlobalCandidates ()

std::auto_ptr < reco::PFCandidateCollection > &	transferPunchThroughCleanedHadronCandidates ()

std::auto_ptr < reco::PFCandidateCollection > &	transferPunchThroughCleanedMuonCandidates ()

virtual	~PFMuonAlgo ()
	destructor More...

Static Public Member Functions
static bool	hasValidTracks (const reco::MuonRef &)

static bool	isGlobalLooseMuon (const reco::PFBlockElement &elt)

static bool	isGlobalLooseMuon (const reco::MuonRef &muonRef)

static bool	isGlobalTightMuon (const reco::PFBlockElement &elt)

static bool	isGlobalTightMuon (const reco::MuonRef &muonRef)

static bool	isIsolatedMuon (const reco::PFBlockElement &elt)

static bool	isIsolatedMuon (const reco::MuonRef &muonRef)

static bool	isLooseMuon (const reco::PFBlockElement &elt)

static bool	isLooseMuon (const reco::MuonRef &muonRef)

static bool	isMuon (const reco::PFBlockElement &elt)

static bool	isMuon (const reco::MuonRef &muonRef)

static bool	isTightMuonPOG (const reco::MuonRef &muonRef)

static bool	isTrackerLooseMuon (const reco::PFBlockElement &elt)

static bool	isTrackerLooseMuon (const reco::MuonRef &muonRef)

static bool	isTrackerTightMuon (const reco::PFBlockElement &elt)

static bool	isTrackerTightMuon (const reco::MuonRef &muonRef)

static void	printMuonProperties (const reco::MuonRef &muonRef)

Private Types
typedef reco::Muon::MuonTrackType	MuonTrackType

typedef reco::Muon::MuonTrackTypePair	MuonTrackTypePair

Private Member Functions
bool	cleanMismeasured (reco::PFCandidate &, unsigned int)

bool	cleanPunchThroughAndFakes (reco::PFCandidate &, reco::PFCandidateCollection *, unsigned int)

void	estimateEventQuantities (const reco::PFCandidateCollection *)

std::pair< double, double >	getMinMaxMET2 (const reco::PFCandidate &)

MuonTrackTypePair	getTrackWithSmallestError (const std::vector< MuonTrackTypePair > &)

std::vector < reco::Muon::MuonTrackTypePair >	goodMuonTracks (const reco::MuonRef &muon, bool includeSA=false)

std::vector < reco::Muon::MuonTrackTypePair >	muonTracks (const reco::MuonRef &muon, bool includeSA=false, double dpt=1e+9)

void	removeDeadCandidates (reco::PFCandidateCollection *, const std::vector< unsigned int > &)

std::vector< MuonTrackTypePair >	tracksPointingAtMET (const std::vector< MuonTrackTypePair > &)

std::vector< MuonTrackTypePair >	tracksWithBetterMET (const std::vector< MuonTrackTypePair > &, const reco::PFCandidate &)

Private Attributes
double	cosmicRejDistance_

double	dzPV_

double	errorCompScale_

double	eventFactorCosmics_

double	eventFractionCleaning_

double	eventFractionRejection_

std::vector< unsigned int >	maskedIndices_

double	maxDPtOPt_

double	metFactorCleaning_

double	metFactorFake_

double	metFactorHighEta_

double	metFactorRejection_

double	metSigForCleaning_

double	metSigForRejection_

double	METX_

double	METY_

int	minPixelHits_

double	minPostCleaningPt_

double	minPunchThroughEnergy_

double	minPunchThroughMomentum_

int	minTrackerHits_

std::auto_ptr < reco::PFCandidateCollection >	pfAddedMuonCandidates_
	the collection of added muon candidates More...

std::auto_ptr < reco::PFCandidateCollection >	pfCleanedTrackerAndGlobalMuonCandidates_
	the collection of tracker/global cleaned muon candidates More...

std::auto_ptr < reco::PFCandidateCollection >	pfCosmicsMuonCleanedCandidates_
	the collection of cosmics cleaned muon candidates More...

std::auto_ptr < reco::PFCandidateCollection >	pfFakeMuonCleanedCandidates_
	the collection of fake cleaned muon candidates More...

std::auto_ptr < reco::PFCandidateCollection >	pfPunchThroughHadronCleanedCandidates_
	the collection of punch-through cleaned neutral hadron candidates More...

std::auto_ptr < reco::PFCandidateCollection >	pfPunchThroughMuonCleanedCandidates_
	the collection of punch-through cleaned muon candidates More...

bool	postCleaning_

double	ptFactorHighEta_

double	punchThroughFactor_

double	punchThroughMETFactor_

double	sumet_

double	sumetPU_

reco::TrackBase::TrackQuality	trackQuality_

const reco::VertexCollection *	vertices_

Detailed Description

Definition at line 12 of file PFMuonAlgo.h.

Member Typedef Documentation

typedef reco::Muon::MuonTrackType PFMuonAlgo::MuonTrackType

private

Definition at line 15 of file PFMuonAlgo.h.

typedef reco::Muon::MuonTrackTypePair PFMuonAlgo::MuonTrackTypePair

private

Definition at line 14 of file PFMuonAlgo.h.

Constructor & Destructor Documentation

PFMuonAlgo::PFMuonAlgo ( )

constructor

Definition at line 17 of file PFMuonAlgo.cc.

                        {
   pfCosmicsMuonCleanedCandidates_ = std::auto_ptr<reco::PFCandidateCollection>(new reco::PFCandidateCollection);
   pfCleanedTrackerAndGlobalMuonCandidates_= std::auto_ptr<reco::PFCandidateCollection>(new reco::PFCandidateCollection);
   pfFakeMuonCleanedCandidates_= std::auto_ptr<reco::PFCandidateCollection>(new reco::PFCandidateCollection);
   pfPunchThroughMuonCleanedCandidates_= std::auto_ptr<reco::PFCandidateCollection>(new reco::PFCandidateCollection);
   pfPunchThroughHadronCleanedCandidates_= std::auto_ptr<reco::PFCandidateCollection>(new reco::PFCandidateCollection);
   pfAddedMuonCandidates_= std::auto_ptr<reco::PFCandidateCollection>(new reco::PFCandidateCollection);
   
 }

virtual PFMuonAlgo::~PFMuonAlgo ( )

inlinevirtual

destructor

Definition at line 27 of file PFMuonAlgo.h.

27 {;}

Member Function Documentation

void PFMuonAlgo::addMissingMuons	(	edm::Handle< reco::MuonCollection >	muons,
		reco::PFCandidateCollection *	cands
	)

Definition at line 938 of file PFMuonAlgo.cc.

References dPhi(), PFRecoTauDiscriminationAgainstElectronDeadECAL_cfi::dR, i, edm::Ref< C, T, F >::isNonnull(), reco::PFCandidate::mu, reco::PFCandidate::muonRef(), mathSSE::sqrt(), reco::PFCandidate::trackRef(), and testEve_cfg::tracks.

Referenced by PFAlgo::reconstructParticles().

                                                                                                       {
   if(!postCleaning_)
     return;
 
 
   if( pfAddedMuonCandidates_.get() )
      pfAddedMuonCandidates_->clear();
   else 
      pfAddedMuonCandidates_.reset( new reco::PFCandidateCollection );
 
 
 
   for ( unsigned imu = 0; imu < muons->size(); ++imu ) {
     reco::MuonRef muonRef( muons, imu );
     bool used = false;
     bool hadron = false;
     for(unsigned i=0; i<cands->size(); i++) {
       const PFCandidate& pfc = cands->at(i);
       if ( !pfc.trackRef().isNonnull() ) continue;
       if ( pfc.trackRef().isNonnull() && pfc.trackRef() == muonRef->track() )
     hadron = true;
       if ( !pfc.muonRef().isNonnull() ) continue;
 
       if ( pfc.muonRef()->innerTrack() == muonRef->innerTrack())
     used = true;
       else {
     // Check if the stand-alone muon is not a spurious copy of an existing muon 
     // (Protection needed for HLT)
     if ( pfc.muonRef()->isStandAloneMuon() && muonRef->isStandAloneMuon() ) { 
       double dEta = pfc.muonRef()->standAloneMuon()->eta() - muonRef->standAloneMuon()->eta();
       double dPhi = pfc.muonRef()->standAloneMuon()->phi() - muonRef->standAloneMuon()->phi();
       double dR = sqrt(dEta*dEta + dPhi*dPhi);
       if ( dR < 0.005 ) { 
         used = true;
       }
     }
       }
   
     if ( used ) break; 
     }
 
     if ( used ||hadron||(!muonRef.isNonnull()) ) continue;
 
 
     TrackMETComparator comparator(METX_,METY_);
     //Low pt dont need to be cleaned
   
     std::vector<reco::Muon::MuonTrackTypePair> tracks  = goodMuonTracks(muonRef,true);
     //If there is at least 1 track choice  try to change the track 
     if(tracks.size()>0) {
 
     //Find tracks that change dramatically MET or Pt
     std::vector<reco::Muon::MuonTrackTypePair> tracksThatChangeMET = tracksPointingAtMET(tracks);
     //From those tracks get the one with smallest MET 
     if (tracksThatChangeMET.size()>0) {
       reco::Muon::MuonTrackTypePair bestTrackType = *std::min_element(tracksThatChangeMET.begin(),tracksThatChangeMET.end(),comparator);
 
       //Make sure it is not cosmic
       if((vertices_->size()==0) ||bestTrackType.first->dz(vertices_->at(0).position())<cosmicRejDistance_){
     
     //make a pfcandidate
     int charge = bestTrackType.first->charge()>0 ? 1 : -1;
     math::XYZTLorentzVector momentum(bestTrackType.first->px(),
                      bestTrackType.first->py(),
                      bestTrackType.first->pz(),
                        sqrt(bestTrackType.first->p()*bestTrackType.first->p()+0.1057*0.1057));
       
     cands->push_back( PFCandidate( charge, 
                       momentum,
                       reco::PFCandidate::mu ) );
 
     changeTrack(cands->back(),bestTrackType);
 
     if (muonRef->track().isNonnull() ) 
       cands->back().setTrackRef( muonRef->track() );
 
     cands->back().setMuonRef(muonRef);
 
 
     pfAddedMuonCandidates_->push_back(cands->back());
 
       }
     }
     }
   }
 }

void PFMuonAlgo::changeTrack	(	reco::PFCandidate &	candidate,
		const MuonTrackTypePair &	track
	)

Definition at line 770 of file PFMuonAlgo.cc.

References reco::Muon::CombinedTrack, reco::Muon::DYT, relval_parameters_module::energy, reco::Muon::InnerTrack, reco::PFCandidate::mu, reco::Muon::Picky, dt_dqm_sourceclient_common_cff::reco, reco::LeafCandidate::setCharge(), reco::PFCandidate::setMuonTrackType(), reco::LeafCandidate::setP4(), reco::PFCandidate::setParticleType(), reco::PFCandidate::setVertexSource(), mathSSE::sqrt(), and reco::Muon::TPFMS.

                                                                                        {
   using namespace reco;
     reco::TrackRef      bestTrack = track.first;
     MuonTrackType trackType = track.second;
     //OK Now redefine the canddiate with that track
     double px = bestTrack->px();
     double py = bestTrack->py();
     double pz = bestTrack->pz();
     double energy = sqrt(bestTrack->p()*bestTrack->p() + 0.1057*0.1057);
 
     candidate.setCharge(bestTrack->charge()>0 ? 1 : -1);
     candidate.setP4(math::XYZTLorentzVector(px,py,pz,energy));
     candidate.setParticleType(reco::PFCandidate::mu);
     //    candidate.setTrackRef( bestTrack );  
     candidate.setMuonTrackType(trackType);
     if(trackType == reco::Muon::InnerTrack)
       candidate.setVertexSource( PFCandidate::kTrkMuonVertex );
     else if(trackType == reco::Muon::CombinedTrack)
       candidate.setVertexSource( PFCandidate::kComMuonVertex );
     else if(trackType == reco::Muon::TPFMS)
       candidate.setVertexSource( PFCandidate::kTPFMSMuonVertex );
     else if(trackType == reco::Muon::Picky)
       candidate.setVertexSource( PFCandidate::kPickyMuonVertex );
     else if(trackType == reco::Muon::DYT)
       candidate.setVertexSource( PFCandidate::kDYTMuonVertex );
   }

bool PFMuonAlgo::cleanMismeasured	(	reco::PFCandidate &	pfc,
		unsigned int	i
	)

private

Definition at line 1046 of file PFMuonAlgo.cc.

References reco::PFCandidate::muonRef(), reco::LeafCandidate::pt(), reco::LeafCandidate::px(), reco::LeafCandidate::py(), dt_dqm_sourceclient_common_cff::reco, mathSSE::sqrt(), and testEve_cfg::tracks.

                                                                      {
   using namespace std;
   using namespace reco;
   bool cleaned=false;
 
   //First define the MET without this guy
   double METNOX = METX_ - pfc.px();
   double METNOY = METY_ - pfc.py();
   double SUMETNO = sumet_  -pfc.pt(); 
   
   TrackMETComparator comparator(METNOX,METNOY);
   //Low pt dont need to be cleaned
   if (pfc.pt()<minPostCleaningPt_)
     return false;
   std::vector<reco::Muon::MuonTrackTypePair> tracks  = goodMuonTracks(pfc.muonRef(),false);
   
 
 
   //If there is more than 1 track choice  try to change the track 
   if(tracks.size()>1) {
     //Find tracks that change dramatically MET or Pt
     std::vector<reco::Muon::MuonTrackTypePair> tracksThatChangeMET = tracksWithBetterMET(tracks,pfc);
     //From those tracks get the one with smallest MET 
     if (tracksThatChangeMET.size()>0) {
       reco::Muon::MuonTrackTypePair bestTrackType = *std::min_element(tracksThatChangeMET.begin(),tracksThatChangeMET.end(),comparator);
       changeTrack(pfc,bestTrackType);
 
       pfCleanedTrackerAndGlobalMuonCandidates_->push_back(pfc);
       //update eventquantities
       METX_ = METNOX+pfc.px();
       METY_ = METNOY+pfc.py();
       sumet_=SUMETNO+pfc.pt();
 
     }      
   }
 
   //Now attempt to kill it 
   if (!(pfc.muonRef()->isGlobalMuon() && pfc.muonRef()->isTrackerMuon())) {
     //define MET significance and SUM ET
     double MET2 = METX_*METX_+METY_*METY_;
     double newMET2 = METNOX*METNOX+METNOY*METNOY;
     double METSig = sqrt(MET2)/sqrt(sumet_-sumetPU_);
     if( METSig>metSigForRejection_)
       if((newMET2 < MET2/metFactorRejection_) &&
      ((SUMETNO-sumetPU_)/(sumet_-sumetPU_)<eventFractionRejection_)) {
        pfFakeMuonCleanedCandidates_->push_back(pfc);
        maskedIndices_.push_back(i);
        METX_ = METNOX;
        METY_ = METNOY;
        sumet_=SUMETNO;
        cleaned=true;
     }
     
   }
     return cleaned;
 
 }

bool PFMuonAlgo::cleanPunchThroughAndFakes	(	reco::PFCandidate &	pfc,
		reco::PFCandidateCollection *	cands,
		unsigned int	imu
	)

private

Definition at line 1157 of file PFMuonAlgo.cc.

References reco::PFCandidate::elementsInBlocks(), reco::LeafCandidate::eta(), reco::PFCandidate::h, reco::PFCandidate::h0, i, edm::Ref< C, T, F >::isNonnull(), edm::Ref< C, T, F >::key(), bookConverter::max, reco::PFCandidate::muonRef(), reco::LeafCandidate::p(), reco::PFCandidate::particleId(), reco::LeafCandidate::pt(), reco::PFCandidate::rawEcalEnergy(), reco::PFCandidate::rawHcalEnergy(), dt_dqm_sourceclient_common_cff::reco, and edm::second().

                                                                                                                 {
   using namespace reco;
 
   bool cleaned=false;
 
   if (pfc.pt()<minPostCleaningPt_)
     return false;
 
 
   double METXNO = METX_-pfc.pt();
   double METYNO = METY_-pfc.pt();
   double MET2NO = METXNO*METXNO+METYNO*METYNO;
   double MET2   = METX_*METX_+METY_*METY_;
   bool fake1=false;
 
   std::pair<double,double> met2 = getMinMaxMET2(pfc);
 
   //Check for Fakes at high pseudorapidity
   if(pfc.muonRef()->standAloneMuon().isNonnull()) 
     fake1 =fabs ( pfc.eta() ) > 2.15 && 
       met2.first<met2.second/2 &&
       MET2NO < MET2/metFactorHighEta_ && 
       pfc.muonRef()->standAloneMuon()->pt() < pfc.pt()/ptFactorHighEta_;
 
   double factor = std::max(2.,2000./(sumet_-pfc.pt()-sumetPU_));
   bool fake2 = ( pfc.pt()/(sumet_-sumetPU_) < 0.25 && MET2NO < MET2/metFactorFake_ && met2.first<met2.second/factor );
 
   bool punchthrough =pfc.p() > minPunchThroughMomentum_ &&  
     pfc.rawHcalEnergy() > minPunchThroughEnergy_ && 
     pfc.rawEcalEnergy()+pfc.rawHcalEnergy() > pfc.p()/punchThroughFactor_ &&
     !isIsolatedMuon(pfc.muonRef()) && MET2NO < MET2/punchThroughMETFactor_;
 
 
   if(fake1 || fake2||punchthrough) {
     // Find the block of the muon
     const PFCandidate::ElementsInBlocks& eleInBlocks = pfc.elementsInBlocks();
     if ( eleInBlocks.size() ) { 
       PFBlockRef blockRefMuon = eleInBlocks[0].first;
       unsigned indexMuon = eleInBlocks[0].second;
       if (eleInBlocks.size()>1)
     indexMuon = eleInBlocks[1].second;
       
       // Check if the muon gave rise to a neutral hadron
       double iHad = 1E9;
       bool hadron = false;
       for ( unsigned i = imu+1; i < cands->size(); ++i ) { 
     const PFCandidate& pfcn = cands->at(i);
         const PFCandidate::ElementsInBlocks& ele = pfcn.elementsInBlocks();
         if ( !ele.size() ) { 
           continue;
         }
         PFBlockRef blockRefHadron = ele[0].first;
         unsigned indexHadron = ele[0].second;
         // We are out of the block -> exit the loop
         if ( blockRefHadron.key() != blockRefMuon.key() ) break;
         // Check that this particle is a neutral hadron
         if ( indexHadron == indexMuon && 
          pfcn.particleId() == reco::PFCandidate::h0 ) {
           iHad = i;
           hadron = true;
         }
         if ( hadron ) break;
       }
       
       if ( hadron ) { 
 
     double rescaleFactor = cands->at(iHad).p()/cands->at(imu).p();
         METX_ -=  cands->at(imu).px() + cands->at(iHad).px();
         METY_ -=  cands->at(imu).py() + cands->at(iHad).py();
         sumet_ -=cands->at(imu).pt();
         cands->at(imu).rescaleMomentum(rescaleFactor);
         maskedIndices_.push_back(iHad);
         pfPunchThroughHadronCleanedCandidates_->push_back(cands->at(iHad));
         cands->at(imu).setParticleType(reco::PFCandidate::h);
         pfPunchThroughMuonCleanedCandidates_->push_back(cands->at(imu));
         METX_ +=  cands->at(imu).px();
         METY_ +=  cands->at(imu).py();    
         sumet_ += cands->at(imu).pt();
 
       } else if ( fake1 || fake2 ) {
         METX_ -=  cands->at(imu).px();
         METY_ -=  cands->at(imu).py();    
         sumet_ -= cands->at(imu).pt();
         maskedIndices_.push_back(imu);
         pfFakeMuonCleanedCandidates_->push_back(cands->at(imu));
         cleaned=true;
       }
     }
   }
   return cleaned;
 }

void PFMuonAlgo::estimateEventQuantities ( const reco::PFCandidateCollection * pfc )

private

Definition at line 808 of file PFMuonAlgo.cc.

References i.

 {
   //SUM ET 
   sumetPU_ = 0.0;
   METX_=0.;
   METY_=0.;
   sumet_=0.0;
   for(reco::PFCandidateCollection::const_iterator i = pfc->begin();i!=pfc->end();++i) {
     sumet_+=i->pt();
     METX_+=i->px();
     METY_+=i->py();
   }
 
 }

std::pair< double, double > PFMuonAlgo::getMinMaxMET2 ( const reco::PFCandidate & pfc )

private

Definition at line 1026 of file PFMuonAlgo.cc.

References i, reco::PFCandidate::muonRef(), funct::pow(), reco::LeafCandidate::px(), reco::LeafCandidate::py(), and testEve_cfg::tracks.

                                                    {
   std::vector<reco::Muon::MuonTrackTypePair> tracks  = goodMuonTracks((pfc.muonRef()),true);
 
   double METXNO = METX_-pfc.px();
   double METYNO = METY_-pfc.py();
   std::vector<double> met2;
   for (unsigned int i=0;i<tracks.size();++i) {
     met2.push_back(pow(METXNO+tracks.at(i).first->px(),2)+pow(METYNO+tracks.at(i).first->py(),2));
   }
 
   //PROTECT for cases of only one track. If there is only one track it will crash .
   //Has never happened but could likely happen!
 
   if(tracks.size()>1)
     return std::make_pair(*std::min_element(met2.begin(),met2.end()),*std::max_element(met2.begin(),met2.end()));
   else
     return std::make_pair(0,0);
 }

reco::Muon::MuonTrackTypePair PFMuonAlgo::getTrackWithSmallestError ( const std::vector< MuonTrackTypePair > & tracks )

private

Definition at line 799 of file PFMuonAlgo.cc.

References MCScenario_CRAFT1_22X::sorter().

                                                                                         {
     TrackPtErrorSorter sorter;
     return *std::min_element(tracks.begin(),tracks.end(),sorter);
 }

std::vector< reco::Muon::MuonTrackTypePair > PFMuonAlgo::goodMuonTracks	(	const reco::MuonRef &	muon,
		bool	includeSA = `false`
	)

private

Definition at line 633 of file PFMuonAlgo.cc.

                                                                                                       {
   return muonTracks(muon,includeSA,maxDPtOPt_);
 }

static bool PFMuonAlgo::hasValidTracks ( const reco::MuonRef & )

static

bool PFMuonAlgo::isGlobalLooseMuon ( const reco::PFBlockElement & elt )

static

Definition at line 198 of file PFMuonAlgo.cc.

References reco::PFBlockElementTrack::muonRef().

                                                              {
 
   const reco::PFBlockElementTrack* eltTrack 
     = dynamic_cast<const reco::PFBlockElementTrack*>(&elt);
 
   assert ( eltTrack );
   reco::MuonRef muonRef = eltTrack->muonRef();
 
   return isGlobalLooseMuon(muonRef);
 
 }

bool PFMuonAlgo::isGlobalLooseMuon ( const reco::MuonRef & muonRef )

static

Definition at line 341 of file PFMuonAlgo.cc.

References muon::isGoodMuon(), edm::Ref< C, T, F >::isNonnull(), min(), query::result, and muon::TMLastStationAngTight.

                                                           {
 
   if ( !muonRef.isNonnull() ) return false;
   if ( !muonRef->isGlobalMuon() ) return false;
   if ( !muonRef->isStandAloneMuon() ) return false;
   
   reco::TrackRef standAloneMu = muonRef->standAloneMuon();
   reco::TrackRef combinedMu = muonRef->combinedMuon();
   reco::TrackRef trackerMu = muonRef->track();
  
   unsigned nMuonHits =
     standAloneMu->hitPattern().numberOfValidMuonDTHits() +
     2*standAloneMu->hitPattern().numberOfValidMuonCSCHits();
     
   bool quality = false;
   
   if ( muonRef->isTrackerMuon() ){
 
     bool result = combinedMu->normalizedChi2() < 100.;
     
     bool laststation =
       muon::isGoodMuon(*muonRef,muon::TMLastStationAngTight);
         
     int nMatches = muonRef->numberOfMatches();    
     
     quality = laststation && nMuonHits > 12 && nMatches > 1;
 
     return result && quality;
     
   }
   else{
 
     // Check the quality of the stand-alone muon : 
     // good chi**2 and large number of hits and good pt error
     if (  nMuonHits <=15  ||
       standAloneMu->normalizedChi2() > 10. || 
       standAloneMu->ptError()/standAloneMu->pt() > 0.20 ) {
       quality = false;
     }
    else { 
       // If the stand-alone muon is good, check the global muon
       if ( combinedMu->normalizedChi2() > standAloneMu->normalizedChi2() ) {
     // If the combined muon is worse than the stand-alone, it 
     // means that either the corresponding tracker track was not 
     // reconstructed, or that the sta muon comes from a late 
     // pion decay (hence with a momentum smaller than the track)
     // Take the stand-alone muon only if its momentum is larger
     // than that of the track
 
     // Note that here we even take the standAlone if it has a smaller pT, in contrast to GlobalTight
     if(standAloneMu->pt() > trackerMu->pt() || combinedMu->normalizedChi2()<5.) quality =  true;
       }
       else { 
     // If the combined muon is better (and good enough), take the 
     // global muon
     if(combinedMu->ptError()/combinedMu->pt() < std::min(0.20,standAloneMu->ptError()/standAloneMu->pt())) 
       quality = true;
     
       }
    }         
   }
   
 
   return quality;
 
 }

bool PFMuonAlgo::isGlobalTightMuon ( const reco::PFBlockElement & elt )

static

Definition at line 185 of file PFMuonAlgo.cc.

References reco::PFBlockElementTrack::muonRef().

                                                              {
 
   const reco::PFBlockElementTrack* eltTrack 
     = dynamic_cast<const reco::PFBlockElementTrack*>(&elt);
 
   assert ( eltTrack );
   reco::MuonRef muonRef = eltTrack->muonRef();
   
   return isGlobalTightMuon(muonRef);
 
 }

bool PFMuonAlgo::isGlobalTightMuon ( const reco::MuonRef & muonRef )

static

Definition at line 250 of file PFMuonAlgo.cc.

References muon::GlobalMuonPromptTight, muon::isGoodMuon(), edm::Ref< C, T, F >::isNonnull(), min(), query::result, and muon::TM2DCompatibilityTight.

                                                           {
 
  if ( !muonRef.isNonnull() ) return false;
 
  if ( !muonRef->isGlobalMuon() ) return false;
  if ( !muonRef->isStandAloneMuon() ) return false;
  
  
  if ( muonRef->isTrackerMuon() ) { 
    
    bool result = muon::isGoodMuon(*muonRef,muon::GlobalMuonPromptTight);
    
    bool isTM2DCompatibilityTight =  muon::isGoodMuon(*muonRef,muon::TM2DCompatibilityTight);   
    int nMatches = muonRef->numberOfMatches();
    bool quality = nMatches > 2 || isTM2DCompatibilityTight;
    
    return result && quality;
    
  } else {
  
    reco::TrackRef standAloneMu = muonRef->standAloneMuon();
    
     // No tracker muon -> Request a perfect stand-alone muon, or an even better global muon
     bool result = false;
       
     // Check the quality of the stand-alone muon : 
     // good chi**2 and large number of hits and good pt error
     if ( ( standAloneMu->hitPattern().numberOfValidMuonDTHits() < 22 &&
        standAloneMu->hitPattern().numberOfValidMuonCSCHits() < 15 ) ||
      standAloneMu->normalizedChi2() > 10. || 
      standAloneMu->ptError()/standAloneMu->pt() > 0.20 ) {
       result = false;
     } else { 
       
       reco::TrackRef combinedMu = muonRef->combinedMuon();
       reco::TrackRef trackerMu = muonRef->track();
             
       // If the stand-alone muon is good, check the global muon
       if ( combinedMu->normalizedChi2() > standAloneMu->normalizedChi2() ) {
     // If the combined muon is worse than the stand-alone, it 
     // means that either the corresponding tracker track was not 
     // reconstructed, or that the sta muon comes from a late 
     // pion decay (hence with a momentum smaller than the track)
     // Take the stand-alone muon only if its momentum is larger
     // than that of the track
     result = standAloneMu->pt() > trackerMu->pt() ;
      } else { 
     // If the combined muon is better (and good enough), take the 
     // global muon
     result = 
       combinedMu->ptError()/combinedMu->pt() < 
       std::min(0.20,standAloneMu->ptError()/standAloneMu->pt());
       }
     }      
 
     return result;    
   }
 
   return false;
 
 }

bool PFMuonAlgo::isIsolatedMuon ( const reco::PFBlockElement & elt )

static

Definition at line 224 of file PFMuonAlgo.cc.

References reco::PFBlockElementTrack::muonRef().

Referenced by GreedyMuonPFCandidateFilter::filter(), and PFAlgo::processBlock().

                                                           {
 
   const reco::PFBlockElementTrack* eltTrack 
     = dynamic_cast<const reco::PFBlockElementTrack*>(&elt);
 
   assert ( eltTrack );
   reco::MuonRef muonRef = eltTrack->muonRef();
 
   return isIsolatedMuon(muonRef);
 
 }

bool PFMuonAlgo::isIsolatedMuon ( const reco::MuonRef & muonRef )

static

Definition at line 432 of file PFMuonAlgo.cc.

References edm::Ref< C, T, F >::isNonnull().

                                                       {
 
 
   if ( !muonRef.isNonnull() ) return false;
   if ( !muonRef->isIsolationValid() ) return false;
   
   // Isolated Muons which are missed by standard cuts are nearly always global+tracker
   if ( !muonRef->isGlobalMuon() ) return false;
 
   // If it's not a tracker muon, only take it if there are valid muon hits
 
   reco::TrackRef standAloneMu = muonRef->standAloneMuon();
 
   if ( !muonRef->isTrackerMuon() ){
     if(standAloneMu->hitPattern().numberOfValidMuonDTHits() == 0 &&
        standAloneMu->hitPattern().numberOfValidMuonCSCHits() ==0) return false;
   }
   
   // for isolation, take the smallest pt available to reject fakes
 
   reco::TrackRef combinedMu = muonRef->combinedMuon();
   double smallestMuPt = combinedMu->pt();
   
   if(standAloneMu->pt()<smallestMuPt) smallestMuPt = standAloneMu->pt();
   
   if(muonRef->isTrackerMuon())
     {
       reco::TrackRef trackerMu = muonRef->track();
       if(trackerMu->pt() < smallestMuPt) smallestMuPt= trackerMu->pt();
     }
      
   double sumPtR03 = muonRef->isolationR03().sumPt;
   double emEtR03 = muonRef->isolationR03().emEt;
   double hadEtR03 = muonRef->isolationR03().hadEt;
   
   double relIso = (sumPtR03 + emEtR03 + hadEtR03)/smallestMuPt;
 
   if(relIso<0.1) return true;
   else return false;
 }

bool PFMuonAlgo::isLooseMuon ( const reco::PFBlockElement & elt )

static

Definition at line 166 of file PFMuonAlgo.cc.

References muon::isLooseMuon(), and reco::PFBlockElementTrack::muonRef().

Referenced by GeneralTracksImporter::importToBlock(), and PFAlgo::processBlock().

                                                        {
 
   const reco::PFBlockElementTrack* eltTrack 
     = dynamic_cast<const reco::PFBlockElementTrack*>(&elt);
 
   assert ( eltTrack );
 
   reco::MuonRef muonRef = eltTrack->muonRef();
 
   return isLooseMuon(muonRef);
 
 }

bool PFMuonAlgo::isLooseMuon ( const reco::MuonRef & muonRef )

static

Definition at line 243 of file PFMuonAlgo.cc.

                                                   {
 
   return isGlobalLooseMuon(muonRef) || isTrackerLooseMuon(muonRef);
 
 }

bool PFMuonAlgo::isMuon ( const reco::PFBlockElement & elt )

static

Definition at line 153 of file PFMuonAlgo.cc.

References reco::isMuon(), and reco::PFBlockElementTrack::muonRef().

Referenced by GeneralTracksImporter::importToBlock(), PFEGammaAlgo::isAMuon(), PFAlgo::processBlock(), and PFElectronAlgo::SetLinks().

                                                   {
 
   const reco::PFBlockElementTrack* eltTrack 
     = dynamic_cast<const reco::PFBlockElementTrack*>(&elt);
 
   assert ( eltTrack );
   reco::MuonRef muonRef = eltTrack->muonRef();
 
   return isMuon(muonRef);
 
 }

bool PFMuonAlgo::isMuon ( const reco::MuonRef & muonRef )

static

Definition at line 237 of file PFMuonAlgo.cc.

                                              {
 
   return isGlobalTightMuon(muonRef) || isTrackerTightMuon(muonRef) || isIsolatedMuon(muonRef);
 }

bool PFMuonAlgo::isTightMuonPOG ( const reco::MuonRef & muonRef )

static

Definition at line 474 of file PFMuonAlgo.cc.

References muon::GlobalMuonPromptTight, and muon::isGoodMuon().

                                                      {
 
   if(!muon::isGoodMuon(*muonRef,muon::GlobalMuonPromptTight)) return false;
 
   if(!muonRef->isTrackerMuon()) return false;
   
   if(muonRef->numberOfMatches()<2) return false;
   
   //const reco::TrackRef& combinedMuon = muonRef->combinedMuon();    
   const reco::TrackRef& combinedMuon = muonRef->globalTrack();    
   
   if(combinedMuon->hitPattern().numberOfValidTrackerHits()<11) return false;
   
   if(combinedMuon->hitPattern().numberOfValidPixelHits()==0) return false;
   
   if(combinedMuon->hitPattern().numberOfValidMuonHits()==0) return false;  
 
   return true;
 
 }

static bool PFMuonAlgo::isTrackerLooseMuon ( const reco::PFBlockElement & elt )

static

bool PFMuonAlgo::isTrackerLooseMuon ( const reco::MuonRef & muonRef )

static

Definition at line 410 of file PFMuonAlgo.cc.

References muon::AllArbitrated, muon::isGoodMuon(), edm::Ref< C, T, F >::isNonnull(), and muon::TMLastStationAngTight.

                                                            {
 
   if ( !muonRef.isNonnull() ) return false;
   if(!muonRef->isTrackerMuon()) return false;
 
   reco::TrackRef trackerMu = muonRef->track();
 
   if(trackerMu->ptError()/trackerMu->pt() > 0.20) return false;
 
   // this doesn't seem to be necessary on the small samples looked at, but keep it around as insurance
   if(trackerMu->pt()>20.) return false;
     
   bool isAllArbitrated = muon::isGoodMuon(*muonRef,muon::AllArbitrated);
   bool isTMLastStationAngTight = muon::isGoodMuon(*muonRef,muon::TMLastStationAngTight);
 
   bool quality = isAllArbitrated && isTMLastStationAngTight;
 
   return quality;
   
 }

bool PFMuonAlgo::isTrackerTightMuon ( const reco::PFBlockElement & elt )

static

Definition at line 211 of file PFMuonAlgo.cc.

References reco::PFBlockElementTrack::muonRef().

                                                               {
 
   const reco::PFBlockElementTrack* eltTrack 
     = dynamic_cast<const reco::PFBlockElementTrack*>(&elt);
 
   assert ( eltTrack );
   reco::MuonRef muonRef = eltTrack->muonRef();
 
   return isTrackerTightMuon(muonRef);
 
 }

bool PFMuonAlgo::isTrackerTightMuon ( const reco::MuonRef & muonRef )

static

Definition at line 313 of file PFMuonAlgo.cc.

References muon::AllArbitrated, reco::TrackBase::hitPattern(), muon::isGoodMuon(), edm::Ref< C, T, F >::isNonnull(), reco::HitPattern::numberOfValidTrackerHits(), and muon::TM2DCompatibilityTight.

                                                            {
 
   if ( !muonRef.isNonnull() ) return false;
     
   if(!muonRef->isTrackerMuon()) return false;
   
   reco::TrackRef trackerMu = muonRef->track();
   const reco::Track& track = *trackerMu;
   
   unsigned nTrackerHits =  track.hitPattern().numberOfValidTrackerHits();
   
   if(nTrackerHits<=12) return false;
   
   bool isAllArbitrated = muon::isGoodMuon(*muonRef,muon::AllArbitrated);
   
   bool isTM2DCompatibilityTight = muon::isGoodMuon(*muonRef,muon::TM2DCompatibilityTight);
   
   if(!isAllArbitrated || !isTM2DCompatibilityTight)  return false;
 
   if((trackerMu->ptError()/trackerMu->pt() > 0.10)){
     //std::cout<<" PT ERROR > 10 % "<< trackerMu->pt() <<std::endl;
     return false;
   }
   return true;
   
 }

std::vector< reco::Muon::MuonTrackTypePair > PFMuonAlgo::muonTracks	(	const reco::MuonRef &	muon,
		bool	includeSA = `false`,
		double	dpt = `1e+9`
	)

private

Definition at line 639 of file PFMuonAlgo.cc.

References reco::Muon::CombinedTrack, reco::Muon::DYT, reco::Muon::InnerTrack, edm::Ref< C, T, F >::isNonnull(), dbtoconf::out, reco::Muon::OuterTrack, reco::Muon::Picky, and reco::Muon::TPFMS.

                                                                                                              {
 
 
 
   std::vector<reco::Muon::MuonTrackTypePair> out;
 
   
   if(muon->globalTrack().isNonnull() && muon->globalTrack()->pt()>0) 
     if(muon->globalTrack()->ptError()/muon->globalTrack()->pt()<dpt)
       out.push_back(std::make_pair(muon->globalTrack(),reco::Muon::CombinedTrack));
 
   if(muon->innerTrack().isNonnull() && muon->innerTrack()->pt()>0) 
     if(muon->innerTrack()->ptError()/muon->innerTrack()->pt()<dpt)//Here Loose!@
       out.push_back(std::make_pair(muon->innerTrack(),reco::Muon::InnerTrack));
 
   bool pickyExists=false; 
   double pickyDpt=99999.; 
   if(muon->pickyTrack().isNonnull() && muon->pickyTrack()->pt()>0) {
     pickyDpt = muon->pickyTrack()->ptError()/muon->pickyTrack()->pt(); 
     if(pickyDpt<dpt) 
       out.push_back(std::make_pair(muon->pickyTrack(),reco::Muon::Picky));
     pickyExists=true;
   }
 
   bool dytExists=false; 
   double dytDpt=99999.; 
   if(muon->dytTrack().isNonnull() && muon->dytTrack()->pt()>0) {
     dytDpt = muon->dytTrack()->ptError()/muon->dytTrack()->pt(); 
     if(dytDpt<dpt) 
       out.push_back(std::make_pair(muon->dytTrack(),reco::Muon::DYT));
     dytExists=true;
   }
 
   //Magic: TPFMS is not a really good track especially under misalignment
   //IT is kind of crap because if mu system is displaced it can make a change
   //So allow TPFMS if there is no picky or the error of tpfms is better than picky
   //AND if there is no DYT or the error of tpfms is better than DYT
   if(muon->tpfmsTrack().isNonnull() && muon->tpfmsTrack()->pt()>0) { 
     double tpfmsDpt = muon->tpfmsTrack()->ptError()/muon->tpfmsTrack()->pt(); 
     if( ( (pickyExists && tpfmsDpt<pickyDpt) || (!pickyExists) ) && 
     ( (dytExists   && tpfmsDpt<dytDpt)   || (!dytExists)   ) && 
     tpfmsDpt<dpt )
       out.push_back(std::make_pair(muon->tpfmsTrack(),reco::Muon::TPFMS));
   }
 
   if(includeSA && muon->outerTrack().isNonnull())
     if(muon->outerTrack()->ptError()/muon->outerTrack()->pt()<dpt)
       out.push_back(std::make_pair(muon->outerTrack(),reco::Muon::OuterTrack));
 
   return out;
 
 }

void PFMuonAlgo::postClean ( reco::PFCandidateCollection * cands )

Definition at line 827 of file PFMuonAlgo.cc.

References reco::PFCandidate::h, i, edm::Ref< C, T, F >::isNonnull(), reco::PFCandidate::mu, reco::PFCandidate::muonRef(), patZpeak::muons, reco::LeafCandidate::p(), reco::PFCandidate::particleId(), reco::LeafCandidate::pt(), reco::LeafCandidate::px(), reco::LeafCandidate::py(), reco::PFCandidate::rawHcalEnergy(), dt_dqm_sourceclient_common_cff::reco, python.multivaluedict::sort(), and reco::PFCandidate::trackRef().

Referenced by PFAlgo::reconstructParticles().

                                                             {
   using namespace std;
   using namespace reco;
   if (!postCleaning_)
     return;
 
   //Initialize vectors
 
   if(pfCosmicsMuonCleanedCandidates_.get() )
     pfCosmicsMuonCleanedCandidates_->clear();
   else 
     pfCosmicsMuonCleanedCandidates_.reset( new reco::PFCandidateCollection );
 
   if(pfCleanedTrackerAndGlobalMuonCandidates_.get() )
     pfCleanedTrackerAndGlobalMuonCandidates_->clear();
   else 
     pfCleanedTrackerAndGlobalMuonCandidates_.reset( new reco::PFCandidateCollection );
 
   if( pfFakeMuonCleanedCandidates_.get() )
      pfFakeMuonCleanedCandidates_->clear();
   else 
      pfFakeMuonCleanedCandidates_.reset( new reco::PFCandidateCollection );
 
 
   if( pfPunchThroughMuonCleanedCandidates_.get() )
      pfPunchThroughMuonCleanedCandidates_->clear();
   else 
      pfPunchThroughMuonCleanedCandidates_.reset( new reco::PFCandidateCollection );
 
   if( pfPunchThroughHadronCleanedCandidates_.get() )
      pfPunchThroughHadronCleanedCandidates_->clear();
   else 
      pfPunchThroughHadronCleanedCandidates_.reset( new reco::PFCandidateCollection );
 
 
 
   pfPunchThroughHadronCleanedCandidates_->clear();
   
   maskedIndices_.clear();
 
   //Estimate MET and SumET
 
   estimateEventQuantities(cands);
 
   
   std::vector<int> muons;
   std::vector<int> cosmics;
   //get the muons
   for(unsigned int i=0;i<cands->size();++i)  {
     const reco::PFCandidate& cand = (*cands)[i];
     if ( cand.particleId() == reco::PFCandidate::mu )
       muons.push_back(i);
     else if ( cand.particleId() == reco::PFCandidate::h && cand.muonRef().isNonnull()) {
     //MET cleaning for muons that are not high purity and became charged hadrons 
       if (cand.pt()>100.0 &&  (!cand.trackRef()->quality(trackQuality_)) && cand.muonRef()->isGlobalMuon() && cand.rawHcalEnergy()<0.05*cand.p()) {
       maskedIndices_.push_back(i);
       pfPunchThroughMuonCleanedCandidates_->push_back(cand);
       }
     }
       
   }
   //Then sort the muon indicess by decsending pt
 
   IndexPtComparator comparator(cands);
   std::sort(muons.begin(),muons.end(),comparator);
 
 
   //first kill cosmics
   double METXCosmics=0;
   double METYCosmics=0;
   double SUMETCosmics=0.0;
 
   for(unsigned int i=0;i<muons.size();++i) {
     const PFCandidate& pfc = cands->at(muons[i]);
     double origin=0.0;
     if(vertices_->size()>0&& vertices_->at(0).isValid() && ! vertices_->at(0).isFake())
       origin = pfc.muonRef()->muonBestTrack()->dxy(vertices_->at(0).position());
 
     if( origin> cosmicRejDistance_) {
       cosmics.push_back(muons[i]);
       METXCosmics +=pfc.px();
       METYCosmics +=pfc.py();
       SUMETCosmics +=pfc.pt();
     }
   }
   double MET2Cosmics = METXCosmics*METXCosmics+METYCosmics*METYCosmics;
 
   if ( SUMETCosmics > (sumet_-sumetPU_)/eventFactorCosmics_ && MET2Cosmics < METX_*METX_+ METY_*METY_)
     for(unsigned int i=0;i<cosmics.size();++i) {
       maskedIndices_.push_back(cosmics[i]);
       pfCosmicsMuonCleanedCandidates_->push_back(cands->at(cosmics[i]));
     }
 
   //Loop on the muons candidates and clean
   for(unsigned int i=0;i<muons.size();++i) {  
     if( cleanMismeasured(cands->at(muons[i]),muons[i]))
       continue;
     cleanPunchThroughAndFakes(cands->at(muons[i]),cands,muons[i]);
   
   }
 
 
 
   //OK Now do the hard job ->remove the candidates that were cleaned 
   removeDeadCandidates(cands,maskedIndices_);
 
 
 
 
 }

void PFMuonAlgo::printMuonProperties ( const reco::MuonRef & muonRef )

static

Definition at line 496 of file PFMuonAlgo.cc.

References runregparse::combined, gather_cfg::cout, delta, reco::TrackBase::hitPattern(), muon::isGoodMuon(), edm::Ref< C, T, F >::isNonnull(), bookConverter::max, min(), reco::HitPattern::numberOfValidTrackerHits(), mathSSE::sqrt(), muon::TM2DCompatibilityLoose, muon::TM2DCompatibilityTight, muon::TMLastStationAngLoose, muon::TMLastStationAngTight, muon::TMLastStationLoose, muon::TMLastStationOptimizedBarrelLowPtLoose, muon::TMLastStationOptimizedBarrelLowPtTight, muon::TMLastStationOptimizedLowPtLoose, muon::TMLastStationOptimizedLowPtTight, muon::TMLastStationTight, muon::TMOneStationLoose, muon::TMOneStationTight, and patCandidatesForDimuonsSequences_cff::tracker.

                                                          {
   
   if ( !muonRef.isNonnull() ) return;
   
   bool isGL = muonRef->isGlobalMuon();
   bool isTR = muonRef->isTrackerMuon();
   bool isST = muonRef->isStandAloneMuon();
 
   std::cout<<" GL: "<<isGL<<" TR: "<<isTR<<" ST: "<<isST<<std::endl;
   std::cout<<" nMatches "<<muonRef->numberOfMatches()<<std::endl;
   
   if ( muonRef->isGlobalMuon() ){
     reco::TrackRef combinedMu = muonRef->combinedMuon();
     std::cout<<" GL,  pt: " << combinedMu->pt() 
     << " +/- " << combinedMu->ptError()/combinedMu->pt() 
          << " chi**2 GBL : " << combinedMu->normalizedChi2()<<std::endl;
     std::cout<< " Total Muon Hits : " << combinedMu->hitPattern().numberOfValidMuonHits()
          << "/" << combinedMu->hitPattern().numberOfLostMuonHits()
     << " DT Hits : " << combinedMu->hitPattern().numberOfValidMuonDTHits()
     << "/" << combinedMu->hitPattern().numberOfLostMuonDTHits()
     << " CSC Hits : " << combinedMu->hitPattern().numberOfValidMuonCSCHits()
     << "/" << combinedMu->hitPattern().numberOfLostMuonCSCHits()
     << " RPC Hits : " << combinedMu->hitPattern().numberOfValidMuonRPCHits()
          << "/" << combinedMu->hitPattern().numberOfLostMuonRPCHits()<<std::endl;
 
     std::cout<<"  # of Valid Tracker Hits "<<combinedMu->hitPattern().numberOfValidTrackerHits()<<std::endl;
     std::cout<<"  # of Valid Pixel Hits "<<combinedMu->hitPattern().numberOfValidPixelHits()<<std::endl;
   }
   if ( muonRef->isStandAloneMuon() ){
     reco::TrackRef standAloneMu = muonRef->standAloneMuon();
     std::cout<<" ST,  pt: " << standAloneMu->pt() 
     << " +/- " << standAloneMu->ptError()/standAloneMu->pt() 
     << " eta : " << standAloneMu->eta()  
     << " DT Hits : " << standAloneMu->hitPattern().numberOfValidMuonDTHits()
     << "/" << standAloneMu->hitPattern().numberOfLostMuonDTHits()
     << " CSC Hits : " << standAloneMu->hitPattern().numberOfValidMuonCSCHits()
     << "/" << standAloneMu->hitPattern().numberOfLostMuonCSCHits()
     << " RPC Hits : " << standAloneMu->hitPattern().numberOfValidMuonRPCHits()
     << "/" << standAloneMu->hitPattern().numberOfLostMuonRPCHits()
          << " chi**2 STA : " << standAloneMu->normalizedChi2()<<std::endl;
       }
 
 
   if ( muonRef->isTrackerMuon() ){
     reco::TrackRef trackerMu = muonRef->track();
     const reco::Track& track = *trackerMu;
     std::cout<<" TR,  pt: " << trackerMu->pt() 
     << " +/- " << trackerMu->ptError()/trackerMu->pt() 
          << " chi**2 TR : " << trackerMu->normalizedChi2()<<std::endl;    
     std::cout<<" nTrackerHits "<<track.hitPattern().numberOfValidTrackerHits()<<std::endl;
     std::cout<< "TMLastStationAngLoose               "
     << muon::isGoodMuon(*muonRef,muon::TMLastStationAngLoose) << std::endl       
     << "TMLastStationAngTight               "
     << muon::isGoodMuon(*muonRef,muon::TMLastStationAngTight) << std::endl          
     << "TMLastStationLoose               "
     << muon::isGoodMuon(*muonRef,muon::TMLastStationLoose) << std::endl       
     << "TMLastStationTight               "
     << muon::isGoodMuon(*muonRef,muon::TMLastStationTight) << std::endl          
     << "TMOneStationLoose                "
     << muon::isGoodMuon(*muonRef,muon::TMOneStationLoose) << std::endl       
     << "TMOneStationTight                "
     << muon::isGoodMuon(*muonRef,muon::TMOneStationTight) << std::endl       
     << "TMLastStationOptimizedLowPtLoose " 
     << muon::isGoodMuon(*muonRef,muon::TMLastStationOptimizedLowPtLoose) << std::endl
     << "TMLastStationOptimizedLowPtTight " 
     << muon::isGoodMuon(*muonRef,muon::TMLastStationOptimizedLowPtTight) << std::endl 
     << "TMLastStationOptimizedBarrelLowPtLoose " 
     << muon::isGoodMuon(*muonRef,muon::TMLastStationOptimizedBarrelLowPtLoose) << std::endl
     << "TMLastStationOptimizedBarrelLowPtTight " 
     << muon::isGoodMuon(*muonRef,muon::TMLastStationOptimizedBarrelLowPtTight) << std::endl 
     << std::endl;
 
   }
 
   std::cout<< "TM2DCompatibilityLoose           "
       << muon::isGoodMuon(*muonRef,muon::TM2DCompatibilityLoose) << std::endl 
       << "TM2DCompatibilityTight           "
        << muon::isGoodMuon(*muonRef,muon::TM2DCompatibilityTight) << std::endl;
 
 
 
   if (      muonRef->isGlobalMuon() &&  muonRef->isTrackerMuon() &&  muonRef->isStandAloneMuon() ){
     reco::TrackRef combinedMu = muonRef->combinedMuon();
     reco::TrackRef trackerMu = muonRef->track();
     reco::TrackRef standAloneMu = muonRef->standAloneMuon();
     
     double sigmaCombined = combinedMu->ptError()/(combinedMu->pt()*combinedMu->pt());    
     double sigmaTracker = trackerMu->ptError()/(trackerMu->pt()*trackerMu->pt());    
     double sigmaStandAlone = standAloneMu->ptError()/(standAloneMu->pt()*standAloneMu->pt());    
     
     bool combined = combinedMu->ptError()/combinedMu->pt() < 0.20;   
     bool tracker = trackerMu->ptError()/trackerMu->pt() < 0.20;      
     bool standAlone = standAloneMu->ptError()/standAloneMu->pt() < 0.20;     
   
     double delta1 =  combined && tracker ?   
       fabs(1./combinedMu->pt() -1./trackerMu->pt())      
       /sqrt(sigmaCombined*sigmaCombined + sigmaTracker*sigmaTracker) : 100.;     
     double delta2 = combined && standAlone ?     
       fabs(1./combinedMu->pt() -1./standAloneMu->pt())   
       /sqrt(sigmaCombined*sigmaCombined + sigmaStandAlone*sigmaStandAlone) : 100.;   
     double delta3 = standAlone && tracker ?      
       fabs(1./standAloneMu->pt() -1./trackerMu->pt())    
       /sqrt(sigmaStandAlone*sigmaStandAlone + sigmaTracker*sigmaTracker) : 100.;     
     
     double delta =   
       standAloneMu->hitPattern().numberOfValidMuonDTHits()+      
       standAloneMu->hitPattern().numberOfValidMuonCSCHits() > 0 ?    
       std::min(delta3,std::min(delta1,delta2)) : std::max(delta3,std::max(delta1,delta2));   
 
     std::cout << "delta = " << delta << " delta1 "<<delta1<<" delta2 "<<delta2<<" delta3 "<<delta3<<std::endl;   
     
     double ratio =   
       combinedMu->ptError()/combinedMu->pt()     
       / (trackerMu->ptError()/trackerMu->pt());      
     //if ( ratio > 2. && delta < 3. ) std::cout << "ALARM ! " << ratio << ", " << delta << std::endl;
     std::cout<<" ratio "<<ratio<<" combined mu pt "<<combinedMu->pt()<<std::endl;
     //bool quality3 =  ( combinedMu->pt() < 50. || ratio < 2. ) && delta <  3.;
 
 
   }
 
     double sumPtR03 = muonRef->isolationR03().sumPt;
     double emEtR03 = muonRef->isolationR03().emEt;
     double hadEtR03 = muonRef->isolationR03().hadEt;    
     double relIsoR03 = (sumPtR03 + emEtR03 + hadEtR03)/muonRef->pt();
     double sumPtR05 = muonRef->isolationR05().sumPt;
     double emEtR05 = muonRef->isolationR05().emEt;
     double hadEtR05 = muonRef->isolationR05().hadEt;    
     double relIsoR05 = (sumPtR05 + emEtR05 + hadEtR05)/muonRef->pt();
     std::cout<<" 0.3 Radion Rel Iso: "<<relIsoR03<<" sumPt "<<sumPtR03<<" emEt "<<emEtR03<<" hadEt "<<hadEtR03<<std::endl;
     std::cout<<" 0.5 Radion Rel Iso: "<<relIsoR05<<" sumPt "<<sumPtR05<<" emEt "<<emEtR05<<" hadEt "<<hadEtR05<<std::endl;
   return;
 
 }

bool PFMuonAlgo::reconstructMuon	(	reco::PFCandidate &	candidate,
		const reco::MuonRef &	muon,
		bool	allowLoose = `false`
	)

Definition at line 701 of file PFMuonAlgo.cc.

References reco::Muon::InnerTrack, muon::isLooseMuon(), reco::isMuon(), edm::Ref< C, T, F >::isNonnull(), dt_dqm_sourceclient_common_cff::reco, reco::PFCandidate::setMuonRef(), and muon::tevOptimized().

Referenced by PFAlgo::reconstructTrack().

                                                                                                      {
     using namespace std;
     using namespace reco;
 
     if (!muon.isNonnull())
       return false;
 
     
 
     bool isMu=false;
 
     if(allowLoose) 
       isMu = isMuon(muon) || isLooseMuon(muon);
     else
       isMu = isMuon(muon);
 
     if( !isMu)
       return false;
 
     //get the valid tracks(without standalone except we allow loose muons)
     //MIKE: Here we need to be careful. If we have a muon inside a dense 
     //jet environment often the track is badly measured. In this case 
     //we should not apply Dpt/Pt<1
  
     std::vector<reco::Muon::MuonTrackTypePair> validTracks = goodMuonTracks(muon);
     if (!allowLoose)
       validTracks = goodMuonTracks(muon);
     else
       validTracks = muonTracks(muon);
 
     if( validTracks.size() ==0)
       return false;
 
 
 
     //check what is the track used.Rerun TuneP
     reco::Muon::MuonTrackTypePair bestTrackPair = muon::tevOptimized(*muon);
     
     TrackRef bestTrack = bestTrackPair.first;
     MuonTrackType trackType = bestTrackPair.second;
 
 
 
     MuonTrackTypePair trackPairWithSmallestError = getTrackWithSmallestError(validTracks);
     TrackRef trackWithSmallestError = trackPairWithSmallestError.first;
     
     if( trackType == reco::Muon::InnerTrack && 
     (!bestTrack->quality(trackQuality_) ||
      bestTrack->ptError()/bestTrack->pt()> errorCompScale_*trackWithSmallestError->ptError()/trackWithSmallestError->pt() )) {
       bestTrack = trackWithSmallestError;
       trackType = trackPairWithSmallestError.second;
     }
     else if (trackType != reco::Muon::InnerTrack &&
          bestTrack->ptError()/bestTrack->pt()> errorCompScale_*trackWithSmallestError->ptError()/trackWithSmallestError->pt())  {
       bestTrack = trackWithSmallestError;
       trackType = trackPairWithSmallestError.second;
       
     }
 
 
     changeTrack(candidate,std::make_pair(bestTrack,trackType));
     candidate.setMuonRef( muon );
 
     return true;
 }

void PFMuonAlgo::removeDeadCandidates	(	reco::PFCandidateCollection *	obj,
		const std::vector< unsigned int > &	indices
	)

private

Definition at line 1250 of file PFMuonAlgo.cc.

References i, and N.

 {
   size_t N = indices.size();
   size_t collSize = obj->size();
 
   for (size_t i = 0 ; i < N ; ++i)
     obj->at(indices.at(i)) = obj->at(collSize-i-1);
 
   obj->resize(collSize - indices.size());
 }

void PFMuonAlgo::setInputsForCleaning ( const reco::VertexCollection * vertices )

Definition at line 1153 of file PFMuonAlgo.cc.

Referenced by PFAlgo::setPFVertexParameters().

                                                                            {
   vertices_ = vertices;
 }

void PFMuonAlgo::setParameters ( const edm::ParameterSet & iConfig )

Definition at line 29 of file PFMuonAlgo.cc.

References edm::ParameterSet::exists(), edm::ParameterSet::getParameter(), reco::TrackBase::qualityByName(), and AlCaHLTBitMon_QueryRunRegistry::string.

Referenced by PFAlgo::setPFMuonAndFakeParameters().

 {
 
   if(iConfig.exists("maxDPtOPt"))
     maxDPtOPt_ = iConfig.getParameter<double>("maxDPtOPt");
   else
     maxDPtOPt_=1.0;
 
   if(iConfig.exists("minTrackerHits"))
     minTrackerHits_ = iConfig.getParameter<int>("minTrackerHits");
   else
     minTrackerHits_ = 8;
 
   if(iConfig.exists("minPixelHits"))
     minPixelHits_ = iConfig.getParameter<int>("minPixelHits");
   else
     minPixelHits_ = 1;
 
   if(iConfig.exists("trackQuality"))
     trackQuality_  = reco::TrackBase::qualityByName(iConfig.getParameter<std::string>("trackQuality"));
   else
     trackQuality_  = reco::TrackBase::qualityByName("highPurity");
 
   if(iConfig.exists("ptErrorScale"))
     errorCompScale_ = iConfig.getParameter<double>("ptErrorScale");
   else
     errorCompScale_ = 4.;
 
   if(iConfig.exists("eventFractionForCleaning"))
     eventFractionCleaning_ = iConfig.getParameter<double>("eventFractionForCleaning");
   else
     eventFractionCleaning_ = 0.75;
 
   if(iConfig.exists("dzPV"))
     dzPV_ = iConfig.getParameter<double>("dzPV");
   else
     dzPV_ = 0.2;
 
   if(iConfig.exists("postMuonCleaning"))
     postCleaning_ = iConfig.getParameter<bool>("postMuonCleaning");
   else
     postCleaning_ = false; //Disable by default (for HLT)
 
   if(iConfig.exists("minPtForPostCleaning"))
     minPostCleaningPt_ = iConfig.getParameter<double>("minPtForPostCleaning");
   else
     minPostCleaningPt_ = 20.;
 
   if(iConfig.exists("eventFactorForCosmics"))
     eventFactorCosmics_ = iConfig.getParameter<double>("eventFactorForCosmics");
   else
     eventFactorCosmics_ = 10.;
 
 
   if(iConfig.exists("metSignificanceForCleaning"))
     metSigForCleaning_ = iConfig.getParameter<double>("metSignificanceForCleaning");
   else
     metSigForCleaning_ = 3.;
 
   if(iConfig.exists("metSignificanceForRejection"))
     metSigForRejection_ = iConfig.getParameter<double>("metSignificanceForRejection");
   else
     metSigForRejection_ = 4.;
 
   if(iConfig.exists("metFactorForCleaning"))
     metFactorCleaning_ = iConfig.getParameter<double>("metFactorForCleaning");
   else
     metFactorCleaning_ = 4.;
 
   if(iConfig.exists("eventFractionForRejection"))
     eventFractionRejection_ = iConfig.getParameter<double>("eventFractionForRejection");
   else
     eventFractionRejection_ = 0.75;
 
   if(iConfig.exists("metFactorForRejection"))
     metFactorRejection_ = iConfig.getParameter<double>("metFactorForRejection");
   else
     metFactorRejection_ =4.;
 
   if(iConfig.exists("metFactorForHighEta"))
     metFactorHighEta_ = iConfig.getParameter<double>("metFactorForHighEta");
   else
     metFactorHighEta_=4;
 
   if(iConfig.exists("ptFactorForHighEta"))
     ptFactorHighEta_ = iConfig.getParameter<double>("ptFactorForHighEta");
   else
     ptFactorHighEta_ = 2.;
 
   if(iConfig.exists("metFactorForFakes"))
     metFactorFake_ = iConfig.getParameter<double>("metFactorForFakes");
   else
     metFactorFake_ = 4.;
 
   if(iConfig.exists("minMomentumForPunchThrough"))
     minPunchThroughMomentum_ = iConfig.getParameter<double>("minMomentumForPunchThrough");
   else
     minPunchThroughMomentum_=100.;
 
   if(iConfig.exists("minEnergyForPunchThrough"))
     minPunchThroughEnergy_ = iConfig.getParameter<double>("minEnergyForPunchThrough");
   else
     minPunchThroughEnergy_ = 100.;
 
   if(iConfig.exists("punchThroughFactor"))
     punchThroughFactor_ = iConfig.getParameter<double>("punchThroughFactor");
   else
     punchThroughFactor_ = 3.;
 
   if(iConfig.exists("punchThroughMETFactor"))
     punchThroughMETFactor_ = iConfig.getParameter<double>("punchThroughMETFactor");
   else
     punchThroughMETFactor_ = 4.;
 
   if(iConfig.exists("cosmicRejectionDistance"))
     cosmicRejDistance_ = iConfig.getParameter<double>("cosmicRejectionDistance");
   else
     cosmicRejDistance_ = 1.0;
 }

std::vector< reco::Muon::MuonTrackTypePair > PFMuonAlgo::tracksPointingAtMET ( const std::vector< MuonTrackTypePair > & tracks )

private

Definition at line 1133 of file PFMuonAlgo.cc.

References i, and funct::pow().

                                                                                   {
   std::vector<reco::Muon::MuonTrackTypePair> outputTracks;
 
 
   double newMET2=0.0;
 
   for( unsigned int i=0;i<tracks.size();++i) {
     //calculate new SUM ET and MET2
     newMET2  = pow(METX_+tracks.at(i).first->px(),2)+pow(METY_+tracks.at(i).first->py(),2);
     
     if(newMET2<(METX_*METX_+METY_*METY_)/metFactorCleaning_)
       outputTracks.push_back(tracks.at(i));
   }
   
   
   return outputTracks;
 }

std::vector< reco::Muon::MuonTrackTypePair > PFMuonAlgo::tracksWithBetterMET	(	const std::vector< MuonTrackTypePair > &	tracks,
		const reco::PFCandidate &	pfc
	)

private

Definition at line 1105 of file PFMuonAlgo.cc.

References i, funct::pow(), reco::LeafCandidate::pt(), reco::LeafCandidate::px(), reco::LeafCandidate::py(), and mathSSE::sqrt().

                                                                                                               {
   std::vector<reco::Muon::MuonTrackTypePair> outputTracks;
 
   double METNOX  = METX_ - pfc.px();
   double METNOY  = METY_ - pfc.py();
   double SUMETNO = sumet_  -pfc.pt(); 
   double MET2 = METX_*METX_+METY_*METY_;
   double newMET2=0.0;
   double newSUMET=0.0;
   double METSIG = sqrt(MET2)/sqrt(sumet_-sumetPU_);
 
 
   if(METSIG>metSigForCleaning_)
   for( unsigned int i=0;i<tracks.size();++i) {
     //calculate new SUM ET and MET2
     newSUMET = SUMETNO+tracks.at(i).first->pt()-sumetPU_;
     newMET2  = pow(METNOX+tracks.at(i).first->px(),2)+pow(METNOY+tracks.at(i).first->py(),2);
     
     if(newSUMET/(sumet_-sumetPU_)>eventFractionCleaning_ &&  newMET2<MET2/metFactorCleaning_)
       outputTracks.push_back(tracks.at(i));
   }
   
   
   return outputTracks;
 }

std::auto_ptr<reco::PFCandidateCollection>& PFMuonAlgo::transferAddedMuonCandidates ( )

inline

Definition at line 101 of file PFMuonAlgo.h.

References pfAddedMuonCandidates_.

                                                                         {
     return pfAddedMuonCandidates_;
   }

std::auto_ptr<reco::PFCandidateCollection>& PFMuonAlgo::transferCleanedCosmicCandidates ( )

inline

Definition at line 81 of file PFMuonAlgo.h.

References pfCosmicsMuonCleanedCandidates_.

                                                                             {
     return pfCosmicsMuonCleanedCandidates_;
   }

std::auto_ptr<reco::PFCandidateCollection>& PFMuonAlgo::transferCleanedFakeCandidates ( )

inline

Definition at line 89 of file PFMuonAlgo.h.

References pfFakeMuonCleanedCandidates_.

                                                                           {
     return pfFakeMuonCleanedCandidates_;
   }

std::auto_ptr<reco::PFCandidateCollection>& PFMuonAlgo::transferCleanedTrackerAndGlobalCandidates ( )

inline

Definition at line 85 of file PFMuonAlgo.h.

References pfCleanedTrackerAndGlobalMuonCandidates_.

                                                                                       {
     return pfCleanedTrackerAndGlobalMuonCandidates_;
   }

std::auto_ptr<reco::PFCandidateCollection>& PFMuonAlgo::transferPunchThroughCleanedHadronCandidates ( )

inline

Definition at line 97 of file PFMuonAlgo.h.

References pfPunchThroughHadronCleanedCandidates_.

                                                                                         {
     return pfPunchThroughHadronCleanedCandidates_;
   }

std::auto_ptr<reco::PFCandidateCollection>& PFMuonAlgo::transferPunchThroughCleanedMuonCandidates ( )

inline

Definition at line 93 of file PFMuonAlgo.h.

References pfPunchThroughMuonCleanedCandidates_.

                                                                                       {
     return pfPunchThroughMuonCleanedCandidates_;
   }

Member Data Documentation

double PFMuonAlgo::cosmicRejDistance_

private

Definition at line 178 of file PFMuonAlgo.h.

double PFMuonAlgo::dzPV_

private

Definition at line 162 of file PFMuonAlgo.h.

double PFMuonAlgo::errorCompScale_

private

Definition at line 159 of file PFMuonAlgo.h.

double PFMuonAlgo::eventFactorCosmics_

private

Definition at line 165 of file PFMuonAlgo.h.

double PFMuonAlgo::eventFractionCleaning_

private

Definition at line 160 of file PFMuonAlgo.h.

double PFMuonAlgo::eventFractionRejection_

private

Definition at line 169 of file PFMuonAlgo.h.

std::vector<unsigned int > PFMuonAlgo::maskedIndices_

private

Definition at line 145 of file PFMuonAlgo.h.

double PFMuonAlgo::maxDPtOPt_

private

Definition at line 154 of file PFMuonAlgo.h.

double PFMuonAlgo::metFactorCleaning_

private

Definition at line 168 of file PFMuonAlgo.h.

double PFMuonAlgo::metFactorFake_

private

Definition at line 173 of file PFMuonAlgo.h.

double PFMuonAlgo::metFactorHighEta_

private

Definition at line 171 of file PFMuonAlgo.h.

double PFMuonAlgo::metFactorRejection_

private

Definition at line 170 of file PFMuonAlgo.h.

double PFMuonAlgo::metSigForCleaning_

private

Definition at line 166 of file PFMuonAlgo.h.

double PFMuonAlgo::metSigForRejection_

private

Definition at line 167 of file PFMuonAlgo.h.

double PFMuonAlgo::METX_

private

Definition at line 181 of file PFMuonAlgo.h.

double PFMuonAlgo::METY_

private

Definition at line 182 of file PFMuonAlgo.h.

int PFMuonAlgo::minPixelHits_

private

Definition at line 156 of file PFMuonAlgo.h.

double PFMuonAlgo::minPostCleaningPt_

private

Definition at line 164 of file PFMuonAlgo.h.

double PFMuonAlgo::minPunchThroughEnergy_

private

Definition at line 175 of file PFMuonAlgo.h.

double PFMuonAlgo::minPunchThroughMomentum_

private

Definition at line 174 of file PFMuonAlgo.h.

int PFMuonAlgo::minTrackerHits_

private

Definition at line 155 of file PFMuonAlgo.h.

std::auto_ptr< reco::PFCandidateCollection > PFMuonAlgo::pfAddedMuonCandidates_

private

the collection of added muon candidates

Definition at line 143 of file PFMuonAlgo.h.

Referenced by transferAddedMuonCandidates().

std::auto_ptr< reco::PFCandidateCollection > PFMuonAlgo::pfCleanedTrackerAndGlobalMuonCandidates_

private

the collection of tracker/global cleaned muon candidates

Definition at line 135 of file PFMuonAlgo.h.

Referenced by transferCleanedTrackerAndGlobalCandidates().

std::auto_ptr< reco::PFCandidateCollection > PFMuonAlgo::pfCosmicsMuonCleanedCandidates_

private

the collection of cosmics cleaned muon candidates

Definition at line 133 of file PFMuonAlgo.h.

Referenced by transferCleanedCosmicCandidates().

std::auto_ptr< reco::PFCandidateCollection > PFMuonAlgo::pfFakeMuonCleanedCandidates_

private

the collection of fake cleaned muon candidates

Definition at line 137 of file PFMuonAlgo.h.

Referenced by transferCleanedFakeCandidates().

std::auto_ptr< reco::PFCandidateCollection > PFMuonAlgo::pfPunchThroughHadronCleanedCandidates_

private

the collection of punch-through cleaned neutral hadron candidates

Definition at line 141 of file PFMuonAlgo.h.

Referenced by transferPunchThroughCleanedHadronCandidates().

std::auto_ptr< reco::PFCandidateCollection > PFMuonAlgo::pfPunchThroughMuonCleanedCandidates_

private

the collection of punch-through cleaned muon candidates

Definition at line 139 of file PFMuonAlgo.h.

Referenced by transferPunchThroughCleanedMuonCandidates().

bool PFMuonAlgo::postCleaning_

private

Definition at line 163 of file PFMuonAlgo.h.

double PFMuonAlgo::ptFactorHighEta_

private

Definition at line 172 of file PFMuonAlgo.h.

double PFMuonAlgo::punchThroughFactor_

private

Definition at line 176 of file PFMuonAlgo.h.

double PFMuonAlgo::punchThroughMETFactor_

private

Definition at line 177 of file PFMuonAlgo.h.

double PFMuonAlgo::sumet_

private

Definition at line 180 of file PFMuonAlgo.h.

double PFMuonAlgo::sumetPU_

private

Definition at line 161 of file PFMuonAlgo.h.

reco::TrackBase::TrackQuality PFMuonAlgo::trackQuality_

private

Definition at line 157 of file PFMuonAlgo.h.

const reco::VertexCollection* PFMuonAlgo::vertices_

private

Definition at line 149 of file PFMuonAlgo.h.

Classes

Public Member Functions

Static Public Member Functions

Private Types

Private Member Functions

Private Attributes

Detailed Description

Member Typedef Documentation

Constructor & Destructor Documentation

Member Function Documentation

Member Data Documentation