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
 
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, 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, 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(), LogDebug, PV3DBase< T, PVType, FrameType >::mag(), MagField, Momentum, GluedGeomDet::monoDet(), 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(), 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  TrajectoryStateOnSurface tsos=itm->updatedState();
365  const TransientTrackingRecHit::ConstRecHitPointer thit=itm->recHit();
366  if((thit->geographicalId().subdetId() == int(StripSubdetector::TIB)) || thit->geographicalId().subdetId()== int(StripSubdetector::TOB)){ //include only barrel
367  const SiStripMatchedRecHit2D* matchedhit=dynamic_cast<const SiStripMatchedRecHit2D*>((*thit).hit());
368  const ProjectedSiStripRecHit2D* phit=dynamic_cast<const ProjectedSiStripRecHit2D*>((*thit).hit());
369  const SiStripRecHit2D* hit=dynamic_cast<const SiStripRecHit2D*>((*thit).hit());
370  if(phit) hit=&(phit->originalHit());
371 
372  LocalVector trackdirection=tsos.localDirection();
373 
374  if(matchedhit){//if matched hit...
375 
376  GluedGeomDet * gdet=(GluedGeomDet *)tracker->idToDet(matchedhit->geographicalId());
377 
378  GlobalVector gtrkdir=gdet->toGlobal(trackdirection);
379 
380  // THIS THE POINTER TO THE MONO HIT OF A MATCHED HIT
381 
382  // top be migrated to the more direct interface of matchedhit
383  cache.push_back(matchedhit->monoHit());
384  const SiStripRecHit2D * monohit = &cache.back();
385  const SiStripRecHit2D::ClusterRef & monocluster=monohit->cluster();
386  const GeomDetUnit * monodet=gdet->monoDet();
387  // this does not exists anymore! either project the matched or use CPE
388  const LocalPoint monoposition = monohit->localPosition();
389 
390  StripSubdetector detid=(StripSubdetector)monohit->geographicalId();
391  id_detector = detid.rawId();
393  const StripTopology& mtopol=(StripTopology&)monodet->topology();
394  pitch_detector = mtopol.localPitch(monoposition);
395  const GlobalPoint monogposition = (monodet->surface()).toGlobal(monoposition);
396  ClSize = (monocluster->amplitudes()).size();
397 
398  const std::vector<uint8_t> amplitudes = monocluster->amplitudes();
399 
400  barycenter = monocluster->barycenter()- 0.5;
401  uint16_t FirstStrip = monocluster->firstStrip();
402  std::vector<uint8_t>::const_iterator idigi;
403  std::vector<uint8_t>::const_iterator begin=amplitudes.begin();
404  nstrip=0;
405  for(idigi=begin; idigi!=amplitudes.end(); idigi++){
406  Amplitudes[nstrip]=*idigi;
407  sumx+=pow(((FirstStrip+idigi-begin)-barycenter),2)*(*idigi);
408  HitCharge+=*idigi;
409  }
411 
412 
413  XGlobal = monogposition.x();
414  YGlobal = monogposition.y();
415  ZGlobal = monogposition.z();
416 
417  Type = detid.subdetId();
418  MonoStereo=detid.stereo();
419 
420  if(detid.subdetId() == int (StripSubdetector::TIB)){
421 
422  Layer = tTopo->tibLayer(detid);
423  bw_fw = tTopo->tibStringInfo(detid)[0];
424  Ext_Int = tTopo->tibStringInfo(detid)[1];
425  }
426  if(detid.subdetId() == int (StripSubdetector::TOB)){
427 
428  Layer = tTopo->tobLayer(detid);
429  bw_fw = tTopo->tobRodInfo(detid)[0];
430  }
431  if(detid.subdetId() == int (StripSubdetector::TID)){
432 
433  Wheel = tTopo->tidWheel(detid);
434  bw_fw = tTopo->tidModuleInfo(detid)[0];
435  }
436  if(detid.subdetId() == int (StripSubdetector::TEC)){
437 
438  Wheel = tTopo->tecWheel(detid);
439  bw_fw = tTopo->tecPetalInfo(detid)[0];
440  }
441 
442 
443  LocalVector monotkdir=monodet->toLocal(gtrkdir);
444 
445  if(monotkdir.z()!=0){
446 
447  // THE LOCAL ANGLE (MONO)
448  float tanangle = monotkdir.x()/monotkdir.z();
449  TanTrackAngleParallel = monotkdir.y()/monotkdir.z();
450  TanTrackAngle = tanangle;
451  detparmap::iterator TheDet=detmap.find(detid.rawId());
452  LocalVector localmagdir;
453  if(TheDet!=detmap.end())localmagdir=TheDet->second->magfield;
454  MagField = localmagdir.mag();
455  if(MagField != 0.){
456  LocalVector monoylocal(0,1,0);
457  float signcorrection = (localmagdir * monoylocal)/(MagField);
458  if(signcorrection!=0)SignCorrection=1/signcorrection;
459  }
460 
461  std::map<const SiStripRecHit2D *,std::pair<float,float>,DetIdLess>::iterator alreadystored=hitangleassociation.find(monohit);
462 
463 
464  if(alreadystored != hitangleassociation.end()){//decide which hit take
465  if(itm->estimate() > alreadystored->second.first){
466  worse_double_hit++;}
467  if(itm->estimate() < alreadystored->second.first){
469  hitangleassociation.insert(std::make_pair(monohit, std::make_pair(itm->estimate(),tanangle)));
470 
471  }}
472  else{
473  hitangleassociation.insert(make_pair(monohit, std::make_pair(itm->estimate(),tanangle)));
474  HitsTree->Fill();
475  hitcounter++;}
476 
477  // THIS THE POINTER TO THE STEREO HIT OF A MATCHED HIT
478 
479  // top be migrated to the more direct interface of matchedhit
480  cache.push_back(matchedhit->stereoHit());
481  const SiStripRecHit2D * stereohit = &cache.back();
482  const SiStripRecHit2D::ClusterRef & stereocluster=stereohit->cluster();
483  const GeomDetUnit * stereodet=gdet->stereoDet();
484  // this does not exists anymore! either project the matched or use CPE
485  const LocalPoint stereoposition = stereohit->localPosition();
486  StripSubdetector detid=(StripSubdetector)stereohit->geographicalId();
487  id_detector = detid.rawId();
489  const StripTopology& stopol=(StripTopology&)stereodet->topology();
490  pitch_detector = stopol.localPitch(stereoposition);
491  const GlobalPoint stereogposition = (stereodet->surface()).toGlobal(stereoposition);
492 
493  ClSize = (stereocluster->amplitudes()).size();
494 
495  const std::vector<uint8_t> amplitudes = stereocluster->amplitudes();
496 
497  barycenter = stereocluster->barycenter()- 0.5;
498  uint16_t FirstStrip = stereocluster->firstStrip();
499  std::vector<uint8_t>::const_iterator idigi;
500  std::vector<uint8_t>::const_iterator begin=amplitudes.begin();
501  nstrip=0;
502  for(idigi=begin; idigi!=amplitudes.end(); idigi++){
503  Amplitudes[nstrip]=*idigi;
504  sumx+=pow(((FirstStrip+idigi-begin)-barycenter),2)*(*idigi);
505  HitCharge+=*idigi;
506  }
508 
509  XGlobal = stereogposition.x();
510  YGlobal = stereogposition.y();
511  ZGlobal = stereogposition.z();
512 
513  Type = detid.subdetId();
514  MonoStereo=detid.stereo();
515 
516  if(detid.subdetId() == int (StripSubdetector::TIB)){
517 
518  Layer = tTopo->tibLayer(detid);
519  bw_fw = tTopo->tibStringInfo(detid)[0];
520  Ext_Int = tTopo->tibStringInfo(detid)[1];
521  }
522  if(detid.subdetId() == int (StripSubdetector::TOB)){
523 
524  Layer = tTopo->tobLayer(detid);
525  bw_fw = tTopo->tobRodInfo(detid)[0];
526  }
527  if(detid.subdetId() == int (StripSubdetector::TID)){
528 
529  Wheel = tTopo->tidWheel(detid);
530  bw_fw = tTopo->tidModuleInfo(detid)[0];
531  }
532  if(detid.subdetId() == int (StripSubdetector::TEC)){
533 
534  Wheel = tTopo->tecWheel(detid);
535  bw_fw = tTopo->tecPetalInfo(detid)[0];
536  }
537 
538 
539  LocalVector stereotkdir=stereodet->toLocal(gtrkdir);
540 
541  if(stereotkdir.z()!=0){
542 
543  // THE LOCAL ANGLE (STEREO)
544  float tanangle = stereotkdir.x()/stereotkdir.z();
545  TanTrackAngleParallel = stereotkdir.y()/stereotkdir.z();
546  TanTrackAngle = tanangle;
547  detparmap::iterator TheDet=detmap.find(detid.rawId());
548  LocalVector localmagdir;
549  if(TheDet!=detmap.end())localmagdir=TheDet->second->magfield;
550  MagField = localmagdir.mag();
551  LocalVector stereoylocal(0,1,0);
552  if(MagField != 0.){
553  float signcorrection = (localmagdir * stereoylocal)/(MagField);
554  if(signcorrection!=0)SignCorrection=1/signcorrection;}
555 
556  std::map<const SiStripRecHit2D *,std::pair<float,float>,DetIdLess>::iterator alreadystored=hitangleassociation.find(stereohit);
557 
558  if(alreadystored != hitangleassociation.end()){//decide which hit take
559  if(itm->estimate() > alreadystored->second.first){
560  worse_double_hit++;}
561  if(itm->estimate() < alreadystored->second.first){
563  hitangleassociation.insert(std::make_pair(stereohit, std::make_pair(itm->estimate(),tanangle)));
564 
565  }}
566  else{
567  hitangleassociation.insert(std::make_pair(stereohit, std::make_pair(itm->estimate(),tanangle)));
568  HitsTree->Fill();
569  hitcounter++;}
570 
571  }
572  }
573  }
574  else if(hit){
575 
576 
577  // hit= POINTER TO THE RECHIT
578 
579  const SiStripRecHit2D::ClusterRef & cluster=hit->cluster();
580 
581  GeomDetUnit * gdet=(GeomDetUnit *)tracker->idToDet(hit->geographicalId());
582  const LocalPoint position = hit->localPosition();
583  StripSubdetector detid=(StripSubdetector)hit->geographicalId();
584  id_detector = detid.rawId();
586  const StripTopology& topol=(StripTopology&)gdet->topology();
587  pitch_detector = topol.localPitch(position);
588  const GlobalPoint gposition = (gdet->surface()).toGlobal(position);
589 
590  ClSize = (cluster->amplitudes()).size();
591 
592  const std::vector<uint8_t> amplitudes = cluster->amplitudes();
593 
594  barycenter = cluster->barycenter()- 0.5;
595  uint16_t FirstStrip = cluster->firstStrip();
596  std::vector<uint8_t>::const_iterator idigi;
597  std::vector<uint8_t>::const_iterator begin=amplitudes.begin();
598  nstrip=0;
599  for(idigi=begin; idigi!=amplitudes.end(); idigi++){
600  Amplitudes[nstrip]=*idigi;
601  sumx+=pow(((FirstStrip+idigi-begin)-barycenter),2)*(*idigi);
602  HitCharge+=*idigi;
603  }
605 
606  XGlobal = gposition.x();
607  YGlobal = gposition.y();
608  ZGlobal = gposition.z();
609 
610  Type = detid.subdetId();
611  MonoStereo=detid.stereo();
612 
613  if(detid.subdetId() == int (StripSubdetector::TIB)){
614 
615  Layer = tTopo->tibLayer(detid);
616  bw_fw = tTopo->tibStringInfo(detid)[0];
617  Ext_Int = tTopo->tibStringInfo(detid)[1];
618  }
619  if(detid.subdetId() == int (StripSubdetector::TOB)){
620 
621  Layer = tTopo->tobLayer(detid);
622  bw_fw = tTopo->tobRodInfo(detid)[0];
623  }
624  if(detid.subdetId() == int (StripSubdetector::TID)){
625 
626  Wheel = tTopo->tidWheel(detid);
627  bw_fw = tTopo->tidModuleInfo(detid)[0];
628  }
629  if(detid.subdetId() == int (StripSubdetector::TEC)){
630 
631  Wheel = tTopo->tecWheel(detid);
632  bw_fw = tTopo->tecPetalInfo(detid)[0];
633  }
634 
635  if(trackdirection.z()!=0){
636 
637  // THE LOCAL ANGLE
638  float tanangle = trackdirection.x()/trackdirection.z();
639  TanTrackAngleParallel = trackdirection.y()/trackdirection.z();
640  TanTrackAngle = tanangle;
641  detparmap::iterator TheDet=detmap.find(detid.rawId());
642  LocalVector localmagdir;
643  if(TheDet!=detmap.end())localmagdir=TheDet->second->magfield;
644  MagField = localmagdir.mag();
645  if(MagField != 0.){
646  LocalVector ylocal(0,1,0);
647  float signcorrection = (localmagdir * ylocal)/(MagField);
648  if(signcorrection!=0)SignCorrection=1/signcorrection;}
649 
650  std::map<const SiStripRecHit2D *,std::pair<float,float>, DetIdLess>::iterator alreadystored=hitangleassociation.find(hit);
651 
652  if(alreadystored != hitangleassociation.end()){//decide which hit take
653  if(itm->estimate() > alreadystored->second.first){
654  worse_double_hit++;}
655  if(itm->estimate() < alreadystored->second.first){
657  hitangleassociation.insert(std::make_pair(hit, std::make_pair(itm->estimate(),tanangle)));
658 
659  }}
660  else{
661  hitangleassociation.insert(std::make_pair(hit,std::make_pair(itm->estimate(), tanangle) ) );
662  HitsTree->Fill();
663  hitcounter++;}
664 
665 
666  }
667  }
668  }
669  }
670  }
671  }
672  std::map<const SiStripRecHit2D *,std::pair<float,float>,DetIdLess>::iterator hitsiter;
673 
674 
675  for(hitsiter=hitangleassociation.begin();hitsiter!=hitangleassociation.end();hitsiter++){
676 
678 
679  const SiStripRecHit2D* hit=hitsiter->first;
680  const SiStripRecHit2D::ClusterRef & cluster=hit->cluster();
681 
682  size=(cluster->amplitudes()).size();
683 
684  StripSubdetector detid=(StripSubdetector)hit->geographicalId();
685 
686  float tangent = hitsiter->second.second;
687 
688  //Sign and XZ plane projection correction applied in TrackLocalAngle (TIB|TOB layers)
689 
690  detparmap::iterator thedet=detmap.find(detid.rawId());
691  LocalVector localmagdir;
692  if(thedet!=detmap.end())localmagdir=thedet->second->magfield;
693  float localmagfield = localmagdir.mag();
694 
695  if(localmagfield != 0.){
696 
697  LocalVector ylocal(0,1,0);
698 
699  float normprojection = (localmagdir * ylocal)/(localmagfield);
700 
701  if(normprojection == 0.)LogDebug("SiStripLAProfileBooker::analyze")<<"Error: YBprojection = 0";
702 
703  else{
704  float signprojcorrection = 1/normprojection;
705  tangent*=signprojcorrection;
706  }
707  }
708 
709  //Filling histograms
710 
711  histomap::iterator thehisto=histos.find(detid.rawId());
712 
713  if(thehisto==histos.end())edm::LogError("SiStripLAProfileBooker::analyze")<<"Error: the profile associated to"<<detid.rawId()<<"does not exist! ";
714  else thehisto->second->Fill(tangent,size);
715 
716  //Summary histograms
718  unsigned int layerid;
719  getlayer(detid,tTopo,name,layerid);
720  histomap::iterator thesummaryhisto=summaryhisto.find(layerid);
721  if(thesummaryhisto==summaryhisto.end())edm::LogError("SiStripLAProfileBooker::analyze")<<"Error: the profile associated to subdet "<<name<<"does not exist! ";
722  else thesummaryhisto->second->Fill(tangent,size);
723 
724  }
725 
726 
727 }
#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
static int position[TOTALCHAMBERS][3]
Definition: ReadPGInfo.cc:509
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
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
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
unsigned int stereo() const
stereo
tuple tracks
Definition: testEve_cfg.py:39
T const * product() const
Definition: ESHandle.h:62
std::vector< Trajectory > TrajectoryCollection
edm::EventID id() const
Definition: EventBase.h:56
#define begin
Definition: vmac.h:30
T x() const
Definition: PV3DBase.h:62
unsigned int tecWheel(const DetId &id) const
const SiStripRecHit2D & originalHit() 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:1187
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 766 of file SiStripLAProfileBooker.cc.

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

