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 148 of file BuilderUtils.cc.

References FWProxyBuilderBase::setupAddElement().

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

149  {
150  TEveBox* eveBox = new TEveBox( "Box" );
151  eveBox->SetDrawFrame( false );
152  eveBox->SetPickable( true );
153  eveBox->SetVertices( &corners[0] );
154 
155  pb->setupAddElement( eveBox, comp );
156  }
FWProxyBuilderBase::setupAddElement
void setupAddElement(TEveElement *el, TEveElement *parent, bool set_color=true) const
Definition: FWProxyBuilderBase.cc:439
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 158 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().

159  {
160  TEveStraightLineSet* container = new TEveStraightLineSet;
161 
162  for( unsigned int iphi = 0; iphi < nLineSegments; ++iphi )
163  {
164  container->AddLine( eta + radius * cos( 2 * M_PI / nLineSegments * iphi ),
165  phi + radius * sin( 2 * M_PI / nLineSegments * iphi ),
166  0.01,
167  eta + radius * cos( 2 * M_PI / nLineSegments * ( iphi + 1 )),
168  phi + radius * sin( 2 * M_PI / nLineSegments * ( iphi + 1 )),
169  0.01 );
170  }
171  pb->setupAddElement( container, comp );
172  }
FWProxyBuilderBase::setupAddElement
void setupAddElement(TEveElement *el, TEveElement *parent, bool set_color=true) const
Definition: FWProxyBuilderBase.cc:439
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 174 of file BuilderUtils.cc.

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

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

175  {
176  TEveScalableStraightLineSet* marker = new TEveScalableStraightLineSet;
177  marker->SetLineWidth( 1 );
178  marker->SetLineStyle( 2 );
179  marker->AddLine( 0, 0, 0, size * cos( phi ), size * sin( phi ), 0 );
180  marker->AddLine( size * 0.9 * cos( phi + 0.03 ), size * 0.9 * sin( phi + 0.03 ), 0, size * cos( phi ), size * sin( phi ), 0 );
181  marker->AddLine( size * 0.9 * cos( phi - 0.03 ), size * 0.9 * sin( phi - 0.03 ), 0, size * cos( phi ), size * sin( phi ), 0 );
182  pb->setupAddElement( marker, comp );
183  }
FWProxyBuilderBase::setupAddElement
void setupAddElement(TEveElement *el, TEveElement *parent, bool set_color=true) const
Definition: FWProxyBuilderBase.cc:439
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 185 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().

186  {
187  double r( 0 );
188  if( theta < pb->context().caloTransAngle() || M_PI - theta < pb->context().caloTransAngle())
189  r = pb->context().caloZ2() / fabs( cos( theta ));
190  else
191  r = pb->context().caloR1() / sin( theta );
192 
193  TEveStraightLineSet* marker = new TEveStraightLineSet;
194  marker->SetLineWidth( 2 );
195  marker->SetLineStyle( 2 );
196  marker->AddLine( r * cos( phi ) * sin( theta ), r * sin( phi ) * sin( theta ), r * cos( theta ),
197  ( r + size ) * cos( phi ) * sin( theta ), ( r + size ) * sin( phi ) * sin( theta ), ( r + size ) * cos( theta ));
198  pb->setupAddElement( marker, comp );
199  }
FWProxyBuilderBase::context
const fireworks::Context & context() const
Definition: FWProxyBuilderBase.cc:513
FWProxyBuilderBase::setupAddElement
void setupAddElement(TEveElement *el, TEveElement *parent, bool set_color=true) const
Definition: FWProxyBuilderBase.cc:439
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:225
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:210
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 201 of file BuilderUtils.cc.

References M_PI, and FWProxyBuilderBase::setupAddElement().

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

202  {
203  TEveStraightLineSet* mainLine = new TEveStraightLineSet;
204  mainLine->AddLine( -5.191, phi, 0.01, 5.191, phi, 0.01 );
205  pb->setupAddElement( mainLine, comp );
206 
207  phi = phi > 0 ? phi - M_PI : phi + M_PI;
208  TEveStraightLineSet* secondLine = new TEveStraightLineSet;
209  secondLine->SetLineStyle( 7 );
210  secondLine->AddLine( -5.191, phi, 0.01, 5.191, phi, 0.01 );
211  pb->setupAddElement( secondLine, comp );
212  }
FWProxyBuilderBase::setupAddElement
void setupAddElement(TEveElement *el, TEveElement *parent, bool set_color=true) const
Definition: FWProxyBuilderBase.cc:439
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:439
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 385 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, edm::EventBase::get(), FWEventItem::getEvent(), FWEventItem::getGeom(), FWGeometry::getParameters(), edm::Ref< C, T, F >::id(), edm::Ref< C, T, F >::isAvailable(), edm::Ref< C, T, F >::isNonnull(), TrackingRecHit::isValid(), fwlog::kDebug, fwlog::kError, localSiStrip(), FWGeometry::localToGlobal(), edm::Handle< T >::product(), DetId::rawId(), reco::Track::recHitsBegin(), reco::Track::recHitsEnd(), and setupAddElement().

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

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

