L1MuonRecoTreeProducer Class Reference

Inheritance diagram for L1MuonRecoTreeProducer:

Public Member Functions
TrajectoryStateOnSurface	cylExtrapTrkSam (reco::TrackRef track, double rho)
void	empty_global ()
void	empty_hlt ()
void	empty_standalone ()
void	empty_tracker ()
	L1MuonRecoTreeProducer (const edm::ParameterSet &)
double	match_trigger (std::vector< int > &trigIndices, const trigger::TriggerObjectCollection &trigObjs, edm::Handle< trigger::TriggerEvent > &triggerEvent, const reco::Muon &mu)
TrajectoryStateOnSurface	surfExtrapTrkSam (reco::TrackRef track, double z)
	~L1MuonRecoTreeProducer () override
Public Member Functions inherited from edm::EDAnalyzer
void	callWhenNewProductsRegistered (std::function< void(BranchDescription const &)> const &func)
	EDAnalyzer ()
SerialTaskQueue *	globalLuminosityBlocksQueue ()
SerialTaskQueue *	globalRunsQueue ()
ModuleDescription const &	moduleDescription () const
std::string	workerType () const
	~EDAnalyzer () override
Public Member Functions inherited from edm::EDConsumerBase
std::vector< ConsumesInfo >	consumesInfo () const
void	convertCurrentProcessAlias (std::string const &processName)
	Convert "@currentProcess" in InputTag process names to the actual current process name. More...
	EDConsumerBase ()
	EDConsumerBase (EDConsumerBase const &)=delete
	EDConsumerBase (EDConsumerBase &&)=default
ESProxyIndex const *	esGetTokenIndices (edm::Transition iTrans) const
std::vector< ESProxyIndex > const &	esGetTokenIndicesVector (edm::Transition iTrans) const
std::vector< ESRecordIndex > const &	esGetTokenRecordIndicesVector (edm::Transition iTrans) const
ProductResolverIndexAndSkipBit	indexFrom (EDGetToken, BranchType, TypeID const &) const
void	itemsMayGet (BranchType, std::vector< ProductResolverIndexAndSkipBit > &) const
void	itemsToGet (BranchType, std::vector< ProductResolverIndexAndSkipBit > &) const
std::vector < ProductResolverIndexAndSkipBit > const &	itemsToGetFrom (BranchType iType) const
void	labelsForToken (EDGetToken iToken, Labels &oLabels) const
void	modulesWhoseProductsAreConsumed (std::array< std::vector< ModuleDescription const * > *, NumBranchTypes > &modulesAll, std::vector< ModuleProcessName > &modulesInPreviousProcesses, ProductRegistry const &preg, std::map< std::string, ModuleDescription const * > const &labelsToDesc, std::string const &processName) const
EDConsumerBase const &	operator= (EDConsumerBase const &)=delete
EDConsumerBase &	operator= (EDConsumerBase &&)=default
bool	registeredToConsume (ProductResolverIndex, bool, BranchType) const
bool	registeredToConsumeMany (TypeID const &, BranchType) const
void	selectInputProcessBlocks (ProductRegistry const &productRegistry, ProcessBlockHelperBase const &processBlockHelperBase)
ProductResolverIndexAndSkipBit	uncheckedIndexFrom (EDGetToken) const
void	updateLookup (BranchType iBranchType, ProductResolverIndexHelper const &, bool iPrefetchMayGet)
void	updateLookup (eventsetup::ESRecordsToProxyIndices const &)
virtual	~EDConsumerBase () noexcept(false)

Public Attributes
L1Analysis::L1AnalysisRecoMuon *	muon
L1Analysis::L1AnalysisRecoMuonDataFormat *	muonData
L1Analysis::L1AnalysisRecoRpcHit *	rpcHit
L1Analysis::L1AnalysisRecoRpcHitDataFormat *	rpcHitData

Private Types
enum	{ GL_MUON = 0, SA_MUON = 1, TR_MUON = 2, TRSA_MUON = 3 }

Private Member Functions
void	analyze (const edm::Event &, const edm::EventSetup &) override
void	beginJob (void) override
void	beginRun (const edm::Run &, const edm::EventSetup &) override
void	endJob () override
void	endRun (const edm::Run &, const edm::EventSetup &) override
FreeTrajectoryState	freeTrajStateMuon (reco::TrackRef track)

Private Attributes
edm::ESHandle< CSCGeometry >	cscGeom
edm::ESHandle< DTGeometry >	dtGeom
edm::Service< TFileService >	fs_
HLTConfigProvider	hltConfig_
std::vector< int >	isoTriggerIndices_
std::vector< std::string >	isoTriggerNames_
unsigned	maxMuon_
unsigned	maxRpcHit_
edm::InputTag	muonTag_
edm::ESHandle< Propagator >	propagatorAlong
edm::ESHandle< Propagator >	propagatorOpposite
edm::ESHandle< RPCGeometry >	rpcGeom
edm::InputTag	rpcHitTag_
bool	runOnPostLS1_
edm::ESHandle< MagneticField >	theBField
edm::ESHandle < GlobalTrackingGeometry >	theTrackingGeometry
TTree *	tree_
std::vector< int >	triggerIndices_
bool	triggerMatching_
double	triggerMaxDeltaR_
std::vector< std::string >	triggerNames_
std::string	triggerProcessLabel_
edm::InputTag	triggerSummaryLabel_

Additional Inherited Members
Public Types inherited from edm::EDAnalyzer
typedef EDAnalyzer	ModuleType
Public Types inherited from edm::EDConsumerBase
typedef ProductLabels	Labels
Static Public Member Functions inherited from edm::EDAnalyzer
static const std::string &	baseType ()
static void	fillDescriptions (ConfigurationDescriptions &descriptions)
static void	prevalidate (ConfigurationDescriptions &)
static bool	wantsGlobalLuminosityBlocks ()
static bool	wantsGlobalRuns ()
static bool	wantsInputProcessBlocks ()
static bool	wantsProcessBlocks ()
static bool	wantsStreamLuminosityBlocks ()
static bool	wantsStreamRuns ()
Protected Member Functions inherited from edm::EDConsumerBase
template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType >	consumes (edm::InputTag const &tag)
template<BranchType B = InEvent>
EDConsumerBaseAdaptor< B >	consumes (edm::InputTag tag) noexcept
EDGetToken	consumes (const TypeToGet &id, edm::InputTag const &tag)
template<BranchType B>
EDGetToken	consumes (TypeToGet const &id, edm::InputTag const &tag)
ConsumesCollector	consumesCollector ()
	Use a ConsumesCollector to gather consumes information from helper functions. More...
template<typename ProductType , BranchType B = InEvent>
void	consumesMany ()
void	consumesMany (const TypeToGet &id)
template<BranchType B>
void	consumesMany (const TypeToGet &id)
template<typename ESProduct , typename ESRecord , Transition Tr = Transition::Event>
auto	esConsumes ()
template<typename ESProduct , typename ESRecord , Transition Tr = Transition::Event>
auto	esConsumes (ESInputTag const &tag)
template<Transition Tr = Transition::Event>
constexpr auto	esConsumes () noexcept
template<Transition Tr = Transition::Event>
auto	esConsumes (ESInputTag tag) noexcept
template<Transition Tr = Transition::Event>
ESGetTokenGeneric	esConsumes (eventsetup::EventSetupRecordKey const &iRecord, eventsetup::DataKey const &iKey)
	Used with EventSetupRecord::doGet. More...
template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType >	mayConsume (edm::InputTag const &tag)
EDGetToken	mayConsume (const TypeToGet &id, edm::InputTag const &tag)
template<BranchType B>
EDGetToken	mayConsume (const TypeToGet &id, edm::InputTag const &tag)
void	resetItemsToGetFrom (BranchType iType)

Detailed Description

Definition at line 102 of file L1MuonRecoTreeProducer.cc.

Member Enumeration Documentation

anonymous enum

private

Enumerator
GL_MUON
SA_MUON
TR_MUON
TRSA_MUON

Definition at line 151 of file L1MuonRecoTreeProducer.cc.

{ GL_MUON = 0, SA_MUON = 1, TR_MUON = 2, TRSA_MUON = 3 };

Constructor & Destructor Documentation

L1MuonRecoTreeProducer::L1MuonRecoTreeProducer ( const edm::ParameterSet &  iConfig )

explicit

Definition at line 188 of file L1MuonRecoTreeProducer.cc.

References fs_, L1Analysis::L1AnalysisRecoRpcHit::getData(), edm::ParameterSet::getParameter(), edm::ParameterSet::getUntrackedParameter(), isoTriggerNames_, TFileService::make(), maxMuon_, maxRpcHit_, muon, muonData, muonTag_, rpcHit, rpcHitData, rpcHitTag_, runOnPostLS1_, AlCaHLTBitMon_QueryRunRegistry::string, tree_, triggerMatching_, triggerMaxDeltaR_, triggerNames_, triggerProcessLabel_, and triggerSummaryLabel_.

                                                                             {
  maxMuon_ = iConfig.getParameter<unsigned int>("maxMuon");
  maxRpcHit_ = iConfig.getParameter<unsigned int>("maxMuon");

  muonTag_ = iConfig.getParameter<edm::InputTag>("muonTag");
  rpcHitTag_ = iConfig.getParameter<edm::InputTag>("rpcHitTag");

  runOnPostLS1_ = iConfig.getParameter<bool>("runOnPostLS1");

  muon = new L1Analysis::L1AnalysisRecoMuon();
  muonData = muon->getData();

  rpcHit = new L1Analysis::L1AnalysisRecoRpcHit();
  rpcHitData = rpcHit->getData();

  // set up output
  tree_ = fs_->make<TTree>("MuonRecoTree", "MuonRecoTree");
  tree_->Branch("Muon", "L1Analysis::L1AnalysisRecoMuonDataFormat", &muonData, 32000, 3);
  tree_->Branch("RpcHit", "L1Analysis::L1AnalysisRecoRpcHitDataFormat", &rpcHitData, 32000, 3);

  //---------------------------------------------------------------------------
  // TRIGGER MATCHING
  // member variables needed for matching, if using reco muons instead of PAT
  //---------------------------------------------------------------------------
  triggerMatching_ = iConfig.getUntrackedParameter<bool>("triggerMatching");
  triggerSummaryLabel_ = iConfig.getParameter<edm::InputTag>("triggerSummaryLabel");
  triggerProcessLabel_ = iConfig.getUntrackedParameter<std::string>("triggerProcessLabel");
  triggerNames_ = iConfig.getParameter<std::vector<std::string> >("triggerNames");
  isoTriggerNames_ = iConfig.getParameter<std::vector<std::string> >("isoTriggerNames");
  triggerMaxDeltaR_ = iConfig.getParameter<double>("triggerMaxDeltaR");
}

L1MuonRecoTreeProducer::~L1MuonRecoTreeProducer ( )

override

Definition at line 220 of file L1MuonRecoTreeProducer.cc.

{}

Member Function Documentation

void L1MuonRecoTreeProducer::analyze ( const edm::Event &  iEvent, const edm::EventSetup &  iSetup  )

overrideprivatevirtual

Implements edm::EDAnalyzer.

Definition at line 398 of file L1MuonRecoTreeProducer.cc.

