CMS 3D CMS Logo
 All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Properties Friends Macros Pages
Namespaces | Classes | Functions | Variables
fireworks Namespace Reference

Namespaces

 expression
 
 geometry
 
 table
 

Classes

class  Context
 
struct  jetScaleMarker
 
class  OptionNode
 
struct  OptionNodePtrCompare
 
struct  OptionNodePtrEqual
 
struct  scaleMarker
 
struct  State
 
class  StateOrdering
 

Functions

bool acceptDataFormatsVersion (TString &n)
 
void addBox (const std::vector< float > &corners, TEveElement *, FWProxyBuilderBase *)
 
void addCircle (double eta, double phi, double radius, const unsigned int nLineSegments, TEveElement *comp, FWProxyBuilderBase *pb)
 
void addDashedArrow (double phi, double size, TEveElement *comp, FWProxyBuilderBase *pb)
 
void addDashedLine (double phi, double theta, double size, TEveElement *comp, FWProxyBuilderBase *pb)
 
void addDoubleLines (double phi, TEveElement *comp, FWProxyBuilderBase *pb)
 
void addRhoZEnergyProjection (FWProxyBuilderBase *, TEveElement *, double r_ecal, double z_ecal, double theta_min, double theta_max, double phi)
 
void addSiStripClusters (const FWEventItem *iItem, const reco::Track &t, class TEveElement *tList, bool addNearbyClusters, bool master)
 
void addStraightLineSegment (TEveStraightLineSet *marker, reco::Candidate const *cand, double scale_factor=2)
 
const TGGC & boldGC ()
 
void createSegment (int detector, bool matchedSegment, float segmentLength, float segmentLimit, float *segmentPosition, float *segmentDirection, float *segmentInnerPoint, float *segmentOuterPoint)
 
const TString datadir ("/src/Fireworks/Core/")
 
void drawEnergyScaledBox3D (const float *corners, float scale, TEveElement *, FWProxyBuilderBase *, bool invert=false)
 
void drawEnergyTower3D (const float *corners, float scale, TEveElement *, FWProxyBuilderBase *, bool reflect=false)
 
void drawEtScaledBox3D (const float *corners, float energy, float maxEnergy, TEveElement *, FWProxyBuilderBase *, bool reflect=false)
 
void drawEtTower3D (const float *corners, float scale, TEveElement *, FWProxyBuilderBase *, bool reflect=false)
 
void energyScaledBox3DCorners (const float *corners, float scale, std::vector< float > &, bool invert=false)
 
void energyTower3DCorners (const float *corners, float scale, std::vector< float > &, bool reflect=false)
 
void etScaledBox3DCorners (const float *corners, float energy, float maxEnergy, std::vector< float > &scaledCorners, bool reflect=false)
 
void etTower3DCorners (const float *corners, float scale, std::vector< float > &, bool reflect=false)
 
const SiStripCluster * extractClusterFromTrackingRecHit (const TrackingRecHit *rh)
 
void getDecomposedVersion (const TString &s, int *out)
 
std::string getLocalTime (const edm::EventBase &event)
 
std::pair< double, double > getPhiRange (const std::vector< double > &phis, double phi)
 
TEveGeoShape * getShape (const char *name, TGeoBBox *shape, Color_t color)
 
std::string getTimeGMT (const edm::EventBase &event)
 
const TGGC & greenGC ()
 
std::string info (const DetId &)
 
std::string info (const std::set< DetId > &)
 
std::string info (const std::vector< DetId > &)
 
void invertBox (std::vector< float > &corners)
 
const TGGC & italicGC ()
 
void localSiStrip (short strip, float *localTop, float *localBottom, const float *pars, unsigned int id)
 
TEveJetCone * makeEveJetCone (const reco::Jet &iData, const fireworks::Context &context)
 
bool makeRhoPhiSuperCluster (FWProxyBuilderBase *, const reco::SuperClusterRef &iCluster, float iPhi, TEveElement &oItemHolder)
 
bool makeRhoZSuperCluster (FWProxyBuilderBase *, const reco::SuperClusterRef &iCluster, float iPhi, TEveElement &oItemHolder)
 
float pixelLocalX (const double mpx, const int m_nrows)
 
float pixelLocalY (const double mpy, const int m_ncols)
 
TEveTrack * prepareCandidate (const reco::Candidate &track, TEveTrackPropagator *propagator)
 
TEveTrack * prepareTrack (const reco::Track &track, TEveTrackPropagator *propagator, const std::vector< TEveVector > &extraRefPoints=std::vector< TEveVector >())
 
void pushNearbyPixelHits (std::vector< TVector3 > &pixelPoints, const FWEventItem &iItem, const reco::Track &t)
 
void pushPixelCluster (std::vector< TVector3 > &pixelPoints, const FWGeometry &geom, DetId id, const SiPixelCluster &c, const float *pars)
 
void pushPixelHits (std::vector< TVector3 > &pixelPoints, const FWEventItem &iItem, const reco::Track &t)
 
const TGGC & redGC ()
 
void setPath (TString &v)
 
void setTrackTypePF (const reco::PFCandidate &pfCand, TAttLine *track)
 
void setupAddElement (TEveElement *el, TEveElement *parent, const FWEventItem *item, bool master, bool color)
 
int * supportedDataFormatsVersion ()
 
std::vector< std::pair< double,
double > > 		thetaBins ()
 

Variables

static const int BIG_PIX_PER_ROC_X = 1
 
static const int BIG_PIX_PER_ROC_Y = 2
 
static const int COLS_PER_ROC = 52
 
static const double MICRON = 1./1000./10.
 
const TString p1 = gSystem->Getenv("CMSSW_BASE") + datadir
 
const TString p2 = gSystem->Getenv("CMSSW_RELEASE_BASE") + datadir
 
static const double PITCHX = 100*MICRON
 
static const double PITCHY = 150*MICRON
 
static const int ROWS_PER_ROC = 80
 
static const std::string subdets [7] = { "UNKNOWN", "PXB", "PXF", "TIB", "TID", "TOB", "TEC" }
 

Function Documentation

bool fireworks::acceptDataFormatsVersion ( TString &  n)

Definition at line 71 of file fwPaths.cc.

References data, getDecomposedVersion(), and supportedDataFormatsVersion().

Referenced by FWFileEntry::openFile().

72 {
73  int data[] = {0, 0, 0};
74  getDecomposedVersion(processConfigurationVersion, data);
75 
76 
77  int* running = supportedDataFormatsVersion();
78  if ((data[0] == 6 && running[0] == 5 && running[1] > 1) ||
79  (data[0] == 5 && data[1] > 1 && running[0] == 6))
80  return true;
81  else
82  return data[0] == running[0];
83 }
fireworks::getDecomposedVersion
void getDecomposedVersion(const TString &s, int *out)
Definition: fwPaths.cc:30
fireworks::supportedDataFormatsVersion
int * supportedDataFormatsVersion()
Definition: fwPaths.cc:43
data
char data[epos_bytes_allocation]
Definition: EPOS_Wrapper.h:82
void fireworks::addBox ( const std::vector< float > &  corners,
TEveElement *  comp,
FWProxyBuilderBase *  pb 
)

Definition at line 146 of file BuilderUtils.cc.

References FWProxyBuilderBase::setupAddElement().

Referenced by drawEnergyScaledBox3D(), drawEnergyTower3D(), drawEtScaledBox3D(), and drawEtTower3D().

147  {
148  TEveBox* eveBox = new TEveBox( "Box" );
149  eveBox->SetDrawFrame( false );
150  eveBox->SetPickable( true );
151  eveBox->SetVertices( &corners[0] );
152 
153  pb->setupAddElement( eveBox, comp );
154  }
FWProxyBuilderBase::setupAddElement
void setupAddElement(TEveElement *el, TEveElement *parent, bool set_color=true) const
Definition: FWProxyBuilderBase.cc:429
AlCaHLTBitMon_QueryRunRegistry.comp
string comp
Definition: AlCaHLTBitMon_QueryRunRegistry.py:248
void fireworks::addCircle ( double  eta,
double  phi,
double  radius,
const unsigned int  nLineSegments,
TEveElement *  comp,
FWProxyBuilderBase *  pb 
)

Definition at line 156 of file BuilderUtils.cc.

References funct::cos(), M_PI, FWProxyBuilderBase::setupAddElement(), and funct::sin().

Referenced by FWJetLegoProxyBuilder::build(), FWL1MuonParticleLegoProxyBuilder::build(), FWL1JetParticleLegoProxyBuilder::build(), FWPFTauProxyBuilder::buildViewType(), and FWCaloTauProxyBuilder::buildViewType().

157  {
158  TEveStraightLineSet* container = new TEveStraightLineSet;
159 
160  for( unsigned int iphi = 0; iphi < nLineSegments; ++iphi )
161  {
162  container->AddLine( eta + radius * cos( 2 * M_PI / nLineSegments * iphi ),
163  phi + radius * sin( 2 * M_PI / nLineSegments * iphi ),
164  0.01,
165  eta + radius * cos( 2 * M_PI / nLineSegments * ( iphi + 1 )),
166  phi + radius * sin( 2 * M_PI / nLineSegments * ( iphi + 1 )),
167  0.01 );
168  }
169  pb->setupAddElement( container, comp );
170  }
FWProxyBuilderBase::setupAddElement
void setupAddElement(TEveElement *el, TEveElement *parent, bool set_color=true) const
Definition: FWProxyBuilderBase.cc:429
funct::sin
Sin< T >::type sin(const T &t)
Definition: Sin.h:22
eta
T eta() const
Definition: Basic3DVectorLD.h:177
AlCaHLTBitMon_QueryRunRegistry.comp
string comp
Definition: AlCaHLTBitMon_QueryRunRegistry.py:248
funct::cos
Cos< T >::type cos(const T &t)
Definition: Cos.h:22
M_PI
#define M_PI
Definition: readTestVector.cc:25
CosmicsPD_Skims.radius
tuple radius
Definition: CosmicsPD_Skims.py:135
phi
Definition: DDAxes.h:10
void fireworks::addDashedArrow ( double  phi,
double  size,
TEveElement *  comp,
FWProxyBuilderBase *  pb 
)

Definition at line 172 of file BuilderUtils.cc.

References funct::cos(), cmsHarvester::marker, FWProxyBuilderBase::setupAddElement(), and funct::sin().

Referenced by FWL1EtMissParticleGlimpseProxyBuilder::build(), and FWMETGlimpseProxyBuilder::build().

173  {
174  TEveScalableStraightLineSet* marker = new TEveScalableStraightLineSet;
175  marker->SetLineWidth( 1 );
176  marker->SetLineStyle( 2 );
177  marker->AddLine( 0, 0, 0, size * cos( phi ), size * sin( phi ), 0 );
178  marker->AddLine( size * 0.9 * cos( phi + 0.03 ), size * 0.9 * sin( phi + 0.03 ), 0, size * cos( phi ), size * sin( phi ), 0 );
179  marker->AddLine( size * 0.9 * cos( phi - 0.03 ), size * 0.9 * sin( phi - 0.03 ), 0, size * cos( phi ), size * sin( phi ), 0 );
180  pb->setupAddElement( marker, comp );
181  }
FWProxyBuilderBase::setupAddElement
void setupAddElement(TEveElement *el, TEveElement *parent, bool set_color=true) const
Definition: FWProxyBuilderBase.cc:429
funct::sin
Sin< T >::type sin(const T &t)
Definition: Sin.h:22
AlCaHLTBitMon_QueryRunRegistry.comp
string comp
Definition: AlCaHLTBitMon_QueryRunRegistry.py:248
funct::cos
Cos< T >::type cos(const T &t)
Definition: Cos.h:22
findQualityFiles.size
tuple size
Write out results.
Definition: findQualityFiles.py:442
phi
Definition: DDAxes.h:10
void fireworks::addDashedLine ( double  phi,
double  theta,
double  size,
TEveElement *  comp,
FWProxyBuilderBase *  pb 
)

Definition at line 183 of file BuilderUtils.cc.

References fireworks::Context::caloR1(), fireworks::Context::caloZ2(), FWProxyBuilderBase::context(), funct::cos(), M_PI, cmsHarvester::marker, alignCSCRings::r, FWProxyBuilderBase::setupAddElement(), and funct::sin().

Referenced by FWL1EmParticleProxyBuilder::build(), FWL1MuonParticleProxyBuilder::build(), FWL1JetParticleProxyBuilder::build(), and FWL1EtMissParticleProxyBuilder::build().

184  {
185  double r( 0 );
186  if( theta < pb->context().caloTransAngle() || M_PI - theta < pb->context().caloTransAngle())
187  r = pb->context().caloZ2() / fabs( cos( theta ));
188  else
189  r = pb->context().caloR1() / sin( theta );
190 
191  TEveStraightLineSet* marker = new TEveStraightLineSet;
192  marker->SetLineWidth( 2 );
193  marker->SetLineStyle( 2 );
194  marker->AddLine( r * cos( phi ) * sin( theta ), r * sin( phi ) * sin( theta ), r * cos( theta ),
195  ( r + size ) * cos( phi ) * sin( theta ), ( r + size ) * sin( phi ) * sin( theta ), ( r + size ) * cos( theta ));
196  pb->setupAddElement( marker, comp );
197  }
FWProxyBuilderBase::context
const fireworks::Context & context() const
Definition: FWProxyBuilderBase.cc:503
FWProxyBuilderBase::setupAddElement
void setupAddElement(TEveElement *el, TEveElement *parent, bool set_color=true) const
Definition: FWProxyBuilderBase.cc:429
funct::sin
Sin< T >::type sin(const T &t)
Definition: Sin.h:22
theta
Geom::Theta< T > theta() const
Definition: Basic3DVectorLD.h:170
fireworks::Context::caloZ2
static float caloZ2(bool offset=true)
Definition: Context.cc:224
AlCaHLTBitMon_QueryRunRegistry.comp
string comp
Definition: AlCaHLTBitMon_QueryRunRegistry.py:248
funct::cos
Cos< T >::type cos(const T &t)
Definition: Cos.h:22
fireworks::Context::caloR1
static float caloR1(bool offset=true)
Definition: Context.cc:209
M_PI
#define M_PI
Definition: readTestVector.cc:25
findQualityFiles.size
tuple size
Write out results.
Definition: findQualityFiles.py:442
phi
Definition: DDAxes.h:10
void fireworks::addDoubleLines ( double  phi,
TEveElement *  comp,
FWProxyBuilderBase *  pb 
)

Definition at line 199 of file BuilderUtils.cc.

References M_PI, and FWProxyBuilderBase::setupAddElement().

Referenced by FWL1EtMissParticleLegoProxyBuilder::build(), and FWMETLegoProxyBuilder::build().

