CMS 3D CMS Logo

 All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Properties Friends Macros Pages
List of all members | Public Member Functions | Static Public Member Functions | Private Member Functions | Private Attributes
FWPFCandidateWithHitsProxyBuilder Class Reference

#include <FWPFCandidateWithHitsProxyBuilder.h>

Inheritance diagram for FWPFCandidateWithHitsProxyBuilder:
FWProxyBuilderBase

Public Member Functions

virtual void build (const FWEventItem *iItem, TEveElementList *product, const FWViewContext *)
 
 FWPFCandidateWithHitsProxyBuilder ()
 
virtual bool havePerViewProduct (FWViewType::EType) const
 
const std::string & purpose () const
 
virtual void scaleProduct (TEveElementList *parent, FWViewType::EType type, const FWViewContext *vc)
 
const std::string & typeName () const
 
const std::string & view () const
 
virtual ~FWPFCandidateWithHitsProxyBuilder ()
 
- Public Member Functions inherited from FWProxyBuilderBase
void build ()
 
virtual bool canHandle (const FWEventItem &)
 
const fireworks::Contextcontext () const
 
TEveElementList * createProduct (FWViewType::EType, const FWViewContext *)
 
 FWProxyBuilderBase ()
 
bool getHaveWindow () const
 
virtual bool haveSingleProduct () const
 
const FWEventItemitem () const
 
virtual void itemBeingDestroyed (const FWEventItem *)
 
void itemChanged (const FWEventItem *)
 
int layer () const
 
void modelChanges (const FWModelIds &)
 
void removePerViewProduct (FWViewType::EType, const FWViewContext *vc)
 
void scaleChanged (const FWViewContext *)
 
void setHaveWindow (bool iFlag)
 
virtual void setInteractionList (FWInteractionList *, const std::string &)
 
virtual void setItem (const FWEventItem *iItem)
 
void setupAddElement (TEveElement *el, TEveElement *parent, bool set_color=true) const
 
void setupElement (TEveElement *el, bool color=true) const
 
virtual bool willHandleInteraction () const
 
virtual ~FWProxyBuilderBase ()
 

Static Public Member Functions

static const std::string & classPurpose ()
 
static const std::string & classRegisterTypeName ()
 
static const std::string & classTypeName ()
 
static const std::string & classView ()
 
- Static Public Member Functions inherited from FWProxyBuilderBase
static bool representsSubPart ()
 
static std::string typeOfBuilder ()
 Used by the plugin system to determine how the proxy uses the data from FWEventItem. More...
 

Private Member Functions

void addHitsForCandidate (const reco::PFCandidate &c, TEveElement *holder, const FWViewContext *vc)
 
 FWPFCandidateWithHitsProxyBuilder (const FWPFCandidateWithHitsProxyBuilder &)
 
const reco::PFRecHitgetHitForDetId (unsigned detId)
 
void initPFRecHitsCollections ()
 
const
FWPFCandidateWithHitsProxyBuilder
operator= (const FWPFCandidateWithHitsProxyBuilder &)
 
void viewContextBoxScale (const float *corners, float scale, bool plotEt, std::vector< float > &scaledCorners, const reco::PFRecHit *)
 

Private Attributes

const reco::PFRecHitCollectionm_collectionHCAL
 

Additional Inherited Members

- Protected Types inherited from FWProxyBuilderBase
typedef std::vector< Product * >
::iterator 
Product_it
 
- Protected Member Functions inherited from FWProxyBuilderBase
virtual void buildViewType (const FWEventItem *iItem, TEveElementList *, FWViewType::EType, const FWViewContext *)
 
virtual void clean ()
 
virtual void cleanLocal ()
 
TEveCompound * createCompound (bool set_color=true, bool propagate_color_to_all_children=false) const
 
 FWProxyBuilderBase (const FWProxyBuilderBase &)
 
void increaseComponentTransparency (unsigned int index, TEveElement *holder, const std::string &name, Char_t transpOffset)
 
