CMS 3D CMS Logo
List of all members | Public Member Functions | Private Attributes
CosmicRegionalSeedGenerator Class Reference

#include <CosmicRegionalSeedGenerator.h>

Inheritance diagram for CosmicRegionalSeedGenerator:
TrackingRegionProducer

Public Member Functions

 CosmicRegionalSeedGenerator (const edm::ParameterSet &conf, edm::ConsumesCollector &&iC)
 
std::vector< std::unique_ptr< TrackingRegion > > regions (const edm::Event &event, const edm::EventSetup &es) const override
 
 ~CosmicRegionalSeedGenerator () override
 
- Public Member Functions inherited from TrackingRegionProducer
virtual ~TrackingRegionProducer ()
 

Private Attributes

float deltaEta_
 
float deltaPhi_
 
double deltaRExclusionSize_
 
bool doJetsExclusionCheck_
 
double jetsPtMin_
 
edm::ESGetToken< MagneticField, IdealMagneticFieldRecord > magfieldToken_
 
edm::EDGetTokenT< MeasurementTrackerEvent > measurementTrackerEventToken_
 
edm::ESGetToken< Propagator, TrackingComponentsRecord > propagatorToken_
 
float ptMin_
 
edm::InputTag recoCaloJetsCollection_
 
edm::EDGetTokenT< reco::CaloJetCollection > recoCaloJetsToken_
 
edm::InputTag recoL2MuonsCollection_
 
edm::EDGetTokenT< reco::RecoChargedCandidateCollection > recoL2MuonsToken_
 
edm::InputTag recoMuonsCollection_
 
edm::EDGetTokenT< reco::MuonCollection > recoMuonsToken_
 
edm::InputTag recoTrackMuonsCollection_
 
edm::EDGetTokenT< reco::TrackCollection > recoTrackMuonsToken_
 
std::string regionBase_
 
float rVertex_
 
float zVertex_
 

Detailed Description

Definition at line 48 of file CosmicRegionalSeedGenerator.h.

Constructor & Destructor Documentation

◆ CosmicRegionalSeedGenerator()

CosmicRegionalSeedGenerator::CosmicRegionalSeedGenerator ( const edm::ParameterSet &  conf,
edm::ConsumesCollector &&  iC 
)
explicit

Definition at line 18 of file CosmicRegionalSeedGenerator.cc.

References deltaEta_, deltaPhi_, deltaRExclusionSize_, doJetsExclusionCheck_, edm::ParameterSet::getParameter(), ProducerED_cfi::InputTag, jetsPtMin_, HLTSiStripMonitoring_cff::MeasurementTrackerEvent, measurementTrackerEventToken_, propagatorToken_, ptMin_, recoCaloJetsCollection_, recoCaloJetsToken_, recoL2MuonsCollection_, recoL2MuonsToken_, recoMuonsCollection_, recoMuonsToken_, recoTrackMuonsCollection_, recoTrackMuonsToken_, regionBase_, rVertex_, AlCaHLTBitMon_QueryRunRegistry::string, and zVertex_.

