[u/mrichter/AliRoot.git] / EVE / EveBase / AliEveHFEditor.cxx

// @(#)root/eve:$Id$
// Main author: Davide Caffarri 2009

/**************************************************************************
 * Copyright(c) 1998-2008, ALICE Experiment at CERN, all rights reserved. *
 * See http://aliceinfo.cern.ch/Offline/AliRoot/License.html for          *
 * full copyright notice.                                                 *
 **************************************************************************/

#include "AliEveHFEditor.h"
#include "AliEveHF.h"

#include "TVirtualPad.h"
#include "TColor.h"
#include <TDatabasePDG.h>

#include <TEveTrack.h>

// Cleanup these includes:
#include "TGLabel.h"
#include "TGButton.h"
#include "TGNumberEntry.h"
#include "TGColorSelect.h"
#include "TGDoubleSlider.h"


//______________________________________________________________________________
// GUI editor for AliEveHF.
//

ClassImp(AliEveHFEditor)

//______________________________________________________________________________
AliEveHFEditor::AliEveHFEditor(const TGWindow *p, Int_t width, Int_t height,
                               UInt_t options, Pixel_t back) :
  TGedFrame(p, width, height, options | kVerticalFrame, back),
  fM(0),
  fnProng(0),
  fInfoLabel0(0),
  fInfoLabel1(0),
  fInfoLabel2(0),
  fXButton(0)
  // Initialize widget pointers to 0
{
  // Constructor.

  MakeTitle("AliEveHF");

  fInfoLabel0 = new TGLabel(this);
  fInfoLabel0->SetTextJustify(kTextLeft);
  AddFrame(fInfoLabel0, new TGLayoutHints(kLHintsLeft|kLHintsExpandX,
                                          8, 0, 2, 0));

  fInfoLabel1 = new TGLabel(this);
  fInfoLabel1->SetTextJustify(kTextLeft);
  AddFrame(fInfoLabel1, new TGLayoutHints(kLHintsLeft|kLHintsExpandX,
                                          8, 0, 2, 0));

  fInfoLabel2 = new TGLabel(this);
  fInfoLabel2->SetTextJustify(kTextLeft);
  AddFrame(fInfoLabel2, new TGLayoutHints(kLHintsLeft|kLHintsExpandX,
                                          8, 0, 2, 0));

  fXButton = new TGTextButton(this, "Detailed View");
  AddFrame(fXButton, new TGLayoutHints(kLHintsLeft|kLHintsExpandX, 1, 1, 0, 0));
  fXButton->Connect("Clicked()", "AliEveHFEditor", this, "DisplayDetailed()");
}

/******************************************************************************/

//______________________________________________________________________________
void AliEveHFEditor::SetModel(TObject* obj)
{
  // Set model object.

  fM = static_cast<AliEveHF*>(obj);

  // Set values of widgets
  fInfoLabel0->SetText(Form("CosPointingAngle = %f",  fM->GetCosPointingAngle()));
  fInfoLabel1->SetText(Form("Pt = %f", fM->GetPt()));

  fInfoLabel2->SetText(Form("PDG Invariant Mass = %f", TDatabasePDG::Instance()->GetParticle(421)->Mass()));

}

/******************************************************************************/

// Implements callback/slot methods

//______________________________________________________________________________
// void AliEveHFEditor::DoXYZZ()
// {
//    // Slot for XYZZ.
//
//    fM->SetXYZZ(fXYZZ->GetValue());
//    Update();
// }

#include <TEveManager.h>
#include <TEveWindow.h>
#include <TEveViewer.h>
#include <TEveScene.h>
#include <TEveGeoNode.h>
#include <TEveProjectionManager.h>

#include <TGLCamera.h>
#include <TGLViewer.h>
#include "TGLCameraOverlay.h"

#include <TLatex.h>
#include <TRootEmbeddedCanvas.h>
#include <TInterpreter.h>

