test
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 Types | Private Attributes
SiStripLAProfileBooker Class Reference

#include <SiStripLAProfileBooker.h>

Inheritance diagram for SiStripLAProfileBooker:
edm::EDAnalyzer edm::EDConsumerBase

Classes

struct  detparameters
 

Public Member Functions

void analyze (const edm::Event &e, const edm::EventSetup &c)
 
void beginRun (const edm::EventSetup &c)
 
void endJob ()
 
void getlayer (const DetId &detid, const TrackerTopology *tTopo, std::string &name, unsigned int &layerid)
 
 SiStripLAProfileBooker (const edm::ParameterSet &conf)
 
 ~SiStripLAProfileBooker ()
 
- Public Member Functions inherited from edm::EDAnalyzer
void callWhenNewProductsRegistered (std::function< void(BranchDescription const &)> const &func)
 
 EDAnalyzer ()
 
ModuleDescription const & moduleDescription () const
 
std::string workerType () const
 
virtual ~EDAnalyzer ()
 
- Public Member Functions inherited from edm::EDConsumerBase
 EDConsumerBase ()
 
ProductHolderIndexAndSkipBit indexFrom (EDGetToken, BranchType, TypeID const &) const
 
void itemsMayGet (BranchType, std::vector< ProductHolderIndexAndSkipBit > &) const
 
void itemsToGet (BranchType, std::vector< ProductHolderIndexAndSkipBit > &) const
 
std::vector
< ProductHolderIndexAndSkipBit >
const & 
itemsToGetFromEvent () const
 
void labelsForToken (EDGetToken iToken, Labels &oLabels) const
 
void modulesDependentUpon (const std::string &iProcessName, std::vector< const char * > &oModuleLabels) const
 
bool registeredToConsume (ProductHolderIndex, bool, BranchType) const
 
bool registeredToConsumeMany (TypeID const &, BranchType) const
 
void updateLookup (BranchType iBranchType, ProductHolderIndexHelper const &)
 
virtual ~EDConsumerBase ()
 

Private Types

typedef std::map< unsigned int,
detparameters * > 
detparmap
 
typedef std::map< unsigned int,
MonitorElement * > 
histomap
 

Private Attributes

uint8_t Amplitudes [100]
 
float barycenter
 
int better_double_hit
 
int bw_fw
 
float chi2norm
 
int ClSize
 
edm::ParameterSet conf_
 
DQMStoredbe_
 
detparmap detmap
 
float EtaTrack
 
TDirectory * Event_Tree
 
int eventcounter
 
int EventCounter
 
int EventNumber
 
TTree * EventTree
 
int Ext_Int
 
TFile * hFile
 
histomap histos
 
float hit_std_dev
 
TDirectory * Hit_Tree
 
int HitCharge
 
int hitcounter
 
int hitcounter_2ndloop
 
int HitNr
 
int HitPerTrack
 
TTree * HitsTree
 
int id_detector
 
int Layer
 
float MagField
 
float Momentum
 
int MonoStereo
 
int nstrip
 
int ParticleCharge
 
float PhiTrack
 
float pitch_detector
 
float pt
 
int RunNumber
 
float SignCorrection
 
int size
 
detparmap summarydetmap
 
histomap summaryhisto
 
float sumx
 
float TanTrackAngle
 
float TanTrackAngleParallel
 
float thick_detector
 
TDirectory * Track_Tree
 
int trackcollsize
 
int TrackCounter
 
const TrackerGeometrytracker
 
TTree * TrackTree
 
int trajsize
 
std::string treename_
 
int Type
 
int Wheel
 
int worse_double_hit
 
float XGlobal
 
float YGlobal
 
float ZGlobal
 

Additional Inherited Members

- Public Types inherited from edm::EDAnalyzer
typedef EDAnalyzer ModuleType
 
- Static Public Member Functions inherited from edm::EDAnalyzer
static const std::string & baseType ()
 
static void fillDescriptions (ConfigurationDescriptions &descriptions)
 
static void prevalidate (ConfigurationDescriptions &)
 
- Protected Member Functions inherited from edm::EDConsumerBase
template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType > consumes (edm::InputTag const &tag)
 