References addBox(), and energyScaledBox3DCorners().

242  {
243  std::vector<float> scaledCorners( 24 );
244  energyScaledBox3DCorners(corners, scale, scaledCorners, invert);
245  addBox( scaledCorners, comp, pb );
246  }
AlCaHLTBitMon_QueryRunRegistry.comp
string comp
Definition: AlCaHLTBitMon_QueryRunRegistry.py:248
HLT_25ns14e33_v1_cff.invert
tuple invert
Definition: HLT_25ns14e33_v1_cff.py:22882
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:215
fireworks::addBox
void addBox(const std::vector< float > &corners, TEveElement *, FWProxyBuilderBase *)
Definition: BuilderUtils.cc:148
void fireworks::drawEnergyTower3D ( const float *  corners,
float  scale,
TEveElement *  comp,
FWProxyBuilderBase *  pb,
bool  reflect = false 
)

Definition at line 325 of file BuilderUtils.cc.

References addBox(), and energyTower3DCorners().

326  {
327  std::vector<float> scaledCorners( 24 );
328  energyTower3DCorners(corners, scale, scaledCorners, reflect);
329  addBox( scaledCorners, comp, pb );
330  }
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:286
fireworks::addBox
void addBox(const std::vector< float > &corners, TEveElement *, FWProxyBuilderBase *)
Definition: BuilderUtils.cc:148
void fireworks::drawEtScaledBox3D ( const float *  corners,
float  energy,
float  maxEnergy,
TEveElement *  comp,
FWProxyBuilderBase *  pb,
bool  reflect = false 
)

Definition at line 278 of file BuilderUtils.cc.

References addBox(), and etScaledBox3DCorners().

279  {
280  std::vector<float> scaledCorners( 24 );
281  etScaledBox3DCorners(corners, energy, maxEnergy, scaledCorners, invert);
282  addBox( scaledCorners, comp, pb );
283  }
AlCaHLTBitMon_QueryRunRegistry.comp
string comp
Definition: AlCaHLTBitMon_QueryRunRegistry.py:248
HLT_25ns14e33_v1_cff.invert
tuple invert
Definition: HLT_25ns14e33_v1_cff.py:22882
fireworks::addBox
void addBox(const std::vector< float > &corners, TEveElement *, FWProxyBuilderBase *)
Definition: BuilderUtils.cc:148
fireworks::etScaledBox3DCorners
void etScaledBox3DCorners(const float *corners, float energy, float maxEnergy, std::vector< float > &scaledCorners, bool reflect=false)
Definition: BuilderUtils.cc:249
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 374 of file BuilderUtils.cc.

References addBox(), and etTower3DCorners().

375  {
376  std::vector<float> scaledCorners( 24 );
377  etTower3DCorners(corners, scale, scaledCorners, reflect);
378  addBox( scaledCorners, comp, pb );
379  }
fireworks::etTower3DCorners
void etTower3DCorners(const float *corners, float scale, std::vector< float > &, bool reflect=false)
Definition: BuilderUtils.cc:334
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:148
void fireworks::energyScaledBox3DCorners ( const float *  corners,
float  scale,
std::vector< float > &  scaledCorners,
bool  invert = false 
)

Definition at line 215 of file BuilderUtils.cc.

References f, i, and invertBox().

Referenced by drawEnergyScaledBox3D().

