CMS 3D CMS Logo
List of all members | Public Member Functions | Public Attributes | Private Types | Private Member Functions | Private Attributes
L1MuonRecoTreeProducer Class Reference
Inheritance diagram for L1MuonRecoTreeProducer:
edm::EDAnalyzer edm::EDConsumerBase

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 &&)=default
 
 EDConsumerBase (EDConsumerBase const &)=delete
 
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 & operator= (EDConsumerBase &&)=default
 
EDConsumerBase const & operator= (EDConsumerBase const &)=delete
 
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
EDGetToken consumes (const TypeToGet &id, edm::InputTag const &tag)
 
template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType > consumes (edm::InputTag const &tag)
 
template<BranchType B = InEvent>
EDConsumerBaseAdaptor< B > consumes (edm::InputTag tag) noexcept
 
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<Transition Tr = Transition::Event>
constexpr auto esConsumes () noexcept
 
template<typename ESProduct , typename ESRecord , Transition Tr = Transition::Event>
auto esConsumes (ESInputTag const &tag)
 
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...
 
EDGetToken mayConsume (const TypeToGet &id, edm::InputTag const &tag)
 
template<BranchType B>
EDGetToken mayConsume (const TypeToGet &id, edm::InputTag const &tag)
 
template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType > mayConsume (edm::InputTag const &tag)
 

Detailed Description

Definition at line 102 of file L1MuonRecoTreeProducer.cc.

Member Enumeration Documentation

◆ anonymous enum

anonymous enum
private
Enumerator
GL_MUON 
SA_MUON 
TR_MUON 
TRSA_MUON 

Definition at line 151 of file L1MuonRecoTreeProducer.cc.

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

Constructor & Destructor Documentation

◆ L1MuonRecoTreeProducer()

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

Definition at line 188 of file L1MuonRecoTreeProducer.cc.

188  {
189  maxMuon_ = iConfig.getParameter<unsigned int>("maxMuon");
190  maxRpcHit_ = iConfig.getParameter<unsigned int>("maxMuon");
191 
192  muonTag_ = iConfig.getParameter<edm::InputTag>("muonTag");
193  rpcHitTag_ = iConfig.getParameter<edm::InputTag>("rpcHitTag");
194 
195  runOnPostLS1_ = iConfig.getParameter<bool>("runOnPostLS1");
196 
197  muon = new L1Analysis::L1AnalysisRecoMuon();
198  muonData = muon->getData();
199 
200  rpcHit = new L1Analysis::L1AnalysisRecoRpcHit();
201  rpcHitData = rpcHit->getData();
202 
203  // set up output
204  tree_ = fs_->make<TTree>("MuonRecoTree", "MuonRecoTree");
205  tree_->Branch("Muon", "L1Analysis::L1AnalysisRecoMuonDataFormat", &muonData, 32000, 3);
206  tree_->Branch("RpcHit", "L1Analysis::L1AnalysisRecoRpcHitDataFormat", &rpcHitData, 32000, 3);
207 
208  //---------------------------------------------------------------------------
209  // TRIGGER MATCHING
210  // member variables needed for matching, if using reco muons instead of PAT
211  //---------------------------------------------------------------------------
212  triggerMatching_ = iConfig.getUntrackedParameter<bool>("triggerMatching");
213  triggerSummaryLabel_ = iConfig.getParameter<edm::InputTag>("triggerSummaryLabel");
214  triggerProcessLabel_ = iConfig.getUntrackedParameter<std::string>("triggerProcessLabel");
215  triggerNames_ = iConfig.getParameter<std::vector<std::string> >("triggerNames");
216  isoTriggerNames_ = iConfig.getParameter<std::vector<std::string> >("isoTriggerNames");
217  triggerMaxDeltaR_ = iConfig.getParameter<double>("triggerMaxDeltaR");
218 }

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

◆ ~L1MuonRecoTreeProducer()

L1MuonRecoTreeProducer::~L1MuonRecoTreeProducer ( )
override

Definition at line 220 of file L1MuonRecoTreeProducer.cc.

220 {}

Member Function Documentation

◆ analyze()

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

Implements edm::EDAnalyzer.

Definition at line 398 of file L1MuonRecoTreeProducer.cc.

