CMS 3D CMS Logo
 All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Properties Friends Macros Pages
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)
 
virtual std::vector
< std::unique_ptr
< TrackingRegion > > 		regions (const edm::Event &event, const edm::EventSetup &es) const override
 
virtual ~CosmicRegionalSeedGenerator ()
 
- Public Member Functions inherited from TrackingRegionProducer
virtual ~TrackingRegionProducer ()
 

Private Attributes

edm::ParameterSet conf_
 
float deltaEta_
 
float deltaPhi_
 
double deltaRExclusionSize_
 
bool doJetsExclusionCheck_
 
double jetsPtMin_
 
edm::EDGetTokenT
< MeasurementTrackerEvent > 		measurementTrackerEventToken_
 
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_
 
edm::ParameterSet regionPSet
 
float rVertex_
 
std::string thePropagatorName_
 
float zVertex_
 

Detailed Description

Definition at line 49 of file CosmicRegionalSeedGenerator.h.

Constructor & Destructor Documentation

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

Definition at line 22 of file CosmicRegionalSeedGenerator.cc.

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

23  :
24  conf_(conf)
25 {
26  edm::LogInfo ("CosmicRegionalSeedGenerator") << "Begin Run:: Constructing CosmicRegionalSeedGenerator";
27 
28  regionPSet = conf_.getParameter<edm::ParameterSet>("RegionPSet");
29  ptMin_ = regionPSet.getParameter<double>("ptMin");
30  rVertex_ = regionPSet.getParameter<double>("rVertex");
31  zVertex_ = regionPSet.getParameter<double>("zVertex");
32  deltaEta_ = regionPSet.getParameter<double>("deltaEtaRegion");
33  deltaPhi_ = regionPSet.getParameter<double>("deltaPhiRegion");
34 
35  edm::ParameterSet toolsPSet = conf_.getParameter<edm::ParameterSet>("ToolsPSet");
36  thePropagatorName_ = toolsPSet.getParameter<std::string>("thePropagatorName");
37  regionBase_ = toolsPSet.getParameter<std::string>("regionBase");
38 
39  edm::ParameterSet collectionsPSet = conf_.getParameter<edm::ParameterSet>("CollectionsPSet");
40  recoMuonsCollection_ = collectionsPSet.getParameter<edm::InputTag>("recoMuonsCollection");
41  recoTrackMuonsCollection_ = collectionsPSet.getParameter<edm::InputTag>("recoTrackMuonsCollection");
42  recoL2MuonsCollection_ = collectionsPSet.getParameter<edm::InputTag>("recoL2MuonsCollection");
43 
44  edm::ParameterSet regionInJetsCheckPSet = conf_.getParameter<edm::ParameterSet>("RegionInJetsCheckPSet");
45  doJetsExclusionCheck_ = regionInJetsCheckPSet.getParameter<bool>("doJetsExclusionCheck");
46  deltaRExclusionSize_ = regionInJetsCheckPSet.getParameter<double>("deltaRExclusionSize");
47  jetsPtMin_ = regionInJetsCheckPSet.getParameter<double>("jetsPtMin");
48  recoCaloJetsCollection_ = regionInJetsCheckPSet.getParameter<edm::InputTag>("recoCaloJetsCollection");
49  recoCaloJetsToken_ = iC.consumes<reco::CaloJetCollection>(recoCaloJetsCollection_);
50  recoMuonsToken_ = iC.consumes<reco::MuonCollection>(recoMuonsCollection_);
51  recoTrackMuonsToken_ = iC.consumes<reco::TrackCollection>(recoTrackMuonsCollection_);
52  recoL2MuonsToken_ = iC.consumes<reco::RecoChargedCandidateCollection>(recoL2MuonsCollection_);
53  measurementTrackerEventToken_ = iC.consumes<MeasurementTrackerEvent>(edm::InputTag("MeasurementTrackerEvent"));
54 
55  edm::LogInfo ("CosmicRegionalSeedGenerator") << "Reco muons collection: " << recoMuonsCollection_ << "\n"
56  << "Reco tracks muons collection: " << recoTrackMuonsCollection_<< "\n"
57  << "Reco L2 muons collection: " << recoL2MuonsCollection_;
58 }
CosmicRegionalSeedGenerator::recoTrackMuonsToken_
edm::EDGetTokenT< reco::TrackCollection > recoTrackMuonsToken_
Definition: CosmicRegionalSeedGenerator.h:77
edm::ConsumesCollector::consumes
EDGetTokenT< ProductType > consumes(edm::InputTag const &tag)
Definition: ConsumesCollector.h:52
edm::ParameterSet::getParameter
T getParameter(std::string const &) const
CosmicRegionalSeedGenerator::recoCaloJetsToken_
edm::EDGetTokenT< reco::CaloJetCollection > recoCaloJetsToken_
Definition: CosmicRegionalSeedGenerator.h:75
reco::TrackCollection
std::vector< Track > TrackCollection
collection of Tracks
Definition: TrackFwd.h:14
AlCaHLTBitMon_QueryRunRegistry.string
string string
Definition: AlCaHLTBitMon_QueryRunRegistry.py:255
reco::MuonCollection
std::vector< Muon > MuonCollection
collection of Muon objects
Definition: MuonFwd.h:9
CosmicRegionalSeedGenerator::recoMuonsToken_
edm::EDGetTokenT< reco::MuonCollection > recoMuonsToken_
Definition: CosmicRegionalSeedGenerator.h:76
CosmicRegionalSeedGenerator::measurementTrackerEventToken_
edm::EDGetTokenT< MeasurementTrackerEvent > measurementTrackerEventToken_
Definition: CosmicRegionalSeedGenerator.h:79
reco::RecoChargedCandidateCollection
std::vector< RecoChargedCandidate > RecoChargedCandidateCollection
collectin of RecoChargedCandidate objects
Definition: RecoChargedCandidateFwd.h:9
MeasurementTrackerEventProducer_cfi.MeasurementTrackerEvent
tuple MeasurementTrackerEvent
Definition: MeasurementTrackerEventProducer_cfi.py:4
HLT_FULL_cff.InputTag
tuple InputTag
Definition: HLT_FULL_cff.py:69759
CosmicRegionalSeedGenerator::recoL2MuonsToken_
edm::EDGetTokenT< reco::RecoChargedCandidateCollection > recoL2MuonsToken_
Definition: CosmicRegionalSeedGenerator.h:78
reco::CaloJetCollection
std::vector< CaloJet > CaloJetCollection
collection of CaloJet objects
Definition: CaloJetCollection.h:15
virtual CosmicRegionalSeedGenerator::~CosmicRegionalSeedGenerator ( )
inlinevirtual