References beam_dqm_sourceclient-live_cfg::beamSpot, makePileupJSON::bx, L1Analysis::L1AnalysisRecoRpcHitDataFormat::bx, L1Analysis::L1AnalysisRecoMuonDataFormat::calo_energy, L1Analysis::L1AnalysisRecoMuonDataFormat::calo_energy3x3, L1Analysis::L1AnalysisRecoMuonDataFormat::ch, L1Analysis::L1AnalysisRecoRpcHitDataFormat::clusterSize, MuonSubdetId::CSC, cscGeom, cylExtrapTrkSam(), DetId::det(), reco::MuonTime::direction(), dtGeom, L1Analysis::L1AnalysisRecoMuonDataFormat::ecal_terr, reco::MuonEnergy::ecal_time, L1Analysis::L1AnalysisRecoMuonDataFormat::ecal_time, reco::MuonEnergy::ecal_timeError, empty_global(), empty_hlt(), empty_standalone(), empty_tracker(), PV3DBase< T, PVType, FrameType >::eta(), L1Analysis::L1AnalysisRecoMuonDataFormat::eta, edm::EventSetup::get(), edm::Event::getByLabel(), GL_MUON, TrajectoryStateOnSurface::globalPosition(), runTauDisplay::gp, L1Analysis::L1AnalysisRecoMuonDataFormat::hcal_terr, reco::MuonEnergy::hcal_time, L1Analysis::L1AnalysisRecoMuonDataFormat::hcal_time, reco::MuonEnergy::hcal_timeError, L1Analysis::L1AnalysisRecoMuonDataFormat::hlt_deltaR, L1Analysis::L1AnalysisRecoMuonDataFormat::hlt_isoDeltaR, L1Analysis::L1AnalysisRecoMuonDataFormat::hlt_isomu, L1Analysis::L1AnalysisRecoMuonDataFormat::hlt_mu, L1Analysis::L1AnalysisRecoMuonDataFormat::howmanytypes, gpuClustering::id, L1Analysis::L1AnalysisRecoMuonDataFormat::imp_point_p, L1Analysis::L1AnalysisRecoMuonDataFormat::imp_point_pt, L1Analysis::L1AnalysisRecoMuonDataFormat::imp_point_x, L1Analysis::L1AnalysisRecoMuonDataFormat::imp_point_y, L1Analysis::L1AnalysisRecoMuonDataFormat::imp_point_z, reco::TransientTrack::impactPointTSCP(), isoTriggerIndices_, TrajectoryStateOnSurface::isValid(), edm::HandleBase::isValid(), TrajectoryStateClosestToPoint::isValid(), edm::InputTag::label(), L1Analysis::L1AnalysisRecoRpcHitDataFormat::layer, RPCDetId::layer(), phase1PixelTopology::layer, CSCSegment::localPosition(), log, PV3DBase< T, PVType, FrameType >::mag(), match_trigger(), maxMuon_, maxRpcHit_, DetId::Muon, muon, muonData, L1Analysis::L1AnalysisRecoRpcHitDataFormat::muonId, muonTag_, L1Analysis::L1AnalysisRecoMuonDataFormat::nMuons, L1Analysis::L1AnalysisRecoMuonDataFormat::normchi2, L1Analysis::L1AnalysisRecoRpcHitDataFormat::nRpcHits, L1Analysis::L1AnalysisRecoMuonDataFormat::numberOfMatchedStations, L1Analysis::L1AnalysisRecoMuonDataFormat::numberOfValidMuonHits, L1Analysis::L1AnalysisRecoMuonDataFormat::p, PV3DBase< T, PVType, FrameType >::perp2(), PV3DBase< T, PVType, FrameType >::phi(), L1Analysis::L1AnalysisRecoMuonDataFormat::phi, L1Analysis::L1AnalysisRecoMuonDataFormat::phi_hb, L1Analysis::L1AnalysisRecoMuonDataFormat::phi_he_n, L1Analysis::L1AnalysisRecoMuonDataFormat::phi_he_p, L1Analysis::L1AnalysisRecoRpcHitDataFormat::phiGlob, Pi, TrajectoryStateClosestToPoint::position(), propagatorAlong, propagatorOpposite, L1Analysis::L1AnalysisRecoMuonDataFormat::pt, L1Analysis::L1AnalysisRecoMuonDataFormat::r_he_n, L1Analysis::L1AnalysisRecoMuonDataFormat::r_he_p, DetId::rawId(), L1Analysis::L1AnalysisRecoMuonDataFormat::rchCSCtype, L1Analysis::L1AnalysisRecoMuonDataFormat::rchEta, L1Analysis::L1AnalysisRecoMuonDataFormat::rchPhi, L1Analysis::L1AnalysisRecoRpcHitDataFormat::region, RPCDetId::region(), HLT_FULL_cff::region, L1Analysis::L1AnalysisRecoRpcHit::Reset(), L1Analysis::L1AnalysisRecoMuon::Reset(), L1Analysis::L1AnalysisRecoRpcHitDataFormat::ring, RPCDetId::ring(), relativeConstraints::ring, L1Analysis::L1AnalysisRecoRpcHitDataFormat::roll, RPCDetId::roll(), MuonSubdetId::RPC, RPCDetId, rpcGeom, rpcHit, rpcHitData, rpcHitTag_, rpcRecHits_cfi::rpcRecHits, findQualityFiles::rr, runOnPostLS1_, L1Analysis::L1AnalysisRecoMuonDataFormat::sa_ch, L1Analysis::L1AnalysisRecoMuonDataFormat::sa_eta, L1Analysis::L1AnalysisRecoMuonDataFormat::sa_imp_point_p, L1Analysis::L1AnalysisRecoMuonDataFormat::sa_imp_point_pt, L1Analysis::L1AnalysisRecoMuonDataFormat::sa_imp_point_x, L1Analysis::L1AnalysisRecoMuonDataFormat::sa_imp_point_y, L1Analysis::L1AnalysisRecoMuonDataFormat::sa_imp_point_z, L1Analysis::L1AnalysisRecoMuonDataFormat::sa_inner_eta, L1Analysis::L1AnalysisRecoMuonDataFormat::sa_inner_phi, L1Analysis::L1AnalysisRecoMuonDataFormat::sa_inner_pt, L1Analysis::L1AnalysisRecoMuonDataFormat::sa_inner_x, L1Analysis::L1AnalysisRecoMuonDataFormat::sa_inner_y, L1Analysis::L1AnalysisRecoMuonDataFormat::sa_inner_z, SA_MUON, L1Analysis::L1AnalysisRecoMuonDataFormat::sa_nChambers, L1Analysis::L1AnalysisRecoMuonDataFormat::sa_nMatches, L1Analysis::L1AnalysisRecoMuonDataFormat::sa_normchi2, L1Analysis::L1AnalysisRecoMuonDataFormat::sa_outer_eta, L1Analysis::L1AnalysisRecoMuonDataFormat::sa_outer_phi, L1Analysis::L1AnalysisRecoMuonDataFormat::sa_outer_pt, L1Analysis::L1AnalysisRecoMuonDataFormat::sa_outer_x, L1Analysis::L1AnalysisRecoMuonDataFormat::sa_outer_y, L1Analysis::L1AnalysisRecoMuonDataFormat::sa_outer_z, L1Analysis::L1AnalysisRecoMuonDataFormat::sa_p, L1Analysis::L1AnalysisRecoMuonDataFormat::sa_phi, L1Analysis::L1AnalysisRecoMuonDataFormat::sa_phi_hb, L1Analysis::L1AnalysisRecoMuonDataFormat::sa_phi_he_n, L1Analysis::L1AnalysisRecoMuonDataFormat::sa_phi_he_p, L1Analysis::L1AnalysisRecoMuonDataFormat::sa_phi_mb1, L1Analysis::L1AnalysisRecoMuonDataFormat::sa_phi_mb2, L1Analysis::L1AnalysisRecoMuonDataFormat::sa_phi_me1_n, L1Analysis::L1AnalysisRecoMuonDataFormat::sa_phi_me1_p, L1Analysis::L1AnalysisRecoMuonDataFormat::sa_phi_me2_n, L1Analysis::L1AnalysisRecoMuonDataFormat::sa_phi_me2_p, L1Analysis::L1AnalysisRecoMuonDataFormat::sa_pseta, L1Analysis::L1AnalysisRecoMuonDataFormat::sa_pt, L1Analysis::L1AnalysisRecoMuonDataFormat::sa_r_he_n, L1Analysis::L1AnalysisRecoMuonDataFormat::sa_r_he_p, L1Analysis::L1AnalysisRecoMuonDataFormat::sa_r_me1_n, L1Analysis::L1AnalysisRecoMuonDataFormat::sa_r_me1_p, L1Analysis::L1AnalysisRecoMuonDataFormat::sa_r_me2_n, L1Analysis::L1AnalysisRecoMuonDataFormat::sa_r_me2_p, L1Analysis::L1AnalysisRecoMuonDataFormat::sa_validhits, L1Analysis::L1AnalysisRecoMuonDataFormat::sa_z_hb, L1Analysis::L1AnalysisRecoMuonDataFormat::sa_z_mb1, L1Analysis::L1AnalysisRecoMuonDataFormat::sa_z_mb2, L1Analysis::L1AnalysisRecoRpcHitDataFormat::sector, RPCDetId::sector(), mathSSE::sqrt(), edm_modernize_messagelogger::stat, L1Analysis::L1AnalysisRecoRpcHitDataFormat::station, relativeConstraints::station, RPCDetId::station(), L1Analysis::L1AnalysisRecoRpcHitDataFormat::strip, DetId::subdetId(), L1Analysis::L1AnalysisRecoRpcHitDataFormat::subsector, RPCDetId::subsector(), surfExtrapTrkSam(), funct::tan(), theBField, theTrackingGeometry, L1Analysis::L1AnalysisRecoMuonDataFormat::time_dir, L1Analysis::L1AnalysisRecoMuonDataFormat::time_inout, L1Analysis::L1AnalysisRecoMuonDataFormat::time_inout_err, L1Analysis::L1AnalysisRecoMuonDataFormat::time_outin, L1Analysis::L1AnalysisRecoMuonDataFormat::time_outin_err, reco::MuonTime::timeAtIpInOut, reco::MuonTime::timeAtIpInOutErr, reco::MuonTime::timeAtIpOutIn, reco::MuonTime::timeAtIpOutInErr, GeomDet::toGlobal(), reco::MuonEnergy::tower, reco::MuonEnergy::towerS9, L1Analysis::L1AnalysisRecoMuonDataFormat::tr_ch, L1Analysis::L1AnalysisRecoMuonDataFormat::tr_d0, L1Analysis::L1AnalysisRecoMuonDataFormat::tr_eta, L1Analysis::L1AnalysisRecoMuonDataFormat::tr_imp_point_p, L1Analysis::L1AnalysisRecoMuonDataFormat::tr_imp_point_pt, L1Analysis::L1AnalysisRecoMuonDataFormat::tr_imp_point_x, L1Analysis::L1AnalysisRecoMuonDataFormat::tr_imp_point_y, L1Analysis::L1AnalysisRecoMuonDataFormat::tr_imp_point_z, TR_MUON, L1Analysis::L1AnalysisRecoMuonDataFormat::tr_normchi2, L1Analysis::L1AnalysisRecoMuonDataFormat::tr_p, L1Analysis::L1AnalysisRecoMuonDataFormat::tr_phi, L1Analysis::L1AnalysisRecoMuonDataFormat::tr_phi_mb1, L1Analysis::L1AnalysisRecoMuonDataFormat::tr_phi_mb2, L1Analysis::L1AnalysisRecoMuonDataFormat::tr_phi_me1_n, L1Analysis::L1AnalysisRecoMuonDataFormat::tr_phi_me1_p, L1Analysis::L1AnalysisRecoMuonDataFormat::tr_phi_me2_n, L1Analysis::L1AnalysisRecoMuonDataFormat::tr_phi_me2_p, L1Analysis::L1AnalysisRecoMuonDataFormat::tr_pt, L1Analysis::L1AnalysisRecoMuonDataFormat::tr_r_me1_n, L1Analysis::L1AnalysisRecoMuonDataFormat::tr_r_me1_p, L1Analysis::L1AnalysisRecoMuonDataFormat::tr_r_me2_n, L1Analysis::L1AnalysisRecoMuonDataFormat::tr_r_me2_p, L1Analysis::L1AnalysisRecoMuonDataFormat::tr_validhits, L1Analysis::L1AnalysisRecoMuonDataFormat::tr_validpixhits, L1Analysis::L1AnalysisRecoMuonDataFormat::tr_z_mb1, L1Analysis::L1AnalysisRecoMuonDataFormat::tr_z_mb2, tree_, triggerIndices_, triggerMatching_, triggerMaxDeltaR_, triggerProcessLabel_, triggerResults, triggerSummaryLabel_, TRSA_MUON, L1Analysis::L1AnalysisRecoMuonDataFormat::type, L1Analysis::L1AnalysisRecoMuonDataFormat::validhits, PV3DBase< T, PVType, FrameType >::x(), L1Analysis::L1AnalysisRecoRpcHitDataFormat::xLoc, PV3DBase< T, PVType, FrameType >::y(), PV3DBase< T, PVType, FrameType >::z(), and L1Analysis::L1AnalysisRecoMuonDataFormat::z_hb.

                                                                                        {
  float pig = TMath::Pi();

  muon->Reset();
  rpcHit->Reset();

  //GP start
  // Get the CSC Geometry
  iSetup.get<MuonGeometryRecord>().get(cscGeom);
  //GP end

  // Get the DT Geometry from the setup
  iSetup.get<MuonGeometryRecord>().get(dtGeom);

  // Get the RPC Geometry from the setup
  iSetup.get<MuonGeometryRecord>().get(rpcGeom);

  //Get the Magnetic field from the setup
  iSetup.get<IdealMagneticFieldRecord>().get(theBField);

  // Get the GlobalTrackingGeometry from the setup
  iSetup.get<GlobalTrackingGeometryRecord>().get(theTrackingGeometry);

  edm::Handle<RPCRecHitCollection> rpcRecHits;

  iEvent.getByLabel(rpcHitTag_, rpcRecHits);
  if (!rpcRecHits.isValid()) {
    edm::LogInfo("L1Prompt") << "can't find RPCRecHitCollection with label " << rpcHitTag_.label();
  } else {
    RPCRecHitCollection::const_iterator recHitIt = rpcRecHits->begin();
    RPCRecHitCollection::const_iterator recHitEnd = rpcRecHits->end();

    int iRpcRecHits = 0;

    for (; recHitIt != recHitEnd; ++recHitIt) {
      if ((unsigned int)iRpcRecHits > maxRpcHit_ - 1)
        continue;

      int cls = recHitIt->clusterSize();
      int firststrip = recHitIt->firstClusterStrip();
      int bx = recHitIt->BunchX();

      RPCDetId rpcId = recHitIt->rpcId();
      int region = rpcId.region();
      int stat = rpcId.station();
      int sect = rpcId.sector();
      int layer = rpcId.layer();
      int subsector = rpcId.subsector();
      int roll = rpcId.roll();
      int ring = rpcId.ring();

      LocalPoint recHitPosLoc = recHitIt->localPosition();
      const BoundPlane &RPCSurface = rpcGeom->roll(rpcId)->surface();
      GlobalPoint recHitPosGlob = RPCSurface.toGlobal(recHitPosLoc);

      float xLoc = recHitPosLoc.x();
      float phiGlob = recHitPosGlob.phi();

      rpcHitData->region.push_back(region);
      rpcHitData->clusterSize.push_back(cls);
      rpcHitData->strip.push_back(firststrip);
      rpcHitData->bx.push_back(bx);
      rpcHitData->xLoc.push_back(xLoc);
      rpcHitData->phiGlob.push_back(phiGlob);
      rpcHitData->station.push_back(stat);
      rpcHitData->sector.push_back(sect);
      rpcHitData->layer.push_back(layer);
      rpcHitData->subsector.push_back(subsector);
      rpcHitData->roll.push_back(roll);
      rpcHitData->ring.push_back(ring);
      rpcHitData->muonId.push_back(-999);  // CB set to invalid now, updated when looking at muon info

      iRpcRecHits++;
    }

    rpcHitData->nRpcHits = iRpcRecHits;
  }

  // Get the muon candidates
  edm::Handle<reco::MuonCollection> mucand;
  iEvent.getByLabel(muonTag_, mucand);
  if (!mucand.isValid()) {
    edm::LogInfo("L1Prompt") << "can't find Muon Collection with label " << muonTag_.label();
    return;
  }

  // Get the beamspot
  edm::Handle<reco::BeamSpot> beamSpot;
  iEvent.getByLabel("offlineBeamSpot", beamSpot);

  // Get the primary vertices
  edm::Handle<std::vector<reco::Vertex> > vertex;
  iEvent.getByLabel(edm::InputTag("offlinePrimaryVertices"), vertex);

  // Get the propagators
  iSetup.get<TrackingComponentsRecord>().get("SmartPropagatorAny", propagatorAlong);
  iSetup.get<TrackingComponentsRecord>().get("SmartPropagatorAnyOpposite", propagatorOpposite);

  for (reco::MuonCollection::const_iterator imu = mucand->begin();
       // for(pat::MuonCollection::const_iterator imu = mucand->begin();
       imu != mucand->end() && (unsigned)muonData->nMuons < maxMuon_;
       imu++) {
    int type = 0;
    if (imu->isGlobalMuon())
      type = type + 1;
    if (imu->isStandAloneMuon())
      type = type + 2;
    if (imu->isTrackerMuon())
      type = type + 4;
    if (imu->isCaloMuon())
      type = type + 8;

    bool isTIGHT =
        (!vertex->empty() && imu->isGlobalMuon() && imu->globalTrack()->normalizedChi2() < 10. &&
         imu->globalTrack()->hitPattern().numberOfValidMuonHits() > 0 && imu->numberOfMatchedStations() > 1 &&
         fabs(imu->innerTrack()->dxy(vertex->at(0).position())) < 0.2 &&
         fabs(imu->innerTrack()->dz(vertex->at(0).position())) < 0.5 &&
         imu->innerTrack()->hitPattern().numberOfValidPixelHits() > 0 &&
         imu->innerTrack()->hitPattern().trackerLayersWithMeasurement() > 5);

    if (isTIGHT)
      type = type + 16;
    if (imu->isPFMuon())
      type = type + 32;

    muonData->howmanytypes.push_back(type);  // muon type is counting calo muons and multiple assignments

    // bool type identifiers are exclusive. is this correct? CB to check
    bool isSA = (!imu->isGlobalMuon() && imu->isStandAloneMuon() && !imu->isTrackerMuon());
    bool isTR = (!imu->isGlobalMuon() && imu->isTrackerMuon() && !imu->isStandAloneMuon());
    bool isGL = (imu->isGlobalMuon());  //&&!(imu->isStandAloneMuon())&&!(imu->isTrackerMuon()));
    bool isTRSA = (!imu->isGlobalMuon() && imu->isStandAloneMuon() && imu->isTrackerMuon());

    // How we fill this. We have 3 blocks of variables: muons_; muons_sa_; muons_tr_
    //   GL   SA   TR      Description               muon_type     Bool id   Filled Vars               Variables to -99999
    //   0    0    0       Muon does not exist         /             /         /                             /
    //   0    0    1       It is a Tracker muon      TR_MUON       isTR      muons_tr_                  muons_sa_,muons_
    //   0    1    0       It is a SA muon           SA_MUON       isSA      muons_sa                   muons_tr,muons_
    //   0    1    1       It is a SA+Tracker muon   TRSA_MUON     isTRSA    muons_sa,muons_tr          muons_
    //   1    0    0       It is a Global only muon  GL_MUON       isGL      muons_,muons_sa,muons_tr        /
    //   1    0    1       Gl w/out SA cannot exist    /            /           /                            /
    //   1    1    0       Gl+SA (no trk-mu match)   GL_MUON       isGL      muons_,muons_sa,muons_tr     none
    //   1    1    1       GL+SA+Tr (all matched)    GL_MUON       isGL      muons_,muons_sa,muons_tr     none

    //---------------------------------------------------------------------
    // TRIGGER MATCHING:
    // if specified the reconstructed muons are matched to a trigger
    //---------------------------------------------------------------------
    if (triggerMatching_) {
      double isoMatchDeltaR = 9999.;
      double matchDeltaR = 9999.;
      int hasIsoTriggered = 0;
      int hasTriggered = 0;

      // first check if the trigger results are valid:
      edm::Handle<edm::TriggerResults> triggerResults;
      iEvent.getByLabel(edm::InputTag("TriggerResults", "", triggerProcessLabel_), triggerResults);

      if (triggerResults.isValid()) {
        edm::Handle<trigger::TriggerEvent> triggerEvent;
        iEvent.getByLabel(triggerSummaryLabel_, triggerEvent);
        if (triggerEvent.isValid()) {
          // get trigger objects:
          const trigger::TriggerObjectCollection triggerObjects = triggerEvent->getObjects();

          matchDeltaR = match_trigger(triggerIndices_, triggerObjects, triggerEvent, (*imu));
          if (matchDeltaR < triggerMaxDeltaR_)
            hasTriggered = 1;

          isoMatchDeltaR = match_trigger(isoTriggerIndices_, triggerObjects, triggerEvent, (*imu));
          if (isoMatchDeltaR < triggerMaxDeltaR_)
            hasIsoTriggered = 1;
        }  // end if (triggerEvent.isValid())
      }    // end if (triggerResults.isValid())

      // fill trigger matching variables:
      muonData->hlt_isomu.push_back(hasIsoTriggered);
      muonData->hlt_mu.push_back(hasTriggered);
      muonData->hlt_isoDeltaR.push_back(isoMatchDeltaR);
      muonData->hlt_deltaR.push_back(matchDeltaR);
    } else {
      empty_hlt();
    }  // end if (triggerMatching_)

    if (isGL || isTR || isSA || isTRSA) {
      muonData->nMuons = muonData->nMuons + 1;
      if (isTR)
        muonData->type.push_back(TR_MUON);
      if (isGL)
        muonData->type.push_back(GL_MUON);
      if (isSA)
        muonData->type.push_back(SA_MUON);
      if (isTRSA)
        muonData->type.push_back(TRSA_MUON);
      if (!isGL)
        empty_global();
      if (!isTR && !isGL && !isTRSA)
        empty_tracker();
      if (!isSA && !isGL && !isTRSA)
        empty_standalone();

      // begin GP
      //---------------------------------------------------------------------
      // RECHIT information in CSC: only for standalone/global muons!
      //---------------------------------------------------------------------
      // An artificial rank to sort RecHits
      // the closer RecHit to key station/layer -> the smaller the rank
      // KK's original idea
      const int lutCSC[4][6] = {
          {26, 24, 22, 21, 23, 25}, {6, 4, 2, 1, 3, 5}, {16, 14, 12, 11, 13, 15}, {36, 34, 32, 31, 33, 35}};

      int globalTypeRCH = -999999;
      // float  localEtaRCH = -999999;
      // float  localPhiRCH = -999999;
      float globalEtaRCH = -999999;
      float globalPhiRCH = -999999;

      if (isSA || isGL) {
        trackingRecHit_iterator hit = imu->outerTrack()->recHitsBegin();
        trackingRecHit_iterator hitEnd = imu->outerTrack()->recHitsEnd();

        for (; hit != hitEnd; ++hit) {
          if (!((*hit)->isValid()))
            continue;

          // Hardware ID of the RecHit (in terms of wire/strips/chambers)
          DetId detid = (*hit)->geographicalId();

          // Interested in muon systems only
          // Look only at CSC Hits (CSC id is 2)
          if (detid.det() != DetId::Muon)
            continue;
          if (detid.subdetId() != MuonSubdetId::CSC)
            continue;

          CSCDetId id(detid.rawId());
          //std::cout << "before Lut id.station = " <<id.station()
          //          << " id.layer() = " << id.layer() << " globalTypeRCH = "
          //          << globalTypeRCH  << std::endl;
          // another sanity check
          if (id.station() < 1)
            continue;

          // Look up some stuff specific to CSCRecHit2D
          // std::cout << " typeid().name() " << typeid(**hit).name() << std::endl;
          // const CSCRecHit2D* CSChit =dynamic_cast<const CSCRecHit2D*>(&**hit);

          const CSCSegment *cscSegment = dynamic_cast<const CSCSegment *>(&**hit);
          //std::cout << "cscSegment = " << cscSegment << std::endl;
          if (cscSegment == nullptr)
            continue;
          // const CSCRecHit2D* CSChit =(CSCRecHit2D*)(&**hit);

          // std::cout << " after CSCRecHit2D, CSChit = "  << CSChit << std::endl;
          // LocalPoint rhitlocal = CSChit->localPosition();
          LocalPoint rhitlocal = cscSegment->localPosition();

          // for debugging purpouses
          //if (printLevel > 0) {
          //std::cout << "!!!!rhitlocal.phi = "<<rhitlocal.phi() << std::endl;
          //std::cout << "rhitlocal.y="<<rhitlocal.y();
          //std::cout << "rhitlocal.z="<<rhitlocal.z();
          //}

          GlobalPoint gp = GlobalPoint(0.0, 0.0, 0.0);

          const CSCChamber *cscchamber = cscGeom->chamber(id);

          if (!cscchamber)
            continue;

          gp = cscchamber->toGlobal(rhitlocal);

          // identify the rechit position
          //int pos = ( ((id.station()-1)*6) + (id.layer()-1) ) + (MAX_CSC_RECHIT*whichMuon);

          // --------------------------------------------------
          // this part has to be deprecated once we are sure of
          // the TMatrixF usage ;)
          // fill the rechits array
          //rchEtaList[pos]      = gp.eta();
          //rchPhiList[pos]      = gp.phi();
          //float phi02PI = gp.phi();
          //if (gp.phi() < 0) phi02PI += (2*PI);

          //rchPhiList_02PI[pos] = phi02PI;
          // --------------------------------------------------

          // --------------------------------------------------
          // See if this hit is closer to the "key" position
          //if( lutCSC[id.station()-1][id.layer()-1]<globalTypeRCH || globalTypeRCH<0 ){
          if (lutCSC[id.station() - 1][3 - 1] < globalTypeRCH || globalTypeRCH < 0) {
            //globalTypeRCH    = lutCSC[id.station()-1][id.layer()-1];
            globalTypeRCH = lutCSC[id.station() - 1][3 - 1];
            // localEtaRCH      = rhitlocal.eta();
            // localPhiRCH      = rhitlocal.phi();
            globalEtaRCH = gp.eta();
            globalPhiRCH = gp.phi();  // phi from -pi to pi
            //std::cout << "globalEtaRCH =  "  <<globalEtaRCH
            //          << " globalPhiRCH = " << globalPhiRCH << std::endl;
            if (globalPhiRCH < 0)
              globalPhiRCH = globalPhiRCH + 2 * pig;  // convert to [0; 2pi]
          }
          // --------------------------------------------------
        }

        hit = imu->outerTrack()->recHitsBegin();

        for (; hit != hitEnd; hit++) {
          if (!((*hit)->isValid()))
            continue;

          DetId detId = (*hit)->geographicalId();

          if (detId.det() != DetId::Muon)
            continue;
          if (detId.subdetId() != MuonSubdetId::RPC)
            continue;

          RPCDetId rpcId = (RPCDetId)(*hit)->geographicalId();

          int region = rpcId.region();
          int stat = rpcId.station();
          int sect = rpcId.sector();
          int layer = rpcId.layer();
          int subsector = rpcId.subsector();
          int roll = rpcId.roll();
          int ring = rpcId.ring();

          float xLoc = (*hit)->localPosition().x();

          for (int iRpcHit = 0; iRpcHit < rpcHitData->nRpcHits; ++iRpcHit) {
            if (region == rpcHitData->region.at(iRpcHit) && stat == rpcHitData->station.at(iRpcHit) &&
                sect == rpcHitData->sector.at(iRpcHit) && layer == rpcHitData->layer.at(iRpcHit) &&
                subsector == rpcHitData->subsector.at(iRpcHit) && roll == rpcHitData->roll.at(iRpcHit) &&
                ring == rpcHitData->ring.at(iRpcHit) && fabs(xLoc - rpcHitData->xLoc.at(iRpcHit)) < 0.01)

              rpcHitData->muonId.at(iRpcHit) = muonData->nMuons;  // CB rpc hit belongs to mu imu
                  // due to cleaning as in 2012 1 rpc hit belogns to 1 mu
          }
        }
      }
      // at the end of the loop, write only the best rechit
      // if(globalPhiRCH > -100.) std::cout << "globalPhiRCH = " << globalPhiRCH
      // << "globalEtaRCH = " << globalEtaRCH << std::endl;
      muonData->rchCSCtype.push_back(globalTypeRCH);
      muonData->rchPhi.push_back(globalPhiRCH);
      muonData->rchEta.push_back(globalEtaRCH);
      // end GP

      if (isGL) {
        // Filling calo energies and calo times
        if (imu->isEnergyValid()) {
          reco::MuonEnergy muon_energy;
          muon_energy = imu->calEnergy();
          muonData->calo_energy.push_back(muon_energy.tower);
          muonData->calo_energy3x3.push_back(muon_energy.towerS9);
          muonData->ecal_time.push_back(muon_energy.ecal_time);
          muonData->ecal_terr.push_back(muon_energy.ecal_timeError);
          muonData->hcal_time.push_back(muon_energy.hcal_time);
          muonData->hcal_terr.push_back(muon_energy.hcal_timeError);
        } else {
          muonData->calo_energy.push_back(-999999);
          muonData->calo_energy3x3.push_back(-999999);
          muonData->ecal_time.push_back(-999999);
          muonData->ecal_terr.push_back(-999999);
          muonData->hcal_time.push_back(-999999);
          muonData->hcal_terr.push_back(-999999);
        }

        // Filling muon time data
        if (imu->isTimeValid()) {
          reco::MuonTime muon_time;
          muon_time = imu->time();
          muonData->time_dir.push_back(muon_time.direction());
          muonData->time_inout.push_back(muon_time.timeAtIpInOut);
          muonData->time_inout_err.push_back(muon_time.timeAtIpInOutErr);
          muonData->time_outin.push_back(muon_time.timeAtIpOutIn);
          muonData->time_outin_err.push_back(muon_time.timeAtIpOutInErr);
        } else {
          muonData->time_dir.push_back(-999999);
          muonData->time_inout.push_back(-999999);
          muonData->time_inout_err.push_back(-999999);
          muonData->time_outin.push_back(-999999);
          muonData->time_outin_err.push_back(-999999);
        }

        // Use the track now, and make a transient track out of it (for extrapolations)
        reco::TrackRef glb_mu = imu->globalTrack();
        reco::TransientTrack ttrack(*glb_mu, &*theBField, theTrackingGeometry);

        // Track quantities
        muonData->ch.push_back(glb_mu->charge());
        muonData->pt.push_back(glb_mu->pt());
        muonData->p.push_back(glb_mu->p());
        muonData->eta.push_back(glb_mu->eta());
        muonData->phi.push_back(glb_mu->phi());
        muonData->normchi2.push_back(glb_mu->normalizedChi2());
        muonData->validhits.push_back(glb_mu->numberOfValidHits());
        muonData->numberOfMatchedStations.push_back(imu->numberOfMatchedStations());
        muonData->numberOfValidMuonHits.push_back(glb_mu->hitPattern().numberOfValidMuonHits());

        // Extrapolation to IP
        if (ttrack.impactPointTSCP().isValid()) {
          muonData->imp_point_x.push_back(ttrack.impactPointTSCP().position().x());
          muonData->imp_point_y.push_back(ttrack.impactPointTSCP().position().y());
          muonData->imp_point_z.push_back(ttrack.impactPointTSCP().position().z());
          muonData->imp_point_p.push_back(sqrt(ttrack.impactPointTSCP().position().mag()));
          muonData->imp_point_pt.push_back(sqrt(ttrack.impactPointTSCP().position().perp2()));
        } else {
          muonData->imp_point_x.push_back(-999999);
          muonData->imp_point_y.push_back(-999999);
          muonData->imp_point_z.push_back(-999999);
          muonData->imp_point_p.push_back(-999999);
          muonData->imp_point_pt.push_back(-999999);
        }

        if (!runOnPostLS1_) {
          //Extrapolation to HB
          TrajectoryStateOnSurface tsos;
          tsos = cylExtrapTrkSam(glb_mu, 235);
          if (tsos.isValid()) {
            double xx = tsos.globalPosition().x();
            double yy = tsos.globalPosition().y();
            double zz = tsos.globalPosition().z();
            muonData->z_hb.push_back(zz);
            double rr = sqrt(xx * xx + yy * yy);
            double cosphi = xx / rr;
            if (yy >= 0)
              muonData->phi_hb.push_back(acos(cosphi));
            else
              muonData->phi_hb.push_back(2 * pig - acos(cosphi));
          } else {
            muonData->phi_hb.push_back(-999999);
            muonData->z_hb.push_back(-999999);
          }

          //Extrapolation to HE+
          tsos = surfExtrapTrkSam(glb_mu, 479);
          if (tsos.isValid()) {
            double xx = tsos.globalPosition().x();
            double yy = tsos.globalPosition().y();
            double rr = sqrt(xx * xx + yy * yy);
            muonData->r_he_p.push_back(rr);
            double cosphi = xx / rr;
            if (yy >= 0)
              muonData->phi_he_p.push_back(acos(cosphi));
            else
              muonData->phi_he_p.push_back(2 * pig - acos(cosphi));
          } else {
            muonData->r_he_p.push_back(-999999);
            muonData->phi_he_p.push_back(-999999);
          }

          //Extrapolation to HE-
          tsos = surfExtrapTrkSam(glb_mu, -479);
          if (tsos.isValid()) {
            double xx = tsos.globalPosition().x();
            double yy = tsos.globalPosition().y();
            double rr = sqrt(xx * xx + yy * yy);
            muonData->r_he_n.push_back(rr);
            double cosphi = xx / rr;
            if (yy >= 0)
              muonData->phi_he_n.push_back(acos(cosphi));
            else
              muonData->phi_he_n.push_back(2 * pig - acos(cosphi));
          } else {
            muonData->r_he_n.push_back(-999999);
            muonData->phi_he_n.push_back(-999999);
          }
        } else {
          muonData->phi_hb.push_back(-999999);
          muonData->z_hb.push_back(-999999);
          muonData->r_he_p.push_back(-999999);
          muonData->phi_he_p.push_back(-999999);
          muonData->r_he_n.push_back(-999999);
          muonData->phi_he_n.push_back(-999999);
        }
      }  // end of IF IS GLOBAL

      // If global muon or tracker muon, fill the tracker track quantities
      if (isTR || isGL || isTRSA) {
        // Take the tracker track and build a transient track out of it
        reco::TrackRef tr_mu = imu->innerTrack();
        reco::TransientTrack ttrack(*tr_mu, &*theBField, theTrackingGeometry);
        // Fill track quantities
        muonData->tr_ch.push_back(tr_mu->charge());
        muonData->tr_pt.push_back(tr_mu->pt());
        muonData->tr_p.push_back(tr_mu->p());
        muonData->tr_eta.push_back(tr_mu->eta());
        muonData->tr_phi.push_back(tr_mu->phi());
        muonData->tr_normchi2.push_back(tr_mu->normalizedChi2());
        muonData->tr_validhits.push_back(tr_mu->numberOfValidHits());
        muonData->tr_validpixhits.push_back(tr_mu->hitPattern().numberOfValidPixelHits());
        // find d0 from vertex position
        edm::Handle<reco::BeamSpot> beamSpot;
        iEvent.getByLabel("offlineBeamSpot", beamSpot);
        muonData->tr_d0.push_back(tr_mu->dxy(beamSpot->position()));

        // Extrapolation to the IP
        if (ttrack.impactPointTSCP().isValid()) {
          muonData->tr_imp_point_x.push_back(ttrack.impactPointTSCP().position().x());
          muonData->tr_imp_point_y.push_back(ttrack.impactPointTSCP().position().y());
          muonData->tr_imp_point_z.push_back(ttrack.impactPointTSCP().position().z());
          muonData->tr_imp_point_p.push_back(sqrt(ttrack.impactPointTSCP().position().mag()));
          muonData->tr_imp_point_pt.push_back(sqrt(ttrack.impactPointTSCP().position().perp2()));
        } else {
          muonData->tr_imp_point_x.push_back(-999999);
          muonData->tr_imp_point_y.push_back(-999999);
          muonData->tr_imp_point_z.push_back(-999999);
          muonData->tr_imp_point_p.push_back(-999999);
          muonData->tr_imp_point_pt.push_back(-999999);
        }

        if (!runOnPostLS1_) {
          TrajectoryStateOnSurface tsos;
          tsos = cylExtrapTrkSam(tr_mu, 410);  // track at MB1 radius - extrapolation
          if (tsos.isValid()) {
            double xx = tsos.globalPosition().x();
            double yy = tsos.globalPosition().y();
            double zz = tsos.globalPosition().z();
            muonData->tr_z_mb1.push_back(zz);
            double rr = sqrt(xx * xx + yy * yy);
            double cosphi = xx / rr;
            if (yy >= 0)
              muonData->tr_phi_mb1.push_back(acos(cosphi));
            else
              muonData->tr_phi_mb1.push_back(2 * pig - acos(cosphi));
          } else {
            muonData->tr_z_mb1.push_back(-999999);
            muonData->tr_phi_mb1.push_back(-999999);
          }

          tsos = cylExtrapTrkSam(tr_mu, 500);  // track at MB2 radius - extrapolation
          if (tsos.isValid()) {
            double xx = tsos.globalPosition().x();
            double yy = tsos.globalPosition().y();
            double zz = tsos.globalPosition().z();
            muonData->tr_z_mb2.push_back(zz);
            double rr = sqrt(xx * xx + yy * yy);
            double cosphi = xx / rr;
            if (yy >= 0)
              muonData->tr_phi_mb2.push_back(acos(cosphi));
            else
              muonData->tr_phi_mb2.push_back(2 * pig - acos(cosphi));
          } else {
            muonData->tr_z_mb2.push_back(-999999);
            muonData->tr_phi_mb2.push_back(-999999);
          }

          tsos = surfExtrapTrkSam(tr_mu, 630);  // track at ME1+ plane - extrapolation
          if (tsos.isValid()) {
            double xx = tsos.globalPosition().x();
            double yy = tsos.globalPosition().y();
            double rr = sqrt(xx * xx + yy * yy);
            muonData->tr_r_me1_p.push_back(rr);
            double cosphi = xx / rr;
            if (yy >= 0)
              muonData->tr_phi_me1_p.push_back(acos(cosphi));
            else
              muonData->tr_phi_me1_p.push_back(2 * pig - acos(cosphi));
          } else {
            muonData->tr_r_me1_p.push_back(-999999);
            muonData->tr_phi_me1_p.push_back(-999999);
          }

          tsos = surfExtrapTrkSam(tr_mu, 790);  // track at ME2+ plane - extrapolation
          if (tsos.isValid()) {
            double xx = tsos.globalPosition().x();
            double yy = tsos.globalPosition().y();
            double rr = sqrt(xx * xx + yy * yy);
            muonData->tr_r_me2_p.push_back(rr);
            double cosphi = xx / rr;
            if (yy >= 0)
              muonData->tr_phi_me2_p.push_back(acos(cosphi));
            else
              muonData->tr_phi_me2_p.push_back(2 * pig - acos(cosphi));
          } else {
            muonData->tr_r_me2_p.push_back(-999999);
            muonData->tr_phi_me2_p.push_back(-999999);
          }

          tsos = surfExtrapTrkSam(tr_mu, -630);  // track at ME1- plane - extrapolation
          if (tsos.isValid()) {
            double xx = tsos.globalPosition().x();
            double yy = tsos.globalPosition().y();
            double rr = sqrt(xx * xx + yy * yy);
            muonData->tr_r_me1_n.push_back(rr);
            double cosphi = xx / rr;
            if (yy >= 0)
              muonData->tr_phi_me1_n.push_back(acos(cosphi));
            else
              muonData->tr_phi_me1_n.push_back(2 * pig - acos(cosphi));
          } else {
            muonData->tr_r_me1_n.push_back(-999999);
            muonData->tr_phi_me1_n.push_back(-999999);
          }

          tsos = surfExtrapTrkSam(tr_mu, -790);  // track at ME2- plane - extrapolation
          if (tsos.isValid()) {
            double xx = tsos.globalPosition().x();
            double yy = tsos.globalPosition().y();
            double rr = sqrt(xx * xx + yy * yy);
            muonData->tr_r_me2_n.push_back(rr);
            double cosphi = xx / rr;
            if (yy >= 0)
              muonData->tr_phi_me2_n.push_back(acos(cosphi));
            else
              muonData->tr_phi_me2_n.push_back(2 * pig - acos(cosphi));
          } else {
            muonData->tr_r_me2_n.push_back(-999999);
            muonData->tr_phi_me2_n.push_back(-999999);
          }
        } else {
          muonData->tr_z_mb1.push_back(-999999);
          muonData->tr_phi_mb1.push_back(-999999);
          muonData->tr_z_mb2.push_back(-999999);
          muonData->tr_phi_mb2.push_back(-999999);
          muonData->tr_r_me1_p.push_back(-999999);
          muonData->tr_phi_me1_p.push_back(-999999);
          muonData->tr_r_me2_p.push_back(-999999);
          muonData->tr_phi_me2_p.push_back(-999999);
          muonData->tr_r_me1_n.push_back(-999999);
          muonData->tr_phi_me1_n.push_back(-999999);
          muonData->tr_r_me2_n.push_back(-999999);
          muonData->tr_phi_me2_n.push_back(-999999);
        }

      }  // end of IF IS TRACKER

      // If global muon or sa muon, fill the sa track quantities
      if (isGL || isSA || isTRSA) {
        muonData->sa_nChambers.push_back(imu->numberOfChambers());
        muonData->sa_nMatches.push_back(imu->numberOfMatches());

        // Take the SA track and build a transient track out of it
        reco::TrackRef sa_mu = imu->outerTrack();
        reco::TransientTrack ttrack(*sa_mu, &*theBField, theTrackingGeometry);

        // Extrapolation to IP
        if (ttrack.impactPointTSCP().isValid()) {
          muonData->sa_imp_point_x.push_back(ttrack.impactPointTSCP().position().x());
          muonData->sa_imp_point_y.push_back(ttrack.impactPointTSCP().position().y());
          muonData->sa_imp_point_z.push_back(ttrack.impactPointTSCP().position().z());
          muonData->sa_imp_point_p.push_back(sqrt(ttrack.impactPointTSCP().position().mag()));
          muonData->sa_imp_point_pt.push_back(sqrt(ttrack.impactPointTSCP().position().perp2()));
        } else {
          muonData->sa_imp_point_x.push_back(-999999);
          muonData->sa_imp_point_y.push_back(-999999);
          muonData->sa_imp_point_z.push_back(-999999);
          muonData->sa_imp_point_p.push_back(-999999);
          muonData->sa_imp_point_pt.push_back(-999999);
        }

        // Extrapolation to MB2

        TrajectoryStateOnSurface tsos;
        tsos = cylExtrapTrkSam(sa_mu, 410);  // track at MB1 radius - extrapolation
        if (tsos.isValid()) {
          double xx = tsos.globalPosition().x();
          double yy = tsos.globalPosition().y();
          double zz = tsos.globalPosition().z();
          muonData->sa_z_mb1.push_back(zz);
          double rr = sqrt(xx * xx + yy * yy);
          double cosphi = xx / rr;
          if (yy >= 0)
            muonData->sa_phi_mb1.push_back(acos(cosphi));
          else
            muonData->sa_phi_mb1.push_back(2 * pig - acos(cosphi));
        } else {
          muonData->sa_z_mb1.push_back(-999999);
          muonData->sa_phi_mb1.push_back(-999999);
        }

        tsos = cylExtrapTrkSam(sa_mu, 500);  // track at MB2 radius - extrapolation
        if (tsos.isValid()) {
          double xx = tsos.globalPosition().x();
          double yy = tsos.globalPosition().y();
          double zz = tsos.globalPosition().z();
          muonData->sa_z_mb2.push_back(zz);
          double rr = sqrt(xx * xx + yy * yy);
          double cosphi = xx / rr;
          if (yy >= 0)
            muonData->sa_phi_mb2.push_back(acos(cosphi));
          else
            muonData->sa_phi_mb2.push_back(2 * pig - acos(cosphi));
          double abspseta = -log(tan(atan(fabs(rr / zz)) / 2.0));
          if (zz >= 0)
            muonData->sa_pseta.push_back(abspseta);
          else
            muonData->sa_pseta.push_back(-abspseta);
        } else {
          muonData->sa_z_mb2.push_back(-999999);
          muonData->sa_phi_mb2.push_back(-999999);
          muonData->sa_pseta.push_back(-999999);
        }

        tsos = surfExtrapTrkSam(sa_mu, 630);  // track at ME1+ plane - extrapolation
        if (tsos.isValid()) {
          double xx = tsos.globalPosition().x();
          double yy = tsos.globalPosition().y();
          double rr = sqrt(xx * xx + yy * yy);
          muonData->sa_r_me1_p.push_back(rr);
          double cosphi = xx / rr;
          if (yy >= 0)
            muonData->sa_phi_me1_p.push_back(acos(cosphi));
          else
            muonData->sa_phi_me1_p.push_back(2 * pig - acos(cosphi));
        } else {
          muonData->sa_r_me1_p.push_back(-999999);
          muonData->sa_phi_me1_p.push_back(-999999);
        }

        // Extrapolation to ME2+
        tsos = surfExtrapTrkSam(sa_mu, 790);
        if (tsos.isValid()) {
          double xx = tsos.globalPosition().x();
          double yy = tsos.globalPosition().y();
          double rr = sqrt(xx * xx + yy * yy);
          muonData->sa_r_me2_p.push_back(rr);
          double cosphi = xx / rr;
          if (yy >= 0)
            muonData->sa_phi_me2_p.push_back(acos(cosphi));
          else
            muonData->sa_phi_me2_p.push_back(2 * pig - acos(cosphi));
        } else {
          muonData->sa_r_me2_p.push_back(-999999);
          muonData->sa_phi_me2_p.push_back(-999999);
        }

        tsos = surfExtrapTrkSam(sa_mu, -630);  // track at ME1- plane - extrapolation
        if (tsos.isValid()) {
          double xx = tsos.globalPosition().x();
          double yy = tsos.globalPosition().y();
          double rr = sqrt(xx * xx + yy * yy);
          muonData->sa_r_me1_n.push_back(rr);
          double cosphi = xx / rr;
          if (yy >= 0)
            muonData->sa_phi_me1_n.push_back(acos(cosphi));
          else
            muonData->sa_phi_me1_n.push_back(2 * pig - acos(cosphi));
        } else {
          muonData->sa_r_me1_n.push_back(-999999);
          muonData->sa_phi_me1_n.push_back(-999999);
        }

        // Extrapolation to ME2-
        tsos = surfExtrapTrkSam(sa_mu, -790);  // track at ME2- disk - extrapolation
        if (tsos.isValid()) {
          double xx = tsos.globalPosition().x();
          double yy = tsos.globalPosition().y();
          double rr = sqrt(xx * xx + yy * yy);
          muonData->sa_r_me2_n.push_back(rr);
          double cosphi = xx / rr;
          if (yy >= 0)
            muonData->sa_phi_me2_n.push_back(acos(cosphi));
          else
            muonData->sa_phi_me2_n.push_back(2 * pig - acos(cosphi));
        } else {
          muonData->sa_r_me2_n.push_back(-999999);
          muonData->sa_phi_me2_n.push_back(-999999);
        }

        // Extrapolation to HB
        tsos = cylExtrapTrkSam(sa_mu, 235);  // track at HB radius - extrapolation
        if (tsos.isValid()) {
          double xx = tsos.globalPosition().x();
          double yy = tsos.globalPosition().y();
          double zz = tsos.globalPosition().z();
          muonData->sa_z_hb.push_back(zz);
          double rr = sqrt(xx * xx + yy * yy);
          double cosphi = xx / rr;
          if (yy >= 0)
            muonData->sa_phi_hb.push_back(acos(cosphi));
          else
            muonData->sa_phi_hb.push_back(2 * pig - acos(cosphi));
        } else {
          muonData->sa_z_hb.push_back(-999999);
          muonData->sa_phi_hb.push_back(-999999);
        }

        // Extrapolation to HE+
        tsos = surfExtrapTrkSam(sa_mu, 479);
        if (tsos.isValid()) {
          double xx = tsos.globalPosition().x();
          double yy = tsos.globalPosition().y();
          double rr = sqrt(xx * xx + yy * yy);
          muonData->sa_r_he_p.push_back(rr);
          double cosphi = xx / rr;
          if (yy >= 0)
            muonData->sa_phi_he_p.push_back(acos(cosphi));
          else
            muonData->sa_phi_he_p.push_back(2 * pig - acos(cosphi));
        } else {
          muonData->sa_r_he_p.push_back(-999999);
          muonData->sa_phi_he_p.push_back(-999999);
        }

        // Extrapolation to HE-
        tsos = surfExtrapTrkSam(sa_mu, -479);
        if (tsos.isValid()) {
          double xx = tsos.globalPosition().x();
          double yy = tsos.globalPosition().y();
          double rr = sqrt(xx * xx + yy * yy);
          muonData->sa_r_he_n.push_back(rr);
          double cosphi = xx / rr;
          if (yy >= 0)
            muonData->sa_phi_he_n.push_back(acos(cosphi));
          else
            muonData->sa_phi_he_n.push_back(2 * pig - acos(cosphi));
        } else {
          muonData->sa_r_he_n.push_back(-999999);
          muonData->sa_phi_he_n.push_back(-999999);
        }

        //  Several track quantities
        muonData->sa_normchi2.push_back(sa_mu->normalizedChi2());
        muonData->sa_validhits.push_back(sa_mu->numberOfValidHits());
        muonData->sa_ch.push_back(sa_mu->charge());
        muonData->sa_pt.push_back(sa_mu->pt());
        muonData->sa_p.push_back(sa_mu->p());
        muonData->sa_eta.push_back(sa_mu->eta());
        muonData->sa_phi.push_back(sa_mu->phi());
        muonData->sa_outer_pt.push_back(sqrt(sa_mu->outerMomentum().Perp2()));
        muonData->sa_inner_pt.push_back(sqrt(sa_mu->innerMomentum().Perp2()));
        muonData->sa_outer_eta.push_back(sa_mu->outerMomentum().Eta());
        muonData->sa_inner_eta.push_back(sa_mu->innerMomentum().Eta());
        muonData->sa_outer_phi.push_back(sa_mu->outerMomentum().Phi());
        muonData->sa_inner_phi.push_back(sa_mu->innerMomentum().Phi());
        muonData->sa_outer_x.push_back(sa_mu->outerPosition().x());
        muonData->sa_outer_y.push_back(sa_mu->outerPosition().y());
        muonData->sa_outer_z.push_back(sa_mu->outerPosition().z());
        muonData->sa_inner_x.push_back(sa_mu->innerPosition().x());
        muonData->sa_inner_y.push_back(sa_mu->innerPosition().y());
        muonData->sa_inner_z.push_back(sa_mu->innerPosition().z());

      }  // end of IF IS STANDALONE

    }  // end of if type 012 found
  }    // end of muon loop

  tree_->Fill();
}