200  {
201  TEveStraightLineSet* mainLine = new TEveStraightLineSet;
202  mainLine->AddLine( -5.191, phi, 0.01, 5.191, phi, 0.01 );
203  pb->setupAddElement( mainLine, comp );
204 
205  phi = phi > 0 ? phi - M_PI : phi + M_PI;
206  TEveStraightLineSet* secondLine = new TEveStraightLineSet;
207  secondLine->SetLineStyle( 7 );
208  secondLine->AddLine( -5.191, phi, 0.01, 5.191, phi, 0.01 );
209  pb->setupAddElement( secondLine, comp );
210  }
FWProxyBuilderBase::setupAddElement
void setupAddElement(TEveElement *el, TEveElement *parent, bool set_color=true) const
Definition: FWProxyBuilderBase.cc:429
AlCaHLTBitMon_QueryRunRegistry.comp
string comp
Definition: AlCaHLTBitMon_QueryRunRegistry.py:248
M_PI
#define M_PI
Definition: readTestVector.cc:25
phi
Definition: DDAxes.h:10
void fireworks::addRhoZEnergyProjection ( FWProxyBuilderBase *  pb,
TEveElement *  container,
double  r_ecal,
double  z_ecal,
double  theta_min,
double  theta_max,
double  phi 
)

Definition at line 60 of file BuilderUtils.cc.

References diffTwoXMLs::r1, diffTwoXMLs::r2, FWProxyBuilderBase::setupAddElement(), edmStreamStallGrapher::t, and funct::tan().

Referenced by FWTauProxyBuilderBase::buildBaseTau(), FWMETProxyBuilder::buildViewType(), and makeRhoZSuperCluster().

64  {
65  TEveGeoManagerHolder gmgr( TEveGeoShape::GetGeoMangeur());
66  double z1 = r_ecal / tan( theta_min );
67  if( z1 > z_ecal ) z1 = z_ecal;
68  if( z1 < -z_ecal ) z1 = -z_ecal;
69  double z2 = r_ecal / tan( theta_max );
70  if( z2 > z_ecal ) z2 = z_ecal;
71  if( z2 < -z_ecal ) z2 = -z_ecal;
72  double r1 = z_ecal * fabs( tan( theta_min ));
73  if( r1 > r_ecal ) r1 = r_ecal;
74  if( phi < 0 ) r1 = -r1;
75  double r2 = z_ecal * fabs( tan( theta_max ));
76  if( r2 > r_ecal ) r2 = r_ecal;
77  if( phi < 0 ) r2 = -r2;
78 
79  if( fabs(r2 - r1) > 1 )
80  {
81  TGeoBBox *sc_box = new TGeoBBox( 0., fabs( r2 - r1 ) / 2, 1 );
82  TEveGeoShape *element = new TEveGeoShape("r-segment");
83  element->SetShape(sc_box);
84  TEveTrans &t = element->RefMainTrans();
85  t(1,4) = 0;
86  t(2,4) = (r2+r1)/2;
87  t(3,4) = fabs(z2)>fabs(z1) ? z2 : z1;
88  pb->setupAddElement(element, container);
89  }
90  if( fabs(z2 - z1) > 1 )
91  {
92  TGeoBBox *sc_box = new TGeoBBox( 0., 1, ( z2 - z1 ) / 2 );
93  TEveGeoShape *element = new TEveGeoShape("z-segment");
94  element->SetShape( sc_box );
95  TEveTrans &t = element->RefMainTrans();
96  t(1,4) = 0;
97  t(2,4) = fabs(r2)>fabs(r1) ? r2 : r1;
98  t(3,4) = (z2+z1)/2;
99  pb->setupAddElement(element, container);
100  }
101  }
FWProxyBuilderBase::setupAddElement
void setupAddElement(TEveElement *el, TEveElement *parent, bool set_color=true) const
Definition: FWProxyBuilderBase.cc:429
funct::tan
Tan< T >::type tan(const T &t)
Definition: Tan.h:22
phi
Definition: DDAxes.h:10
void fireworks::addSiStripClusters ( const FWEventItem *  iItem,
const reco::Track &  t,
class TEveElement *  tList,
bool  addNearbyClusters,
bool  master 
)

Definition at line 381 of file TrackUtils.cc.

References edmNew::DetSet< T >::begin(), SiStripRecHit1D::cluster(), SiStripRecHit2D::cluster(), FWGeometry::contains(), edmNew::DetSet< T >::end(), edmNew::DetSetVector< T >::end(), extractClusterFromTrackingRecHit(), edmNew::DetSetVector< T >::find(), SiStripCluster::firstStrip(), fwLog, TrackingRecHit::geographicalId(), relativeConstraints::geom, FWEventItem::getGeom(), FWGeometry::getParameters(), edm::Ref< C, T, F >::isAvailable(), edm::Ref< C, T, F >::isNonnull(), TrackingRecHit::isValid(), fwlog::kDebug, fwlog::kError, localSiStrip(), FWGeometry::localToGlobal(), edm::Ref< C, T, F >::product(), DetId::rawId(), reco::Track::recHitsBegin(), reco::Track::recHitsEnd(), and setupAddElement().

Referenced by FWTrackHitsDetailView::addHits(), FWConvTrackHitsDetailView::addHits(), and FWTracksRecHitsProxyBuilder::build().

382 {
383  // master is true if the product is for proxy builder
384  const FWGeometry *geom = iItem->getGeom();
385 
386  const edmNew::DetSetVector<SiStripCluster> * allClusters = 0;
387  if( addNearbyClusters )
388  {
389  for( trackingRecHit_iterator it = t.recHitsBegin(), itEnd = t.recHitsEnd(); it != itEnd; ++it )
390  {
391  if( typeid( **it ) == typeid( SiStripRecHit2D ))
392  {
393  const SiStripRecHit2D &hit = static_cast<const SiStripRecHit2D &>( **it );
394  if( hit.cluster().isNonnull() && hit.cluster().isAvailable())
395  {
396  allClusters = hit.cluster().product();
397  break;
398  }
399  }
400  else if( typeid( **it ) == typeid( SiStripRecHit1D ))
401  {
402  const SiStripRecHit1D &hit = static_cast<const SiStripRecHit1D &>( **it );
403  if( hit.cluster().isNonnull() && hit.cluster().isAvailable())
404  {
405  allClusters = hit.cluster().product();
406  break;
407  }
408  }
409  }
410  }
411 
412  for( trackingRecHit_iterator it = t.recHitsBegin(), itEnd = t.recHitsEnd(); it != itEnd; ++it )
413  {
414  unsigned int rawid = (*it)->geographicalId();
415  if( ! geom->contains( rawid ))
416  {
417  fwLog( fwlog::kError )
418  << "failed to get geometry of SiStripCluster with detid: "
419  << rawid << std::endl;
420 
421  continue;
422  }
423 
424  const float* pars = geom->getParameters( rawid );
425 
426  // -- get phi from SiStripHit
427  TrackingRecHitRef rechitRef = *it;
428  const TrackingRecHit *rechit = &( *rechitRef );
429  const SiStripCluster *cluster = extractClusterFromTrackingRecHit( rechit );
430 
431  if( cluster )
432  {
433  if( allClusters != 0 )
434  {
435  const edmNew::DetSet<SiStripCluster> & clustersOnThisDet = (*allClusters)[rechit->geographicalId().rawId()];
436 
437  for( edmNew::DetSet<SiStripCluster>::const_iterator itc = clustersOnThisDet.begin(), edc = clustersOnThisDet.end(); itc != edc; ++itc )
438  {
439 
440  TEveStraightLineSet *scposition = new TEveStraightLineSet;
441  scposition->SetDepthTest( false );
442  scposition->SetPickable( kTRUE );
443 
444  short firststrip = itc->firstStrip();
445 
446  if( &*itc == cluster )
447  {
448  scposition->SetTitle( Form( "Exact SiStripCluster from TrackingRecHit, first strip %d", firststrip ));
449  scposition->SetLineColor( kGreen );
450  }
451  else
452  {
453  scposition->SetTitle( Form( "SiStripCluster, first strip %d", firststrip ));
454  scposition->SetLineColor( kRed );
455  }
456 
457  float localTop[3] = { 0.0, 0.0, 0.0 };
458  float localBottom[3] = { 0.0, 0.0, 0.0 };
459 
460  fireworks::localSiStrip( firststrip, localTop, localBottom, pars, rawid );
461 
462  float globalTop[3];
463  float globalBottom[3];
464  geom->localToGlobal( rawid, localTop, globalTop, localBottom, globalBottom );
465 
466  scposition->AddLine( globalTop[0], globalTop[1], globalTop[2],
467  globalBottom[0], globalBottom[1], globalBottom[2] );
468 
469  setupAddElement( scposition, tList, iItem, master, false );
470  }
471  }
472  else
473  {
474  short firststrip = cluster->firstStrip();
475  TEveStraightLineSet *scposition = new TEveStraightLineSet;
476  scposition->SetDepthTest( false );
477  scposition->SetPickable( kTRUE );
478  scposition->SetTitle( Form( "SiStripCluster, first strip %d", firststrip ));
479 
480  float localTop[3] = { 0.0, 0.0, 0.0 };
481  float localBottom[3] = { 0.0, 0.0, 0.0 };
482 
483  fireworks::localSiStrip( firststrip, localTop, localBottom, pars, rawid );
484 
485  float globalTop[3];
486  float globalBottom[3];
487  geom->localToGlobal( rawid, localTop, globalTop, localBottom, globalBottom );
488 
489  scposition->AddLine( globalTop[0], globalTop[1], globalTop[2],
490  globalBottom[0], globalBottom[1], globalBottom[2] );
491 
492  setupAddElement( scposition, tList, iItem, master, true );
493  }
494  }
495  else if( !rechit->isValid() && ( rawid != 0 )) // lost hit
496  {
497  if( allClusters != 0 )
498  {
499  edmNew::DetSetVector<SiStripCluster>::const_iterator itds = allClusters->find( rawid );
500  if( itds != allClusters->end())
501  {
502  const edmNew::DetSet<SiStripCluster> & clustersOnThisDet = *itds;
503  for( edmNew::DetSet<SiStripCluster>::const_iterator itc = clustersOnThisDet.begin(), edc = clustersOnThisDet.end(); itc != edc; ++itc )
504  {
505  short firststrip = itc->firstStrip();
506 
507  TEveStraightLineSet *scposition = new TEveStraightLineSet;
508  scposition->SetDepthTest( false );
509  scposition->SetPickable( kTRUE );
510  scposition->SetTitle( Form( "Lost SiStripCluster, first strip %d", firststrip ));
511 
512  float localTop[3] = { 0.0, 0.0, 0.0 };
513  float localBottom[3] = { 0.0, 0.0, 0.0 };
514 
515  fireworks::localSiStrip( firststrip, localTop, localBottom, pars, rawid );
516 
517  float globalTop[3];
518  float globalBottom[3];
519  geom->localToGlobal( rawid, localTop, globalTop, localBottom, globalBottom );
520 
521  scposition->AddLine( globalTop[0], globalTop[1], globalTop[2],
522  globalBottom[0], globalBottom[1], globalBottom[2] );
523 
524 
525  setupAddElement( scposition, tList, iItem, master, false );
526  scposition->SetLineColor( kRed );
527  }
528  }
529  }
530  }
531  else
532  {
533  fwLog( fwlog::kDebug )
534  << "*ANOTHER* option possible: valid=" << rechit->isValid()
535  << ", rawid=" << rawid << std::endl;
536  }
537  }
538 }
SiStripRecHit1D::cluster
ClusterRef cluster() const
Definition: SiStripRecHit1D.h:28
edmNew::DetSetVector::const_iterator
boost::transform_iterator< IterHelp, const_IdIter > const_iterator
Definition: DetSetVectorNew.h:118
edmNew::DetSetVector::end
const_iterator end(bool update=false) const
Definition: DetSetVectorNew.h:331
FWGeometry::getParameters
const float * getParameters(unsigned int id) const
Definition: FWGeometry.cc:293
SiStripCluster::firstStrip
uint16_t firstStrip() const
Definition: SiStripCluster.h:35
edmNew::DetSet::const_iterator
data_type const * const_iterator
Definition: DetSetNew.h:30
edm::Ref::isAvailable
bool isAvailable() const
Definition: Ref.h:276
fireworks::localSiStrip
void localSiStrip(short strip, float *localTop, float *localBottom, const float *pars, unsigned int id)
Definition: TrackUtils.cc:291
DetId::rawId
uint32_t rawId() const
get the raw id
Definition: DetId.h:43
edm::Ref::isNonnull
bool isNonnull() const
Checks for non-null.
Definition: Ref.h:250
autoMagneticFieldProducer_cfi.master
tuple master
Definition: autoMagneticFieldProducer_cfi.py:99
fireworks::extractClusterFromTrackingRecHit
const SiStripCluster * extractClusterFromTrackingRecHit(const TrackingRecHit *rh)
Definition: TrackUtils.cc:358
FWGeometry::localToGlobal
void localToGlobal(unsigned int id, const float *local, float *global, bool translatep=true) const
Definition: FWGeometry.cc:325
SiStripRecHit2D::cluster
ClusterRef cluster() const
Definition: SiStripRecHit2D.h:29
FWGeometry::contains
bool contains(unsigned int id) const
Definition: FWGeometry.h:109
reco::Track::recHitsBegin
trackingRecHit_iterator recHitsBegin() const
Iterator to first hit on the track.
Definition: Track.h:63
fwLog
#define fwLog(_level_)
Definition: fwLog.h:50
TrackingRecHit::isValid
bool isValid() const
Definition: TrackingRecHit.h:143
edmNew::DetSetVector::find
const_iterator find(id_type i, bool update=false) const
Definition: DetSetVectorNew.h:310
edmNew::DetSet::end
iterator end()
Definition: DetSetNew.h:70
TrackingRecHit::geographicalId
DetId geographicalId() const
Definition: TrackingRecHit.h:120
edm::Ref::product
C const * product() const
Accessor for product collection.
Definition: Ref.h:601
fireworks::setupAddElement
void setupAddElement(TEveElement *el, TEveElement *parent, const FWEventItem *item, bool master, bool color)
Definition: TrackUtils.cc:336
FWEventItem::getGeom
const FWGeometry * getGeom() const
Definition: FWEventItem.cc:681
reco::Track::recHitsEnd
trackingRecHit_iterator recHitsEnd() const
Iterator to last hit on the track.
Definition: Track.h:65
edmNew::DetSet::begin
iterator begin()
Definition: DetSetNew.h:67
void fireworks::addStraightLineSegment ( TEveStraightLineSet *  marker,
reco::Candidate const *  cand,
double  scale_factor = 2 
)

Definition at line 23 of file CandidateUtils.cc.

References funct::cos(), phi, reco::Candidate::phi(), reco::Candidate::pt(), funct::sin(), findQualityFiles::size, reco::Candidate::theta(), and theta().

