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

#include <ElectronSeedAnalyzer.h>

Inheritance diagram for ElectronSeedAnalyzer:
edm::one::EDAnalyzer<> edm::one::EDAnalyzerBase edm::EDConsumerBase

Public Member Functions

void analyze (const edm::Event &, const edm::EventSetup &) override
 
void beginJob () override
 
 ElectronSeedAnalyzer (const edm::ParameterSet &conf)
 
void endJob () override
 
 ~ElectronSeedAnalyzer () override
 
- Public Member Functions inherited from edm::one::EDAnalyzer<>
 EDAnalyzer ()=default
 
SerialTaskQueue * globalLuminosityBlocksQueue () final
 
SerialTaskQueue * globalRunsQueue () final
 
bool wantsGlobalLuminosityBlocks () const final
 
bool wantsGlobalRuns () const final
 
- Public Member Functions inherited from edm::one::EDAnalyzerBase
void callWhenNewProductsRegistered (std::function< void(BranchDescription const &)> const &func)
 
 EDAnalyzerBase ()
 
ModuleDescription const & moduleDescription () const
 
bool wantsStreamLuminosityBlocks () const
 
bool wantsStreamRuns () const
 
 ~EDAnalyzerBase () override
 
- Public Member Functions inherited from edm::EDConsumerBase
std::vector< ConsumesInfo > consumesInfo () const
 
void convertCurrentProcessAlias (std::string const &processName)
 Convert "@currentProcess" in InputTag process names to the actual current process name. More...
 
 EDConsumerBase ()
 
 EDConsumerBase (EDConsumerBase &&)=default
 
 EDConsumerBase (EDConsumerBase const &)=delete
 
ESProxyIndex const * esGetTokenIndices (edm::Transition iTrans) const
 
ProductResolverIndexAndSkipBit indexFrom (EDGetToken, BranchType, TypeID const &) const
 
void itemsMayGet (BranchType, std::vector< ProductResolverIndexAndSkipBit > &) const
 
void itemsToGet (BranchType, std::vector< ProductResolverIndexAndSkipBit > &) const
 
std::vector< ProductResolverIndexAndSkipBit > const & itemsToGetFrom (BranchType iType) const
 
void labelsForToken (EDGetToken iToken, Labels &oLabels) const
 
void modulesWhoseProductsAreConsumed (std::vector< ModuleDescription const * > &modules, ProductRegistry const &preg, std::map< std::string, ModuleDescription const * > const &labelsToDesc, std::string const &processName) const
 
EDConsumerBase & operator= (EDConsumerBase &&)=default
 
EDConsumerBase const & operator= (EDConsumerBase const &)=delete
 
bool registeredToConsume (ProductResolverIndex, bool, BranchType) const
 
bool registeredToConsumeMany (TypeID const &, BranchType) const
 
ProductResolverIndexAndSkipBit uncheckedIndexFrom (EDGetToken) const
 
void updateLookup (BranchType iBranchType, ProductResolverIndexHelper const &, bool iPrefetchMayGet)
 
void updateLookup (eventsetup::ESRecordsToProxyIndices const &)
 
virtual ~EDConsumerBase () noexcept(false)
 

Private Attributes

edm::InputTag beamSpot_
 
TH1F * histecaldriveneMCmatched_
 
TH1F * histecaldrivenetaMCmatched_
 
TH1F * histecaldrivenptMCmatched_
 
TH1F * histeclu_
 
TH1F * histeffeta_
 
TH1F * histeffpt_
 
TH1F * histeMC_
 
TH1F * histeMCmatched_
 
TH1F * histeoverp_
 
TH1F * histeta_
 
TH1F * histetaclu_
 
TH1F * histetaMC_
 
TH1F * histetaMCmatched_
 
TH1F * histetclu_
 
TFile * histfile_
 
TH1I * histnbclus_
 
TH1I * histnbseeds_
 
TH1I * histnrseeds_
 
TH1F * histp_
 
TH1F * histpt_
 
TH1F * histptMC_
 
TH1F * histptMCmatched_
 
TH1F * histq_
 
TH1F * histtrackerdriveneMCmatched_
 
TH1F * histtrackerdrivenetaMCmatched_
 
TH1F * histtrackerdrivenptMCmatched_
 
edm::InputTag inputCollection_
 
float mcEnergy [10]
 
float mcEta [10]
 
float mcPhi [10]
 
float mcPt [10]
 
float mcQ [10]
 
float seedDphi1 [10]
 
float seedDphi2 [10]
 
float seedDrz1 [10]
 
float seedDrz2 [10]
 
float seedEta [10]
 
int seedLayer1 [10]
 
int seedLayer2 [10]
 
float seedMomentum [10]
 
float seedPhi [10]
 
float seedPhi1 [10]
 
float seedPhi2 [10]
 
float seedPt [10]
 
float seedQ [10]
 
float seedRz1 [10]
 
float seedRz2 [10]
 
int seedSide1 [10]
 
int seedSide2 [10]
 
int seedSubdet1 [10]
 
int seedSubdet2 [10]
 
float superclusterEnergy [10]
 
float superclusterEt [10]
 
float superclusterEta [10]
 
float superclusterPhi [10]
 
TrajectoryStateTransform transformer_
 
TTree * tree_
 

Additional Inherited Members

- Public Types inherited from edm::one::EDAnalyzerBase
typedef EDAnalyzerBase ModuleType
 
- Public Types inherited from edm::EDConsumerBase
typedef ProductLabels Labels
 
- Static Public Member Functions inherited from edm::one::EDAnalyzerBase
static const std::string & baseType ()
 
static void fillDescriptions (ConfigurationDescriptions &descriptions)
 
static void prevalidate (ConfigurationDescriptions &descriptions)
 
- Protected Member Functions inherited from edm::EDConsumerBase
EDGetToken consumes (const TypeToGet &id, edm::InputTag const &tag)
 
template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType > consumes (edm::InputTag const &tag)
 
template<BranchType B>
EDGetToken consumes (TypeToGet const &id, edm::InputTag const &tag)
 
ConsumesCollector consumesCollector ()
 Use a ConsumesCollector to gather consumes information from helper functions. More...
 
template<typename ProductType , BranchType B = InEvent>
void consumesMany ()
 
void consumesMany (const TypeToGet &id)
 
template<BranchType B>
void consumesMany (const TypeToGet &id)
 
template<typename ESProduct , typename ESRecord , Transition Tr = Transition::Event>
auto esConsumes ()
 
template<typename ESProduct , typename ESRecord , Transition Tr = Transition::Event>
auto esConsumes (ESInputTag const &tag)
 
EDGetToken mayConsume (const TypeToGet &id, edm::InputTag const &tag)
 
template<BranchType B>
EDGetToken mayConsume (const TypeToGet &id, edm::InputTag const &tag)
 
template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType > mayConsume (edm::InputTag const &tag)
 

Detailed Description

Definition at line 33 of file ElectronSeedAnalyzer.h.

Constructor & Destructor Documentation

◆ ElectronSeedAnalyzer()

ElectronSeedAnalyzer::ElectronSeedAnalyzer ( const edm::ParameterSet &  conf)
explicit

Definition at line 56 of file ElectronSeedAnalyzer.cc.

57  : beamSpot_(conf.getParameter<edm::InputTag>("beamSpot")) {
58  inputCollection_ = conf.getParameter<edm::InputTag>("inputCollection");
59  histfile_ = new TFile("electronpixelseeds.root", "RECREATE");
60 }

References edm::ParameterSet::getParameter(), histfile_, and inputCollection_.

◆ ~ElectronSeedAnalyzer()

ElectronSeedAnalyzer::~ElectronSeedAnalyzer ( )
override

Definition at line 173 of file ElectronSeedAnalyzer.cc.

173  {
174  // do anything here that needs to be done at desctruction time
175  // (e.g. close files, deallocate resources etc.)
176  //tree_->Print();
177  histfile_->Write();
178  histeMC_->Write();
179  histfile_->Close();
180 }

References histeMC_, and histfile_.

Member Function Documentation

◆ analyze()

void ElectronSeedAnalyzer::analyze ( const edm::Event &  e,
const edm::EventSetup &  iSetup 
)
overridevirtual

Implements edm::one::EDAnalyzerBase.

Definition at line 182 of file ElectronSeedAnalyzer.cc.