398  {
399  float pig = TMath::Pi();
400 
401  muon->Reset();
402  rpcHit->Reset();
403 
404  //GP start
405  // Get the CSC Geometry
406  iSetup.get<MuonGeometryRecord>().get(cscGeom);
407  //GP end
408 
409  // Get the DT Geometry from the setup
410  iSetup.get<MuonGeometryRecord>().get(dtGeom);
411 
412  // Get the RPC Geometry from the setup
413  iSetup.get<MuonGeometryRecord>().get(rpcGeom);
414 
415  //Get the Magnetic field from the setup
416  iSetup.get<IdealMagneticFieldRecord>().get(theBField);
417 
418  // Get the GlobalTrackingGeometry from the setup
419  iSetup.get<GlobalTrackingGeometryRecord>().get(theTrackingGeometry);
420 
421  edm::Handle<RPCRecHitCollection> rpcRecHits;
422 
423  iEvent.getByLabel(rpcHitTag_, rpcRecHits);
424  if (!rpcRecHits.isValid()) {
425  edm::LogInfo("L1Prompt") << "can't find RPCRecHitCollection with label " << rpcHitTag_.label();
426  } else {
427  RPCRecHitCollection::const_iterator recHitIt = rpcRecHits->begin();
428  RPCRecHitCollection::const_iterator recHitEnd = rpcRecHits->end();
429 
430  int iRpcRecHits = 0;
431 
432  for (; recHitIt != recHitEnd; ++recHitIt) {
433  if ((unsigned int)iRpcRecHits > maxRpcHit_ - 1)
434  continue;
435 
436  int cls = recHitIt->clusterSize();
437  int firststrip = recHitIt->firstClusterStrip();
438  int bx = recHitIt->BunchX();
439 
440  RPCDetId rpcId = recHitIt->rpcId();
441  int region = rpcId.region();
442  int stat = rpcId.station();
443  int sect = rpcId.sector();
444  int layer = rpcId.layer();
445  int subsector = rpcId.subsector();
446  int roll = rpcId.roll();
447  int ring = rpcId.ring();
448 
449  LocalPoint recHitPosLoc = recHitIt->localPosition();
450  const BoundPlane &RPCSurface = rpcGeom->roll(rpcId)->surface();
451  GlobalPoint recHitPosGlob = RPCSurface.toGlobal(recHitPosLoc);
452 
453  float xLoc = recHitPosLoc.x();
454  float phiGlob = recHitPosGlob.phi();
455 
456  rpcHitData->region.push_back(region);
457  rpcHitData->clusterSize.push_back(cls);
458  rpcHitData->strip.push_back(firststrip);
459  rpcHitData->bx.push_back(bx);
460  rpcHitData->xLoc.push_back(xLoc);
461  rpcHitData->phiGlob.push_back(phiGlob);
462  rpcHitData->station.push_back(stat);
463  rpcHitData->sector.push_back(sect);
464  rpcHitData->layer.push_back(layer);
465  rpcHitData->subsector.push_back(subsector);
466  rpcHitData->roll.push_back(roll);
467  rpcHitData->ring.push_back(ring);
468  rpcHitData->muonId.push_back(-999); // CB set to invalid now, updated when looking at muon info
469 
470  iRpcRecHits++;
471  }
472 
473  rpcHitData->nRpcHits = iRpcRecHits;
474  }
475 
476  // Get the muon candidates
477  edm::Handle<reco::MuonCollection> mucand;
478  iEvent.getByLabel(muonTag_, mucand);
479  if (!mucand.isValid()) {
480  edm::LogInfo("L1Prompt") << "can't find Muon Collection with label " << muonTag_.label();
481  return;
482  }
483 
484  // Get the beamspot
485  edm::Handle<reco::BeamSpot> beamSpot;
486  iEvent.getByLabel("offlineBeamSpot", beamSpot);
487 
488  // Get the primary vertices
489  edm::Handle<std::vector<reco::Vertex> > vertex;
490  iEvent.getByLabel(edm::InputTag("offlinePrimaryVertices"), vertex);
491 
492  // Get the propagators
493  iSetup.get<TrackingComponentsRecord>().get("SmartPropagatorAny", propagatorAlong);
494  iSetup.get<TrackingComponentsRecord>().get("SmartPropagatorAnyOpposite", propagatorOpposite);
495 
496  for (reco::MuonCollection::const_iterator imu = mucand->begin();
497  // for(pat::MuonCollection::const_iterator imu = mucand->begin();
498  imu != mucand->end() && (unsigned)muonData->nMuons < maxMuon_;
499  imu++) {
500  int type = 0;
501  if (imu->isGlobalMuon())
502  type = type + 1;
503  if (imu->isStandAloneMuon())
504  type = type + 2;
505  if (imu->isTrackerMuon())
506  type = type + 4;
507  if (imu->isCaloMuon())
508  type = type + 8;
509 
510  bool isTIGHT =
511  (!vertex->empty() && imu->isGlobalMuon() && imu->globalTrack()->normalizedChi2() < 10. &&
512  imu->globalTrack()->hitPattern().numberOfValidMuonHits() > 0 && imu->numberOfMatchedStations() > 1 &&
513  fabs(imu->innerTrack()->dxy(vertex->at(0).position())) < 0.2 &&
514  fabs(imu->innerTrack()->dz(vertex->at(0).position())) < 0.5 &&
515  imu->innerTrack()->hitPattern().numberOfValidPixelHits() > 0 &&
516  imu->innerTrack()->hitPattern().trackerLayersWithMeasurement() > 5);
517 
518  if (isTIGHT)
519  type = type + 16;
520  if (imu->isPFMuon())
521  type = type + 32;
522 
523  muonData->howmanytypes.push_back(type); // muon type is counting calo muons and multiple assignments
524 
525  // bool type identifiers are exclusive. is this correct? CB to check
526  bool isSA = (!imu->isGlobalMuon() && imu->isStandAloneMuon() && !imu->isTrackerMuon());
527  bool isTR = (!imu->isGlobalMuon() && imu->isTrackerMuon() && !imu->isStandAloneMuon());
528  bool isGL = (imu->isGlobalMuon()); //&&!(imu->isStandAloneMuon())&&!(imu->isTrackerMuon()));
529  bool isTRSA = (!imu->isGlobalMuon() && imu->isStandAloneMuon() && imu->isTrackerMuon());
530 
531  // How we fill this. We have 3 blocks of variables: muons_; muons_sa_; muons_tr_
532  // GL SA TR Description muon_type Bool id Filled Vars Variables to -99999
533  // 0 0 0 Muon does not exist / / / /
534  // 0 0 1 It is a Tracker muon TR_MUON isTR muons_tr_ muons_sa_,muons_
535  // 0 1 0 It is a SA muon SA_MUON isSA muons_sa muons_tr,muons_
536  // 0 1 1 It is a SA+Tracker muon TRSA_MUON isTRSA muons_sa,muons_tr muons_
537  // 1 0 0 It is a Global only muon GL_MUON isGL muons_,muons_sa,muons_tr /
538  // 1 0 1 Gl w/out SA cannot exist / / / /
539  // 1 1 0 Gl+SA (no trk-mu match) GL_MUON isGL muons_,muons_sa,muons_tr none
540  // 1 1 1 GL+SA+Tr (all matched) GL_MUON isGL muons_,muons_sa,muons_tr none
541 
542  //---------------------------------------------------------------------
543  // TRIGGER MATCHING:
544  // if specified the reconstructed muons are matched to a trigger
545  //---------------------------------------------------------------------
546  if (triggerMatching_) {
547  double isoMatchDeltaR = 9999.;
548  double matchDeltaR = 9999.;
549  int hasIsoTriggered = 0;
550  int hasTriggered = 0;
551 
552  // first check if the trigger results are valid:
553  edm::Handle<edm::TriggerResults> triggerResults;
554  iEvent.getByLabel(edm::InputTag("TriggerResults", "", triggerProcessLabel_), triggerResults);
555 
556  if (triggerResults.isValid()) {
557  edm::Handle<trigger::TriggerEvent> triggerEvent;
558  iEvent.getByLabel(triggerSummaryLabel_, triggerEvent);
559  if (triggerEvent.isValid()) {
560  // get trigger objects:
561  const trigger::TriggerObjectCollection triggerObjects = triggerEvent->getObjects();
562 
563  matchDeltaR = match_trigger(triggerIndices_, triggerObjects, triggerEvent, (*imu));
564  if (matchDeltaR < triggerMaxDeltaR_)
565  hasTriggered = 1;
566 
567  isoMatchDeltaR = match_trigger(isoTriggerIndices_, triggerObjects, triggerEvent, (*imu));
568  if (isoMatchDeltaR < triggerMaxDeltaR_)
569  hasIsoTriggered = 1;
570  } // end if (triggerEvent.isValid())
571  } // end if (triggerResults.isValid())
572 
573  // fill trigger matching variables:
574  muonData->hlt_isomu.push_back(hasIsoTriggered);
575  muonData->hlt_mu.push_back(hasTriggered);
576  muonData->hlt_isoDeltaR.push_back(isoMatchDeltaR);
577  muonData->hlt_deltaR.push_back(matchDeltaR);
578  } else {
579  empty_hlt();
580  } // end if (triggerMatching_)
581 
582  if (isGL || isTR || isSA || isTRSA) {
583  muonData->nMuons = muonData->nMuons + 1;
584  if (isTR)
585  muonData->type.push_back(TR_MUON);
586  if (isGL)
587  muonData->type.push_back(GL_MUON);
588  if (isSA)
589  muonData->type.push_back(SA_MUON);
590  if (isTRSA)
591  muonData->type.push_back(TRSA_MUON);
592  if (!isGL)
593  empty_global();
594  if (!isTR && !isGL && !isTRSA)
595  empty_tracker();
596  if (!isSA && !isGL && !isTRSA)
597  empty_standalone();
598 
599  // begin GP
600  //---------------------------------------------------------------------
601  // RECHIT information in CSC: only for standalone/global muons!
602  //---------------------------------------------------------------------
603  // An artificial rank to sort RecHits
604  // the closer RecHit to key station/layer -> the smaller the rank
605  // KK's original idea
606  const int lutCSC[4][6] = {
607  {26, 24, 22, 21, 23, 25}, {6, 4, 2, 1, 3, 5}, {16, 14, 12, 11, 13, 15}, {36, 34, 32, 31, 33, 35}};
608 
609  int globalTypeRCH = -999999;
610  // float localEtaRCH = -999999;
611  // float localPhiRCH = -999999;
612  float globalEtaRCH = -999999;
613  float globalPhiRCH = -999999;
614 
615  if (isSA || isGL) {
616  trackingRecHit_iterator hit = imu->outerTrack()->recHitsBegin();
617  trackingRecHit_iterator hitEnd = imu->outerTrack()->recHitsEnd();
618 
619  for (; hit != hitEnd; ++hit) {
620  if (!((*hit)->isValid()))
621  continue;
622 
623  // Hardware ID of the RecHit (in terms of wire/strips/chambers)
624  DetId detid = (*hit)->geographicalId();
625 
626  // Interested in muon systems only
627  // Look only at CSC Hits (CSC id is 2)
628  if (detid.det() != DetId::Muon)
629  continue;
630  if (detid.subdetId() != MuonSubdetId::CSC)
631  continue;
632 
633  CSCDetId id(detid.rawId());
634  //std::cout << "before Lut id.station = " <<id.station()
635  // << " id.layer() = " << id.layer() << " globalTypeRCH = "
636  // << globalTypeRCH << std::endl;
637  // another sanity check
638  if (id.station() < 1)
639  continue;
640 
641  // Look up some stuff specific to CSCRecHit2D
642  // std::cout << " typeid().name() " << typeid(**hit).name() << std::endl;
643  // const CSCRecHit2D* CSChit =dynamic_cast<const CSCRecHit2D*>(&**hit);
644 
645  const CSCSegment *cscSegment = dynamic_cast<const CSCSegment *>(&**hit);
646  //std::cout << "cscSegment = " << cscSegment << std::endl;
647  if (cscSegment == nullptr)
648  continue;
649  // const CSCRecHit2D* CSChit =(CSCRecHit2D*)(&**hit);
650 
651  // std::cout << " after CSCRecHit2D, CSChit = " << CSChit << std::endl;
652  // LocalPoint rhitlocal = CSChit->localPosition();
653  LocalPoint rhitlocal = cscSegment->localPosition();
654 
655  // for debugging purpouses
656  //if (printLevel > 0) {
657  //std::cout << "!!!!rhitlocal.phi = "<<rhitlocal.phi() << std::endl;
658  //std::cout << "rhitlocal.y="<<rhitlocal.y();
659  //std::cout << "rhitlocal.z="<<rhitlocal.z();
660  //}
661 
662  GlobalPoint gp = GlobalPoint(0.0, 0.0, 0.0);
663 
664  const CSCChamber *cscchamber = cscGeom->chamber(id);
665 
666  if (!cscchamber)
667  continue;
668 
669  gp = cscchamber->toGlobal(rhitlocal);
670 
671  // identify the rechit position
672  //int pos = ( ((id.station()-1)*6) + (id.layer()-1) ) + (MAX_CSC_RECHIT*whichMuon);
673 
674  // --------------------------------------------------
675  // this part has to be deprecated once we are sure of
676  // the TMatrixF usage ;)
677  // fill the rechits array
678  //rchEtaList[pos] = gp.eta();
679  //rchPhiList[pos] = gp.phi();
680  //float phi02PI = gp.phi();
681  //if (gp.phi() < 0) phi02PI += (2*PI);
682 
683  //rchPhiList_02PI[pos] = phi02PI;
684  // --------------------------------------------------
685 
686  // --------------------------------------------------
687  // See if this hit is closer to the "key" position
688  //if( lutCSC[id.station()-1][id.layer()-1]<globalTypeRCH || globalTypeRCH<0 ){
689  if (lutCSC[id.station() - 1][3 - 1] < globalTypeRCH || globalTypeRCH < 0) {
690  //globalTypeRCH = lutCSC[id.station()-1][id.layer()-1];
691  globalTypeRCH = lutCSC[id.station() - 1][3 - 1];
692  // localEtaRCH = rhitlocal.eta();
693  // localPhiRCH = rhitlocal.phi();
694  globalEtaRCH = gp.eta();
695  globalPhiRCH = gp.phi(); // phi from -pi to pi
696  //std::cout << "globalEtaRCH = " <<globalEtaRCH
697  // << " globalPhiRCH = " << globalPhiRCH << std::endl;
698  if (globalPhiRCH < 0)
699  globalPhiRCH = globalPhiRCH + 2 * pig; // convert to [0; 2pi]
700  }
701  // --------------------------------------------------
702  }
703 
704  hit = imu->outerTrack()->recHitsBegin();
705 
706  for (; hit != hitEnd; hit++) {
707  if (!((*hit)->isValid()))
708  continue;
709 
710  DetId detId = (*hit)->geographicalId();
711 
712  if (detId.det() != DetId::Muon)
713  continue;
714  if (detId.subdetId() != MuonSubdetId::RPC)
715  continue;
716 
717  RPCDetId rpcId = (RPCDetId)(*hit)->geographicalId();
718 
719  int region = rpcId.region();
720  int stat = rpcId.station();
721  int sect = rpcId.sector();
722  int layer = rpcId.layer();
723  int subsector = rpcId.subsector();
724  int roll = rpcId.roll();
725  int ring = rpcId.ring();
726 
727  float xLoc = (*hit)->localPosition().x();
728 
729  for (int iRpcHit = 0; iRpcHit < rpcHitData->nRpcHits; ++iRpcHit) {
730  if (region == rpcHitData->region.at(iRpcHit) && stat == rpcHitData->station.at(iRpcHit) &&
731  sect == rpcHitData->sector.at(iRpcHit) && layer == rpcHitData->layer.at(iRpcHit) &&
732  subsector == rpcHitData->subsector.at(iRpcHit) && roll == rpcHitData->roll.at(iRpcHit) &&
733  ring == rpcHitData->ring.at(iRpcHit) && fabs(xLoc - rpcHitData->xLoc.at(iRpcHit)) < 0.01)
734 
735  rpcHitData->muonId.at(iRpcHit) = muonData->nMuons; // CB rpc hit belongs to mu imu
736  // due to cleaning as in 2012 1 rpc hit belogns to 1 mu
737  }
738  }
739  }
740  // at the end of the loop, write only the best rechit
741  // if(globalPhiRCH > -100.) std::cout << "globalPhiRCH = " << globalPhiRCH
742  // << "globalEtaRCH = " << globalEtaRCH << std::endl;
743  muonData->rchCSCtype.push_back(globalTypeRCH);
744  muonData->rchPhi.push_back(globalPhiRCH);
745  muonData->rchEta.push_back(globalEtaRCH);
746  // end GP
747 
748  if (isGL) {
749  // Filling calo energies and calo times
750  if (imu->isEnergyValid()) {
751  reco::MuonEnergy muon_energy;
752  muon_energy = imu->calEnergy();
753  muonData->calo_energy.push_back(muon_energy.tower);
754  muonData->calo_energy3x3.push_back(muon_energy.towerS9);
755  muonData->ecal_time.push_back(muon_energy.ecal_time);
756  muonData->ecal_terr.push_back(muon_energy.ecal_timeError);
757  muonData->hcal_time.push_back(muon_energy.hcal_time);
758  muonData->hcal_terr.push_back(muon_energy.hcal_timeError);
759  } else {
760  muonData->calo_energy.push_back(-999999);
761  muonData->calo_energy3x3.push_back(-999999);
762  muonData->ecal_time.push_back(-999999);
763  muonData->ecal_terr.push_back(-999999);
764  muonData->hcal_time.push_back(-999999);
765  muonData->hcal_terr.push_back(-999999);
766  }
767 
768  // Filling muon time data
769  if (imu->isTimeValid()) {
770  reco::MuonTime muon_time;
771  muon_time = imu->time();
772  muonData->time_dir.push_back(muon_time.direction());
773  muonData->time_inout.push_back(muon_time.timeAtIpInOut);
774  muonData->time_inout_err.push_back(muon_time.timeAtIpInOutErr);
775  muonData->time_outin.push_back(muon_time.timeAtIpOutIn);
776  muonData->time_outin_err.push_back(muon_time.timeAtIpOutInErr);
777  } else {
778  muonData->time_dir.push_back(-999999);
779  muonData->time_inout.push_back(-999999);
780  muonData->time_inout_err.push_back(-999999);
781  muonData->time_outin.push_back(-999999);
782  muonData->time_outin_err.push_back(-999999);
783  }
784 
785  // Use the track now, and make a transient track out of it (for extrapolations)
786  reco::TrackRef glb_mu = imu->globalTrack();
787  reco::TransientTrack ttrack(*glb_mu, &*theBField, theTrackingGeometry);
788 
789  // Track quantities
790  muonData->ch.push_back(glb_mu->charge());
791  muonData->pt.push_back(glb_mu->pt());
792  muonData->p.push_back(glb_mu->p());
793  muonData->eta.push_back(glb_mu->eta());
794  muonData->phi.push_back(glb_mu->phi());
795  muonData->normchi2.push_back(glb_mu->normalizedChi2());
796  muonData->validhits.push_back(glb_mu->numberOfValidHits());
797  muonData->numberOfMatchedStations.push_back(imu->numberOfMatchedStations());
798  muonData->numberOfValidMuonHits.push_back(glb_mu->hitPattern().numberOfValidMuonHits());
799 
800  // Extrapolation to IP
801  if (ttrack.impactPointTSCP().isValid()) {
802  muonData->imp_point_x.push_back(ttrack.impactPointTSCP().position().x());
803  muonData->imp_point_y.push_back(ttrack.impactPointTSCP().position().y());
804  muonData->imp_point_z.push_back(ttrack.impactPointTSCP().position().z());
805  muonData->imp_point_p.push_back(sqrt(ttrack.impactPointTSCP().position().mag()));
806  muonData->imp_point_pt.push_back(sqrt(ttrack.impactPointTSCP().position().perp2()));
807  } else {
808  muonData->imp_point_x.push_back(-999999);
809  muonData->imp_point_y.push_back(-999999);
810  muonData->imp_point_z.push_back(-999999);
811  muonData->imp_point_p.push_back(-999999);
812  muonData->imp_point_pt.push_back(-999999);
813  }
814 
815  if (!runOnPostLS1_) {
816  //Extrapolation to HB
817  TrajectoryStateOnSurface tsos;
818  tsos = cylExtrapTrkSam(glb_mu, 235);
819  if (tsos.isValid()) {
820  double xx = tsos.globalPosition().x();
821  double yy = tsos.globalPosition().y();
822  double zz = tsos.globalPosition().z();
823  muonData->z_hb.push_back(zz);
824  double rr = sqrt(xx * xx + yy * yy);
825  double cosphi = xx / rr;
826  if (yy >= 0)
827  muonData->phi_hb.push_back(acos(cosphi));
828  else
829  muonData->phi_hb.push_back(2 * pig - acos(cosphi));
830  } else {
831  muonData->phi_hb.push_back(-999999);
832  muonData->z_hb.push_back(-999999);
833  }
834 
835  //Extrapolation to HE+
836  tsos = surfExtrapTrkSam(glb_mu, 479);
837  if (tsos.isValid()) {
838  double xx = tsos.globalPosition().x();
839  double yy = tsos.globalPosition().y();
840  double rr = sqrt(xx * xx + yy * yy);
841  muonData->r_he_p.push_back(rr);
842  double cosphi = xx / rr;
843  if (yy >= 0)
844  muonData->phi_he_p.push_back(acos(cosphi));
845  else
846  muonData->phi_he_p.push_back(2 * pig - acos(cosphi));
847  } else {
848  muonData->r_he_p.push_back(-999999);
849  muonData->phi_he_p.push_back(-999999);
850  }
851 
852  //Extrapolation to HE-
853  tsos = surfExtrapTrkSam(glb_mu, -479);
854  if (tsos.isValid()) {
855  double xx = tsos.globalPosition().x();
856  double yy = tsos.globalPosition().y();
857  double rr = sqrt(xx * xx + yy * yy);
858  muonData->r_he_n.push_back(rr);
859  double cosphi = xx / rr;
860  if (yy >= 0)
861  muonData->phi_he_n.push_back(acos(cosphi));
862  else
863  muonData->phi_he_n.push_back(2 * pig - acos(cosphi));
864  } else {
865  muonData->r_he_n.push_back(-999999);
866  muonData->phi_he_n.push_back(-999999);
867  }
868  } else {
869  muonData->phi_hb.push_back(-999999);
870  muonData->z_hb.push_back(-999999);
871  muonData->r_he_p.push_back(-999999);
872  muonData->phi_he_p.push_back(-999999);
873  muonData->r_he_n.push_back(-999999);
874  muonData->phi_he_n.push_back(-999999);
875  }
876  } // end of IF IS GLOBAL
877 
878  // If global muon or tracker muon, fill the tracker track quantities
879  if (isTR || isGL || isTRSA) {
880  // Take the tracker track and build a transient track out of it
881  reco::TrackRef tr_mu = imu->innerTrack();
882  reco::TransientTrack ttrack(*tr_mu, &*theBField, theTrackingGeometry);
883  // Fill track quantities
884  muonData->tr_ch.push_back(tr_mu->charge());
885  muonData->tr_pt.push_back(tr_mu->pt());
886  muonData->tr_p.push_back(tr_mu->p());
887  muonData->tr_eta.push_back(tr_mu->eta());
888  muonData->tr_phi.push_back(tr_mu->phi());
889  muonData->tr_normchi2.push_back(tr_mu->normalizedChi2());
890  muonData->tr_validhits.push_back(tr_mu->numberOfValidHits());
891  muonData->tr_validpixhits.push_back(tr_mu->hitPattern().numberOfValidPixelHits());
892  // find d0 from vertex position
894  edm::Handle<reco::BeamSpot> beamSpot;
895  iEvent.getByLabel("offlineBeamSpot", beamSpot);
896  muonData->tr_d0.push_back(tr_mu->dxy(beamSpot->position()));
897 
898  // Extrapolation to the IP
899  if (ttrack.impactPointTSCP().isValid()) {
900  muonData->tr_imp_point_x.push_back(ttrack.impactPointTSCP().position().x());
901  muonData->tr_imp_point_y.push_back(ttrack.impactPointTSCP().position().y());
902  muonData->tr_imp_point_z.push_back(ttrack.impactPointTSCP().position().z());
903  muonData->tr_imp_point_p.push_back(sqrt(ttrack.impactPointTSCP().position().mag()));
904  muonData->tr_imp_point_pt.push_back(sqrt(ttrack.impactPointTSCP().position().perp2()));
905  } else {
906  muonData->tr_imp_point_x.push_back(-999999);
907  muonData->tr_imp_point_y.push_back(-999999);
908  muonData->tr_imp_point_z.push_back(-999999);
909  muonData->tr_imp_point_p.push_back(-999999);
910  muonData->tr_imp_point_pt.push_back(-999999);
911  }
912 
913  if (!runOnPostLS1_) {
914  TrajectoryStateOnSurface tsos;
915  tsos = cylExtrapTrkSam(tr_mu, 410); // track at MB1 radius - extrapolation
916  if (tsos.isValid()) {
917  double xx = tsos.globalPosition().x();
918  double yy = tsos.globalPosition().y();
919  double zz = tsos.globalPosition().z();
920  muonData->tr_z_mb1.push_back(zz);
921  double rr = sqrt(xx * xx + yy * yy);
922  double cosphi = xx / rr;
923  if (yy >= 0)
924  muonData->tr_phi_mb1.push_back(acos(cosphi));
925  else
926  muonData->tr_phi_mb1.push_back(2 * pig - acos(cosphi));
927  } else {
928  muonData->tr_z_mb1.push_back(-999999);
929  muonData->tr_phi_mb1.push_back(-999999);
930  }
931 
932  tsos = cylExtrapTrkSam(tr_mu, 500); // track at MB2 radius - extrapolation
933  if (tsos.isValid()) {
934  double xx = tsos.globalPosition().x();
935  double yy = tsos.globalPosition().y();
936  double zz = tsos.globalPosition().z();
937  muonData->tr_z_mb2.push_back(zz);
938  double rr = sqrt(xx * xx + yy * yy);
939  double cosphi = xx / rr;
940  if (yy >= 0)
941  muonData->tr_phi_mb2.push_back(acos(cosphi));
942  else
943  muonData->tr_phi_mb2.push_back(2 * pig - acos(cosphi));
944  } else {
945  muonData->tr_z_mb2.push_back(-999999);
946  muonData->tr_phi_mb2.push_back(-999999);
947  }
948 
949  tsos = surfExtrapTrkSam(tr_mu, 630); // track at ME1+ plane - extrapolation
950  if (tsos.isValid()) {
951  double xx = tsos.globalPosition().x();
952  double yy = tsos.globalPosition().y();
953  double rr = sqrt(xx * xx + yy * yy);
954  muonData->tr_r_me1_p.push_back(rr);
955  double cosphi = xx / rr;
956  if (yy >= 0)
957  muonData->tr_phi_me1_p.push_back(acos(cosphi));
958  else
959  muonData->tr_phi_me1_p.push_back(2 * pig - acos(cosphi));
960  } else {
961  muonData->tr_r_me1_p.push_back(-999999);
962  muonData->tr_phi_me1_p.push_back(-999999);
963  }
964 
965  tsos = surfExtrapTrkSam(tr_mu, 790); // track at ME2+ plane - extrapolation
966  if (tsos.isValid()) {
967  double xx = tsos.globalPosition().x();
968  double yy = tsos.globalPosition().y();
969  double rr = sqrt(xx * xx + yy * yy);
970  muonData->tr_r_me2_p.push_back(rr);
971  double cosphi = xx / rr;
972  if (yy >= 0)
973  muonData->tr_phi_me2_p.push_back(acos(cosphi));
974  else
975  muonData->tr_phi_me2_p.push_back(2 * pig - acos(cosphi));
976  } else {
977  muonData->tr_r_me2_p.push_back(-999999);
978  muonData->tr_phi_me2_p.push_back(-999999);
979  }
980 
981  tsos = surfExtrapTrkSam(tr_mu, -630); // track at ME1- plane - extrapolation
982  if (tsos.isValid()) {
983  double xx = tsos.globalPosition().x();
984  double yy = tsos.globalPosition().y();
985  double rr = sqrt(xx * xx + yy * yy);
986  muonData->tr_r_me1_n.push_back(rr);
987  double cosphi = xx / rr;
988  if (yy >= 0)
989  muonData->tr_phi_me1_n.push_back(acos(cosphi));
990  else
991  muonData->tr_phi_me1_n.push_back(2 * pig - acos(cosphi));
992  } else {
993  muonData->tr_r_me1_n.push_back(-999999);
994  muonData->tr_phi_me1_n.push_back(-999999);
995  }
996 
997  tsos = surfExtrapTrkSam(tr_mu, -790); // track at ME2- plane - extrapolation
998  if (tsos.isValid()) {
999  double xx = tsos.globalPosition().x();
1000  double yy = tsos.globalPosition().y();
1001  double rr = sqrt(xx * xx + yy * yy);
1002  muonData->tr_r_me2_n.push_back(rr);
1003  double cosphi = xx / rr;
1004  if (yy >= 0)
1005  muonData->tr_phi_me2_n.push_back(acos(cosphi));
1006  else
1007  muonData->tr_phi_me2_n.push_back(2 * pig - acos(cosphi));
1008  } else {
1009  muonData->tr_r_me2_n.push_back(-999999);
1010  muonData->tr_phi_me2_n.push_back(-999999);
1011  }
1012  } else {
1013  muonData->tr_z_mb1.push_back(-999999);
1014  muonData->tr_phi_mb1.push_back(-999999);
1015  muonData->tr_z_mb2.push_back(-999999);
1016  muonData->tr_phi_mb2.push_back(-999999);
1017  muonData->tr_r_me1_p.push_back(-999999);
1018  muonData->tr_phi_me1_p.push_back(-999999);
1019  muonData->tr_r_me2_p.push_back(-999999);
1020  muonData->tr_phi_me2_p.push_back(-999999);
1021  muonData->tr_r_me1_n.push_back(-999999);
1022  muonData->tr_phi_me1_n.push_back(-999999);
1023  muonData->tr_r_me2_n.push_back(-999999);
1024  muonData->tr_phi_me2_n.push_back(-999999);
1025  }
1026 
1027  } // end of IF IS TRACKER
1028 
1029  // If global muon or sa muon, fill the sa track quantities
1030  if (isGL || isSA || isTRSA) {
1031  muonData->sa_nChambers.push_back(imu->numberOfChambers());
1032  muonData->sa_nMatches.push_back(imu->numberOfMatches());
1033 
1034  // Take the SA track and build a transient track out of it
1035  reco::TrackRef sa_mu = imu->outerTrack();
1036  reco::TransientTrack ttrack(*sa_mu, &*theBField, theTrackingGeometry);
1037 
1038  // Extrapolation to IP
1039  if (ttrack.impactPointTSCP().isValid()) {
1040  muonData->sa_imp_point_x.push_back(ttrack.impactPointTSCP().position().x());
1041  muonData->sa_imp_point_y.push_back(ttrack.impactPointTSCP().position().y());
1042  muonData->sa_imp_point_z.push_back(ttrack.impactPointTSCP().position().z());
1043  muonData->sa_imp_point_p.push_back(sqrt(ttrack.impactPointTSCP().position().mag()));
1044  muonData->sa_imp_point_pt.push_back(sqrt(ttrack.impactPointTSCP().position().perp2()));
1045  } else {
1046  muonData->sa_imp_point_x.push_back(-999999);
1047  muonData->sa_imp_point_y.push_back(-999999);
1048  muonData->sa_imp_point_z.push_back(-999999);
1049  muonData->sa_imp_point_p.push_back(-999999);
1050  muonData->sa_imp_point_pt.push_back(-999999);
1051  }
1052 
1053  // Extrapolation to MB2
1054 
1055  TrajectoryStateOnSurface tsos;
1056  tsos = cylExtrapTrkSam(sa_mu, 410); // track at MB1 radius - extrapolation
1057  if (tsos.isValid()) {
1058  double xx = tsos.globalPosition().x();
1059  double yy = tsos.globalPosition().y();
1060  double zz = tsos.globalPosition().z();
1061  muonData->sa_z_mb1.push_back(zz);
1062  double rr = sqrt(xx * xx + yy * yy);
1063  double cosphi = xx / rr;
1064  if (yy >= 0)
1065  muonData->sa_phi_mb1.push_back(acos(cosphi));
1066  else
1067  muonData->sa_phi_mb1.push_back(2 * pig - acos(cosphi));
1068  } else {
1069  muonData->sa_z_mb1.push_back(-999999);
1070  muonData->sa_phi_mb1.push_back(-999999);
1071  }
1072 
1073  tsos = cylExtrapTrkSam(sa_mu, 500); // track at MB2 radius - extrapolation
1074  if (tsos.isValid()) {
1075  double xx = tsos.globalPosition().x();
1076  double yy = tsos.globalPosition().y();
1077  double zz = tsos.globalPosition().z();
1078  muonData->sa_z_mb2.push_back(zz);
1079  double rr = sqrt(xx * xx + yy * yy);
1080  double cosphi = xx / rr;
1081  if (yy >= 0)
1082  muonData->sa_phi_mb2.push_back(acos(cosphi));
1083  else
1084  muonData->sa_phi_mb2.push_back(2 * pig - acos(cosphi));
1085  double abspseta = -log(tan(atan(fabs(rr / zz)) / 2.0));
1086  if (zz >= 0)
1087  muonData->sa_pseta.push_back(abspseta);
1088  else
1089  muonData->sa_pseta.push_back(-abspseta);
1090  } else {
1091  muonData->sa_z_mb2.push_back(-999999);
1092  muonData->sa_phi_mb2.push_back(-999999);
1093  muonData->sa_pseta.push_back(-999999);
1094  }
1095 
1096  tsos = surfExtrapTrkSam(sa_mu, 630); // track at ME1+ plane - extrapolation
1097  if (tsos.isValid()) {
1098  double xx = tsos.globalPosition().x();
1099  double yy = tsos.globalPosition().y();
1100  double rr = sqrt(xx * xx + yy * yy);
1101  muonData->sa_r_me1_p.push_back(rr);
1102  double cosphi = xx / rr;
1103  if (yy >= 0)
1104  muonData->sa_phi_me1_p.push_back(acos(cosphi));
1105  else
1106  muonData->sa_phi_me1_p.push_back(2 * pig - acos(cosphi));
1107  } else {
1108  muonData->sa_r_me1_p.push_back(-999999);
1109  muonData->sa_phi_me1_p.push_back(-999999);
1110  }
1111 
1112  // Extrapolation to ME2+
1113  tsos = surfExtrapTrkSam(sa_mu, 790);
1114  if (tsos.isValid()) {
1115  double xx = tsos.globalPosition().x();
1116  double yy = tsos.globalPosition().y();
1117  double rr = sqrt(xx * xx + yy * yy);
1118  muonData->sa_r_me2_p.push_back(rr);
1119  double cosphi = xx / rr;
1120  if (yy >= 0)
1121  muonData->sa_phi_me2_p.push_back(acos(cosphi));
1122  else
1123  muonData->sa_phi_me2_p.push_back(2 * pig - acos(cosphi));
1124  } else {
1125  muonData->sa_r_me2_p.push_back(-999999);
1126  muonData->sa_phi_me2_p.push_back(-999999);
1127  }
1128 
1129  tsos = surfExtrapTrkSam(sa_mu, -630); // track at ME1- plane - extrapolation
1130  if (tsos.isValid()) {
1131  double xx = tsos.globalPosition().x();
1132  double yy = tsos.globalPosition().y();
1133  double rr = sqrt(xx * xx + yy * yy);
1134  muonData->sa_r_me1_n.push_back(rr);
1135  double cosphi = xx / rr;
1136  if (yy >= 0)
1137  muonData->sa_phi_me1_n.push_back(acos(cosphi));
1138  else
1139  muonData->sa_phi_me1_n.push_back(2 * pig - acos(cosphi));
1140  } else {
1141  muonData->sa_r_me1_n.push_back(-999999);
1142  muonData->sa_phi_me1_n.push_back(-999999);
1143  }
1144 
1145  // Extrapolation to ME2-
1146  tsos = surfExtrapTrkSam(sa_mu, -790); // track at ME2- disk - extrapolation
1147  if (tsos.isValid()) {
1148  double xx = tsos.globalPosition().x();
1149  double yy = tsos.globalPosition().y();
1150  double rr = sqrt(xx * xx + yy * yy);
1151  muonData->sa_r_me2_n.push_back(rr);
1152  double cosphi = xx / rr;
1153  if (yy >= 0)
1154  muonData->sa_phi_me2_n.push_back(acos(cosphi));
1155  else
1156  muonData->sa_phi_me2_n.push_back(2 * pig - acos(cosphi));
1157  } else {
1158  muonData->sa_r_me2_n.push_back(-999999);
1159  muonData->sa_phi_me2_n.push_back(-999999);
1160  }
1161 
1162  // Extrapolation to HB
1163  tsos = cylExtrapTrkSam(sa_mu, 235); // track at HB radius - extrapolation
1164  if (tsos.isValid()) {
1165  double xx = tsos.globalPosition().x();
1166  double yy = tsos.globalPosition().y();
1167  double zz = tsos.globalPosition().z();
1168  muonData->sa_z_hb.push_back(zz);
1169  double rr = sqrt(xx * xx + yy * yy);
1170  double cosphi = xx / rr;
1171  if (yy >= 0)
1172  muonData->sa_phi_hb.push_back(acos(cosphi));
1173  else
1174  muonData->sa_phi_hb.push_back(2 * pig - acos(cosphi));
1175  } else {
1176  muonData->sa_z_hb.push_back(-999999);
1177  muonData->sa_phi_hb.push_back(-999999);
1178  }
1179 
1180  // Extrapolation to HE+
1181  tsos = surfExtrapTrkSam(sa_mu, 479);
1182  if (tsos.isValid()) {
1183  double xx = tsos.globalPosition().x();
1184  double yy = tsos.globalPosition().y();
1185  double rr = sqrt(xx * xx + yy * yy);
1186  muonData->sa_r_he_p.push_back(rr);
1187  double cosphi = xx / rr;
1188  if (yy >= 0)
1189  muonData->sa_phi_he_p.push_back(acos(cosphi));
1190  else
1191  muonData->sa_phi_he_p.push_back(2 * pig - acos(cosphi));
1192  } else {
1193  muonData->sa_r_he_p.push_back(-999999);
1194  muonData->sa_phi_he_p.push_back(-999999);
1195  }
1196 
1197  // Extrapolation to HE-
1198  tsos = surfExtrapTrkSam(sa_mu, -479);
1199  if (tsos.isValid()) {
1200  double xx = tsos.globalPosition().x();
1201  double yy = tsos.globalPosition().y();
1202  double rr = sqrt(xx * xx + yy * yy);
1203  muonData->sa_r_he_n.push_back(rr);
1204  double cosphi = xx / rr;
1205  if (yy >= 0)
1206  muonData->sa_phi_he_n.push_back(acos(cosphi));
1207  else
1208  muonData->sa_phi_he_n.push_back(2 * pig - acos(cosphi));
1209  } else {
1210  muonData->sa_r_he_n.push_back(-999999);
1211  muonData->sa_phi_he_n.push_back(-999999);
1212  }
1213 
1214  // Several track quantities
1215  muonData->sa_normchi2.push_back(sa_mu->normalizedChi2());
1216  muonData->sa_validhits.push_back(sa_mu->numberOfValidHits());
1217  muonData->sa_ch.push_back(sa_mu->charge());
1218  muonData->sa_pt.push_back(sa_mu->pt());
1219  muonData->sa_p.push_back(sa_mu->p());
1220  muonData->sa_eta.push_back(sa_mu->eta());
1221  muonData->sa_phi.push_back(sa_mu->phi());
1222  muonData->sa_outer_pt.push_back(sqrt(sa_mu->outerMomentum().Perp2()));
1223  muonData->sa_inner_pt.push_back(sqrt(sa_mu->innerMomentum().Perp2()));
1224  muonData->sa_outer_eta.push_back(sa_mu->outerMomentum().Eta());
1225  muonData->sa_inner_eta.push_back(sa_mu->innerMomentum().Eta());
1226  muonData->sa_outer_phi.push_back(sa_mu->outerMomentum().Phi());
1227  muonData->sa_inner_phi.push_back(sa_mu->innerMomentum().Phi());
1228  muonData->sa_outer_x.push_back(sa_mu->outerPosition().x());
1229  muonData->sa_outer_y.push_back(sa_mu->outerPosition().y());
1230  muonData->sa_outer_z.push_back(sa_mu->outerPosition().z());
1231  muonData->sa_inner_x.push_back(sa_mu->innerPosition().x());
1232  muonData->sa_inner_y.push_back(sa_mu->innerPosition().y());
1233  muonData->sa_inner_z.push_back(sa_mu->innerPosition().z());
1234 
1235  } // end of IF IS STANDALONE
1236 
1237  } // end of if type 012 found
1238  } // end of muon loop
1239 
1240  tree_->Fill();
1241 }