Referenced by FWMuonGlimpseProxyBuilder::build(), and FWElectronGlimpseProxyBuilder::build().

26 {
27  double phi = cand->phi();
28  double theta = cand->theta();
29  double size = cand->pt() * scale_factor;
30  marker->AddLine( 0, 0, 0, size * cos(phi)*sin(theta), size *sin(phi)*sin(theta), size*cos(theta));
31 }
funct::sin
Sin< T >::type sin(const T &t)
Definition: Sin.h:22
theta
Geom::Theta< T > theta() const
Definition: Basic3DVectorLD.h:170
funct::cos
Cos< T >::type cos(const T &t)
Definition: Cos.h:22
findQualityFiles.size
tuple size
Write out results.
Definition: findQualityFiles.py:442
phi
Definition: DDAxes.h:10
const TGGC & fireworks::boldGC ( )

Definition at line 12 of file GlobalContexts.cc.

References fetchall_from_DQM_v2::pool.

Referenced by FWPSetTableManager::cellRenderer().

13 {
14  static TGGC s_boldGC(*gClient->GetResourcePool()->GetFrameGC());
15 
16  TGFontPool *pool = gClient->GetFontPool();
17  //TGFont *font = pool->FindFontByHandle(s_boldGC.GetFont());
18  //FontAttributes_t attributes = font->GetFontAttributes();
19 
20  /*
21  This doesn't seem to work:
22  attributes.fWeight = 1;
23  TGFont *newFont = pool->GetFont(attributes.fFamily, 9,
24  attributes.fWeight, attributes.fSlant);
25 
26  But this does:
27  */
28 
29  TGFont* newFont = pool->GetFont("-*-helvetica-bold-r-*-*-12-*-*-*-*-*-iso8859-1");
30 
31  if ( ! newFont )
32  return s_boldGC;
33 
34  s_boldGC.SetFont(newFont->GetFontHandle());
35 
36  return s_boldGC;
37 }
fetchall_from_DQM_v2.pool
tuple pool
Definition: fetchall_from_DQM_v2.py:139
void fireworks::createSegment ( int  detector,
bool  matchedSegment,
float  segmentLength,
float  segmentLimit,
float *  segmentPosition,
float *  segmentDirection,
float *  segmentInnerPoint,
float *  segmentOuterPoint 
)

Definition at line 10 of file SegmentUtils.cc.

References funct::cos(), MuonSubdetId::CSC, MuonSubdetId::DT, fwLog, fwlog::kWarning, mag(), mathSSE::sqrt(), funct::tan(), and theta().

18  {
19  if( detector == MuonSubdetId::CSC )
20  {
21  if( matchedSegment )
22  {
23  segmentOuterPoint[0] = segmentPosition[0] + segmentLength*segmentDirection[0];
24  segmentOuterPoint[1] = segmentPosition[1] + segmentLength*segmentDirection[1];
25  segmentOuterPoint[2] = segmentLength;
26 
27  if( fabs(segmentOuterPoint[1]) > segmentLimit )
28  segmentOuterPoint[1] = segmentLimit * ( segmentOuterPoint[1] / fabs( segmentOuterPoint[1] ));
29 
30  segmentInnerPoint[0] = segmentPosition[0] - segmentLength*segmentDirection[0];
31  segmentInnerPoint[1] = segmentPosition[1] - segmentLength*segmentDirection[1];
32  segmentInnerPoint[2] = -segmentLength;
33 
34  if( fabs(segmentInnerPoint[1]) > segmentLimit )
35  segmentInnerPoint[1] = segmentLimit * ( segmentInnerPoint[1] / fabs( segmentInnerPoint[1] ));
36 
37  return;
38  }
39  else
40  {
41  segmentOuterPoint[0] = segmentPosition[0] + segmentDirection[0] * ( segmentPosition[2] / segmentDirection[2] );
42  segmentOuterPoint[1] = segmentPosition[1] + segmentDirection[1] * ( segmentPosition[2] / segmentDirection[2] );
43  segmentOuterPoint[2] = segmentPosition[2];
44 
45  if( fabs( segmentOuterPoint[1] ) > segmentLength )
46  segmentOuterPoint[1] = segmentLength * ( segmentOuterPoint[1]/fabs( segmentOuterPoint[1] ));
47 
48  segmentInnerPoint[0] = segmentPosition[0] - segmentDirection[0] * ( segmentPosition[2] / segmentDirection[2] );
49  segmentInnerPoint[1] = segmentPosition[1] - segmentDirection[1] * ( segmentPosition[2] / segmentDirection[2] );
50  segmentInnerPoint[2] = -segmentPosition[2];
51 
52  if( fabs( segmentInnerPoint[1] ) > segmentLength )
53  segmentInnerPoint[1] = segmentLength * ( segmentInnerPoint[1] / fabs( segmentInnerPoint[1] ));
54 
55  return;
56  }
57  }
58 
59  if ( detector == MuonSubdetId::DT )
60  {
61  if( matchedSegment )
62  {
63  segmentOuterPoint[0] = segmentPosition[0] + segmentLength*segmentDirection[0];
64  segmentOuterPoint[1] = segmentPosition[1] + segmentLength*segmentDirection[1];
65  segmentOuterPoint[2] = segmentLength;
66 
67  segmentInnerPoint[0] = segmentPosition[0] - segmentLength*segmentDirection[0];
68  segmentInnerPoint[1] = segmentPosition[1] - segmentLength*segmentDirection[1];
69  segmentInnerPoint[2] = -segmentLength;
70 
71  return;
72  }
73  else
74  {
75  double mag = sqrt( segmentDirection[0] * segmentDirection[0] +
76  segmentDirection[1] * segmentDirection[1] +
77  segmentDirection[2] * segmentDirection[2] );
78 
79  double theta = atan2( sqrt( segmentDirection[0] * segmentDirection[0]
80  + segmentDirection[1] * segmentDirection[1] ), segmentDirection[2] );
81 
82  double newSegmentLength = segmentLength / cos( theta );
83 
84  segmentInnerPoint[0] = segmentPosition[0] + ( segmentDirection[0] / mag ) * newSegmentLength;
85  segmentInnerPoint[1] = segmentPosition[1] + ( segmentDirection[1] / mag ) * newSegmentLength;
86  segmentInnerPoint[2] = segmentPosition[2] + ( segmentDirection[2] / mag ) * newSegmentLength;
87 
88  segmentOuterPoint[0] = segmentPosition[0] - ( segmentDirection[0] / mag ) * newSegmentLength;
89  segmentOuterPoint[1] = segmentPosition[1] - ( segmentDirection[1] / mag ) * newSegmentLength;
90  segmentOuterPoint[2] = segmentPosition[2] - ( segmentDirection[2] / mag ) * newSegmentLength;
91 
92  if( fabs(segmentOuterPoint[0]) > segmentLimit )
93  {
94  segmentOuterPoint[0] = segmentLimit * ( segmentOuterPoint[0]/fabs( segmentOuterPoint[0] ));
95  segmentOuterPoint[1] = ( segmentOuterPoint[1]/fabs( segmentOuterPoint[1] )) * tan( theta );
96  }
97 
98  return;
99  }
100  }
101 
102  fwLog( fwlog::kWarning ) << "MuonSubdetId: " << detector << std::endl;
103  return;
104  }
mag
T mag() const
The vector magnitude. Equivalent to sqrt(vec.mag2())
Definition: Basic3DVectorLD.h:147
theta
Geom::Theta< T > theta() const
Definition: Basic3DVectorLD.h:170
MuonSubdetId::CSC
static const int CSC
Definition: MuonSubdetId.h:13
mathSSE::sqrt
T sqrt(T t)
Definition: SSEVec.h:48
funct::cos
Cos< T >::type cos(const T &t)
Definition: Cos.h:22
funct::tan
Tan< T >::type tan(const T &t)
Definition: Tan.h:22
fwLog
#define fwLog(_level_)
Definition: fwLog.h:50
MuonSubdetId::DT
static const int DT
Definition: MuonSubdetId.h:12
const TString fireworks::datadir ( "/src/Fireworks/Core/"  )
void fireworks::drawEnergyScaledBox3D ( const float *  corners,
float  scale,
TEveElement *  comp,
FWProxyBuilderBase *  pb,
bool  invert = false 
)

Definition at line 239 of file BuilderUtils.cc.

References addBox(), and energyScaledBox3DCorners().

240  {
241  std::vector<float> scaledCorners( 24 );
242  energyScaledBox3DCorners(corners, scale, scaledCorners, invert);
243  addBox( scaledCorners, comp, pb );
244  }
AlCaHLTBitMon_QueryRunRegistry.comp
string comp
Definition: AlCaHLTBitMon_QueryRunRegistry.py:248
pileupReCalc_HLTpaths.scale
int scale
Definition: pileupReCalc_HLTpaths.py:111
fireworks::energyScaledBox3DCorners
void energyScaledBox3DCorners(const float *corners, float scale, std::vector< float > &, bool invert=false)
Definition: BuilderUtils.cc:213
fireworks::addBox
void addBox(const std::vector< float > &corners, TEveElement *, FWProxyBuilderBase *)
Definition: BuilderUtils.cc:146
void fireworks::drawEnergyTower3D ( const float *  corners,
float  scale,
TEveElement *  comp,
FWProxyBuilderBase *  pb,
bool  reflect = false 
)

Definition at line 323 of file BuilderUtils.cc.

References addBox(), and energyTower3DCorners().

324  {
325  std::vector<float> scaledCorners( 24 );
326  energyTower3DCorners(corners, scale, scaledCorners, reflect);
327  addBox( scaledCorners, comp, pb );
328  }
AlCaHLTBitMon_QueryRunRegistry.comp
string comp
Definition: AlCaHLTBitMon_QueryRunRegistry.py:248
pileupReCalc_HLTpaths.scale
int scale
Definition: pileupReCalc_HLTpaths.py:111
fireworks::energyTower3DCorners
void energyTower3DCorners(const float *corners, float scale, std::vector< float > &, bool reflect=false)
Definition: BuilderUtils.cc:284
fireworks::addBox
void addBox(const std::vector< float > &corners, TEveElement *, FWProxyBuilderBase *)
Definition: BuilderUtils.cc:146
void fireworks::drawEtScaledBox3D ( const float *  corners,
float  energy,
float  maxEnergy,
TEveElement *  comp,
FWProxyBuilderBase *  pb,
bool  reflect = false 
)

Definition at line 276 of file BuilderUtils.cc.

References addBox(), and etScaledBox3DCorners().

277  {
278  std::vector<float> scaledCorners( 24 );
279  etScaledBox3DCorners(corners, energy, maxEnergy, scaledCorners, invert);
280  addBox( scaledCorners, comp, pb );
281  }
AlCaHLTBitMon_QueryRunRegistry.comp
string comp
Definition: AlCaHLTBitMon_QueryRunRegistry.py:248
fireworks::addBox
void addBox(const std::vector< float > &corners, TEveElement *, FWProxyBuilderBase *)
Definition: BuilderUtils.cc:146
fireworks::etScaledBox3DCorners
void etScaledBox3DCorners(const float *corners, float energy, float maxEnergy, std::vector< float > &scaledCorners, bool reflect=false)
Definition: BuilderUtils.cc:247
relval_parameters_module.energy
string energy
Definition: relval_parameters_module.py:29
void fireworks::drawEtTower3D ( const float *  corners,
float  scale,
TEveElement *  comp,
FWProxyBuilderBase *  pb,
bool  reflect = false 
)

Definition at line 372 of file BuilderUtils.cc.

References addBox(), and etTower3DCorners().

373  {
374  std::vector<float> scaledCorners( 24 );
375  etTower3DCorners(corners, scale, scaledCorners, reflect);
376  addBox( scaledCorners, comp, pb );
377  }
fireworks::etTower3DCorners
void etTower3DCorners(const float *corners, float scale, std::vector< float > &, bool reflect=false)
Definition: BuilderUtils.cc:332
AlCaHLTBitMon_QueryRunRegistry.comp
string comp
Definition: AlCaHLTBitMon_QueryRunRegistry.py:248
pileupReCalc_HLTpaths.scale
int scale
Definition: pileupReCalc_HLTpaths.py:111
fireworks::addBox
void addBox(const std::vector< float > &corners, TEveElement *, FWProxyBuilderBase *)
Definition: BuilderUtils.cc:146
void fireworks::energyScaledBox3DCorners ( const float *  corners,
float  scale,
std::vector< float > &  scaledCorners,
bool  invert = false 
)

Definition at line 213 of file BuilderUtils.cc.

References f, i, and invertBox().

Referenced by drawEnergyScaledBox3D().

214  {
215  std::vector<float> centre( 3, 0 );
216 
217  for( unsigned int i = 0; i < 24; i += 3 )
218  {
219  centre[0] += corners[i];
220  centre[1] += corners[i + 1];
221  centre[2] += corners[i + 2];
222  }
223 
224  for( unsigned int i = 0; i < 3; ++i )
225  centre[i] *= 1.0f / 8.0f;
226 
227  // Coordinates for a scaled version of the original box
228  for( unsigned int i = 0; i < 24; i += 3 )
229  {
230  scaledCorners[i] = centre[0] + ( corners[i] - centre[0] ) * scale;
231  scaledCorners[i + 1] = centre[1] + ( corners[i + 1] - centre[1] ) * scale;
232  scaledCorners[i + 2] = centre[2] + ( corners[i + 2] - centre[2] ) * scale;
233  }
234 
235  if( invert )
236  invertBox( scaledCorners );
237  }
i
int i
Definition: DBlmapReader.cc:9
pileupReCalc_HLTpaths.scale
int scale
Definition: pileupReCalc_HLTpaths.py:111
fireworks::invertBox
void invertBox(std::vector< float > &corners)
Definition: BuilderUtils.cc:127
f
double f[11][100]
Definition: MuScleFitUtils.cc:78
void fireworks::energyTower3DCorners ( const float *  corners,
float  scale,
std::vector< float > &  scaledCorners,
bool  reflect = false 
)

Definition at line 284 of file BuilderUtils.cc.

References diffTreeTool::diff, i, and pileupReCalc_HLTpaths::scale.

Referenced by FWPRCaloTowerProxyBuilder::build(), FWCaloClusterProxyBuilder::build(), FWPCaloHitProxyBuilder::build(), and drawEnergyTower3D().

