CandidateSimMuonMatcher.cc File Reference

#include "L1Trigger/L1TMuonOverlapPhase1/interface/Tools/CandidateSimMuonMatcher.h"
#include "L1Trigger/L1TMuon/interface/MicroGMTConfiguration.h"
#include "FWCore/MessageLogger/interface/MessageLogger.h"
#include "FWCore/Framework/interface/Event.h"
#include "FWCore/Framework/interface/EventSetup.h"
#include "DataFormats/GeometrySurface/interface/BoundCylinder.h"
#include "DataFormats/GeometrySurface/interface/SimpleCylinderBounds.h"
#include "DataFormats/GeometrySurface/interface/BoundDisk.h"
#include "DataFormats/GeometrySurface/interface/SimpleDiskBounds.h"
#include "DataFormats/MuonDetId/interface/RPCDetId.h"
#include "TrackingTools/TrajectoryParametrization/interface/GlobalTrajectoryParameters.h"
#include "TFile.h"
#include "TH1D.h"

Go to the source code of this file.

Functions
double	foldPhi (double phi)
double	hwGmtPhiToGlobalPhi (int phi)
float	normal_pdf (float x, float m, float s)
bool	simTrackIsMuonInOmtf (const SimTrack &simTrack)
bool	simTrackIsMuonInOmtfBx0 (const SimTrack &simTrack)
bool	trackingParticleIsMuonInOmtfBx0 (const TrackingParticle &trackingParticle)
bool	trackingParticleIsMuonInOmtfEvent0 (const TrackingParticle &trackingParticle)

Function Documentation

double foldPhi

(

double

phi

)

Definition at line 33 of file CandidateSimMuonMatcher.cc.

References M_PI, and PVValHelper::phi.

Referenced by CandidateSimMuonMatcher::match().

                           {
  if (phi > M_PI)
    return (phi - 2 * M_PI);
  else if (phi < -M_PI)
    return (phi + 2 * M_PI);

  return phi;
}

double hwGmtPhiToGlobalPhi

(

int

phi

)

Definition at line 28 of file CandidateSimMuonMatcher.cc.

References M_PI.

Referenced by CandidateSimMuonMatcher::ghostBust(), and CandidateSimMuonMatcher::match().

                                    {
  double phiGmtUnit = 2. * M_PI / 576.;
  return phi * phiGmtUnit;
}

float normal_pdf	(	float	x,
		float	m,
		float	s
	)

Definition at line 342 of file CandidateSimMuonMatcher.cc.

References a, funct::exp(), and visualization-live-secondInstance_cfg::m.

Referenced by CandidateSimMuonMatcher::match().

                                            {
  static const float inv_sqrt_2pi = 0.3989422804014327;
  float a = (x - m) / s;

  return inv_sqrt_2pi / s * std::exp(-0.5 * a * a);
}

bool simTrackIsMuonInOmtf

(

const SimTrack &

simTrack

)

Definition at line 87 of file CandidateSimMuonMatcher.cc.

References funct::abs(), LogTrace, CoreSimTrack::momentum(), and CoreSimTrack::type().

Referenced by simTrackIsMuonInOmtfBx0().

                                                    {
  if (abs(simTrack.type()) == 13 || abs(simTrack.type()) == 1000015) {  //|| tpPtr->pt() > 20 //todo 1000015 is stau
    //only muons
  } else
    return false;

  //in the overlap, the propagation of muons with pt less then ~3.2 fails - the actual threshold depends slightly on eta,
  if (simTrack.momentum().pt() < 2.5)
    return false;

  LogTrace("l1tOmtfEventPrint") << "simTrackIsMuonInOmtf, simTrack type " << std::setw(3) << simTrack.type() << " pt "
                                << std::setw(9) << simTrack.momentum().pt() << " eta " << std::setw(9)
                                << simTrack.momentum().eta() << " phi " << std::setw(9) << simTrack.momentum().phi()
                                << std::endl;

  //higher margin for matching must be used than actual OMTF region (i.e. 0.82 - 1.24),
  //otherwise many candidates are marked as ghosts
  //if( (fabs(simTrack.momentum().eta()) >= 0.82 ) && (fabs(simTrack.momentum().eta()) <= 1.24) ) {
  if ((fabs(simTrack.momentum().eta()) >= 0.72) && (fabs(simTrack.momentum().eta()) <= 1.3)) {
  } else
    return false;

  return true;
}

bool simTrackIsMuonInOmtfBx0

(

const SimTrack &

simTrack

)