virtual void localModelChanges (const FWModelId &iId, TEveElement *iCompound, FWViewType::EType viewType, const FWViewContext *vc)
 
virtual void modelChanges (const FWModelIds &, Product *)
 
const FWProxyBuilderBaseoperator= (const FWProxyBuilderBase &)
 
virtual bool visibilityModelChanges (const FWModelId &, TEveElement *, FWViewType::EType, const FWViewContext *)
 
- Protected Attributes inherited from FWProxyBuilderBase
std::vector< Product * > m_products
 

Detailed Description

Definition at line 12 of file FWPFCandidateWithHitsProxyBuilder.h.

Constructor & Destructor Documentation

FWPFCandidateWithHitsProxyBuilder::FWPFCandidateWithHitsProxyBuilder ( )
inline

Definition at line 15 of file FWPFCandidateWithHitsProxyBuilder.h.

15 {}
virtual FWPFCandidateWithHitsProxyBuilder::~FWPFCandidateWithHitsProxyBuilder ( )
inlinevirtual

Definition at line 16 of file FWPFCandidateWithHitsProxyBuilder.h.

16 {}
FWPFCandidateWithHitsProxyBuilder::FWPFCandidateWithHitsProxyBuilder ( const FWPFCandidateWithHitsProxyBuilder )
private

Member Function Documentation

void FWPFCandidateWithHitsProxyBuilder::addHitsForCandidate ( const reco::PFCandidate c,
TEveElement *  holder,
const FWViewContext vc 
)
private

Definition at line 247 of file FWPFCandidateWithHitsProxyBuilder.cc.

References FWProxyBuilderBase::context(), reco::PFCandidate::elementsInBlocks(), reco::PFRecHit::energy(), fwLog, edm::Ref< C, T, F >::get(), FWGeometry::getCorners(), FWViewContext::getEnergyScale(), fireworks::Context::getGeom(), getHitForDetId(), FWViewEnergyScale::getPlotEt(), FWViewEnergyScale::getScaleFactor3D(), edm::Ref< C, T, F >::id(), edm::Ref< C, T, F >::key(), fwlog::kInfo, fwlog::kWarning, fireworks::lastIdx, m_collectionHCAL, pileupReCalc_HLTpaths::scale, FWProxyBuilderBase::setupAddElement(), and viewContextBoxScale().

Referenced by build().