285  {
286  for( int i = 0; i < 24; ++i )
287  scaledCorners[i] = corners[i];
288  // Coordinates of a front face scaled
289  if( reflect )
290  {
291  // We know, that an ES rechit geometry in -Z needs correction.
292  // The back face is actually its front face.
293  for( unsigned int i = 0; i < 12; i += 3 )
294  {
295  TEveVector diff( corners[i] - corners[i + 12], corners[i + 1] - corners[i + 13], corners[i + 2] - corners[i + 14] );
296  diff.Normalize();
297  diff *= scale;
298 
299  scaledCorners[i] = corners[i] + diff.fX;
300  scaledCorners[i + 1] = corners[i + 1] + diff.fY;
301  scaledCorners[i + 2] = corners[i + 2] + diff.fZ;
302  }
303  }
304  else
305  {
306  for( unsigned int i = 0; i < 12; i += 3 )
307  {
308  TEveVector diff( corners[i + 12] - corners[i], corners[i + 13] - corners[i + 1], corners[i + 14] - corners[i + 2] );
309  diff.Normalize();
310  diff *= scale;
311 
312  scaledCorners[i] = corners[i + 12];
313  scaledCorners[i + 1] = corners[i + 13];
314  scaledCorners[i + 2] = corners[i + 14];
315 
316  scaledCorners[i + 12] = corners[i + 12] + diff.fX;
317  scaledCorners[i + 13] = corners[i + 13] + diff.fY;
318  scaledCorners[i + 14] = corners[i + 14] + diff.fZ;
319  }
320  }
321  }
i
int i
Definition: DBlmapReader.cc:9
pileupReCalc_HLTpaths.scale
int scale
Definition: pileupReCalc_HLTpaths.py:111
void fireworks::etScaledBox3DCorners ( const float *  corners,
float  energy,
float  maxEnergy,
std::vector< float > &  scaledCorners,
bool  reflect = false 
)

Definition at line 247 of file BuilderUtils.cc.

References trackerHits::c, f, i, invertBox(), pileupReCalc_HLTpaths::scale, and funct::sin().

Referenced by drawEtScaledBox3D().

248  {
249  std::vector<float> centre( 3, 0 );
250 
251  for( unsigned int i = 0; i < 24; i += 3 )
252  {
253  centre[0] += corners[i];
254  centre[1] += corners[i + 1];
255  centre[2] += corners[i + 2];
256  }
257 
258  for( unsigned int i = 0; i < 3; ++i )
259  centre[i] *= 1.0f / 8.0f;
260 
261  TEveVector c( centre[0], centre[1], centre[2] );
262  float scale = energy / maxEnergy * sin( c.Theta());
263 
264  // Coordinates for a scaled version of the original box
265  for( unsigned int i = 0; i < 24; i += 3 )
266  {
267  scaledCorners[i] = centre[0] + ( corners[i] - centre[0] ) * scale;
268  scaledCorners[i + 1] = centre[1] + ( corners[i + 1] - centre[1] ) * scale;
269  scaledCorners[i + 2] = centre[2] + ( corners[i + 2] - centre[2] ) * scale;
270  }
271 
272  if( invert )
273  invertBox( scaledCorners );
274  }
i
int i
Definition: DBlmapReader.cc:9
funct::sin
Sin< T >::type sin(const T &t)
Definition: Sin.h:22
pileupReCalc_HLTpaths.scale
int scale
Definition: pileupReCalc_HLTpaths.py:111
fireworks::invertBox
void invertBox(std::vector< float > &corners)
Definition: BuilderUtils.cc:127
f
double f[11][100]
Definition: MuScleFitUtils.cc:78
relval_parameters_module.energy
string energy
Definition: relval_parameters_module.py:29
void fireworks::etTower3DCorners ( const float *  corners,
float  scale,
std::vector< float > &  scaledCorners,
bool  reflect = false 
)

Definition at line 332 of file BuilderUtils.cc.

References diffTreeTool::diff, i, and funct::sin().

Referenced by drawEtTower3D().

333  {
334  for( int i = 0; i < 24; ++i )
335  scaledCorners[i] = corners[i];
336  // Coordinates of a front face scaled
337  if( reflect )
338  {
339  // We know, that an ES rechit geometry in -Z needs correction.
340  // The back face is actually its front face.
341  for( unsigned int i = 0; i < 12; i += 3 )
342  {
343  TEveVector diff( corners[i] - corners[i + 12], corners[i + 1] - corners[i + 13], corners[i + 2] - corners[i + 14] );
344  diff.Normalize();
345  diff *= ( scale * sin( diff.Theta()));
346 
347  scaledCorners[i] = corners[i] + diff.fX;
348  scaledCorners[i + 1] = corners[i + 1] + diff.fY;
349  scaledCorners[i + 2] = corners[i + 2] + diff.fZ;
350  }
351  }
352  else
353  {
354  for( unsigned int i = 0; i < 12; i += 3 )
355  {
356  TEveVector diff( corners[i + 12] - corners[i], corners[i + 13] - corners[i + 1], corners[i + 14] - corners[i + 2] );
357  diff.Normalize();
358  diff *= ( scale * sin( diff.Theta()));
359 
360  scaledCorners[i] = corners[i + 12];
361  scaledCorners[i + 1] = corners[i + 13];
362  scaledCorners[i + 2] = corners[i + 14];
363 
364  scaledCorners[i + 12] = corners[i + 12] + diff.fX;
365  scaledCorners[i + 13] = corners[i + 13] + diff.fY;
366  scaledCorners[i + 14] = corners[i + 14] + diff.fZ;
367  }
368  }
369  }
i
int i
Definition: DBlmapReader.cc:9
funct::sin
Sin< T >::type sin(const T &t)
Definition: Sin.h:22
pileupReCalc_HLTpaths.scale
int scale
Definition: pileupReCalc_HLTpaths.py:111
const SiStripCluster * fireworks::extractClusterFromTrackingRecHit ( const TrackingRecHit *  rh)

Definition at line 358 of file TrackUtils.cc.

References fwLog, and fwlog::kDebug.

Referenced by addSiStripClusters().

359 {
360  const SiStripCluster* cluster = 0;
361 
362  if( const SiStripRecHit2D* hit2D = dynamic_cast<const SiStripRecHit2D*>( rechit ))
363  {
364  fwLog( fwlog::kDebug ) << "hit 2D ";
365 
366  cluster = hit2D->cluster().get();
367  }
368  if( cluster == 0 )
369  {
370  if( const SiStripRecHit1D* hit1D = dynamic_cast<const SiStripRecHit1D*>( rechit ))
371  {
372  fwLog( fwlog::kDebug ) << "hit 1D ";
373 
374  cluster = hit1D->cluster().get();
375  }
376  }
377  return cluster;
378 }
fwLog
#define fwLog(_level_)
Definition: fwLog.h:50
void fireworks::getDecomposedVersion ( const TString &  s,
int *  out 
)

Definition at line 30 of file fwPaths.cc.

Referenced by acceptDataFormatsVersion(), FWFileEntry::openFile(), and supportedDataFormatsVersion().

31 {
32  TPMERegexp re("CMSSW_(\\d+)_(\\d+)_");
33  re.Match(s);
34  if (re.NMatches() > 2)
35  {
36  out[0] = re[1].Atoi();
37  out[1] = re[2].Atoi();
38  }
39 }
dbtoconf.out
tuple out
Definition: dbtoconf.py:99
std::string fireworks::getLocalTime ( const edm::EventBase &  event)

Definition at line 113 of file BuilderUtils.cc.

References AlCaHLTBitMon_QueryRunRegistry::string, edmStreamStallGrapher::t, runonSM::text, edm::EventBase::time(), and edm::Timestamp::value().

Referenced by CmsShowMainFrame::loadEvent(), and FWEventAnnotation::updateOverlayText().

114  {
115  time_t t( event.time().value() >> 32 );
116  std::string text( asctime( localtime( &t ) ) );
117  size_t pos = text.find('\n');
118  if( pos != std::string::npos ) text = text.substr( 0, pos );
119  text += " ";
120  if( daylight )
121  text += tzname[1];
122  else
123  text += tzname[0];
124  return text;
125  }
AlCaHLTBitMon_QueryRunRegistry.string
string string
Definition: AlCaHLTBitMon_QueryRunRegistry.py:255
runonSM.text
tuple text
Definition: runonSM.py:42
edm::Timestamp::value
TimeValue_t value() const
Definition: Timestamp.h:56
edm::EventBase::time
edm::Timestamp time() const
Definition: EventBase.h:57
std::pair< double, double > fireworks::getPhiRange ( const std::vector< double > &  phis,
double  phi 
)

Definition at line 21 of file BuilderUtils.cc.

References i, M_PI, bookConverter::max, and min().

Referenced by FWTauProxyBuilderBase::buildBaseTau(), and makeRhoPhiSuperCluster().

22  {
23  double min = 100;
24  double max = -100;
25 
26  for( std::vector<double>::const_iterator i = phis.begin();
27  i != phis.end(); ++i )
28  {
29  double aphi = *i;
30  // make phi continuous around jet phi
31  if( aphi - phi > M_PI ) aphi -= 2*M_PI;
32  if( phi - aphi > M_PI ) aphi += 2*M_PI;
33  if( aphi > max ) max = aphi;
34  if( aphi < min ) min = aphi;
35  }
36 
37  if( min > max ) return std::pair<double,double>( 0, 0 );
38 
39  return std::pair<double,double>( min, max );
40  }
i
int i
Definition: DBlmapReader.cc:9
min
T min(T a, T b)
Definition: MathUtil.h:58
M_PI
#define M_PI
Definition: readTestVector.cc:25
phi
Definition: DDAxes.h:10
TEveGeoShape * fireworks::getShape ( const char *  name,
TGeoBBox *  shape,
Color_t  color 
)

Definition at line 42 of file BuilderUtils.cc.

References trackerHits::c.

Referenced by FWTauProxyBuilderBase::buildBaseTau(), FWMETProxyBuilder::buildViewType(), FWGlimpseView::FWGlimpseView(), and makeRhoPhiSuperCluster().

45  {
46  TEveGeoShape* egs = new TEveGeoShape( name );
47  TColor* c = gROOT->GetColor( color );
48  Float_t rgba[4] = { 1, 0, 0, 1 };
49  if( c )
50  {
51  rgba[0] = c->GetRed();
52  rgba[1] = c->GetGreen();
53  rgba[2] = c->GetBlue();
54  }
55  egs->SetMainColorRGB( rgba[0], rgba[1], rgba[2] );
56  egs->SetShape( shape );
57  return egs;
58  }
create_public_lumi_plots.color
list color
Definition: create_public_lumi_plots.py:1800
std::string fireworks::getTimeGMT ( const edm::EventBase &  event)

Definition at line 103 of file BuilderUtils.cc.

References AlCaHLTBitMon_QueryRunRegistry::string, edmStreamStallGrapher::t, runonSM::text, edm::EventBase::time(), and edm::Timestamp::value().

104  {
105  time_t t( event.time().value() >> 32 );
106  std::string text( asctime( gmtime( &t ) ) );
107  size_t pos = text.find( '\n' );
108  if( pos != std::string::npos ) text = text.substr( 0, pos );
109  text += " GMT";
110  return text;
111  }
AlCaHLTBitMon_QueryRunRegistry.string
string string
Definition: AlCaHLTBitMon_QueryRunRegistry.py:255
runonSM.text
tuple text
Definition: runonSM.py:42
edm::Timestamp::value
TimeValue_t value() const
Definition: Timestamp.h:56
edm::EventBase::time
edm::Timestamp time() const
Definition: EventBase.h:57
const TGGC & fireworks::greenGC ( )

Definition at line 39 of file GlobalContexts.cc.

40 {
41  static TGGC s_greenGC(*gClient->GetResourcePool()->GetFrameGC());
42  s_greenGC.SetForeground(gVirtualX->GetPixel(kGreen-5));
43  return s_greenGC;
44 }
std::string fireworks::info ( const DetId &  id)

Definition at line 680 of file TrackUtils.cc.

References DetId::Calo, CSCDetId::chamber(), MuonSubdetId::CSC, HcalDetId::depth(), MuonSubdetId::DT, DetId::Ecal, EcalBarrel, EcalEndcap, EcalLaserPnDiode, EcalPreshower, EcalTriggerTower, CSCDetId::endcap(), DetId::Hcal, HcalBarrel, HcalEmpty, HcalEndcap, HcalForward, HcalOther, HcalOuter, HcalTriggerTower, EEDetId::ic(), CaloTowerDetId::ieta(), HcalDetId::ieta(), EBDetId::ieta(), CaloTowerDetId::iphi(), HcalDetId::iphi(), EBDetId::iphi(), EEDetId::iquadrant(), EEDetId::isc(), EEDetId::ix(), EEDetId::iy(), PXBDetId::layer(), TOBDetId::layer(), TIBDetId::layer(), CSCDetId::layer(), RPCDetId::layer(), DetId::Muon, PixelSubdetector::PixelBarrel, PixelSubdetector::PixelEndcap, RPCDetId::region(), RPCDetId::ring(), CSCDetId::ring(), RPCDetId::roll(), MuonSubdetId::RPC, DTChamberId::sector(), RPCDetId::sector(), DTChamberId::station(), RPCDetId::station(), CSCDetId::station(), HcalDetId::subdet(), RPCDetId::subsector(), StripSubdetector::TEC, StripSubdetector::TIB, StripSubdetector::TID, StripSubdetector::TOB, EBDetId::tower_ieta(), EBDetId::tower_iphi(), DetId::Tracker, DTChamberId::wheel(), TIDDetId::wheel(), and TECDetId::wheel().

Referenced by FWTrackHitsDetailView::addModules(), FWConvTrackHitsDetailView::addModules(), FWTracksModulesProxyBuilder::build(), and info().