19  : magfieldToken_(iC.esConsumes()) {
20  edm::LogInfo("CosmicRegionalSeedGenerator") << "Begin Run:: Constructing CosmicRegionalSeedGenerator";
21 
22  edm::ParameterSet regionPSet = conf.getParameter<edm::ParameterSet>("RegionPSet");
23  ptMin_ = regionPSet.getParameter<double>("ptMin");
24  rVertex_ = regionPSet.getParameter<double>("rVertex");
25  zVertex_ = regionPSet.getParameter<double>("zVertex");
26  deltaEta_ = regionPSet.getParameter<double>("deltaEtaRegion");
27  deltaPhi_ = regionPSet.getParameter<double>("deltaPhiRegion");
28 
29  edm::ParameterSet toolsPSet = conf.getParameter<edm::ParameterSet>("ToolsPSet");
30  propagatorToken_ = iC.esConsumes(edm::ESInputTag("", toolsPSet.getParameter<std::string>("thePropagatorName")));
31  regionBase_ = toolsPSet.getParameter<std::string>("regionBase");
32 
33  edm::ParameterSet collectionsPSet = conf.getParameter<edm::ParameterSet>("CollectionsPSet");
34  recoMuonsCollection_ = collectionsPSet.getParameter<edm::InputTag>("recoMuonsCollection");
35  recoTrackMuonsCollection_ = collectionsPSet.getParameter<edm::InputTag>("recoTrackMuonsCollection");
36  recoL2MuonsCollection_ = collectionsPSet.getParameter<edm::InputTag>("recoL2MuonsCollection");
37 
38  edm::ParameterSet regionInJetsCheckPSet = conf.getParameter<edm::ParameterSet>("RegionInJetsCheckPSet");
39  doJetsExclusionCheck_ = regionInJetsCheckPSet.getParameter<bool>("doJetsExclusionCheck");
40  deltaRExclusionSize_ = regionInJetsCheckPSet.getParameter<double>("deltaRExclusionSize");
41  jetsPtMin_ = regionInJetsCheckPSet.getParameter<double>("jetsPtMin");
42  recoCaloJetsCollection_ = regionInJetsCheckPSet.getParameter<edm::InputTag>("recoCaloJetsCollection");
43  recoCaloJetsToken_ = iC.consumes<reco::CaloJetCollection>(recoCaloJetsCollection_);
44  recoMuonsToken_ = iC.consumes<reco::MuonCollection>(recoMuonsCollection_);
45  recoTrackMuonsToken_ = iC.consumes<reco::TrackCollection>(recoTrackMuonsCollection_);
46  recoL2MuonsToken_ = iC.consumes<reco::RecoChargedCandidateCollection>(recoL2MuonsCollection_);
47  measurementTrackerEventToken_ = iC.consumes<MeasurementTrackerEvent>(edm::InputTag("MeasurementTrackerEvent"));
48 
49  edm::LogInfo("CosmicRegionalSeedGenerator") << "Reco muons collection: " << recoMuonsCollection_ << "\n"
50  << "Reco tracks muons collection: " << recoTrackMuonsCollection_ << "\n"
51  << "Reco L2 muons collection: " << recoL2MuonsCollection_;
52 }
CosmicRegionalSeedGenerator::recoTrackMuonsToken_
edm::EDGetTokenT< reco::TrackCollection > recoTrackMuonsToken_
Definition: CosmicRegionalSeedGenerator.h:74
edm::ConsumesCollector::consumes
EDGetTokenT< ProductType > consumes(edm::InputTag const &tag)
Definition: ConsumesCollector.h:55
ProducerED_cfi.InputTag
InputTag
Definition: ProducerED_cfi.py:5
edm::ParameterSet::getParameter
T getParameter(std::string const &) const
Definition: ParameterSet.h:307
CosmicRegionalSeedGenerator::recoCaloJetsToken_
edm::EDGetTokenT< reco::CaloJetCollection > recoCaloJetsToken_
Definition: CosmicRegionalSeedGenerator.h:72
reco::TrackCollection
std::vector< Track > TrackCollection
collection of Tracks
Definition: TrackFwd.h:14
CosmicRegionalSeedGenerator::magfieldToken_
edm::ESGetToken< MagneticField, IdealMagneticFieldRecord > magfieldToken_
Definition: CosmicRegionalSeedGenerator.h:67
AlCaHLTBitMon_QueryRunRegistry.string
string string
Definition: AlCaHLTBitMon_QueryRunRegistry.py:256
reco::MuonCollection
std::vector< Muon > MuonCollection
collection of Muon objects
Definition: MuonFwd.h:9
CosmicRegionalSeedGenerator::propagatorToken_
edm::ESGetToken< Propagator, TrackingComponentsRecord > propagatorToken_
Definition: CosmicRegionalSeedGenerator.h:66
CosmicRegionalSeedGenerator::recoMuonsToken_
edm::EDGetTokenT< reco::MuonCollection > recoMuonsToken_
Definition: CosmicRegionalSeedGenerator.h:73
CosmicRegionalSeedGenerator::measurementTrackerEventToken_
edm::EDGetTokenT< MeasurementTrackerEvent > measurementTrackerEventToken_
Definition: CosmicRegionalSeedGenerator.h:76
edm::LogInfo
Log< level::Info, false > LogInfo
Definition: MessageLogger.h:131
reco::RecoChargedCandidateCollection
std::vector< RecoChargedCandidate > RecoChargedCandidateCollection
collectin of RecoChargedCandidate objects
Definition: RecoChargedCandidateFwd.h:9
CosmicRegionalSeedGenerator::recoL2MuonsToken_
edm::EDGetTokenT< reco::RecoChargedCandidateCollection > recoL2MuonsToken_
Definition: CosmicRegionalSeedGenerator.h:75
HLTSiStripMonitoring_cff.MeasurementTrackerEvent
MeasurementTrackerEvent
Definition: HLTSiStripMonitoring_cff.py:222
reco::CaloJetCollection
std::vector< CaloJet > CaloJetCollection
collection of CaloJet objects
Definition: CaloJetCollection.h:15

◆ ~CosmicRegionalSeedGenerator()

CosmicRegionalSeedGenerator::~CosmicRegionalSeedGenerator ( )
inlineoverride

Definition at line 52 of file CosmicRegionalSeedGenerator.h.

52 {}

Member Function Documentation

◆ regions()

std::vector< std::unique_ptr< TrackingRegion > > CosmicRegionalSeedGenerator::regions ( const edm::Event &  event,
const edm::EventSetup &  es 
) const
overridevirtual

Implements TrackingRegionProducer.

Definition at line 54 of file CosmicRegionalSeedGenerator.cc.