182  {
183  //Retrieve tracker topology from geometry
184  edm::ESHandle<TrackerTopology> tTopo;
185  iSetup.get<TrackerTopologyRcd>().get(tTopo);
186 
187  edm::ESHandle<TrackerGeometry> pDD;
188  edm::ESHandle<MagneticField> theMagField;
189  iSetup.get<TrackerDigiGeometryRecord>().get(pDD);
190  iSetup.get<IdealMagneticFieldRecord>().get(theMagField);
191 
192  // rereads the seeds for test purposes
193  typedef edm::OwnVector<TrackingRecHit> recHitContainer;
194  typedef recHitContainer::const_iterator const_iterator;
195  typedef std::pair<const_iterator, const_iterator> range;
196 
197  // get beam spot
198  edm::Handle<reco::BeamSpot> theBeamSpot;
199  e.getByLabel(beamSpot_, theBeamSpot);
200 
201  // get seeds
202 
203  edm::Handle<ElectronSeedCollection> elSeeds;
204  e.getByLabel(inputCollection_, elSeeds);
205  edm::LogInfo("") << "\n\n =================> Treating event " << e.id() << " Number of seeds "
206  << elSeeds.product()->size();
207  int is = 0;
208 
209  float mass = .000511; // electron propagation
210  PropagatorWithMaterial *prop1stLayer = new PropagatorWithMaterial(oppositeToMomentum, mass, &(*theMagField));
211  PropagatorWithMaterial *prop2ndLayer = new PropagatorWithMaterial(alongMomentum, mass, &(*theMagField));
212 
213  float dphi1 = 0., dphi2 = 0., drz1 = 0., drz2 = 0.;
214  float phi1 = 0., phi2 = 0., rz1 = 0., rz2 = 0.;
215 
216  for (ElectronSeedCollection::const_iterator MyS = (*elSeeds).begin(); MyS != (*elSeeds).end(); ++MyS) {
217  LogDebug("") << "\nSeed nr " << is << ": ";
218  range r = (*MyS).recHits();
219  LogDebug("") << " Number of RecHits= " << (*MyS).nHits();
220  const GeomDet *det1 = nullptr;
221  const GeomDet *det2 = nullptr;
222 
223  TrajectorySeed::const_iterator it = r.first;
224  DetId id1 = (*it).geographicalId();
225  det1 = pDD->idToDet(id1);
226  LogDebug("") << " First hit local x,y,z " << (*it).localPosition() << " det " << id1.det() << " subdet "
227  << id1.subdetId();
228  LogDebug("") << " First hit global " << det1->toGlobal((*it).localPosition());
229  //std::cout <<" First hit local x,y,z "<<(*it).localPosition()<<" det "<<id1.det()<<" subdet "<<id1.subdetId()<< std::endl;
230  //std::cout <<" First hit global "<<det1->toGlobal((*it).localPosition())<< std::endl;
231  it++;
232  DetId id2 = (*it).geographicalId();
233  det2 = pDD->idToDet(id2);
234  LogDebug("") << " Second hit local x,y,z " << (*it).localPosition() << " det " << id2.det() << " subdet "
235  << id2.subdetId();
236  LogDebug("") << " Second hit global " << det2->toGlobal((*it).localPosition());
237  //std::cout <<" Second hit local x,y,z "<<(*it).localPosition()<<" det "<<id2.det()<<" subdet "<<id2.subdetId()<< std::endl;
238  //std::cout <<" Second hit global "<<det2->toGlobal((*it).localPosition()) << std::endl;
239 
240  // state on last det
241  const GeomDet *det = nullptr;
242  for (TrackingRecHitCollection::const_iterator rhits = r.first; rhits != r.second; rhits++)
243  det = pDD->idToDet(((*rhits)).geographicalId());
244  TrajectoryStateOnSurface t =
245  trajectoryStateTransform::transientState((*MyS).startingState(), &(det->surface()), &(*theMagField));
246 
247  // debug
248 
249  LogDebug("") << " ElectronSeed outermost state position: " << t.globalPosition();
250  LogDebug("") << " ElectronSeed outermost state momentum: " << t.globalMomentum();
251  edm::RefToBase<CaloCluster> caloCluster = (*MyS).caloCluster();
252  if (caloCluster.isNull())
253  continue;
254  edm::Ref<SuperClusterCollection> theClus = caloCluster.castTo<SuperClusterRef>();
255  LogDebug("") << " ElectronSeed superCluster energy: " << theClus->energy() << ", position: " << theClus->position();
256  LogDebug("") << " ElectronSeed outermost state Pt: " << t.globalMomentum().perp();
257  LogDebug("") << " ElectronSeed supercluster Et: "
258  << theClus->energy() * sin(2. * atan(exp(-theClus->position().eta())));
259  LogDebug("") << " ElectronSeed outermost momentum direction eta: " << t.globalMomentum().eta();
260  LogDebug("") << " ElectronSeed supercluster eta: " << theClus->position().eta();
261  LogDebug("") << " ElectronSeed seed charge: " << (*MyS).getCharge();
262  LogDebug("") << " ElectronSeed E/p: " << theClus->energy() / t.globalMomentum().mag();
263 
264  // retreive SC and compute distances between hit position and prediction the same
265  // way as in the PixelHitMatcher
266 
267  // inputs are charge, cluster position, vertex position, cluster energy and B field
268  int charge = int((*MyS).getCharge());
269  GlobalPoint xmeas(theClus->position().x(), theClus->position().y(), theClus->position().z());
270  GlobalPoint vprim(theBeamSpot->position().x(), theBeamSpot->position().y(), theBeamSpot->position().z());
271  float energy = theClus->energy();
272 
273  auto fts = trackingTools::ftsFromVertexToPoint(*theMagField, xmeas, vprim, energy, charge);
274  //std::cout << "[PixelHitMatcher::compatibleSeeds] fts position, momentum " <<
275  // fts.parameters().position() << " " << fts.parameters().momentum() << std::endl;
276 
277  PerpendicularBoundPlaneBuilder bpb;
278  TrajectoryStateOnSurface tsos(fts, *bpb(fts.position(), fts.momentum()));
279 
280  // TrajectorySeed::range r=(*seeds.product())[i].recHits();
281  // TrajectorySeed::range r=(*seeds)[i].recHits();
282 
283  // first Hit
284  it = r.first;
285  DetId id = (*it).geographicalId();
286  const GeomDet *geomdet = pDD->idToDet((*it).geographicalId());
287  LocalPoint lp = (*it).localPosition();
288  GlobalPoint hitPos = geomdet->surface().toGlobal(lp);
289 
290  TrajectoryStateOnSurface tsos1;
291  tsos1 = prop1stLayer->propagate(tsos, geomdet->surface());
292 
293  if (tsos1.isValid()) {
294  //UB add test on phidiff
295  float SCl_phi = xmeas.phi();
296  float localDphi = SCl_phi - hitPos.phi();
297  if (localDphi > CLHEP::pi)
298  localDphi -= (2 * CLHEP::pi);
299  if (localDphi < -CLHEP::pi)
300  localDphi += (2 * CLHEP::pi);
301  if (std::abs(localDphi) > 2.5)
302  continue;
303 
304  phi1 = hitPos.phi();
305  dphi1 = hitPos.phi() - tsos1.globalPosition().phi();
306  rz1 = hitPos.perp();
307  drz1 = hitPos.perp() - tsos1.globalPosition().perp();
308  if (id.subdetId() % 2 == 1) {
309  drz1 = hitPos.z() - tsos1.globalPosition().z();
310  rz1 = hitPos.z();
311  }
312 
313  // now second Hit
314  it++;
315  DetId id2 = (*it).geographicalId();
316  const GeomDet *geomdet2 = pDD->idToDet((*it).geographicalId());
317  TrajectoryStateOnSurface tsos2;
318 
319  // compute the z vertex from the cluster point and the found pixel hit
320  double pxHit1z = hitPos.z();
321  double pxHit1x = hitPos.x();
322  double pxHit1y = hitPos.y();
323  double r1diff = (pxHit1x - vprim.x()) * (pxHit1x - vprim.x()) + (pxHit1y - vprim.y()) * (pxHit1y - vprim.y());
324  r1diff = sqrt(r1diff);
325  double r2diff = (xmeas.x() - pxHit1x) * (xmeas.x() - pxHit1x) + (xmeas.y() - pxHit1y) * (xmeas.y() - pxHit1y);
326  r2diff = sqrt(r2diff);
327  double zVertexPred = pxHit1z - r1diff * (xmeas.z() - pxHit1z) / r2diff;
328 
329  GlobalPoint vertexPred(vprim.x(), vprim.y(), zVertexPred);
330 
331  auto fts2 = trackingTools::ftsFromVertexToPoint(*theMagField, hitPos, vertexPred, energy, charge);
332  tsos2 = prop2ndLayer->propagate(fts2, geomdet2->surface());
333 
334  if (tsos2.isValid()) {
335  LocalPoint lp2 = (*it).localPosition();
336  GlobalPoint hitPos2 = geomdet2->surface().toGlobal(lp2);
337  phi2 = hitPos2.phi();
338  dphi2 = hitPos2.phi() - tsos2.globalPosition().phi();
339  rz2 = hitPos2.perp();
340  drz2 = hitPos2.perp() - tsos2.globalPosition().perp();
341  if (id2.subdetId() % 2 == 1) {
342  rz2 = hitPos2.z();
343  drz2 = hitPos2.z() - tsos2.globalPosition().z();
344  }
345  }
346  }
347 
348  // fill the tree and histos
349 
350  histpt_->Fill(t.globalMomentum().perp());
351  histetclu_->Fill(theClus->energy() * sin(2. * atan(exp(-theClus->position().eta()))));
352  histeta_->Fill(t.globalMomentum().eta());
353  histetaclu_->Fill(theClus->position().eta());
354  histq_->Fill((*MyS).getCharge());
355  histeoverp_->Fill(theClus->energy() / t.globalMomentum().mag());
356 
357  if (is < 10) {
358  superclusterEnergy[is] = theClus->energy();
359  superclusterEta[is] = theClus->position().eta();
360  superclusterPhi[is] = theClus->position().phi();
361  superclusterEt[is] = theClus->energy() * sin(2. * atan(exp(-theClus->position().eta())));
362  seedMomentum[is] = t.globalMomentum().mag();
363  seedEta[is] = t.globalMomentum().eta();
364  seedPhi[is] = t.globalMomentum().phi();
365  seedPt[is] = t.globalMomentum().perp();
366  seedQ[is] = (*MyS).getCharge();
367  seedSubdet1[is] = id1.subdetId();
368  seedLayer1[is] = tTopo->layer(id1);
369  seedSide1[is] = tTopo->side(id1);
370  seedPhi1[is] = phi1;
371  seedRz1[is] = rz1;
372  seedDphi1[is] = dphi1;
373  seedDrz1[is] = drz1;
374  seedSubdet2[is] = id2.subdetId();
375  seedLayer2[is] = tTopo->layer(id2);
376  seedSide2[is] = tTopo->side(id2);
377  seedDphi2[is] = dphi2;
378  seedDrz2[is] = drz2;
379  seedPhi2[is] = phi2;
380  seedRz2[is] = rz2;
381  }
382 
383  is++;
384  }
385 
386  histnbseeds_->Fill(elSeeds.product()->size());
387 
388  // get input clusters
389 
390  edm::Handle<SuperClusterCollection> clusters;
391  //CC to be changed according to supercluster input
392  e.getByLabel("correctedHybridSuperClusters", clusters);
393  histnbclus_->Fill(clusters.product()->size());
394  if (!clusters.product()->empty())
395  histnrseeds_->Fill(elSeeds.product()->size());
396  // get MC information
397 
398  edm::Handle<edm::HepMCProduct> HepMCEvt;
399  // this one is empty branch in current test files
400  //e.getByLabel("generatorSmeared", "", HepMCEvt);
401  //e.getByLabel("source", "", HepMCEvt);
402  e.getByLabel("generatorSmeared", "", HepMCEvt);
403 
404  const HepMC::GenEvent *MCEvt = HepMCEvt->GetEvent();
405  HepMC::GenParticle *genPc = nullptr;
406  HepMC::FourVector pAssSim;
407  int ip = 0;
408  for (HepMC::GenEvent::particle_const_iterator partIter = MCEvt->particles_begin(); partIter != MCEvt->particles_end();
409  ++partIter) {
410  for (HepMC::GenEvent::vertex_const_iterator vertIter = MCEvt->vertices_begin(); vertIter != MCEvt->vertices_end();
411  ++vertIter) {
412  // CLHEP::HepLorentzVector creation = (*partIter)->CreationVertex();
413  HepMC::GenVertex *creation = (*partIter)->production_vertex();
414  // CLHEP::HepLorentzVector momentum = (*partIter)->Momentum();
415  HepMC::FourVector momentum = (*partIter)->momentum();
416  // HepPDT::ParticleID id = (*partIter)->particleID(); // electrons and positrons are 11 and -11
417  int id = (*partIter)->pdg_id(); // electrons and positrons are 11 and -11
418  LogDebug("") << "MC particle id " << id << ", creationVertex " << (*creation) << " cm, initialMomentum "
419  << momentum.rho() << " GeV/c" << std::endl;
420 
421  if (id == 11 || id == -11) {
422  // single primary electrons or electrons from Zs or Ws
423  HepMC::GenParticle *mother = nullptr;
424  if ((*partIter)->production_vertex()) {
425  if ((*partIter)->production_vertex()->particles_begin(HepMC::parents) !=
426  (*partIter)->production_vertex()->particles_end(HepMC::parents))
427  mother = *((*partIter)->production_vertex()->particles_begin(HepMC::parents));
428  }
429  if (((mother == nullptr) || ((mother != nullptr) && (mother->pdg_id() == 23)) ||
430  ((mother != nullptr) && (mother->pdg_id() == 32)) ||
431  ((mother != nullptr) && (std::abs(mother->pdg_id()) == 24)))) {
432  genPc = (*partIter);
433  pAssSim = genPc->momentum();
434 
435  // EWK fiducial
436  //if (pAssSim.perp()> 100. || std::abs(pAssSim.eta())> 2.5) continue;
437  //if (pAssSim.perp()< 20. || (std::abs(pAssSim.eta())> 1.4442 && std::abs(pAssSim.eta())< 1.56) || std::abs(pAssSim.eta())> 2.5) continue;
438  // reconstruction fiducial
439  //if (pAssSim.perp()< 5. || std::abs(pAssSim.eta())> 2.5) continue;
440  if (std::abs(pAssSim.eta()) > 2.5)
441  continue;
442 
443  histptMC_->Fill(pAssSim.perp());
444  histetaMC_->Fill(pAssSim.eta());
445  histeMC_->Fill(pAssSim.rho());
446 
447  // looking for the best matching gsf electron
448  bool okSeedFound = false;
449  double seedOkRatio = 999999.;
450 
451  // find best matched seed
452  reco::ElectronSeed bestElectronSeed;
453  for (ElectronSeedCollection::const_iterator gsfIter = (*elSeeds).begin(); gsfIter != (*elSeeds).end();
454  ++gsfIter) {
455  range r = gsfIter->recHits();
456  const GeomDet *det = nullptr;
457  for (TrackingRecHitCollection::const_iterator rhits = r.first; rhits != r.second; rhits++)
458  det = pDD->idToDet(((*rhits)).geographicalId());
459  TrajectoryStateOnSurface t =
460  trajectoryStateTransform::transientState(gsfIter->startingState(), &(det->surface()), &(*theMagField));
461 
462  float eta = t.globalMomentum().eta();
463  float phi = t.globalMomentum().phi();
464  float p = t.globalMomentum().mag();
465  double dphi = phi - pAssSim.phi();
466  if (std::abs(dphi) > CLHEP::pi)
467  dphi = dphi < 0 ? (CLHEP::twopi) + dphi : dphi - CLHEP::twopi;
468  double deltaR = sqrt(std::pow((eta - pAssSim.eta()), 2) + std::pow(dphi, 2));
469  if (deltaR < 0.15) {
470  // if ( deltaR < 0.3 ){
471  //if ( (genPc->pdg_id() == 11) && (gsfIter->charge() < 0.) || (genPc->pdg_id() == -11) &&
472  //(gsfIter->charge() > 0.) ){
473  double tmpSeedRatio = p / pAssSim.t();
474  if (std::abs(tmpSeedRatio - 1) < std::abs(seedOkRatio - 1)) {
475  seedOkRatio = tmpSeedRatio;
476  bestElectronSeed = *gsfIter;
477  okSeedFound = true;
478  }
479  //}
480  }
481  } // loop over rec ele to look for the best one
482 
483  // analysis when the mc track is found
484  if (okSeedFound) {
485  histptMCmatched_->Fill(pAssSim.perp());
486  histetaMCmatched_->Fill(pAssSim.eta());
487  histeMCmatched_->Fill(pAssSim.rho());
488  if (ip < 10) {
489  mcEnergy[ip] = pAssSim.rho();
490  mcEta[ip] = pAssSim.eta();
491  mcPhi[ip] = pAssSim.phi();
492  mcPt[ip] = pAssSim.perp();
493  mcQ[ip] = ((id == 11) ? -1. : +1.);
494  }
495  }
496 
497  // efficiency for ecal driven only
498  okSeedFound = false;
499  seedOkRatio = 999999.;
500 
501  // find best matched seed
502  for (ElectronSeedCollection::const_iterator gsfIter = (*elSeeds).begin(); gsfIter != (*elSeeds).end();
503  ++gsfIter) {
504  range r = gsfIter->recHits();
505  const GeomDet *det = nullptr;
506  for (TrackingRecHitCollection::const_iterator rhits = r.first; rhits != r.second; rhits++)
507  det = pDD->idToDet(((*rhits)).geographicalId());
508  TrajectoryStateOnSurface t =
509  trajectoryStateTransform::transientState(gsfIter->startingState(), &(det->surface()), &(*theMagField));
510 
511  float eta = t.globalMomentum().eta();
512  float phi = t.globalMomentum().phi();
513  float p = t.globalMomentum().mag();
514  double dphi = phi - pAssSim.phi();
515  if (std::abs(dphi) > CLHEP::pi)
516  dphi = dphi < 0 ? (CLHEP::twopi) + dphi : dphi - CLHEP::twopi;
517  double deltaR = sqrt(std::pow((eta - pAssSim.eta()), 2) + std::pow(dphi, 2));
518  if (gsfIter->isEcalDriven()) {
519  if (deltaR < 0.15) {
520  // if ( deltaR < 0.3 ){
521  //if ( (genPc->pdg_id() == 11) && (gsfIter->charge() < 0.) || (genPc->pdg_id() == -11) &&
522  //(gsfIter->charge() > 0.) ){
523  double tmpSeedRatio = p / pAssSim.t();
524  if (std::abs(tmpSeedRatio - 1) < std::abs(seedOkRatio - 1)) {
525  seedOkRatio = tmpSeedRatio;
526  bestElectronSeed = *gsfIter;
527  okSeedFound = true;
528  }
529  //}
530  }
531  } // end if ecal driven
532  } // loop over rec ele to look for the best one
533 
534  // analysis when the mc track is found
535  if (okSeedFound) {
536  histecaldrivenptMCmatched_->Fill(pAssSim.perp());
537  histecaldrivenetaMCmatched_->Fill(pAssSim.eta());
538  histecaldriveneMCmatched_->Fill(pAssSim.rho());
539  }
540 
541  // efficiency for tracker driven only
542  okSeedFound = false;
543  seedOkRatio = 999999.;
544 
545  // find best matched seed
546  for (ElectronSeedCollection::const_iterator gsfIter = (*elSeeds).begin(); gsfIter != (*elSeeds).end();
547  ++gsfIter) {
548  range r = gsfIter->recHits();
549  const GeomDet *det = nullptr;
550  for (TrackingRecHitCollection::const_iterator rhits = r.first; rhits != r.second; rhits++)
551  det = pDD->idToDet(((*rhits)).geographicalId());
552  TrajectoryStateOnSurface t =
553  trajectoryStateTransform::transientState(gsfIter->startingState(), &(det->surface()), &(*theMagField));
554 
555  float eta = t.globalMomentum().eta();
556  float phi = t.globalMomentum().phi();
557  float p = t.globalMomentum().mag();
558  double dphi = phi - pAssSim.phi();
559  if (std::abs(dphi) > CLHEP::pi)
560  dphi = dphi < 0 ? (CLHEP::twopi) + dphi : dphi - CLHEP::twopi;
561  double deltaR = sqrt(std::pow((eta - pAssSim.eta()), 2) + std::pow(dphi, 2));
562  if (gsfIter->isTrackerDriven()) {
563  if (deltaR < 0.15) {
564  // if ( deltaR < 0.3 ){
565  //if ( (genPc->pdg_id() == 11) && (gsfIter->charge() < 0.) || (genPc->pdg_id() == -11) &&
566  //(gsfIter->charge() > 0.) ){
567  double tmpSeedRatio = p / pAssSim.t();
568  if (std::abs(tmpSeedRatio - 1) < std::abs(seedOkRatio - 1)) {
569  seedOkRatio = tmpSeedRatio;
570  bestElectronSeed = *gsfIter;
571  okSeedFound = true;
572  }
573  //}
574  }
575  } // end if ecal driven
576  } // loop over rec ele to look for the best one
577 
578  // analysis when the mc track is found
579  if (okSeedFound) {
580  histtrackerdrivenptMCmatched_->Fill(pAssSim.perp());
581  histtrackerdrivenetaMCmatched_->Fill(pAssSim.eta());
582  histtrackerdriveneMCmatched_->Fill(pAssSim.rho());
583  }
584 
585  } // end if mother W or Z
586 
587  } // end if gen part is electron
588 
589  } // end loop on vertices
590 
591  ip++;
592 
593  } // end loop on gen particles
594 
595  //tree_->Fill();
596 }