680  {
681  std::ostringstream oss;
682 
683  oss << "DetId: " << id.rawId() << "\n";
684 
685  switch ( id.det() ) {
686 
687  case DetId::Tracker:
688  switch ( id.subdetId() ) {
689  case StripSubdetector::TIB:
690  {
691  oss <<"TIB "<<TIBDetId(id).layer();
692  }
693  break;
694  case StripSubdetector::TOB:
695  {
696  oss <<"TOB "<<TOBDetId(id).layer();
697  }
698  break;
699  case StripSubdetector::TEC:
700  {
701  oss <<"TEC "<<TECDetId(id).wheel();
702  }
703  break;
704  case StripSubdetector::TID:
705  {
706  oss <<"TID "<<TIDDetId(id).wheel();
707  }
708  break;
709  case (int) PixelSubdetector::PixelBarrel:
710  {
711  oss <<"PixBarrel "<< PXBDetId(id).layer();
712  }
713  break;
714  case (int) PixelSubdetector::PixelEndcap:
715  {
716  oss <<"PixEndcap "<< PXBDetId(id).layer();
717  }
718  break;
719  }
720  break;
721 
722  case DetId::Muon:
723  switch ( id.subdetId() ) {
724  case MuonSubdetId::DT:
725  {
726  DTChamberId detId(id.rawId());
727  oss << "DT chamber (wheel, station, sector): "
728  << detId.wheel() << ", "
729  << detId.station() << ", "
730  << detId.sector();
731  }
732  break;
733  case MuonSubdetId::CSC:
734  {
735  CSCDetId detId(id.rawId());
736  oss << "CSC chamber (endcap, station, ring, chamber, layer): "
737  << detId.endcap() << ", "
738  << detId.station() << ", "
739  << detId.ring() << ", "
740  << detId.chamber() << ", "
741  << detId.layer();
742  }
743  break;
744  case MuonSubdetId::RPC:
745  {
746  RPCDetId detId(id.rawId());
747  oss << "RPC chamber ";
748  switch ( detId.region() ) {
749  case 0:
750  oss << "/ barrel / (wheel, station, sector, layer, subsector, roll): "
751  << detId.ring() << ", "
752  << detId.station() << ", "
753  << detId.sector() << ", "
754  << detId.layer() << ", "
755  << detId.subsector() << ", "
756  << detId.roll();
757  break;
758  case 1:
759  oss << "/ forward endcap / (wheel, station, sector, layer, subsector, roll): "
760  << detId.ring() << ", "
761  << detId.station() << ", "
762  << detId.sector() << ", "
763  << detId.layer() << ", "
764  << detId.subsector() << ", "
765  << detId.roll();
766  break;
767  case -1:
768  oss << "/ backward endcap / (wheel, station, sector, layer, subsector, roll): "
769  << detId.ring() << ", "
770  << detId.station() << ", "
771  << detId.sector() << ", "
772  << detId.layer() << ", "
773  << detId.subsector() << ", "
774  << detId.roll();
775  break;
776  }
777  }
778  break;
779  }
780  break;
781 
782  case DetId::Calo:
783  {
784  CaloTowerDetId detId( id.rawId() );
785  oss << "CaloTower (ieta, iphi): "
786  << detId.ieta() << ", "
787  << detId.iphi();
788  }
789  break;
790 
791  case DetId::Ecal:
792  switch ( id.subdetId() ) {
793  case EcalBarrel:
794  {
795  EBDetId detId(id);
796  oss << "EcalBarrel (ieta, iphi, tower_ieta, tower_iphi): "
797  << detId.ieta() << ", "
798  << detId.iphi() << ", "
799  << detId.tower_ieta() << ", "
800  << detId.tower_iphi();
801  }
802  break;
803  case EcalEndcap:
804  {
805  EEDetId detId(id);
806  oss << "EcalEndcap (ix, iy, SuperCrystal, crystal, quadrant): "
807  << detId.ix() << ", "
808  << detId.iy() << ", "
809  << detId.isc() << ", "
810  << detId.ic() << ", "
811  << detId.iquadrant();
812  }
813  break;
814  case EcalPreshower:
815  oss << "EcalPreshower";
816  break;
817  case EcalTriggerTower:
818  oss << "EcalTriggerTower";
819  break;
820  case EcalLaserPnDiode:
821  oss << "EcalLaserPnDiode";
822  break;
823  }
824  break;
825 
826  case DetId::Hcal:
827  {
828  HcalDetId detId(id);
829  switch ( detId.subdet() ) {
830  case HcalEmpty:
831  oss << "HcalEmpty ";
832  break;
833  case HcalBarrel:
834  oss << "HcalBarrel ";
835  break;
836  case HcalEndcap:
837  oss << "HcalEndcap ";
838  break;
839  case HcalOuter:
840  oss << "HcalOuter ";
841  break;
842  case HcalForward:
843  oss << "HcalForward ";
844  break;
845  case HcalTriggerTower:
846  oss << "HcalTriggerTower ";
847  break;
848  case HcalOther:
849  oss << "HcalOther ";
850  break;
851  }
852  oss << "(ieta, iphi, depth):"
853  << detId.ieta() << ", "
854  << detId.iphi() << ", "
855  << detId.depth();
856  }
857  break;
858  default :;
859  }
860  return oss.str();
861 }
TOBDetId::layer
unsigned int layer() const
layer id
Definition: TOBDetId.h:39
PXBDetId::layer
unsigned int layer() const
layer id
Definition: PXBDetId.h:35
MuonSubdetId::CSC
static const int CSC
Definition: MuonSubdetId.h:13
TECDetId::wheel
unsigned int wheel() const
wheel id
Definition: TECDetId.h:52
TIBDetId::layer
unsigned int layer() const
layer id
Definition: TIBDetId.h:41
MuonSubdetId::RPC
static const int RPC
Definition: MuonSubdetId.h:14
MuonSubdetId::DT
static const int DT
Definition: MuonSubdetId.h:12
TIDDetId::wheel
unsigned int wheel() const
wheel id
Definition: TIDDetId.h:50
std::string fireworks::info ( const std::set< DetId > &  idSet)

Definition at line 864 of file TrackUtils.cc.

References info(), AlCaHLTBitMon_QueryRunRegistry::string, and runonSM::text.

864  {
865  std::string text;
866  for(std::set<DetId>::const_iterator id = idSet.begin(), idEnd = idSet.end(); id != idEnd; ++id)
867  {
868  text += info(*id);
869  text += "\n";
870  }
871  return text;
872 }
info
static const TGPicture * info(bool iBackgroundIsBlack)
Definition: FWCollectionSummaryWidget.cc:170
AlCaHLTBitMon_QueryRunRegistry.string
string string
Definition: AlCaHLTBitMon_QueryRunRegistry.py:255
runonSM.text
tuple text
Definition: runonSM.py:42
std::string fireworks::info ( const std::vector< DetId > &  idSet)

Definition at line 875 of file TrackUtils.cc.

References info(), AlCaHLTBitMon_QueryRunRegistry::string, and runonSM::text.

875  {
876  std::string text;
877  for(std::vector<DetId>::const_iterator id = idSet.begin(), idEnd = idSet.end(); id != idEnd; ++id)
878  {
879  text += info(*id);
880  text += "\n";
881  }
882  return text;
883 }
info
static const TGPicture * info(bool iBackgroundIsBlack)
Definition: FWCollectionSummaryWidget.cc:170
AlCaHLTBitMon_QueryRunRegistry.string
string string
Definition: AlCaHLTBitMon_QueryRunRegistry.py:255
runonSM.text
tuple text
Definition: runonSM.py:42
void fireworks::invertBox ( std::vector< float > &  corners)

Definition at line 127 of file BuilderUtils.cc.

References std::swap().

Referenced by energyScaledBox3DCorners(), etScaledBox3DCorners(), and FWCaloRecHitDigitSetProxyBuilder::viewContextBoxScale().

128  {
129  std::swap( corners[0], corners[9] );
130  std::swap( corners[1], corners[10] );
131  std::swap( corners[2], corners[11] );
132 
133  std::swap( corners[3], corners[6] );
134  std::swap( corners[4], corners[7] );
135  std::swap( corners[5], corners[8] );
136 
137  std::swap( corners[12], corners[21] );
138  std::swap( corners[13], corners[22] );
139  std::swap( corners[14], corners[23] );
140 
141  std::swap( corners[15], corners[18] );
142  std::swap( corners[16], corners[19] );
143  std::swap( corners[17], corners[20] );
144  }
std::swap
void swap(edm::DataFrameContainer &lhs, edm::DataFrameContainer &rhs)
Definition: DataFrameContainer.h:223
const TGGC & fireworks::italicGC ( )

Definition at line 53 of file GlobalContexts.cc.

References asciidump::attributes, and fetchall_from_DQM_v2::pool.

Referenced by FWPSetTableManager::cellRenderer().

54 {
55  static TGGC s_italicGC(*gClient->GetResourcePool()->GetFrameGC());
56 
57  TGFontPool *pool = gClient->GetFontPool();
58  TGFont *font = pool->FindFontByHandle(s_italicGC.GetFont());
59  FontAttributes_t attributes = font->GetFontAttributes();
60 
61  attributes.fSlant = 1;
62  TGFont *newFont = pool->GetFont(attributes.fFamily, 9,
63  attributes.fWeight, attributes.fSlant);
64 
65  s_italicGC.SetFont(newFont->GetFontHandle());
66 
67  return s_italicGC;
68 }
fetchall_from_DQM_v2.pool
tuple pool
Definition: fetchall_from_DQM_v2.py:139
asciidump.attributes
list attributes
Definition: asciidump.py:415
void fireworks::localSiStrip ( short  strip,
float *  localTop,
float *  localBottom,
const float *  pars,
unsigned int  id 
)

Definition at line 291 of file TrackUtils.cc.

References delta, fwLog, fwlog::kError, evf::evtn::offset(), funct::tan(), and ecaldqm::topology().

Referenced by addSiStripClusters(), FWSiStripClusterProxyBuilder::build(), and FWSiStripDigiProxyBuilder::build().

292 {
293  Float_t topology = pars[0];
294  Float_t halfStripLength = pars[2] * 0.5;
295 
296  Double_t localCenter[3] = { 0.0, 0.0, 0.0 };
297  localTop[1] = halfStripLength;
298  localBottom[1] = -halfStripLength;
299 
300  if( topology == 1 ) // RadialStripTopology
301  {
302  // stripAngle = phiOfOneEdge + strip * angularWidth
303  // localY = originToIntersection * tan( stripAngle )
304  Float_t stripAngle = tan( pars[5] + strip * pars[6] );
305  Float_t delta = halfStripLength * stripAngle;
306  localCenter[0] = pars[4] * stripAngle;
307  localTop[0] = localCenter[0] + delta;
308  localBottom[0] = localCenter[0] - delta;
309  }
310  else if( topology == 2 ) // RectangularStripTopology
311  {
312  // offset = -numberOfStrips/2. * pitch
313  // localY = strip * pitch + offset
314  Float_t offset = -pars[1] * 0.5 * pars[3];
315  localCenter[0] = strip * pars[3] + offset;
316  localTop[0] = localCenter[0];
317  localBottom[0] = localCenter[0];
318  }
319  else if( topology == 3 ) // TrapezoidalStripTopology
320  {
321  fwLog( fwlog::kError )
322  << "did not expect TrapezoidalStripTopology of "
323  << id << std::endl;
324  }
325  else if( pars[0] == 0 ) // StripTopology
326  {
327  fwLog( fwlog::kError )
328  << "did not find StripTopology of "
329  << id << std::endl;
330  }
331 }
delta
dbl * delta
Definition: mlp_gen.cc:36
ecaldqm::topology
CaloTopology const * topology(0)
funct::tan
Tan< T >::type tan(const T &t)
Definition: Tan.h:22
evf::evtn::offset
unsigned int offset(bool)
Definition: GlobalEventNumber.cc:66
fwLog
#define fwLog(_level_)
Definition: fwLog.h:50
TEveJetCone * fireworks::makeEveJetCone ( const reco::Jet &  iData,
const fireworks::Context &  context 
)

Definition at line 7 of file makeEveJetCone.cc.

References trackerHits::c, fireworks::Context::caloMaxEta(), fireworks::Context::caloR1(), fireworks::Context::caloR2(), fireworks::Context::caloTransAngle(), fireworks::Context::caloZ1(), fireworks::Context::caloZ2(), reco::LeafCandidate::eta(), reco::Jet::etaetaMoment(), reco::Jet::getJetConstituents(), metsig::jet, reco::LeafCandidate::phi(), reco::Jet::phiphiMoment(), mathSSE::sqrt(), funct::tan(), and reco::LeafCandidate::vertex().

Referenced by FWTauProxyBuilderBase::buildBaseTau(), and FWJetProxyBuilder::requestCommon().

8 {
9  TEveJetCone* jet = new TEveJetCone();
10  jet->SetApex(TEveVector(iData.vertex().x(),iData.vertex().y(),iData.vertex().z()));
11 
12  reco::Jet::Constituents c = iData.getJetConstituents();
13  bool haveData = true;
14  for ( reco::Jet::Constituents::const_iterator itr = c.begin(); itr != c.end(); ++itr )
15  {
16  if ( !itr->isAvailable() ) {
17  haveData = false;
18  break;
19  }
20  }
21 
22  double eta_size = 0.2;
23  double phi_size = 0.2;
24  if ( haveData ){
25  eta_size = sqrt(iData.etaetaMoment());
26  phi_size = sqrt(iData.phiphiMoment());
27  }
28 
29  static const float offr = 5;
30  static const float offz = offr/tan(context.caloTransAngle());
31  if (iData.eta() < context.caloMaxEta())
32  jet->SetCylinder(context.caloR1(false) -offr, context.caloZ1(false)-offz);
33  else
34  jet->SetCylinder(context.caloR2(false) -offr, context.caloZ2(false)-offz);
35 
36 
37  jet-> AddEllipticCone(iData.eta(), iData.phi(), eta_size, phi_size);
38  jet->SetPickable(kTRUE);
39  return jet;
40 }
fireworks::Context::caloTransAngle
static float caloTransAngle()
Definition: Context.cc:234
reco::LeafCandidate::vertex
virtual const Point & vertex() const
vertex position (overwritten by PF...)
Definition: LeafCandidate.h:225
reco::LeafCandidate::phi
virtual float phi() const
momentum azimuthal angle
Definition: LeafCandidate.h:180
fireworks::Context::caloZ2
static float caloZ2(bool offset=true)
Definition: Context.cc:224
reco::Jet::Constituents
std::vector< Constituent > Constituents
Definition: Jet.h:23
fireworks::Context::caloR2
static float caloR2(bool offset=true)
Definition: Context.cc:214
reco::LeafCandidate::eta
virtual float eta() const
momentum pseudorapidity
Definition: LeafCandidate.h:184
mathSSE::sqrt
T sqrt(T t)
Definition: SSEVec.h:48
fireworks::Context::caloR1
static float caloR1(bool offset=true)
Definition: Context.cc:209
funct::tan
Tan< T >::type tan(const T &t)
Definition: Tan.h:22
reco::Jet::etaetaMoment
float etaetaMoment() const
eta-eta second moment, ET weighted
Definition: Jet.cc:223
reco::Jet::phiphiMoment
float phiphiMoment() const
phi-phi second moment, ET weighted
Definition: Jet.cc:240
fireworks::Context::caloMaxEta
static double caloMaxEta()
Definition: Context.cc:239
fireworks::Context::caloZ1
static float caloZ1(bool offset=true)
Definition: Context.cc:219
reco::Jet::getJetConstituents
virtual Constituents getJetConstituents() const
list of constituents
Definition: Jet.cc:350
bool fireworks::makeRhoPhiSuperCluster ( FWProxyBuilderBase *  pb,
const reco::SuperClusterRef &  iCluster,
float  iPhi,
TEveElement &  oItemHolder 
)

Definition at line 28 of file makeSuperCluster.cc.

References fireworks::Context::caloR1(), FWDisplayProperties::color(), FWProxyBuilderBase::context(), FWEventItem::defaultDisplayProperties(), end, FWGeometry::getCorners(), fireworks::Context::getGeom(), getPhiRange(), getShape(), i, edm::Ref< C, T, F >::isAvailable(), FWProxyBuilderBase::item(), M_PI, colinearityKinematic::Phi, alignCSCRings::r, and FWProxyBuilderBase::setupAddElement().