References funct::abs(), ALCARECOTkAlJpsiMuMu_cff::charge, kinem::delta_R(), deltaEta_, deltaPhi_, deltaRExclusionSize_, doJetsExclusionCheck_, PV3DBase< T, PVType, FrameType >::eta(), edm::EventSetup::getData(), createfilelist::int, edm::EDGetTokenT< T >::isUninitialized(), edm::HandleBase::isValid(), metsig::jet, jetsPtMin_, HLT_2024v11_cff::L2Muon, LogDebug, mag(), volumeBasedMagneticField_160812_cfi::magfield, magfieldToken_, HLTSiStripMonitoring_cff::measurementTracker, HLTSiStripMonitoring_cff::MeasurementTrackerEvent, measurementTrackerEventToken_, SurfaceOrientation::outer, PV3DBase< T, PVType, FrameType >::perp(), PV3DBase< T, PVType, FrameType >::phi(), Pi, edm::Handle< T >::product(), TrackCandidateProducer_cfi::propagator, propagatorToken_, ptMin_, recoCaloJetsToken_, recoL2MuonsCollection_, recoL2MuonsToken_, recoMuonsCollection_, recoMuonsToken_, recoTrackMuonsCollection_, recoTrackMuonsToken_, regionBase_, mps_fire::result, rVertex_, mathSSE::sqrt(), Vector3DBase< T, FrameTag >::unit(), findQualityFiles::v, PV3DBase< T, PVType, FrameType >::z(), and zVertex_.