References pwdgSkimBPark_cfi::beamSpot, l1GtPatternGenerator_cfi::bx, L1Analysis::L1AnalysisRecoRpcHitDataFormat::bx, L1Analysis::L1AnalysisRecoMuonDataFormat::calo_energy, L1Analysis::L1AnalysisRecoMuonDataFormat::calo_energy3x3, L1Analysis::L1AnalysisRecoMuonDataFormat::ch, CSCGeometry::chamber(), L1Analysis::L1AnalysisRecoRpcHitDataFormat::clusterSize, MuonSubdetId::CSC, cscGeom, cscRecHitPSet::cscSegment, 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(), L1Analysis::L1AnalysisRecoMuonDataFormat::eta, edm::EventSetup::get(), get, 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, triggerObjects_cff::id, iEvent, 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, dqm-mbProfile::log, PV3DBase< T, PVType, FrameType >::mag(), match_trigger(), maxMuon_, maxRpcHit_, DetId::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::L1AnalysisRecoRpcHitDataFormat::ring, RPCDetId::ring(), relativeConstraints::ring, L1Analysis::L1AnalysisRecoRpcHitDataFormat::roll, RPCGeometry::roll(), RPCDetId::roll(), MuonSubdetId::RPC, rpcGeom, rpcHit, rpcHitData, rpcHitTag_, dtTriggerPhase2PrimitiveDigis_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(), GeomDet::surface(), 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_, PDWG_DiPhoton_SD_cff::triggerEvent, triggerIndices_, triggerMatching_, triggerMaxDeltaR_, triggerMatchMonitor_cfi::triggerObjects, triggerProcessLabel_, triggerResults, triggerSummaryLabel_, TRSA_MUON, L1Analysis::L1AnalysisRecoMuonDataFormat::type, L1Analysis::L1AnalysisRecoMuonDataFormat::validhits, bphysicsOniaDQM_cfi::vertex, PV3DBase< T, PVType, FrameType >::x(), L1Analysis::L1AnalysisRecoRpcHitDataFormat::xLoc, geometryCSVtoXML::xx, PV3DBase< T, PVType, FrameType >::y(), geometryCSVtoXML::yy, PV3DBase< T, PVType, FrameType >::z(), L1Analysis::L1AnalysisRecoMuonDataFormat::z_hb, and geometryCSVtoXML::zz.