Referenced by FWPhotonProxyBuilder::buildViewType(), and FWElectronProxyBuilder::buildViewType().

32 {
33  if( !iCluster.isAvailable()) return false;
34  TEveGeoManagerHolder gmgr( TEveGeoShape::GetGeoMangeur());
35 
36  std::vector< std::pair<DetId, float> > detids = iCluster->hitsAndFractions();
37  std::vector<double> phis;
38  for( std::vector<std::pair<DetId, float> >::const_iterator id = detids.begin(), end = detids.end(); id != end; ++id )
39  {
40  const float* corners = pb->context().getGeom()->getCorners( id->first.rawId());
41  if( corners != 0 )
42  {
43  std::vector<float> centre( 3, 0 );
44 
45  for( unsigned int i = 0; i < 24; i += 3 )
46  {
47  centre[0] += corners[i];
48  centre[1] += corners[i + 1];
49  centre[2] += corners[i + 2];
50  }
51 
52  phis.push_back( TEveVector( centre[0], centre[1], centre[2] ).Phi());
53  }
54  }
55  std::pair<double,double> phiRange = fireworks::getPhiRange( phis, iPhi );
56  const double r = pb->context().caloR1();
57  TGeoBBox *sc_box = new TGeoTubeSeg( r - 2, r , 1,
58  phiRange.first * 180 / M_PI - 0.5,
59  phiRange.second * 180 / M_PI + 0.5 ); // 0.5 is roughly half size of a crystal
60  TEveGeoShape *sc = fireworks::getShape( "supercluster", sc_box, pb->item()->defaultDisplayProperties().color());
61  sc->SetPickable( kTRUE );
62  pb->setupAddElement( sc, &oItemHolder );
63  return true;
64 }
FWProxyBuilderBase::context
const fireworks::Context & context() const
Definition: FWProxyBuilderBase.cc:503
i
int i
Definition: DBlmapReader.cc:9
FWEventItem::defaultDisplayProperties
const FWDisplayProperties & defaultDisplayProperties() const
Definition: FWEventItem.cc:451
fireworks::getPhiRange
std::pair< double, double > getPhiRange(const std::vector< double > &phis, double phi)
Definition: BuilderUtils.cc:21
fireworks::Context::getGeom
const FWGeometry * getGeom() const
Definition: Context.h:83
FWProxyBuilderBase::setupAddElement
void setupAddElement(TEveElement *el, TEveElement *parent, bool set_color=true) const
Definition: FWProxyBuilderBase.cc:429
edm::Ref::isAvailable
bool isAvailable() const
Definition: Ref.h:276
FWProxyBuilderBase::item
const FWEventItem * item() const
Definition: FWProxyBuilderBase.h:67
FWDisplayProperties::color
Color_t color() const
Definition: FWDisplayProperties.h:51
fireworks::Context::caloR1
static float caloR1(bool offset=true)
Definition: Context.cc:209
end
#define end
Definition: vmac.h:37
M_PI
#define M_PI
Definition: readTestVector.cc:25
FWGeometry::getCorners
const float * getCorners(unsigned int id) const
Definition: FWGeometry.cc:277
fireworks::getShape
TEveGeoShape * getShape(const char *name, TGeoBBox *shape, Color_t color)
Definition: BuilderUtils.cc:42
bool fireworks::makeRhoZSuperCluster ( FWProxyBuilderBase *  pb,
const reco::SuperClusterRef &  iCluster,
float  iPhi,
TEveElement &  oItemHolder 
)

Definition at line 66 of file makeSuperCluster.cc.

References addRhoZEnergyProjection(), Reference_intrackfit_cff::barrel, fireworks::Context::caloR1(), fireworks::Context::caloR2(), fireworks::Context::caloTransAngle(), fireworks::Context::caloZ1(), fireworks::Context::caloZ2(), FWProxyBuilderBase::context(), end, FWGeometry::getCorners(), fireworks::Context::getGeom(), i, edm::Ref< C, T, F >::isAvailable(), Pi, and theta().

Referenced by FWPhotonProxyBuilder::buildViewType(), and FWElectronProxyBuilder::buildViewType().

70 {
71  if( !iCluster.isAvailable()) return false;
72  TEveGeoManagerHolder gmgr( TEveGeoShape::GetGeoMangeur());
73  double theta_max = 0;
74  double theta_min = 10;
75  std::vector<std::pair<DetId, float> > detids = iCluster->hitsAndFractions();
76  for( std::vector<std::pair<DetId, float> >::const_iterator id = detids.begin(), end = detids.end(); id != end; ++id )
77  {
78  const float* corners = pb->context().getGeom()->getCorners( id->first.rawId());
79  if( corners != 0 )
80  {
81  std::vector<float> centre( 3, 0 );
82 
83  for( unsigned int i = 0; i < 24; i += 3 )
84  {
85  centre[0] += corners[i];
86  centre[1] += corners[i + 1];
87  centre[2] += corners[i + 2];
88  }
89 
90  double theta = TEveVector( centre[0], centre[1], centre[2] ).Theta();
91  if( theta > theta_max ) theta_max = theta;
92  if( theta < theta_min ) theta_min = theta;
93  }
94  }
95  // expand theta range by the size of a crystal to avoid segments of zero length
96  bool barrel = true;
97  if ((theta_max > 0 && theta_max < pb->context().caloTransAngle()) ||
98  ( theta_min > (TMath::Pi() -pb->context().caloTransAngle())) )
99  {
100  barrel = false;
101  }
102 
103  double z_ecal = barrel ? pb->context().caloZ1() : pb->context().caloZ2();
104  double r_ecal = barrel ? pb->context().caloR1() : pb->context().caloR2();
105 
106  fireworks::addRhoZEnergyProjection( pb, &oItemHolder, r_ecal-1, z_ecal-1,
107  theta_min - 0.003, theta_max + 0.003,
108  iPhi );
109 
110  return true;
111 }
Pi
const double Pi
Definition: CosmicMuonParameters.h:18
FWProxyBuilderBase::context
const fireworks::Context & context() const
Definition: FWProxyBuilderBase.cc:503
i
int i
Definition: DBlmapReader.cc:9
Reference_intrackfit_cff.barrel
barrel
Definition: Reference_intrackfit_cff.py:37
fireworks::Context::caloTransAngle
static float caloTransAngle()
Definition: Context.cc:234
fireworks::Context::getGeom
const FWGeometry * getGeom() const
Definition: Context.h:83
theta
Geom::Theta< T > theta() const
Definition: Basic3DVectorLD.h:170
fireworks::Context::caloZ2
static float caloZ2(bool offset=true)
Definition: Context.cc:224
edm::Ref::isAvailable
bool isAvailable() const
Definition: Ref.h:276
fireworks::Context::caloR2
static float caloR2(bool offset=true)
Definition: Context.cc:214
fireworks::Context::caloR1
static float caloR1(bool offset=true)
Definition: Context.cc:209
end
#define end
Definition: vmac.h:37
fireworks::addRhoZEnergyProjection
void addRhoZEnergyProjection(FWProxyBuilderBase *, TEveElement *, double r_ecal, double z_ecal, double theta_min, double theta_max, double phi)
Definition: BuilderUtils.cc:60
FWGeometry::getCorners
const float * getCorners(unsigned int id) const
Definition: FWGeometry.cc:277
fireworks::Context::caloZ1
static float caloZ1(bool offset=true)
Definition: Context.cc:219
float fireworks::pixelLocalX ( const double  mpx,
const int  m_nrows 
)

Definition at line 164 of file TrackUtils.cc.

References BIG_PIX_PER_ROC_X, PITCHX, and ROWS_PER_ROC.

Referenced by FWSiPixelClusterProxyBuilder::build(), FWSiPixelDigiProxyBuilder::build(), pushPixelCluster(), and ValidateGeometry::validatePixelTopology().

165 {
166  // Calculate the edge of the active sensor with respect to the center,
167  // that is simply the half-size.
168  // Take into account large pixels
169  const double xoffset = -( nrows + BIG_PIX_PER_ROC_X * nrows / ROWS_PER_ROC ) / 2. * PITCHX;
170 
171  // Measurement to local transformation for X coordinate
172  // X coordinate is in the ROC row number direction
173  // Copy from RectangularPixelTopology::localX implementation
174  int binoffx = int(mpx); // truncate to int
175  double fractionX = mpx - binoffx; // find the fraction
176  double local_PITCHX = PITCHX; // defaultpitch
177  if( binoffx > 80 ) { // ROC 1 - handles x on edge cluster
178  binoffx = binoffx + 2;
179  } else if( binoffx == 80 ) { // ROC 1
180  binoffx = binoffx+1;
181  local_PITCHX = 2 * PITCHX;
182  } else if( binoffx == 79 ) { // ROC 0
183  binoffx = binoffx + 0;
184  local_PITCHX = 2 * PITCHX;
185  } else if( binoffx >= 0 ) { // ROC 0
186  binoffx = binoffx + 0;
187  }
188 
189  // The final position in local coordinates
190  double lpX = double( binoffx * PITCHX ) + fractionX * local_PITCHX + xoffset;
191 
192  return lpX;
193 }
fireworks::ROWS_PER_ROC
static const int ROWS_PER_ROC
Definition: TrackUtils.cc:62
fireworks::PITCHX
static const double PITCHX
Definition: TrackUtils.cc:58
fireworks::BIG_PIX_PER_ROC_X
static const int BIG_PIX_PER_ROC_X
Definition: TrackUtils.cc:63
float fireworks::pixelLocalY ( const double  mpy,
const int  m_ncols 
)

Definition at line 197 of file TrackUtils.cc.

References BIG_PIX_PER_ROC_Y, COLS_PER_ROC, and PITCHY.

Referenced by FWSiPixelClusterProxyBuilder::build(), FWSiPixelDigiProxyBuilder::build(), pushPixelCluster(), and ValidateGeometry::validatePixelTopology().

198 {
199  // Calculate the edge of the active sensor with respect to the center,
200  // that is simply the half-size.
201  // Take into account large pixels
202  double yoffset = -( ncols + BIG_PIX_PER_ROC_Y * ncols / COLS_PER_ROC ) / 2. * PITCHY;
203 
204  // Measurement to local transformation for Y coordinate
205  // Y is in the ROC column number direction
206  // Copy from RectangularPixelTopology::localY implementation
207  int binoffy = int( mpy ); // truncate to int
208  double fractionY = mpy - binoffy; // find the fraction
209  double local_PITCHY = PITCHY; // defaultpitch
210 
211  if( binoffy>416 ) { // ROC 8, not real ROC
212  binoffy = binoffy+17;
213  } else if( binoffy == 416 ) { // ROC 8
214  binoffy = binoffy + 16;
215  local_PITCHY = 2 * PITCHY;
216  } else if( binoffy == 415 ) { // ROC 7, last big pixel
217  binoffy = binoffy + 15;
218  local_PITCHY = 2 * PITCHY;
219  } else if( binoffy > 364 ) { // ROC 7
220  binoffy = binoffy + 15;
221  } else if( binoffy == 364 ) { // ROC 7
222  binoffy = binoffy + 14;
223  local_PITCHY = 2 * PITCHY;
224  } else if( binoffy == 363 ) { // ROC 6
225  binoffy = binoffy + 13;
226  local_PITCHY = 2 * PITCHY;
227  } else if( binoffy > 312 ) { // ROC 6
228  binoffy = binoffy + 13;
229  } else if( binoffy == 312 ) { // ROC 6
230  binoffy = binoffy + 12;
231  local_PITCHY = 2 * PITCHY;
232  } else if( binoffy == 311 ) { // ROC 5
233  binoffy = binoffy + 11;
234  local_PITCHY = 2 * PITCHY;
235  } else if( binoffy > 260 ) { // ROC 5
236  binoffy = binoffy + 11;
237  } else if( binoffy == 260 ) { // ROC 5
238  binoffy = binoffy + 10;
239  local_PITCHY = 2 * PITCHY;
240  } else if( binoffy == 259 ) { // ROC 4
241  binoffy = binoffy + 9;
242  local_PITCHY = 2 * PITCHY;
243  } else if( binoffy > 208 ) { // ROC 4
244  binoffy = binoffy + 9;
245  } else if(binoffy == 208 ) { // ROC 4
246  binoffy = binoffy + 8;
247  local_PITCHY = 2 * PITCHY;
248  } else if( binoffy == 207 ) { // ROC 3
249  binoffy = binoffy + 7;
250  local_PITCHY = 2 * PITCHY;
251  } else if( binoffy > 156 ) { // ROC 3
252  binoffy = binoffy + 7;
253  } else if( binoffy == 156 ) { // ROC 3
254  binoffy = binoffy + 6;
255  local_PITCHY = 2 * PITCHY;
256  } else if( binoffy == 155 ) { // ROC 2
257  binoffy = binoffy + 5;
258  local_PITCHY = 2 * PITCHY;
259  } else if( binoffy > 104 ) { // ROC 2
260  binoffy = binoffy + 5;
261  } else if( binoffy == 104 ) { // ROC 2
262  binoffy = binoffy + 4;
263  local_PITCHY = 2 * PITCHY;
264  } else if( binoffy == 103 ) { // ROC 1
265  binoffy = binoffy + 3;
266  local_PITCHY = 2 * PITCHY;
267  } else if( binoffy > 52 ) { // ROC 1
268  binoffy = binoffy + 3;
269  } else if( binoffy == 52 ) { // ROC 1
270  binoffy = binoffy + 2;
271  local_PITCHY = 2 * PITCHY;
272  } else if( binoffy == 51 ) { // ROC 0
273  binoffy = binoffy + 1;
274  local_PITCHY = 2 * PITCHY;
275  } else if( binoffy > 0 ) { // ROC 0
276  binoffy=binoffy + 1;
277  } else if( binoffy == 0 ) { // ROC 0
278  binoffy = binoffy + 0;
279  local_PITCHY = 2 * PITCHY;
280  }
281 
282  // The final position in local coordinates
283  double lpY = double( binoffy * PITCHY ) + fractionY * local_PITCHY + yoffset;
284 
285  return lpY;
286 }
fireworks::COLS_PER_ROC
static const int COLS_PER_ROC
Definition: TrackUtils.cc:61
fireworks::BIG_PIX_PER_ROC_Y
static const int BIG_PIX_PER_ROC_Y
Definition: TrackUtils.cc:60
fireworks::PITCHY
static const double PITCHY
Definition: TrackUtils.cc:59
TEveTrack * fireworks::prepareCandidate ( const reco::Candidate &  track,
TEveTrackPropagator *  propagator 
)

Definition at line 10 of file CandidateUtils.cc.

References reco::Candidate::charge(), reco::Candidate::px(), reco::Candidate::py(), reco::Candidate::pz(), edmStreamStallGrapher::t, and reco::Candidate::vertex().

