CMS 3D CMS Logo
 All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Properties Friends Macros Pages
List of all members | Classes | Public Member Functions | Private Member Functions | Private Attributes
EvtPlaneProducer Class Reference

#include <RecoHI/EvtPlaneProducer/src/EvtPlaneProducer.cc>

Inheritance diagram for EvtPlaneProducer:
edm::EDProducer edm::ProducerBase edm::ProductRegistryHelper

Classes

class  GenPlane
 

Public Member Functions

 EvtPlaneProducer (const edm::ParameterSet &)
 
 ~EvtPlaneProducer ()
 
- Public Member Functions inherited from edm::EDProducer
 EDProducer ()
 
virtual ~EDProducer ()
 
- Public Member Functions inherited from edm::ProducerBase
 ProducerBase ()
 
void registerProducts (ProducerBase *, ProductRegistry *, ModuleDescription const &)
 
boost::function< void(const
BranchDescription &)> 		registrationCallback () const
 used by the fwk to register list of products More...
 
virtual ~ProducerBase ()
 

Private Member Functions

virtual void beginJob ()
 
virtual void endJob ()
 
virtual void produce (edm::Event &, const edm::EventSetup &)
 

Private Attributes

edm::InputTag caloCollection_
 
double chi2_
 
double dzerr_
 
double maxet_
 
double maxpt_
 
double maxvtx_
 
double minet_
 
double minpt_
 
double minvtx_
 
GenPlane * rp [NumEPNames]
 
bool storeNames_
 
edm::InputTag trackCollection_
 
bool useECAL_
 
bool useHCAL_
 
bool useTrack_
 
bool useTrackPtWeight_
 
edm::InputTag vtxCollection_
 

Additional Inherited Members

- Public Types inherited from edm::EDProducer
typedef EDProducer ModuleType
 
typedef WorkerT< EDProducer > WorkerType
 
- Public Types inherited from edm::ProducerBase
typedef
ProductRegistryHelper::TypeLabelList 		TypeLabelList
 
- Static Public Member Functions inherited from edm::EDProducer
static const std::string & baseType ()
 
static void fillDescriptions (ConfigurationDescriptions &descriptions)
 
static void prevalidate (ConfigurationDescriptions &descriptions)
 
- Protected Member Functions inherited from edm::EDProducer
CurrentProcessingContext const * currentContext () const
 
- Protected Member Functions inherited from edm::ProducerBase
template<class TProducer , class TMethod >
void callWhenNewProductsRegistered (TProducer *iProd, TMethod iMethod)
 

Detailed Description

Description: <one line="" class="" summary>="">

Implementation: <Notes on="" implementation>="">

Definition at line 68 of file EvtPlaneProducer.cc.

Constructor & Destructor Documentation

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

Definition at line 168 of file EvtPlaneProducer.cc.

References data, hi::EPEtaMax1, hi::EPEtaMax2, hi::EPEtaMin1, hi::EPEtaMin2, hi::EPNames, hi::EPOrder, edm::ParameterSet::getParameter(), edm::ParameterSet::getUntrackedParameter(), i, and hi::NumEPNames.