216  {
217  std::vector<float> centre( 3, 0 );
218 
219  for( unsigned int i = 0; i < 24; i += 3 )
220  {
221  centre[0] += corners[i];
222  centre[1] += corners[i + 1];
223  centre[2] += corners[i + 2];
224  }
225 
226  for( unsigned int i = 0; i < 3; ++i )
227  centre[i] *= 1.0f / 8.0f;
228 
229  // Coordinates for a scaled version of the original box
230  for( unsigned int i = 0; i < 24; i += 3 )
231  {
232  scaledCorners[i] = centre[0] + ( corners[i] - centre[0] ) * scale;
233  scaledCorners[i + 1] = centre[1] + ( corners[i + 1] - centre[1] ) * scale;
234  scaledCorners[i + 2] = centre[2] + ( corners[i + 2] - centre[2] ) * scale;
235  }
236 
237  if( invert )
238  invertBox( scaledCorners );
239  }
i
int i
Definition: DBlmapReader.cc:9
HLT_25ns14e33_v1_cff.invert
tuple invert
Definition: HLT_25ns14e33_v1_cff.py:22882
pileupReCalc_HLTpaths.scale
int scale
Definition: pileupReCalc_HLTpaths.py:111
fireworks::invertBox
void invertBox(std::vector< float > &corners)
Definition: BuilderUtils.cc:129
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 286 of file BuilderUtils.cc.

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

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

287  {
288  for( int i = 0; i < 24; ++i )
289  scaledCorners[i] = corners[i];
290  // Coordinates of a front face scaled
291  if( reflect )
292  {
293  // We know, that an ES rechit geometry in -Z needs correction.
294  // The back face is actually its front face.
295  for( unsigned int i = 0; i < 12; i += 3 )
296  {
297  TEveVector diff( corners[i] - corners[i + 12], corners[i + 1] - corners[i + 13], corners[i + 2] - corners[i + 14] );
298  diff.Normalize();
299  diff *= scale;
300 
301  scaledCorners[i] = corners[i] + diff.fX;
302  scaledCorners[i + 1] = corners[i + 1] + diff.fY;
303  scaledCorners[i + 2] = corners[i + 2] + diff.fZ;
304  }
305  }
306  else
307  {
308  for( unsigned int i = 0; i < 12; i += 3 )
309  {
310  TEveVector diff( corners[i + 12] - corners[i], corners[i + 13] - corners[i + 1], corners[i + 14] - corners[i + 2] );
311  diff.Normalize();
312  diff *= scale;
313 
314  scaledCorners[i] = corners[i + 12];
315  scaledCorners[i + 1] = corners[i + 13];
316  scaledCorners[i + 2] = corners[i + 14];
317 
318  scaledCorners[i + 12] = corners[i + 12] + diff.fX;
319  scaledCorners[i + 13] = corners[i + 13] + diff.fY;
320  scaledCorners[i + 14] = corners[i + 14] + diff.fZ;
321  }
322  }
323  }
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 249 of file BuilderUtils.cc.

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

Referenced by drawEtScaledBox3D().

250  {
251  std::vector<float> centre( 3, 0 );
252 
253  for( unsigned int i = 0; i < 24; i += 3 )
254  {
255  centre[0] += corners[i];
256  centre[1] += corners[i + 1];
257  centre[2] += corners[i + 2];
258  }
259 
260  for( unsigned int i = 0; i < 3; ++i )
261  centre[i] *= 1.0f / 8.0f;
262 
263  TEveVector c( centre[0], centre[1], centre[2] );
264  float scale = energy / maxEnergy * sin( c.Theta());
265 
266  // Coordinates for a scaled version of the original box
267  for( unsigned int i = 0; i < 24; i += 3 )
268  {
269  scaledCorners[i] = centre[0] + ( corners[i] - centre[0] ) * scale;
270  scaledCorners[i + 1] = centre[1] + ( corners[i + 1] - centre[1] ) * scale;
271  scaledCorners[i + 2] = centre[2] + ( corners[i + 2] - centre[2] ) * scale;
272  }
273 
274  if( invert )
275  invertBox( scaledCorners );
276  }
i
int i
Definition: DBlmapReader.cc:9
funct::sin
Sin< T >::type sin(const T &t)
Definition: Sin.h:22
HLT_25ns14e33_v1_cff.invert
tuple invert
Definition: HLT_25ns14e33_v1_cff.py:22882
pileupReCalc_HLTpaths.scale
int scale
Definition: pileupReCalc_HLTpaths.py:111
fireworks::invertBox
void invertBox(std::vector< float > &corners)
Definition: BuilderUtils.cc:129
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 334 of file BuilderUtils.cc.

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

Referenced by drawEtTower3D().