void AliEveHFEditor::DisplayDetailed()
{
  TEveWindowSlot *slot = TEveWindow::CreateWindowMainFrame();
  TEveWindowPack *pack = slot->MakePack();
  pack->SetShowTitleBar(kFALSE);
  pack->SetHorizontal();

  //
  // This part is for getting the different objects to display
  //
  char displayInfo[100] = {0};
  snprintf(displayInfo,100,"pt = %.3f",fM->GetPt());
  TEveLine *lhfTransverseMomentumDirection = new TEveLine(displayInfo);
  lhfTransverseMomentumDirection->SetLineColor(kOrange+8);
  lhfTransverseMomentumDirection->SetLineWidth(2);
  lhfTransverseMomentumDirection->SetLineStyle(2);
  lhfTransverseMomentumDirection->SetLineWidth(2);
  Float_t scalePt = 100.; // this needs to be available as a ruler
  lhfTransverseMomentumDirection->SetPoint(0,fM->fRecDecayHF.fX, fM->fRecDecayHF.fY, fM->fRecDecayHF.fZ);
  lhfTransverseMomentumDirection->SetPoint(1,scalePt*fM->fRecDecayMomHF.fX, scalePt*fM->fRecDecayMomHF.fY,0);

  TEvePointSet *pvlocation = new TEvePointSet("PV location");
  pvlocation->SetNextPoint(fM->fRecBirthHF.fX, fM->fRecBirthHF.fY, fM->fRecBirthHF.fZ);
  pvlocation->SetTitle("pv location");
  pvlocation->SetMarkerStyle(4);
  pvlocation->SetMarkerSize(2.5);
  pvlocation->SetMarkerColor(7);

  TEvePointSet *hflocation = new TEvePointSet("HF location");
  hflocation->SetNextPoint(fM->fRecDecayHF.fX, fM->fRecDecayHF.fY, fM->fRecDecayHF.fZ);
  hflocation->SetTitle("HF location");
  hflocation->SetMarkerStyle(4);
  hflocation->SetMarkerSize(2.5);
  hflocation->SetMarkerColor(kOrange+8);

  TEveUtil::LoadMacro("clusters_from_index.C");

  TEveTrack *negTrack =  fM->GetNegTrack();
  TEveTrack *posTrack =  fM->GetPosTrack();


  char macroWithIndex[100] = {0};
  Int_t daughterIndex = 0;
  TEvePointSet *negDaughterCluster = 0;
  TEvePointSet *posDaughterCluster = 0;

  daughterIndex = negTrack->GetIndex();
  snprintf(macroWithIndex,100,"clusters_from_index(%d)",daughterIndex);
  Long_t negResult = gInterpreter->ProcessLine(macroWithIndex);
  if (negResult) {
    negDaughterCluster = reinterpret_cast<TEvePointSet*>(negResult);
    if (negDaughterCluster){
      negDaughterCluster->SetMarkerStyle(4);
      negDaughterCluster->SetMarkerSize(1.5);
      negDaughterCluster->SetMarkerColor(kBlue+3);
    }
  }
  else
  {
    Warning("DisplayDetailed", "Import of negative daughter's clusters failed.");
  }

  daughterIndex = posTrack->GetIndex();
  snprintf(macroWithIndex,100,"clusters_from_index(%d)",daughterIndex);
  Long_t posResult = gInterpreter->ProcessLine(macroWithIndex);
  if (posResult) {
    posDaughterCluster = reinterpret_cast<TEvePointSet*>(posResult);
    if (posDaughterCluster){
      posDaughterCluster->SetMarkerStyle(4);
      posDaughterCluster->SetMarkerSize(1.5);
      posDaughterCluster->SetMarkerColor(kRed+3);
    }
  }
  else
  {
    Warning("DisplayDetailed", "Import of positive daughter's clusters failed.");
  }

  //
  // This part is for the bending plane view
  //
  pack->NewSlot()->MakeCurrent();
  TEveViewer *bpViewer = gEve->SpawnNewViewer("HF bending plane View");
  TEveScene  *bpScene  = gEve->SpawnNewScene("HF bending plane Scene");

  TEveUtil::LoadMacro("geom_gentle.C");
  Long_t result = gInterpreter->ProcessLine("geom_gentle_rphi()");
  if (result)
  {
    TEveGeoShape *geomRPhi = reinterpret_cast<TEveGeoShape*>(result);
    geomRPhi->IncDenyDestroy();
    TEveProjectionManager *projMgr = new TEveProjectionManager(TEveProjection::kPT_RPhi);
    projMgr->ImportElements(geomRPhi, bpScene);
  }
  else
  {
    Warning("DisplayDetailed", "Import of R-Phi geometry failed.");
  }

  bpViewer->AddScene(bpScene);
  bpScene->AddElement(fM);
  bpScene->AddElement(lhfTransverseMomentumDirection);
  bpScene->AddElement(pvlocation);
  bpScene->AddElement(hflocation);

  if (negDaughterCluster) bpScene->AddElement(negDaughterCluster);
  if (posDaughterCluster) bpScene->AddElement(posDaughterCluster);

  // This is the to-do list for the bending plane:
  // 1. fix the view to orthographic XOY (no rotation allowed but moving the center ok) ->done!
  // 2. show axis and tickles along X and Y ->done!
  //       -> note for the projection the cartesian scales are not very useful
  //       -> propose a phi and R scale which rotates with a reset at 0;
  //       -> propose a transformation for an eta scale (keep the z one);
  // 3. show the center, the main vertex and the detectors for this view ->done!
  // 4. show V0 direction in the bending plane with arrow length proportional to pT ->done!
  // 5. show angles with respect to axis (phi angle) ->almost.
  // 6. show clusters in the ITS and in the TPC associated with the daughter tracks
  //       -> include a radius cut for plotting only ITS and TPC ->done!
  bpViewer->GetGLViewer()->SetCurrentCamera(TGLViewer::kCameraOrthoXOY);
  bpViewer->GetGLViewer()->ResetCamerasAfterNextUpdate();
  TGLViewer *lbpGLViewer = bpViewer->GetGLViewer();
  TGLCameraOverlay* co = lbpGLViewer->GetCameraOverlay();
  co->SetShowOrthographic(true); //(false);
  co->SetOrthographicMode(TGLCameraOverlay::kAxis); // ::kPlaneIntersect or ::kBar
  // end of the bending plane part

  //
  // This part is for the decay plane view
  //
  pack->NewSlot()->MakeCurrent();
  TEveViewer *dpViewer = gEve->SpawnNewViewer("HF decay plane View");
  TEveScene  *dpScene  = gEve->SpawnNewScene("HF decay plane Scene");

  dpViewer->AddScene(dpScene);

  result = gInterpreter->ProcessLine("geom_gentle(kFALSE)");
  if (result)
  {
    TEveGeoShape *geom = reinterpret_cast<TEveGeoShape*>(result);
    geom->IncDenyDestroy();
    geom->FindChild("TRD+TOF")->SetRnrState(kFALSE);
    geom->FindChild("PHOS")   ->SetRnrState(kFALSE);
    geom->FindChild("HMPID")  ->SetRnrState(kFALSE);
    dpScene->AddElement(geom);
  }
  else
  {
    Warning("DisplayDetailed", "Import of 3D geometry failed.");
  }

  dpScene->AddElement(fM);
  dpScene->AddElement(lhfTransverseMomentumDirection);
  dpScene->AddElement(pvlocation);
  dpScene->AddElement(hflocation);
  if (negDaughterCluster) dpScene->AddElement(negDaughterCluster);
  if (posDaughterCluster) dpScene->AddElement(posDaughterCluster);

  // This is the to-do list for the decay plane:
  // 1. fix the view to decay plane (no rotation allowed but moving the center ok)
  // 2. show V0 direction with a vertical arrow length proportional to pT -> done!
  // 3. show the center, the main vertex and the detectors for this view -> done!
  // 4. show the x,y and z axis and the different angles
  //       -> this needs a referential object that we can move around
  //          or fix to a selected point (origin being the default)
  // 5. draw the dca between daughters and the extrapolation to the main vertex.
  //       -> this is an issue since we only store the distance: check with J.Belikov
  // 6. show clusters in the ITS and in the TPC associated with the daughter tracks
  //       -> include a radius cut for plotting only ITS and TPC ->done!
  dpViewer->GetGLViewer()->ResetCamerasAfterNextUpdate();
  TGLCamera& dpCam = dpViewer->GetGLViewer()->CurrentCamera();
  dpCam.SetExternalCenter(kTRUE);
  dpCam.SetCenterVec(fM->fRecDecayHF.fX, fM->fRecDecayHF.fY, fM->fRecDecayHF.fZ);
  dpCam.RotateRad(0,-TMath::Pi()/2.); // RotateRad rotates in radians (hRotate,vRotate)
  //  Here rotate the _view_ (not the camera) by (fM->GetPhi() - TMath::Pi()/2.)

  // In the end maybe truck and rotate properly...
  //  dpCam.Truck(0,200);// Truck center wrt the view panel: (x=0 pixels, y pixels)
  //  dpCam.Rotate(0,50,0,0); // Rotate in pixels (xDelta,yDelta)

  // end of the decay plane part

  //
  // This part is for displaying the information
  //
  slot = pack->NewSlot();

  TEveWindowFrame *frame = slot->MakeFrame(new TRootEmbeddedCanvas());
  frame->SetElementName("Details");

  // Print and show detailed information about the HF
  // Calculation of the invariant mass with the max prob PID hypothesis first
  // pseudorapidity, phi angle, pt, radius, dcas
  char info[100] = {0};

  snprintf(info,100,"Cos Pointing angle = %.3f ",fM->GetCosPointingAngle());
  TLatex* ltx = new TLatex(0.05, 0.9, info);
  ltx->SetTextSize(0.08);
  ltx->DrawLatex(0.05, 0.8, info);

  snprintf(info,100,"p_{T} = %.3f [GeV/c]",fM->GetPt());
  ltx->DrawLatex(0.05, 0.7, info);

  snprintf(info,100,"#eta = - ln( tan(#theta/2) ) = %.3f",fM->GetEta());
  ltx->DrawLatex(0.05, 0.6, info);

  snprintf(info,100,"D^{0} inv. mass = %.3f", fM->GetInvariantMassPart());
  ltx->DrawLatex(0.05, 0.5, info);

  snprintf(info,100,"D^{0}bar inv. mass = %.3f", fM->GetInvariantMassAntiPart());
  ltx->DrawLatex(0.05, 0.4, info);

  gEve->Redraw3D();
}
Commit	Line	Data
fa3634bd	1	// @(#)root/eve:$Id$
	2	// Main author: Davide Caffarri 2009
	3
	4	/**************************************************************************
	5	* Copyright(c) 1998-2008, ALICE Experiment at CERN, all rights reserved. *
	6	* See http://aliceinfo.cern.ch/Offline/AliRoot/License.html for *
	7	* full copyright notice. *
	8	**************************************************************************/
	9
	10	#include "AliEveHFEditor.h"
	11	#include "AliEveHF.h"
	12
	13	#include "TVirtualPad.h"
	14	#include "TColor.h"