void L1MuonRecoTreeProducer::beginJob ( void  )

overrideprivatevirtual

Reimplemented from edm::EDAnalyzer.

Definition at line 1337 of file L1MuonRecoTreeProducer.cc.

{}

void L1MuonRecoTreeProducer::beginRun ( const edm::Run &  run, const edm::EventSetup &  eventSetup  )

overrideprivatevirtual

Reimplemented from edm::EDAnalyzer.

Definition at line 1282 of file L1MuonRecoTreeProducer.cc.

References cms::cuda::assert(), gather_cfg::cout, hltConfig_, HLTConfigProvider::init(), isoTriggerIndices_, isoTriggerNames_, HLTConfigProvider::size(), triggerIndices_, triggerMatching_, HLTConfigProvider::triggerName(), triggerNames_, and triggerProcessLabel_.

                                                                                        {
  // Prepare for trigger matching for each new run:
  // Look up triggetIndices in the HLT config for the different paths
  if (triggerMatching_) {
    bool changed = true;
    if (!hltConfig_.init(run, eventSetup, triggerProcessLabel_, changed)) {
      // if you can't initialize hlt configuration, crash!
      std::cout << "Error: didn't find process" << triggerProcessLabel_ << std::endl;
      assert(false);
    }

    bool enableWildcard = true;
    for (size_t iTrig = 0; iTrig < triggerNames_.size(); ++iTrig) {
      // prepare for regular expression (with wildcards) functionality:
      TString tNameTmp = TString(triggerNames_[iTrig]);
      TRegexp tNamePattern = TRegexp(tNameTmp, enableWildcard);
      int tIndex = -1;
      // find the trigger index:
      for (unsigned ipath = 0; ipath < hltConfig_.size(); ++ipath) {
        // use TString since it provides reg exp functionality:
        TString tmpName = TString(hltConfig_.triggerName(ipath));
        if (tmpName.Contains(tNamePattern)) {
          tIndex = int(ipath);
          triggerIndices_.push_back(tIndex);
        }
      }
      if (tIndex < 0) {  // if can't find trigger path at all, give warning:
        std::cout << "Warning: Could not find trigger" << triggerNames_[iTrig] << std::endl;
        //assert(false);
      }
    }  // end for triggerNames
    for (size_t iTrig = 0; iTrig < isoTriggerNames_.size(); ++iTrig) {
      // prepare for regular expression functionality:
      TString tNameTmp = TString(isoTriggerNames_[iTrig]);
      TRegexp tNamePattern = TRegexp(tNameTmp, enableWildcard);
      int tIndex = -1;
      // find the trigger index:
      for (unsigned ipath = 0; ipath < hltConfig_.size(); ++ipath) {
        // use TString since it provides reg exp functionality:
        TString tmpName = TString(hltConfig_.triggerName(ipath));
        if (tmpName.Contains(tNamePattern)) {
          tIndex = int(ipath);
          isoTriggerIndices_.push_back(tIndex);
        }
      }
      if (tIndex < 0) {  // if can't find trigger path at all, give warning:
        std::cout << "Warning: Could not find trigger" << isoTriggerNames_[iTrig] << std::endl;
        //assert(false);
      }
    }  // end for isoTriggerNames
  }    // end if (triggerMatching_)
}