Referenced by FWElectronLegoProxyBuilder::build(), FWCandidateProxyBuilder::build(), FWGenParticleProxyBuilder::build(), FWMuonBuilder::buildMuon(), and FWElectronProxyBuilder::requestCommon().

12 {
13  TEveRecTrack t;
14  t.fBeta = 1.;
15  t.fP = TEveVector( track.px(), track.py(), track.pz() );
16  t.fV = TEveVector( track.vertex().x(), track.vertex().y(), track.vertex().z() );
17  t.fSign = track.charge();
18  TEveTrack* trk = new TEveTrack(&t, propagator);
19  return trk;
20 }
reco::Candidate::pz
virtual double pz() const =0
z coordinate of momentum vector
reco::Candidate::vertex
virtual const Point & vertex() const =0
vertex position
reco::Candidate::charge
virtual int charge() const =0
electric charge
reco::Candidate::py
virtual double py() const =0
y coordinate of momentum vector
reco::Candidate::px
virtual double px() const =0
x coordinate of momentum vector
TEveTrack * fireworks::prepareTrack ( const reco::Track &  track,
TEveTrackPropagator *  propagator,
const std::vector< TEveVector > &  extraRefPoints = std::vector<TEveVector>() 
)

Definition at line 69 of file TrackUtils.cc.

References reco::TrackBase::charge(), reco::Track::extra(), f, i, reco::Track::innerMomentum(), reco::Track::innerOk(), reco::Track::innerPosition(), edm::Ref< C, T, F >::isAvailable(), j, reco::Track::outerMomentum(), reco::Track::outerOk(), reco::Track::outerPosition(), AlCaHLTBitMon_ParallelJobs::p, point, position, reco::TrackBase::pt(), reco::TrackBase::px(), reco::TrackBase::py(), reco::TrackBase::pz(), python.multivaluedict::sort(), edmStreamStallGrapher::t, reco::btau::trackMomentum, findQualityFiles::v, TrackValidation_HighPurity_cff::valid, reco::TrackBase::vx(), reco::TrackBase::vy(), and reco::TrackBase::vz().

Referenced by FWTauProxyBuilderBase::addConstituentTracks(), FWTauProxyBuilderBase::addLeadTrack(), DummyEvelyser::analyze(), FWElectronLegoProxyBuilder::build(), FWTrackProxyBuilder::build(), FWTrackHitsDetailView::build(), FWConvTrackHitsDetailView::build(), FWMuonBuilder::buildMuon(), and FWElectronProxyBuilder::requestCommon().

72 {
73  // To make use of all available information, we have to order states
74  // properly first. Propagator should take care of y=0 transition.
75 
76  std::vector<State> refStates;
77  TEveVector trackMomentum(track.px(), track.py(), track.pz());
78  refStates.push_back(State(TEveVector(track.vx(), track.vy(), track.vz()),
79  trackMomentum));
80  if( track.extra().isAvailable() ) {
81  if (track.innerOk()) {
82  const reco::TrackBase::Point &v = track.innerPosition();
83  const reco::TrackBase::Vector &p = track.innerMomentum();
84  refStates.push_back(State(TEveVector(v.x(), v.y(), v.z()), TEveVector(p.x(), p.y(), p.z())));
85  }
86  if (track.outerOk()) {
87  const reco::TrackBase::Point &v = track.outerPosition();
88  const reco::TrackBase::Vector &p = track.outerMomentum();
89  refStates.push_back(State(TEveVector(v.x(), v.y(), v.z()), TEveVector(p.x(), p.y(), p.z())));
90  }
91  }
92  for( std::vector<TEveVector>::const_iterator point = extraRefPoints.begin(), pointEnd = extraRefPoints.end();
93  point != pointEnd; ++point )
94  refStates.push_back(State(*point));
95  if( track.pt()>1 )
96  std::sort( refStates.begin(), refStates.end(), StateOrdering(trackMomentum) );
97 
98  // * if the first state has non-zero momentum use it as a starting point
99  // and all other points as PathMarks to follow
100  // * if the first state has only position, try the last state. If it has
101  // momentum we propagate backword, if not, we look for the first one
102  // on left that has momentum and ignore all earlier.
103  //
104 
105  TEveRecTrack t;
106  t.fBeta = 1.;
107  t.fSign = track.charge();
108 
109  if( refStates.front().valid ) {
110  t.fV = refStates.front().position;
111  t.fP = refStates.front().momentum;
112  TEveTrack* trk = new TEveTrack( &t, propagator );
113  for( unsigned int i(1); i<refStates.size()-1; ++i) {
114  if( refStates[i].valid )
115  trk->AddPathMark( TEvePathMark( TEvePathMark::kReference, refStates[i].position, refStates[i].momentum ) );
116  else
117  trk->AddPathMark( TEvePathMark( TEvePathMark::kDaughter, refStates[i].position ) );
118  }
119  if( refStates.size()>1 ) {
120  trk->AddPathMark( TEvePathMark( TEvePathMark::kDecay, refStates.back().position ) );
121  }
122  return trk;
123  }
124 
125  if( refStates.back().valid ) {
126  t.fSign = (-1)*track.charge();
127  t.fV = refStates.back().position;
128  t.fP = refStates.back().momentum * (-1.0f);
129  TEveTrack* trk = new TEveTrack( &t, propagator );
130  unsigned int i( refStates.size()-1 );
131  for(; i>0; --i) {
132  if ( refStates[i].valid )
133  trk->AddPathMark( TEvePathMark( TEvePathMark::kReference, refStates[i].position, refStates[i].momentum*(-1.0f) ) );
134  else
135  trk->AddPathMark( TEvePathMark( TEvePathMark::kDaughter, refStates[i].position ) );
136  }
137  if ( refStates.size()>1 ) {
138  trk->AddPathMark( TEvePathMark( TEvePathMark::kDecay, refStates.front().position ) );
139  }
140  return trk;
141  }
142 
143  unsigned int i(0);
144  while( i<refStates.size() && !refStates[i].valid ) ++i;
145  assert( i < refStates.size() );
146 
147  t.fV = refStates[i].position;
148  t.fP = refStates[i].momentum;
149  TEveTrack* trk = new TEveTrack( &t, propagator );
150  for( unsigned int j(i+1); j<refStates.size()-1; ++j ) {
151  if( refStates[i].valid )
152  trk->AddPathMark( TEvePathMark( TEvePathMark::kReference, refStates[i].position, refStates[i].momentum ) );
153  else
154  trk->AddPathMark( TEvePathMark( TEvePathMark::kDaughter, refStates[i].position ) );
155  }
156  if ( i < refStates.size() ) {
157  trk->AddPathMark( TEvePathMark( TEvePathMark::kDecay, refStates.back().position ) );
158  }
159  return trk;
160 }
TrackValidation_HighPurity_cff.valid
tuple valid
Definition: TrackValidation_HighPurity_cff.py:20
i
int i
Definition: DBlmapReader.cc:9
reco::Track::extra
const TrackExtraRef & extra() const
reference to &quot;extra&quot; object
Definition: Track.h:97
reco::Track::innerOk
bool innerOk() const
return true if the innermost hit is valid
Definition: Track.h:40
reco::TrackBase::px
double px() const
x coordinate of momentum vector
Definition: TrackBase.h:131
edm::Ref::isAvailable
bool isAvailable() const
Definition: Ref.h:276
reco::Track::outerPosition
const math::XYZPoint & outerPosition() const
position of the outermost hit
Definition: Track.h:47
reco::Track::innerPosition
const math::XYZPoint & innerPosition() const
position of the innermost hit
Definition: Track.h:42
reco::TrackBase::pt
double pt() const
track transverse momentum
Definition: TrackBase.h:129
j
int j
Definition: DBlmapReader.cc:9
f
double f[11][100]
Definition: MuScleFitUtils.cc:78
reco::TrackBase::Point
math::XYZPoint Point
point in the space
Definition: TrackBase.h:74
reco::TrackBase::pz
double pz() const
z coordinate of momentum vector
Definition: TrackBase.h:135
reco::TrackBase::vz
double vz() const
z coordinate of the reference point on track
Definition: TrackBase.h:145
reco::Track::outerMomentum
const math::XYZVector & outerMomentum() const
momentum vector at the outermost hit position
Definition: Track.h:49
reco::Track::outerOk
bool outerOk() const
return true if the outermost hit is valid
Definition: Track.h:38
AlCaHLTBitMon_ParallelJobs.p
tuple p
Definition: AlCaHLTBitMon_ParallelJobs.py:152
reco::TrackBase::vy
double vy() const
y coordinate of the reference point on track
Definition: TrackBase.h:143
reco::Track::innerMomentum
const math::XYZVector & innerMomentum() const
momentum vector at the innermost hit position
Definition: Track.h:45
position
static int position[264][3]
Definition: ReadPGInfo.cc:509
reco::TrackBase::charge
int charge() const
track electric charge
Definition: TrackBase.h:111
reco::TrackBase::Vector
math::XYZVector Vector
spatial vector
Definition: TrackBase.h:72
point
*vegas h *****************************************************used in the default bin number in original ***version of VEGAS is ***a higher bin number might help to derive a more precise ***grade subtle point
Definition: invegas.h:5
reco::TrackBase::py
double py() const
y coordinate of momentum vector
Definition: TrackBase.h:133
reco::TrackBase::vx
double vx() const
x coordinate of the reference point on track
Definition: TrackBase.h:141
void fireworks::pushNearbyPixelHits ( std::vector< TVector3 > &  pixelPoints,
const FWEventItem &  iItem,
const reco::Track &  t 
)

Definition at line 543 of file TrackUtils.cc.

References edmNew::DetSet< T >::begin(), SiPixelRecHit::cluster(), FWGeometry::contains(), edmNew::DetSet< T >::end(), edmNew::DetSetVector< T >::end(), edmNew::DetSetVector< T >::find(), fwLog, relativeConstraints::geom, FWEventItem::getGeom(), FWGeometry::getParameters(), edm::Ref< C, T, F >::isAvailable(), edm::Ref< C, T, F >::isNonnull(), fwlog::kError, PixelSubdetector::PixelBarrel, PixelSubdetector::PixelEndcap, edm::Ref< C, T, F >::product(), pushPixelCluster(), reco::Track::recHitsBegin(), and reco::Track::recHitsEnd().

Referenced by FWTrackHitsDetailView::addHits(), and FWConvTrackHitsDetailView::addHits().

544 {
545  const edmNew::DetSetVector<SiPixelCluster> * allClusters = 0;
546  for( trackingRecHit_iterator it = t.recHitsBegin(), itEnd = t.recHitsEnd(); it != itEnd; ++it)
547  {
548  if( typeid(**it) == typeid( SiPixelRecHit ))
549  {
550  const SiPixelRecHit &hit = static_cast<const SiPixelRecHit &>(**it);
551  if( hit.cluster().isNonnull() && hit.cluster().isAvailable())
552  {
553  allClusters = hit.cluster().product();
554  break;
555  }
556  }
557  }
558  if( allClusters == 0 ) return;
559 
560  const FWGeometry *geom = iItem.getGeom();
561 
562  for( trackingRecHit_iterator it = t.recHitsBegin(), itEnd = t.recHitsEnd(); it != itEnd; ++it )
563  {
564  const TrackingRecHit* rh = &(**it);
565 
566  DetId id = (*it)->geographicalId();
567  if( ! geom->contains( id ))
568  {
569  fwLog( fwlog::kError )
570  << "failed to get geometry of Tracker Det with raw id: "
571  << id.rawId() << std::endl;
572 
573  continue;
574  }
575 
576  // -- in which detector are we?
577  unsigned int subdet = (unsigned int)id.subdetId();
578  if(( subdet != PixelSubdetector::PixelBarrel ) && ( subdet != PixelSubdetector::PixelEndcap )) continue;
579 
580  const SiPixelCluster* hitCluster = 0;
581  if( const SiPixelRecHit* pixel = dynamic_cast<const SiPixelRecHit*>( rh ))
582  hitCluster = pixel->cluster().get();
583  edmNew::DetSetVector<SiPixelCluster>::const_iterator itds = allClusters->find(id.rawId());
584  if( itds != allClusters->end())
585  {
586  const edmNew::DetSet<SiPixelCluster> & clustersOnThisDet = *itds;
587  for( edmNew::DetSet<SiPixelCluster>::const_iterator itc = clustersOnThisDet.begin(), edc = clustersOnThisDet.end(); itc != edc; ++itc )
588  {
589  if( &*itc != hitCluster )
590  pushPixelCluster( pixelPoints, *geom, id, *itc, geom->getParameters( id ));
591  }
592  }
593  }
594 }
edmNew::DetSetVector::const_iterator
boost::transform_iterator< IterHelp, const_IdIter > const_iterator
Definition: DetSetVectorNew.h:118
edmNew::DetSetVector::end
const_iterator end(bool update=false) const
Definition: DetSetVectorNew.h:331
FWGeometry::getParameters
const float * getParameters(unsigned int id) const
Definition: FWGeometry.cc:293
edmNew::DetSet::const_iterator
data_type const * const_iterator
Definition: DetSetNew.h:30
edm::Ref::isAvailable
bool isAvailable() const
Definition: Ref.h:276
edm::Ref::isNonnull
bool isNonnull() const
Checks for non-null.
Definition: Ref.h:250
FWGeometry::contains
bool contains(unsigned int id) const
Definition: FWGeometry.h:109
reco::Track::recHitsBegin
trackingRecHit_iterator recHitsBegin() const
Iterator to first hit on the track.
Definition: Track.h:63
DetId
Definition: DetId.h:18
fwLog
#define fwLog(_level_)
Definition: fwLog.h:50
SiPixelRecHit::cluster
ClusterRef cluster() const
Definition: SiPixelRecHit.h:49
edmNew::DetSetVector::find
const_iterator find(id_type i, bool update=false) const
Definition: DetSetVectorNew.h:310
SiPixelCluster
Pixel cluster – collection of neighboring pixels above threshold.
Definition: SiPixelCluster.h:33
edmNew::DetSet::end
iterator end()
Definition: DetSetNew.h:70
edm::Ref::product
C const * product() const
Accessor for product collection.
Definition: Ref.h:601
fireworks::pushPixelCluster
void pushPixelCluster(std::vector< TVector3 > &pixelPoints, const FWGeometry &geom, DetId id, const SiPixelCluster &c, const float *pars)
Definition: TrackUtils.cc:646
FWEventItem::getGeom
const FWGeometry * getGeom() const
Definition: FWEventItem.cc:681
SiPixelRecHit
Our base class.
Definition: SiPixelRecHit.h:23
reco::Track::recHitsEnd
trackingRecHit_iterator recHitsEnd() const
Iterator to last hit on the track.
Definition: Track.h:65
edmNew::DetSet::begin
iterator begin()
Definition: DetSetNew.h:67
void fireworks::pushPixelCluster ( std::vector< TVector3 > &  pixelPoints,
const FWGeometry &  geom,
DetId  id,
const SiPixelCluster &  c,
const float *  pars 
)