References funct::abs(), alongMomentum, beamSpot_, edm::RefToBase< T >::castTo(), ALCARECOTkAlJpsiMuMu_cff::charge, bsc_activity_cfg::clusters, PbPb_ZMuSkimMuonDPG_cff::deltaR, MillePedeFileConverter_cfg::e, HCALHighEnergyHPDFilter_cfi::energy, PVValHelper::eta, JetChargeProducer_cfi::exp, trackingTools::ftsFromVertexToPoint(), GenParticle::GenParticle, edm::EventSetup::get(), get, edm::HepMCProduct::GetEvent(), TrajectoryStateOnSurface::globalPosition(), histecaldriveneMCmatched_, histecaldrivenetaMCmatched_, histecaldrivenptMCmatched_, histeMC_, histeMCmatched_, histeoverp_, histeta_, histetaclu_, histetaMC_, histetaMCmatched_, histetclu_, histnbclus_, histnbseeds_, histnrseeds_, histpt_, histptMC_, histptMCmatched_, histq_, histtrackerdriveneMCmatched_, histtrackerdrivenetaMCmatched_, histtrackerdrivenptMCmatched_, globals_cff::id1, globals_cff::id2, TrackerGeometry::idToDet(), inputCollection_, createfilelist::int, edm::RefToBase< T >::isNull(), TrajectoryStateOnSurface::isValid(), TrackerTopology::layer(), LogDebug, EgHLTOffHistBins_cfi::mass, mcEnergy, mcEta, mcPhi, mcPt, mcQ, oppositeToMomentum, AlCaHLTBitMon_ParallelJobs::p, parents, PV3DBase< T, PVType, FrameType >::perp(), phi, PV3DBase< T, PVType, FrameType >::phi(), pi, reco::BeamSpot::position(), funct::pow(), edm::Handle< T >::product(), PropagatorWithMaterial::propagate(), alignCSCRings::r, FastTimerService_cff::range, seedDphi1, seedDphi2, seedDrz1, seedDrz2, seedEta, seedLayer1, seedLayer2, seedMomentum, seedPhi, seedPhi1, seedPhi2, seedPt, seedQ, seedRz1, seedRz2, seedSide1, seedSide2, seedSubdet1, seedSubdet2, TrackerTopology::side(), funct::sin(), mathSSE::sqrt(), superclusterEnergy, superclusterEt, superclusterEta, superclusterPhi, GeomDet::surface(), OrderedSet::t, GeomDet::toGlobal(), Surface::toGlobal(), trajectoryStateTransform::transientState(), PV3DBase< T, PVType, FrameType >::x(), PV3DBase< T, PVType, FrameType >::y(), and PV3DBase< T, PVType, FrameType >::z().