TrajectoryStateOnSurface L1MuonRecoTreeProducer::cylExtrapTrkSam	(	reco::TrackRef	track,
		double	rho
	)

Definition at line 1244 of file L1MuonRecoTreeProducer.cc.

References Cylinder::build(), freeTrajStateMuon(), TrajectoryStateOnSurface::isValid(), propagatorAlong, propagatorOpposite, and makeMuonMisalignmentScenario::rot.

Referenced by analyze().

                                                                                               {
  Cylinder::PositionType pos(0, 0, 0);
  Cylinder::RotationType rot;
  Cylinder::CylinderPointer myCylinder = Cylinder::build(pos, rot, rho);

  FreeTrajectoryState recoStart = freeTrajStateMuon(track);
  TrajectoryStateOnSurface recoProp;
  recoProp = propagatorAlong->propagate(recoStart, *myCylinder);
  if (!recoProp.isValid()) {
    recoProp = propagatorOpposite->propagate(recoStart, *myCylinder);
  }
  return recoProp;
}

void L1MuonRecoTreeProducer::empty_global

(

)

Definition at line 263 of file L1MuonRecoTreeProducer.cc.

References L1Analysis::L1AnalysisRecoMuonDataFormat::calo_energy, L1Analysis::L1AnalysisRecoMuonDataFormat::calo_energy3x3, L1Analysis::L1AnalysisRecoMuonDataFormat::ch, L1Analysis::L1AnalysisRecoMuonDataFormat::ecal_terr, L1Analysis::L1AnalysisRecoMuonDataFormat::ecal_time, L1Analysis::L1AnalysisRecoMuonDataFormat::eta, L1Analysis::L1AnalysisRecoMuonDataFormat::hcal_terr, L1Analysis::L1AnalysisRecoMuonDataFormat::hcal_time, L1Analysis::L1AnalysisRecoMuonDataFormat::hlt_deltaR, L1Analysis::L1AnalysisRecoMuonDataFormat::hlt_isoDeltaR, L1Analysis::L1AnalysisRecoMuonDataFormat::hlt_isomu, L1Analysis::L1AnalysisRecoMuonDataFormat::hlt_mu, L1Analysis::L1AnalysisRecoMuonDataFormat::imp_point_p, L1Analysis::L1AnalysisRecoMuonDataFormat::imp_point_pt, L1Analysis::L1AnalysisRecoMuonDataFormat::imp_point_x, L1Analysis::L1AnalysisRecoMuonDataFormat::imp_point_y, L1Analysis::L1AnalysisRecoMuonDataFormat::imp_point_z, muonData, L1Analysis::L1AnalysisRecoMuonDataFormat::normchi2, L1Analysis::L1AnalysisRecoMuonDataFormat::numberOfMatchedStations, L1Analysis::L1AnalysisRecoMuonDataFormat::numberOfValidMuonHits, L1Analysis::L1AnalysisRecoMuonDataFormat::p, L1Analysis::L1AnalysisRecoMuonDataFormat::phi, L1Analysis::L1AnalysisRecoMuonDataFormat::phi_hb, L1Analysis::L1AnalysisRecoMuonDataFormat::phi_he_n, L1Analysis::L1AnalysisRecoMuonDataFormat::phi_he_p, L1Analysis::L1AnalysisRecoMuonDataFormat::pt, L1Analysis::L1AnalysisRecoMuonDataFormat::r_he_n, L1Analysis::L1AnalysisRecoMuonDataFormat::r_he_p, L1Analysis::L1AnalysisRecoMuonDataFormat::time_dir, L1Analysis::L1AnalysisRecoMuonDataFormat::time_inout, L1Analysis::L1AnalysisRecoMuonDataFormat::time_inout_err, L1Analysis::L1AnalysisRecoMuonDataFormat::time_outin, L1Analysis::L1AnalysisRecoMuonDataFormat::time_outin_err, L1Analysis::L1AnalysisRecoMuonDataFormat::validhits, and L1Analysis::L1AnalysisRecoMuonDataFormat::z_hb.