Definition at line 55 of file CosmicRegionalSeedGenerator.h.

55 {}

Member Function Documentation

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

Implements TrackingRegionProducer.

Definition at line 60 of file CosmicRegionalSeedGenerator.cc.

References funct::abs(), RecoTauCleanerPlugins::charge, kinem::delta_R(), deltaEta_, deltaPhi_, deltaRExclusionSize_, doJetsExclusionCheck_, PV3DBase< T, PVType, FrameType >::eta(), edm::EventSetup::get(), edm::Ref< C, T, F >::get(), edm::EDGetTokenT< T >::isUninitialized(), edm::HandleBase::isValid(), metsig::jet, jetsPtMin_, LogDebug, mag(), HLT_FULL_cff::measurementTracker, MeasurementTrackerEventProducer_cfi::MeasurementTrackerEvent, measurementTrackerEventToken_, SurfaceOrientation::outer, PV3DBase< T, PVType, FrameType >::perp(), PV3DBase< T, PVType, FrameType >::phi(), Pi, edm::Handle< T >::product(), edm::ESHandle< class >::product(), ptMin_, recoCaloJetsToken_, recoL2MuonsCollection_, recoL2MuonsToken_, recoMuonsCollection_, recoMuonsToken_, recoTrackMuonsCollection_, recoTrackMuonsToken_, regionBase_, regionPSet, query::result, rVertex_, mathSSE::sqrt(), thePropagatorName_, Vector3DBase< T, FrameTag >::unit(), findQualityFiles::v, PV3DBase< T, PVType, FrameType >::z(), and zVertex_.

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

Member Data Documentation

edm::ParameterSet CosmicRegionalSeedGenerator::conf_
private