◆ beginJob()

void ElectronSeedAnalyzer::beginJob ( void  )
overridevirtual

Reimplemented from edm::one::EDAnalyzerBase.

Definition at line 62 of file ElectronSeedAnalyzer.cc.

62  {
63  histfile_->cd();
64  tree_ = new TTree("ElectronSeeds", "ElectronSeed validation ntuple");
65  tree_->Branch("mcEnergy", mcEnergy, "mcEnergy[10]/F");
66  tree_->Branch("mcEta", mcEta, "mcEta[10]/F");
67  tree_->Branch("mcPhi", mcPhi, "mcPhi[10]/F");
68  tree_->Branch("mcPt", mcPt, "mcPt[10]/F");
69  tree_->Branch("mcQ", mcQ, "mcQ[10]/F");
70  tree_->Branch("superclusterEnergy", superclusterEnergy, "superclusterEnergy[10]/F");
71  tree_->Branch("superclusterEta", superclusterEta, "superclusterEta[10]/F");
72  tree_->Branch("superclusterPhi", superclusterPhi, "superclusterPhi[10]/F");
73  tree_->Branch("superclusterEt", superclusterEt, "superclusterEt[10]/F");
74  tree_->Branch("seedMomentum", seedMomentum, "seedMomentum[10]/F");
75  tree_->Branch("seedEta", seedEta, "seedEta[10]/F");
76  tree_->Branch("seedPhi", seedPhi, "seedPhi[10]/F");
77  tree_->Branch("seedPt", seedPt, "seedPt[10]/F");
78  tree_->Branch("seedQ", seedQ, "seedQ[10]/F");
79  tree_->Branch("seedSubdet1", seedSubdet1, "seedSubdet1[10]/I");
80  tree_->Branch("seedLayer1", seedLayer1, "seedLayer1[10]/I");
81  tree_->Branch("seedSide1", seedSide1, "seedSide1[10]/I");
82  tree_->Branch("seedPhi1", seedPhi1, "seedPhi1[10]/F");
83  tree_->Branch("seedDphi1", seedDphi1, "seedDphi1[10]/F");
84  tree_->Branch("seedDrz1", seedDrz1, "seedDrz1[10]/F");
85  tree_->Branch("seedRz1", seedRz1, "seedRz1[10]/F");
86  tree_->Branch("seedSubdet2", seedSubdet2, "seedSubdet2[10]/I");
87  tree_->Branch("seedLayer2", seedLayer2, "seedLayer2[10]/I");
88  tree_->Branch("seedSide2", seedSide2, "seedSide2[10]/I");
89  tree_->Branch("seedPhi2", seedPhi2, "seedPhi2[10]/F");
90  tree_->Branch("seedDphi2", seedDphi2, "seedDphi2[10]/F");
91  tree_->Branch("seedRz2", seedRz2, "seedRz2[10]/F");
92  tree_->Branch("seedDrz2", seedDrz2, "seedDrz2[10]/F");
93  histeMC_ = new TH1F("eMC", "MC particle energy", 100, 0., 100.);
94  histeMCmatched_ = new TH1F("eMCmatched", "matched MC particle energy", 100, 0., 100.);
95  histecaldriveneMCmatched_ =
96  new TH1F("ecaldriveneMCmatched", "matched MC particle energy, ecal driven", 100, 0., 100.);
97  histtrackerdriveneMCmatched_ =
98  new TH1F("trackerdriveneMCmatched", "matched MC particle energy, tracker driven", 100, 0., 100.);
99  histp_ = new TH1F("p", "seed p", 100, 0., 100.);
100  histeclu_ = new TH1F("clus energy", "supercluster energy", 100, 0., 100.);
101  histpt_ = new TH1F("pt", "seed pt", 100, 0., 100.);
102  histptMC_ = new TH1F("ptMC", "MC particle pt", 100, 0., 100.);
103  histptMCmatched_ = new TH1F("ptMCmatched", "matched MC particle pt", 100, 0., 100.);
104  histecaldrivenptMCmatched_ = new TH1F("ecaldrivenptMCmatched", "matched MC particle pt, ecal driven", 100, 0., 100.);
105  histtrackerdrivenptMCmatched_ =
106  new TH1F("trackerdrivenptMCmatched", "matched MC particle pt, tracker driven", 100, 0., 100.);
107  histetclu_ = new TH1F("Et", "supercluster Et", 100, 0., 100.);
108  histeffpt_ = new TH1F("pt eff", "seed effciency vs pt", 100, 0., 100.);
109  histeta_ = new TH1F("seed eta", "seed eta", 100, -2.5, 2.5);
110  histetaMC_ = new TH1F("etaMC", "MC particle eta", 100, -2.5, 2.5);
111  histetaMCmatched_ = new TH1F("etaMCmatched", "matched MC particle eta", 100, -2.5, 2.5);
112  histecaldrivenetaMCmatched_ =
113  new TH1F("ecaldrivenetaMCmatched", "matched MC particle eta, ecal driven", 100, -2.5, 2.5);
114  histtrackerdrivenetaMCmatched_ =
115  new TH1F("trackerdrivenetaMCmatched", "matched MC particle eta, tracker driven", 100, -2.5, 2.5);
116  histetaclu_ = new TH1F("clus eta", "supercluster eta", 100, -2.5, 2.5);
117  histeffeta_ = new TH1F("eta eff", "seed effciency vs eta", 100, -2.5, 2.5);
118  histq_ = new TH1F("q", "seed charge", 100, -2.5, 2.5);
119  histeoverp_ = new TH1F("E/p", "seed E/p", 100, 0., 10.);
120  histnbseeds_ = new TH1I("nrs", "Nr of seeds ", 50, 0., 25.);
121  histnbclus_ = new TH1I("nrclus", "Nr of superclusters ", 50, 0., 25.);
122  histnrseeds_ = new TH1I("ns", "Nr of seeds if clusters", 50, 0., 25.);
123 }