248 {
249  reco::PFCandidate::ElementsInBlocks eleInBlocks = cand.elementsInBlocks();
250 
251  TEveBoxSet* boxset = 0;
252  TEveStraightLineSet* lineset = 0;
253 
254  for(unsigned elIdx=0; elIdx<eleInBlocks.size(); elIdx++)
255  {
256  // unsigned ieTrack = 0;
257  // unsigned ieECAL = 0;
258  unsigned ieHCAL = 0;
259 
260  reco::PFBlockRef blockRef = eleInBlocks[elIdx].first;
261  unsigned indexInBlock = eleInBlocks[elIdx].second;
262  edm::Ptr<reco::PFBlock> myBlock(blockRef.id(),blockRef.get(), blockRef.key());
263  /*
264  if (myBlock->elements()[indexInBlock].type() == 1)
265  ieTrack = indexInBlock;
266  if (myBlock->elements()[indexInBlock].type() == 4)
267  ieECAL = indexInBlock;
268  */
269  if (myBlock->elements()[indexInBlock].type() == 5)
270  ieHCAL = indexInBlock;
271 
272 
273  std::vector<float> scaledCorners(24);
274  float scale = vc->getEnergyScale()->getScaleFactor3D()/50;
275  if (ieHCAL && m_collectionHCAL) {
276  reco::PFClusterRef hcalclusterRef=myBlock->elements()[ieHCAL].clusterRef();
277  edm::Ptr<reco::PFCluster> myCluster(hcalclusterRef.id(),hcalclusterRef.get(), hcalclusterRef.key());
278  if (myCluster.get())
279  {
280  const std::vector< std::pair<DetId, float> > & hitsandfracs = myCluster->hitsAndFractions();
281 
282  if (!boxset)
283  {
284  boxset = new TEveBoxSet();
285  boxset->Reset(TEveBoxSet::kBT_FreeBox, true, hitsandfracs.size());
286  boxset->SetAntiFlick(true);
287  boxset->SetAlwaysSecSelect(1);
288  boxset->SetPickable(1);
289  boxset->SetTooltipCBFoo(boxset_tooltip_callback);
290  }
291 
292  if (!lineset)
293  {
294  lineset = new TEveStraightLineSet();
295  }
296 
297  bool hitsFound = false;
298  for ( int ihandf=0, lastIdx=(int)(hitsandfracs.size()); ihandf<lastIdx; ihandf++)
299  {
300  unsigned int hitDetId = hitsandfracs[ihandf].first;
301  const float* corners = context().getGeom()->getCorners(hitDetId);
302  const reco::PFRecHit* hit = getHitForDetId(hitDetId);
303  if (hit)
304  {
305  viewContextBoxScale( corners, hit->energy()*scale, vc->getEnergyScale()->getPlotEt(), scaledCorners, hit);
306  boxset->AddBox( &scaledCorners[0]);
307  // setup last box
308  boxset->DigitColor(holder->GetMainColor());
309  boxset->DigitUserData((void*)hit);
310  addBoxAsLines(lineset, &scaledCorners[0]);
311  hitsFound = true;
312  }
313  /*
314  // AMT: don't add lines if hit is not found becuse of unconsistency of scaling.
315  else
316  {
317  addBoxAsLines(lineset, corners);
318  }
319  */
320  }
321  if (!hitsFound)
322  fwLog(fwlog::kWarning) << Form("Can't find matching hits with for HCAL block %d in %s collection. Number of hits %d.\n", elIdx, cname.c_str(), (int)hitsandfracs.size());
323 
324 
325  }
326  else
327  {
328  fwLog(fwlog::kInfo) << "empty cluster \n";
329  }
330  }
331  } // endloop cand.elementsInBlocks();
332 
333 
334  if (boxset) {
335  boxset->RefitPlex();
336  setupAddElement(boxset, holder);
337  }
338 
339  if (lineset) {
340  setupAddElement(lineset, holder);
341  }
342 }
const fireworks::Context & context() const
float getScaleFactor3D() const
const FWGeometry * getGeom() const
Definition: Context.h:83
void setupAddElement(TEveElement *el, TEveElement *parent, bool set_color=true) const
FWViewEnergyScale * getEnergyScale() const
key_type key() const
Accessor for product key.
Definition: Ref.h:264
std::vector< ElementInBlock > ElementsInBlocks
Definition: PFCandidate.h:386
ProductID id() const
Accessor for product ID.
Definition: Ref.h:258
Particle flow rechit (rechit + geometry and topology information). See clustering algorithm in PFClus...
Definition: PFRecHit.h:35
T const * get() const
Returns C++ pointer to the item.
Definition: Ref.h:244
const float * getCorners(unsigned int id) const
Definition: FWGeometry.cc:280
#define fwLog(_level_)
Definition: fwLog.h:50
void viewContextBoxScale(const float *corners, float scale, bool plotEt, std::vector< float > &scaledCorners, const reco::PFRecHit *)
const reco::PFRecHitCollection * m_collectionHCAL
const reco::PFRecHit * getHitForDetId(unsigned detId)
double energy() const
rechit energy
Definition: PFRecHit.h:112
unsigned int lastIdx
bool getPlotEt() const
void FWPFCandidateWithHitsProxyBuilder::build ( const FWEventItem iItem,
TEveElementList *  product,
const FWViewContext vc 
)
virtual

Reimplemented from FWProxyBuilderBase.

Definition at line 79 of file FWPFCandidateWithHitsProxyBuilder.cc.