Referenced by analyze().

                                          {
  muonData->ch.push_back(-999999);
  muonData->pt.push_back(-999999);
  muonData->p.push_back(-999999);
  muonData->eta.push_back(-999999);
  muonData->phi.push_back(-999999);
  muonData->normchi2.push_back(-999999);
  muonData->validhits.push_back(-999999);
  muonData->numberOfMatchedStations.push_back(-999999);
  muonData->numberOfValidMuonHits.push_back(-999999);
  muonData->imp_point_x.push_back(-999999);
  muonData->imp_point_y.push_back(-999999);
  muonData->imp_point_z.push_back(-999999);
  muonData->imp_point_p.push_back(-999999);
  muonData->imp_point_pt.push_back(-999999);
  muonData->phi_hb.push_back(-999999);
  muonData->z_hb.push_back(-999999);
  muonData->r_he_p.push_back(-999999);
  muonData->phi_he_p.push_back(-999999);
  muonData->r_he_n.push_back(-999999);
  muonData->phi_he_n.push_back(-999999);
  muonData->calo_energy.push_back(-999999);
  muonData->calo_energy3x3.push_back(-999999);
  muonData->ecal_time.push_back(-999999);
  muonData->ecal_terr.push_back(-999999);
  muonData->hcal_time.push_back(-999999);
  muonData->hcal_terr.push_back(-999999);
  muonData->time_dir.push_back(-999999);
  muonData->time_inout.push_back(-999999);
  muonData->time_inout_err.push_back(-999999);
  muonData->time_outin.push_back(-999999);
  muonData->time_outin_err.push_back(-999999);

  muonData->hlt_isomu.push_back(-999999);
  muonData->hlt_mu.push_back(-999999);
  muonData->hlt_isoDeltaR.push_back(-999999);
  muonData->hlt_deltaR.push_back(-999999);
}