55  {
56  std::vector<std::unique_ptr<TrackingRegion>> result;
57 
58  const MeasurementTrackerEvent* measurementTracker = nullptr;
59  if (!measurementTrackerEventToken_.isUninitialized()) {
60  edm::Handle<MeasurementTrackerEvent> hmte;
61  event.getByToken(measurementTrackerEventToken_, hmte);
62  measurementTracker = hmte.product();
63  }
64 
65  //get the propagator
66  auto const& propagator = es.getData(propagatorToken_);
67  auto const& magfield = es.getData(magfieldToken_);
68 
69  //________________________________________
70  //
71  //Seeding on Sta muon (MC && Datas)
72  //________________________________________
73 
74  if (regionBase_ == "seedOnStaMuon" || regionBase_.empty()) {
75  LogDebug("CosmicRegionalSeedGenerator") << "Seeding on stand alone muons ";
76 
77  //get collections
78  //+++++++++++++++
79 
80  //get the muon collection
81  edm::Handle<reco::MuonCollection> muonsHandle;
82  event.getByToken(recoMuonsToken_, muonsHandle);
83  if (!muonsHandle.isValid()) {
84  edm::LogError("CollectionNotFound") << "Error::No reco muons collection (" << recoMuonsCollection_
85  << ") in the event - Please verify the name of the muon collection";
86  return result;
87  }
88 
89  LogDebug("CosmicRegionalSeedGenerator") << "Muons collection size = " << muonsHandle->size();
90 
91  //get the jet collection
92  edm::Handle<CaloJetCollection> caloJetsHandle;
93  event.getByToken(recoCaloJetsToken_, caloJetsHandle);
94 
95  //definition of the region
96  //+++++++++++++++++++++++++
97 
98  for (reco::MuonCollection::const_iterator staMuon = muonsHandle->begin(); staMuon != muonsHandle->end();
99  ++staMuon) {
100  //select sta muons
101  if (!staMuon->isStandAloneMuon()) {
102  LogDebug("CosmicRegionalSeedGenerator") << "This muon is not a stand alone muon";
103  continue;
104  }
105 
106  //bit 25 as a coverage -1.4 < eta < 1.4
107  if (abs(staMuon->standAloneMuon()->eta()) > 1.5)
108  continue;
109 
110  //debug
111  LogDebug("CosmicRegionalSeedGenerator") << "Muon stand alone found in the collection - in muons chambers: \n "
112  << "Position = " << staMuon->standAloneMuon()->outerPosition() << "\n "
113  << "Momentum = " << staMuon->standAloneMuon()->outerMomentum() << "\n "
114  << "Eta = " << staMuon->standAloneMuon()->eta() << "\n "
115  << "Phi = " << staMuon->standAloneMuon()->phi();
116 
117  //initial position, momentum, charge
118 
119  GlobalPoint initialRegionPosition(staMuon->standAloneMuon()->referencePoint().x(),
120  staMuon->standAloneMuon()->referencePoint().y(),
121  staMuon->standAloneMuon()->referencePoint().z());
122  GlobalVector initialRegionMomentum(staMuon->standAloneMuon()->momentum().x(),
123  staMuon->standAloneMuon()->momentum().y(),
124  staMuon->standAloneMuon()->momentum().z());
125  int charge = (int)staMuon->standAloneMuon()->charge();
126 
127  LogDebug("CosmicRegionalSeedGenerator") << "Initial region - Reference point of the sta muon: \n "
128  << "Position = " << initialRegionPosition << "\n "
129  << "Momentum = " << initialRegionMomentum << "\n "
130  << "Eta = " << initialRegionPosition.eta() << "\n "
131  << "Phi = " << initialRegionPosition.phi() << "\n "
132  << "Charge = " << charge;
133 
134  //propagation on the last layers of TOB
135  if (staMuon->standAloneMuon()->outerPosition().y() > 0)
136  initialRegionMomentum *= -1;
137  GlobalTrajectoryParameters glb_parameters(
138  initialRegionPosition, initialRegionMomentum, charge, propagator.magneticField());
139  FreeTrajectoryState fts(glb_parameters);
140  StateOnTrackerBound onBounds(&propagator);
141  TrajectoryStateOnSurface outer = onBounds(fts);
142 
143  if (!outer.isValid()) {
144  //edm::LogError("FailedPropagation") << "Trajectory state on surface not valid" ;
145  LogDebug("CosmicRegionalSeedGenerator") << "Trajectory state on surface not valid";
146  continue;
147  }
148 
149  //final position & momentum
150  GlobalPoint regionPosition = outer.globalPosition();
151  GlobalVector regionMom = outer.globalMomentum();
152 
153  LogDebug("CosmicRegionalSeedGenerator")
154  << "Region after propagation: \n "
155  << "Position = " << outer.globalPosition() << "\n "
156  << "Momentum = " << outer.globalMomentum() << "\n "
157  << "R = " << regionPosition.perp() << " ---- z = " << regionPosition.z() << "\n "
158  << "Eta = " << outer.globalPosition().eta() << "\n "
159  << "Phi = " << outer.globalPosition().phi();
160 
161  //step back
162  double stepBack = 1;
163  GlobalPoint center = regionPosition + stepBack * regionMom.unit();
164  GlobalVector v = stepBack * regionMom.unit();
165  LogDebug("CosmicRegionalSeedGenerator") << "Step back vector = " << v << "\n";
166 
167  //exclude region built in jets
168  if (doJetsExclusionCheck_) {
169  double delta_R_min = 1000.;
170  for (CaloJetCollection::const_iterator jet = caloJetsHandle->begin(); jet != caloJetsHandle->end(); jet++) {
171  if (jet->pt() < jetsPtMin_)
172  continue;
173 
174  double deta = center.eta() - jet->eta();
175  double dphi = fabs(center.phi() - jet->phi());
176  if (dphi > TMath::Pi())
177  dphi = 2 * TMath::Pi() - dphi;
178 
179  double delta_R = sqrt(deta * deta + dphi * dphi);
180  if (delta_R < delta_R_min)
181  delta_R_min = delta_R;
182 
183  } //end loop on jets
184 
185  if (delta_R_min < deltaRExclusionSize_) {
186  LogDebug("CosmicRegionalSeedGenerator") << "Region built too close from a jet";
187  continue;
188  }
189  } //end if doJetsExclusionCheck
190 
191  //definition of the region
192 
193  result.push_back(std::make_unique<CosmicTrackingRegion>(
194  (-1) * regionMom, center, ptMin_, rVertex_, zVertex_, deltaEta_, deltaPhi_, magfield, measurementTracker));
195 
196  LogDebug("CosmicRegionalSeedGenerator")
197  << "Final CosmicTrackingRegion \n "
198  << "Position = " << center << "\n "
199  << "Direction = " << result.back()->direction() << "\n "
200  << "Distance from the region on the layer = " << (regionPosition - center).mag() << "\n "
201  << "Eta = " << center.eta() << "\n "
202  << "Phi = " << center.phi();
203 
204  } //end loop on muons
205 
206  } //end if SeedOnStaMuon
207 
208  //________________________________________
209  //
210  //Seeding on cosmic muons (MC && Datas)
211  //________________________________________
212 
213  if (regionBase_ == "seedOnCosmicMuon") {
214  LogDebug("CosmicRegionalSeedGenerator") << "Seeding on cosmic muons tracks";
215 
216  //get collections
217  //+++++++++++++++
218 
219  //get the muon collection
220  edm::Handle<reco::TrackCollection> cosmicMuonsHandle;
221  event.getByToken(recoTrackMuonsToken_, cosmicMuonsHandle);
222  if (!cosmicMuonsHandle.isValid()) {
223  edm::LogError("CollectionNotFound")
224  << "Error::No cosmic muons collection (" << recoTrackMuonsCollection_
225  << ") in the event - Please verify the name of the muon reco track collection";
226  return result;
227  }
228 
229  LogDebug("CosmicRegionalSeedGenerator") << "Cosmic muons tracks collection size = " << cosmicMuonsHandle->size();
230 
231  //get the jet collection
232  edm::Handle<CaloJetCollection> caloJetsHandle;
233  event.getByToken(recoCaloJetsToken_, caloJetsHandle);
234 
235  //definition of the region
236  //+++++++++++++++++++++++++
237 
238  for (reco::TrackCollection::const_iterator cosmicMuon = cosmicMuonsHandle->begin();
239  cosmicMuon != cosmicMuonsHandle->end();
240  ++cosmicMuon) {
241  //bit 25 as a coverage -1.4 < eta < 1.4
242  if (abs(cosmicMuon->eta()) > 1.5)
243  continue;
244 
245  //initial position, momentum, charge
246  GlobalPoint initialRegionPosition(
247  cosmicMuon->referencePoint().x(), cosmicMuon->referencePoint().y(), cosmicMuon->referencePoint().z());
248  GlobalVector initialRegionMomentum(
249  cosmicMuon->momentum().x(), cosmicMuon->momentum().y(), cosmicMuon->momentum().z());
250  int charge = (int)cosmicMuon->charge();
251 
252  LogDebug("CosmicRegionalSeedGenerator") << "Position and momentum of the muon track in the muon chambers: \n "
253  << "x = " << cosmicMuon->outerPosition().x() << "\n "
254  << "y = " << cosmicMuon->outerPosition().y() << "\n "
255  << "y = " << cosmicMuon->pt() << "\n "
256  << "Initial region - Reference point of the cosmic muon track: \n "
257  << "Position = " << initialRegionPosition << "\n "
258  << "Momentum = " << initialRegionMomentum << "\n "
259  << "Eta = " << initialRegionPosition.eta() << "\n "
260  << "Phi = " << initialRegionPosition.phi() << "\n "
261  << "Charge = " << charge;
262 
263  //propagation on the last layers of TOB
264  if (cosmicMuon->outerPosition().y() > 0 && cosmicMuon->momentum().y() < 0)
265  initialRegionMomentum *= -1;
266  GlobalTrajectoryParameters glb_parameters(
267  initialRegionPosition, initialRegionMomentum, charge, propagator.magneticField());
268  FreeTrajectoryState fts(glb_parameters);
269  StateOnTrackerBound onBounds(&propagator);
270  TrajectoryStateOnSurface outer = onBounds(fts);
271 
272  if (!outer.isValid()) {
273  //edm::LogError("FailedPropagation") << "Trajectory state on surface not valid" ;
274  LogDebug("CosmicRegionalSeedGenerator") << "Trajectory state on surface not valid";
275  continue;
276  }
277 
278  //final position & momentum
279  GlobalPoint regionPosition = outer.globalPosition();
280  GlobalVector regionMom = outer.globalMomentum();
281 
282  LogDebug("CosmicRegionalSeedGenerator")
283  << "Region after propagation: \n "
284  << "Position = " << outer.globalPosition() << "\n "
285  << "Momentum = " << outer.globalMomentum() << "\n "
286  << "R = " << regionPosition.perp() << " ---- z = " << regionPosition.z() << "\n "
287  << "Eta = " << outer.globalPosition().eta() << "\n "
288  << "Phi = " << outer.globalPosition().phi();
289 
290  //step back
291  double stepBack = 1;
292  GlobalPoint center = regionPosition + stepBack * regionMom.unit();
293  GlobalVector v = stepBack * regionMom.unit();
294  LogDebug("CosmicRegionalSeedGenerator") << "Step back vector = " << v << "\n";
295 
296  //exclude region built in jets
297  if (doJetsExclusionCheck_) {
298  double delta_R_min = 1000.;
299  for (CaloJetCollection::const_iterator jet = caloJetsHandle->begin(); jet != caloJetsHandle->end(); jet++) {
300  if (jet->pt() < jetsPtMin_)
301  continue;
302 
303  double deta = center.eta() - jet->eta();
304  double dphi = fabs(center.phi() - jet->phi());
305  if (dphi > TMath::Pi())
306  dphi = 2 * TMath::Pi() - dphi;
307 
308  double delta_R = sqrt(deta * deta + dphi * dphi);
309  if (delta_R < delta_R_min)
310  delta_R_min = delta_R;
311 
312  } //end loop on jets
313 
314  if (delta_R_min < deltaRExclusionSize_) {
315  LogDebug("CosmicRegionalSeedGenerator") << "Region built too close from a jet";
316  continue;
317  }
318  } // end if doJetsExclusionCheck
319 
320  //definition of the region
321  result.push_back(std::make_unique<CosmicTrackingRegion>(
322  (-1) * regionMom, center, ptMin_, rVertex_, zVertex_, deltaEta_, deltaPhi_, magfield, measurementTracker));
323 
324  LogDebug("CosmicRegionalSeedGenerator")
325  << "Final CosmicTrackingRegion \n "
326  << "Position = " << center << "\n "
327  << "Direction = " << result.back()->direction() << "\n "
328  << "Distance from the region on the layer = " << (regionPosition - center).mag() << "\n "
329  << "Eta = " << center.eta() << "\n "
330  << "Phi = " << center.phi();
331 
332  } //end loop on muons
333 
334  } //end if SeedOnCosmicMuon
335 
336  //________________________________________
337  //
338  //Seeding on L2 muons (MC && Datas)
339  //________________________________________
340 
341  if (regionBase_ == "seedOnL2Muon") {
342  LogDebug("CosmicRegionalSeedGenerator") << "Seeding on L2 muons";
343 
344  //get collections
345  //+++++++++++++++
346 
347  //get the muon collection
348  edm::Handle<reco::RecoChargedCandidateCollection> L2MuonsHandle;
349  event.getByToken(recoL2MuonsToken_, L2MuonsHandle);
350 
351  if (!L2MuonsHandle.isValid()) {
352  edm::LogError("CollectionNotFound") << "Error::No L2 muons collection (" << recoL2MuonsCollection_
353  << ") in the event - Please verify the name of the L2 muon collection";
354  return result;
355  }
356 
357  LogDebug("CosmicRegionalSeedGenerator") << "L2 muons collection size = " << L2MuonsHandle->size();
358 
359  //definition of the region
360  //+++++++++++++++++++++++++
361 
362  for (reco::RecoChargedCandidateCollection::const_iterator L2Muon = L2MuonsHandle->begin();
363  L2Muon != L2MuonsHandle->end();
364  ++L2Muon) {
365  reco::TrackRef tkL2Muon = L2Muon->get<reco::TrackRef>();
366 
367  //bit 25 as a coverage -1.4 < eta < 1.4
368  if (abs(tkL2Muon->eta()) > 1.5)
369  continue;
370 
371  //initial position, momentum, charge
372  GlobalPoint initialRegionPosition(
373  tkL2Muon->referencePoint().x(), tkL2Muon->referencePoint().y(), tkL2Muon->referencePoint().z());
374  GlobalVector initialRegionMomentum(tkL2Muon->momentum().x(), tkL2Muon->momentum().y(), tkL2Muon->momentum().z());
375  int charge = (int)tkL2Muon->charge();
376 
377  LogDebug("CosmicRegionalSeedGenerator") << "Position and momentum of the L2 muon track in the muon chambers: \n "
378  << "x = " << tkL2Muon->outerPosition().x() << "\n "
379  << "y = " << tkL2Muon->outerPosition().y() << "\n "
380  << "y = " << tkL2Muon->pt() << "\n "
381  << "Initial region - Reference point of the L2 muon track: \n "
382  << "Position = " << initialRegionPosition << "\n "
383  << "Momentum = " << initialRegionMomentum << "\n "
384  << "Eta = " << initialRegionPosition.eta() << "\n "
385  << "Phi = " << initialRegionPosition.phi() << "\n "
386  << "Charge = " << charge;
387 
388  //seeding only in the bottom
389  if (tkL2Muon->outerPosition().y() > 0) {
390  LogDebug("CosmicRegionalSeedGenerator") << "L2 muon in the TOP --- Region not created";
391  return result;
392  }
393 
394  GlobalTrajectoryParameters glb_parameters(
395  initialRegionPosition, initialRegionMomentum, charge, propagator.magneticField());
396  FreeTrajectoryState fts(glb_parameters);
397  StateOnTrackerBound onBounds(&propagator);
398  TrajectoryStateOnSurface outer = onBounds(fts);
399 
400  if (!outer.isValid()) {
401  //edm::LogError("FailedPropagation") << "Trajectory state on surface not valid" ;
402  LogDebug("CosmicRegionalSeedGenerator") << "Trajectory state on surface not valid";
403  continue;
404  }
405 
406  //final position & momentum
407  GlobalPoint regionPosition = outer.globalPosition();
408  GlobalVector regionMom = outer.globalMomentum();
409 
410  LogDebug("CosmicRegionalSeedGenerator")
411  << "Region after propagation: \n "
412  << "Position = " << outer.globalPosition() << "\n "
413  << "Momentum = " << outer.globalMomentum() << "\n "
414  << "R = " << regionPosition.perp() << " ---- z = " << regionPosition.z() << "\n "
415  << "Eta = " << outer.globalPosition().eta() << "\n "
416  << "Phi = " << outer.globalPosition().phi();
417 
418  //step back
419  double stepBack = 1;
420  GlobalPoint center = regionPosition + stepBack * regionMom.unit();
421  GlobalVector v = stepBack * regionMom.unit();
422  LogDebug("CosmicRegionalSeedGenerator") << "Step back vector = " << v << "\n";
423 
424  //definition of the region
425  result.push_back(std::make_unique<CosmicTrackingRegion>(
426  (-1) * regionMom, center, ptMin_, rVertex_, zVertex_, deltaEta_, deltaPhi_, magfield, measurementTracker));
427 
428  LogDebug("CosmicRegionalSeedGenerator")
429  << "Final L2TrackingRegion \n "
430  << "Position = " << center << "\n "
431  << "Direction = " << result.back()->direction() << "\n "
432  << "Distance from the region on the layer = " << (regionPosition - center).mag() << "\n "
433  << "Eta = " << center.eta() << "\n "
434  << "Phi = " << center.phi();
435 
436  } //end loop on muons
437 
438  } //end if SeedOnL2Muon
439 
440  return result;
441 }
CosmicRegionalSeedGenerator::recoTrackMuonsToken_
edm::EDGetTokenT< reco::TrackCollection > recoTrackMuonsToken_
Definition: CosmicRegionalSeedGenerator.h:74
Pi
const double Pi
Definition: CosmicMuonParameters.h:18
PV3DBase::perp
T perp() const
Definition: PV3DBase.h:69
CosmicRegionalSeedGenerator::recoCaloJetsToken_
edm::EDGetTokenT< reco::CaloJetCollection > recoCaloJetsToken_
Definition: CosmicRegionalSeedGenerator.h:72
edm::EventSetup::getData
T const & getData(const ESGetToken< T, R > &iToken) const noexcept(false)
Definition: EventSetup.h:119
mps_fire.result
result
Definition: mps_fire.py:311
PV3DBase::z
T z() const
Definition: PV3DBase.h:61
HLTSiStripMonitoring_cff.measurementTracker
measurementTracker
Definition: HLTSiStripMonitoring_cff.py:178
TrackCandidateProducer_cfi.propagator
propagator
Definition: TrackCandidateProducer_cfi.py:17
PV3DBase::phi
Geom::Phi< T > phi() const
Definition: PV3DBase.h:66
PV3DBase::eta
T eta() const
Definition: PV3DBase.h:73
edm::Handle::product
T const * product() const
Definition: Handle.h:70
edm::EDGetTokenT::isUninitialized
constexpr bool isUninitialized() const noexcept
Definition: EDGetToken.h:104
CosmicRegionalSeedGenerator::magfieldToken_
edm::ESGetToken< MagneticField, IdealMagneticFieldRecord > magfieldToken_
Definition: CosmicRegionalSeedGenerator.h:67
edm::LogError
Log< level::Error, false > LogError
Definition: MessageLogger.h:129
HLT_2024v11_cff.L2Muon
L2Muon
Definition: HLT_2024v11_cff.py:8391
CosmicRegionalSeedGenerator::propagatorToken_
edm::ESGetToken< Propagator, TrackingComponentsRecord > propagatorToken_
Definition: CosmicRegionalSeedGenerator.h:66
mathSSE::sqrt
T sqrt(T t)
Definition: SSEVec.h:19
funct::abs
Abs< T >::type abs(const T &t)
Definition: Abs.h:22
CosmicRegionalSeedGenerator::recoMuonsToken_
edm::EDGetTokenT< reco::MuonCollection > recoMuonsToken_
Definition: CosmicRegionalSeedGenerator.h:73
kinem::delta_R
double delta_R(double eta1, double phi1, double eta2, double phi2)
Definition: AnglesUtil.h:106
ALCARECOTkAlJpsiMuMu_cff.charge
charge
Definition: ALCARECOTkAlJpsiMuMu_cff.py:47
CosmicRegionalSeedGenerator::measurementTrackerEventToken_
edm::EDGetTokenT< MeasurementTrackerEvent > measurementTrackerEventToken_
Definition: CosmicRegionalSeedGenerator.h:76
mag
T mag() const
The vector magnitude. Equivalent to sqrt(vec.mag2())
Definition: Basic3DVectorLD.h:130
edm::HandleBase::isValid
bool isValid() const
Definition: HandleBase.h:70
Vector3DBase::unit
Vector3DBase unit() const
Definition: Vector3DBase.h:54
CosmicRegionalSeedGenerator::recoL2MuonsToken_
edm::EDGetTokenT< reco::RecoChargedCandidateCollection > recoL2MuonsToken_
Definition: CosmicRegionalSeedGenerator.h:75
HLTSiStripMonitoring_cff.MeasurementTrackerEvent
MeasurementTrackerEvent
Definition: HLTSiStripMonitoring_cff.py:222
volumeBasedMagneticField_160812_cfi.magfield
magfield
Definition: volumeBasedMagneticField_160812_cfi.py:11
LogDebug
#define LogDebug(id)
Definition: MessageLogger.h:241