335  {
336  for( int i = 0; i < 24; ++i )
337  scaledCorners[i] = corners[i];
338  // Coordinates of a front face scaled
339  if( reflect )
340  {
341  // We know, that an ES rechit geometry in -Z needs correction.
342  // The back face is actually its front face.
343  for( unsigned int i = 0; i < 12; i += 3 )
344  {
345  TEveVector diff( corners[i] - corners[i + 12], corners[i + 1] - corners[i + 13], corners[i + 2] - corners[i + 14] );
346  diff.Normalize();
347  diff *= ( scale * sin( diff.Theta()));
348 
349  scaledCorners[i] = corners[i] + diff.fX;
350  scaledCorners[i + 1] = corners[i + 1] + diff.fY;
351  scaledCorners[i + 2] = corners[i + 2] + diff.fZ;
352  }
353  }
354  else
355  {
356  for( unsigned int i = 0; i < 12; i += 3 )
357  {
358  TEveVector diff( corners[i + 12] - corners[i], corners[i + 13] - corners[i + 1], corners[i + 14] - corners[i + 2] );
359  diff.Normalize();
360  diff *= ( scale * sin( diff.Theta()));
361 
362  scaledCorners[i] = corners[i + 12];
363  scaledCorners[i + 1] = corners[i + 13];
364  scaledCorners[i + 2] = corners[i + 14];
365 
366  scaledCorners[i + 12] = corners[i + 12] + diff.fX;
367  scaledCorners[i + 13] = corners[i + 13] + diff.fY;
368  scaledCorners[i + 14] = corners[i + 14] + diff.fZ;
369  }
370  }
371  }
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 362 of file TrackUtils.cc.

References fwLog, and fwlog::kDebug.

Referenced by addSiStripClusters(), and FWTrajectorySeedProxyBuilder::build().

363 {
364  const SiStripCluster* cluster = 0;
365 
366  if( const SiStripRecHit2D* hit2D = dynamic_cast<const SiStripRecHit2D*>( rechit ))
367  {
368  fwLog( fwlog::kDebug ) << "hit 2D ";
369 
370  cluster = hit2D->cluster().get();
371  }
372  if( cluster == 0 )
373  {
374  if( const SiStripRecHit1D* hit1D = dynamic_cast<const SiStripRecHit1D*>( rechit ))
375  {
376  fwLog( fwlog::kDebug ) << "hit 1D ";
377 
378  cluster = hit1D->cluster().get();
379  }
380  }
381  return cluster;
382 }
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 
116  time_t t( event.time().value() >> 32 );
117  struct tm * xx = localtime( &t );
118  std::string text( asctime(xx) );
119  size_t pos = text.find('\n');
120  if( pos != std::string::npos ) text = text.substr( 0, pos );
121  text += " ";
122  if( xx->tm_isdst )
123  text += tzname[1];
124  else
125  text += tzname[0];
126  return text;
127  }
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:61
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 EnergyCorrector::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:61
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 690 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().

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

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

874  {
875  std::string text;
876  for(std::set<DetId>::const_iterator id = idSet.begin(), idEnd = idSet.end(); id != idEnd; ++id)
877  {
878  text += info(*id);
879  text += "\n";
880  }
881  return text;
882 }
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 885 of file TrackUtils.cc.

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

885  {
886  std::string text;
887  for(std::vector<DetId>::const_iterator id = idSet.begin(), idEnd = idSet.end(); id != idEnd; ++id)
888  {
889  text += info(*id);
890  text += "\n";
891  }
892  return text;
893 }
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 129 of file BuilderUtils.cc.

References std::swap().

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

130  {
131  std::swap( corners[0], corners[9] );
132  std::swap( corners[1], corners[10] );
133  std::swap( corners[2], corners[11] );
134 
135  std::swap( corners[3], corners[6] );
136  std::swap( corners[4], corners[7] );
137  std::swap( corners[5], corners[8] );
138 
139  std::swap( corners[12], corners[21] );
140  std::swap( corners[13], corners[22] );
141  std::swap( corners[14], corners[23] );
142 
143  std::swap( corners[15], corners[18] );
144  std::swap( corners[16], corners[19] );
145  std::swap( corners[17], corners[20] );
146  }
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 295 of file TrackUtils.cc.

References delta, fwLog, fwlog::kError, hltrates_dqm_sourceclient-live_cfg::offset, funct::tan(), and ecaldqm::topology().

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

