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 "boost/dynamic_bitset.hpp"
#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	simTrackIsMuonInBx0 (const SimTrack &simTrack)
bool	simTrackIsMuonInOmtf (const SimTrack &simTrack)
bool	simTrackIsMuonInOmtfBx0 (const SimTrack &simTrack)
bool	trackingParticleIsMuonInBx0 (const TrackingParticle &trackingParticle)
bool	trackingParticleIsMuonInOmtfBx0 (const TrackingParticle &trackingParticle)
bool	trackingParticleIsMuonInOmtfEvent0 (const TrackingParticle &trackingParticle)

Function Documentation

foldPhi()

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

hwGmtPhiToGlobalPhi()

double hwGmtPhiToGlobalPhi

(

int

phi

)

Definition at line 28 of file CandidateSimMuonMatcher.cc.

References M_PI, and PVValHelper::phi.

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

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

normal_pdf()

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

Definition at line 376 of file CandidateSimMuonMatcher.cc.

References a, JetChargeProducer_cfi::exp, visualization-live-secondInstance_cfg::m, alignCSCRings::s, and x.

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

simTrackIsMuonInBx0()

bool simTrackIsMuonInBx0

(

const SimTrack &

simTrack

)

Definition at line 120 of file CandidateSimMuonMatcher.cc.

References funct::abs(), and cscDigiValidation_cfi::simTrack.

Referenced by CandidateSimMuonMatcher::observeEventEnd().

                                                   {
  if (std::abs(simTrack.type()) == 13 ||
      std::abs(simTrack.type()) == 1000015) {  // 1000015 is stau, todo use other selection (e.g. pt>20) if needed
    //only muons
    if (simTrack.eventId().bunchCrossing() == 0)
      return true;
  }
  return false;
}

simTrackIsMuonInOmtf()

bool simTrackIsMuonInOmtf

(

const SimTrack &

simTrack

)

Definition at line 87 of file CandidateSimMuonMatcher.cc.

References funct::abs(), LogTrace, and cscDigiValidation_cfi::simTrack.

Referenced by simTrackIsMuonInOmtfBx0().

                                                    {
  if (std::abs(simTrack.type()) == 13 ||
      std::abs(simTrack.type()) == 1000015) {  // 1000015 is stau, todo use other selection (e.g. pt>20) 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 (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;

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

  return true;
}

simTrackIsMuonInOmtfBx0()

bool simTrackIsMuonInOmtfBx0

(

const SimTrack &

simTrack

)

Definition at line 113 of file CandidateSimMuonMatcher.cc.

References cscDigiValidation_cfi::simTrack, and simTrackIsMuonInOmtf().

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

  return simTrackIsMuonInOmtf(simTrack);
}

trackingParticleIsMuonInBx0()

bool trackingParticleIsMuonInBx0

(

const TrackingParticle &

trackingParticle

)

Definition at line 130 of file CandidateSimMuonMatcher.cc.

References funct::abs(), EncodedEventId::bunchCrossing(), TrackingParticle::eventId(), and TrackingParticle::pdgId().

Referenced by CandidateSimMuonMatcher::observeEventEnd().

                                                                           {
  if (std::abs(trackingParticle.pdgId()) == 13 ||
      std::abs(trackingParticle.pdgId()) ==
          1000015) {  // 1000015 is stau, todo use other selection (e.g. pt>20) if needed
    //only muons
    if (trackingParticle.eventId().bunchCrossing() == 0)
      return true;
  }
  return false;
}

trackingParticleIsMuonInOmtfBx0()

bool trackingParticleIsMuonInOmtfBx0

(

const TrackingParticle &

trackingParticle

)

Definition at line 141 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 trackingParticleIsMuonInOmtfEvent0().

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

  if (std::abs(trackingParticle.pdgId()) == 13 || std::abs(trackingParticle.pdgId()) == 1000015) {
    // 1000015 is stau, todo use other selection (e.g. pt>20) 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();

  //some margin for matching must be used on top of actual OMTF region,
  //i.e. (0.82-1.24)=>(0.72-1.3),
  //otherwise many candidates are marked as ghosts
  if ((std::abs(trackingParticle.momentum().eta()) >= 0.72) && (std::abs(trackingParticle.momentum().eta()) <= 1.3)) {
  } else
    return false;

  return true;
}

trackingParticleIsMuonInOmtfEvent0()

bool trackingParticleIsMuonInOmtfEvent0

(

const TrackingParticle &

trackingParticle

)

Definition at line 182 of file CandidateSimMuonMatcher.cc.

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

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

  return trackingParticleIsMuonInOmtfBx0(trackingParticle);
}