References addHitsForCandidate(), objects.IsoTrackAnalyzer::candidates, reco::LeafCandidate::charge(), FWDisplayProperties::color(), AlCaHLTBitMon_QueryRunRegistry::comp, FWProxyBuilderBase::context(), FWProxyBuilderBase::createCompound(), FWEventItem::defaultDisplayProperties(), FWEventItem::get(), customizeTrackingMonitorSeedNumber::idx, initPFRecHitsCollections(), reco::LeafCandidate::px(), reco::LeafCandidate::py(), reco::LeafCandidate::pz(), fireworks::setTrackTypePF(), FWProxyBuilderBase::setupAddElement(), lumiQTWidget::t, and reco::PFCandidate::vertex().

80 {
81  // init PFCandiate collection
83  iItem->get( candidates );
84  if( candidates == 0 ) return;
85 
86 
87  Int_t idx = 0;
89  for( reco::PFCandidateCollection::const_iterator it = candidates->begin(), itEnd = candidates->end(); it != itEnd; ++it, ++idx)
90  {
91  TEveCompound* comp = createCompound();
92  setupAddElement( comp, product );
93  // printf("products size %d/%d \n", (int)iItem->size(), product->NumChildren());
94 
95  const reco::PFCandidate& cand = *it;
96 
97  // track
98  {
99  TEveRecTrack t;
100  t.fBeta = 1.;
101  t.fP = TEveVector( cand.px(), cand.py(), cand.pz() );
102  t.fV = TEveVector( cand.vertex().x(), cand.vertex().y(), cand.vertex().z() );
103  t.fSign = cand.charge();
104  TEveTrack* trk = new TEveTrack(&t, context().getTrackPropagator() );
105  trk->MakeTrack();
106  fireworks::setTrackTypePF( cand, trk );
107  setupAddElement( trk, comp);
108  }
109  // hits
110  {
111  comp->SetMainColor(iItem->defaultDisplayProperties().color());
112  addHitsForCandidate(cand, comp, vc);
113  }
114 
115  }
116 }
const fireworks::Context & context() const
const FWDisplayProperties & defaultDisplayProperties() const
Definition: FWEventItem.cc:453
void setupAddElement(TEveElement *el, TEveElement *parent, bool set_color=true) const
void get(const T *&oData) const
Definition: FWEventItem.h:85
Color_t color() const
virtual int charge() const final
electric charge
Definition: LeafCandidate.h:91
virtual const Point & vertex() const
vertex position (overwritten by PF...)
Definition: PFCandidate.cc:647
virtual double py() const final
y coordinate of momentum vector
virtual double pz() const final
z coordinate of momentum vector
std::vector< reco::PFCandidate > PFCandidateCollection
collection of PFCandidates
tuple idx
DEBUGGING if hasattr(process,&quot;trackMonIterativeTracking2012&quot;): print &quot;trackMonIterativeTracking2012 D...
void setTrackTypePF(const reco::PFCandidate &pfCand, TAttLine *track)
Particle reconstructed by the particle flow algorithm.
Definition: PFCandidate.h:39
TEveCompound * createCompound(bool set_color=true, bool propagate_color_to_all_children=false) const
void addHitsForCandidate(const reco::PFCandidate &c, TEveElement *holder, const FWViewContext *vc)
virtual double px() const final
x coordinate of momentum vector
const std::string & FWPFCandidateWithHitsProxyBuilder::classPurpose ( )
static

Definition at line 344 of file FWPFCandidateWithHitsProxyBuilder.cc.

const std::string & FWPFCandidateWithHitsProxyBuilder::classRegisterTypeName ( )
static

Definition at line 344 of file FWPFCandidateWithHitsProxyBuilder.cc.

const std::string & FWPFCandidateWithHitsProxyBuilder::classTypeName ( )
static

Definition at line 344 of file FWPFCandidateWithHitsProxyBuilder.cc.

const std::string & FWPFCandidateWithHitsProxyBuilder::classView ( )
static

Definition at line 344 of file FWPFCandidateWithHitsProxyBuilder.cc.

const reco::PFRecHit * FWPFCandidateWithHitsProxyBuilder::getHitForDetId ( unsigned  detId)
private

Definition at line 174 of file FWPFCandidateWithHitsProxyBuilder.cc.

References m_collectionHCAL.