◆ beginJob()

void L1MuonRecoTreeProducer::beginJob ( void  )
overrideprivatevirtual

Reimplemented from edm::EDAnalyzer.

Definition at line 1337 of file L1MuonRecoTreeProducer.cc.

1337 {}

◆ beginRun()

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

Reimplemented from edm::EDAnalyzer.

Definition at line 1282 of file L1MuonRecoTreeProducer.cc.

1282  {
1283  // Prepare for trigger matching for each new run:
1284  // Look up triggetIndices in the HLT config for the different paths
1285  if (triggerMatching_) {
1286  bool changed = true;
1287  if (!hltConfig_.init(run, eventSetup, triggerProcessLabel_, changed)) {
1288  // if you can't initialize hlt configuration, crash!
1289  std::cout << "Error: didn't find process" << triggerProcessLabel_ << std::endl;
1290  assert(false);
1291  }
1292 
1293  bool enableWildcard = true;
1294  for (size_t iTrig = 0; iTrig < triggerNames_.size(); ++iTrig) {
1295  // prepare for regular expression (with wildcards) functionality:
1296  TString tNameTmp = TString(triggerNames_[iTrig]);
1297  TRegexp tNamePattern = TRegexp(tNameTmp, enableWildcard);
1298  int tIndex = -1;
1299  // find the trigger index:
1300  for (unsigned ipath = 0; ipath < hltConfig_.size(); ++ipath) {
1301  // use TString since it provides reg exp functionality:
1302  TString tmpName = TString(hltConfig_.triggerName(ipath));
1303  if (tmpName.Contains(tNamePattern)) {
1304  tIndex = int(ipath);
1305  triggerIndices_.push_back(tIndex);
1306  }
1307  }
1308  if (tIndex < 0) { // if can't find trigger path at all, give warning:
1309  std::cout << "Warning: Could not find trigger" << triggerNames_[iTrig] << std::endl;
1310  //assert(false);
1311  }
1312  } // end for triggerNames
1313  for (size_t iTrig = 0; iTrig < isoTriggerNames_.size(); ++iTrig) {
1314  // prepare for regular expression functionality:
1315  TString tNameTmp = TString(isoTriggerNames_[iTrig]);
1316  TRegexp tNamePattern = TRegexp(tNameTmp, enableWildcard);
1317  int tIndex = -1;
1318  // find the trigger index:
1319  for (unsigned ipath = 0; ipath < hltConfig_.size(); ++ipath) {
1320  // use TString since it provides reg exp functionality:
1321  TString tmpName = TString(hltConfig_.triggerName(ipath));
1322  if (tmpName.Contains(tNamePattern)) {
1323  tIndex = int(ipath);
1324  isoTriggerIndices_.push_back(tIndex);
1325  }
1326  }
1327  if (tIndex < 0) { // if can't find trigger path at all, give warning:
1328  std::cout << "Warning: Could not find trigger" << isoTriggerNames_[iTrig] << std::endl;
1329  //assert(false);
1330  }
1331  } // end for isoTriggerNames
1332  } // end if (triggerMatching_)
1333 }

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

◆ cylExtrapTrkSam()

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

Definition at line 1244 of file L1MuonRecoTreeProducer.cc.

1244  {
1245  Cylinder::PositionType pos(0, 0, 0);
1246  Cylinder::RotationType rot;
1247  Cylinder::CylinderPointer myCylinder = Cylinder::build(pos, rot, rho);
1248 
1249  FreeTrajectoryState recoStart = freeTrajStateMuon(track);
1250  TrajectoryStateOnSurface recoProp;
1251  recoProp = propagatorAlong->propagate(recoStart, *myCylinder);
1252  if (!recoProp.isValid()) {
1253  recoProp = propagatorOpposite->propagate(recoStart, *myCylinder);
1254  }
1255  return recoProp;
1256 }

References Cylinder::build(), freeTrajStateMuon(), TrajectoryStateOnSurface::isValid(), Propagator::propagate(), propagatorAlong, propagatorOpposite, rho, makeMuonMisalignmentScenario::rot, and HLT_FULL_cff::track.

Referenced by analyze().

◆ empty_global()

void L1MuonRecoTreeProducer::empty_global ( )

Definition at line 263 of file L1MuonRecoTreeProducer.cc.

263  {
264  muonData->ch.push_back(-999999);
265  muonData->pt.push_back(-999999);
266  muonData->p.push_back(-999999);
267  muonData->eta.push_back(-999999);
268  muonData->phi.push_back(-999999);
269  muonData->normchi2.push_back(-999999);
270  muonData->validhits.push_back(-999999);
271  muonData->numberOfMatchedStations.push_back(-999999);
272  muonData->numberOfValidMuonHits.push_back(-999999);
273  muonData->imp_point_x.push_back(-999999);
274  muonData->imp_point_y.push_back(-999999);
275  muonData->imp_point_z.push_back(-999999);
276  muonData->imp_point_p.push_back(-999999);
277  muonData->imp_point_pt.push_back(-999999);
278  muonData->phi_hb.push_back(-999999);
279  muonData->z_hb.push_back(-999999);
280  muonData->r_he_p.push_back(-999999);
281  muonData->phi_he_p.push_back(-999999);
282  muonData->r_he_n.push_back(-999999);
283  muonData->phi_he_n.push_back(-999999);
284  muonData->calo_energy.push_back(-999999);
285  muonData->calo_energy3x3.push_back(-999999);
286  muonData->ecal_time.push_back(-999999);
287  muonData->ecal_terr.push_back(-999999);
288  muonData->hcal_time.push_back(-999999);
289  muonData->hcal_terr.push_back(-999999);
290  muonData->time_dir.push_back(-999999);
291  muonData->time_inout.push_back(-999999);
292  muonData->time_inout_err.push_back(-999999);
293  muonData->time_outin.push_back(-999999);
294  muonData->time_outin_err.push_back(-999999);
295 
296  muonData->hlt_isomu.push_back(-999999);
297  muonData->hlt_mu.push_back(-999999);
298  muonData->hlt_isoDeltaR.push_back(-999999);
299  muonData->hlt_deltaR.push_back(-999999);
300 }

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().

◆ empty_hlt()

void L1MuonRecoTreeProducer::empty_hlt ( )

Definition at line 388 of file L1MuonRecoTreeProducer.cc.

388  {
389  muonData->hlt_isomu.push_back(-999999);
390  muonData->hlt_mu.push_back(-999999);
391  muonData->hlt_isoDeltaR.push_back(-9999999);
392  muonData->hlt_deltaR.push_back(-999999);
393 }

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

Referenced by analyze().

◆ empty_standalone()

void L1MuonRecoTreeProducer::empty_standalone ( )

Definition at line 333 of file L1MuonRecoTreeProducer.cc.