Definition at line 646 of file TrackUtils.cc.

References cuy::col, fwLog, fwlog::kDebug, FWGeometry::localToGlobal(), SiPixelCluster::minPixelCol(), SiPixelCluster::minPixelRow(), pixelLocalX(), and pixelLocalY().

Referenced by pushNearbyPixelHits(), and pushPixelHits().

647 {
648  double row = c.minPixelRow();
649  double col = c.minPixelCol();
650  float lx = 0.;
651  float ly = 0.;
652 
653  int nrows = (int)pars[0];
654  int ncols = (int)pars[1];
655  lx = pixelLocalX( row, nrows );
656  ly = pixelLocalY( col, ncols );
657 
658  fwLog( fwlog::kDebug )
659  << ", row: " << row << ", col: " << col
660  << ", lx: " << lx << ", ly: " << ly ;
661 
662  float local[3] = { lx, ly, 0. };
663  float global[3];
664  geom.localToGlobal( id, local, global );
665  TVector3 pb( global[0], global[1], global[2] );
666  pixelPoints.push_back( pb );
667 
668  fwLog( fwlog::kDebug )
669  << " x: " << pb.X()
670  << ", y: " << pb.Y()
671  << " z: " << pb.Z()
672  << " eta: " << pb.Eta()
673  << ", phi: " << pb.Phi()
674  << " rho: " << pb.Pt() << std::endl;
675 }
SiPixelCluster::minPixelCol
int minPixelCol() const
Definition: SiPixelCluster.h:155
fireworks::pixelLocalY
float pixelLocalY(const double mpy, const int m_ncols)
Definition: TrackUtils.cc:197
SiPixelCluster::minPixelRow
int minPixelRow() const
Definition: SiPixelCluster.h:153
FWGeometry::localToGlobal
void localToGlobal(unsigned int id, const float *local, float *global, bool translatep=true) const
Definition: FWGeometry.cc:325
fwLog
#define fwLog(_level_)
Definition: fwLog.h:50
cuy.col
int col
Definition: cuy.py:1008
fireworks::pixelLocalX
float pixelLocalX(const double mpx, const int m_nrows)
Definition: TrackUtils.cc:164
void fireworks::pushPixelHits ( std::vector< TVector3 > &  pixelPoints,
const FWEventItem &  iItem,
const reco::Track &  t 
)

Definition at line 599 of file TrackUtils.cc.

References trackerHits::c, FWGeometry::contains(), reco::TrackBase::dz(), reco::TrackBase::eta(), fwLog, relativeConstraints::geom, FWEventItem::getGeom(), FWGeometry::getParameters(), fwlog::kDebug, fwlog::kError, PixelSubdetector::PixelBarrel, PixelSubdetector::PixelEndcap, pushPixelCluster(), reco::Track::recHitsBegin(), reco::Track::recHitsEnd(), subdets, and reco::TrackBase::vz().

Referenced by FWTrackHitsDetailView::addHits(), FWConvTrackHitsDetailView::addHits(), and FWTracksRecHitsProxyBuilder::build().

600 {
601  /*
602  * -- return for each Pixel Hit a 3D point
603  */
604  const FWGeometry *geom = iItem.getGeom();
605 
606  double dz = t.dz();
607  double vz = t.vz();
608  double etaT = t.eta();
609 
610  fwLog( fwlog::kDebug ) << "Track eta: " << etaT << ", vz: " << vz << ", dz: " << dz
611  << std::endl;
612 
613  int cnt = 0;
614  for( trackingRecHit_iterator it = t.recHitsBegin(), itEnd = t.recHitsEnd(); it != itEnd; ++it )
615  {
616  const TrackingRecHit* rh = &(**it);
617  // -- get position of center of wafer, assuming (0,0,0) is the center
618  DetId id = (*it)->geographicalId();
619  if( ! geom->contains( id ))
620  {
621  fwLog( fwlog::kError )
622  << "failed to get geometry of Tracker Det with raw id: "
623  << id.rawId() << std::endl;
624 
625  continue;
626  }
627 
628  // -- in which detector are we?
629  unsigned int subdet = (unsigned int)id.subdetId();
630 
631  if(( subdet == PixelSubdetector::PixelBarrel ) || ( subdet == PixelSubdetector::PixelEndcap ))
632  {
633  fwLog( fwlog::kDebug ) << cnt++ << " -- "
634  << subdets[subdet];
635 
636  if( const SiPixelRecHit* pixel = dynamic_cast<const SiPixelRecHit*>( rh ))
637  {
638  const SiPixelCluster& c = *( pixel->cluster());
639  pushPixelCluster( pixelPoints, *geom, id, c, geom->getParameters( id ));
640  }
641  }
642  }
643 }
FWGeometry::getParameters
const float * getParameters(unsigned int id) const
Definition: FWGeometry.cc:293
reco::TrackBase::eta
double eta() const
pseudorapidity of momentum vector
Definition: TrackBase.h:139
FWGeometry::contains
bool contains(unsigned int id) const
Definition: FWGeometry.h:109
reco::Track::recHitsBegin
trackingRecHit_iterator recHitsBegin() const
Iterator to first hit on the track.
Definition: Track.h:63
subdets
char const * subdets[11]
Definition: HcalDetDiagNoiseMonitor.cc:76
reco::TrackBase::dz
double dz() const
dz parameter (= dsz/cos(lambda)). This is the track z0 w.r.t (0,0,0) only if the refPoint is close to...
Definition: TrackBase.h:125
reco::TrackBase::vz
double vz() const
z coordinate of the reference point on track
Definition: TrackBase.h:145
DetId
Definition: DetId.h:18
fwLog
#define fwLog(_level_)
Definition: fwLog.h:50
SiPixelCluster
Pixel cluster – collection of neighboring pixels above threshold.
Definition: SiPixelCluster.h:33
fireworks::pushPixelCluster
void pushPixelCluster(std::vector< TVector3 > &pixelPoints, const FWGeometry &geom, DetId id, const SiPixelCluster &c, const float *pars)
Definition: TrackUtils.cc:646
FWEventItem::getGeom
const FWGeometry * getGeom() const
Definition: FWEventItem.cc:681
SiPixelRecHit
Our base class.
Definition: SiPixelRecHit.h:23
reco::Track::recHitsEnd
trackingRecHit_iterator recHitsEnd() const
Iterator to last hit on the track.
Definition: Track.h:65
const TGGC & fireworks::redGC ( )

Definition at line 46 of file GlobalContexts.cc.

47 {
48  static TGGC s_redGC(*gClient->GetResourcePool()->GetFrameGC());
49  s_redGC.SetForeground(gVirtualX->GetPixel(kRed-5));
50  return s_redGC;
51 }
void fireworks::setPath ( TString &  v)

Definition at line 15 of file fwPaths.cc.

References fwLog, fwlog::kError, p1, and p2.

Referenced by CmsShowHelpPopup::CmsShowHelpPopup(), main(), objMon< T >::monitorDenominator(), CmsShowMainBase::sendVersionInfo(), and CmsShowMainFrame::showFWorksInfo().

16 {
17  if (gSystem->AccessPathName(p1 + v) == kFALSE)
18  {
19  v.Prepend(p1);
20  return;
21  }
22 
23  v.Prepend(p2);
24  if (gSystem->AccessPathName(v))
25  fwLog(fwlog::kError) << "Can't access path " << v << std::endl;
26 }
p2
double p2[4]
Definition: TauolaWrapper.h:90
fwLog
#define fwLog(_level_)
Definition: fwLog.h:50
p1
double p1[4]
Definition: TauolaWrapper.h:89
void fireworks::setTrackTypePF ( const reco::PFCandidate &  pfCand,
TAttLine *  track 
)

Definition at line 17 of file setTrackTypePF.cc.

References dt_dqm_sourceclient_common_cff::reco.

Referenced by FWPFCandidateWithHitsProxyBuilder::build(), FWPFPatJetLegoProxyBuilder< T >::build(), FWPFPatJet3DProxyBuilder< T >::build(), FWPFCandidate3DProxyBuilder::build(), and FWPFCandidatesLegoProxyBuilder::build().

18  {
19 
20  using namespace reco;
21 
22  // switch (pfCand.particleId() ) {
23  // case PFCandidate::e: track->SetLineStyle(5); break;
24  // case PFCandidate::mu: track->SetLineStyle(6); break;
25  // case PFCandidate::h0: track->SetLineStyle(3); break;
26  // case PFCandidate::gamma: track->SetLineStyle(7); break;
27  // default: break;
28  // }
29  }
dt_dqm_sourceclient_common_cff.reco
tuple reco
Definition: dt_dqm_sourceclient_common_cff.py:107
void fireworks::setupAddElement ( TEveElement *  el,
TEveElement *  parent,
const FWEventItem *  item,
bool  master,
bool  color 
)

Definition at line 336 of file TrackUtils.cc.

References FWDisplayProperties::color(), FWEventItem::defaultDisplayProperties(), and FWDisplayProperties::transparency().

Referenced by addSiStripClusters(), and FWDTRecHitProxyBuilder::buildViewType().

337 {
338  if (master)
339  {
340  el->CSCTakeAnyParentAsMaster();
341  el->SetPickable(true);
342  }
343 
344  if (color)
345  {
346  el->CSCApplyMainColorToMatchingChildren();
347  el->CSCApplyMainTransparencyToMatchingChildren();
348  el->SetMainColor(item->defaultDisplayProperties().color());
349  assert((item->defaultDisplayProperties().transparency() >= 0)
350  && (item->defaultDisplayProperties().transparency() <= 100));
351  el->SetMainTransparency(item->defaultDisplayProperties().transparency());
352  }
353  parent->AddElement(el);
354 }
FWEventItem::defaultDisplayProperties
const FWDisplayProperties & defaultDisplayProperties() const
Definition: FWEventItem.cc:451
dbtoconf.parent
list parent
Definition: dbtoconf.py:74
FWDisplayProperties::color
Color_t color() const
Definition: FWDisplayProperties.h:51
FWDisplayProperties::transparency
Char_t transparency() const
Definition: FWDisplayProperties.h:55
autoMagneticFieldProducer_cfi.master
tuple master
Definition: autoMagneticFieldProducer_cfi.py:99
create_public_lumi_plots.color
list color
Definition: create_public_lumi_plots.py:1800
int * fireworks::supportedDataFormatsVersion ( )

Definition at line 43 of file fwPaths.cc.

References getDecomposedVersion(), and findQualityFiles::v.

Referenced by acceptDataFormatsVersion(), CmsShowSearchFiles::CmsShowSearchFiles(), FWFileEntry::openFile(), and CmsShowSearchFiles::readInfo().

44 {
45  static int mm[] = {0, 0, 0};
46 
47  if (!mm[0]) {
48  TString v;
49  if (gSystem->Getenv("CMSSW_VERSION"))
50  {
51  v = gSystem->Getenv("CMSSW_VERSION");
52 
53  }
54  else
55  {
56  TString versionFileName = gSystem->Getenv("CMSSW_BASE");
57  versionFileName += "/src/Fireworks/Core/data/version.txt";
58  std::ifstream fs(versionFileName);
59  TString infoText;
60  infoText.ReadLine(fs); infoText.ReadLine(fs);
61  fs.close();
62  v = &infoText[13];
63  }
64 
65  getDecomposedVersion(v, &mm[0]);
66  }
67 
68  return &mm[0];
69 }
fireworks::getDecomposedVersion
void getDecomposedVersion(const TString &s, int *out)
Definition: fwPaths.cc:30
std::vector< std::pair< double, double > > fireworks::thetaBins ( )

Definition at line 48 of file thetaBins.cc.

References create_public_lumi_plots::exp, i, fw3dlego::xbins, and fw3dlego::xbins_n.

Referenced by FWCaloTauProxyBuilder::buildViewType().

49  {
50  const int n_bins = fw3dlego::xbins_n - 1;
51  std::vector<std::pair<double,double> > thetaBins(n_bins);
52  for (int i = 0; i < n_bins; ++i )
53  {
54  thetaBins[i].first = 2*atan( exp(-fw3dlego::xbins[i]) );
55  thetaBins[i].second = 2*atan( exp(-fw3dlego::xbins[i+1]) );
56  }
57  return thetaBins;
58  }
i
int i
Definition: DBlmapReader.cc:9
create_public_lumi_plots.exp
tuple exp
Definition: create_public_lumi_plots.py:1087
fw3dlego::xbins
const double xbins[]
Definition: fw3dlego_xbins.cc:17
fireworks::thetaBins
std::vector< std::pair< double, double > > thetaBins()
Definition: thetaBins.cc:48
fw3dlego::xbins_n
const int xbins_n
Definition: fw3dlego_xbins.cc:16

Variable Documentation

const int fireworks::BIG_PIX_PER_ROC_X = 1
static

Definition at line 63 of file TrackUtils.cc.

Referenced by PixelTopologyBuilder::build(), TrackerGeomBuilderFromGeometricDet::build(), and pixelLocalX().

const int fireworks::BIG_PIX_PER_ROC_Y = 2
static

Definition at line 60 of file TrackUtils.cc.

Referenced by PixelTopologyBuilder::build(), TrackerGeomBuilderFromGeometricDet::build(), and pixelLocalY().

const int fireworks::COLS_PER_ROC = 52
static

Definition at line 61 of file TrackUtils.cc.

Referenced by pixelLocalY(), and RectangularPixelTopology::RectangularPixelTopology().

const double fireworks::MICRON = 1./1000./10.
static

Definition at line 57 of file TrackUtils.cc.

const TString fireworks::p1 = gSystem->Getenv("CMSSW_BASE") + datadir

Definition at line 12 of file fwPaths.cc.

Referenced by setPath().

const TString fireworks::p2 = gSystem->Getenv("CMSSW_RELEASE_BASE") + datadir

Definition at line 13 of file fwPaths.cc.

Referenced by setPath().

const double fireworks::PITCHX = 100*MICRON
static

Definition at line 58 of file TrackUtils.cc.

Referenced by pixelLocalX().

const double fireworks::PITCHY = 150*MICRON
static

Definition at line 59 of file TrackUtils.cc.

Referenced by pixelLocalY().

const int fireworks::ROWS_PER_ROC = 80
static

Definition at line 62 of file TrackUtils.cc.

Referenced by pixelLocalX(), and RectangularPixelTopology::RectangularPixelTopology().

const std::string fireworks::subdets[7] = { "UNKNOWN", "PXB", "PXF", "TIB", "TID", "TOB", "TEC" }
static

Definition at line 66 of file TrackUtils.cc.

Referenced by pushPixelHits().

Generated for CMSSW Reference Manual by doxygen 1.8.5