Member Data Documentation

◆ deltaEta_

float CosmicRegionalSeedGenerator::deltaEta_
private

Definition at line 61 of file CosmicRegionalSeedGenerator.h.

Referenced by CosmicRegionalSeedGenerator(), and regions().

◆ deltaPhi_

float CosmicRegionalSeedGenerator::deltaPhi_
private

Definition at line 62 of file CosmicRegionalSeedGenerator.h.

Referenced by CosmicRegionalSeedGenerator(), and regions().

◆ deltaRExclusionSize_

double CosmicRegionalSeedGenerator::deltaRExclusionSize_
private

Definition at line 79 of file CosmicRegionalSeedGenerator.h.

Referenced by CosmicRegionalSeedGenerator(), and regions().

◆ doJetsExclusionCheck_

bool CosmicRegionalSeedGenerator::doJetsExclusionCheck_
private

Definition at line 78 of file CosmicRegionalSeedGenerator.h.

Referenced by CosmicRegionalSeedGenerator(), and regions().

◆ jetsPtMin_

double CosmicRegionalSeedGenerator::jetsPtMin_
private

Definition at line 80 of file CosmicRegionalSeedGenerator.h.

Referenced by CosmicRegionalSeedGenerator(), and regions().

◆ magfieldToken_

edm::ESGetToken<MagneticField, IdealMagneticFieldRecord> CosmicRegionalSeedGenerator::magfieldToken_
private

