CMS 3D CMS Logo

FWGeoTopNode.cc
Go to the documentation of this file.
1 // -*- C++ -*-
2 //
3 // Package: Core
4 // Class : FWGeoTopNode
5 //
6 // Implementation:
7 // [Notes on implementation]
8 //
9 // Original Author: Matevz Tadel, Alja Mrak Tadel
10 // Created: Thu Jun 23 01:24:51 CEST 2011
11 //
12 
13 // system include files
14 
15 // user include files
16 
17 #include "TEveTrans.h"
18 #include "TEveViewer.h"
19 #include "TEveManager.h"
20 #include "TEveUtil.h"
21 
22 
23 #include "TROOT.h"
24 #include "TBuffer3D.h"
25 #include "TBuffer3DTypes.h"
26 #include "TVirtualViewer3D.h"
27 #include "TColor.h"
28 #include "TGLScenePad.h"
29 #include "TGLPhysicalShape.h"
30 #include "TGLSelectRecord.h"
31 #include "TGLViewer.h"
32 #include "TGLWidget.h"
33 
34 #include "TGeoShape.h"
35 #include "TGeoVolume.h"
36 #include "TGeoNode.h"
37 #include "TGeoShapeAssembly.h"
38 #include "TGeoCompositeShape.h"
39 #include "TGeoBoolNode.h"
40 #include "TGeoManager.h"
41 #include "TGeoMatrix.h"
42 #include "TVirtualGeoPainter.h"
43 
49 
50 TGLViewer* FWGeoTopNode::s_pickedViewer = 0;
52 
53 UInt_t FWGeoTopNode::phyID(int tableIdx)
54 {
55 return UInt_t(tableIdx + 2);
56 
57 }
58 
59 int FWGeoTopNode::tableIdx(TGLPhysicalShape* ps)
60 {
61  return ps->ID() - 2;
62 }
63 //______________________________________________________________________________
64 void FWGeoTopNode::EraseFromSet(std::set<TGLPhysicalShape*>& sset, TGLPhysicalShape* id)
65 {
66  sset.erase(id);
67  SetStateOf(id);
68 }
69 
70 //______________________________________________________________________________
71 void FWGeoTopNode::ClearSet(std::set<TGLPhysicalShape*>& sset)
72 {
73  while (!sset.empty())
74  {
75  TGLPhysicalShape *id = *sset.begin();
76  sset.erase(id);
77  SetStateOf(id);
78  }
79 }
80 //______________________________________________________________________________
81 void FWGeoTopNode::SetStateOf(TGLPhysicalShape* id)
82 {
84 
85  if (fSted.find(id) != fSted.end())
86  {
87  id->Select(1);
89  }
90  else if (fHted.find(id) != fHted.end())
91  {
92  id->Select(3);
94  }
95  else
96  {
97  id->Select(0);
100  }
101 
102 }
103 
104 //______________________________________________________________________________
105 void FWGeoTopNode::ProcessSelection(TGLSelectRecord& rec, std::set<TGLPhysicalShape*>& sset, TGLPhysicalShape* id)
106 {
107  // printf("FWGeoTopNode::ProcessSelection ===============================\n");
108 
109  m_scene->BeginUpdate();
110 
111  if (sset.empty())
112  {
113  if (id)
114  {
115  sset.insert(id);
116  rec.SetSecSelResult(TGLSelectRecord::kEnteringSelection);
117  }
118  }
119  else
120  {
121  if (id)
122  {
123  if (rec.GetMultiple())
124  {
125  if (sset.find(id) == sset.end())
126  {
127  sset.insert(id);
128  rec.SetSecSelResult(TGLSelectRecord::kModifyingInternalSelection);
129  }
130  else
131  {
132  EraseFromSet(sset, id);
133  if (sset.empty())
134  rec.SetSecSelResult(TGLSelectRecord::kLeavingSelection);
135  else
136  rec.SetSecSelResult(TGLSelectRecord::kModifyingInternalSelection);
137  }
138  }
139  else
140  {
141  if (sset.size() != 1 || sset.find(id) == sset.end())
142  {
143  ClearSet(sset);
144  sset.insert(id);
145  rec.SetSecSelResult(TGLSelectRecord::kModifyingInternalSelection);
146  }
147  }
148  }
149  else
150  {
151  if (!rec.GetMultiple())
152  {
153  ClearSet(sset);
154  rec.SetSecSelResult(TGLSelectRecord::kLeavingSelection);
155  }
156  }
157  }
158 
159  if (id)
160  {
161  SetStateOf(id);
162  }
163 
164  if (rec.GetSecSelResult() != TGLSelectRecord::kNone)
165  {
166  m_scene->EndUpdate(kTRUE, kFALSE, kTRUE);
167  gEve->Redraw3D();
168 
170  }
171  else
172  {
173  m_scene->EndUpdate(kFALSE, kFALSE, kFALSE);
174  }
175 }
176 
177 //______________________________________________________________________________
179 {
180  // printf("FWGeoTopNode::selectPhysicalFromTable
181 
182  TGLPhysicalShape* ps = m_scene->FindPhysical(phyID(tableIndex));
183  if (ps) {
184  fSted.insert(ps);
185  ps->Select(1);
186  // printf("selectPhysicalFromTable found physical \n");
187  return true;
188  }
189  else if ( tableManager()->refEntries().at(tableIndex).testBit(FWGeometryTableManagerBase::kVisNodeSelf))
190  fwLog(fwlog::kInfo) << "Selected entry not drawn in GL viewer. \n" ;
191  return false;
192 }
193 
194 //______________________________________________________________________________
196 {
197  for (std::set<TGLPhysicalShape*>::iterator it = fSted.begin(); it != fSted.end(); ++it)
198  {
199  printf("FWGeoTopNode::printSelected %s \n", tableManager()->refEntries().at(tableIdx(*it)).name() );
200  }
201 }
202 
203 //______________________________________________________________________________
204 
206 {
207  // Note: if object would be rendered, this would return fSted.begin().
208 
209  if (fSted.size() <= 1)
210  {
211  int cnt = 0;
212  for (FWGeometryTableManagerBase::Entries_i i = tableManager()->refEntries().begin(); i != tableManager()->refEntries().end(); ++i, ++cnt)
213  {
214  if (i->testBit(FWGeometryTableManagerBase::kSelected)) return cnt;
215  }
216  }
217  return -1;
218 }
219 
220 //______________________________________________________________________________
222 {
223  // Fill bounding-box information. Virtual from TAttBBox.
224 
225  BBoxZero(1.0f);
226 }
227 
228 //______________________________________________________________________________
229 void FWGeoTopNode::setupBuffMtx(TBuffer3D& buff, const TGeoHMatrix& mat)
230 {
231  const Double_t *r = mat.GetRotationMatrix();
232  const Double_t *t = mat.GetTranslation();
233  const Double_t *s = mat.GetScale();
234  Double_t *m = buff.fLocalMaster;
235  m[0] = r[0]*s[0]; m[1] = r[1]*s[1]; m[2] = r[2]*s[2]; m[3] = 0;
236  m[4] = r[3]*s[0]; m[5] = r[4]*s[1]; m[6] = r[5]*s[2]; m[7] = 0;
237  m[8] = r[6]*s[0]; m[9] = r[7]*s[1]; m[10] = r[8]*s[2]; m[11] = 0;
238  m[12] = t[0]; m[13] = t[1]; m[15] = t[2]; m[15] = 1;
239 
240  buff.fLocalFrame = kTRUE;
241 }
242 
243 // ______________________________________________________________________
244 void FWGeoTopNode::paintShape( Int_t tableIndex, const TGeoHMatrix& nm, bool volumeColor, bool isParentNode)
245 {
246  static const TEveException eh("FWGeoTopNode::paintShape ");
247 
248  // printf("paint sahpe id %d\n", tableIndex );
249 
251  UChar_t transparency = wrapTransparency(data, isParentNode);
252  // printf("trans %d \n", transparency );
253  if (transparency >= 100) return;
254 
255  if (data.m_node->GetVolume()->IsAssembly()) return;
256  TGeoShape* shape = data.m_node->GetVolume()->GetShape();
257 
258  TGeoCompositeShape* compositeShape = dynamic_cast<TGeoCompositeShape*>(shape);
259  if (compositeShape)
260  {
261  // m_scene->fNextCompositeID = phyID(tableIndex);
262 
263  Double_t halfLengths[3] = { compositeShape->GetDX(), compositeShape->GetDY(), compositeShape->GetDZ() };
264 
265  TBuffer3D buff(TBuffer3DTypes::kComposite);
266  buff.fID = data.m_node->GetVolume();
267  buff.fColor = volumeColor ? data.m_node->GetVolume()->GetLineColor() : data.m_color ;
268  buff.fTransparency = transparency;// data.m_node->GetVolume()->GetTransparency();
269 
270  nm.GetHomogenousMatrix(buff.fLocalMaster);
271  buff.fLocalFrame = kTRUE; // Always enforce local frame (no geo manager).
272  buff.SetAABoundingBox(compositeShape->GetOrigin(), halfLengths);
273  buff.SetSectionsValid(TBuffer3D::kCore|TBuffer3D::kBoundingBox);
274 
275  Bool_t paintComponents = kTRUE;
276  // Start a composite shape, identified by this buffer
277  if (TBuffer3D::GetCSLevel() == 0) {
278  paintComponents = m_scene->OpenCompositeWithPhyID(phyID(tableIndex), buff);
279  }
280 
281  TBuffer3D::IncCSLevel();
282 
283  // Paint the boolean node - will add more buffers to viewer
284  TGeoHMatrix xxx;
285  TGeoMatrix *gst = TGeoShape::GetTransform();
286  TGeoShape::SetTransform(&xxx);
287 
288  if (paintComponents) compositeShape->GetBoolNode()->Paint("");
289  TGeoShape::SetTransform(gst);
290  // Close the composite shape
291  if (TBuffer3D::DecCSLevel() == 0)
292  gPad->GetViewer3D()->CloseComposite();
293 
294 
295  // m_scene->fNextCompositeID = 0;
296  }
297  else
298  {
299  TBuffer3D& buff = (TBuffer3D&) shape->GetBuffer3D (TBuffer3D::kCore, kFALSE);
300  setupBuffMtx(buff, nm);
301  buff.fID = data.m_node->GetVolume();
302  buff.fColor = volumeColor ? data.m_node->GetVolume()->GetLineColor() : data.m_color ;
303  buff.fTransparency = transparency;// data.m_node->GetVolume()->GetTransparency();
304 
305 
306  nm.GetHomogenousMatrix(buff.fLocalMaster);
307  buff.fLocalFrame = kTRUE; // Always enforce local frame (no geo manager).
308 
309  Int_t sections = TBuffer3D::kBoundingBox | TBuffer3D::kShapeSpecific;
310  shape->GetBuffer3D(sections, kTRUE);
311 
312  Int_t reqSec = gPad->GetViewer3D()->AddObject(phyID(tableIndex), buff);
313 
314  if (reqSec != TBuffer3D::kNone) {
315  // This shouldn't happen, but I suspect it does sometimes.
316  if (reqSec & TBuffer3D::kCore)
317  Warning(eh, "Core section required again for shape='%s'. This shouldn't happen.", GetName());
318  shape->GetBuffer3D(reqSec, kTRUE);
319  reqSec = gPad->GetViewer3D()->AddObject(phyID(tableIndex), buff);
320  }
321 
322  if (reqSec != TBuffer3D::kNone)
323  Warning(eh, "Extra section required: reqSec=%d, shape=%s.", reqSec, GetName());
324  }
325 }
326 
327 // ______________________________________________________________________
328 void FWGeoTopNode::Paint(Option_t* opt)
329 {
330  static const TEveException eh("FWGeoTopNode::Paint ");
331 
332  TBuffer3D buff(TBuffer3DTypes::kGeneric);
333 
334  // Section kCore
335  buff.fID = this;
336  buff.fColor = GetMainColor();
337  buff.fTransparency = GetMainTransparency();
338  if (HasMainTrans()) RefMainTrans().SetBuffer3D(buff);
339 
340  buff.SetSectionsValid(TBuffer3D::kCore);
341 
342  Int_t reqSections = gPad->GetViewer3D()->AddObject(1, buff);
343  if (reqSections != TBuffer3D::kNone)
344  {
345  Warning(eh, "IsA='%s'. Viewer3D requires more sections (%d). Only direct-rendering supported.",
346  ClassName(), reqSections);
347  }
348 }
349 
350 // ______________________________________________________________________
352 {
353  ClearSet(fSted);
354  for (FWGeometryTableManagerBase::Entries_i i = tableManager()->refEntries().begin(); i != tableManager()->refEntries().end(); ++i)
356 }
357 
358 // ______________________________________________________________________
360 {
361  ClearSet(fHted);
362  for (FWGeometryTableManagerBase::Entries_i i = tableManager()->refEntries().begin(); i != tableManager()->refEntries().end(); ++i)
364 }
365 
366 
367 // ______________________________________________________________________
368 FWPopupMenu* FWGeoTopNode::setPopupMenu(int iX, int iY, TGLViewer* v, bool overlap)
369 {
370  if (getFirstSelectedTableIndex() < 0)
371  {
372  if (fSted.empty()) fwLog(fwlog::kInfo) << "No menu -- no node/entry selected \n";
373  return 0;
374  }
375 
376  FWPopupMenu* nodePopup = new FWPopupMenu();
377 
378  nodePopup->AddEntry("Set As Top Node", kSetTopNode);
379  nodePopup->AddEntry("Set As Top Node and Reset Camera", kSetTopNodeCam);
380  nodePopup->AddSeparator();
381  if (v) {
382  nodePopup->AddEntry("Rnr Off", kVisSelfOff);
383  }
384  nodePopup->AddEntry("Turn Render On For Children", kVisChldOn);
385  nodePopup->AddEntry("Turn Render Off For Children", kVisChldOff);
386  nodePopup->AddEntry("Apply Color To Children", kApplyChldCol);
387  nodePopup->AddEntry("Apply Color Recursively", kApplyChldColRec);
388  nodePopup->AddSeparator();
389 
390  if (overlap)
391  nodePopup->AddEntry("Print Overlap", kPrintOverlap);
392  nodePopup->AddEntry("Print Path", kPrintPath);
393  nodePopup->AddEntry("Print Shape", kPrintShape);
394  nodePopup->AddEntry("Print Material", kPrintMaterial);
395 
396  nodePopup->AddSeparator();
397  if (v) {
398  Window_t wdummy;
399  Int_t x,y;
400  gVirtualX->TranslateCoordinates(gClient->GetDefaultRoot()->GetId(), v->GetGLWidget()->GetId(), iX, iY, x, y, wdummy);
401  TGLVector3 pnt(x, y, 0.5*v->GetSelRec().GetMinZ());
402  v->CurrentCamera().WindowToViewport(pnt);
403  s_pickedCamera3DCenter = v->CurrentCamera().ViewportToWorld(pnt);
404  // s_pickedCamera3DCenter.Dump();
405  s_pickedViewer = v;
406 
407  nodePopup->AddEntry("Set Camera Center", kCamera);
408  }
409 
410  nodePopup->PlaceMenu(iX, iY,true,true);
411  return nodePopup;
412 }
413 
414 
416 {
417  if (isParentNode)
418  {
419  return TMath::Max((Char_t)browser()->getMinParentTransparency(), data.m_transparency) * browser()->getParentTransparencyFactor();
420  }
421  else
422  {
423  return TMath::Max((Char_t)browser()->getMinLeafTransparency(), data.m_transparency) * browser()->getLeafTransparencyFactor();
424  }
425 }
std::set< TGLPhysicalShape * > fSted
Definition: FWGeoTopNode.h:74
virtual FWGeometryTableManagerBase * tableManager()
Definition: FWGeoTopNode.h:70
void ClearSet(std::set< TGLPhysicalShape * > &sset)
Definition: FWGeoTopNode.cc:71
bool OpenCompositeWithPhyID(UInt_t phyID, const TBuffer3D &buffer)
virtual void Paint(Option_t *option="")
void printSelected()
void dataChanged()
Classes which inherit from FWTableManagerBase must call this when their underlying data changes...
int getFirstSelectedTableIndex()
bool overlap(const reco::Muon &muon1, const reco::Muon &muon2, double pullX=1.0, double pullY=1.0, bool checkAdjacentChambers=false)
static const G4AffineTransform & GetTransform(const G4TouchableHistory *touchable, int depth)
FWGeoTopNodeGLScene * m_scene
Definition: FWGeoTopNode.h:68
UChar_t wrapTransparency(FWGeometryTableManagerBase::NodeInfo &data, bool leafNode)
double f[11][100]
virtual void ComputeBBox()
T Max(T a, T b)
Definition: MathUtil.h:44
FWPopupMenu * setPopupMenu(int iX, int iY, TGLViewer *v, bool)
#define fwLog(_level_)
Definition: fwLog.h:50
virtual FWGeometryTableViewBase * browser()
Definition: FWGeoTopNode.h:71
bool selectPhysicalFromTable(int)
static TGLViewer * s_pickedViewer
Definition: FWGeoTopNode.h:87
std::set< TGLPhysicalShape * > fHted
Definition: FWGeoTopNode.h:73
void paintShape(Int_t idx, const TGeoHMatrix &nm, bool volumeColor, bool parentNode)
#define begin
Definition: vmac.h:30
char data[epos_bytes_allocation]
Definition: EPOS_Wrapper.h:82
virtual void UnHighlighted()
void EraseFromSet(std::set< TGLPhysicalShape * > &sset, TGLPhysicalShape *id)
Definition: FWGeoTopNode.cc:64
void ProcessSelection(TGLSelectRecord &rec, std::set< TGLPhysicalShape * > &sset, TGLPhysicalShape *id)
void SetStateOf(TGLPhysicalShape *id)
Definition: FWGeoTopNode.cc:81
static int tableIdx(TGLPhysicalShape *ps)
Definition: FWGeoTopNode.cc:59
virtual void UnSelected()
static TGLVector3 s_pickedCamera3DCenter
Definition: FWGeoTopNode.h:86
void setupBuffMtx(TBuffer3D &buff, const TGeoHMatrix &mat)
static UInt_t phyID(int tableIdx)
Definition: FWGeoTopNode.cc:53