References histecaldriveneMCmatched_, histecaldrivenetaMCmatched_, histecaldrivenptMCmatched_, histeclu_, histeffeta_, histeffpt_, histeMC_, histeMCmatched_, histeoverp_, histeta_, histetaclu_, histetaMC_, histetaMCmatched_, histetclu_, histfile_, histnbclus_, histnbseeds_, histnrseeds_, histp_, histpt_, histptMC_, histptMCmatched_, histq_, histtrackerdriveneMCmatched_, histtrackerdrivenetaMCmatched_, histtrackerdrivenptMCmatched_, mcEnergy, mcEta, mcPhi, mcPt, mcQ, seedDphi1, seedDphi2, seedDrz1, seedDrz2, seedEta, seedLayer1, seedLayer2, seedMomentum, seedPhi, seedPhi1, seedPhi2, seedPt, seedQ, seedRz1, seedRz2, seedSide1, seedSide2, seedSubdet1, seedSubdet2, superclusterEnergy, superclusterEt, superclusterEta, superclusterPhi, and tree_.

◆ endJob()

void ElectronSeedAnalyzer::endJob ( void  )
overridevirtual

Reimplemented from edm::one::EDAnalyzerBase.

Definition at line 125 of file ElectronSeedAnalyzer.cc.

125  {
126  histfile_->cd();
127  tree_->Print();
128  tree_->Write();
129 
130  // efficiency vs eta
131  TH1F *histetaEff = (TH1F *)histetaMCmatched_->Clone("histetaEff");
132  histetaEff->Reset();
133  histetaEff->Divide(histetaMCmatched_, histeta_, 1, 1, "b");
134  histetaEff->Print();
135  histetaEff->GetXaxis()->SetTitle("#eta");
136  histetaEff->GetYaxis()->SetTitle("Efficiency");
137 
138  // efficiency vs pt
139  TH1F *histptEff = (TH1F *)histptMCmatched_->Clone("histotEff");
140  histptEff->Reset();
141  histptEff->Divide(histptMCmatched_, histpt_, 1, 1, "b");
142  histptEff->Print();
143  histptEff->GetXaxis()->SetTitle("p_{T}");
144  histptEff->GetYaxis()->SetTitle("Efficiency");
145 
146  histeMCmatched_->Write();
147  histecaldriveneMCmatched_->Write();
148  histtrackerdriveneMCmatched_->Write();
149  histeMC_->Write();
150  histp_->Write();
151  histeclu_->Write();
152  histpt_->Write();
153  histptMCmatched_->Write();
154  histecaldrivenptMCmatched_->Write();
155  histtrackerdrivenptMCmatched_->Write();
156  histptMC_->Write();
157  histetclu_->Write();
158  histeffpt_->Write();
159  histeta_->Write();
160  histetaMCmatched_->Write();
161  histecaldrivenetaMCmatched_->Write();
162  histtrackerdrivenetaMCmatched_->Write();
163  histetaMC_->Write();
164  histetaclu_->Write();
165  histeffeta_->Write();
166  histq_->Write();
167  histeoverp_->Write();
168  histnbseeds_->Write();
169  histnbclus_->Write();
170  histnrseeds_->Write();
171 }

References histecaldriveneMCmatched_, histecaldrivenetaMCmatched_, histecaldrivenptMCmatched_, histeclu_, histeffeta_, histeffpt_, histeMC_, histeMCmatched_, histeoverp_, histeta_, histetaclu_, histetaMC_, histetaMCmatched_, histetclu_, histfile_, histnbclus_, histnbseeds_, histnrseeds_, histp_, histpt_, histptMC_, histptMCmatched_, histq_, histtrackerdriveneMCmatched_, histtrackerdrivenetaMCmatched_, histtrackerdrivenptMCmatched_, and tree_.

Member Data Documentation

◆ beamSpot_

edm::InputTag ElectronSeedAnalyzer::beamSpot_
private

Definition at line 81 of file ElectronSeedAnalyzer.h.

Referenced by analyze().

◆ histecaldriveneMCmatched_

TH1F* ElectronSeedAnalyzer::histecaldriveneMCmatched_
private

Definition at line 56 of file ElectronSeedAnalyzer.h.

Referenced by analyze(), beginJob(), and endJob().

◆ histecaldrivenetaMCmatched_