Definition at line 67 of file CosmicRegionalSeedGenerator.h.

Referenced by regions().

◆ measurementTrackerEventToken_

edm::EDGetTokenT<MeasurementTrackerEvent> CosmicRegionalSeedGenerator::measurementTrackerEventToken_
private

Definition at line 76 of file CosmicRegionalSeedGenerator.h.

Referenced by CosmicRegionalSeedGenerator(), and regions().

◆ propagatorToken_

edm::ESGetToken<Propagator, TrackingComponentsRecord> CosmicRegionalSeedGenerator::propagatorToken_
private

Definition at line 66 of file CosmicRegionalSeedGenerator.h.

Referenced by CosmicRegionalSeedGenerator(), and regions().

◆ ptMin_

float CosmicRegionalSeedGenerator::ptMin_
private

Definition at line 58 of file CosmicRegionalSeedGenerator.h.

Referenced by CosmicRegionalSeedGenerator(), and regions().

◆ recoCaloJetsCollection_

edm::InputTag CosmicRegionalSeedGenerator::recoCaloJetsCollection_
private

Definition at line 81 of file CosmicRegionalSeedGenerator.h.

Referenced by CosmicRegionalSeedGenerator().

◆ recoCaloJetsToken_

edm::EDGetTokenT<reco::CaloJetCollection> CosmicRegionalSeedGenerator::recoCaloJetsToken_
private