void L1MuonRecoTreeProducer::empty_hlt

(

)

Definition at line 388 of file L1MuonRecoTreeProducer.cc.

References L1Analysis::L1AnalysisRecoMuonDataFormat::hlt_deltaR, L1Analysis::L1AnalysisRecoMuonDataFormat::hlt_isoDeltaR, L1Analysis::L1AnalysisRecoMuonDataFormat::hlt_isomu, L1Analysis::L1AnalysisRecoMuonDataFormat::hlt_mu, and muonData.

Referenced by analyze().

                                       {
  muonData->hlt_isomu.push_back(-999999);
  muonData->hlt_mu.push_back(-999999);
  muonData->hlt_isoDeltaR.push_back(-9999999);
  muonData->hlt_deltaR.push_back(-999999);
}

void L1MuonRecoTreeProducer::empty_standalone

(

)

Definition at line 333 of file L1MuonRecoTreeProducer.cc.

References muonData, L1Analysis::L1AnalysisRecoMuonDataFormat::sa_ch, L1Analysis::L1AnalysisRecoMuonDataFormat::sa_eta, L1Analysis::L1AnalysisRecoMuonDataFormat::sa_imp_point_p, L1Analysis::L1AnalysisRecoMuonDataFormat::sa_imp_point_pt, L1Analysis::L1AnalysisRecoMuonDataFormat::sa_imp_point_x, L1Analysis::L1AnalysisRecoMuonDataFormat::sa_imp_point_y, L1Analysis::L1AnalysisRecoMuonDataFormat::sa_imp_point_z, L1Analysis::L1AnalysisRecoMuonDataFormat::sa_inner_eta, L1Analysis::L1AnalysisRecoMuonDataFormat::sa_inner_phi, L1Analysis::L1AnalysisRecoMuonDataFormat::sa_inner_pt, L1Analysis::L1AnalysisRecoMuonDataFormat::sa_inner_x, L1Analysis::L1AnalysisRecoMuonDataFormat::sa_inner_y, L1Analysis::L1AnalysisRecoMuonDataFormat::sa_inner_z, L1Analysis::L1AnalysisRecoMuonDataFormat::sa_nChambers, L1Analysis::L1AnalysisRecoMuonDataFormat::sa_nMatches, L1Analysis::L1AnalysisRecoMuonDataFormat::sa_normchi2, L1Analysis::L1AnalysisRecoMuonDataFormat::sa_outer_eta, L1Analysis::L1AnalysisRecoMuonDataFormat::sa_outer_phi, L1Analysis::L1AnalysisRecoMuonDataFormat::sa_outer_pt, L1Analysis::L1AnalysisRecoMuonDataFormat::sa_outer_x, L1Analysis::L1AnalysisRecoMuonDataFormat::sa_outer_y, L1Analysis::L1AnalysisRecoMuonDataFormat::sa_outer_z, L1Analysis::L1AnalysisRecoMuonDataFormat::sa_p, L1Analysis::L1AnalysisRecoMuonDataFormat::sa_phi, L1Analysis::L1AnalysisRecoMuonDataFormat::sa_phi_hb, L1Analysis::L1AnalysisRecoMuonDataFormat::sa_phi_he_n, L1Analysis::L1AnalysisRecoMuonDataFormat::sa_phi_he_p, L1Analysis::L1AnalysisRecoMuonDataFormat::sa_phi_mb1, L1Analysis::L1AnalysisRecoMuonDataFormat::sa_phi_mb2, L1Analysis::L1AnalysisRecoMuonDataFormat::sa_phi_me1_n, L1Analysis::L1AnalysisRecoMuonDataFormat::sa_phi_me1_p, L1Analysis::L1AnalysisRecoMuonDataFormat::sa_phi_me2_n, L1Analysis::L1AnalysisRecoMuonDataFormat::sa_phi_me2_p, L1Analysis::L1AnalysisRecoMuonDataFormat::sa_pseta, L1Analysis::L1AnalysisRecoMuonDataFormat::sa_pt, L1Analysis::L1AnalysisRecoMuonDataFormat::sa_r_he_n, L1Analysis::L1AnalysisRecoMuonDataFormat::sa_r_he_p, L1Analysis::L1AnalysisRecoMuonDataFormat::sa_r_me1_n, L1Analysis::L1AnalysisRecoMuonDataFormat::sa_r_me1_p, L1Analysis::L1AnalysisRecoMuonDataFormat::sa_r_me2_n, L1Analysis::L1AnalysisRecoMuonDataFormat::sa_r_me2_p, L1Analysis::L1AnalysisRecoMuonDataFormat::sa_validhits, L1Analysis::L1AnalysisRecoMuonDataFormat::sa_z_hb, L1Analysis::L1AnalysisRecoMuonDataFormat::sa_z_mb1, and L1Analysis::L1AnalysisRecoMuonDataFormat::sa_z_mb2.