169 {
170 
171 
172  //register your products
173  vtxCollection_ = iConfig.getParameter<edm::InputTag>("vtxCollection_");
174  caloCollection_ = iConfig.getParameter<edm::InputTag>("caloCollection_");
175  trackCollection_ = iConfig.getParameter<edm::InputTag>("trackCollection_");
176  useECAL_ = iConfig.getUntrackedParameter<bool>("useECAL_",true);
177  useHCAL_ = iConfig.getUntrackedParameter<bool>("useHCAL_",true);
178  useTrack_ = iConfig.getUntrackedParameter<bool>("useTrack",true);
179  useTrackPtWeight_ = iConfig.getUntrackedParameter<bool>("useTrackPtWeight_",true);
180  minet_ = iConfig.getUntrackedParameter<double>("minet_",0.2);
181  maxet_ = iConfig.getUntrackedParameter<double>("maxet_",500.);
182  minpt_ = iConfig.getUntrackedParameter<double>("minpt_",0.3);
183  maxpt_ = iConfig.getUntrackedParameter<double>("maxpt_",2.5);
184  minvtx_ = iConfig.getUntrackedParameter<double>("minvtx_",-50.);
185  maxvtx_ = iConfig.getUntrackedParameter<double>("maxvtx_",50.);
186  dzerr_ = iConfig.getUntrackedParameter<double>("dzerr_",10.);
187  chi2_ = iConfig.getUntrackedParameter<double>("chi2_",40.);
188 
189  storeNames_ = 1;
190  produces<reco::EvtPlaneCollection>("recoLevel");
191  for(int i = 0; i<NumEPNames; i++ ) {
192  rp[i] = new GenPlane(EPNames[i].data(),EPEtaMin1[i],EPEtaMax1[i],EPEtaMin2[i],EPEtaMax2[i],EPOrder[i]);
193  }
194 }
edm::ParameterSet::getParameter
T getParameter(std::string const &) const
edm::ParameterSet::getUntrackedParameter
T getUntrackedParameter(std::string const &, T const &) const
hi::EPOrder
const int EPOrder[]
Definition: HiEvtPlaneList.h:120
i
int i
Definition: DBlmapReader.cc:9
hi::EPEtaMin2
const double EPEtaMin2[]
Definition: HiEvtPlaneList.h:167
hi::EPNames
const std::string EPNames[]
Definition: HiEvtPlaneList.h:103
hi::EPEtaMin1
const double EPEtaMin1[]
Definition: HiEvtPlaneList.h:136
EvtPlaneProducer::rp
GenPlane * rp[NumEPNames]
Definition: EvtPlaneProducer.cc:130
data
char data[epos_bytes_allocation]
Definition: EPOS_Wrapper.h:82
EvtPlaneProducer::caloCollection_
edm::InputTag caloCollection_
Definition: EvtPlaneProducer.cc:138
EvtPlaneProducer::vtxCollection_
edm::InputTag vtxCollection_
Definition: EvtPlaneProducer.cc:137
hi::EPEtaMax2
const double EPEtaMax2[]
Definition: HiEvtPlaneList.h:182
hi::NumEPNames
static const int NumEPNames
Definition: HiEvtPlaneList.h:198
EvtPlaneProducer::trackCollection_
edm::InputTag trackCollection_
Definition: EvtPlaneProducer.cc:139
hi::EPEtaMax1
const double EPEtaMax1[]
Definition: HiEvtPlaneList.h:152
EvtPlaneProducer::~EvtPlaneProducer ( )

Definition at line 197 of file EvtPlaneProducer.cc.

198 {
199 
200  // do anything here that needs to be done at desctruction time
201  // (e.g. close files, deallocate resources etc.)
202 
203 }

Member Function Documentation

void EvtPlaneProducer::beginJob ( void  )
privatevirtual

Reimplemented from edm::EDProducer.

Definition at line 444 of file EvtPlaneProducer.cc.

445 {
446 }
void EvtPlaneProducer::endJob ( void  )
privatevirtual

Reimplemented from edm::EDProducer.

Definition at line 450 of file EvtPlaneProducer.cc.

450  {
451 }
void EvtPlaneProducer::produce ( edm::Event &  iEvent,
const edm::EventSetup &  iSetup 
)
privatevirtual

Implements edm::EDProducer.

Definition at line 212 of file EvtPlaneProducer.cc.

References trackerHits::c, HiggsCouplings::c5, funct::cos(), hi::EPNames, hi::EPTracksMid3, hi::EPTracksMid4, hi::EPTracksMid5, hi::EPTracksMid6, hi::EPTracksNeg3, hi::EPTracksNeg4, hi::EPTracksNeg5, hi::EPTracksNeg6, hi::EPTracksPos3, hi::EPTracksPos4, hi::EPTracksPos5, hi::EPTracksPos6, hi::etCaloHFM, hi::etCaloHFP, hi::etEcal, hi::etEcalM, hi::etEcalP, hi::etHcal, hi::etHcalM, hi::etHcalP, hi::etHF, hi::etHF3, hi::etHF4, hi::etHF5, hi::etHF6, hi::etHFm, hi::etHFm3, hi::etHFm4, hi::etHFm5, hi::etHFm6, hi::etHFp, hi::etHFp3, hi::etHFp4, hi::etHFp5, hi::etHFp6, hi::EvtPlaneFromTracksMidEta, hi::EvtPTracksNegEtaGap, hi::EvtPTracksPosEtaGap, edm::Event::getByLabel(), i, edm::HandleBase::isValid(), j, hi::NumEPNames, edm::Handle< T >::product(), edm::Event::put(), dt_dqm_sourceclient_common_cff::reco, reset(), alignCSCRings::s, indexGen::s2, funct::sin(), mathSSE::sqrt(), testEve_cfg::tracks, and w().