c180f65d	15	#include <TDatabasePDG.h>
fa3634bd	16
	17	#include <TEveTrack.h>
	18
	19	// Cleanup these includes:
	20	#include "TGLabel.h"
	21	#include "TGButton.h"
	22	#include "TGNumberEntry.h"
	23	#include "TGColorSelect.h"
	24	#include "TGDoubleSlider.h"
	25
	26
	27	//______________________________________________________________________________
	28	// GUI editor for AliEveHF.
	29	//
	30
	31	ClassImp(AliEveHFEditor)
	32
	33	//______________________________________________________________________________
	34	AliEveHFEditor::AliEveHFEditor(const TGWindow *p, Int_t width, Int_t height,
	35	UInt_t options, Pixel_t back) :
	36	TGedFrame(p, width, height, options \| kVerticalFrame, back),
	37	fM(0),
	38	fnProng(0),
	39	fInfoLabel0(0),
	40	fInfoLabel1(0),
	41	fInfoLabel2(0),
	42	fXButton(0)
	43	// Initialize widget pointers to 0
	44	{
	45	// Constructor.
	46
	47	MakeTitle("AliEveHF");
	48
	49	fInfoLabel0 = new TGLabel(this);
	50	fInfoLabel0->SetTextJustify(kTextLeft);
	51	AddFrame(fInfoLabel0, new TGLayoutHints(kLHintsLeft\|kLHintsExpandX,
	52	8, 0, 2, 0));
	53
	54	fInfoLabel1 = new TGLabel(this);
	55	fInfoLabel1->SetTextJustify(kTextLeft);
	56	AddFrame(fInfoLabel1, new TGLayoutHints(kLHintsLeft\|kLHintsExpandX,
	57	8, 0, 2, 0));
	58
	59	fInfoLabel2 = new TGLabel(this);
	60	fInfoLabel2->SetTextJustify(kTextLeft);
	61	AddFrame(fInfoLabel2, new TGLayoutHints(kLHintsLeft\|kLHintsExpandX,
	62	8, 0, 2, 0));
	63
	64	fXButton = new TGTextButton(this, "Detailed View");
	65	AddFrame(fXButton, new TGLayoutHints(kLHintsLeft\|kLHintsExpandX, 1, 1, 0, 0));
	66	fXButton->Connect("Clicked()", "AliEveHFEditor", this, "DisplayDetailed()");
	67	}
	68
	69	/******************************************************************************/
	70
	71	//______________________________________________________________________________
	72	void AliEveHFEditor::SetModel(TObject* obj)
	73	{
	74	// Set model object.
	75
04a6d27f	76	fM = static_cast<AliEveHF*>(obj);
fa3634bd	77
	78	// Set values of widgets
	79	fInfoLabel0->SetText(Form("CosPointingAngle = %f", fM->GetCosPointingAngle()));
	80	fInfoLabel1->SetText(Form("Pt = %f", fM->GetPt()));
f8d5a389	81
fa3634bd	82	fInfoLabel2->SetText(Form("PDG Invariant Mass = %f", TDatabasePDG::Instance()->GetParticle(421)->Mass()));
f8d5a389	83
fa3634bd	84	}
	85
	86	/******************************************************************************/
	87
	88	// Implements callback/slot methods
	89
	90	//______________________________________________________________________________
	91	// void AliEveHFEditor::DoXYZZ()
	92	// {
	93	// // Slot for XYZZ.
	94	//
	95	// fM->SetXYZZ(fXYZZ->GetValue());
	96	// Update();
	97	// }
	98
	99	#include <TEveManager.h>
	100	#include <TEveWindow.h>
	101	#include <TEveViewer.h>
	102	#include <TEveScene.h>
	103	#include <TEveGeoNode.h>
	104	#include <TEveProjectionManager.h>
	105
	106	#include <TGLCamera.h>
	107	#include <TGLViewer.h>
	108	#include "TGLCameraOverlay.h"
	109
	110	#include <TLatex.h>
	111	#include <TRootEmbeddedCanvas.h>
	112	#include <TInterpreter.h>
	113
	114	void AliEveHFEditor::DisplayDetailed()
	115	{
	116	TEveWindowSlot *slot = TEveWindow::CreateWindowMainFrame();
	117	TEveWindowPack *pack = slot->MakePack();
	118	pack->SetShowTitleBar(kFALSE);
	119	pack->SetHorizontal();
	120
	121	//
	122	// This part is for getting the different objects to display
	123	//
	124	char displayInfo[100] = {0};
452b3436	125	snprintf(displayInfo,100,"pt = %.3f",fM->GetPt());
fa3634bd	126	TEveLine *lhfTransverseMomentumDirection = new TEveLine(displayInfo);
	127	lhfTransverseMomentumDirection->SetLineColor(kOrange+8);
	128	lhfTransverseMomentumDirection->SetLineWidth(2);
	129	lhfTransverseMomentumDirection->SetLineStyle(2);
	130	lhfTransverseMomentumDirection->SetLineWidth(2);
	131	Float_t scalePt = 100.; // this needs to be available as a ruler
	132	lhfTransverseMomentumDirection->SetPoint(0,fM->fRecDecayHF.fX, fM->fRecDecayHF.fY, fM->fRecDecayHF.fZ);
12cc265f	133	lhfTransverseMomentumDirection->SetPoint(1,scalePtfM->fRecDecayMomHF.fX, scalePtfM->fRecDecayMomHF.fY,0);
fa3634bd	134
	135	TEvePointSet *pvlocation = new TEvePointSet("PV location");
	136	pvlocation->SetNextPoint(fM->fRecBirthHF.fX, fM->fRecBirthHF.fY, fM->fRecBirthHF.fZ);
	137	pvlocation->SetTitle("pv location");
	138	pvlocation->SetMarkerStyle(4);
	139	pvlocation->SetMarkerSize(2.5);
	140	pvlocation->SetMarkerColor(7);
	141
	142	TEvePointSet *hflocation = new TEvePointSet("HF location");
	143	hflocation->SetNextPoint(fM->fRecDecayHF.fX, fM->fRecDecayHF.fY, fM->fRecDecayHF.fZ);
	144	hflocation->SetTitle("HF location");
	145	hflocation->SetMarkerStyle(4);
	146	hflocation->SetMarkerSize(2.5);
	147	hflocation->SetMarkerColor(kOrange+8);
	148
	149	TEveUtil::LoadMacro("clusters_from_index.C");
	150
	151	TEveTrack *negTrack = fM->GetNegTrack();
	152	TEveTrack *posTrack = fM->GetPosTrack();
	153
	154
	155	char macroWithIndex[100] = {0};
	156	Int_t daughterIndex = 0;
	157	TEvePointSet *negDaughterCluster = 0;
	158	TEvePointSet *posDaughterCluster = 0;
f8d5a389	159
fa3634bd	160	daughterIndex = negTrack->GetIndex();
452b3436	161	snprintf(macroWithIndex,100,"clusters_from_index(%d)",daughterIndex);
fa3634bd	162	Long_t negResult = gInterpreter->ProcessLine(macroWithIndex);
	163	if (negResult) {
	164	negDaughterCluster = reinterpret_cast<TEvePointSet*>(negResult);
	165	if (negDaughterCluster){
f8d5a389	166	negDaughterCluster->SetMarkerStyle(4);
f8d5a389	167	negDaughterCluster->SetMarkerSize(1.5);
fa3634bd	168	negDaughterCluster->SetMarkerColor(kBlue+3);
	169	}
	170	}
	171	else
	172	{
	173	Warning("DisplayDetailed", "Import of negative daughter's clusters failed.");
	174	}
	175
	176	daughterIndex = posTrack->GetIndex();
452b3436	177	snprintf(macroWithIndex,100,"clusters_from_index(%d)",daughterIndex);
fa3634bd	178	Long_t posResult = gInterpreter->ProcessLine(macroWithIndex);
	179	if (posResult) {
	180	posDaughterCluster = reinterpret_cast<TEvePointSet*>(posResult);
	181	if (posDaughterCluster){
f8d5a389	182	posDaughterCluster->SetMarkerStyle(4);
f8d5a389	183	posDaughterCluster->SetMarkerSize(1.5);
fa3634bd	184	posDaughterCluster->SetMarkerColor(kRed+3);
	185	}
	186	}
	187	else
	188	{
	189	Warning("DisplayDetailed", "Import of positive daughter's clusters failed.");
	190	}
	191
	192	//
	193	// This part is for the bending plane view
	194	//
	195	pack->NewSlot()->MakeCurrent();
	196	TEveViewer *bpViewer = gEve->SpawnNewViewer("HF bending plane View");
	197	TEveScene *bpScene = gEve->SpawnNewScene("HF bending plane Scene");
	198
	199	TEveUtil::LoadMacro("geom_gentle.C");
	200	Long_t result = gInterpreter->ProcessLine("geom_gentle_rphi()");
	201	if (result)
	202	{
	203	TEveGeoShape geomRPhi = reinterpret_cast<TEveGeoShape>(result);
	204	geomRPhi->IncDenyDestroy();
edd2b192	205	TEveProjectionManager *projMgr = new TEveProjectionManager(TEveProjection::kPT_RPhi);
fa3634bd	206	projMgr->ImportElements(geomRPhi, bpScene);
	207	}
	208	else
	209	{
	210	Warning("DisplayDetailed", "Import of R-Phi geometry failed.");
	211	}
	212
	213	bpViewer->AddScene(bpScene);
	214	bpScene->AddElement(fM);
	215	bpScene->AddElement(lhfTransverseMomentumDirection);
	216	bpScene->AddElement(pvlocation);
	217	bpScene->AddElement(hflocation);
	218
	219	if (negDaughterCluster) bpScene->AddElement(negDaughterCluster);
	220	if (posDaughterCluster) bpScene->AddElement(posDaughterCluster);
	221
	222	// This is the to-do list for the bending plane:
f8d5a389	223	// 1. fix the view to orthographic XOY (no rotation allowed but moving the center ok) ->done!
fa3634bd	224	// 2. show axis and tickles along X and Y ->done!
	225	// -> note for the projection the cartesian scales are not very useful
	226	// -> propose a phi and R scale which rotates with a reset at 0;
	227	// -> propose a transformation for an eta scale (keep the z one);
	228	// 3. show the center, the main vertex and the detectors for this view ->done!
	229	// 4. show V0 direction in the bending plane with arrow length proportional to pT ->done!
	230	// 5. show angles with respect to axis (phi angle) ->almost.
	231	// 6. show clusters in the ITS and in the TPC associated with the daughter tracks
	232	// -> include a radius cut for plotting only ITS and TPC ->done!
	233	bpViewer->GetGLViewer()->SetCurrentCamera(TGLViewer::kCameraOrthoXOY);
	234	bpViewer->GetGLViewer()->ResetCamerasAfterNextUpdate();
	235	TGLViewer *lbpGLViewer = bpViewer->GetGLViewer();
	236	TGLCameraOverlay* co = lbpGLViewer->GetCameraOverlay();
	237	co->SetShowOrthographic(true); //(false);
	238	co->SetOrthographicMode(TGLCameraOverlay::kAxis); // ::kPlaneIntersect or ::kBar
	239	// end of the bending plane part
	240
	241	//
	242	// This part is for the decay plane view
	243	//
	244	pack->NewSlot()->MakeCurrent();
	245	TEveViewer *dpViewer = gEve->SpawnNewViewer("HF decay plane View");
	246	TEveScene *dpScene = gEve->SpawnNewScene("HF decay plane Scene");
	247
	248	dpViewer->AddScene(dpScene);
	249
	250	result = gInterpreter->ProcessLine("geom_gentle(kFALSE)");
	251	if (result)
	252	{
	253	TEveGeoShape geom = reinterpret_cast<TEveGeoShape>(result);
	254	geom->IncDenyDestroy();
	255	geom->FindChild("TRD+TOF")->SetRnrState(kFALSE);
	256	geom->FindChild("PHOS") ->SetRnrState(kFALSE);
	257	geom->FindChild("HMPID") ->SetRnrState(kFALSE);
	258	dpScene->AddElement(geom);
	259	}
	260	else
	261	{
	262	Warning("DisplayDetailed", "Import of 3D geometry failed.");
	263	}
	264
	265	dpScene->AddElement(fM);
	266	dpScene->AddElement(lhfTransverseMomentumDirection);
	267	dpScene->AddElement(pvlocation);
	268	dpScene->AddElement(hflocation);
	269	if (negDaughterCluster) dpScene->AddElement(negDaughterCluster);
	270	if (posDaughterCluster) dpScene->AddElement(posDaughterCluster);
	271
	272	// This is the to-do list for the decay plane:
	273	// 1. fix the view to decay plane (no rotation allowed but moving the center ok)
	274	// 2. show V0 direction with a vertical arrow length proportional to pT -> done!
	275	// 3. show the center, the main vertex and the detectors for this view -> done!
	276	// 4. show the x,y and z axis and the different angles
	277	// -> this needs a referential object that we can move around
	278	// or fix to a selected point (origin being the default)
	279	// 5. draw the dca between daughters and the extrapolation to the main vertex.
	280	// -> this is an issue since we only store the distance: check with J.Belikov
	281	// 6. show clusters in the ITS and in the TPC associated with the daughter tracks
	282	// -> include a radius cut for plotting only ITS and TPC ->done!
	283	dpViewer->GetGLViewer()->ResetCamerasAfterNextUpdate();
	284	TGLCamera& dpCam = dpViewer->GetGLViewer()->CurrentCamera();
	285	dpCam.SetExternalCenter(kTRUE);
	286	dpCam.SetCenterVec(fM->fRecDecayHF.fX, fM->fRecDecayHF.fY, fM->fRecDecayHF.fZ);
	287	dpCam.RotateRad(0,-TMath::Pi()/2.); // RotateRad rotates in radians (hRotate,vRotate)
288	// Here rotate the _view_ (not the camera) by (fM->GetPhi() - TMath::Pi()/2.)
289
290	// In the end maybe truck and rotate properly...
291	// dpCam.Truck(0,200);// Truck center wrt the view panel: (x=0 pixels, y pixels)
292	// dpCam.Rotate(0,50,0,0); // Rotate in pixels (xDelta,yDelta)
293
294	// end of the decay plane part
295
296	//
297	// This part is for displaying the information
298	//
299	slot = pack->NewSlot();
f8d5a389	300
fa3634bd	301	TEveWindowFrame *frame = slot->MakeFrame(new TRootEmbeddedCanvas());
	302	frame->SetElementName("Details");
	303
	304	// Print and show detailed information about the HF
	305	// Calculation of the invariant mass with the max prob PID hypothesis first
	306	// pseudorapidity, phi angle, pt, radius, dcas
	307	char info[100] = {0};
	308
452b3436	309	snprintf(info,100,"Cos Pointing angle = %.3f ",fM->GetCosPointingAngle());
fa3634bd	310	TLatex* ltx = new TLatex(0.05, 0.9, info);
	311	ltx->SetTextSize(0.08);
	312	ltx->DrawLatex(0.05, 0.8, info);
f8d5a389	313
452b3436	314	snprintf(info,100,"p_{T} = %.3f [GeV/c]",fM->GetPt());
fa3634bd	315	ltx->DrawLatex(0.05, 0.7, info);
fa3634bd	316
452b3436	317	snprintf(info,100,"#eta = - ln( tan(#theta/2) ) = %.3f",fM->GetEta());
fa3634bd	318	ltx->DrawLatex(0.05, 0.6, info);
fa3634bd	319
452b3436	320	snprintf(info,100,"D^{0} inv. mass = %.3f", fM->GetInvariantMassPart());
f8d5a389	321	ltx->DrawLatex(0.05, 0.5, info);
fa3634bd	322
452b3436	323	snprintf(info,100,"D^{0}bar inv. mass = %.3f", fM->GetInvariantMassAntiPart());
f8d5a389	324	ltx->DrawLatex(0.05, 0.4, info);
fa3634bd	325
	326	gEve->Redraw3D();
	327	}