296 {
297  Float_t topology = pars[0];
298  Float_t halfStripLength = pars[2] * 0.5;
299 
300  Double_t localCenter[3] = { 0.0, 0.0, 0.0 };
301  localTop[1] = halfStripLength;
302  localBottom[1] = -halfStripLength;
303 
304  if( topology == 1 ) // RadialStripTopology
305  {
306  // stripAngle = phiOfOneEdge + strip * angularWidth
307  // localY = originToIntersection * tan( stripAngle )
308  Float_t stripAngle = tan( pars[5] + strip * pars[6] );
309  Float_t delta = halfStripLength * stripAngle;
310  localCenter[0] = pars[4] * stripAngle;
311  localTop[0] = localCenter[0] + delta;
312  localBottom[0] = localCenter[0] - delta;
313  }
314  else if( topology == 2 ) // RectangularStripTopology
315  {
316  // offset = -numberOfStrips/2. * pitch
317  // localY = strip * pitch + offset
318  Float_t offset = -pars[1] * 0.5 * pars[3];
319  localCenter[0] = strip * pars[3] + offset;
320  localTop[0] = localCenter[0];
321  localBottom[0] = localCenter[0];
322  }
323  else if( topology == 3 ) // TrapezoidalStripTopology
324  {
325  fwLog( fwlog::kError )
326  << "did not expect TrapezoidalStripTopology of "
327  << id << std::endl;
328  }
329  else if( pars[0] == 0 ) // StripTopology
330  {
331  fwLog( fwlog::kError )
332  << "did not find StripTopology of "
333  << id << std::endl;
334  }
335 }
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
fwLog
#define fwLog(_level_)
Definition: fwLog.h:50
hltrates_dqm_sourceclient-live_cfg.offset
tuple offset
Definition: hltrates_dqm_sourceclient-live_cfg.py:81
TEveJetCone * fireworks::makeEveJetCone ( const reco::Jet &  iData,
const fireworks::Context &  context 
)

Definition at line 7 of file makeEveJetCone.cc.

References EnergyCorrector::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:235
reco::LeafCandidate::vertex
virtual const Point & vertex() const
vertex position (overwritten by PF...)
Definition: LeafCandidate.h:150
fireworks::Context::caloZ2
static float caloZ2(bool offset=true)
Definition: Context.cc:225
reco::Jet::Constituents
std::vector< Constituent > Constituents
Definition: Jet.h:23
reco::Jet::etaetaMoment
float etaetaMoment() const
eta-eta second moment, ET weighted
reco::Jet::getJetConstituents
virtual Constituents getJetConstituents() const
list of constituents
reco::LeafCandidate::eta
virtual double eta() const
momentum pseudorapidity
Definition: LeafCandidate.h:137
fireworks::Context::caloR2
static float caloR2(bool offset=true)
Definition: Context.cc:215
mathSSE::sqrt
T sqrt(T t)
Definition: SSEVec.h:48
fireworks::Context::caloR1
static float caloR1(bool offset=true)
Definition: Context.cc:210
funct::tan
Tan< T >::type tan(const T &t)
Definition: Tan.h:22
reco::Jet::phiphiMoment
float phiphiMoment() const
phi-phi second moment, ET weighted
fireworks::Context::caloMaxEta
static double caloMaxEta()
Definition: Context.cc:240
fireworks::Context::caloZ1
static float caloZ1(bool offset=true)
Definition: Context.cc:220
reco::LeafCandidate::phi
virtual double phi() const
momentum azimuthal angle
Definition: LeafCandidate.h:133
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  if (detids.empty()) return false;
38  std::vector<double> phis;
39  for( std::vector<std::pair<DetId, float> >::const_iterator id = detids.begin(), end = detids.end(); id != end; ++id )
40  {
41  const float* corners = pb->context().getGeom()->getCorners( id->first.rawId());
42  if( corners != 0 )
43  {
44  std::vector<float> centre( 3, 0 );
45 
46  for( unsigned int i = 0; i < 24; i += 3 )
47  {
48  centre[0] += corners[i];
49  centre[1] += corners[i + 1];
50  centre[2] += corners[i + 2];
51  }
52 
53  phis.push_back( TEveVector( centre[0], centre[1], centre[2] ).Phi());
54  }
55  }
56  std::pair<double,double> phiRange = fireworks::getPhiRange( phis, iPhi );
57  const double r = pb->context().caloR1();
58  TGeoBBox *sc_box = new TGeoTubeSeg( r - 2, r , 1,
59  phiRange.first * 180 / M_PI - 0.5,
60  phiRange.second * 180 / M_PI + 0.5 ); // 0.5 is roughly half size of a crystal
61  TEveGeoShape *sc = fireworks::getShape( "supercluster", sc_box, pb->item()->defaultDisplayProperties().color());
62  sc->SetPickable( kTRUE );
63  pb->setupAddElement( sc, &oItemHolder );
64  return true;
65 }
edm::Ref::isAvailable
bool isAvailable() const
Definition: Ref.h:614
FWProxyBuilderBase::context
const fireworks::Context & context() const
Definition: FWProxyBuilderBase.cc:513
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:439
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:210
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:280
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 67 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().