Definition at line 112 of file CandidateSimMuonMatcher.cc.

References EncodedEventId::bunchCrossing(), CoreSimTrack::eventId(), and simTrackIsMuonInOmtf().

Referenced by CandidateSimMuonMatcher::observeEventEnd().

                                                       {
  if (simTrack.eventId().bunchCrossing() != 0)
    return false;

  return simTrackIsMuonInOmtf(simTrack);
}

bool trackingParticleIsMuonInOmtfBx0

(

const TrackingParticle &

trackingParticle

)

Definition at line 119 of file CandidateSimMuonMatcher.cc.

References funct::abs(), EncodedEventId::bunchCrossing(), EncodedEventId::event(), TrackingParticle::eventId(), TrackingParticle::g4Tracks(), edm::Ref< C, T, F >::isNonnull(), LogTrace, TrackingParticle::momentum(), TrackingParticle::parentVertex(), TrackingParticle::pdgId(), and TrackingParticle::pt().

Referenced by CandidateSimMuonMatcher::observeEventEnd(), and trackingParticleIsMuonInOmtfEvent0().

                                                                               {
  /*if(trackingParticle.eventId().event() != 0)
      LogTrace("l1tOmtfEventPrint") <<"trackingParticleIsMuonInOmtfBx0, pdgId "<<std::setw(3)<<trackingParticle.pdgId()<<" pt "<<std::setw(9)<<trackingParticle.pt()
              <<" eta "<<std::setw(9)<<trackingParticle.momentum().eta()<<" phi "<<std::setw(9)<<trackingParticle.momentum().phi()<<" event "<<trackingParticle.eventId().event()
              <<" bx "<<trackingParticle.eventId().bunchCrossing()<<" eventNot0"<<std::endl;*/

  if (trackingParticle.eventId().bunchCrossing() != 0)
    return false;

  if (abs(trackingParticle.pdgId()) == 13 || abs(trackingParticle.pdgId()) == 1000015) {
    // 1000015 is stau, todo use other selection if needed
    //only muons
  } else
    return false;

  //in the overlap, the propagation of muons with pt less then ~3.2 fails, the actual threshold depends slightly on eta,
  if (trackingParticle.pt() < 2.5)
    return false;

  if (trackingParticle.parentVertex().isNonnull())
    LogTrace("l1tOmtfEventPrint") << "trackingParticleIsMuonInOmtfBx0, pdgId " << std::setw(3)
                                  << trackingParticle.pdgId() << " pt " << std::setw(9) << trackingParticle.pt()
                                  << " eta " << std::setw(9) << trackingParticle.momentum().eta() << " phi "
                                  << std::setw(9) << trackingParticle.momentum().phi() << " event "
                                  << trackingParticle.eventId().event() << " trackId "
                                  << trackingParticle.g4Tracks().at(0).trackId() << " parentVertex Rho "
                                  << trackingParticle.parentVertex()->position().Rho() << " eta "
                                  << trackingParticle.parentVertex()->position().eta() << " phi "
                                  << trackingParticle.parentVertex()->position().phi() << std::endl;
  else
    LogTrace("l1tOmtfEventPrint") << "trackingParticleIsMuonInOmtfBx0, pdgId " << std::setw(3)
                                  << trackingParticle.pdgId() << " pt " << std::setw(9) << trackingParticle.pt()
                                  << " eta " << std::setw(9) << trackingParticle.momentum().eta() << " phi "
                                  << std::setw(9) << trackingParticle.momentum().phi() << " trackId "
                                  << trackingParticle.g4Tracks().at(0).trackId();

  //higher margin for matching must be used than actual OMTF region (i.e. 0.82 - 1.24),
  //otherwise many candidates are marked as ghosts
  //if( (fabs(simTrack.momentum().eta()) >= 0.82 ) && (fabs(trackingParticle.momentum().eta()) <= 1.24) ) {
  if ((fabs(trackingParticle.momentum().eta()) >= 0.72) && (fabs(trackingParticle.momentum().eta()) <= 1.3)) {
  } else
    return false;

  return true;
}

bool trackingParticleIsMuonInOmtfEvent0

(

const TrackingParticle &

trackingParticle

)

Definition at line 165 of file CandidateSimMuonMatcher.cc.

References EncodedEventId::event(), TrackingParticle::eventId(), and trackingParticleIsMuonInOmtfBx0().

                                                                                  {
  if (trackingParticle.eventId().event() != 0)
    return false;

  return trackingParticleIsMuonInOmtfBx0(trackingParticle);
}