TH1F* ElectronSeedAnalyzer::histecaldrivenetaMCmatched_
private

Definition at line 70 of file ElectronSeedAnalyzer.h.

Referenced by analyze(), beginJob(), and endJob().

◆ histecaldrivenptMCmatched_

TH1F* ElectronSeedAnalyzer::histecaldrivenptMCmatched_
private

Definition at line 63 of file ElectronSeedAnalyzer.h.

Referenced by analyze(), beginJob(), and endJob().

◆ histeclu_

TH1F* ElectronSeedAnalyzer::histeclu_
private

Definition at line 59 of file ElectronSeedAnalyzer.h.

Referenced by beginJob(), and endJob().

◆ histeffeta_

TH1F* ElectronSeedAnalyzer::histeffeta_
private

Definition at line 73 of file ElectronSeedAnalyzer.h.

Referenced by beginJob(), and endJob().

◆ histeffpt_

TH1F* ElectronSeedAnalyzer::histeffpt_
private

Definition at line 66 of file ElectronSeedAnalyzer.h.

Referenced by beginJob(), and endJob().

◆ histeMC_

TH1F* ElectronSeedAnalyzer::histeMC_
private

Definition at line 54 of file ElectronSeedAnalyzer.h.

Referenced by analyze(), beginJob(), endJob(), and ~ElectronSeedAnalyzer().

◆ histeMCmatched_

TH1F* ElectronSeedAnalyzer::histeMCmatched_
private

Definition at line 55 of file ElectronSeedAnalyzer.h.

Referenced by analyze(), beginJob(), and endJob().

◆ histeoverp_

TH1F* ElectronSeedAnalyzer::histeoverp_
private

Definition at line 75 of file ElectronSeedAnalyzer.h.

Referenced by analyze(), beginJob(), and endJob().

◆ histeta_

TH1F* ElectronSeedAnalyzer::histeta_
private

Definition at line 67 of file ElectronSeedAnalyzer.h.

Referenced by analyze(), beginJob(), and endJob().

◆ histetaclu_

TH1F* ElectronSeedAnalyzer::histetaclu_
private

Definition at line 72 of file ElectronSeedAnalyzer.h.

Referenced by analyze(), beginJob(), and endJob().

◆ histetaMC_

TH1F* ElectronSeedAnalyzer::histetaMC_
private

Definition at line 68 of file ElectronSeedAnalyzer.h.

Referenced by analyze(), beginJob(), and endJob().

◆ histetaMCmatched_

TH1F* ElectronSeedAnalyzer::histetaMCmatched_
private

Definition at line 69 of file ElectronSeedAnalyzer.h.

Referenced by analyze(), beginJob(), and endJob().

◆ histetclu_

TH1F* ElectronSeedAnalyzer::histetclu_
private

Definition at line 65 of file ElectronSeedAnalyzer.h.

Referenced by analyze(), beginJob(), and endJob().

◆ histfile_

TFile* ElectronSeedAnalyzer::histfile_
private

Definition at line 44 of file ElectronSeedAnalyzer.h.

Referenced by beginJob(), ElectronSeedAnalyzer(), endJob(), and ~ElectronSeedAnalyzer().

◆ histnbclus_

TH1I* ElectronSeedAnalyzer::histnbclus_
private

Definition at line 78 of file ElectronSeedAnalyzer.h.

Referenced by analyze(), beginJob(), and endJob().

◆ histnbseeds_

TH1I* ElectronSeedAnalyzer::histnbseeds_
private

Definition at line 77 of file ElectronSeedAnalyzer.h.

Referenced by analyze(), beginJob(), and endJob().

◆ histnrseeds_

TH1I* ElectronSeedAnalyzer::histnrseeds_
private

Definition at line 76 of file ElectronSeedAnalyzer.h.

Referenced by analyze(), beginJob(), and endJob().

◆ histp_

TH1F* ElectronSeedAnalyzer::histp_
private

Definition at line 58 of file ElectronSeedAnalyzer.h.

Referenced by beginJob(), and endJob().

◆ histpt_

TH1F* ElectronSeedAnalyzer::histpt_
private

Definition at line 60 of file ElectronSeedAnalyzer.h.

Referenced by analyze(), beginJob(), and endJob().

◆ histptMC_

TH1F* ElectronSeedAnalyzer::histptMC_
private

Definition at line 61 of file ElectronSeedAnalyzer.h.

Referenced by analyze(), beginJob(), and endJob().

◆ histptMCmatched_

TH1F* ElectronSeedAnalyzer::histptMCmatched_
private

Definition at line 62 of file ElectronSeedAnalyzer.h.

Referenced by analyze(), beginJob(), and endJob().

◆ histq_

TH1F* ElectronSeedAnalyzer::histq_
private

Definition at line 74 of file ElectronSeedAnalyzer.h.

Referenced by analyze(), beginJob(), and endJob().

◆ histtrackerdriveneMCmatched_

TH1F* ElectronSeedAnalyzer::histtrackerdriveneMCmatched_
private

Definition at line 57 of file ElectronSeedAnalyzer.h.

Referenced by analyze(), beginJob(), and endJob().

◆ histtrackerdrivenetaMCmatched_

TH1F* ElectronSeedAnalyzer::histtrackerdrivenetaMCmatched_
private

Definition at line 71 of file ElectronSeedAnalyzer.h.

Referenced by analyze(), beginJob(), and endJob().

◆ histtrackerdrivenptMCmatched_

TH1F* ElectronSeedAnalyzer::histtrackerdrivenptMCmatched_
private

Definition at line 64 of file ElectronSeedAnalyzer.h.

Referenced by analyze(), beginJob(), and endJob().

◆ inputCollection_

edm::InputTag ElectronSeedAnalyzer::inputCollection_
private

Definition at line 80 of file ElectronSeedAnalyzer.h.

Referenced by analyze(), and ElectronSeedAnalyzer().

◆ mcEnergy

float ElectronSeedAnalyzer::mcEnergy[10]
private

Definition at line 46 of file ElectronSeedAnalyzer.h.

Referenced by analyze(), and beginJob().

◆ mcEta

float ElectronSeedAnalyzer::mcEta[10]
private

Definition at line 46 of file ElectronSeedAnalyzer.h.

Referenced by analyze(), and beginJob().

◆ mcPhi

float ElectronSeedAnalyzer::mcPhi[10]
private

Definition at line 46 of file ElectronSeedAnalyzer.h.

Referenced by analyze(), and beginJob().

◆ mcPt

float ElectronSeedAnalyzer::mcPt[10]
private

Definition at line 46 of file ElectronSeedAnalyzer.h.

Referenced by analyze(), and beginJob().

◆ mcQ

float ElectronSeedAnalyzer::mcQ[10]
private

Definition at line 46 of file ElectronSeedAnalyzer.h.

Referenced by analyze(), and beginJob().

◆ seedDphi1

float ElectronSeedAnalyzer::seedDphi1[10]
private

Definition at line 52 of file ElectronSeedAnalyzer.h.

Referenced by analyze(), and beginJob().

◆ seedDphi2

float ElectronSeedAnalyzer::seedDphi2[10]
private

Definition at line 52 of file ElectronSeedAnalyzer.h.

Referenced by analyze(), and beginJob().

◆ seedDrz1

float ElectronSeedAnalyzer::seedDrz1[10]
private

Definition at line 52 of file ElectronSeedAnalyzer.h.

Referenced by analyze(), and beginJob().

◆ seedDrz2

float ElectronSeedAnalyzer::seedDrz2[10]
private

Definition at line 52 of file ElectronSeedAnalyzer.h.

Referenced by analyze(), and beginJob().

◆ seedEta

float ElectronSeedAnalyzer::seedEta[10]
private

Definition at line 48 of file ElectronSeedAnalyzer.h.

Referenced by analyze(), and beginJob().

◆ seedLayer1

int ElectronSeedAnalyzer::seedLayer1[10]
private

Definition at line 50 of file ElectronSeedAnalyzer.h.

Referenced by analyze(), and beginJob().

◆ seedLayer2

int ElectronSeedAnalyzer::seedLayer2[10]
private

Definition at line 50 of file ElectronSeedAnalyzer.h.

Referenced by analyze(), and beginJob().

◆ seedMomentum

float ElectronSeedAnalyzer::seedMomentum[10]
private

Definition at line 48 of file ElectronSeedAnalyzer.h.

Referenced by analyze(), and beginJob().

◆ seedPhi

float ElectronSeedAnalyzer::seedPhi[10]
private

Definition at line 48 of file ElectronSeedAnalyzer.h.

Referenced by analyze(), and beginJob().

◆ seedPhi1