71 {
72  if( !iCluster.isAvailable()) return false;
73  TEveGeoManagerHolder gmgr( TEveGeoShape::GetGeoMangeur());
74  double theta_max = 0;
75  double theta_min = 10;
76  std::vector<std::pair<DetId, float> > detids = iCluster->hitsAndFractions();
77  if (detids.empty()) return false;
78  for( std::vector<std::pair<DetId, float> >::const_iterator id = detids.begin(), end = detids.end(); id != end; ++id )
79  {
80  const float* corners = pb->context().getGeom()->getCorners( id->first.rawId());
81  if( corners != 0 )
82  {
83  std::vector<float> centre( 3, 0 );
84 
85  for( unsigned int i = 0; i < 24; i += 3 )
86  {
87  centre[0] += corners[i];
88  centre[1] += corners[i + 1];
89  centre[2] += corners[i + 2];
90  }
91 
92  double theta = TEveVector( centre[0], centre[1], centre[2] ).Theta();
93  if( theta > theta_max ) theta_max = theta;
94  if( theta < theta_min ) theta_min = theta;
95  }
96  }
97  // expand theta range by the size of a crystal to avoid segments of zero length
98  bool barrel = true;
99  if ((theta_max > 0 && theta_max < pb->context().caloTransAngle()) ||
100  ( theta_min > (TMath::Pi() -pb->context().caloTransAngle())) )
101  {
102  barrel = false;
103  }
104 
105  double z_ecal = barrel ? pb->context().caloZ1() : pb->context().caloZ2();
106  double r_ecal = barrel ? pb->context().caloR1() : pb->context().caloR2();
107 
108  fireworks::addRhoZEnergyProjection( pb, &oItemHolder, r_ecal-1, z_ecal-1,
109  theta_min - 0.003, theta_max + 0.003,
110  iPhi );
111 
112  return true;
113 }
edm::Ref::isAvailable
bool isAvailable() const
Definition: Ref.h:614
Pi
const double Pi
Definition: CosmicMuonParameters.h:18
FWProxyBuilderBase::context
const fireworks::Context & context() const
Definition: FWProxyBuilderBase.cc:513
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:235
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:225
fireworks::Context::caloR2
static float caloR2(bool offset=true)
Definition: Context.cc:215
fireworks::Context::caloR1
static float caloR1(bool offset=true)
Definition: Context.cc:210
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:280
fireworks::Context::caloZ1
static float caloZ1(bool offset=true)
Definition: Context.cc:220
float fireworks::pixelLocalX ( const double  mpx,
const int  m_nrows 
)

Definition at line 168 of file TrackUtils.cc.

References BIG_PIX_PER_ROC_X, PITCHX, and ROWS_PER_ROC.

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

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

Definition at line 201 of file TrackUtils.cc.

References BIG_PIX_PER_ROC_Y, COLS_PER_ROC, and PITCHY.

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

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

Definition at line 73 of file TrackUtils.cc.

References assert(), 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().

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

Definition at line 551 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, edm::EventBase::get(), FWEventItem::getEvent(), FWEventItem::getGeom(), FWGeometry::getParameters(), edm::Ref< C, T, F >::id(), edm::Ref< C, T, F >::isAvailable(), edm::Ref< C, T, F >::isNonnull(), fwlog::kError, PixelSubdetector::PixelBarrel, PixelSubdetector::PixelEndcap, edm::Handle< T >::product(), pushPixelCluster(), reco::Track::recHitsBegin(), and reco::Track::recHitsEnd().

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