333  {
334  muonData->sa_imp_point_x.push_back(-999999);
335  muonData->sa_imp_point_y.push_back(-999999);
336  muonData->sa_imp_point_z.push_back(-999999);
337  muonData->sa_imp_point_p.push_back(-999999);
338  muonData->sa_imp_point_pt.push_back(-999999);
339  muonData->sa_z_mb2.push_back(-999999);
340  muonData->sa_phi_mb2.push_back(-999999);
341  muonData->sa_pseta.push_back(-999999);
342 
343  muonData->sa_z_hb.push_back(-999999);
344  muonData->sa_phi_hb.push_back(-999999);
345 
346  muonData->sa_r_he_p.push_back(-999999);
347  muonData->sa_phi_he_p.push_back(-999999);
348  muonData->sa_r_he_n.push_back(-999999);
349  muonData->sa_phi_he_n.push_back(-999999);
350 
351  muonData->sa_normchi2.push_back(-999999);
352  muonData->sa_validhits.push_back(-999999);
353  muonData->sa_ch.push_back(-999999);
354  muonData->sa_pt.push_back(-999999);
355  muonData->sa_p.push_back(-999999);
356  muonData->sa_eta.push_back(-999999);
357  muonData->sa_phi.push_back(-999999);
358  muonData->sa_outer_pt.push_back(-999999);
359  muonData->sa_inner_pt.push_back(-999999);
360  muonData->sa_outer_eta.push_back(-999999);
361  muonData->sa_inner_eta.push_back(-999999);
362  muonData->sa_outer_phi.push_back(-999999);
363  muonData->sa_inner_phi.push_back(-999999);
364  muonData->sa_outer_x.push_back(-999999);
365  muonData->sa_outer_y.push_back(-999999);
366  muonData->sa_outer_z.push_back(-999999);
367  muonData->sa_inner_x.push_back(-999999);
368  muonData->sa_inner_y.push_back(-999999);
369  muonData->sa_inner_z.push_back(-999999);
370 
371  muonData->sa_r_me2_p.push_back(-999999);
372  muonData->sa_phi_me2_p.push_back(-999999);
373 
374  muonData->sa_r_me2_n.push_back(-999999);
375  muonData->sa_phi_me2_n.push_back(-999999);
376 
377  muonData->sa_z_mb1.push_back(-999999);
378  muonData->sa_phi_mb1.push_back(-999999);
379  muonData->sa_r_me1_p.push_back(-999999);
380  muonData->sa_phi_me1_p.push_back(-999999);
381  muonData->sa_r_me1_n.push_back(-999999);
382  muonData->sa_phi_me1_n.push_back(-999999);
383 
384  muonData->sa_nChambers.push_back(-999);
385  muonData->sa_nMatches.push_back(-999);
386 }

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().

◆ empty_tracker()

void L1MuonRecoTreeProducer::empty_tracker ( )

Definition at line 302 of file L1MuonRecoTreeProducer.cc.

302  {
303  muonData->tr_ch.push_back(-999999);
304  muonData->tr_pt.push_back(-999999);
305  muonData->tr_p.push_back(-999999);
306  muonData->tr_eta.push_back(-999999);
307  muonData->tr_phi.push_back(-999999);
308  muonData->tr_normchi2.push_back(-999999);
309  muonData->tr_validhits.push_back(-999999);
310  muonData->tr_validpixhits.push_back(-999999);
311  muonData->tr_d0.push_back(-999999);
312  muonData->tr_imp_point_x.push_back(-999999);
313  muonData->tr_imp_point_y.push_back(-999999);
314  muonData->tr_imp_point_z.push_back(-999999);
315  muonData->tr_imp_point_p.push_back(-999999);
316  muonData->tr_imp_point_pt.push_back(-999999);
317 
318  muonData->tr_z_mb2.push_back(-999999);
319  muonData->tr_phi_mb2.push_back(-999999);
320  muonData->tr_r_me2_p.push_back(-999999);
321  muonData->tr_phi_me2_p.push_back(-999999);
322  muonData->tr_r_me2_n.push_back(-999999);
323  muonData->tr_phi_me2_n.push_back(-999999);
324 
325  muonData->tr_z_mb1.push_back(-999999);
326  muonData->tr_phi_mb1.push_back(-999999);
327  muonData->tr_r_me1_p.push_back(-999999);
328  muonData->tr_phi_me1_p.push_back(-999999);
329  muonData->tr_r_me1_n.push_back(-999999);
330  muonData->tr_phi_me1_n.push_back(-999999);
331 }

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().

◆ endJob()

void L1MuonRecoTreeProducer::endJob ( void  )
overrideprivatevirtual

Reimplemented from edm::EDAnalyzer.

Definition at line 1340 of file L1MuonRecoTreeProducer.cc.

1340 { delete muon; }

References muon.

◆ endRun()

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

Reimplemented from edm::EDAnalyzer.

Definition at line 1335 of file L1MuonRecoTreeProducer.cc.

1335 {}

◆ freeTrajStateMuon()

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

Definition at line 1273 of file L1MuonRecoTreeProducer.cc.

1273  {
1274  GlobalPoint innerPoint(track->innerPosition().x(), track->innerPosition().y(), track->innerPosition().z());
1275  GlobalVector innerVec(track->innerMomentum().x(), track->innerMomentum().y(), track->innerMomentum().z());
1276 
1277  FreeTrajectoryState recoStart(innerPoint, innerVec, track->charge(), &*theBField);
1278 
1279  return recoStart;
1280 }

References theBField, and HLT_FULL_cff::track.

Referenced by cylExtrapTrkSam(), and surfExtrapTrkSam().

◆ match_trigger()

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.

229  {
230  double matchDeltaR = 9999;
231 
232  for (size_t iTrigIndex = 0; iTrigIndex < trigIndices.size(); ++iTrigIndex) {
233  int triggerIndex = trigIndices[iTrigIndex];
234  const std::vector<std::string> moduleLabels(hltConfig_.moduleLabels(triggerIndex));
235  // find index of the last module:
236  const unsigned moduleIndex = hltConfig_.size(triggerIndex) - 2;
237  // find index of HLT trigger name:
238  const unsigned hltFilterIndex =
239  triggerEvent->filterIndex(edm::InputTag(moduleLabels[moduleIndex], "", triggerProcessLabel_));
240 
241  if (hltFilterIndex < triggerEvent->sizeFilters()) {
242  const trigger::Keys triggerKeys(triggerEvent->filterKeys(hltFilterIndex));
243  const trigger::Vids triggerVids(triggerEvent->filterIds(hltFilterIndex));
244 
245  const unsigned nTriggers = triggerVids.size();
246  for (size_t iTrig = 0; iTrig < nTriggers; ++iTrig) {
247  // loop over all trigger objects:
248  const trigger::TriggerObject trigObject = trigObjs[triggerKeys[iTrig]];
249 
250  double dRtmp = deltaR(mu, trigObject);
251 
252  if (dRtmp < matchDeltaR) {
253  matchDeltaR = dRtmp;
254  }
255 
256  } // loop over different trigger objects
257  } // if trigger is in event (should apply hltFilter with used trigger...)
258  } // loop over muon candidates
259 
260  return matchDeltaR;
261 }

References PbPb_ZMuSkimMuonDPG_cff::deltaR, hltConfig_, HLTConfigProvider::moduleLabels(), amptDefaultParameters_cff::mu, HLTConfigProvider::size(), PDWG_DiPhoton_SD_cff::triggerEvent, triggerProcessLabel_, and TriggerAnalyzer::trigObjs.

Referenced by analyze().

◆ surfExtrapTrkSam()

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

Definition at line 1259 of file L1MuonRecoTreeProducer.cc.

1259  {
1260  Plane::PositionType pos(0, 0, z);
1261  Plane::RotationType rot;
1262  Plane::PlanePointer myPlane = Plane::build(pos, rot);
1263 
1264  FreeTrajectoryState recoStart = freeTrajStateMuon(track);
1265  TrajectoryStateOnSurface recoProp;
1266  recoProp = propagatorAlong->propagate(recoStart, *myPlane);
1267  if (!recoProp.isValid()) {
1268  recoProp = propagatorOpposite->propagate(recoStart, *myPlane);
1269  }
1270  return recoProp;
1271 }

References Plane::build(), freeTrajStateMuon(), TrajectoryStateOnSurface::isValid(), Propagator::propagate(), propagatorAlong, propagatorOpposite, makeMuonMisalignmentScenario::rot, HLT_FULL_cff::track, and z.

Referenced by analyze().

Member Data Documentation

◆ cscGeom

edm::ESHandle<CSCGeometry> L1MuonRecoTreeProducer::cscGeom
private

Definition at line 154 of file L1MuonRecoTreeProducer.cc.

Referenced by analyze().

◆ dtGeom

edm::ESHandle<DTGeometry> L1MuonRecoTreeProducer::dtGeom
private

Definition at line 158 of file L1MuonRecoTreeProducer.cc.

Referenced by analyze().

◆ fs_

edm::Service<TFileService> L1MuonRecoTreeProducer::fs_
private

Definition at line 176 of file L1MuonRecoTreeProducer.cc.

Referenced by L1MuonRecoTreeProducer().

◆ hltConfig_

HLTConfigProvider L1MuonRecoTreeProducer::hltConfig_
private

Definition at line 149 of file L1MuonRecoTreeProducer.cc.

Referenced by beginRun(), and match_trigger().

◆ isoTriggerIndices_

std::vector<int> L1MuonRecoTreeProducer::isoTriggerIndices_
private

Definition at line 146 of file L1MuonRecoTreeProducer.cc.

Referenced by analyze(), and beginRun().

◆ isoTriggerNames_

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

Definition at line 143 of file L1MuonRecoTreeProducer.cc.

Referenced by beginRun(), and L1MuonRecoTreeProducer().

◆ maxMuon_

unsigned L1MuonRecoTreeProducer::maxMuon_
private

Definition at line 133 of file L1MuonRecoTreeProducer.cc.

Referenced by analyze(), and L1MuonRecoTreeProducer().

◆ maxRpcHit_

unsigned L1MuonRecoTreeProducer::maxRpcHit_
private

Definition at line 134 of file L1MuonRecoTreeProducer.cc.

Referenced by analyze(), and L1MuonRecoTreeProducer().

◆ muon

L1Analysis::L1AnalysisRecoMuon* L1MuonRecoTreeProducer::muon

Definition at line 126 of file L1MuonRecoTreeProducer.cc.

Referenced by endJob().

◆ muonData

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().

◆ muonTag_

edm::InputTag L1MuonRecoTreeProducer::muonTag_
private

Definition at line 184 of file L1MuonRecoTreeProducer.cc.

Referenced by analyze(), and L1MuonRecoTreeProducer().

◆ propagatorAlong

edm::ESHandle<Propagator> L1MuonRecoTreeProducer::propagatorAlong
private

Definition at line 170 of file L1MuonRecoTreeProducer.cc.

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

◆ propagatorOpposite

edm::ESHandle<Propagator> L1MuonRecoTreeProducer::propagatorOpposite
private

Definition at line 171 of file L1MuonRecoTreeProducer.cc.

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

◆ rpcGeom

edm::ESHandle<RPCGeometry> L1MuonRecoTreeProducer::rpcGeom
private

Definition at line 161 of file L1MuonRecoTreeProducer.cc.

Referenced by analyze().

◆ rpcHit

L1Analysis::L1AnalysisRecoRpcHit* L1MuonRecoTreeProducer::rpcHit

Definition at line 129 of file L1MuonRecoTreeProducer.cc.

Referenced by analyze(), and L1MuonRecoTreeProducer().

◆ rpcHitData

L1Analysis::L1AnalysisRecoRpcHitDataFormat* L1MuonRecoTreeProducer::rpcHitData

Definition at line 130 of file L1MuonRecoTreeProducer.cc.

Referenced by analyze(), and L1MuonRecoTreeProducer().

◆ rpcHitTag_

edm::InputTag L1MuonRecoTreeProducer::rpcHitTag_
private

Definition at line 185 of file L1MuonRecoTreeProducer.cc.

Referenced by analyze(), and L1MuonRecoTreeProducer().

◆ runOnPostLS1_

bool L1MuonRecoTreeProducer::runOnPostLS1_
private

Definition at line 181 of file L1MuonRecoTreeProducer.cc.

Referenced by analyze(), and L1MuonRecoTreeProducer().

◆ theBField

edm::ESHandle<MagneticField> L1MuonRecoTreeProducer::theBField
private

Definition at line 164 of file L1MuonRecoTreeProducer.cc.

Referenced by analyze(), and freeTrajStateMuon().

◆ theTrackingGeometry

edm::ESHandle<GlobalTrackingGeometry> L1MuonRecoTreeProducer::theTrackingGeometry
private

Definition at line 167 of file L1MuonRecoTreeProducer.cc.

Referenced by analyze().

◆ tree_

TTree* L1MuonRecoTreeProducer::tree_
private

Definition at line 179 of file L1MuonRecoTreeProducer.cc.

Referenced by analyze(), and L1MuonRecoTreeProducer().

◆ triggerIndices_

std::vector<int> L1MuonRecoTreeProducer::triggerIndices_
private

Definition at line 147 of file L1MuonRecoTreeProducer.cc.

Referenced by analyze(), and beginRun().

◆ triggerMatching_

bool L1MuonRecoTreeProducer::triggerMatching_
private

Definition at line 140 of file L1MuonRecoTreeProducer.cc.

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

◆ triggerMaxDeltaR_

double L1MuonRecoTreeProducer::triggerMaxDeltaR_
private

Definition at line 148 of file L1MuonRecoTreeProducer.cc.

Referenced by analyze(), and L1MuonRecoTreeProducer().

◆ triggerNames_

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

Definition at line 144 of file L1MuonRecoTreeProducer.cc.

Referenced by beginRun(), and L1MuonRecoTreeProducer().

◆ triggerProcessLabel_

std::string L1MuonRecoTreeProducer::triggerProcessLabel_
private

Definition at line 142 of file L1MuonRecoTreeProducer.cc.

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

◆ triggerSummaryLabel_

edm::InputTag L1MuonRecoTreeProducer::triggerSummaryLabel_
private

Definition at line 141 of file L1MuonRecoTreeProducer.cc.