float ElectronSeedAnalyzer::seedPhi1[10]
private

Definition at line 53 of file ElectronSeedAnalyzer.h.

Referenced by analyze(), and beginJob().

◆ seedPhi2

float ElectronSeedAnalyzer::seedPhi2[10]
private

Definition at line 53 of file ElectronSeedAnalyzer.h.

Referenced by analyze(), and beginJob().

◆ seedPt

float ElectronSeedAnalyzer::seedPt[10]
private

Definition at line 48 of file ElectronSeedAnalyzer.h.

Referenced by analyze(), and beginJob().

◆ seedQ

float ElectronSeedAnalyzer::seedQ[10]
private

Definition at line 48 of file ElectronSeedAnalyzer.h.

Referenced by analyze(), and beginJob().

◆ seedRz1

float ElectronSeedAnalyzer::seedRz1[10]
private

Definition at line 53 of file ElectronSeedAnalyzer.h.

Referenced by analyze(), and beginJob().

◆ seedRz2

float ElectronSeedAnalyzer::seedRz2[10]
private

Definition at line 53 of file ElectronSeedAnalyzer.h.

Referenced by analyze(), and beginJob().

◆ seedSide1

int ElectronSeedAnalyzer::seedSide1[10]
private

Definition at line 51 of file ElectronSeedAnalyzer.h.

Referenced by analyze(), and beginJob().

◆ seedSide2

int ElectronSeedAnalyzer::seedSide2[10]
private

Definition at line 51 of file ElectronSeedAnalyzer.h.

Referenced by analyze(), and beginJob().

◆ seedSubdet1

int ElectronSeedAnalyzer::seedSubdet1[10]
private

Definition at line 49 of file ElectronSeedAnalyzer.h.

Referenced by analyze(), and beginJob().

◆ seedSubdet2

int ElectronSeedAnalyzer::seedSubdet2[10]
private

Definition at line 49 of file ElectronSeedAnalyzer.h.

Referenced by analyze(), and beginJob().

◆ superclusterEnergy

float ElectronSeedAnalyzer::superclusterEnergy[10]
private

Definition at line 47 of file ElectronSeedAnalyzer.h.

Referenced by analyze(), and beginJob().

◆ superclusterEt

float ElectronSeedAnalyzer::superclusterEt[10]
private

Definition at line 47 of file ElectronSeedAnalyzer.h.

Referenced by analyze(), and beginJob().

◆ superclusterEta

float ElectronSeedAnalyzer::superclusterEta[10]
private

Definition at line 47 of file ElectronSeedAnalyzer.h.

Referenced by analyze(), and beginJob().

◆ superclusterPhi

float ElectronSeedAnalyzer::superclusterPhi[10]
private

Definition at line 47 of file ElectronSeedAnalyzer.h.

Referenced by analyze(), and beginJob().

◆ transformer_

TrajectoryStateTransform ElectronSeedAnalyzer::transformer_
private

Definition at line 42 of file ElectronSeedAnalyzer.h.

◆ tree_

TTree* ElectronSeedAnalyzer::tree_
private

Definition at line 45 of file ElectronSeedAnalyzer.h.

Referenced by beginJob(), and endJob().