Referenced by analyze().

                                              {
  muonData->sa_imp_point_x.push_back(-999999);
  muonData->sa_imp_point_y.push_back(-999999);
  muonData->sa_imp_point_z.push_back(-999999);
  muonData->sa_imp_point_p.push_back(-999999);
  muonData->sa_imp_point_pt.push_back(-999999);
  muonData->sa_z_mb2.push_back(-999999);
  muonData->sa_phi_mb2.push_back(-999999);
  muonData->sa_pseta.push_back(-999999);

  muonData->sa_z_hb.push_back(-999999);
  muonData->sa_phi_hb.push_back(-999999);

  muonData->sa_r_he_p.push_back(-999999);
  muonData->sa_phi_he_p.push_back(-999999);
  muonData->sa_r_he_n.push_back(-999999);
  muonData->sa_phi_he_n.push_back(-999999);

  muonData->sa_normchi2.push_back(-999999);
  muonData->sa_validhits.push_back(-999999);
  muonData->sa_ch.push_back(-999999);
  muonData->sa_pt.push_back(-999999);
  muonData->sa_p.push_back(-999999);
  muonData->sa_eta.push_back(-999999);
  muonData->sa_phi.push_back(-999999);
  muonData->sa_outer_pt.push_back(-999999);
  muonData->sa_inner_pt.push_back(-999999);
  muonData->sa_outer_eta.push_back(-999999);
  muonData->sa_inner_eta.push_back(-999999);
  muonData->sa_outer_phi.push_back(-999999);
  muonData->sa_inner_phi.push_back(-999999);
  muonData->sa_outer_x.push_back(-999999);
  muonData->sa_outer_y.push_back(-999999);
  muonData->sa_outer_z.push_back(-999999);
  muonData->sa_inner_x.push_back(-999999);
  muonData->sa_inner_y.push_back(-999999);
  muonData->sa_inner_z.push_back(-999999);

  muonData->sa_r_me2_p.push_back(-999999);
  muonData->sa_phi_me2_p.push_back(-999999);

  muonData->sa_r_me2_n.push_back(-999999);
  muonData->sa_phi_me2_n.push_back(-999999);

  muonData->sa_z_mb1.push_back(-999999);
  muonData->sa_phi_mb1.push_back(-999999);
  muonData->sa_r_me1_p.push_back(-999999);
  muonData->sa_phi_me1_p.push_back(-999999);
  muonData->sa_r_me1_n.push_back(-999999);
  muonData->sa_phi_me1_n.push_back(-999999);

  muonData->sa_nChambers.push_back(-999);
  muonData->sa_nMatches.push_back(-999);
}

void L1MuonRecoTreeProducer::empty_tracker

(

)

Definition at line 302 of file L1MuonRecoTreeProducer.cc.

References muonData, L1Analysis::L1AnalysisRecoMuonDataFormat::tr_ch, L1Analysis::L1AnalysisRecoMuonDataFormat::tr_d0, L1Analysis::L1AnalysisRecoMuonDataFormat::tr_eta, L1Analysis::L1AnalysisRecoMuonDataFormat::tr_imp_point_p, L1Analysis::L1AnalysisRecoMuonDataFormat::tr_imp_point_pt, L1Analysis::L1AnalysisRecoMuonDataFormat::tr_imp_point_x, L1Analysis::L1AnalysisRecoMuonDataFormat::tr_imp_point_y, L1Analysis::L1AnalysisRecoMuonDataFormat::tr_imp_point_z, L1Analysis::L1AnalysisRecoMuonDataFormat::tr_normchi2, L1Analysis::L1AnalysisRecoMuonDataFormat::tr_p, L1Analysis::L1AnalysisRecoMuonDataFormat::tr_phi, L1Analysis::L1AnalysisRecoMuonDataFormat::tr_phi_mb1, L1Analysis::L1AnalysisRecoMuonDataFormat::tr_phi_mb2, L1Analysis::L1AnalysisRecoMuonDataFormat::tr_phi_me1_n, L1Analysis::L1AnalysisRecoMuonDataFormat::tr_phi_me1_p, L1Analysis::L1AnalysisRecoMuonDataFormat::tr_phi_me2_n, L1Analysis::L1AnalysisRecoMuonDataFormat::tr_phi_me2_p, L1Analysis::L1AnalysisRecoMuonDataFormat::tr_pt, L1Analysis::L1AnalysisRecoMuonDataFormat::tr_r_me1_n, L1Analysis::L1AnalysisRecoMuonDataFormat::tr_r_me1_p, L1Analysis::L1AnalysisRecoMuonDataFormat::tr_r_me2_n, L1Analysis::L1AnalysisRecoMuonDataFormat::tr_r_me2_p, L1Analysis::L1AnalysisRecoMuonDataFormat::tr_validhits, L1Analysis::L1AnalysisRecoMuonDataFormat::tr_validpixhits, L1Analysis::L1AnalysisRecoMuonDataFormat::tr_z_mb1, and L1Analysis::L1AnalysisRecoMuonDataFormat::tr_z_mb2.