Referenced by analyze(), and L1MuonRecoTreeProducer().

L1Analysis::L1AnalysisRecoMuonDataFormat::sa_nChambers
std::vector< int > sa_nChambers
Definition: L1AnalysisRecoMuonDataFormat.h:266
Vector3DBase
Definition: Vector3DBase.h:8
L1Analysis::L1AnalysisRecoMuonDataFormat::sa_outer_y
std::vector< double > sa_outer_y
Definition: L1AnalysisRecoMuonDataFormat.h:225
TkRotation< float >
reco::MuonEnergy
Definition: MuonEnergy.h:17
reco::MuonEnergy::ecal_time
float ecal_time
Calorimeter timing.
Definition: MuonEnergy.h:47
L1Analysis::L1AnalysisRecoMuonDataFormat::tr_phi_me1_n
std::vector< double > tr_phi_me1_n
Definition: L1AnalysisRecoMuonDataFormat.h:205
MuonSubdetId::CSC
static constexpr int CSC
Definition: MuonSubdetId.h:12
L1MuonRecoTreeProducer::isoTriggerNames_
std::vector< std::string > isoTriggerNames_
Definition: L1MuonRecoTreeProducer.cc:143
geometryCSVtoXML.zz
zz
Definition: geometryCSVtoXML.py:19
L1MuonRecoTreeProducer::runOnPostLS1_
bool runOnPostLS1_
Definition: L1MuonRecoTreeProducer.cc:181
pwdgSkimBPark_cfi.beamSpot
beamSpot
Definition: pwdgSkimBPark_cfi.py:5
HLT_FULL_cff.track
track
Definition: HLT_FULL_cff.py:11724
L1Analysis::L1AnalysisRecoMuonDataFormat::sa_imp_point_y
std::vector< double > sa_imp_point_y
Definition: L1AnalysisRecoMuonDataFormat.h:231
L1Analysis::L1AnalysisRecoMuonDataFormat::hcal_terr
std::vector< double > hcal_terr
Definition: L1AnalysisRecoMuonDataFormat.h:258
L1Analysis::L1AnalysisRecoMuonDataFormat::time_inout
std::vector< double > time_inout
Definition: L1AnalysisRecoMuonDataFormat.h:261
RPCDetId::station
int station() const
Definition: RPCDetId.h:78
RPCDetId::region
int region() const
Region id: 0 for Barrel, +/-1 For +/- Endcap.
Definition: RPCDetId.h:53
L1MuonRecoTreeProducer::propagatorAlong
edm::ESHandle< Propagator > propagatorAlong
Definition: L1MuonRecoTreeProducer.cc:170
L1Analysis::L1AnalysisRecoMuonDataFormat::imp_point_x
std::vector< double > imp_point_x
Definition: L1AnalysisRecoMuonDataFormat.h:165
L1Analysis::L1AnalysisRecoMuonDataFormat::z_hb
std::vector< double > z_hb
Definition: L1AnalysisRecoMuonDataFormat.h:171
RPCDetId::subsector
int subsector() const
SubSector id : some sectors are divided along the phi direction in subsectors (from 1 to 4 in Barrel,...
Definition: RPCDetId.h:88
L1Analysis::L1AnalysisRecoMuonDataFormat::tr_imp_point_z
std::vector< double > tr_imp_point_z
Definition: L1AnalysisRecoMuonDataFormat.h:189
L1Analysis::L1AnalysisRecoMuonDataFormat::imp_point_z
std::vector< double > imp_point_z
Definition: L1AnalysisRecoMuonDataFormat.h:167
triggerMatchMonitor_cfi.triggerObjects
triggerObjects
Definition: triggerMatchMonitor_cfi.py:15
trigger::Vids
std::vector< int > Vids
Definition: TriggerTypeDefs.h:21
L1Analysis::L1AnalysisRecoRpcHitDataFormat::bx
std::vector< int > bx
Definition: L1AnalysisRecoRpcHitDataFormat.h:43
L1Analysis::L1AnalysisRecoMuonDataFormat::calo_energy
std::vector< double > calo_energy
Definition: L1AnalysisRecoMuonDataFormat.h:253
L1Analysis::L1AnalysisRecoMuonDataFormat::sa_pt
std::vector< double > sa_pt
Definition: L1AnalysisRecoMuonDataFormat.h:214
L1Analysis::L1AnalysisRecoMuonDataFormat::sa_phi_he_n
std::vector< double > sa_phi_he_n
Definition: L1AnalysisRecoMuonDataFormat.h:240
muon
Definition: MuonCocktails.h:17
L1MuonRecoTreeProducer::isoTriggerIndices_
std::vector< int > isoTriggerIndices_
Definition: L1MuonRecoTreeProducer.cc:146
L1MuonRecoTreeProducer::TRSA_MUON
Definition: L1MuonRecoTreeProducer.cc:151
PV3DBase::x
T x() const
Definition: PV3DBase.h:59
L1MuonRecoTreeProducer::triggerMaxDeltaR_
double triggerMaxDeltaR_
Definition: L1MuonRecoTreeProducer.cc:148
L1Analysis::L1AnalysisRecoMuonDataFormat::sa_phi_me2_n
std::vector< double > sa_phi_me2_n
Definition: L1AnalysisRecoMuonDataFormat.h:244
amptDefaultParameters_cff.mu
mu
Definition: amptDefaultParameters_cff.py:16
L1Analysis::L1AnalysisRecoMuonDataFormat::rchCSCtype
std::vector< int > rchCSCtype
Definition: L1AnalysisRecoMuonDataFormat.h:270
L1Analysis::L1AnalysisRecoMuonDataFormat::sa_imp_point_z
std::vector< double > sa_imp_point_z
Definition: L1AnalysisRecoMuonDataFormat.h:232
TriggerAnalyzer.trigObjs
trigObjs
Definition: TriggerAnalyzer.py:89
findQualityFiles.rr
string rr
Definition: findQualityFiles.py:185
relativeConstraints.station
station
Definition: relativeConstraints.py:67
L1Analysis::L1AnalysisRecoMuonDataFormat::tr_phi_me1_p
std::vector< double > tr_phi_me1_p
Definition: L1AnalysisRecoMuonDataFormat.h:203
L1Analysis::L1AnalysisRecoMuon
Definition: L1AnalysisRecoMuon.h:15
L1MuonRecoTreeProducer::match_trigger
double match_trigger(std::vector< int > &trigIndices, const trigger::TriggerObjectCollection &trigObjs, edm::Handle< trigger::TriggerEvent > &triggerEvent, const reco::Muon &mu)
Definition: L1MuonRecoTreeProducer.cc:226
L1Analysis::L1AnalysisRecoMuonDataFormat::rchPhi
std::vector< double > rchPhi
Definition: L1AnalysisRecoMuonDataFormat.h:271
DetId::det
constexpr Detector det() const
get the detector field from this detid
Definition: DetId.h:46
L1Analysis::L1AnalysisRecoRpcHitDataFormat::sector
std::vector< int > sector
Definition: L1AnalysisRecoRpcHitDataFormat.h:49
L1Analysis::L1AnalysisRecoMuonDataFormat::hlt_isoDeltaR
std::vector< double > hlt_isoDeltaR
Definition: L1AnalysisRecoMuonDataFormat.h:277
RPCDetId
Definition: RPCDetId.h:16
TrajectoryStateOnSurface::globalPosition
GlobalPoint globalPosition() const
Definition: TrajectoryStateOnSurface.h:65
gather_cfg.cout
cout
Definition: gather_cfg.py:144
L1Analysis::L1AnalysisRecoMuonDataFormat::sa_ch
std::vector< double > sa_ch
Definition: L1AnalysisRecoMuonDataFormat.h:213
pos
Definition: PixelAliasList.h:18
L1Analysis::L1AnalysisRecoMuonDataFormat::r_he_p
std::vector< double > r_he_p
Definition: L1AnalysisRecoMuonDataFormat.h:172
L1Analysis::L1AnalysisRecoMuonDataFormat::tr_phi_me2_p
std::vector< double > tr_phi_me2_p
Definition: L1AnalysisRecoMuonDataFormat.h:196
l1GtPatternGenerator_cfi.bx
bx
Definition: l1GtPatternGenerator_cfi.py:18
L1Analysis::L1AnalysisRecoMuonDataFormat::sa_outer_eta
std::vector< double > sa_outer_eta
Definition: L1AnalysisRecoMuonDataFormat.h:220
triggerResults
static const std::string triggerResults
Definition: EdmProvDump.cc:45
L1Analysis::L1AnalysisRecoMuonDataFormat::sa_imp_point_x
std::vector< double > sa_imp_point_x
Definition: L1AnalysisRecoMuonDataFormat.h:230
L1MuonRecoTreeProducer::rpcHitTag_
edm::InputTag rpcHitTag_
Definition: L1MuonRecoTreeProducer.cc:185
cms::cuda::assert
assert(be >=bs)
L1Analysis::L1AnalysisRecoMuonDataFormat::numberOfValidMuonHits
std::vector< double > numberOfValidMuonHits
Definition: L1AnalysisRecoMuonDataFormat.h:163
GlobalTrackingGeometryRecord
Definition: GlobalTrackingGeometryRecord.h:17
RPCGeometry::roll
const RPCRoll * roll(RPCDetId id) const
Return a roll given its id.
Definition: RPCGeometry.cc:50
L1MuonRecoTreeProducer::muonData
L1Analysis::L1AnalysisRecoMuonDataFormat * muonData
Definition: L1MuonRecoTreeProducer.cc:127
L1Analysis::L1AnalysisRecoRpcHitDataFormat::region
std::vector< int > region
Definition: L1AnalysisRecoRpcHitDataFormat.h:40
reco::MuonTime::timeAtIpOutInErr
float timeAtIpOutInErr
Definition: MuonTime.h:17
edm::ParameterSet::getUntrackedParameter
T getUntrackedParameter(std::string const &, T const &) const
ReferenceCountingPointer< Cylinder >
L1Analysis::L1AnalysisRecoMuonDataFormat::sa_phi_mb1
std::vector< double > sa_phi_mb1
Definition: L1AnalysisRecoMuonDataFormat.h:247
L1Analysis::L1AnalysisRecoMuonDataFormat::sa_z_mb1
std::vector< double > sa_z_mb1
Definition: L1AnalysisRecoMuonDataFormat.h:246
reco::MuonEnergy::hcal_time
float hcal_time
Definition: MuonEnergy.h:49
edm::LogInfo
Log< level::Info, false > LogInfo
Definition: MessageLogger.h:125
L1MuonRecoTreeProducer::fs_
edm::Service< TFileService > fs_
Definition: L1MuonRecoTreeProducer.cc:176
L1Analysis::L1AnalysisRecoMuonDataFormat::sa_phi_me2_p
std::vector< double > sa_phi_me2_p
Definition: L1AnalysisRecoMuonDataFormat.h:243
edm::Handle< RPCRecHitCollection >
reco::MuonEnergy::tower
float tower
Definition: MuonEnergy.h:20
L1MuonRecoTreeProducer::tree_
TTree * tree_
Definition: L1MuonRecoTreeProducer.cc:179
L1Analysis::L1AnalysisRecoMuonDataFormat::eta
std::vector< double > eta
Definition: L1AnalysisRecoMuonDataFormat.h:159
L1Analysis::L1AnalysisRecoRpcHitDataFormat::phiGlob
std::vector< float > phiGlob
Definition: L1AnalysisRecoRpcHitDataFormat.h:46
L1Analysis::L1AnalysisRecoMuonDataFormat::time_dir
std::vector< double > time_dir
Definition: L1AnalysisRecoMuonDataFormat.h:260
L1MuonRecoTreeProducer::rpcHitData
L1Analysis::L1AnalysisRecoRpcHitDataFormat * rpcHitData
Definition: L1MuonRecoTreeProducer.cc:130
L1Analysis::L1AnalysisRecoRpcHitDataFormat::clusterSize
std::vector< int > clusterSize
Definition: L1AnalysisRecoRpcHitDataFormat.h:41
edm::Ref< TrackCollection >
L1MuonRecoTreeProducer::SA_MUON
Definition: L1MuonRecoTreeProducer.cc:151
L1MuonRecoTreeProducer::dtGeom
edm::ESHandle< DTGeometry > dtGeom
Definition: L1MuonRecoTreeProducer.cc:158
reco::MuonTime
Definition: MuonTime.h:5
L1Analysis::L1AnalysisRecoMuonDataFormat::tr_imp_point_pt
std::vector< double > tr_imp_point_pt
Definition: L1AnalysisRecoMuonDataFormat.h:191
PV3DBase::z
T z() const
Definition: PV3DBase.h:61
L1Analysis::L1AnalysisRecoMuonDataFormat::sa_outer_pt
std::vector< double > sa_outer_pt
Definition: L1AnalysisRecoMuonDataFormat.h:218
IdealMagneticFieldRecord
Definition: IdealMagneticFieldRecord.h:11
edm::InputTag::label
std::string const & label() const
Definition: InputTag.h:36
HLTConfigProvider::moduleLabels
const std::vector< std::string > & moduleLabels(unsigned int trigger) const
label(s) of module(s) on a trigger path
Definition: HLTConfigProvider.h:76
DetId
Definition: DetId.h:17
GeomDet::surface
const Plane & surface() const
The nominal surface of the GeomDet.
Definition: GeomDet.h:37
TrajectoryStateOnSurface
Definition: TrajectoryStateOnSurface.h:16
L1MuonRecoTreeProducer::maxRpcHit_
unsigned maxRpcHit_
Definition: L1MuonRecoTreeProducer.cc:134
L1Analysis::L1AnalysisRecoMuonDataFormat::sa_inner_z
std::vector< double > sa_inner_z
Definition: L1AnalysisRecoMuonDataFormat.h:229
edm::EventSetup::get
T get() const
Definition: EventSetup.h:87
L1Analysis::L1AnalysisRecoMuonDataFormat::sa_inner_eta
std::vector< double > sa_inner_eta
Definition: L1AnalysisRecoMuonDataFormat.h:221
L1Analysis::L1AnalysisRecoRpcHit
Definition: L1AnalysisRecoRpcHit.h:15
L1Analysis::L1AnalysisRecoMuonDataFormat::r_he_n
std::vector< double > r_he_n
Definition: L1AnalysisRecoMuonDataFormat.h:173
L1Analysis::L1AnalysisRecoMuonDataFormat::sa_outer_z
std::vector< double > sa_outer_z
Definition: L1AnalysisRecoMuonDataFormat.h:226
L1Analysis::L1AnalysisRecoMuonDataFormat::tr_phi_mb1
std::vector< double > tr_phi_mb1
Definition: L1AnalysisRecoMuonDataFormat.h:201
cscRecHitPSet.cscSegment
cscSegment
Definition: cscRecHitPSet.py:10
L1MuonRecoTreeProducer::cscGeom
edm::ESHandle< CSCGeometry > cscGeom
Definition: L1MuonRecoTreeProducer.cc:154
L1Analysis::L1AnalysisRecoMuonDataFormat::tr_z_mb1
std::vector< double > tr_z_mb1
Definition: L1AnalysisRecoMuonDataFormat.h:200
mathSSE::sqrt
T sqrt(T t)
Definition: SSEVec.h:19
L1MuonRecoTreeProducer::triggerIndices_
std::vector< int > triggerIndices_
Definition: L1MuonRecoTreeProducer.cc:147
DDAxes::z
L1Analysis::L1AnalysisRecoMuonDataFormat::ecal_time
std::vector< double > ecal_time
Definition: L1AnalysisRecoMuonDataFormat.h:255
L1Analysis::L1AnalysisRecoMuonDataFormat::sa_p
std::vector< double > sa_p
Definition: L1AnalysisRecoMuonDataFormat.h:215
CSCChamber
Definition: CSCChamber.h:22
L1Analysis::L1AnalysisRecoMuonDataFormat::sa_inner_y
std::vector< double > sa_inner_y
Definition: L1AnalysisRecoMuonDataFormat.h:228
L1Analysis::L1AnalysisRecoRpcHit::getData
L1AnalysisRecoRpcHitDataFormat * getData()
Definition: L1AnalysisRecoRpcHit.h:23
L1Analysis::L1AnalysisRecoMuonDataFormat::sa_phi_me1_n
std::vector< double > sa_phi_me1_n
Definition: L1AnalysisRecoMuonDataFormat.h:251
L1Analysis::L1AnalysisRecoMuonDataFormat::time_outin
std::vector< double > time_outin
Definition: L1AnalysisRecoMuonDataFormat.h:263
L1Analysis::L1AnalysisRecoMuonDataFormat::sa_phi_hb
std::vector< double > sa_phi_hb
Definition: L1AnalysisRecoMuonDataFormat.h:235
L1Analysis::L1AnalysisRecoRpcHitDataFormat::station
std::vector< int > station
Definition: L1AnalysisRecoRpcHitDataFormat.h:48
L1MuonRecoTreeProducer::TR_MUON
Definition: L1MuonRecoTreeProducer.cc:151
GlobalPoint
Global3DPoint GlobalPoint
Definition: GlobalPoint.h:10
Point3DBase< float, LocalTag >
L1Analysis::L1AnalysisRecoMuonDataFormat::sa_r_he_p
std::vector< double > sa_r_he_p
Definition: L1AnalysisRecoMuonDataFormat.h:237
L1MuonRecoTreeProducer::GL_MUON
Definition: L1MuonRecoTreeProducer.cc:151
L1Analysis::L1AnalysisRecoMuonDataFormat::sa_validhits
std::vector< double > sa_validhits
Definition: L1AnalysisRecoMuonDataFormat.h:212
CSCSegment
Definition: CSCSegment.h:21
L1Analysis::L1AnalysisRecoMuonDataFormat::sa_imp_point_p
std::vector< double > sa_imp_point_p
Definition: L1AnalysisRecoMuonDataFormat.h:233
trigger::TriggerObject
Single trigger physics object (e.g., an isolated muon)
Definition: TriggerObject.h:21
L1Analysis::L1AnalysisRecoMuonDataFormat::tr_phi_mb2
std::vector< double > tr_phi_mb2
Definition: L1AnalysisRecoMuonDataFormat.h:194
L1Analysis::L1AnalysisRecoMuonDataFormat::p
std::vector< double > p
Definition: L1AnalysisRecoMuonDataFormat.h:158
PbPb_ZMuSkimMuonDPG_cff.deltaR
deltaR
Definition: PbPb_ZMuSkimMuonDPG_cff.py:63
edm::OwnVector::const_iterator
Definition: OwnVector.h:41
DetId::subdetId
constexpr int subdetId() const
get the contents of the subdetector field (not cast into any detector's numbering enum)
Definition: DetId.h:48
L1Analysis::L1AnalysisRecoRpcHit::Reset
void Reset()
Definition: L1AnalysisRecoRpcHit.h:20
DDAxes::rho
phase1PixelTopology::layer
constexpr std::array< uint8_t, layerIndexSize > layer
Definition: phase1PixelTopology.h:99
L1Analysis::L1AnalysisRecoMuonDataFormat::tr_validpixhits
std::vector< double > tr_validpixhits
Definition: L1AnalysisRecoMuonDataFormat.h:184
L1Analysis::L1AnalysisRecoMuonDataFormat::sa_phi
std::vector< double > sa_phi
Definition: L1AnalysisRecoMuonDataFormat.h:217
HLTConfigProvider::size
unsigned int size() const
number of trigger paths in trigger table
Definition: HLTConfigProvider.h:63
L1MuonRecoTreeProducer::empty_standalone
void empty_standalone()
Definition: L1MuonRecoTreeProducer.cc:333
L1MuonRecoTreeProducer::triggerSummaryLabel_
edm::InputTag triggerSummaryLabel_
Definition: L1MuonRecoTreeProducer.cc:141
Plane::build
static PlanePointer build(Args &&... args)
Definition: Plane.h:33
L1MuonRecoTreeProducer::hltConfig_
HLTConfigProvider hltConfig_
Definition: L1MuonRecoTreeProducer.cc:149
L1MuonRecoTreeProducer::rpcHit
L1Analysis::L1AnalysisRecoRpcHit * rpcHit
Definition: L1MuonRecoTreeProducer.cc:129
L1Analysis::L1AnalysisRecoMuonDataFormat::sa_r_he_n
std::vector< double > sa_r_he_n
Definition: L1AnalysisRecoMuonDataFormat.h:238
L1Analysis::L1AnalysisRecoMuonDataFormat::phi
std::vector< double > phi
Definition: L1AnalysisRecoMuonDataFormat.h:160
runTauDisplay.gp
gp
Definition: runTauDisplay.py:431
RPCDetId::roll
int roll() const
Definition: RPCDetId.h:92
L1MuonRecoTreeProducer::surfExtrapTrkSam
TrajectoryStateOnSurface surfExtrapTrkSam(reco::TrackRef track, double z)
Definition: L1MuonRecoTreeProducer.cc:1259
L1Analysis::L1AnalysisRecoMuonDataFormat::ch
std::vector< double > ch
Definition: L1AnalysisRecoMuonDataFormat.h:156
L1MuonRecoTreeProducer::empty_tracker
void empty_tracker()
Definition: L1MuonRecoTreeProducer.cc:302
L1Analysis::L1AnalysisRecoMuonDataFormat::sa_inner_x
std::vector< double > sa_inner_x
Definition: L1AnalysisRecoMuonDataFormat.h:227
bphysicsOniaDQM_cfi.vertex
vertex
Definition: bphysicsOniaDQM_cfi.py:7
GeomDet::toGlobal
GlobalPoint toGlobal(const Local2DPoint &lp) const
Conversion to the global R.F. from the R.F. of the GeomDet.
Definition: GeomDet.h:49
L1MuonRecoTreeProducer::freeTrajStateMuon
FreeTrajectoryState freeTrajStateMuon(reco::TrackRef track)
Definition: L1MuonRecoTreeProducer.cc:1273
L1Analysis::L1AnalysisRecoMuonDataFormat::nMuons
int nMuons
Definition: L1AnalysisRecoMuonDataFormat.h:150
reco::MuonTime::timeAtIpInOutErr
float timeAtIpInOutErr
Definition: MuonTime.h:14
L1Analysis::L1AnalysisRecoMuonDataFormat::calo_energy3x3
std::vector< double > calo_energy3x3
Definition: L1AnalysisRecoMuonDataFormat.h:254
geometryCSVtoXML.yy
yy
Definition: geometryCSVtoXML.py:19
reco::MuonEnergy::hcal_timeError
float hcal_timeError
Definition: MuonEnergy.h:50
reco::MuonTime::timeAtIpOutIn
float timeAtIpOutIn
b) particle is moving from outside in
Definition: MuonTime.h:16
L1Analysis::L1AnalysisRecoMuonDataFormat::sa_outer_x
std::vector< double > sa_outer_x
Definition: L1AnalysisRecoMuonDataFormat.h:224
L1Analysis::L1AnalysisRecoMuonDataFormat::tr_z_mb2
std::vector< double > tr_z_mb2
Definition: L1AnalysisRecoMuonDataFormat.h:193
type
type
Definition: SiPixelVCal_PayloadInspector.cc:39
L1Analysis::L1AnalysisRecoMuonDataFormat::howmanytypes
std::vector< int > howmanytypes
Definition: L1AnalysisRecoMuonDataFormat.h:154
trigger::TriggerObjectCollection
std::vector< TriggerObject > TriggerObjectCollection
collection of trigger physics objects (e.g., all isolated muons)
Definition: TriggerObject.h:75
edm::RangeMap::const_iterator
C::const_iterator const_iterator
constant access iterator type
Definition: RangeMap.h:43
L1MuonRecoTreeProducer::muonTag_
edm::InputTag muonTag_
Definition: L1MuonRecoTreeProducer.cc:184
CSCDetId
Definition: CSCDetId.h:26
L1MuonRecoTreeProducer::cylExtrapTrkSam
TrajectoryStateOnSurface cylExtrapTrkSam(reco::TrackRef track, double rho)
Definition: L1MuonRecoTreeProducer.cc:1244
L1Analysis::L1AnalysisRecoMuonDataFormat::imp_point_pt
std::vector< double > imp_point_pt
Definition: L1AnalysisRecoMuonDataFormat.h:169
PV3DBase::y
T y() const
Definition: PV3DBase.h:60
funct::tan
Tan< T >::type tan(const T &t)
Definition: Tan.h:22
Cylinder::build
static CylinderPointer build(const PositionType &pos, const RotationType &rot, Scalar radius, Bounds *bounds=nullptr)
Definition: Cylinder.h:45
L1Analysis::L1AnalysisRecoMuonDataFormat::phi_he_n
std::vector< double > phi_he_n
Definition: L1AnalysisRecoMuonDataFormat.h:175
HLT_FULL_cff.region
region
Definition: HLT_FULL_cff.py:88286
createfilelist.int
int
Definition: createfilelist.py:10
L1Analysis::L1AnalysisRecoRpcHitDataFormat::nRpcHits
int nRpcHits
Definition: L1AnalysisRecoRpcHitDataFormat.h:38
iEvent
int iEvent
Definition: GenABIO.cc:224
reco::MuonEnergy::towerS9
float towerS9
total energy in 3x3 tower shape
Definition: MuonEnergy.h:22
L1Analysis::L1AnalysisRecoMuonDataFormat::tr_p
std::vector< double > tr_p
Definition: L1AnalysisRecoMuonDataFormat.h:180
trigger::Keys
std::vector< size_type > Keys
Definition: TriggerTypeDefs.h:19
L1Analysis::L1AnalysisRecoMuonDataFormat::tr_r_me1_n
std::vector< double > tr_r_me1_n
Definition: L1AnalysisRecoMuonDataFormat.h:204
L1Analysis::L1AnalysisRecoMuonDataFormat::hcal_time
std::vector< double > hcal_time
Definition: L1AnalysisRecoMuonDataFormat.h:257
L1Analysis::L1AnalysisRecoMuonDataFormat::sa_inner_phi
std::vector< double > sa_inner_phi
Definition: L1AnalysisRecoMuonDataFormat.h:223
L1Analysis::L1AnalysisRecoMuonDataFormat::tr_d0
std::vector< double > tr_d0
Definition: L1AnalysisRecoMuonDataFormat.h:186
Propagator::propagate
TrajectoryStateOnSurface propagate(STA const &state, SUR const &surface) const
Definition: Propagator.h:50
L1Analysis::L1AnalysisRecoMuonDataFormat::tr_normchi2
std::vector< double > tr_normchi2
Definition: L1AnalysisRecoMuonDataFormat.h:185
L1Analysis::L1AnalysisRecoMuonDataFormat::sa_imp_point_pt
std::vector< double > sa_imp_point_pt
Definition: L1AnalysisRecoMuonDataFormat.h:234
L1Analysis::L1AnalysisRecoMuonDataFormat::sa_r_me2_p
std::vector< double > sa_r_me2_p
Definition: L1AnalysisRecoMuonDataFormat.h:241
get
#define get
PDWG_DiPhoton_SD_cff.triggerEvent
triggerEvent
Definition: PDWG_DiPhoton_SD_cff.py:39
L1MuonRecoTreeProducer::theTrackingGeometry
edm::ESHandle< GlobalTrackingGeometry > theTrackingGeometry
Definition: L1MuonRecoTreeProducer.cc:167
TFileService::make
T * make(const Args &...args) const
make new ROOT object
Definition: TFileService.h:64
AlCaHLTBitMon_QueryRunRegistry.string
string string
Definition: AlCaHLTBitMon_QueryRunRegistry.py:256
L1Analysis::L1AnalysisRecoRpcHitDataFormat::xLoc
std::vector< float > xLoc
Definition: L1AnalysisRecoRpcHitDataFormat.h:45
L1MuonRecoTreeProducer::muon
L1Analysis::L1AnalysisRecoMuon * muon
Definition: L1MuonRecoTreeProducer.cc:126
L1Analysis::L1AnalysisRecoMuonDataFormat::normchi2
std::vector< double > normchi2
Definition: L1AnalysisRecoMuonDataFormat.h:164
L1Analysis::L1AnalysisRecoRpcHitDataFormat::layer
std::vector< int > layer
Definition: L1AnalysisRecoRpcHitDataFormat.h:50
RPCDetId::ring
int ring() const
Definition: RPCDetId.h:59
reco::MuonEnergy::ecal_timeError
float ecal_timeError
Definition: MuonEnergy.h:48
L1Analysis::L1AnalysisRecoRpcHitDataFormat::muonId
std::vector< int > muonId
Definition: L1AnalysisRecoRpcHitDataFormat.h:54
L1Analysis::L1AnalysisRecoMuonDataFormat::rchEta
std::vector< double > rchEta
Definition: L1AnalysisRecoMuonDataFormat.h:272
DetId::rawId
constexpr uint32_t rawId() const
get the raw id
Definition: DetId.h:57
writedatasetfile.run
run
Definition: writedatasetfile.py:27
L1Analysis::L1AnalysisRecoMuonDataFormat::sa_eta
std::vector< double > sa_eta
Definition: L1AnalysisRecoMuonDataFormat.h:216
L1Analysis::L1AnalysisRecoMuonDataFormat::tr_pt
std::vector< double > tr_pt
Definition: L1AnalysisRecoMuonDataFormat.h:179
reco::TransientTrack
Definition: TransientTrack.h:19
L1Analysis::L1AnalysisRecoMuonDataFormat::tr_r_me2_n
std::vector< double > tr_r_me2_n
Definition: L1AnalysisRecoMuonDataFormat.h:197
HLTConfigProvider::init
bool init(const edm::Run &iRun, const edm::EventSetup &iSetup, const std::string &processName, bool &changed)
d'tor
Definition: HLTConfigProvider.cc:36
FreeTrajectoryState
Definition: FreeTrajectoryState.h:27
L1Analysis::L1AnalysisRecoMuonDataFormat::sa_normchi2
std::vector< double > sa_normchi2
Definition: L1AnalysisRecoMuonDataFormat.h:211
L1Analysis::L1AnalysisRecoMuonDataFormat::sa_r_me1_p
std::vector< double > sa_r_me1_p
Definition: L1AnalysisRecoMuonDataFormat.h:248
dtTriggerPhase2PrimitiveDigis_cfi.rpcRecHits
rpcRecHits
Definition: dtTriggerPhase2PrimitiveDigis_cfi.py:46
L1MuonRecoTreeProducer::rpcGeom
edm::ESHandle< RPCGeometry > rpcGeom
Definition: L1MuonRecoTreeProducer.cc:161
MuonSubdetId::RPC
static constexpr int RPC
Definition: MuonSubdetId.h:13
triggerObjects_cff.id
id
Definition: triggerObjects_cff.py:29
L1Analysis::L1AnalysisRecoMuonDataFormat::sa_phi_he_p
std::vector< double > sa_phi_he_p
Definition: L1AnalysisRecoMuonDataFormat.h:239
L1Analysis::L1AnalysisRecoMuonDataFormat::sa_pseta
std::vector< double > sa_pseta
Definition: L1AnalysisRecoMuonDataFormat.h:210
L1Analysis::L1AnalysisRecoMuonDataFormat::sa_r_me2_n
std::vector< double > sa_r_me2_n
Definition: L1AnalysisRecoMuonDataFormat.h:242
L1Analysis::L1AnalysisRecoMuonDataFormat::tr_r_me2_p
std::vector< double > tr_r_me2_p
Definition: L1AnalysisRecoMuonDataFormat.h:195
L1Analysis::L1AnalysisRecoMuonDataFormat::imp_point_p
std::vector< double > imp_point_p
Definition: L1AnalysisRecoMuonDataFormat.h:168
L1Analysis::L1AnalysisRecoMuonDataFormat::sa_nMatches
std::vector< int > sa_nMatches
Definition: L1AnalysisRecoMuonDataFormat.h:267
BoundPlane
RPCDetId::sector
int sector() const
Sector id: the group of chambers at same phi (and increasing r)
Definition: RPCDetId.h:81
L1Analysis::L1AnalysisRecoMuonDataFormat::tr_imp_point_y
std::vector< double > tr_imp_point_y
Definition: L1AnalysisRecoMuonDataFormat.h:188
L1Analysis::L1AnalysisRecoMuonDataFormat::validhits
std::vector< double > validhits
Definition: L1AnalysisRecoMuonDataFormat.h:161
L1Analysis::L1AnalysisRecoMuonDataFormat::sa_outer_phi
std::vector< double > sa_outer_phi
Definition: L1AnalysisRecoMuonDataFormat.h:222
relativeConstraints.ring
ring
Definition: relativeConstraints.py:68
L1Analysis::L1AnalysisRecoMuonDataFormat::tr_imp_point_p
std::vector< double > tr_imp_point_p
Definition: L1AnalysisRecoMuonDataFormat.h:190
makeMuonMisalignmentScenario.rot
rot
Definition: makeMuonMisalignmentScenario.py:322
L1Analysis::L1AnalysisRecoMuonDataFormat::sa_z_mb2
std::vector< double > sa_z_mb2
Definition: L1AnalysisRecoMuonDataFormat.h:209
L1Analysis::L1AnalysisRecoMuonDataFormat::ecal_terr
std::vector< double > ecal_terr
Definition: L1AnalysisRecoMuonDataFormat.h:256
edm::ParameterSet::getParameter
T getParameter(std::string const &) const
Definition: ParameterSet.h:303
L1Analysis::L1AnalysisRecoMuonDataFormat::phi_he_p
std::vector< double > phi_he_p
Definition: L1AnalysisRecoMuonDataFormat.h:174
L1Analysis::L1AnalysisRecoRpcHitDataFormat::strip
std::vector< int > strip
Definition: L1AnalysisRecoRpcHitDataFormat.h:42
dqm-mbProfile.log
log
Definition: dqm-mbProfile.py:17
L1Analysis::L1AnalysisRecoMuonDataFormat::sa_z_hb
std::vector< double > sa_z_hb
Definition: L1AnalysisRecoMuonDataFormat.h:236
Pi
const double Pi
Definition: CosmicMuonParameters.h:18
L1Analysis::L1AnalysisRecoMuonDataFormat::tr_r_me1_p
std::vector< double > tr_r_me1_p
Definition: L1AnalysisRecoMuonDataFormat.h:202
L1MuonRecoTreeProducer::maxMuon_
unsigned maxMuon_
Definition: L1MuonRecoTreeProducer.cc:133
HLTConfigProvider::triggerName
const std::string & triggerName(unsigned int triggerIndex) const
Definition: HLTConfigProvider.h:70
reco::MuonTime::direction
Direction direction() const
direction estimation based on time dispersion
Definition: MuonTime.h:20
L1Analysis::L1AnalysisRecoMuonDataFormat::sa_phi_me1_p
std::vector< double > sa_phi_me1_p
Definition: L1AnalysisRecoMuonDataFormat.h:249
L1Analysis::L1AnalysisRecoMuonDataFormat::sa_phi_mb2
std::vector< double > sa_phi_mb2
Definition: L1AnalysisRecoMuonDataFormat.h:208
L1Analysis::L1AnalysisRecoMuonDataFormat::imp_point_y
std::vector< double > imp_point_y
Definition: L1AnalysisRecoMuonDataFormat.h:166
L1Analysis::L1AnalysisRecoRpcHitDataFormat::roll
std::vector< int > roll
Definition: L1AnalysisRecoRpcHitDataFormat.h:52
DetId::Muon
Definition: DetId.h:26
L1Analysis::L1AnalysisRecoMuonDataFormat::sa_inner_pt
std::vector< double > sa_inner_pt
Definition: L1AnalysisRecoMuonDataFormat.h:219
L1MuonRecoTreeProducer::empty_hlt
void empty_hlt()
Definition: L1MuonRecoTreeProducer.cc:388
L1MuonRecoTreeProducer::theBField
edm::ESHandle< MagneticField > theBField
Definition: L1MuonRecoTreeProducer.cc:164
L1MuonRecoTreeProducer::triggerMatching_
bool triggerMatching_
Definition: L1MuonRecoTreeProducer.cc:140
L1Analysis::L1AnalysisRecoMuonDataFormat::numberOfMatchedStations
std::vector< double > numberOfMatchedStations
Definition: L1AnalysisRecoMuonDataFormat.h:162
reco::MuonTime::timeAtIpInOut
float timeAtIpInOut
Definition: MuonTime.h:13
L1Analysis::L1AnalysisRecoMuonDataFormat::time_inout_err
std::vector< double > time_inout_err
Definition: L1AnalysisRecoMuonDataFormat.h:262
L1Analysis::L1AnalysisRecoMuonDataFormat::time_outin_err
std::vector< double > time_outin_err
Definition: L1AnalysisRecoMuonDataFormat.h:264
L1Analysis::L1AnalysisRecoMuonDataFormat::tr_phi
std::vector< double > tr_phi
Definition: L1AnalysisRecoMuonDataFormat.h:182
edm_modernize_messagelogger.stat
stat
Definition: edm_modernize_messagelogger.py:27
edm::HandleBase::isValid
bool isValid() const
Definition: HandleBase.h:70
L1Analysis::L1AnalysisRecoMuonDataFormat::phi_hb
std::vector< double > phi_hb
Definition: L1AnalysisRecoMuonDataFormat.h:170
L1MuonRecoTreeProducer::triggerProcessLabel_
std::string triggerProcessLabel_
Definition: L1MuonRecoTreeProducer.cc:142
RPCDetId::layer
int layer() const
Definition: RPCDetId.h:85
L1Analysis::L1AnalysisRecoMuonDataFormat::sa_r_me1_n
std::vector< double > sa_r_me1_n
Definition: L1AnalysisRecoMuonDataFormat.h:250
L1Analysis::L1AnalysisRecoMuonDataFormat::tr_phi_me2_n
std::vector< double > tr_phi_me2_n
Definition: L1AnalysisRecoMuonDataFormat.h:198
L1Analysis::L1AnalysisRecoMuonDataFormat::pt
std::vector< double > pt
Definition: L1AnalysisRecoMuonDataFormat.h:157
L1Analysis::L1AnalysisRecoMuonDataFormat::hlt_mu
std::vector< int > hlt_mu
Definition: L1AnalysisRecoMuonDataFormat.h:276
L1Analysis::L1AnalysisRecoMuonDataFormat::type
std::vector< int > type
Definition: L1AnalysisRecoMuonDataFormat.h:153
MuonGeometryRecord
Definition: MuonGeometryRecord.h:34
L1Analysis::L1AnalysisRecoMuonDataFormat::hlt_deltaR
std::vector< double > hlt_deltaR
Definition: L1AnalysisRecoMuonDataFormat.h:278
L1MuonRecoTreeProducer::triggerNames_
std::vector< std::string > triggerNames_
Definition: L1MuonRecoTreeProducer.cc:144
L1MuonRecoTreeProducer::propagatorOpposite
edm::ESHandle< Propagator > propagatorOpposite
Definition: L1MuonRecoTreeProducer.cc:171
edm::InputTag
Definition: InputTag.h:15
PV3DBase::phi
Geom::Phi< T > phi() const
Definition: PV3DBase.h:66
L1MuonRecoTreeProducer::empty_global
void empty_global()
Definition: L1MuonRecoTreeProducer.cc:263
CSCGeometry::chamber
const CSCChamber * chamber(CSCDetId id) const
Return the chamber corresponding to given DetId.
Definition: CSCGeometry.cc:100
L1Analysis::L1AnalysisRecoRpcHitDataFormat::ring
std::vector< int > ring
Definition: L1AnalysisRecoRpcHitDataFormat.h:53
geometryCSVtoXML.xx
xx
Definition: geometryCSVtoXML.py:19
TrajectoryStateOnSurface::isValid
bool isValid() const
Definition: TrajectoryStateOnSurface.h:54
hit
Definition: SiStripHitEffFromCalibTree.cc:88
L1Analysis::L1AnalysisRecoMuonDataFormat::tr_ch
std::vector< double > tr_ch
Definition: L1AnalysisRecoMuonDataFormat.h:178
L1Analysis::L1AnalysisRecoMuonDataFormat::tr_imp_point_x
std::vector< double > tr_imp_point_x
Definition: L1AnalysisRecoMuonDataFormat.h:187
L1Analysis::L1AnalysisRecoMuonDataFormat::hlt_isomu
std::vector< int > hlt_isomu
Definition: L1AnalysisRecoMuonDataFormat.h:275
L1Analysis::L1AnalysisRecoRpcHitDataFormat::subsector
std::vector< int > subsector
Definition: L1AnalysisRecoRpcHitDataFormat.h:51
L1Analysis::L1AnalysisRecoMuonDataFormat::tr_eta
std::vector< double > tr_eta
Definition: L1AnalysisRecoMuonDataFormat.h:181
L1Analysis::L1AnalysisRecoMuonDataFormat::tr_validhits
std::vector< double > tr_validhits
Definition: L1AnalysisRecoMuonDataFormat.h:183
TrackingComponentsRecord
Definition: TrackingComponentsRecord.h:12
Generated for CMSSW Reference Manual by doxygen 1.8.16