Referenced by addHitsForCandidate().

175 {
176 
177  for (reco::PFRecHitCollection::const_iterator it = m_collectionHCAL->begin(); it != m_collectionHCAL->end(); ++it)
178  {
179 
180  if ( it->detId() == candIdx)
181  {
182  return &(*it);
183  }
184  }
185  return 0;
186 }
const reco::PFRecHitCollection * m_collectionHCAL
virtual bool FWPFCandidateWithHitsProxyBuilder::havePerViewProduct ( FWViewType::EType  ) const
inlinevirtual

Reimplemented from FWProxyBuilderBase.

Definition at line 20 of file FWPFCandidateWithHitsProxyBuilder.h.

20 { return true; }
void FWPFCandidateWithHitsProxyBuilder::initPFRecHitsCollections ( )
private

Definition at line 119 of file FWPFCandidateWithHitsProxyBuilder.cc.

References fwLog, edm::EventBase::getByLabel(), FWEventItem::getEvent(), edm::HandleBase::isValid(), FWProxyBuilderBase::item(), fwlog::kError, m_collectionHCAL, FWEventItem::name(), and GlobalPosition_Frontier_DevDB_cff::tag.

Referenced by build().

120 {
121  // ref hcal collections
123 
124 
125  m_collectionHCAL =0;
126  try
127  {
128  // edm::InputTag tag("particleFlowRecHitHCAL");
129  edm::InputTag tag(cname);
130  item()->getEvent()->getByLabel(tag, handle_hits);
131  if (handle_hits.isValid())
132  {
133  m_collectionHCAL = &*handle_hits;
134  }
135  else
136  {
137  fwLog(fwlog::kError) <<"FWPFCandidateWithHitsProxyBuilder, item " << item()->name() <<": Failed to access collection with name " << cname << "." << std::endl;
138  }
139  }
140  catch (...)
141  {
142  fwLog(fwlog::kError) <<"FWPFCandidateWithHitsProxyBuilder, item " << item()->name() <<": Failed to access collection with name " << cname << "." << std::endl;
143  }
144 }
const std::string & name() const
Definition: FWEventItem.cc:502
const FWEventItem * item() const
bool isValid() const
Definition: HandleBase.h:75
const edm::EventBase * getEvent() const
Definition: FWEventItem.h:148
#define fwLog(_level_)
Definition: fwLog.h:50
const reco::PFRecHitCollection * m_collectionHCAL
bool getByLabel(InputTag const &, Handle< T > &) const
Definition: EventBase.h:90
const FWPFCandidateWithHitsProxyBuilder& FWPFCandidateWithHitsProxyBuilder::operator= ( const FWPFCandidateWithHitsProxyBuilder )
private
const std::string& FWPFCandidateWithHitsProxyBuilder::purpose ( ) const
inline
void FWPFCandidateWithHitsProxyBuilder::scaleProduct ( TEveElementList *  parent,
FWViewType::EType  type,
const FWViewContext vc 
)
virtual

Reimplemented from FWProxyBuilderBase.

Definition at line 189 of file FWPFCandidateWithHitsProxyBuilder.cc.

References reco::PFRecHit::detId(), reco::PFRecHit::energy(), FWGeometry::getCorners(), FWViewContext::getEnergyScale(), FWEventItem::getGeom(), FWViewEnergyScale::getPlotEt(), FWViewEnergyScale::getScaleFactor3D(), i, FWProxyBuilderBase::item(), AlCaHLTBitMon_ParallelJobs::p, pileupReCalc_HLTpaths::scale, and viewContextBoxScale().