EDGetToken consumes (const TypeToGet &id, 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 ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType > mayConsume (edm::InputTag const &tag)
 
EDGetToken mayConsume (const TypeToGet &id, edm::InputTag const &tag)
 
template<BranchType B>
EDGetToken mayConsume (const TypeToGet &id, edm::InputTag const &tag)
 

Detailed Description

Definition at line 30 of file SiStripLAProfileBooker.h.

Member Typedef Documentation

typedef std::map<unsigned int, detparameters*> SiStripLAProfileBooker::detparmap
private

Definition at line 49 of file SiStripLAProfileBooker.h.

typedef std::map<unsigned int, MonitorElement*> SiStripLAProfileBooker::histomap
private

Definition at line 50 of file SiStripLAProfileBooker.h.

Constructor & Destructor Documentation

SiStripLAProfileBooker::SiStripLAProfileBooker ( const edm::ParameterSet conf)
explicit

Definition at line 64 of file SiStripLAProfileBooker.cc.

64  :
65  conf_(conf)
66 {
67 }
SiStripLAProfileBooker::~SiStripLAProfileBooker ( )

Definition at line 261 of file SiStripLAProfileBooker.cc.

References detmap, hFile, and summarydetmap.

261  {
262  detparmap::iterator detpariter;
263  for( detpariter=detmap.begin(); detpariter!=detmap.end();++detpariter)delete detpariter->second;
264  for( detpariter=summarydetmap.begin(); detpariter!=summarydetmap.end();++detpariter)delete detpariter->second;
265  delete hFile;
266 }

Member Function Documentation

void SiStripLAProfileBooker::analyze ( const edm::Event e,
const edm::EventSetup c 
)
virtual

Implements edm::EDAnalyzer.

Definition at line 270 of file SiStripLAProfileBooker.cc.

References Amplitudes, barycenter, begin, better_double_hit, Surface::bounds(), bw_fw, chi2norm, ClSize, SiStripRecHit2D::cluster(), conf_, edm::AssociationMap< edm::OneToOne< std::vector< Trajectory >, reco::TrackCollection, unsigned short > >::const_iterator, cond::rpcobgas::detid, detmap, EtaTrack, edm::EventID::event(), eventcounter, EventNumber, EventTree, Ext_Int, TrackingRecHit::geographicalId(), edm::EventSetup::get(), edm::Event::getByLabel(), getlayer(), edm::ParameterSet::getParameter(), histos, hit_std_dev, HitCharge, hitcounter, hitcounter_2ndloop, HitNr, HitPerTrack, HitsTree, i, edm::EventBase::id(), id_detector, TrackerGeometry::idToDet(), Layer, TrajectoryStateOnSurface::localDirection(), StripTopology::localPitch(), BaseTrackerRecHit::localPosition(), LogDebug, PV3DBase< T, PVType, FrameType >::mag(), MagField, Momentum, GluedGeomDet::monoDet(), SiStripMatchedRecHit2D::monoHit(), MonoStereo, mergeVDriftHistosByStation::name, nstrip, ProjectedSiStripRecHit2D::originalHit(), ParticleCharge, PhiTrack, pitch_detector, position, funct::pow(), edm::ESHandle< class >::product(), pt, DetId::rawId(), edm::EventID::run(), RunNumber, SignCorrection, size, GeomDet::specificSurface(), mathSSE::sqrt(), StripSubdetector::stereo(), GluedGeomDet::stereoDet(), SiStripMatchedRecHit2D::stereoHit(), AlCaHLTBitMon_QueryRunRegistry::string, DetId::subdetId(), summaryhisto, sumx, GeomDet::surface(), TanTrackAngle, TanTrackAngleParallel, StripSubdetector::TEC, TrackerTopology::tecPetalInfo(), TrackerTopology::tecWheel(), thick_detector, Bounds::thickness(), StripSubdetector::TIB, TrackerTopology::tibLayer(), TrackerTopology::tibStringInfo(), StripSubdetector::TID, TrackerTopology::tidModuleInfo(), TrackerTopology::tidWheel(), StripSubdetector::TOB, TrackerTopology::tobLayer(), TrackerTopology::tobRodInfo(), GeomDet::toGlobal(), GeomDet::toLocal(), GeomDetUnit::topology(), trackcollsize, tracker, testEve_cfg::tracks, TrackTree, trajsize, Type, Wheel, worse_double_hit, PV3DBase< T, PVType, FrameType >::x(), XGlobal, PV3DBase< T, PVType, FrameType >::y(), YGlobal, PV3DBase< T, PVType, FrameType >::z(), and ZGlobal.

271 {
272  //Retrieve tracker topology from geometry
273  edm::ESHandle<TrackerTopology> tTopoHandle;
274  es.get<IdealGeometryRecord>().get(tTopoHandle);
275  const TrackerTopology* const tTopo = tTopoHandle.product();
276 
277  RunNumber = e.id().run();
278  EventNumber = e.id().event();
279 
280  eventcounter++;
281 
282  EventTree->Fill();
283 
284  //Analysis of Trajectory-RecHits
285 
286  edm::InputTag TkTag = conf_.getParameter<edm::InputTag>("Tracks");
287 
288  edm::Handle<reco::TrackCollection> trackCollection;
289  e.getByLabel(TkTag,trackCollection);
290 
292  e.getByLabel(TkTag,TrajectoryCollection);
293 
295  e.getByLabel(TkTag, TrajTrackMap);
296 
297  const reco::TrackCollection *tracks=trackCollection.product();
298 
299  // FIXME this has to be changed to use pointers to clusters...
300  std::map<const SiStripRecHit2D*,std::pair<float,float>,DetIdLess> hitangleassociation;
301  std::list<SiStripRecHit2D> cache; // ugly, inefficient, effective in making the above working
302 
303  trackcollsize=tracks->size();
304  trajsize=TrajectoryCollection->size();
305 
306  edm::LogInfo("SiStripLAProfileBooker::analyze") <<" Number of tracks in event = "<<trackcollsize<<"\n";
307  edm::LogInfo("SiStripLAProfileBooker::analyze") <<" Number of trajectories in event = "<<trajsize<<"\n";
308 
310 
311  for(TrajTrackIter = TrajTrackMap->begin(); TrajTrackIter!= TrajTrackMap->end(); TrajTrackIter++){ //loop on trajectories
312 
313  if(TrajTrackIter->key->foundHits()>=5){
314 
315  TrackTree->Fill();
316 
317  ParticleCharge = -99;
318  Momentum = -99;
319  pt = -99;
320  chi2norm = -99;
321  HitPerTrack = -99;
322  EtaTrack = -99;
323  PhiTrack = -99;
324 
325  ParticleCharge = TrajTrackIter->val->charge();
326  pt = TrajTrackIter->val->pt();
327  Momentum = TrajTrackIter->val->p();
328  chi2norm = TrajTrackIter->val->normalizedChi2();
329  EtaTrack = TrajTrackIter->val->eta();
330  PhiTrack = TrajTrackIter->val->phi();
331  HitPerTrack = TrajTrackIter->key->foundHits();
332 
333  std::vector<TrajectoryMeasurement> TMeas=TrajTrackIter->key->measurements();
334  std::vector<TrajectoryMeasurement>::iterator itm;
335 
336  for (itm=TMeas.begin();itm!=TMeas.end();itm++){ //loop on hits
337 
338  int i;
339  for(i=0;i<100;i++){Amplitudes[i]=0;}
340 
341  TanTrackAngle = -99;
343  ClSize = -99;
344  HitCharge = 0;
345  Type = -99;
346  Layer = -99;
347  Wheel = -99;
348  bw_fw = -99;
349  Ext_Int = -99;
350  MonoStereo = -99;
351  MagField = -99;
352  SignCorrection = -99;
353  XGlobal = -99;
354  YGlobal = -99;
355  ZGlobal = -99;
356  barycenter = -99;
357  hit_std_dev = -99;
358  sumx = 0;
359  id_detector=-1;
360  thick_detector=-1;
361  pitch_detector=-1;
362  HitNr = 1;
363 
364  SiStripRecHit2D lhit;
365  TrajectoryStateOnSurface tsos=itm->updatedState();
366  const TransientTrackingRecHit::ConstRecHitPointer thit=itm->recHit();
367  if((thit->geographicalId().subdetId() == int(StripSubdetector::TIB)) || thit->geographicalId().subdetId()== int(StripSubdetector::TOB)){ //include only barrel
368  const SiStripMatchedRecHit2D* matchedhit=dynamic_cast<const SiStripMatchedRecHit2D*>((*thit).hit());
369  const ProjectedSiStripRecHit2D* phit=dynamic_cast<const ProjectedSiStripRecHit2D*>((*thit).hit());
370  const SiStripRecHit2D* hit=dynamic_cast<const SiStripRecHit2D*>((*thit).hit());
371  if(phit) {lhit = phit->originalHit(); hit = &lhit;}
372 
373  LocalVector trackdirection=tsos.localDirection();
374 
375  if(matchedhit){//if matched hit...
376 
377  GluedGeomDet * gdet=(GluedGeomDet *)tracker->idToDet(matchedhit->geographicalId());
378 
379  GlobalVector gtrkdir=gdet->toGlobal(trackdirection);
380 
381  // THIS THE POINTER TO THE MONO HIT OF A MATCHED HIT
382 
383  // top be migrated to the more direct interface of matchedhit
384  cache.push_back(matchedhit->monoHit());
385  const SiStripRecHit2D * monohit = &cache.back();
386  const SiStripRecHit2D::ClusterRef & monocluster=monohit->cluster();
387  const GeomDetUnit * monodet=gdet->monoDet();
388  // this does not exists anymore! either project the matched or use CPE
389  const LocalPoint monoposition = monohit->localPosition();
390 
392  id_detector = detid.rawId();
394  const StripTopology& mtopol=(StripTopology&)monodet->topology();
395  pitch_detector = mtopol.localPitch(monoposition);
396  const GlobalPoint monogposition = (monodet->surface()).toGlobal(monoposition);
397  ClSize = (monocluster->amplitudes()).size();
398 
399  const std::vector<uint8_t> amplitudes = monocluster->amplitudes();
400 
401  barycenter = monocluster->barycenter()- 0.5;
402  uint16_t FirstStrip = monocluster->firstStrip();
403  std::vector<uint8_t>::const_iterator idigi;
404  std::vector<uint8_t>::const_iterator begin=amplitudes.begin();
405  nstrip=0;
406  for(idigi=begin; idigi!=amplitudes.end(); idigi++){
407  Amplitudes[nstrip]=*idigi;
408  sumx+=pow(((FirstStrip+idigi-begin)-barycenter),2)*(*idigi);
409  HitCharge+=*idigi;
410  }
412 
413 
414  XGlobal = monogposition.x();
415  YGlobal = monogposition.y();
416  ZGlobal = monogposition.z();
417 
418  Type = detid.subdetId();
419  MonoStereo=detid.stereo();
420 
421  if(detid.subdetId() == int (StripSubdetector::TIB)){
422 
423  Layer = tTopo->tibLayer(detid);
424  bw_fw = tTopo->tibStringInfo(detid)[0];
425  Ext_Int = tTopo->tibStringInfo(detid)[1];
426  }
427  if(detid.subdetId() == int (StripSubdetector::TOB)){
428 
429  Layer = tTopo->tobLayer(detid);
430  bw_fw = tTopo->tobRodInfo(detid)[0];
431  }
432  if(detid.subdetId() == int (StripSubdetector::TID)){
433 
434  Wheel = tTopo->tidWheel(detid);
435  bw_fw = tTopo->tidModuleInfo(detid)[0];
436  }
437  if(detid.subdetId() == int (StripSubdetector::TEC)){
438 
439  Wheel = tTopo->tecWheel(detid);
440  bw_fw = tTopo->tecPetalInfo(detid)[0];
441  }
442 
443 
444  LocalVector monotkdir=monodet->toLocal(gtrkdir);
445 
446  if(monotkdir.z()!=0){
447 
448  // THE LOCAL ANGLE (MONO)
449  float tanangle = monotkdir.x()/monotkdir.z();
450  TanTrackAngleParallel = monotkdir.y()/monotkdir.z();
451  TanTrackAngle = tanangle;
452  detparmap::iterator TheDet=detmap.find(detid.rawId());
453  LocalVector localmagdir;
454  if(TheDet!=detmap.end())localmagdir=TheDet->second->magfield;
455  MagField = localmagdir.mag();
456  if(MagField != 0.){
457  LocalVector monoylocal(0,1,0);
458  float signcorrection = (localmagdir * monoylocal)/(MagField);
459  if(signcorrection!=0)SignCorrection=1/signcorrection;
460  }
461 
462  std::map<const SiStripRecHit2D *,std::pair<float,float>,DetIdLess>::iterator alreadystored=hitangleassociation.find(monohit);
463 
464 
465  if(alreadystored != hitangleassociation.end()){//decide which hit take
466  if(itm->estimate() > alreadystored->second.first){
467  worse_double_hit++;}
468  if(itm->estimate() < alreadystored->second.first){
470  hitangleassociation.insert(std::make_pair(monohit, std::make_pair(itm->estimate(),tanangle)));
471 
472  }}
473  else{
474  hitangleassociation.insert(make_pair(monohit, std::make_pair(itm->estimate(),tanangle)));
475  HitsTree->Fill();
476  hitcounter++;}
477 
478  // THIS THE POINTER TO THE STEREO HIT OF A MATCHED HIT
479 
480  // top be migrated to the more direct interface of matchedhit
481  cache.push_back(matchedhit->stereoHit());
482  const SiStripRecHit2D * stereohit = &cache.back();
483  const SiStripRecHit2D::ClusterRef & stereocluster=stereohit->cluster();
484  const GeomDetUnit * stereodet=gdet->stereoDet();
485  // this does not exists anymore! either project the matched or use CPE
486  const LocalPoint stereoposition = stereohit->localPosition();
488  id_detector = detid.rawId();
490  const StripTopology& stopol=(StripTopology&)stereodet->topology();
491  pitch_detector = stopol.localPitch(stereoposition);
492  const GlobalPoint stereogposition = (stereodet->surface()).toGlobal(stereoposition);
493 
494  ClSize = (stereocluster->amplitudes()).size();
495 
496  const std::vector<uint8_t> amplitudes = stereocluster->amplitudes();
497 
498  barycenter = stereocluster->barycenter()- 0.5;
499  uint16_t FirstStrip = stereocluster->firstStrip();
500  std::vector<uint8_t>::const_iterator idigi;
501  std::vector<uint8_t>::const_iterator begin=amplitudes.begin();
502  nstrip=0;
503  for(idigi=begin; idigi!=amplitudes.end(); idigi++){
504  Amplitudes[nstrip]=*idigi;
505  sumx+=pow(((FirstStrip+idigi-begin)-barycenter),2)*(*idigi);
506  HitCharge+=*idigi;
507  }
509 
510  XGlobal = stereogposition.x();
511  YGlobal = stereogposition.y();
512  ZGlobal = stereogposition.z();
513 
514  Type = detid.subdetId();
515  MonoStereo=detid.stereo();
516 
517  if(detid.subdetId() == int (StripSubdetector::TIB)){
518 
519  Layer = tTopo->tibLayer(detid);
520  bw_fw = tTopo->tibStringInfo(detid)[0];
521  Ext_Int = tTopo->tibStringInfo(detid)[1];
522  }
523  if(detid.subdetId() == int (StripSubdetector::TOB)){
524 
525  Layer = tTopo->tobLayer(detid);
526  bw_fw = tTopo->tobRodInfo(detid)[0];
527  }
528  if(detid.subdetId() == int (StripSubdetector::TID)){
529 
530  Wheel = tTopo->tidWheel(detid);
531  bw_fw = tTopo->tidModuleInfo(detid)[0];
532  }
533  if(detid.subdetId() == int (StripSubdetector::TEC)){
534 
535  Wheel = tTopo->tecWheel(detid);
536  bw_fw = tTopo->tecPetalInfo(detid)[0];
537  }
538 
539 
540  LocalVector stereotkdir=stereodet->toLocal(gtrkdir);
541 
542  if(stereotkdir.z()!=0){
543 
544  // THE LOCAL ANGLE (STEREO)
545  float tanangle = stereotkdir.x()/stereotkdir.z();
546  TanTrackAngleParallel = stereotkdir.y()/stereotkdir.z();
547  TanTrackAngle = tanangle;
548  detparmap::iterator TheDet=detmap.find(detid.rawId());
549  LocalVector localmagdir;
550  if(TheDet!=detmap.end())localmagdir=TheDet->second->magfield;
551  MagField = localmagdir.mag();
552  LocalVector stereoylocal(0,1,0);
553  if(MagField != 0.){
554  float signcorrection = (localmagdir * stereoylocal)/(MagField);
555  if(signcorrection!=0)SignCorrection=1/signcorrection;}
556 
557  std::map<const SiStripRecHit2D *,std::pair<float,float>,DetIdLess>::iterator alreadystored=hitangleassociation.find(stereohit);
558 
559  if(alreadystored != hitangleassociation.end()){//decide which hit take
560  if(itm->estimate() > alreadystored->second.first){
561  worse_double_hit++;}
562  if(itm->estimate() < alreadystored->second.first){
564  hitangleassociation.insert(std::make_pair(stereohit, std::make_pair(itm->estimate(),tanangle)));
565 
566  }}
567  else{
568  hitangleassociation.insert(std::make_pair(stereohit, std::make_pair(itm->estimate(),tanangle)));
569  HitsTree->Fill();
570  hitcounter++;}
571 
572  }
573  }
574  }
575  else if(hit){
576 
577 
578  // hit= POINTER TO THE RECHIT
579 
580  const SiStripRecHit2D::ClusterRef & cluster=hit->cluster();
581 
583  const LocalPoint position = hit->localPosition();
585  id_detector = detid.rawId();
587  const StripTopology& topol=(StripTopology&)gdet->topology();
588  pitch_detector = topol.localPitch(position);
589  const GlobalPoint gposition = (gdet->surface()).toGlobal(position);
590 
591  ClSize = (cluster->amplitudes()).size();
592 
593  const std::vector<uint8_t> amplitudes = cluster->amplitudes();
594 
595  barycenter = cluster->barycenter()- 0.5;
596  uint16_t FirstStrip = cluster->firstStrip();
597  std::vector<uint8_t>::const_iterator idigi;
598  std::vector<uint8_t>::const_iterator begin=amplitudes.begin();
599  nstrip=0;
600  for(idigi=begin; idigi!=amplitudes.end(); idigi++){
601  Amplitudes[nstrip]=*idigi;
602  sumx+=pow(((FirstStrip+idigi-begin)-barycenter),2)*(*idigi);
603  HitCharge+=*idigi;
604  }
606 
607  XGlobal = gposition.x();
608  YGlobal = gposition.y();
609  ZGlobal = gposition.z();
610 
611  Type = detid.subdetId();
612  MonoStereo=detid.stereo();
613 
614  if(detid.subdetId() == int (StripSubdetector::TIB)){
615 
616  Layer = tTopo->tibLayer(detid);
617  bw_fw = tTopo->tibStringInfo(detid)[0];
618  Ext_Int = tTopo->tibStringInfo(detid)[1];
619  }
620  if(detid.subdetId() == int (StripSubdetector::TOB)){
621 
622  Layer = tTopo->tobLayer(detid);
623  bw_fw = tTopo->tobRodInfo(detid)[0];
624  }
625  if(detid.subdetId() == int (StripSubdetector::TID)){
626 
627  Wheel = tTopo->tidWheel(detid);
628  bw_fw = tTopo->tidModuleInfo(detid)[0];
629  }
630  if(detid.subdetId() == int (StripSubdetector::TEC)){
631 
632  Wheel = tTopo->tecWheel(detid);
633  bw_fw = tTopo->tecPetalInfo(detid)[0];
634  }
635 
636  if(trackdirection.z()!=0){
637 
638  // THE LOCAL ANGLE
639  float tanangle = trackdirection.x()/trackdirection.z();
640  TanTrackAngleParallel = trackdirection.y()/trackdirection.z();
641  TanTrackAngle = tanangle;
642  detparmap::iterator TheDet=detmap.find(detid.rawId());
643  LocalVector localmagdir;
644  if(TheDet!=detmap.end())localmagdir=TheDet->second->magfield;
645  MagField = localmagdir.mag();
646  if(MagField != 0.){
647  LocalVector ylocal(0,1,0);
648  float signcorrection = (localmagdir * ylocal)/(MagField);
649  if(signcorrection!=0)SignCorrection=1/signcorrection;}
650 
651  std::map<const SiStripRecHit2D *,std::pair<float,float>, DetIdLess>::iterator alreadystored=hitangleassociation.find(hit);
652 
653  if(alreadystored != hitangleassociation.end()){//decide which hit take
654  if(itm->estimate() > alreadystored->second.first){
655  worse_double_hit++;}
656  if(itm->estimate() < alreadystored->second.first){
658  hitangleassociation.insert(std::make_pair(hit, std::make_pair(itm->estimate(),tanangle)));
659 
660  }}
661  else{
662  hitangleassociation.insert(std::make_pair(hit,std::make_pair(itm->estimate(), tanangle) ) );
663  HitsTree->Fill();
664  hitcounter++;}
665 
666 
667  }
668  }
669  }
670  }
671  }
672  }
673  std::map<const SiStripRecHit2D *,std::pair<float,float>,DetIdLess>::iterator hitsiter;
674 
675 
676  for(hitsiter=hitangleassociation.begin();hitsiter!=hitangleassociation.end();hitsiter++){
677 
679 
680  const SiStripRecHit2D* hit=hitsiter->first;
681  const SiStripRecHit2D::ClusterRef & cluster=hit->cluster();
682 
683  size=(cluster->amplitudes()).size();
684 
686 
687  float tangent = hitsiter->second.second;
688 
689  //Sign and XZ plane projection correction applied in TrackLocalAngle (TIB|TOB layers)
690 
691  detparmap::iterator thedet=detmap.find(detid.rawId());
692  LocalVector localmagdir;
693  if(thedet!=detmap.end())localmagdir=thedet->second->magfield;
694  float localmagfield = localmagdir.mag();
695 
696  if(localmagfield != 0.){
697 
698  LocalVector ylocal(0,1,0);
699 
700  float normprojection = (localmagdir * ylocal)/(localmagfield);
701 
702  if(normprojection == 0.)LogDebug("SiStripLAProfileBooker::analyze")<<"Error: YBprojection = 0";
703 
704  else{
705  float signprojcorrection = 1/normprojection;
706  tangent*=signprojcorrection;
707  }
708  }
709 
710  //Filling histograms
711 
712  histomap::iterator thehisto=histos.find(detid.rawId());
713 
714  if(thehisto==histos.end())edm::LogError("SiStripLAProfileBooker::analyze")<<"Error: the profile associated to"<<detid.rawId()<<"does not exist! ";
715  else thehisto->second->Fill(tangent,size);
716 
717  //Summary histograms
719  unsigned int layerid;
720  getlayer(detid,tTopo,name,layerid);
721  histomap::iterator thesummaryhisto=summaryhisto.find(layerid);
722  if(thesummaryhisto==summaryhisto.end())edm::LogError("SiStripLAProfileBooker::analyze")<<"Error: the profile associated to subdet "<<name<<"does not exist! ";
723  else thesummaryhisto->second->Fill(tangent,size);
724 
725  }
726 
727 
728 }
#define LogDebug(id)
RunNumber_t run() const
Definition: EventID.h:42
T getParameter(std::string const &) const
EventNumber_t event() const
Definition: EventID.h:44
int i
Definition: DBlmapReader.cc:9
unsigned int tibLayer(const DetId &id) const
const GeomDetUnit * monoDet() const
Definition: GluedGeomDet.h:20
std::vector< unsigned int > tidModuleInfo(const DetId &id) const
LocalVector localDirection() const
GlobalPoint toGlobal(const Local2DPoint &lp) const
Conversion to the global R.F. from the R.F. of the GeomDet.
Definition: GeomDet.h:47
std::vector< Track > TrackCollection
collection of Tracks
Definition: TrackFwd.h:10
T y() const
Definition: PV3DBase.h:63
LocalPoint toLocal(const GlobalPoint &gp) const
Conversion to the R.F. of the GeomDet.
Definition: GeomDet.h:62
const Bounds & bounds() const
Definition: Surface.h:128
unsigned int tidWheel(const DetId &id) const
std::vector< unsigned int > tibStringInfo(const DetId &id) const
const Plane & surface() const
The nominal surface of the GeomDet.
Definition: GeomDet.h:35
uint32_t rawId() const
get the raw id
Definition: DetId.h:43
virtual float thickness() const =0
virtual const Topology & topology() const =0
virtual float localPitch(const LocalPoint &) const =0
std::vector< unsigned int > tecPetalInfo(const DetId &id) const
T mag() const
Definition: PV3DBase.h:67
std::vector< unsigned int > tobRodInfo(const DetId &id) const
std::shared_ptr< TrackingRecHit const > ConstRecHitPointer
void getlayer(const DetId &detid, const TrackerTopology *tTopo, std::string &name, unsigned int &layerid)
T sqrt(T t)
Definition: SSEVec.h:48
T z() const
Definition: PV3DBase.h:64
ClusterRef cluster() const
int subdetId() const
get the contents of the subdetector field (not cast into any detector&#39;s numbering enum) ...
Definition: DetId.h:37
virtual const GeomDet * idToDet(DetId) const
bool getByLabel(InputTag const &tag, Handle< PROD > &result) const
Definition: Event.h:390
SiStripRecHit2D originalHit() const
unsigned int stereo() const
stereo
SiStripRecHit2D stereoHit() const
tuple tracks
Definition: testEve_cfg.py:39
T const * product() const
Definition: ESHandle.h:62
std::vector< Trajectory > TrajectoryCollection
SiStripRecHit2D monoHit() const
edm::EventID id() const
Definition: EventBase.h:56
#define begin
Definition: vmac.h:30
static int position[264][3]
Definition: ReadPGInfo.cc:509
DetId geographicalId() const
T x() const
Definition: PV3DBase.h:62
virtual LocalPoint localPosition() const
unsigned int tecWheel(const DetId &id) const
Power< A, B >::type pow(const A &a, const B &b)
Definition: Power.h:40
const GeomDetUnit * stereoDet() const
Definition: GluedGeomDet.h:21
const Plane & specificSurface() const
Same as surface(), kept for backward compatibility.
Definition: GeomDet.h:38
unsigned int tobLayer(const DetId &id) const
const TrackerGeometry * tracker
void SiStripLAProfileBooker::beginRun ( const edm::EventSetup c)

Definition at line 71 of file SiStripLAProfileBooker.cc.

References Amplitudes, better_double_hit, DQMStore::bookProfile(), Surface::bounds(), bw_fw, chi2norm, ClSize, conf_, SiStripHistoId::createHistoId(), dbe_, detmap, EtaTrack, Event_Tree, eventcounter, EventCounter, EventNumber, EventTree, Ext_Int, edm::EventSetup::get(), getlayer(), edm::ParameterSet::getParameter(), edm::ParameterSet::getUntrackedParameter(), hFile, histos, hit_std_dev, Hit_Tree, HitCharge, hitcounter, hitcounter_2ndloop, HitNr, HitPerTrack, HitsTree, id_detector, TrackerGeometry::idToDet(), MagneticField::inTesla(), edm::InputTag::label(), Layer, StripTopology::localPitch(), SiStripLAProfileBooker::detparameters::magfield, MagField, HLT_ES_cff::magfield, Momentum, MonoStereo, mergeVDriftHistosByStation::name, cppFunctionSkipper::operator, AlCaHLTBitMon_ParallelJobs::p, ParticleCharge, PhiTrack, SiStripLAProfileBooker::detparameters::pitch, pitch_detector, edm::ESHandle< class >::product(), pt, DetId::rawId(), RunNumber, MonitorElement::setAxisTitle(), SiStripFolderOrganizer::setDetectorFolder(), SiStripFolderOrganizer::setSiStripFolder(), SignCorrection, GeomDet::specificSurface(), AlCaHLTBitMon_QueryRunRegistry::string, DetId::subdetId(), summarydetmap, summaryhisto, GeomDet::surface(), TanTrackAngle, TanTrackAngleParallel, StripSubdetector::TEC, thick_detector, Bounds::thickness(), SiStripLAProfileBooker::detparameters::thickness, StripSubdetector::TIB, StripSubdetector::TID, StripSubdetector::TOB, toLocal(), GeomDetUnit::topology(), Track_Tree, trackcollsize, TrackCounter, tracker, TrackTree, trajsize, Type, Wheel, worse_double_hit, XGlobal, YGlobal, and ZGlobal.

71  {
72 
73  //Retrieve tracker topology from geometry
75  c.get<IdealGeometryRecord>().get(tTopoHandle);
76  const TrackerTopology* const tTopo = tTopoHandle.product();
77 
78  //get magnetic field and geometry from ES
80  c.get<IdealMagneticFieldRecord>().get(esmagfield);
81  const MagneticField * magfield=&(*esmagfield);
82 
84  c.get<TrackerDigiGeometryRecord>().get(estracker);
85  tracker=&(*estracker);
86 
87  std::vector<uint32_t> activeDets;
88  edm::ESHandle<SiStripDetCabling> tkmechstruct=0;
89  if (conf_.getParameter<bool>("UseStripCablingDB")){
90  c.get<SiStripDetCablingRcd>().get(tkmechstruct);
91  activeDets.clear();
92  tkmechstruct->addActiveDetectorsRawIds(activeDets);
93  }
94  else {
95  const TrackerGeometry::DetIdContainer& Id = estracker->detIds();
96  TrackerGeometry::DetIdContainer::const_iterator Iditer;
97  activeDets.clear();
98  for(Iditer=Id.begin();Iditer!=Id.end();Iditer++){
99  activeDets.push_back(Iditer->rawId());
100  }
101  }
102 
103  edm::InputTag TkTag = conf_.getParameter<edm::InputTag>("Tracks");
104  //Get Ids;
105  double ModuleRangeMin=conf_.getParameter<double>("ModuleXMin");
106  double ModuleRangeMax=conf_.getParameter<double>("ModuleXMax");
107  double TIBRangeMin=conf_.getParameter<double>("TIBXMin");
108  double TIBRangeMax=conf_.getParameter<double>("TIBXMax");
109  double TOBRangeMin=conf_.getParameter<double>("TOBXMin");
110  double TOBRangeMax=conf_.getParameter<double>("TOBXMax");
111  int TIB_bin=conf_.getParameter<int>("TIB_bin");
112  int TOB_bin=conf_.getParameter<int>("TOB_bin");
113  int SUM_bin=conf_.getParameter<int>("SUM_bin");
114 
115  hFile = new TFile (conf_.getUntrackedParameter<std::string>("treeName").c_str(), "RECREATE" );
116 
117  Hit_Tree = hFile->mkdir("Hit_Tree");
118  Track_Tree = hFile->mkdir("Track_Tree");
119  Event_Tree = hFile->mkdir("Event_Tree");
120 
121  HitsTree = new TTree("HitsTree", "HitsTree");
122 
123  HitsTree->Branch("RunNumber", &RunNumber, "RunNumber/I");
124  HitsTree->Branch("EventNumber", &EventNumber, "EventNumber/I");
125  HitsTree->Branch("TanTrackAngle", &TanTrackAngle, "TanTrackAngle/F");
126  HitsTree->Branch("TanTrackAngleParallel", &TanTrackAngleParallel, "TanTrackAngleParallel/F");
127  HitsTree->Branch("ClSize", &ClSize, "ClSize/I");
128  HitsTree->Branch("HitCharge", &HitCharge, "HitCharge/I");
129  HitsTree->Branch("Hit_Std_Dev", &hit_std_dev, "hit_std_dev/F");
130  HitsTree->Branch("Type", &Type, "Type/I");
131  HitsTree->Branch("Layer", &Layer, "Layer/I");
132  HitsTree->Branch("Wheel", &Wheel, "Wheel/I");
133  HitsTree->Branch("bw_fw", &bw_fw, "bw_fw/I");
134  HitsTree->Branch("Ext_Int", &Ext_Int, "Ext_Int/I");
135  HitsTree->Branch("MonoStereo", &MonoStereo, "MonoStereo/I");
136  HitsTree->Branch("MagField", &MagField, "MagField/F");
137  HitsTree->Branch("SignCorrection", &SignCorrection, "SignCorrection/F");
138  HitsTree->Branch("XGlobal", &XGlobal, "XGlobal/F");
139  HitsTree->Branch("YGlobal", &YGlobal, "YGlobal/F");
140  HitsTree->Branch("ZGlobal", &ZGlobal, "ZGlobal/F");
141  HitsTree->Branch("ParticleCharge", &ParticleCharge, "ParticleCharge/I");
142  HitsTree->Branch("Momentum", &Momentum, "Momentum/F");
143  HitsTree->Branch("pt", &pt, "pt/F");
144  HitsTree->Branch("chi2norm", &chi2norm, "chi2norm/F");
145  HitsTree->Branch("EtaTrack", &EtaTrack, "EtaTrack/F");
146  HitsTree->Branch("PhiTrack", &PhiTrack, "PhiTrack/F");
147  HitsTree->Branch("TrajSize", &trajsize, "trajsize/I");
148  HitsTree->Branch("HitNr", &HitNr, "HitNr/I");
149  HitsTree->Branch("HitPerTrack", &HitPerTrack, "HitPerTrack/I");
150  HitsTree->Branch("id_detector", &id_detector, "id_detector/I");
151  HitsTree->Branch("thick_detector", &thick_detector, "thick_detector/F");
152  HitsTree->Branch("pitch_detector", &pitch_detector, "pitch_detector/F");
153  HitsTree->Branch("Amplitudes", Amplitudes, "Amplitudes[ClSize]/I");
154 
155  HitsTree->SetDirectory(Hit_Tree);
156 
157  TrackTree = new TTree("TrackTree", "TrackTree");
158 
159  TrackTree->Branch("TrackCounter", &TrackCounter, "TrackCounter/I");
160 
161  TrackTree->SetDirectory(Track_Tree);
162 
163  EventTree = new TTree("EventTree", "EventTree");
164 
165  EventTree->Branch("EventCounter", &EventCounter, "EventCounter/I");
166 
167  EventTree->SetDirectory(Event_Tree);
168 
169 
170  // use SistripHistoId for producing histogram id (and title)
171  SiStripHistoId hidmanager;
172 
173  // create SiStripFolderOrganizer
174  SiStripFolderOrganizer folder_organizer;
175 
177 
178  //get all detids
179 
180  for(std::vector<uint32_t>::const_iterator Id = activeDets.begin(); Id!=activeDets.end(); Id++){
181 
182  // for(Iditer=Id.begin();Iditer!=Id.end();Iditer++){ //loop on detids
183  DetId Iditero=DetId(*Id);
184  DetId *Iditer=&Iditero;
185  if((Iditer->subdetId() == int(StripSubdetector::TIB)) || (Iditer->subdetId() == int(StripSubdetector::TOB))){ //include only barrel
186 
187  int module_bin = 0;
188  if(Iditer->subdetId() == int(StripSubdetector::TIB)){
189  module_bin = TIB_bin;
190  }else{
191  module_bin = TOB_bin;
192  }
193 
194  // create a TProfile for each module
195  StripSubdetector subid(*Iditer);
196  std::string hid;
197  //Mono single sided detectors
198  LocalPoint p;
199  const GeomDetUnit * stripdet=dynamic_cast<const GeomDetUnit*>(tracker->idToDet(subid));
200  if(stripdet==0)continue;
201  const StripTopology& topol=(const StripTopology&)stripdet->topology();
202  float thickness=stripdet->specificSurface().bounds().thickness();
203 
204  folder_organizer.setDetectorFolder(Iditer->rawId(), tTopo);
205  hid = hidmanager.createHistoId(TkTag.label().c_str(),"det",Iditer->rawId());
206  MonitorElement * profile=dbe_->bookProfile(hid,hid,module_bin,ModuleRangeMin,ModuleRangeMax,20,0,5,"");
207  detparameters *param=new detparameters;
208  histos[Iditer->rawId()] = profile;
209  detmap[Iditer->rawId()] = param;
210  param->thickness = thickness*10000;
211  param->pitch = topol.localPitch(p)*10000;
212 
213  const GlobalPoint globalp = (stripdet->surface()).toGlobal(p);
214  GlobalVector globalmagdir = magfield->inTesla(globalp);
215  param->magfield=(stripdet->surface()).toLocal(globalmagdir);
216 
217  profile->setAxisTitle("tan(#theta_{t})",1);
218  profile->setAxisTitle("Cluster size",2);
219 
220  // create a summary histo if it does not exist already
222  unsigned int layerid;
223  getlayer(subid,tTopo,name,layerid);
224  name+=TkTag.label().c_str();
225  if(summaryhisto.find(layerid)==(summaryhisto.end())){
226  folder_organizer.setSiStripFolder();
227  MonitorElement * summaryprofile=0;
228  if (subid.subdetId()==int (StripSubdetector::TIB)||subid.subdetId()==int (StripSubdetector::TID))
229  summaryprofile=dbe_->bookProfile(name,name,SUM_bin,TIBRangeMin,TIBRangeMax,20,0,5,"");
230  else if (subid.subdetId()==int (StripSubdetector::TOB)||subid.subdetId()==int (StripSubdetector::TEC))
231  summaryprofile=dbe_->bookProfile(name,name,SUM_bin,TOBRangeMin,TOBRangeMax,20,0,5,"");
232  if(summaryprofile){
233  detparameters *summaryparam=new detparameters;
234  summaryhisto[layerid] = summaryprofile;
235  summarydetmap[layerid] = summaryparam;
236  summaryparam->thickness = thickness*10000;
237  summaryparam->pitch = topol.localPitch(p)*10000;
238  summaryprofile->setAxisTitle("tan(#theta_{t})",1);
239  summaryprofile->setAxisTitle("Cluster size",2);
240  }
241  }
242 
243  }
244  }
245 
246  trackcollsize = 0;
247  trajsize = 0;
248  RunNumber = 0;
249  EventNumber = 0;
250  hitcounter = 0;
251  hitcounter_2ndloop = 0;
252  worse_double_hit = 0;
253  better_double_hit = 0;
254  eventcounter = 0;
255 
256  EventCounter = 1;
257  TrackCounter = 1;
258 
259 }
T getParameter(std::string const &) const
T getUntrackedParameter(std::string const &, T const &) const
virtual GlobalVector inTesla(const GlobalPoint &gp) const =0
Field value ad specified global point, in Tesla.
tuple magfield
Definition: HLT_ES_cff.py:2311
const Bounds & bounds() const
Definition: Surface.h:128
LocalVector toLocal(const reco::Track::Vector &v, const Surface &s)
const Plane & surface() const
The nominal surface of the GeomDet.
Definition: GeomDet.h:35
uint32_t rawId() const
get the raw id
Definition: DetId.h:43
virtual float thickness() const =0
virtual const Topology & topology() const =0
virtual float localPitch(const LocalPoint &) const =0
void getlayer(const DetId &detid, const TrackerTopology *tTopo, std::string &name, unsigned int &layerid)
MonitorElement * bookProfile(const char *name, const char *title, int nchX, double lowX, double highX, int nchY, double lowY, double highY, const char *option="s")
Definition: DQMStore.cc:1268
int subdetId() const
get the contents of the subdetector field (not cast into any detector&#39;s numbering enum) ...
Definition: DetId.h:37
virtual const GeomDet * idToDet(DetId) const
Definition: DetId.h:18
const T & get() const
Definition: EventSetup.h:55
T const * product() const
Definition: ESHandle.h:62
std::string const & label() const
Definition: InputTag.h:42
void setAxisTitle(const std::string &title, int axis=1)
set x-, y- or z-axis title (axis=1, 2, 3 respectively)
const Plane & specificSurface() const
Same as surface(), kept for backward compatibility.
Definition: GeomDet.h:38
std::vector< DetId > DetIdContainer
const TrackerGeometry * tracker
void SiStripLAProfileBooker::endJob ( void  )
virtual

Reimplemented from edm::EDAnalyzer.

Definition at line 767 of file SiStripLAProfileBooker.cc.

References conf_, dbe_, edm::ParameterSet::getUntrackedParameter(), hFile, DQMStore::save(), and AlCaHLTBitMon_QueryRunRegistry::string.

767  {
768 
769  std::string outputFile_ =conf_.getUntrackedParameter<std::string>("fileName", "LorentzAngle.root");
770  dbe_->save(outputFile_);
771 
772  hFile->Write();
773  hFile->Close();
774 }
T getUntrackedParameter(std::string const &, T const &) const
void save(const std::string &filename, const std::string &path="", const std::string &pattern="", const std::string &rewrite="", const uint32_t run=0, const uint32_t lumi=0, SaveReferenceTag ref=SaveWithReference, int minStatus=dqm::qstatus::STATUS_OK, const std::string &fileupdate="RECREATE", const bool resetMEsAfterWriting=false)
Definition: DQMStore.cc:2540
void SiStripLAProfileBooker::getlayer ( const DetId detid,
const TrackerTopology tTopo,
std::string &  name,
unsigned int &  layerid 
)

Definition at line 733 of file SiStripLAProfileBooker.cc.

References DetId::subdetId(), StripSubdetector::TEC, TrackerTopology::tecRing(), StripSubdetector::TIB, TrackerTopology::tibLayer(), StripSubdetector::TID, TrackerTopology::tidRing(), StripSubdetector::TOB, and TrackerTopology::tobLayer().

Referenced by analyze(), and beginRun().

733  {
734  int layer=0;
735  std::stringstream layernum;
736 
737  if(detid.subdetId() == int (StripSubdetector::TIB)){
738 
739  name+="TIB_Layer_";
740  layer = tTopo->tibLayer(detid);
741  }
742 
743  else if(detid.subdetId() == int (StripSubdetector::TID)){
744 
745  name+="TID_Ring_";
746  layer = tTopo->tidRing(detid);
747  }
748 
749  else if(detid.subdetId() == int (StripSubdetector::TOB)){
750 
751  name+="TOB_Layer_";
752  layer = tTopo->tobLayer(detid);
753 
754  }
755 
756  else if(detid.subdetId() == int (StripSubdetector::TEC)){
757 
758  name+="TEC_Ring_";
759  layer = tTopo->tecRing(detid);
760  }
761  layernum<<layer;
762  name+=layernum.str();
763  layerid=detid.subdetId()*10+layer;
764 
765 }
unsigned int tibLayer(const DetId &id) const
unsigned int tidRing(const DetId &id) const
unsigned int tecRing(const DetId &id) const
ring id
int subdetId() const
get the contents of the subdetector field (not cast into any detector&#39;s numbering enum) ...
Definition: DetId.h:37
unsigned int tobLayer(const DetId &id) const

Member Data Documentation

uint8_t SiStripLAProfileBooker::Amplitudes[100]
private

Definition at line 57 of file SiStripLAProfileBooker.h.

Referenced by analyze(), and beginRun().

float SiStripLAProfileBooker::barycenter
private

Definition at line 53 of file SiStripLAProfileBooker.h.

Referenced by analyze().

int SiStripLAProfileBooker::better_double_hit
private

Definition at line 54 of file SiStripLAProfileBooker.h.

Referenced by analyze(), and beginRun().

int SiStripLAProfileBooker::bw_fw
private

Definition at line 52 of file SiStripLAProfileBooker.h.

Referenced by analyze(), and beginRun().

float SiStripLAProfileBooker::chi2norm
private

Definition at line 53 of file SiStripLAProfileBooker.h.

Referenced by analyze(), and beginRun().

int SiStripLAProfileBooker::ClSize
private

Definition at line 52 of file SiStripLAProfileBooker.h.

Referenced by analyze(), and beginRun().

edm::ParameterSet SiStripLAProfileBooker::conf_
private

Definition at line 71 of file SiStripLAProfileBooker.h.

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

DQMStore* SiStripLAProfileBooker::dbe_
private

Definition at line 67 of file SiStripLAProfileBooker.h.

Referenced by beginRun(), and endJob().

detparmap SiStripLAProfileBooker::detmap
private

Definition at line 69 of file SiStripLAProfileBooker.h.

Referenced by analyze(), beginRun(), and ~SiStripLAProfileBooker().

float SiStripLAProfileBooker::EtaTrack
private

Definition at line 53 of file SiStripLAProfileBooker.h.

Referenced by analyze(), and beginRun().

TDirectory * SiStripLAProfileBooker::Event_Tree
private

Definition at line 62 of file SiStripLAProfileBooker.h.

Referenced by beginRun().

int SiStripLAProfileBooker::eventcounter
private

Definition at line 54 of file SiStripLAProfileBooker.h.

Referenced by analyze(), and beginRun().

int SiStripLAProfileBooker::EventCounter
private

Definition at line 55 of file SiStripLAProfileBooker.h.

Referenced by beginRun().

int SiStripLAProfileBooker::EventNumber
private

Definition at line 52 of file SiStripLAProfileBooker.h.

Referenced by analyze(), and beginRun().

TTree * SiStripLAProfileBooker::EventTree
private

Definition at line 59 of file SiStripLAProfileBooker.h.

Referenced by analyze(), and beginRun().

int SiStripLAProfileBooker::Ext_Int
private

Definition at line 52 of file SiStripLAProfileBooker.h.

Referenced by analyze(), and beginRun().

TFile* SiStripLAProfileBooker::hFile
private

Definition at line 60 of file SiStripLAProfileBooker.h.

Referenced by beginRun(), endJob(), and ~SiStripLAProfileBooker().

histomap SiStripLAProfileBooker::histos
private
float SiStripLAProfileBooker::hit_std_dev
private

Definition at line 53 of file SiStripLAProfileBooker.h.

Referenced by analyze(), and beginRun().

TDirectory* SiStripLAProfileBooker::Hit_Tree
private

Definition at line 62 of file SiStripLAProfileBooker.h.

Referenced by beginRun().

int SiStripLAProfileBooker::HitCharge
private

Definition at line 52 of file SiStripLAProfileBooker.h.

Referenced by analyze(), and beginRun().

int SiStripLAProfileBooker::hitcounter
private

Definition at line 54 of file SiStripLAProfileBooker.h.

Referenced by analyze(), and beginRun().

int SiStripLAProfileBooker::hitcounter_2ndloop
private

Definition at line 54 of file SiStripLAProfileBooker.h.

Referenced by analyze(), and beginRun().

int SiStripLAProfileBooker::HitNr
private

Definition at line 54 of file SiStripLAProfileBooker.h.

Referenced by analyze(), and beginRun().

int SiStripLAProfileBooker::HitPerTrack
private

Definition at line 54 of file SiStripLAProfileBooker.h.

Referenced by analyze(), and beginRun().

TTree* SiStripLAProfileBooker::HitsTree
private

Definition at line 59 of file SiStripLAProfileBooker.h.

Referenced by analyze(), and beginRun().

int SiStripLAProfileBooker::id_detector
private

Definition at line 55 of file SiStripLAProfileBooker.h.

Referenced by analyze(), and beginRun().

int SiStripLAProfileBooker::Layer
private

Definition at line 52 of file SiStripLAProfileBooker.h.

Referenced by analyze(), and beginRun().

float SiStripLAProfileBooker::MagField
private

Definition at line 53 of file SiStripLAProfileBooker.h.

Referenced by analyze(), and beginRun().

float SiStripLAProfileBooker::Momentum
private

Definition at line 53 of file SiStripLAProfileBooker.h.

Referenced by analyze(), and beginRun().

int SiStripLAProfileBooker::MonoStereo
private

Definition at line 52 of file SiStripLAProfileBooker.h.

Referenced by analyze(), and beginRun().

int SiStripLAProfileBooker::nstrip
private

Definition at line 54 of file SiStripLAProfileBooker.h.

Referenced by analyze().

int SiStripLAProfileBooker::ParticleCharge
private

Definition at line 52 of file SiStripLAProfileBooker.h.

Referenced by analyze(), and beginRun().

float SiStripLAProfileBooker::PhiTrack
private

Definition at line 53 of file SiStripLAProfileBooker.h.

Referenced by analyze(), and beginRun().

float SiStripLAProfileBooker::pitch_detector
private

Definition at line 56 of file SiStripLAProfileBooker.h.

Referenced by analyze(), and beginRun().

float SiStripLAProfileBooker::pt
private

Definition at line 53 of file SiStripLAProfileBooker.h.

Referenced by analyze(), and beginRun().

int SiStripLAProfileBooker::RunNumber
private

Definition at line 52 of file SiStripLAProfileBooker.h.

Referenced by analyze(), beginRun(), and PopulateDB.PopulateDB::run().

float SiStripLAProfileBooker::SignCorrection
private

Definition at line 53 of file SiStripLAProfileBooker.h.

Referenced by analyze(), and beginRun().

int SiStripLAProfileBooker::size
private
detparmap SiStripLAProfileBooker::summarydetmap
private

Definition at line 70 of file SiStripLAProfileBooker.h.

Referenced by beginRun(), and ~SiStripLAProfileBooker().

histomap SiStripLAProfileBooker::summaryhisto
private

Definition at line 65 of file SiStripLAProfileBooker.h.

Referenced by analyze(), and beginRun().

float SiStripLAProfileBooker::sumx
private

Definition at line 53 of file SiStripLAProfileBooker.h.

Referenced by analyze().

float SiStripLAProfileBooker::TanTrackAngle
private

Definition at line 53 of file SiStripLAProfileBooker.h.

Referenced by analyze(), and beginRun().

float SiStripLAProfileBooker::TanTrackAngleParallel
private

Definition at line 53 of file SiStripLAProfileBooker.h.

Referenced by analyze(), and beginRun().

float SiStripLAProfileBooker::thick_detector
private

Definition at line 56 of file SiStripLAProfileBooker.h.

Referenced by analyze(), and beginRun().

TDirectory * SiStripLAProfileBooker::Track_Tree
private

Definition at line 62 of file SiStripLAProfileBooker.h.

Referenced by beginRun().

int SiStripLAProfileBooker::trackcollsize
private

Definition at line 52 of file SiStripLAProfileBooker.h.

Referenced by analyze(), and beginRun().

int SiStripLAProfileBooker::TrackCounter
private

Definition at line 55 of file SiStripLAProfileBooker.h.

Referenced by beginRun().

const TrackerGeometry* SiStripLAProfileBooker::tracker
private

Definition at line 74 of file SiStripLAProfileBooker.h.

Referenced by analyze(), and beginRun().

TTree * SiStripLAProfileBooker::TrackTree
private

Definition at line 59 of file SiStripLAProfileBooker.h.

Referenced by analyze(), and beginRun().

int SiStripLAProfileBooker::trajsize
private

Definition at line 52 of file SiStripLAProfileBooker.h.

Referenced by analyze(), and beginRun().

std::string SiStripLAProfileBooker::treename_
private

Definition at line 72 of file SiStripLAProfileBooker.h.

int SiStripLAProfileBooker::Type
private

Definition at line 52 of file SiStripLAProfileBooker.h.

Referenced by analyze(), beginRun(), and BeamSpotObj.BeamSpot::Reset().

int SiStripLAProfileBooker::Wheel
private

Definition at line 52 of file SiStripLAProfileBooker.h.

Referenced by analyze(), and beginRun().

int SiStripLAProfileBooker::worse_double_hit
private

Definition at line 54 of file SiStripLAProfileBooker.h.

Referenced by analyze(), and beginRun().

float SiStripLAProfileBooker::XGlobal
private

Definition at line 53 of file SiStripLAProfileBooker.h.

Referenced by analyze(), and beginRun().

float SiStripLAProfileBooker::YGlobal
private

Definition at line 53 of file SiStripLAProfileBooker.h.

Referenced by analyze(), and beginRun().

float SiStripLAProfileBooker::ZGlobal
private

Definition at line 53 of file SiStripLAProfileBooker.h.

Referenced by analyze(), and beginRun().