Definition at line 60 of file CosmicRegionalSeedGenerator.h.

Referenced by CosmicRegionalSeedGenerator().

float CosmicRegionalSeedGenerator::deltaEta_
private

Definition at line 66 of file CosmicRegionalSeedGenerator.h.

Referenced by CosmicRegionalSeedGenerator(), and regions().

float CosmicRegionalSeedGenerator::deltaPhi_
private

Definition at line 67 of file CosmicRegionalSeedGenerator.h.

Referenced by CosmicRegionalSeedGenerator(), and regions().

double CosmicRegionalSeedGenerator::deltaRExclusionSize_
private

Definition at line 83 of file CosmicRegionalSeedGenerator.h.

Referenced by CosmicRegionalSeedGenerator(), and regions().

bool CosmicRegionalSeedGenerator::doJetsExclusionCheck_
private

Definition at line 82 of file CosmicRegionalSeedGenerator.h.

Referenced by CosmicRegionalSeedGenerator(), and regions().

double CosmicRegionalSeedGenerator::jetsPtMin_
private

Definition at line 84 of file CosmicRegionalSeedGenerator.h.

Referenced by CosmicRegionalSeedGenerator(), and regions().

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

Definition at line 79 of file CosmicRegionalSeedGenerator.h.

Referenced by CosmicRegionalSeedGenerator(), and regions().

float CosmicRegionalSeedGenerator::ptMin_
private

Definition at line 63 of file CosmicRegionalSeedGenerator.h.

Referenced by CosmicRegionalSeedGenerator(), and regions().

edm::InputTag CosmicRegionalSeedGenerator::recoCaloJetsCollection_
private

Definition at line 85 of file CosmicRegionalSeedGenerator.h.

Referenced by CosmicRegionalSeedGenerator().

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

Definition at line 75 of file CosmicRegionalSeedGenerator.h.

Referenced by CosmicRegionalSeedGenerator(), and regions().

edm::InputTag CosmicRegionalSeedGenerator::recoL2MuonsCollection_
private

Definition at line 74 of file CosmicRegionalSeedGenerator.h.

Referenced by CosmicRegionalSeedGenerator(), and regions().

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

Definition at line 78 of file CosmicRegionalSeedGenerator.h.

Referenced by CosmicRegionalSeedGenerator(), and regions().

edm::InputTag CosmicRegionalSeedGenerator::recoMuonsCollection_
private

Definition at line 72 of file CosmicRegionalSeedGenerator.h.

Referenced by CosmicRegionalSeedGenerator(), and regions().

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

Definition at line 76 of file CosmicRegionalSeedGenerator.h.

Referenced by CosmicRegionalSeedGenerator(), and regions().

edm::InputTag CosmicRegionalSeedGenerator::recoTrackMuonsCollection_
private

Definition at line 73 of file CosmicRegionalSeedGenerator.h.

Referenced by CosmicRegionalSeedGenerator(), and regions().

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

Definition at line 77 of file CosmicRegionalSeedGenerator.h.

Referenced by CosmicRegionalSeedGenerator(), and regions().

std::string CosmicRegionalSeedGenerator::regionBase_
private

Definition at line 70 of file CosmicRegionalSeedGenerator.h.

Referenced by CosmicRegionalSeedGenerator(), and regions().

edm::ParameterSet CosmicRegionalSeedGenerator::regionPSet
private

Definition at line 61 of file CosmicRegionalSeedGenerator.h.

Referenced by CosmicRegionalSeedGenerator(), and regions().

float CosmicRegionalSeedGenerator::rVertex_
private

Definition at line 64 of file CosmicRegionalSeedGenerator.h.

Referenced by CosmicRegionalSeedGenerator(), and regions().

std::string CosmicRegionalSeedGenerator::thePropagatorName_
private

Definition at line 69 of file CosmicRegionalSeedGenerator.h.

Referenced by CosmicRegionalSeedGenerator(), and regions().

float CosmicRegionalSeedGenerator::zVertex_
private

Definition at line 65 of file CosmicRegionalSeedGenerator.h.

Referenced by CosmicRegionalSeedGenerator(), and regions().

Generated for CMSSW Reference Manual by doxygen 1.8.5