Definition at line 72 of file CosmicRegionalSeedGenerator.h.

Referenced by CosmicRegionalSeedGenerator(), and regions().

◆ recoL2MuonsCollection_

edm::InputTag CosmicRegionalSeedGenerator::recoL2MuonsCollection_
private

Definition at line 71 of file CosmicRegionalSeedGenerator.h.

Referenced by CosmicRegionalSeedGenerator(), and regions().

◆ recoL2MuonsToken_

edm::EDGetTokenT<reco::RecoChargedCandidateCollection> CosmicRegionalSeedGenerator::recoL2MuonsToken_
private

Definition at line 75 of file CosmicRegionalSeedGenerator.h.

Referenced by CosmicRegionalSeedGenerator(), and regions().

◆ recoMuonsCollection_

edm::InputTag CosmicRegionalSeedGenerator::recoMuonsCollection_
private

Definition at line 69 of file CosmicRegionalSeedGenerator.h.

Referenced by CosmicRegionalSeedGenerator(), and regions().

◆ recoMuonsToken_

edm::EDGetTokenT<reco::MuonCollection> CosmicRegionalSeedGenerator::recoMuonsToken_
private

Definition at line 73 of file CosmicRegionalSeedGenerator.h.

Referenced by CosmicRegionalSeedGenerator(), and regions().