TrackerGeometry::idToDet
const TrackerGeomDet * idToDet(DetId) const override
Definition: TrackerGeometry.cc:193
ElectronSeedAnalyzer::histecaldrivenetaMCmatched_
TH1F * histecaldrivenetaMCmatched_
Definition: ElectronSeedAnalyzer.h:70
edm::RefToBase::isNull
bool isNull() const
Checks for null.
Definition: RefToBase.h:295
ElectronSeedAnalyzer::superclusterEt
float superclusterEt[10]
Definition: ElectronSeedAnalyzer.h:47
FastTimerService_cff.range
range
Definition: FastTimerService_cff.py:34
TrackerTopology::side
unsigned int side(const DetId &id) const
Definition: TrackerTopology.cc:28
ElectronSeedAnalyzer::histtrackerdriveneMCmatched_
TH1F * histtrackerdriveneMCmatched_
Definition: ElectronSeedAnalyzer.h:57
ElectronSeedAnalyzer::histnbseeds_
TH1I * histnbseeds_
Definition: ElectronSeedAnalyzer.h:77
GeomDet
Definition: GeomDet.h:27
ElectronSeedAnalyzer::histecaldriveneMCmatched_
TH1F * histecaldriveneMCmatched_
Definition: ElectronSeedAnalyzer.h:56
edm::Handle::product
T const * product() const
Definition: Handle.h:70
PropagatorWithMaterial::propagate
TrajectoryStateOnSurface propagate(STA const &state, SUR const &surface) const
Definition: Propagator.h:50
ElectronSeedAnalyzer::histfile_
TFile * histfile_
Definition: ElectronSeedAnalyzer.h:44
ElectronSeedAnalyzer::seedDrz1
float seedDrz1[10]
Definition: ElectronSeedAnalyzer.h:52
PV3DBase::x
T x() const
Definition: PV3DBase.h:59
AlCaHLTBitMon_ParallelJobs.p
p
Definition: AlCaHLTBitMon_ParallelJobs.py:153
globals_cff.id1
id1
Definition: globals_cff.py:32
TrajectoryStateOnSurface::globalPosition
GlobalPoint globalPosition() const
Definition: TrajectoryStateOnSurface.h:65
edm::LogInfo
Definition: MessageLogger.h:254
ElectronSeedAnalyzer::histpt_
TH1F * histpt_
Definition: ElectronSeedAnalyzer.h:60
ElectronSeedAnalyzer::histeta_
TH1F * histeta_
Definition: ElectronSeedAnalyzer.h:67
ElectronSeedAnalyzer::superclusterEta
float superclusterEta[10]
Definition: ElectronSeedAnalyzer.h:47
ElectronSeedAnalyzer::seedMomentum
float seedMomentum[10]
Definition: ElectronSeedAnalyzer.h:48
oppositeToMomentum
Definition: PropagationDirection.h:4
TrajectorySeed::const_iterator
recHitContainer::const_iterator const_iterator
Definition: TrajectorySeed.h:20
TrackerTopology::layer
unsigned int layer(const DetId &id) const
Definition: TrackerTopology.cc:47
ElectronSeedAnalyzer::histtrackerdrivenetaMCmatched_
TH1F * histtrackerdrivenetaMCmatched_
Definition: ElectronSeedAnalyzer.h:71
ElectronSeedAnalyzer::mcEta
float mcEta[10]
Definition: ElectronSeedAnalyzer.h:46
ElectronSeedAnalyzer::histeMCmatched_
TH1F * histeMCmatched_
Definition: ElectronSeedAnalyzer.h:55
ElectronSeedAnalyzer::seedSide2
int seedSide2[10]
Definition: ElectronSeedAnalyzer.h:51
ElectronSeedAnalyzer::histeffpt_
TH1F * histeffpt_
Definition: ElectronSeedAnalyzer.h:66
edm::Handle< reco::BeamSpot >
ElectronSeedAnalyzer::seedSubdet1
int seedSubdet1[10]
Definition: ElectronSeedAnalyzer.h:49
ElectronSeedAnalyzer::seedPhi2
float seedPhi2[10]
Definition: ElectronSeedAnalyzer.h:53
ElectronSeedAnalyzer::histeffeta_
TH1F * histeffeta_
Definition: ElectronSeedAnalyzer.h:73
HepMC::GenEvent
Definition: hepmc_rootio.cc:9
ElectronSeedAnalyzer::histtrackerdrivenptMCmatched_
TH1F * histtrackerdrivenptMCmatched_
Definition: ElectronSeedAnalyzer.h:64
edm::Ref< SuperClusterCollection >
funct::sin
Sin< T >::type sin(const T &t)
Definition: Sin.h:22
PV3DBase::z
T z() const
Definition: PV3DBase.h:61
IdealMagneticFieldRecord
Definition: IdealMagneticFieldRecord.h:11
DetId
Definition: DetId.h:17
ElectronSeedAnalyzer::histnrseeds_
TH1I * histnrseeds_
Definition: ElectronSeedAnalyzer.h:76
ElectronSeedAnalyzer::histnbclus_
TH1I * histnbclus_
Definition: ElectronSeedAnalyzer.h:78
GeomDet::surface
const Plane & surface() const
The nominal surface of the GeomDet.
Definition: GeomDet.h:37
TrajectoryStateOnSurface
Definition: TrajectoryStateOnSurface.h:16
PropagatorWithMaterial
Definition: PropagatorWithMaterial.h:25
reco::ElectronSeed
Definition: ElectronSeed.h:51
trackingTools::ftsFromVertexToPoint
FreeTrajectoryState ftsFromVertexToPoint(MagneticField const &magField, GlobalPoint const &xmeas, GlobalPoint const &xvert, float momentum, TrackCharge charge)
Definition: ftsFromVertexToPoint.cc:17
ElectronSeedAnalyzer::seedDphi1
float seedDphi1[10]
Definition: ElectronSeedAnalyzer.h:52
edm::EventSetup::get
T get() const
Definition: EventSetup.h:73
ElectronSeedAnalyzer::seedLayer2
int seedLayer2[10]
Definition: ElectronSeedAnalyzer.h:50
PVValHelper::eta
Definition: PVValidationHelpers.h:69
TrackerDigiGeometryRecord
Definition: TrackerDigiGeometryRecord.h:15
mathSSE::sqrt
T sqrt(T t)
Definition: SSEVec.h:19
Surface::toGlobal
GlobalPoint toGlobal(const Point2DBase< Scalar, LocalTag > lp) const
Definition: Surface.h:79
edm::ESHandle< TrackerTopology >
HCALHighEnergyHPDFilter_cfi.energy
energy
Definition: HCALHighEnergyHPDFilter_cfi.py:5
ElectronSeedAnalyzer::mcPhi
float mcPhi[10]
Definition: ElectronSeedAnalyzer.h:46
ElectronSeedAnalyzer::seedPhi1
float seedPhi1[10]
Definition: ElectronSeedAnalyzer.h:53
reco::BeamSpot::position
const Point & position() const
position
Definition: BeamSpot.h:59
Point3DBase< float, GlobalTag >
ElectronSeedAnalyzer::seedPhi
float seedPhi[10]
Definition: ElectronSeedAnalyzer.h:48
OrderedSet.t
t
Definition: OrderedSet.py:90
PbPb_ZMuSkimMuonDPG_cff.deltaR
deltaR
Definition: PbPb_ZMuSkimMuonDPG_cff.py:63
ElectronSeedAnalyzer::histetaMC_
TH1F * histetaMC_
Definition: ElectronSeedAnalyzer.h:68
ElectronSeedAnalyzer::histeclu_
TH1F * histeclu_
Definition: ElectronSeedAnalyzer.h:59
bsc_activity_cfg.clusters
clusters
Definition: bsc_activity_cfg.py:36
ALCARECOTkAlJpsiMuMu_cff.charge
charge
Definition: ALCARECOTkAlJpsiMuMu_cff.py:47
GeomDet::toGlobal
GlobalPoint toGlobal(const Local2DPoint &lp) const
Conversion to the global R.F. from the R.F. of the GeomDet.
Definition: GeomDet.h:49
ElectronSeedAnalyzer::seedDrz2
float seedDrz2[10]
Definition: ElectronSeedAnalyzer.h:52
LogDebug
#define LogDebug(id)
Definition: MessageLogger.h:670
ElectronSeedAnalyzer::histeMC_
TH1F * histeMC_
Definition: ElectronSeedAnalyzer.h:54
ElectronSeedAnalyzer::tree_
TTree * tree_
Definition: ElectronSeedAnalyzer.h:45
ElectronSeedAnalyzer::histeoverp_
TH1F * histeoverp_
Definition: ElectronSeedAnalyzer.h:75
ElectronSeedAnalyzer::histecaldrivenptMCmatched_
TH1F * histecaldrivenptMCmatched_
Definition: ElectronSeedAnalyzer.h:63
ElectronSeedAnalyzer::mcEnergy
float mcEnergy[10]
Definition: ElectronSeedAnalyzer.h:46
ElectronSeedAnalyzer::histetaclu_
TH1F * histetaclu_
Definition: ElectronSeedAnalyzer.h:72
PV3DBase::y
T y() const
Definition: PV3DBase.h:60
createfilelist.int
int
Definition: createfilelist.py:10
trajectoryStateTransform::transientState
TrajectoryStateOnSurface transientState(const PTrajectoryStateOnDet &ts, const Surface *surface, const MagneticField *field)
Definition: TrajectoryStateTransform.cc:35
ElectronSeedAnalyzer::beamSpot_
edm::InputTag beamSpot_
Definition: ElectronSeedAnalyzer.h:81
ElectronSeedAnalyzer::histetaMCmatched_
TH1F * histetaMCmatched_
Definition: ElectronSeedAnalyzer.h:69
ElectronSeedAnalyzer::histetclu_
TH1F * histetclu_
Definition: ElectronSeedAnalyzer.h:65
edm::RefToBase::castTo
REF castTo() const
Definition: RefToBase.h:257
edm::HepMCProduct::GetEvent
const HepMC::GenEvent * GetEvent() const
Definition: HepMCProduct.h:34
get
#define get
ElectronSeedAnalyzer::histq_
TH1F * histq_
Definition: ElectronSeedAnalyzer.h:74
alignCSCRings.r
r
Definition: alignCSCRings.py:93
ElectronSeedAnalyzer::seedEta
float seedEta[10]
Definition: ElectronSeedAnalyzer.h:48
ElectronSeedAnalyzer::seedRz2
float seedRz2[10]
Definition: ElectronSeedAnalyzer.h:53
DDAxes::phi
edm::ParameterSet::getParameter
T getParameter(std::string const &) const
ElectronSeedAnalyzer::seedRz1
float seedRz1[10]
Definition: ElectronSeedAnalyzer.h:53
ElectronSeedAnalyzer::mcQ
float mcQ[10]
Definition: ElectronSeedAnalyzer.h:46
GenParticle.GenParticle
GenParticle
Definition: GenParticle.py:18
ElectronSeedAnalyzer::seedDphi2
float seedDphi2[10]
Definition: ElectronSeedAnalyzer.h:52
ElectronSeedAnalyzer::seedLayer1
int seedLayer1[10]
Definition: ElectronSeedAnalyzer.h:50
ElectronSeedAnalyzer::seedPt
float seedPt[10]
Definition: ElectronSeedAnalyzer.h:48
ElectronSeedAnalyzer::seedSide1
int seedSide1[10]
Definition: ElectronSeedAnalyzer.h:51
PerpendicularBoundPlaneBuilder
Definition: PerpendicularBoundPlaneBuilder.h:11
ElectronSeedAnalyzer::seedSubdet2
int seedSubdet2[10]
Definition: ElectronSeedAnalyzer.h:49
ElectronSeedAnalyzer::superclusterPhi
float superclusterPhi[10]
Definition: ElectronSeedAnalyzer.h:47
EgHLTOffHistBins_cfi.mass
mass
Definition: EgHLTOffHistBins_cfi.py:34
edm::RefToBase< CaloCluster >
ElectronSeedAnalyzer::histptMC_
TH1F * histptMC_
Definition: ElectronSeedAnalyzer.h:61
funct::pow
Power< A, B >::type pow(const A &a, const B &b)
Definition: Power.h:30
ElectronSeedAnalyzer::mcPt
float mcPt[10]
Definition: ElectronSeedAnalyzer.h:46
pi
const Double_t pi
Definition: trackSplitPlot.h:36
funct::abs
Abs< T >::type abs(const T &t)
Definition: Abs.h:22
ElectronSeedAnalyzer::superclusterEnergy
float superclusterEnergy[10]
Definition: ElectronSeedAnalyzer.h:47
ElectronSeedAnalyzer::histptMCmatched_
TH1F * histptMCmatched_
Definition: ElectronSeedAnalyzer.h:62
TrackerTopologyRcd
Definition: TrackerTopologyRcd.h:10
globals_cff.id2
id2
Definition: globals_cff.py:33
ElectronSeedAnalyzer::histp_
TH1F * histp_
Definition: ElectronSeedAnalyzer.h:58
JetChargeProducer_cfi.exp
exp
Definition: JetChargeProducer_cfi.py:6
parents
TPRegexp parents
Definition: eve_filter.cc:21
ElectronSeedAnalyzer::inputCollection_
edm::InputTag inputCollection_
Definition: ElectronSeedAnalyzer.h:80
PV3DBase::perp
T perp() const
Definition: PV3DBase.h:69
ElectronSeedAnalyzer::seedQ
float seedQ[10]
Definition: ElectronSeedAnalyzer.h:48
edm::InputTag
Definition: InputTag.h:15
alongMomentum
Definition: PropagationDirection.h:4
PV3DBase::phi
Geom::Phi< T > phi() const
Definition: PV3DBase.h:66
TrajectoryStateOnSurface::isValid
bool isValid() const
Definition: TrajectoryStateOnSurface.h:54
edm::OwnVector< TrackingRecHit >
MillePedeFileConverter_cfg.e
e
Definition: MillePedeFileConverter_cfg.py:37
Generated for CMSSW Reference Manual by doxygen 1.8.16