CMS 3D CMS Logo
 All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Properties Friends Macros Pages
HITRegionalPixelSeedGenerator.h
Go to the documentation of this file.
1 #ifndef HITRegionalPixelSeedGenerator_h
2 #define HITRegionalPixelSeedGenerator_h
3 
4 //
5 // Class: HITRegionalPixelSeedGenerator
6 
7 
8 #include "FWCore/Framework/interface/EDProducer.h"
9 #include "FWCore/Framework/interface/Event.h"
10 #include "DataFormats/Common/interface/Handle.h"
11 #include "FWCore/Framework/interface/EventSetup.h"
12 #include "FWCore/Utilities/interface/InputTag.h"
13 #include "RecoTracker/TkTrackingRegions/interface/TrackingRegionProducer.h"
14 #include "RecoTracker/TkTrackingRegions/interface/GlobalTrackingRegion.h"
15 #include "RecoTracker/TkTrackingRegions/interface/RectangularEtaPhiTrackingRegion.h"
16 #include "DataFormats/VertexReco/interface/Vertex.h"
17 #include "DataFormats/VertexReco/interface/VertexFwd.h"
18 #include "DataFormats/TrackReco/interface/Track.h"
19 #include "DataFormats/JetReco/interface/CaloJetCollection.h"
20 #include "DataFormats/Math/interface/Vector3D.h"
21 #include "DataFormats/L1Trigger/interface/L1JetParticle.h"
22 #include "DataFormats/L1Trigger/interface/L1JetParticleFwd.h"
23 #include "DataFormats/HLTReco/interface/TriggerFilterObjectWithRefs.h"
24 
25 // Math
26 #include "Math/GenVector/VectorUtil.h"
27 #include "Math/GenVector/PxPyPzE4D.h"
28 
29 #include "FWCore/ParameterSet/interface/ParameterSet.h"
30 
31 #include "DataFormats/Common/interface/Ref.h"
32 #include "DataFormats/JetReco/interface/Jet.h"
33 
34 class HITRegionalPixelSeedGenerator : public TrackingRegionProducer {
35  public:
36 
37  explicit HITRegionalPixelSeedGenerator(const edm::ParameterSet& conf_)
38  {
39  edm::LogInfo ("HITRegionalPixelSeedGenerator")<<"Enter the HITRegionalPixelSeedGenerator";
40 
41  edm::ParameterSet regionPSet = conf_.getParameter<edm::ParameterSet>("RegionPSet");
42 
43  ptmin=regionPSet.getParameter<double>("ptMin");
44  originradius=regionPSet.getParameter<double>("originRadius");
45  halflength=regionPSet.getParameter<double>("originHalfLength");
46  vertexSrc=regionPSet.getParameter<std::string>("vertexSrc");
47  etaCenter_=regionPSet.getParameter<double>("etaCenter");
48  phiCenter_=regionPSet.getParameter<double>("phiCenter");
49  deltaTrackEta = regionPSet.getParameter<double>("deltaEtaTrackRegion");
50  deltaTrackPhi = regionPSet.getParameter<double>("deltaPhiTrackRegion");
51  deltaL1JetEta = regionPSet.getParameter<double>("deltaEtaL1JetRegion");
52  deltaL1JetPhi = regionPSet.getParameter<double>("deltaPhiL1JetRegion");
53  tracksrc_ = regionPSet.getParameter<edm::InputTag>("trackSrc");
54  isoTracksrc_ = regionPSet.getParameter<edm::InputTag>("isoTrackSrc");
55  l1jetsrc_ = regionPSet.getParameter<edm::InputTag>("l1tjetSrc");
56  usejets_ = regionPSet.getParameter<bool>("useL1Jets");
57  usetracks_ = regionPSet.getParameter<bool>("useTracks");
58  useIsoTracks_ = regionPSet.getParameter<bool>("useIsoTracks");
59  fixedReg_ = regionPSet.getParameter<bool>("fixedReg");
60  }
61 
62  virtual ~HITRegionalPixelSeedGenerator() {}
63 
64 
65  virtual std::vector<TrackingRegion* > regions(const edm::Event& e, const edm::EventSetup& es) const
66  {
67  std::vector<TrackingRegion* > result;
68  float originz =0.;
69 
70  double deltaZVertex = halflength;
71 
72 
73 
74  if (usetracks_)
75  {
76  edm::Handle<reco::TrackCollection> tracks;
77  e.getByLabel(tracksrc_, tracks);
78 
79  edm::Handle<reco::VertexCollection> vertices;
80  e.getByLabel(vertexSrc,vertices);
81  const reco::VertexCollection vertCollection = *(vertices.product());
82  reco::VertexCollection::const_iterator ci = vertCollection.begin();
83 
84  if(vertCollection.size() > 0)
85  {
86  originz = ci->z();
87  }
88  else
89  {
90  deltaZVertex = 15.;
91  }
92 
93  GlobalVector globalVector(0,0,1);
94  if(tracks->size() == 0) return result;
95 
96  reco::TrackCollection::const_iterator itr = tracks->begin();
97  for(;itr != tracks->end();itr++)
98  {
99 
100  GlobalVector ptrVec((itr)->px(),(itr)->py(),(itr)->pz());
101  globalVector = ptrVec;
102 
103 
104  RectangularEtaPhiTrackingRegion* etaphiRegion = new RectangularEtaPhiTrackingRegion(globalVector,
105  GlobalPoint(0,0,originz),
106  ptmin,
107  originradius,
108  deltaZVertex,
109  deltaTrackEta,
110  deltaTrackPhi);
111  result.push_back(etaphiRegion);
112 
113  }
114  }
115 
116  if (useIsoTracks_)
117  {
118  edm::Handle<trigger::TriggerFilterObjectWithRefs> isotracks;
119  e.getByLabel(isoTracksrc_, isotracks);
120 
121  std::vector< edm::Ref<reco::IsolatedPixelTrackCandidateCollection> > isoPixTrackRefs;
122 
123  isotracks->getObjects(trigger::TriggerTrack, isoPixTrackRefs);
124 
125  edm::Handle<reco::VertexCollection> vertices;
126  e.getByLabel(vertexSrc,vertices);
127  const reco::VertexCollection vertCollection = *(vertices.product());
128  reco::VertexCollection::const_iterator ci = vertCollection.begin();
129 
130  if(vertCollection.size() > 0)
131  {
132  originz = ci->z();
133  }
134  else
135  {
136  deltaZVertex = 15.;
137  }
138 
139  GlobalVector globalVector(0,0,1);
140  if(isoPixTrackRefs.size() == 0) return result;
141 
142  for(uint32_t p=0; p<isoPixTrackRefs.size(); p++)
143  {
144  GlobalVector ptrVec((isoPixTrackRefs[p]->track())->px(),(isoPixTrackRefs[p]->track())->py(),(isoPixTrackRefs[p]->track())->pz());
145  globalVector = ptrVec;
146 
147  RectangularEtaPhiTrackingRegion* etaphiRegion = new RectangularEtaPhiTrackingRegion(globalVector,
148  GlobalPoint(0,0,originz),
149  ptmin,
150  originradius,
151  deltaZVertex,
152  deltaTrackEta,
153  deltaTrackPhi);
154  result.push_back(etaphiRegion);
155  }
156  }
157 
158  if (usejets_)
159  {
160  edm::Handle<l1extra::L1JetParticleCollection> jets;
161  e.getByLabel(l1jetsrc_, jets);
162 
163  edm::Handle<reco::VertexCollection> vertices;
164  e.getByLabel(vertexSrc,vertices);
165  const reco::VertexCollection vertCollection = *(vertices.product());
166  reco::VertexCollection::const_iterator ci = vertCollection.begin();
167  if(vertCollection.size() > 0)
168  {
169  originz = ci->z();
170  }
171  else
172  {
173  deltaZVertex = 15.;
174  }
175 
176  GlobalVector globalVector(0,0,1);
177  if(jets->size() == 0) return result;
178 
179  for (l1extra::L1JetParticleCollection::const_iterator iJet = jets->begin(); iJet != jets->end(); iJet++)
180  {
181  GlobalVector jetVector(iJet->p4().x(), iJet->p4().y(), iJet->p4().z());
182  GlobalPoint vertex(0, 0, originz);
183 
184  RectangularEtaPhiTrackingRegion* etaphiRegion = new RectangularEtaPhiTrackingRegion( jetVector,
185  vertex,
186  ptmin,
187  originradius,
188  deltaZVertex,
189  deltaL1JetEta,
190  deltaL1JetPhi );
191  result.push_back(etaphiRegion);
192  }
193  }
194  if (fixedReg_)
195  {
196  GlobalVector fixedVector(cos(phiCenter_)*sin(2*atan(exp(-etaCenter_))), sin(phiCenter_)*sin(2*atan(exp(-etaCenter_))), cos(2*atan(exp(-etaCenter_))));
197  GlobalPoint vertex(0, 0, originz);
198 
199  edm::Handle<reco::VertexCollection> vertices;
200  e.getByLabel(vertexSrc,vertices);
201  const reco::VertexCollection vertCollection = *(vertices.product());
202  reco::VertexCollection::const_iterator ci = vertCollection.begin();
203  if(vertCollection.size() > 0)
204  {
205  originz = ci->z();
206  }
207  else
208  {
209  deltaZVertex = 15.;
210  }
211 
212  RectangularEtaPhiTrackingRegion* etaphiRegion = new RectangularEtaPhiTrackingRegion( fixedVector,
213  vertex,
214  ptmin,
215  originradius,
216  deltaZVertex,
217  deltaL1JetEta,
218  deltaL1JetPhi );
219  result.push_back(etaphiRegion);
220  }
221 
222  return result;
223  }
224 
225 
226  private:
227  edm::ParameterSet conf_;
228 
229  float ptmin;
230  float originradius;
231  float halflength;
232  double etaCenter_;
233  double phiCenter_;
234  float deltaTrackEta;
235  float deltaTrackPhi;
236  float deltaL1JetEta;
237  float deltaL1JetPhi;
238  edm::InputTag tracksrc_;
239  edm::InputTag isoTracksrc_;
240  std::string vertexSrc;
241  edm::InputTag l1jetsrc_;
242  bool usejets_;
243  bool usetracks_;
244  bool fixedReg_;
245  bool useIsoTracks_;
246 
247 };
248 
249 #endif
250 
251 
252 
edm::ParameterSet::getParameter
T getParameter(std::string const &) const
create_public_lumi_plots.exp
tuple exp
Definition: create_public_lumi_plots.py:1093
funct::sin
Sin< T >::type sin(const T &t)
Definition: Sin.h:22
GlobalPoint
Global3DPoint GlobalPoint
Definition: GlobalPoint.h:10
reco::VertexCollection
std::vector< Vertex > VertexCollection
collection of Vertex objects
Definition: VertexFwd.h:9
fwrapper::jets
vector< PseudoJet > jets
Definition: fastjetfortran_madfks.cc:43
query.result
tuple result
Definition: query.py:137
funct::cos
Cos< T >::type cos(const T &t)
Definition: Cos.h:22
edm::Event::getByLabel
bool getByLabel(InputTag const &tag, Handle< PROD > &result) const
Definition: Event.h:356
HITRegionalPixelSeedGenerator::regions
virtual std::vector< TrackingRegion * > regions(const edm::Event &e, const edm::EventSetup &es) const
Definition: HITRegionalPixelSeedGenerator.h:65
testEve_cfg.tracks
tuple tracks
Definition: testEve_cfg.py:39
HITRegionalPixelSeedGenerator::HITRegionalPixelSeedGenerator
HITRegionalPixelSeedGenerator(const edm::ParameterSet &conf_)
Definition: HITRegionalPixelSeedGenerator.h:37
edm::Handle::product
T const * product() const
Definition: Handle.h:74
AlCaHLTBitMon_ParallelJobs.p
tuple p
Definition: AlCaHLTBitMon_ParallelJobs.py:152
Generated for CMSSW Reference Manual by doxygen 1.8.5