◆ recoTrackMuonsCollection_

edm::InputTag CosmicRegionalSeedGenerator::recoTrackMuonsCollection_
private

Definition at line 70 of file CosmicRegionalSeedGenerator.h.

Referenced by CosmicRegionalSeedGenerator(), and regions().

◆ recoTrackMuonsToken_

edm::EDGetTokenT<reco::TrackCollection> CosmicRegionalSeedGenerator::recoTrackMuonsToken_
private

Definition at line 74 of file CosmicRegionalSeedGenerator.h.

Referenced by CosmicRegionalSeedGenerator(), and regions().

◆ regionBase_

std::string CosmicRegionalSeedGenerator::regionBase_
private

Definition at line 64 of file CosmicRegionalSeedGenerator.h.

Referenced by CosmicRegionalSeedGenerator(), and regions().

◆ rVertex_

float CosmicRegionalSeedGenerator::rVertex_
private

Definition at line 59 of file CosmicRegionalSeedGenerator.h.

Referenced by CosmicRegionalSeedGenerator(), and regions().

◆ zVertex_

float CosmicRegionalSeedGenerator::zVertex_
private

Definition at line 60 of file CosmicRegionalSeedGenerator.h.

Referenced by CosmicRegionalSeedGenerator(), and regions().

Generated for CMSSW Reference Manual by doxygen 1.8.14