213 {
214  using namespace edm;
215  using namespace reco;
216 
217  int vs_sell;
218  float vzr_sell;
219  //
220  //Get Vertex
221  //
222  edm::Handle<reco::VertexCollection> vertexCollection3;
223  iEvent.getByLabel(vtxCollection_,vertexCollection3);
224  const reco::VertexCollection * vertices3 = vertexCollection3.product();
225  vs_sell = vertices3->size();
226  if(vs_sell>0) {
227  vzr_sell = vertices3->begin()->z();
228  } else
229  vzr_sell = -999.9;
230  //
231  for(int i = 0; i<NumEPNames; i++) rp[i]->reset();
232  if(vzr_sell>minvtx_ && vzr_sell<maxvtx_) {
233  //calorimetry part
234 
235  double tower_eta, tower_phi;
236  double tower_energyet, tower_energyet_e, tower_energyet_h;
237  double s1, s2, s13,s23,s14,s24,s15,s25,s16,s26;
238  Handle<CaloTowerCollection> calotower;
239  iEvent.getByLabel(caloCollection_,calotower);
240 
241  if(calotower.isValid()){
242 
243  for (CaloTowerCollection::const_iterator j = calotower->begin();j !=calotower->end(); j++) {
244  tower_eta = j->eta();
245  tower_phi = j->phi();
246  tower_energyet_e = j->emEt();
247  tower_energyet_h = j->hadEt();
248  tower_energyet = tower_energyet_e + tower_energyet_h;
249 
250  s1 = sin(2.*tower_phi);
251  s2 = cos(2.*tower_phi);
252  s13 = sin(3.*tower_phi);
253  s23 = cos(3.*tower_phi);
254  s14 = sin(4.*tower_phi);
255  s24 = cos(4.*tower_phi);
256  s15 = sin(5.*tower_phi);
257  s25 = cos(5.*tower_phi);
258  s16 = sin(6.*tower_phi);
259  s26 = cos(6.*tower_phi);
260 
261  if(tower_energyet<minet_) continue;
262  if(tower_energyet>maxet_) continue;;
263 
264  rp[etEcal ]->addParticle (tower_energyet_e, s1, s2, tower_eta);
265  rp[etEcalP ]->addParticle (tower_energyet_e, s1, s2, tower_eta);
266  rp[etEcalM ]->addParticle (tower_energyet_e, s1, s2, tower_eta);
267  rp[etHcal ]->addParticle (tower_energyet_h, s1, s2, tower_eta);
268  rp[etHcalP ]->addParticle (tower_energyet_h, s1, s2, tower_eta);
269  rp[etHcalM ]->addParticle (tower_energyet_h, s1, s2, tower_eta);
270  rp[etHF ]->addParticle (tower_energyet, s1, s2, tower_eta);
271  rp[etHFp ]->addParticle (tower_energyet, s1, s2, tower_eta);
272  rp[etHFm ]->addParticle (tower_energyet, s1, s2, tower_eta);
273  rp[etCaloHFP ]->addParticle (tower_energyet, s1, s2, tower_eta);
274  rp[etCaloHFM ]->addParticle (tower_energyet, s1, s2, tower_eta);
275 
276  rp[etHF3 ]->addParticle (tower_energyet, s13, s23, tower_eta);
277  rp[etHFp3 ]->addParticle (tower_energyet, s13, s23, tower_eta);
278  rp[etHFm3 ]->addParticle (tower_energyet, s13, s23, tower_eta);
279 
280  rp[etHF4 ]->addParticle (tower_energyet, s14, s24, tower_eta);
281  rp[etHFp4 ]->addParticle (tower_energyet, s14, s24, tower_eta);
282  rp[etHFm4 ]->addParticle (tower_energyet, s14, s24, tower_eta);
283 
284  rp[etHF5 ]->addParticle (tower_energyet, s15, s25, tower_eta);
285  rp[etHFp5 ]->addParticle (tower_energyet, s15, s25, tower_eta);
286  rp[etHFm5 ]->addParticle (tower_energyet, s15, s25, tower_eta);
287 
288  rp[etHF6 ]->addParticle (tower_energyet, s16, s26, tower_eta);
289  rp[etHFp6 ]->addParticle (tower_energyet, s16, s26, tower_eta);
290  rp[etHFm6 ]->addParticle (tower_energyet, s16, s26, tower_eta);
291  }
292  }
293  //Tracking part
294 
295  double track_eta;
296  double track_phi;
297  double track_pt;
298 #ifdef TRACKCOLLECTION
299  Handle<reco::TrackCollection> tracks;
300  iEvent.getByLabel(trackCollection_, tracks);
301 
302  if(tracks.isValid()){
303  for(reco::TrackCollection::const_iterator j = tracks->begin(); j != tracks->end(); j++){
304 
305  //Find possible collections with command line: edmDumpEventContent *.root
306 #endif
307 #ifdef RECOCHARGEDCANDIDATECOLLECTION
308  edm::Handle<reco::RecoChargedCandidateCollection> trackCollection;
309  iEvent.getByLabel("allMergedPtSplit12Tracks",trackCollection);
310  // iEvent.getByLabel("hiGoodMergedTracks",trackCollection);
311  if(trackCollection.isValid()){
312 
313  const reco::RecoChargedCandidateCollection * tracks = trackCollection.product();
314  for(reco::RecoChargedCandidateCollection::const_iterator j = tracks->begin(); j != tracks->end(); j++){
315 #endif
316  edm::Handle<reco::VertexCollection> vertex;
317  iEvent.getByLabel(vtxCollection_, vertex);
318 
319 // find the vertex point and error
320 
321  math::XYZPoint vtxPoint(0.0,0.0,0.0);
322  double vzErr =0.0, vxErr=0.0, vyErr=0.0;
323  if(vertex->size()>0) {
324  vtxPoint=vertex->begin()->position();
325  vzErr=vertex->begin()->zError();
326  vxErr=vertex->begin()->xError();
327  vyErr=vertex->begin()->yError();
328  }
329  bool accepted = true;
330  bool isPixel = false;
331  // determine if the track is a pixel track
332  if ( j->numberOfValidHits() < 7 ) isPixel = true;
333 
334  // determine the vertex significance
335  double d0=0.0, dz=0.0, d0sigma=0.0, dzsigma=0.0;
336  d0 = -1.*j->dxy(vtxPoint);
337  dz = j->dz(vtxPoint);
338  d0sigma = sqrt(j->d0Error()*j->d0Error()+vxErr*vyErr);
339  dzsigma = sqrt(j->dzError()*j->dzError()+vzErr*vzErr);
340 
341  // cuts for pixel tracks
342  if( isPixel )
343  {
344  // dz significance cut
345  if ( fabs(dz/dzsigma) > dzerr_ ) accepted = false;
346 
347  // chi2/ndof cut
348  if ( j->normalizedChi2() > chi2_ ) accepted = false;
349  }
350 
351  // cuts for full tracks
352  if ( ! isPixel)
353  {
354  // dz and d0 significance cuts
355  if ( fabs(dz/dzsigma) > 3 ) accepted = false;
356  if ( fabs(d0/d0sigma) > 3 ) accepted = false;
357 
358  // pt resolution cut
359  if ( j->ptError()/j->pt() > 0.05 ) accepted = false;
360 
361  // number of valid hits cut
362  if ( j->numberOfValidHits() < 12 ) accepted = false;
363  }
364  if( accepted ) {
365  track_eta = j->eta();
366  track_phi = j->phi();
367  track_pt = j->pt();
368  double s =sin(2*track_phi);
369  double c =cos(2*track_phi);
370  double s3 =sin(3*track_phi);
371  double c3 =cos(3*track_phi);
372  double s4 =sin(4*track_phi);
373  double c4 =cos(4*track_phi);
374  double s5 =sin(5*track_phi);
375  double c5 =cos(5*track_phi);
376  double s6 =sin(6*track_phi);
377  double c6 =cos(6*track_phi);
378  double w = 1;
379  if(useTrackPtWeight_) {
380  w = track_pt;
381  if(w>2.5) w=2.0; //v2 starts decreasing above ~2.5 GeV/c
382  }
383  if(track_pt<minpt_) continue;
384  if(track_pt>maxpt_) continue;
385  rp[EvtPlaneFromTracksMidEta]->addParticle(w,s,c,track_eta);
386  rp[EvtPTracksPosEtaGap]->addParticle(w,s,c,track_eta);
387  rp[EvtPTracksNegEtaGap]->addParticle(w,s,c,track_eta);
388  rp[EPTracksMid3]->addParticle(w,s3,c3,track_eta);
389  rp[EPTracksPos3]->addParticle(w,s3,c3,track_eta);
390  rp[EPTracksNeg3]->addParticle(w,s3,c3,track_eta);
391  rp[EPTracksMid4]->addParticle(w,s4,c4,track_eta);
392  rp[EPTracksPos4]->addParticle(w,s4,c4,track_eta);
393  rp[EPTracksNeg4]->addParticle(w,s4,c4,track_eta);
394  rp[EPTracksMid5]->addParticle(w,s5,c5,track_eta);
395  rp[EPTracksPos5]->addParticle(w,s5,c5,track_eta);
396  rp[EPTracksNeg5]->addParticle(w,s5,c5,track_eta);
397  rp[EPTracksMid6]->addParticle(w,s6,c6,track_eta);
398  rp[EPTracksPos6]->addParticle(w,s6,c6,track_eta);
399  rp[EPTracksNeg6]->addParticle(w,s6,c6,track_eta);
400 
401  }
402 
403  }
404  }
405  }
406  std::auto_ptr<EvtPlaneCollection> evtplaneOutput(new EvtPlaneCollection);
407 
408  EvtPlane *ep[NumEPNames];
409 
410  double ang=-10;
411  double sv = 0;
412  double cv = 0;
413 
414  for(int i = 0; i<NumEPNames; i++) {
415  rp[i]->getAngle(ang,sv,cv);
416  if(storeNames_) ep[i] = new EvtPlane(ang,sv,cv,EPNames[i]);
417  else ep[i] = new EvtPlane(ang,sv,cv,"");
418  }
419  if(useTrack_) {
420  for(int i = 0; i<9; i++) {
421  evtplaneOutput->push_back(*ep[i]);
422  }
423  }
424  for(int i = 9; i<NumEPNames; i++) {
425  if(useECAL_ && !useHCAL_) {
426  if(EPNames[i].rfind("Ecal")!=string::npos) {
427  evtplaneOutput->push_back(*ep[i]);
428  }
429  } else if (useHCAL_ && !useECAL_) {
430  if(EPNames[i].rfind("Hcal")!=string::npos) {
431  evtplaneOutput->push_back(*ep[i]);
432  }
433  }else if (useECAL_ && useHCAL_) {
434  evtplaneOutput->push_back(*ep[i]);
435  }
436  }
437 
438  iEvent.put(evtplaneOutput, "recoLevel");
439  // storeNames_ = 0;
440 }
EvtPlaneProducer::GenPlane::getAngle
double getAngle(double &ang, double &sv, double &cv)
Definition: EvtPlaneProducer.cc:101
i
int i
Definition: DBlmapReader.cc:9
funct::sin
Sin< T >::type sin(const T &t)
Definition: Sin.h:22
edm::SortedCollection< CaloTower >::const_iterator
std::vector< CaloTower >::const_iterator const_iterator
Definition: SortedCollection.h:80
reco::VertexCollection
std::vector< Vertex > VertexCollection
collection of Vertex objects
Definition: VertexFwd.h:9
dt_dqm_sourceclient_common_cff.reco
tuple reco
Definition: dt_dqm_sourceclient_common_cff.py:105
hi::EPNames
const std::string EPNames[]
Definition: HiEvtPlaneList.h:103
indexGen.s2
tuple s2
Definition: indexGen.py:106
EvtPlaneProducer::rp
GenPlane * rp[NumEPNames]
Definition: EvtPlaneProducer.cc:130
edm::Event::put
OrphanHandle< PROD > put(std::auto_ptr< PROD > product)
Put a new product.
Definition: Event.h:85
mathSSE::sqrt
T sqrt(T t)
Definition: SSEVec.h:46
funct::cos
Cos< T >::type cos(const T &t)
Definition: Cos.h:22
reco::EvtPlaneCollection
std::vector< EvtPlane > EvtPlaneCollection
Definition: EvtPlane.h:34
j
int j
Definition: DBlmapReader.cc:9
edm::HandleBase::isValid
bool isValid() const
Definition: HandleBase.h:76
edm::Event::getByLabel
bool getByLabel(InputTag const &tag, Handle< PROD > &result) const
Definition: Event.h:356
reco::RecoChargedCandidateCollection
std::vector< RecoChargedCandidate > RecoChargedCandidateCollection
collectin of RecoChargedCandidate objects
Definition: RecoChargedCandidateFwd.h:9
EvtPlaneProducer::GenPlane::addParticle
void addParticle(double w, double s, double c, double eta)
Definition: EvtPlaneProducer.cc:88
testEve_cfg.tracks
tuple tracks
Definition: testEve_cfg.py:39
math::XYZPoint
XYZPointD XYZPoint
point in space with cartesian internal representation
Definition: Point3D.h:13
edm::Handle::product
T const * product() const
Definition: Handle.h:74
EvtPlaneProducer::caloCollection_
edm::InputTag caloCollection_
Definition: EvtPlaneProducer.cc:138
EvtPlaneProducer::vtxCollection_
edm::InputTag vtxCollection_
Definition: EvtPlaneProducer.cc:137
hi::NumEPNames
static const int NumEPNames
Definition: HiEvtPlaneList.h:198
EvtPlaneProducer::trackCollection_
edm::InputTag trackCollection_
Definition: EvtPlaneProducer.cc:139
reset
void reset(double vett[256])
Definition: TPedValues.cc:11
w
T w() const
Definition: newBasic3DVector.h:234

Member Data Documentation

edm::InputTag EvtPlaneProducer::caloCollection_
private

Definition at line 138 of file EvtPlaneProducer.cc.

double EvtPlaneProducer::chi2_
private

Definition at line 151 of file EvtPlaneProducer.cc.

double EvtPlaneProducer::dzerr_
private

Definition at line 150 of file EvtPlaneProducer.cc.

double EvtPlaneProducer::maxet_
private

Definition at line 145 of file EvtPlaneProducer.cc.

double EvtPlaneProducer::maxpt_
private

Definition at line 147 of file EvtPlaneProducer.cc.

double EvtPlaneProducer::maxvtx_
private

Definition at line 149 of file EvtPlaneProducer.cc.

double EvtPlaneProducer::minet_
private

Definition at line 144 of file EvtPlaneProducer.cc.

double EvtPlaneProducer::minpt_
private

Definition at line 146 of file EvtPlaneProducer.cc.

double EvtPlaneProducer::minvtx_
private

Definition at line 148 of file EvtPlaneProducer.cc.

GenPlane* EvtPlaneProducer::rp[NumEPNames]
private

Definition at line 130 of file EvtPlaneProducer.cc.

bool EvtPlaneProducer::storeNames_
private

Definition at line 153 of file EvtPlaneProducer.cc.

edm::InputTag EvtPlaneProducer::trackCollection_
private

Definition at line 139 of file EvtPlaneProducer.cc.

bool EvtPlaneProducer::useECAL_
private

Definition at line 140 of file EvtPlaneProducer.cc.

bool EvtPlaneProducer::useHCAL_
private

Definition at line 141 of file EvtPlaneProducer.cc.

bool EvtPlaneProducer::useTrack_
private

Definition at line 142 of file EvtPlaneProducer.cc.

bool EvtPlaneProducer::useTrackPtWeight_
private

Definition at line 143 of file EvtPlaneProducer.cc.

edm::InputTag EvtPlaneProducer::vtxCollection_
private

Definition at line 137 of file EvtPlaneProducer.cc.

Generated for CMSSW Reference Manual by doxygen 1.8.5