PFMuonAlgo Class Reference

#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.

{;}

Member Function Documentation

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

Definition at line 948 of file PFMuonAlgo.cc.

References DeDxDiscriminatorTools::charge(), 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 754 of file PFMuonAlgo.cc.

References reco::Muon::CombinedTrack, relval_parameters_module::energy, reco::Muon::InnerTrack, reco::PFCandidate::mu, reco::Muon::Picky, L1Trigger_dataformats::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.13957*0.13957);

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

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

private

Definition at line 1056 of file PFMuonAlgo.cc.

References reco::PFCandidate::muonRef(), reco::LeafCandidate::pt(), reco::LeafCandidate::px(), reco::LeafCandidate::py(), L1Trigger_dataformats::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 1167 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(), max(), reco::PFCandidate::muonRef(), reco::LeafCandidate::p(), reco::PFCandidate::particleId(), reco::LeafCandidate::pt(), reco::PFCandidate::rawEcalEnergy(), reco::PFCandidate::rawHcalEnergy(), and L1Trigger_dataformats::reco.

                                                                                                                {
  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;
      for ( unsigned iele = 1; iele < eleInBlocks.size(); ++iele ) { 
    indexMuon = eleInBlocks[iele].second;
    break;
      }
      
      // 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 790 of file PFMuonAlgo.cc.

References abs, and i.

{
  //SUM ET from PU
  sumetPU_ = 0.0;
  METX_=0.;
  METY_=0.;
  for (unsigned short i=1 ;i<vertices_->size();++i ) {
    if ( !vertices_->at(i).isValid() || vertices_->at(i).isFake() ) continue; 
    vertices_->at(i);
    for ( reco::Vertex::trackRef_iterator itr = vertices_->at(i).tracks_begin();
      itr <  vertices_->at(i).tracks_end(); ++itr ) { 
      sumetPU_ += (*itr)->pt();
    }
  }
  sumetPU_ /= 0.65;
  //SUM ET and MET
  sumet_=0.0;
  double METXCh=0.0;
  double METYCh=0.0;
  double METXNeut=0.0;
  double METYNeut=0.0;


  for(reco::PFCandidateCollection::const_iterator i = pfc->begin();i!=pfc->end();++i) {
    sumet_+=i->pt();

    if (vertices_->size()>0 && vertices_->at(0).isValid()&& !vertices_->at(0).isFake()) {
      //If charged and from PV or muon
      if( (i->charge() !=0 && i->trackRef().isNonnull() && vertices_->size()>0&& i->trackRef()->dz(vertices_->at(0).position())<dzPV_)||(abs(i->pdgId())==13)) {
    METXCh+=i->px();
    METYCh+=i->py();
      }
      //If charged and not from PV(assume there is a neutral balancing it)
      else if( i->charge() !=0 && i->trackRef().isNonnull() && i->trackRef()->dz(vertices_->at(0).position())>dzPV_) {
    METXNeut-=i->px();
    METYNeut-=i->py();
      }
      //Neutral
      else if( !(i->charge() !=0 && i->trackRef().isNonnull())) {
    METXNeut+=i->px();
    METYNeut+=i->py();
      }
    } //else if we dont have a vertex make standard PFMET
    else {
    METXCh+=i->px();
    METYCh+=i->py();
    }
    METX_ = (METXCh+METXNeut);
    METY_ = (METYCh+METYNeut);
  }

}

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

private

Definition at line 1036 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 781 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 PFAlgo::processBlock(), and PFBlockAlgo::setInput().

                                                       {


  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 PFAlgo::processBlock(), PFBlockAlgo::setInput(), PFElectronAlgo::SetLinks(), and PFEGammaAlgo::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::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()) 
    if(muon->globalTrack()->ptError()/muon->globalTrack()->pt()<dpt)
      out.push_back(std::make_pair(muon->globalTrack(),reco::Muon::CombinedTrack));

  if(muon->innerTrack().isNonnull()) 
    if(muon->innerTrack()->ptError()/muon->innerTrack()->pt()<dpt)//Here Loose!@
      out.push_back(std::make_pair(muon->innerTrack(),reco::Muon::InnerTrack));

  bool pickyExists=false; 
  if(muon->pickyTrack().isNonnull()) {
    if(muon->pickyTrack()->ptError()/muon->pickyTrack()->pt()<dpt) 
      out.push_back(std::make_pair(muon->pickyTrack(),reco::Muon::Picky));
    pickyExists=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
  if(muon->tpfmsTrack().isNonnull() && ((pickyExists && muon->tpfmsTrack()->ptError()/muon->tpfmsTrack()->pt()<muon->pickyTrack()->ptError()/muon->pickyTrack()->pt())||(!pickyExists)) ) 
    if(muon->tpfmsTrack()->ptError()/muon->tpfmsTrack()->pt()<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 846 of file PFMuonAlgo.cc.

References i, reco::PFCandidate::mu, reco::PFCandidate::muonRef(), patZpeak::muons, reco::LeafCandidate::pt(), reco::LeafCandidate::px(), reco::LeafCandidate::py(), L1Trigger_dataformats::reco, and python.multivaluedict::sort().

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) 
    if ( cands->at(i).particleId() == reco::PFCandidate::mu )
      muons.push_back(i);

  //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)
      pfCosmicsMuonCleanedCandidates_->push_back(cands->at(muons[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(), 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 685 of file PFMuonAlgo.cc.

References reco::Muon::InnerTrack, muon::isLooseMuon(), reco::isMuon(), edm::Ref< C, T, F >::isNonnull(), L1Trigger_dataformats::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 1262 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 1163 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 1143 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 1115 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.