190 {
191  std::vector<float> scaledCorners(24);
192 
193  float scale = vc->getEnergyScale()->getScaleFactor3D()/50;
194  for (TEveElement::List_i i=parent->BeginChildren(); i!=parent->EndChildren(); ++i)
195  {
196  if ((*i)->NumChildren() > 1)
197  {
198  TEveElement::List_i xx = (*i)->BeginChildren(); ++xx;
199  TEveBoxSet* boxset = dynamic_cast<TEveBoxSet*>(*xx);
200  ++xx;
201  TEveStraightLineSet* lineset = dynamic_cast<TEveStraightLineSet*>(*xx);
202  TEveChunkManager::iterator li(lineset->GetLinePlex());
203  li.next();
204 
205 
206  TEveChunkManager* plex = boxset->GetPlex();
207  if (plex->N())
208  {
209  for (int atomIdx=0; atomIdx < plex->Size(); ++atomIdx)
210  {
211 
212  TEveBoxSet::BFreeBox_t* atom = (TEveBoxSet::BFreeBox_t*)boxset->GetPlex()->Atom(atomIdx);
213  reco::PFRecHit* hit = (reco::PFRecHit*)boxset->GetUserData(atomIdx);
214  const float* corners = item()->getGeom()->getCorners(hit->detId());
215  viewContextBoxScale(corners, hit->energy()*scale, vc->getEnergyScale()->getPlotEt(), scaledCorners, hit);
216  memcpy(atom->fVertices, &scaledCorners[0], sizeof(atom->fVertices));
217 
218  editBoxInLineSet(li, &scaledCorners[0]);
219  }
220 
221  for (TEveProjectable::ProjList_i p = lineset->BeginProjecteds(); p != lineset->EndProjecteds(); ++p)
222  {
223  TEveStraightLineSetProjected* projLineSet = (TEveStraightLineSetProjected*)(*p);
224  projLineSet->UpdateProjection();
225  }
226  }
227  }
228  }
229 }
float getScaleFactor3D() const
int i
Definition: DBlmapReader.cc:9
unsigned detId() const
rechit detId
Definition: PFRecHit.h:106
FWViewEnergyScale * getEnergyScale() const
const FWEventItem * item() const
Particle flow rechit (rechit + geometry and topology information). See clustering algorithm in PFClus...
Definition: PFRecHit.h:35
const float * getCorners(unsigned int id) const
Definition: FWGeometry.cc:280
void viewContextBoxScale(const float *corners, float scale, bool plotEt, std::vector< float > &scaledCorners, const reco::PFRecHit *)
double energy() const
rechit energy
Definition: PFRecHit.h:112
const FWGeometry * getGeom() const
Definition: FWEventItem.cc:683
bool getPlotEt() const
const std::string& FWPFCandidateWithHitsProxyBuilder::typeName ( ) const
inline
const std::string& FWPFCandidateWithHitsProxyBuilder::view ( ) const
inline
void FWPFCandidateWithHitsProxyBuilder::viewContextBoxScale ( const float *  corners,
float  scale,
bool  plotEt,
std::vector< float > &  scaledCorners,
const reco::PFRecHit  
)
private

Definition at line 147 of file FWPFCandidateWithHitsProxyBuilder.cc.

References f, and i.

Referenced by addHitsForCandidate(), and scaleProduct().

148 {
149  static TEveVector vtmp;
150  vtmp.Set(0.f, 0.f, 0.f);
151  for( unsigned int i = 0; i < 24; i += 3 )
152  {
153  vtmp[0] += corners[i];
154  vtmp[1] += corners[i + 1];
155  vtmp[2] += corners[i + 2];
156  }
157  vtmp *= 1.f/8.f;
158 
159  if (plotEt)
160  {
161  scale *= vtmp.Perp()/vtmp.Mag();
162  }
163 
164  // Coordinates for a scaled version of the original box
165  for( unsigned int i = 0; i < 24; i += 3 )
166  {
167  scaledCorners[i] = vtmp[0] + ( corners[i] - vtmp[0] ) * scale;
168  scaledCorners[i + 1] = vtmp[1] + ( corners[i + 1] - vtmp[1] ) * scale;
169  scaledCorners[i + 2] = vtmp[2] + ( corners[i + 2] - vtmp[2] ) * scale;
170  }
171 }
int i
Definition: DBlmapReader.cc:9
double f[11][100]

Member Data Documentation

const reco::PFRecHitCollection* FWPFCandidateWithHitsProxyBuilder::m_collectionHCAL
private