766  {
767 
768  std::string outputFile_ =conf_.getUntrackedParameter<std::string>("fileName", "LorentzAngle.root");
769  dbe_->save(outputFile_);
770 
771  hFile->Write();
772  hFile->Close();
773 }
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, SaveReferenceTag ref=SaveWithReference, int minStatus=dqm::qstatus::STATUS_OK, const std::string &fileupdate="RECREATE")
Definition: DQMStore.cc:2297
void SiStripLAProfileBooker::getlayer ( const DetId detid,
const TrackerTopology tTopo,
std::string &  name,
unsigned int &  layerid 
)

Definition at line 732 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().

732  {
733  int layer=0;
734  std::stringstream layernum;
735 
736  if(detid.subdetId() == int (StripSubdetector::TIB)){
737 
738  name+="TIB_Layer_";
739  layer = tTopo->tibLayer(detid);
740  }
741 
742  else if(detid.subdetId() == int (StripSubdetector::TID)){
743 
744  name+="TID_Ring_";
745  layer = tTopo->tidRing(detid);
746  }
747 
748  else if(detid.subdetId() == int (StripSubdetector::TOB)){
749 
750  name+="TOB_Layer_";
751  layer = tTopo->tobLayer(detid);
752 
753  }
754 
755  else if(detid.subdetId() == int (StripSubdetector::TEC)){
756 
757  name+="TEC_Ring_";
758  layer = tTopo->tecRing(detid);
759  }
760  layernum<<layer;
761  name+=layernum.str();
762  layerid=detid.subdetId()*10+layer;
763 
764 }
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

Definition at line 54 of file SiStripLAProfileBooker.h.

Referenced by cuy.FindIssue::__init__(), and analyze().

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().