552 {
553  const edmNew::DetSetVector<SiPixelCluster> * allClusters = 0;
554  for( trackingRecHit_iterator it = t.recHitsBegin(), itEnd = t.recHitsEnd(); it != itEnd; ++it)
555  {
556  if( typeid(**it) == typeid( SiPixelRecHit ))
557  {
558  const SiPixelRecHit &hit = static_cast<const SiPixelRecHit &>(**it);
559  if( hit.cluster().isNonnull() && hit.cluster().isAvailable())
560  {
561  edm::Handle<edmNew::DetSetVector<SiPixelCluster> > allClustersHandle;
562  iItem.getEvent()->get(hit.cluster().id(), allClustersHandle);
563  allClusters = allClustersHandle.product();
564  break;
565  }
566  }
567  }
568  if( allClusters == 0 ) return;
569 
570  const FWGeometry *geom = iItem.getGeom();
571 
572  for( trackingRecHit_iterator it = t.recHitsBegin(), itEnd = t.recHitsEnd(); it != itEnd; ++it )
573  {
574  const TrackingRecHit* rh = &(**it);
575 
576  DetId id = (*it)->geographicalId();
577  if( ! geom->contains( id ))
578  {
579  fwLog( fwlog::kError )
580  << "failed to get geometry of Tracker Det with raw id: "
581  << id.rawId() << std::endl;
582 
583  continue;
584  }
585 
586  // -- in which detector are we?
587  unsigned int subdet = (unsigned int)id.subdetId();
588  if(( subdet != PixelSubdetector::PixelBarrel ) && ( subdet != PixelSubdetector::PixelEndcap )) continue;
589 
590  const SiPixelCluster* hitCluster = 0;
591  if( const SiPixelRecHit* pixel = dynamic_cast<const SiPixelRecHit*>( rh ))
592  hitCluster = pixel->cluster().get();
593  edmNew::DetSetVector<SiPixelCluster>::const_iterator itds = allClusters->find(id.rawId());
594  if( itds != allClusters->end())
595  {
596  const edmNew::DetSet<SiPixelCluster> & clustersOnThisDet = *itds;
597  for( edmNew::DetSet<SiPixelCluster>::const_iterator itc = clustersOnThisDet.begin(), edc = clustersOnThisDet.end(); itc != edc; ++itc )
598  {
599  if( &*itc != hitCluster )
600  pushPixelCluster( pixelPoints, *geom, id, *itc, geom->getParameters( id ));
601  }
602  }
603  }
604 }
edm::Ref::isAvailable
bool isAvailable() const
Definition: Ref.h:614
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:339
edm::Ref::isNonnull
bool isNonnull() const
Checks for non-null.
Definition: Ref.h:250
FWGeometry::getParameters
const float * getParameters(unsigned int id) const
Definition: FWGeometry.cc:296
edm::EventBase::get
bool get(ProductID const &, Handle< T > &) const
Definition: EventBase.h:104
edmNew::DetSet::const_iterator
data_type const * const_iterator
Definition: DetSetNew.h:30
edm::Ref::id
ProductID id() const
Accessor for product ID.
Definition: Ref.h:256
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:104
FWEventItem::getEvent
const edm::EventBase * getEvent() const
Definition: FWEventItem.h:146
DetId
Definition: DetId.h:18
edm::Handle::product
T const * product() const
Definition: Handle.h:81
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:318
SiPixelCluster
Pixel cluster – collection of neighboring pixels above threshold.
Definition: SiPixelCluster.h:33
edmNew::DetSet::end
iterator end()
Definition: DetSetNew.h:70
fireworks::pushPixelCluster
void pushPixelCluster(std::vector< TVector3 > &pixelPoints, const FWGeometry &geom, DetId id, const SiPixelCluster &c, const float *pars)
Definition: TrackUtils.cc:656
FWEventItem::getGeom
const FWGeometry * getGeom() const
Definition: FWEventItem.cc:681
trackingRecHit_iterator
TrackingRecHitCollection::base::const_iterator trackingRecHit_iterator
iterator over a vector of reference to TrackingRecHit in the same collection
Definition: TrackingRecHitFwd.h:19
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:109
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 656 of file TrackUtils.cc.

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

Referenced by pushNearbyPixelHits(), and pushPixelHits().

657 {
658  double row = c.minPixelRow();
659  double col = c.minPixelCol();
660  float lx = 0.;
661  float ly = 0.;
662 
663  int nrows = (int)pars[0];
664  int ncols = (int)pars[1];
665  lx = pixelLocalX( row, nrows );
666  ly = pixelLocalY( col, ncols );
667 
668  fwLog( fwlog::kDebug )
669  << ", row: " << row << ", col: " << col
670  << ", lx: " << lx << ", ly: " << ly ;
671 
672  float local[3] = { lx, ly, 0. };
673  float global[3];
674  geom.localToGlobal( id, local, global );
675  TVector3 pb( global[0], global[1], global[2] );
676  pixelPoints.push_back( pb );
677 
678  fwLog( fwlog::kDebug )
679  << " x: " << pb.X()
680  << ", y: " << pb.Y()
681  << " z: " << pb.Z()
682  << " eta: " << pb.Eta()
683  << ", phi: " << pb.Phi()
684  << " rho: " << pb.Pt() << std::endl;
685 }
SiPixelCluster::minPixelCol
int minPixelCol() const
Definition: SiPixelCluster.h:155
fireworks::pixelLocalY
float pixelLocalY(const double mpy, const int m_ncols)
Definition: TrackUtils.cc:201
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:328
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:168
void fireworks::pushPixelHits ( std::vector< TVector3 > &  pixelPoints,
const FWEventItem &  iItem,
const reco::Track &  t 
)