Referenced by analyze().

                                           {
  muonData->tr_ch.push_back(-999999);
  muonData->tr_pt.push_back(-999999);
  muonData->tr_p.push_back(-999999);
  muonData->tr_eta.push_back(-999999);
  muonData->tr_phi.push_back(-999999);
  muonData->tr_normchi2.push_back(-999999);
  muonData->tr_validhits.push_back(-999999);
  muonData->tr_validpixhits.push_back(-999999);
  muonData->tr_d0.push_back(-999999);
  muonData->tr_imp_point_x.push_back(-999999);
  muonData->tr_imp_point_y.push_back(-999999);
  muonData->tr_imp_point_z.push_back(-999999);
  muonData->tr_imp_point_p.push_back(-999999);
  muonData->tr_imp_point_pt.push_back(-999999);

  muonData->tr_z_mb2.push_back(-999999);
  muonData->tr_phi_mb2.push_back(-999999);
  muonData->tr_r_me2_p.push_back(-999999);
  muonData->tr_phi_me2_p.push_back(-999999);
  muonData->tr_r_me2_n.push_back(-999999);
  muonData->tr_phi_me2_n.push_back(-999999);

  muonData->tr_z_mb1.push_back(-999999);
  muonData->tr_phi_mb1.push_back(-999999);
  muonData->tr_r_me1_p.push_back(-999999);
  muonData->tr_phi_me1_p.push_back(-999999);
  muonData->tr_r_me1_n.push_back(-999999);
  muonData->tr_phi_me1_n.push_back(-999999);
}

void L1MuonRecoTreeProducer::endJob ( void  )

overrideprivatevirtual

Reimplemented from edm::EDAnalyzer.

Definition at line 1340 of file L1MuonRecoTreeProducer.cc.

References muon.

{ delete muon; }

void L1MuonRecoTreeProducer::endRun ( const edm::Run &  run, const edm::EventSetup &  eventSetup  )

overrideprivatevirtual

Reimplemented from edm::EDAnalyzer.

Definition at line 1335 of file L1MuonRecoTreeProducer.cc.

{}

FreeTrajectoryState L1MuonRecoTreeProducer::freeTrajStateMuon ( reco::TrackRef  track )

private

Definition at line 1273 of file L1MuonRecoTreeProducer.cc.

References theBField.

Referenced by cylExtrapTrkSam(), and surfExtrapTrkSam().

                                                                                {
  GlobalPoint innerPoint(track->innerPosition().x(), track->innerPosition().y(), track->innerPosition().z());
  GlobalVector innerVec(track->innerMomentum().x(), track->innerMomentum().y(), track->innerMomentum().z());

  FreeTrajectoryState recoStart(innerPoint, innerVec, track->charge(), &*theBField);

  return recoStart;
}

double L1MuonRecoTreeProducer::match_trigger	(	std::vector< int > &	trigIndices,
		const trigger::TriggerObjectCollection &	trigObjs,
		edm::Handle< trigger::TriggerEvent > &	triggerEvent,
		const reco::Muon &	mu
	)

Definition at line 226 of file L1MuonRecoTreeProducer.cc.

References HLT_FULL_cff::deltaR, hltConfig_, HLTConfigProvider::moduleLabels(), HLTConfigProvider::size(), and triggerProcessLabel_.

Referenced by analyze().

                                                                 {
  double matchDeltaR = 9999;

  for (size_t iTrigIndex = 0; iTrigIndex < trigIndices.size(); ++iTrigIndex) {
    int triggerIndex = trigIndices[iTrigIndex];
    const std::vector<std::string> moduleLabels(hltConfig_.moduleLabels(triggerIndex));
    // find index of the last module:
    const unsigned moduleIndex = hltConfig_.size(triggerIndex) - 2;
    // find index of HLT trigger name:
    const unsigned hltFilterIndex =
        triggerEvent->filterIndex(edm::InputTag(moduleLabels[moduleIndex], "", triggerProcessLabel_));

    if (hltFilterIndex < triggerEvent->sizeFilters()) {
      const trigger::Keys triggerKeys(triggerEvent->filterKeys(hltFilterIndex));
      const trigger::Vids triggerVids(triggerEvent->filterIds(hltFilterIndex));

      const unsigned nTriggers = triggerVids.size();
      for (size_t iTrig = 0; iTrig < nTriggers; ++iTrig) {
        // loop over all trigger objects:
        const trigger::TriggerObject trigObject = trigObjs[triggerKeys[iTrig]];

        double dRtmp = deltaR(mu, trigObject);

        if (dRtmp < matchDeltaR) {
          matchDeltaR = dRtmp;
        }

      }  // loop over different trigger objects
    }    // if trigger is in event (should apply hltFilter with used trigger...)
  }      // loop over muon candidates

  return matchDeltaR;
}

TrajectoryStateOnSurface L1MuonRecoTreeProducer::surfExtrapTrkSam	(	reco::TrackRef	track,
		double	z
	)

Definition at line 1259 of file L1MuonRecoTreeProducer.cc.

References Plane::build(), freeTrajStateMuon(), TrajectoryStateOnSurface::isValid(), propagatorAlong, propagatorOpposite, and makeMuonMisalignmentScenario::rot.

Referenced by analyze().

                                                                                              {
  Plane::PositionType pos(0, 0, z);
  Plane::RotationType rot;
  Plane::PlanePointer myPlane = Plane::build(pos, rot);

  FreeTrajectoryState recoStart = freeTrajStateMuon(track);
  TrajectoryStateOnSurface recoProp;
  recoProp = propagatorAlong->propagate(recoStart, *myPlane);
  if (!recoProp.isValid()) {
    recoProp = propagatorOpposite->propagate(recoStart, *myPlane);
  }
  return recoProp;
}

Member Data Documentation

edm::ESHandle<CSCGeometry> L1MuonRecoTreeProducer::cscGeom

private

Definition at line 154 of file L1MuonRecoTreeProducer.cc.

Referenced by analyze().

edm::ESHandle<DTGeometry> L1MuonRecoTreeProducer::dtGeom

private

Definition at line 158 of file L1MuonRecoTreeProducer.cc.

Referenced by analyze().

edm::Service<TFileService> L1MuonRecoTreeProducer::fs_

private

Definition at line 176 of file L1MuonRecoTreeProducer.cc.

Referenced by L1MuonRecoTreeProducer().

HLTConfigProvider L1MuonRecoTreeProducer::hltConfig_

private

Definition at line 149 of file L1MuonRecoTreeProducer.cc.

Referenced by beginRun(), and match_trigger().

std::vector<int> L1MuonRecoTreeProducer::isoTriggerIndices_

private

Definition at line 146 of file L1MuonRecoTreeProducer.cc.

Referenced by analyze(), and beginRun().

std::vector<std::string> L1MuonRecoTreeProducer::isoTriggerNames_

private

Definition at line 143 of file L1MuonRecoTreeProducer.cc.

Referenced by beginRun(), and L1MuonRecoTreeProducer().

unsigned L1MuonRecoTreeProducer::maxMuon_

private

Definition at line 133 of file L1MuonRecoTreeProducer.cc.

Referenced by analyze(), and L1MuonRecoTreeProducer().

unsigned L1MuonRecoTreeProducer::maxRpcHit_

private

Definition at line 134 of file L1MuonRecoTreeProducer.cc.

Referenced by analyze(), and L1MuonRecoTreeProducer().

L1Analysis::L1AnalysisRecoMuon* L1MuonRecoTreeProducer::muon

Definition at line 126 of file L1MuonRecoTreeProducer.cc.

Referenced by analyze(), endJob(), and L1MuonRecoTreeProducer().

L1Analysis::L1AnalysisRecoMuonDataFormat* L1MuonRecoTreeProducer::muonData

Definition at line 127 of file L1MuonRecoTreeProducer.cc.

Referenced by analyze(), empty_global(), empty_hlt(), empty_standalone(), empty_tracker(), and L1MuonRecoTreeProducer().

edm::InputTag L1MuonRecoTreeProducer::muonTag_

private

Definition at line 184 of file L1MuonRecoTreeProducer.cc.

Referenced by analyze(), and L1MuonRecoTreeProducer().

edm::ESHandle<Propagator> L1MuonRecoTreeProducer::propagatorAlong

private

Definition at line 170 of file L1MuonRecoTreeProducer.cc.

Referenced by analyze(), cylExtrapTrkSam(), and surfExtrapTrkSam().

edm::ESHandle<Propagator> L1MuonRecoTreeProducer::propagatorOpposite

private

Definition at line 171 of file L1MuonRecoTreeProducer.cc.

Referenced by analyze(), cylExtrapTrkSam(), and surfExtrapTrkSam().

edm::ESHandle<RPCGeometry> L1MuonRecoTreeProducer::rpcGeom

private

Definition at line 161 of file L1MuonRecoTreeProducer.cc.

Referenced by analyze().

L1Analysis::L1AnalysisRecoRpcHit* L1MuonRecoTreeProducer::rpcHit

Definition at line 129 of file L1MuonRecoTreeProducer.cc.

Referenced by analyze(), and L1MuonRecoTreeProducer().

L1Analysis::L1AnalysisRecoRpcHitDataFormat* L1MuonRecoTreeProducer::rpcHitData

Definition at line 130 of file L1MuonRecoTreeProducer.cc.

Referenced by analyze(), and L1MuonRecoTreeProducer().

edm::InputTag L1MuonRecoTreeProducer::rpcHitTag_

private

Definition at line 185 of file L1MuonRecoTreeProducer.cc.

Referenced by analyze(), and L1MuonRecoTreeProducer().

bool L1MuonRecoTreeProducer::runOnPostLS1_

private

Definition at line 181 of file L1MuonRecoTreeProducer.cc.

Referenced by analyze(), and L1MuonRecoTreeProducer().

edm::ESHandle<MagneticField> L1MuonRecoTreeProducer::theBField

private

Definition at line 164 of file L1MuonRecoTreeProducer.cc.

Referenced by analyze(), and freeTrajStateMuon().

edm::ESHandle<GlobalTrackingGeometry> L1MuonRecoTreeProducer::theTrackingGeometry

private

Definition at line 167 of file L1MuonRecoTreeProducer.cc.

Referenced by analyze().

TTree* L1MuonRecoTreeProducer::tree_

private

Definition at line 179 of file L1MuonRecoTreeProducer.cc.

Referenced by analyze(), and L1MuonRecoTreeProducer().

std::vector<int> L1MuonRecoTreeProducer::triggerIndices_

private

Definition at line 147 of file L1MuonRecoTreeProducer.cc.

Referenced by analyze(), and beginRun().

bool L1MuonRecoTreeProducer::triggerMatching_

private

Definition at line 140 of file L1MuonRecoTreeProducer.cc.

Referenced by analyze(), beginRun(), and L1MuonRecoTreeProducer().

double L1MuonRecoTreeProducer::triggerMaxDeltaR_

private

Definition at line 148 of file L1MuonRecoTreeProducer.cc.

Referenced by analyze(), and L1MuonRecoTreeProducer().

std::vector<std::string> L1MuonRecoTreeProducer::triggerNames_

private

Definition at line 144 of file L1MuonRecoTreeProducer.cc.

Referenced by beginRun(), and L1MuonRecoTreeProducer().

std::string L1MuonRecoTreeProducer::triggerProcessLabel_

private

Definition at line 142 of file L1MuonRecoTreeProducer.cc.

Referenced by analyze(), beginRun(), L1MuonRecoTreeProducer(), and match_trigger().

edm::InputTag L1MuonRecoTreeProducer::triggerSummaryLabel_

private

Definition at line 141 of file L1MuonRecoTreeProducer.cc.

Referenced by analyze(), and L1MuonRecoTreeProducer().