Definition at line 609 of file TrackUtils.cc.

References EnergyCorrector::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().

610 {
611  /*
612  * -- return for each Pixel Hit a 3D point
613  */
614  const FWGeometry *geom = iItem.getGeom();
615 
616  double dz = t.dz();
617  double vz = t.vz();
618  double etaT = t.eta();
619 
620  fwLog( fwlog::kDebug ) << "Track eta: " << etaT << ", vz: " << vz << ", dz: " << dz
621  << std::endl;
622 
623  int cnt = 0;
624  for( trackingRecHit_iterator it = t.recHitsBegin(), itEnd = t.recHitsEnd(); it != itEnd; ++it )
625  {
626  const TrackingRecHit* rh = &(**it);
627  // -- get position of center of wafer, assuming (0,0,0) is the center
628  DetId id = (*it)->geographicalId();
629  if( ! geom->contains( id ))
630  {
631  fwLog( fwlog::kError )
632  << "failed to get geometry of Tracker Det with raw id: "
633  << id.rawId() << std::endl;
634 
635  continue;
636  }
637 
638  // -- in which detector are we?
639  unsigned int subdet = (unsigned int)id.subdetId();
640 
641  if(( subdet == PixelSubdetector::PixelBarrel ) || ( subdet == PixelSubdetector::PixelEndcap ))
642  {
643  fwLog( fwlog::kDebug ) << cnt++ << " -- "
644  << subdets[subdet];
645 
646  if( const SiPixelRecHit* pixel = dynamic_cast<const SiPixelRecHit*>( rh ))
647  {
648  const SiPixelCluster& c = *( pixel->cluster());
649  pushPixelCluster( pixelPoints, *geom, id, c, geom->getParameters( id ));
650  }
651  }
652  }
653 }
FWGeometry::getParameters
const float * getParameters(unsigned int id) const
Definition: FWGeometry.cc:296
reco::TrackBase::eta
double eta() const
pseudorapidity of momentum vector
Definition: TrackBase.h:604
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:104
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:562
reco::TrackBase::vz
double vz() const
z coordinate of the reference point on track
Definition: TrackBase.h:622
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:656
FWEventItem::getGeom
const FWGeometry * getGeom() const
Definition: FWEventItem.cc:681
trackingRecHit_iterator
TrackingRecHitCollection::base::const_iterator trackingRecHit_iterator
iterator over a vector of reference to TrackingRecHit in the same collection
Definition: TrackingRecHitFwd.h:19
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:109
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(), 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:106
void fireworks::setupAddElement ( TEveElement *  el,
TEveElement *  parent,
const FWEventItem *  item,
bool  master,
bool  color 
)

Definition at line 340 of file TrackUtils.cc.

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

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

341 {
342  if (master)
343  {
344  el->CSCTakeAnyParentAsMaster();
345  el->SetPickable(true);
346  }
347 
348  if (color)
349  {
350  el->CSCApplyMainColorToMatchingChildren();
351  el->CSCApplyMainTransparencyToMatchingChildren();
352  el->SetMainColor(item->defaultDisplayProperties().color());
353  assert((item->defaultDisplayProperties().transparency() >= 0)
354  && (item->defaultDisplayProperties().transparency() <= 100));
355  el->SetMainTransparency(item->defaultDisplayProperties().transparency());
356  }
357  parent->AddElement(el);
358 }
FWEventItem::defaultDisplayProperties
const FWDisplayProperties & defaultDisplayProperties() const
Definition: FWEventItem.cc:451
dbtoconf.parent
list parent
Definition: dbtoconf.py:74
assert
assert(m_qm.get())
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 67 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 64 of file TrackUtils.cc.

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

const int fireworks::COLS_PER_ROC = 52
static

Definition at line 65 of file TrackUtils.cc.

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

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

Definition at line 61 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 62 of file TrackUtils.cc.

Referenced by pixelLocalX().

const double fireworks::PITCHY = 150*MICRON
static

Definition at line 63 of file TrackUtils.cc.

Referenced by pixelLocalY().

const int fireworks::ROWS_PER_ROC = 80
static

Definition at line 66 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 70 of file TrackUtils.cc.

Referenced by pushPixelHits().

Generated for CMSSW Reference Manual by doxygen 1.8.5