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

// $Id$
// Main authors: Matevz Tadel & Alja Mrak-Tadel: 2006, 2007

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

//------------------------------------------------------------------------
//
// AliEveCascade class
//
//------------------------------------------------------------------------

//#include <TPolyLine3D.h>
//#include <TColor.h>

//#include <vector>

#include <TDatabasePDG.h>
#include <TEveTrackPropagator.h>
#include <TEveManager.h>

#include "AliEveCascade.h"


ClassImp(AliEveCascade)

//______________________________________________________________________________
AliEveCascade::AliEveCascade() :
  TEvePointSet(),

  fRecBirthV(),
  fRecDecayV(),
  fRecDecayP(),
  fRecDecayV0(),
  fBacTrack(0),
  fNegTrack(0),
  fPosTrack(0),
  fRnrStyle(0),
  fPointingCurve(0),
  fV0Path(0),
  fESDIndex(-1),
  fDaughterDCA(999),
  fChi2Cascade(-1),
  fLambdaP(0.,0.,0.),
  fBachP(0.,0.,0.)
{
  // Default constructor.

  // Override from TEveElement.
  fPickable = kTRUE;
  fMainColorPtr = &fMarkerColor;
}

//______________________________________________________________________________
AliEveCascade::AliEveCascade(TEveRecTrack* tBac, TEveRecTrack* tNeg, TEveRecTrack* tPos, TEveRecV0* v0, TEveRecCascade* cascade, TEveTrackPropagator* rs) :
  TEvePointSet(),

  fRecBirthV(cascade->fCascadeBirth),
  fRecDecayV(cascade->fCascadeVCa),
  fRecDecayP(cascade->fPBac + v0->fPNeg + v0->fPPos),
  fRecDecayV0(v0->fVCa),
  fBacTrack(new AliEveTrack(tBac, rs)),
  fNegTrack(new AliEveTrack(tNeg, rs)),
  fPosTrack(new AliEveTrack(tPos, rs)),
  

  fRnrStyle(rs),
  fPointingCurve(new TEveLine("Pointing Curve")),
  fV0Path(new TEveLine("V0 Path")),
  fESDIndex(-1),
  fDaughterDCA(999), // DCA between Xi daughters is properly set in esd_make_cascade via the setter
  fChi2Cascade(-1),
  fLambdaP(0.,0.,0.), // Momemtum of Lambda is properly set in esd_make_cascade via the setter
  fBachP(0.,0.,0.)    // Momemtum of Lambda is properly set in esd_make_cascade via the setter
{
  // Constructor with full Cascade specification.

  // Override from TEveElement.
  fPickable = kTRUE;
  fMainColorPtr = &fMarkerColor;

  fMarkerStyle = 2;
  fMarkerColor = kMagenta+2;
  fMarkerSize  = 1;

  fPointingCurve->SetTitle("Xi pointing curve");
  fPointingCurve->SetLineColor(fMarkerColor);
  fPointingCurve->SetLineStyle(9);
  fPointingCurve->SetLineWidth(2);  
  fPointingCurve->IncDenyDestroy();
  AddElement(fPointingCurve);

  fV0Path->SetTitle("V0 path");
  fV0Path->SetLineColor(fMarkerColor);
  fV0Path->SetLineStyle(3);
  fV0Path->SetLineWidth(2);
  fV0Path->IncDenyDestroy();
  AddElement(fV0Path);

  fBacTrack->SetLineColor(kMagenta);
  fBacTrack->SetStdTitle();
  fBacTrack->IncDenyDestroy();
  AddElement(fBacTrack);
  
  fNegTrack->SetLineColor(kCyan+2);  // in V0 = kCyan
  fNegTrack->SetStdTitle();
  fNegTrack->IncDenyDestroy();
  AddElement(fNegTrack);
  
  fPosTrack->SetLineColor(kRed+2);  // in Eve V0 = kRed
  fPosTrack->SetStdTitle();
  fPosTrack->IncDenyDestroy();
  AddElement(fPosTrack);  
}

//______________________________________________________________________________
AliEveCascade::~AliEveCascade()
{
  // Destructor. Dereferences bachelor track and pointing-line objects.

  fBacTrack->DecDenyDestroy();
  fNegTrack->DecDenyDestroy();
  fPosTrack->DecDenyDestroy();
  fPointingCurve->DecDenyDestroy();
  fV0Path->DecDenyDestroy();
}

//______________________________________________________________________________
void AliEveCascade::MakeCascade()
{
  // Set all dependant components for drawing.

  SetPoint(0, fRecDecayV.fX, fRecDecayV.fY, fRecDecayV.fZ);

  fBacTrack->MakeTrack();  
  fNegTrack->MakeTrack();
  fPosTrack->MakeTrack();

  fPointingCurve->SetPoint(0, fRecBirthV.fX, fRecBirthV.fY, fRecBirthV.fZ);
  fPointingCurve->SetPoint(1, fRecDecayV.fX, fRecDecayV.fY, fRecDecayV.fZ);

  fV0Path->SetPoint(0, fRecDecayV.fX, fRecDecayV.fY, fRecDecayV.fZ);
  fV0Path->SetPoint(1, fRecDecayV0.fX, fRecDecayV0.fY, fRecDecayV0.fZ);
}


//______________________________________________________________________________
Double_t AliEveCascade::GetInvMass(Int_t cascadePdgCodeHyp) const
{
  // Returns Invariant Mass assuming the masses of the daughter particles
	TEveVector lBachMomentum = fBacTrack->GetMomentum();
  // Does not work properly because momenta at the primary vertex for V0 daughters (same for AliEVEv0) + Bach
	
	Double_t lBachMass = 0.0;
	if(cascadePdgCodeHyp == kXiMinus || cascadePdgCodeHyp == -kXiMinus )
		lBachMass = TDatabasePDG::Instance()->GetParticle("pi-")->Mass();
	else if(cascadePdgCodeHyp == kOmegaMinus || cascadePdgCodeHyp == -kOmegaMinus )
		lBachMass = TDatabasePDG::Instance()->GetParticle("K-")->Mass();
	
	Double_t lLambdaMass = TDatabasePDG::Instance()->GetParticle("Lambda0")->Mass();
	
	printf("\n Cascade : check the mass of the bachelor (%.5f) - the Lambda (%.5f)\n",lBachMass, lLambdaMass);
	
	Double_t eBach   = TMath::Sqrt(lBachMass*lBachMass     + lBachMomentum.Mag2());
	//Double_t eBach   = TMath::Sqrt(lBachMass  *lBachMass   + fBachP.Mag2());
	Double_t eLambda = TMath::Sqrt(lLambdaMass*lLambdaMass + fLambdaP.Mag2());
	
	TVector3 lCascadeMom(fBachP + fLambdaP);
	/*
	printf("Ptot(Lambda) : %f GeV/c \n", fLambdaP.Mag());
	printf("Ptot(Bach)   : %f GeV/c \n", lBachMomentum.Mag());
	printf("Ptot(BachESD): %f GeV/c \n", fBachP.Mag());
	printf("Ptot(Casc)   : %f GeV/c \n", fRecDecayP.Mag());
	printf("Ptot(CascESD): %f GeV/c \n", lCascadeMom.Mag());
	*/
	//return TMath::Sqrt( (eBach+eLambda)*(eBach+eLambda) - fRecDecayP.Mag2() );	// numerically instable ...
	//return TMath::Sqrt( ((eBach+eLambda) - lCascadeMom.Mag()) * ((eBach+eLambda) + lCascadeMom.Mag()) );
	return TMath::Sqrt( ((eBach+eLambda) - fRecDecayP.Mag()) * ((eBach+eLambda) + fRecDecayP.Mag()) );
}


//------------------------------------------------------------------------
//
// AliEveCascadeList class
//
//------------------------------------------------------------------------

ClassImp(AliEveCascadeList)

//______________________________________________________________________________
AliEveCascadeList::AliEveCascadeList() :
  TEveElementList(),
  fTitle(),
  fRnrStyle(0),
  fRnrDaughters(kTRUE),
  fRnrCascadevtx(kTRUE),
  fRnrCascadepath(kTRUE),
  fBacColor(0),
  fMinRCut(0),
  fMaxRCut(250),
  fMinDaughterDCA(0),
  fMaxDaughterDCA(1),
  fMinPt(0),
  fMaxPt(20),
  fInvMassHyp(kXiMinus),
  fMinInvariantMass(1.0),
  fMaxInvariantMass(6.0)
{
  // Default constructor.

  fChildClass = AliEveCascade::Class(); // override member from base TEveElementList
}

//______________________________________________________________________________
AliEveCascadeList::AliEveCascadeList(TEveTrackPropagator* rs) :
  TEveElementList(),
  fTitle(),
  fRnrStyle(rs),
  fRnrDaughters(kTRUE),
  fRnrCascadevtx(kTRUE),
  fRnrCascadepath(kTRUE),
  fBacColor(0),
  fMinRCut(0),
  fMaxRCut(250),
  fMinDaughterDCA(0),
  fMaxDaughterDCA(1),
  fMinPt(0),
  fMaxPt(20),
  fInvMassHyp(kXiMinus),
  fMinInvariantMass(1.0),
  fMaxInvariantMass(6.0)
{
  // Constructor with given track-propagator..

  fChildClass = AliEveCascade::Class(); // override member from base TEveElementList

  Init();
}

//______________________________________________________________________________
AliEveCascadeList::AliEveCascadeList(const Text_t* name, TEveTrackPropagator* rs) :
  TEveElementList(),
  fTitle(),
  fRnrStyle(rs),
  fRnrDaughters(kTRUE),
  fRnrCascadevtx(kTRUE),
  fRnrCascadepath(kTRUE),
  fBacColor(0),
  fMinRCut(0),
  fMaxRCut(100),
  fMinDaughterDCA(0),
  fMaxDaughterDCA(1),
  fMinPt(0),
  fMaxPt(20),
  fInvMassHyp(kXiMinus),
  fMinInvariantMass(1.0),
  fMaxInvariantMass(6.0)
{
  // Standard constructor.

  fChildClass = AliEveCascade::Class(); // override member from base TEveElementList

  Init();
  SetName(name);
}

//______________________________________________________________________________
void AliEveCascadeList::Init()
{
  // Initialize members needed for drawing operations.

  if (fRnrStyle== 0) fRnrStyle = new TEveTrackPropagator;
}

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

//______________________________________________________________________________
void AliEveCascadeList::MakeCascades()
{
  // Call MakeCascade() for all elements.

  for(List_i i=fChildren.begin(); i!=fChildren.end(); ++i) {
    ((AliEveCascade*)(*i))->MakeCascade();
  }
  gEve->Redraw3D();
}

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

//______________________________________________________________________________
void AliEveCascadeList::FilterByRadius(Float_t minR, Float_t maxR)
{
  // Select visibility of elements based on their axial radius.

  fMinRCut = minR;
  fMaxRCut = maxR;

  for(List_i i = fChildren.begin(); i != fChildren.end(); ++i)
  {
    AliEveCascade* cascade = (AliEveCascade*) *i;
    Float_t  rad = cascade->GetRadius();
    Bool_t  show = rad >= fMinRCut && rad <= fMaxRCut;
    cascade->SetRnrState(show);
  }
  ElementChanged();
  gEve->Redraw3D();
}

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

//______________________________________________________________________________
void AliEveCascadeList::FilterByDaughterDCA(Float_t minDaughterDCA, Float_t maxDaughterDCA)
{
  // Select visibility of elements based on the DCA between daughters.

  fMinDaughterDCA = minDaughterDCA;
  fMaxDaughterDCA = maxDaughterDCA;

  for(List_i i = fChildren.begin(); i != fChildren.end(); ++i)
  {
    AliEveCascade* cascade = (AliEveCascade*) *i;
    Float_t  dca = cascade->GetDaughterDCA();
    Bool_t  show = dca >= fMinDaughterDCA && dca <= fMaxDaughterDCA;
    cascade->SetRnrState(show);
  }
  ElementChanged();
  gEve->Redraw3D();
}

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

//______________________________________________________________________________
void AliEveCascadeList::FilterByPt(Float_t minPt, Float_t maxPt)
{
  // Select visibility of elements based on the Cascade pt.

  fMinPt = minPt;
  fMaxPt = maxPt;

  for(List_i i = fChildren.begin(); i != fChildren.end(); ++i)
  {
    AliEveCascade* cascade = (AliEveCascade*) *i;
    Float_t  pt = cascade->GetPt();
    Bool_t  show = pt >= fMinPt && pt <= fMaxPt;
    cascade->SetRnrState(show);
  }
  ElementChanged();
  gEve->Redraw3D();
}

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

//______________________________________________________________________________
void AliEveCascadeList::FilterByInvariantMass(Float_t minInvariantMass, Float_t maxInvariantMass, Int_t cascadePdgCodeHyp)
{
  // Select visibility of elements based on the Cascade invariant mass.

	fMinInvariantMass = minInvariantMass;
	fMaxInvariantMass = maxInvariantMass;

	for(List_i i = fChildren.begin(); i != fChildren.end(); ++i)
	{
		AliEveCascade* cascade = (AliEveCascade*) *i;
		Float_t  invMass = (Float_t)cascade->GetInvMass(cascadePdgCodeHyp);
		Bool_t   show    = invMass >= fMinInvariantMass && invMass <= fMaxInvariantMass;
		cascade->SetRnrState(show);
	}
	ElementChanged();
	gEve->Redraw3D();
}
Commit	Line	Data
2991d609	1	// $Id$
	2	// Main authors: Matevz Tadel & Alja Mrak-Tadel: 2006, 2007
	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
9851a663	10	//------------------------------------------------------------------------
	11	//
	12	// AliEveCascade class
	13	//
	14	//------------------------------------------------------------------------
	15
	16	//#include <TPolyLine3D.h>
	17	//#include <TColor.h>
	18
	19	//#include <vector>
2991d609	20
c180f65d	21	#include <TDatabasePDG.h>
2991d609	22	#include <TEveTrackPropagator.h>
	23	#include <TEveManager.h>
	24
9851a663	25	#include "AliEveCascade.h"
2991d609	26
2991d609	27
2991d609	28
2991d609	29
	30	ClassImp(AliEveCascade)
	31
	32	//______________________________________________________________________________
	33	AliEveCascade::AliEveCascade() :
	34	TEvePointSet(),
	35
	36	fRecBirthV(),
	37	fRecDecayV(),
	38	fRecDecayP(),
73ce6137	39	fRecDecayV0(),
2991d609	40	fBacTrack(0),
9851a663	41	fNegTrack(0),
9851a663	42	fPosTrack(0),
2991d609	43	fRnrStyle(0),
73ce6137	44	fPointingCurve(0),
73ce6137	45	fV0Path(0),
2991d609	46	fESDIndex(-1),
2991d609	47	fDaughterDCA(999),
9851a663	48	fChi2Cascade(-1),
	49	fLambdaP(0.,0.,0.),
	50	fBachP(0.,0.,0.)
2991d609	51	{
	52	// Default constructor.
	53
	54	// Override from TEveElement.
	55	fPickable = kTRUE;
	56	fMainColorPtr = &fMarkerColor;
	57	}
	58
	59	//______________________________________________________________________________
9851a663	60	AliEveCascade::AliEveCascade(TEveRecTrack* tBac, TEveRecTrack* tNeg, TEveRecTrack* tPos, TEveRecV0* v0, TEveRecCascade* cascade, TEveTrackPropagator* rs) :
2991d609	61	TEvePointSet(),
	62
	63	fRecBirthV(cascade->fCascadeBirth),
	64	fRecDecayV(cascade->fCascadeVCa),
73ce6137	65	fRecDecayP(cascade->fPBac + v0->fPNeg + v0->fPPos),
73ce6137	66	fRecDecayV0(v0->fVCa),
0e33c639	67	fBacTrack(new AliEveTrack(tBac, rs)),
9851a663	68	fNegTrack(new AliEveTrack(tNeg, rs)),
	69	fPosTrack(new AliEveTrack(tPos, rs)),
	70
2991d609	71
2991d609	72	fRnrStyle(rs),
73ce6137	73	fPointingCurve(new TEveLine("Pointing Curve")),
73ce6137	74	fV0Path(new TEveLine("V0 Path")),
2991d609	75	fESDIndex(-1),
9851a663	76	fDaughterDCA(999), // DCA between Xi daughters is properly set in esd_make_cascade via the setter
	77	fChi2Cascade(-1),
	78	fLambdaP(0.,0.,0.), // Momemtum of Lambda is properly set in esd_make_cascade via the setter
	79	fBachP(0.,0.,0.) // Momemtum of Lambda is properly set in esd_make_cascade via the setter
2991d609	80	{
	81	// Constructor with full Cascade specification.
	82
	83	// Override from TEveElement.
	84	fPickable = kTRUE;
	85	fMainColorPtr = &fMarkerColor;
	86
	87	fMarkerStyle = 2;
9851a663	88	fMarkerColor = kMagenta+2;
2991d609	89	fMarkerSize = 1;
2991d609	90
9851a663	91	fPointingCurve->SetTitle("Xi pointing curve");
73ce6137	92	fPointingCurve->SetLineColor(fMarkerColor);
9851a663	93	fPointingCurve->SetLineStyle(9);
9851a663	94	fPointingCurve->SetLineWidth(2);
73ce6137	95	fPointingCurve->IncDenyDestroy();
	96	AddElement(fPointingCurve);
	97
9851a663	98	fV0Path->SetTitle("V0 path");
73ce6137	99	fV0Path->SetLineColor(fMarkerColor);
	100	fV0Path->SetLineStyle(3);
	101	fV0Path->SetLineWidth(2);
	102	fV0Path->IncDenyDestroy();
	103	AddElement(fV0Path);
2991d609	104
9851a663	105	fBacTrack->SetLineColor(kMagenta);
2991d609	106	fBacTrack->SetStdTitle();
2991d609	107	fBacTrack->IncDenyDestroy();
2991d609	108	AddElement(fBacTrack);
9851a663	109
	110	fNegTrack->SetLineColor(kCyan+2); // in V0 = kCyan
	111	fNegTrack->SetStdTitle();
	112	fNegTrack->IncDenyDestroy();
	113	AddElement(fNegTrack);
	114
	115	fPosTrack->SetLineColor(kRed+2); // in Eve V0 = kRed
	116	fPosTrack->SetStdTitle();
	117	fPosTrack->IncDenyDestroy();
	118	AddElement(fPosTrack);
2991d609	119	}
	120
	121	//______________________________________________________________________________
	122	AliEveCascade::~AliEveCascade()
	123	{
	124	// Destructor. Dereferences bachelor track and pointing-line objects.
	125
	126	fBacTrack->DecDenyDestroy();
9851a663	127	fNegTrack->DecDenyDestroy();
9851a663	128	fPosTrack->DecDenyDestroy();
73ce6137	129	fPointingCurve->DecDenyDestroy();
73ce6137	130	fV0Path->DecDenyDestroy();
2991d609	131	}
	132
	133	//______________________________________________________________________________
	134	void AliEveCascade::MakeCascade()
	135	{
	136	// Set all dependant components for drawing.
	137
	138	SetPoint(0, fRecDecayV.fX, fRecDecayV.fY, fRecDecayV.fZ);
	139
9851a663	140	fBacTrack->MakeTrack();
	141	fNegTrack->MakeTrack();
	142	fPosTrack->MakeTrack();
2991d609	143
73ce6137	144	fPointingCurve->SetPoint(0, fRecBirthV.fX, fRecBirthV.fY, fRecBirthV.fZ);
	145	fPointingCurve->SetPoint(1, fRecDecayV.fX, fRecDecayV.fY, fRecDecayV.fZ);
	146
	147	fV0Path->SetPoint(0, fRecDecayV.fX, fRecDecayV.fY, fRecDecayV.fZ);
	148	fV0Path->SetPoint(1, fRecDecayV0.fX, fRecDecayV0.fY, fRecDecayV0.fZ);
2991d609	149	}
	150
	151
9851a663	152
	153	//______________________________________________________________________________
	154	Double_t AliEveCascade::GetInvMass(Int_t cascadePdgCodeHyp) const
	155	{
	156	// Returns Invariant Mass assuming the masses of the daughter particles
	157	TEveVector lBachMomentum = fBacTrack->GetMomentum();
	158	// Does not work properly because momenta at the primary vertex for V0 daughters (same for AliEVEv0) + Bach
	159
	160	Double_t lBachMass = 0.0;
	161	if(cascadePdgCodeHyp == kXiMinus \|\| cascadePdgCodeHyp == -kXiMinus )
	162	lBachMass = TDatabasePDG::Instance()->GetParticle("pi-")->Mass();
	163	else if(cascadePdgCodeHyp == kOmegaMinus \|\| cascadePdgCodeHyp == -kOmegaMinus )
	164	lBachMass = TDatabasePDG::Instance()->GetParticle("K-")->Mass();
	165
	166	Double_t lLambdaMass = TDatabasePDG::Instance()->GetParticle("Lambda0")->Mass();
	167
	168	printf("\n Cascade : check the mass of the bachelor (%.5f) - the Lambda (%.5f)\n",lBachMass, lLambdaMass);
	169
	170	Double_t eBach = TMath::Sqrt(lBachMass*lBachMass + lBachMomentum.Mag2());
	171	//Double_t eBach = TMath::Sqrt(lBachMass *lBachMass + fBachP.Mag2());
	172	Double_t eLambda = TMath::Sqrt(lLambdaMass*lLambdaMass + fLambdaP.Mag2());
	173
	174	TVector3 lCascadeMom(fBachP + fLambdaP);
	175	/*
	176	printf("Ptot(Lambda) : %f GeV/c \n", fLambdaP.Mag());
	177	printf("Ptot(Bach) : %f GeV/c \n", lBachMomentum.Mag());
	178	printf("Ptot(BachESD): %f GeV/c \n", fBachP.Mag());
	179	printf("Ptot(Casc) : %f GeV/c \n", fRecDecayP.Mag());
	180	printf("Ptot(CascESD): %f GeV/c \n", lCascadeMom.Mag());
	181	*/
	182	//return TMath::Sqrt( (eBach+eLambda)*(eBach+eLambda) - fRecDecayP.Mag2() ); // numerically instable ...
	183	//return TMath::Sqrt( ((eBach+eLambda) - lCascadeMom.Mag()) * ((eBach+eLambda) + lCascadeMom.Mag()) );
	184	return TMath::Sqrt( ((eBach+eLambda) - fRecDecayP.Mag()) * ((eBach+eLambda) + fRecDecayP.Mag()) );
	185	}
	186
	187
	188
	189	//------------------------------------------------------------------------
	190	//
	191	// AliEveCascadeList class
	192	//
	193	//------------------------------------------------------------------------
2991d609	194
	195	ClassImp(AliEveCascadeList)
	196
	197	//______________________________________________________________________________
	198	AliEveCascadeList::AliEveCascadeList() :
	199	TEveElementList(),
	200	fTitle(),
	201	fRnrStyle(0),
	202	fRnrDaughters(kTRUE),
	203	fRnrCascadevtx(kTRUE),
	204	fRnrCascadepath(kTRUE),
	205	fBacColor(0),
	206	fMinRCut(0),
	207	fMaxRCut(250),
	208	fMinDaughterDCA(0),
	209	fMaxDaughterDCA(1),
	210	fMinPt(0),
9851a663	211	fMaxPt(20),
	212	fInvMassHyp(kXiMinus),
	213	fMinInvariantMass(1.0),
	214	fMaxInvariantMass(6.0)
2991d609	215	{
	216	// Default constructor.
	217
	218	fChildClass = AliEveCascade::Class(); // override member from base TEveElementList
	219	}
	220
	221	//______________________________________________________________________________
	222	AliEveCascadeList::AliEveCascadeList(TEveTrackPropagator* rs) :
	223	TEveElementList(),
	224	fTitle(),
	225	fRnrStyle(rs),
	226	fRnrDaughters(kTRUE),
	227	fRnrCascadevtx(kTRUE),
	228	fRnrCascadepath(kTRUE),
	229	fBacColor(0),
	230	fMinRCut(0),
	231	fMaxRCut(250),
	232	fMinDaughterDCA(0),
	233	fMaxDaughterDCA(1),
	234	fMinPt(0),
9851a663	235	fMaxPt(20),
	236	fInvMassHyp(kXiMinus),
	237	fMinInvariantMass(1.0),
	238	fMaxInvariantMass(6.0)
2991d609	239	{
	240	// Constructor with given track-propagator..
	241
	242	fChildClass = AliEveCascade::Class(); // override member from base TEveElementList
	243
	244	Init();
	245	}
	246
	247	//______________________________________________________________________________
	248	AliEveCascadeList::AliEveCascadeList(const Text_t* name, TEveTrackPropagator* rs) :
	249	TEveElementList(),
	250	fTitle(),
	251	fRnrStyle(rs),
	252	fRnrDaughters(kTRUE),
	253	fRnrCascadevtx(kTRUE),
	254	fRnrCascadepath(kTRUE),
	255	fBacColor(0),
	256	fMinRCut(0),
	257	fMaxRCut(100),
	258	fMinDaughterDCA(0),
	259	fMaxDaughterDCA(1),
	260	fMinPt(0),
9851a663	261	fMaxPt(20),
	262	fInvMassHyp(kXiMinus),
	263	fMinInvariantMass(1.0),
	264	fMaxInvariantMass(6.0)
2991d609	265	{
	266	// Standard constructor.
	267
	268	fChildClass = AliEveCascade::Class(); // override member from base TEveElementList
	269
	270	Init();
	271	SetName(name);
	272	}
	273
	274	//______________________________________________________________________________
	275	void AliEveCascadeList::Init()
	276	{
	277	// Initialize members needed for drawing operations.
	278
	279	if (fRnrStyle== 0) fRnrStyle = new TEveTrackPropagator;
	280	}
	281
	282	/******************************************************************************/
	283
	284	//______________________________________________________________________________
	285	void AliEveCascadeList::MakeCascades()
	286	{
	287	// Call MakeCascade() for all elements.
	288
	289	for(List_i i=fChildren.begin(); i!=fChildren.end(); ++i) {
	290	((AliEveCascade)(i))->MakeCascade();
	291	}
	292	gEve->Redraw3D();
	293	}
	294
	295	/******************************************************************************/
	296
	297	//______________________________________________________________________________
	298	void AliEveCascadeList::FilterByRadius(Float_t minR, Float_t maxR)
	299	{
	300	// Select visibility of elements based on their axial radius.
	301
	302	fMinRCut = minR;
	303	fMaxRCut = maxR;
	304
	305	for(List_i i = fChildren.begin(); i != fChildren.end(); ++i)
	306	{
	307	AliEveCascade* cascade = (AliEveCascade) i;
	308	Float_t rad = cascade->GetRadius();
	309	Bool_t show = rad >= fMinRCut && rad <= fMaxRCut;
	310	cascade->SetRnrState(show);
	311	}
	312	ElementChanged();
	313	gEve->Redraw3D();
	314	}
	315
	316	/******************************************************************************/
	317
	318	//______________________________________________________________________________
	319	void AliEveCascadeList::FilterByDaughterDCA(Float_t minDaughterDCA, Float_t maxDaughterDCA)
	320	{
	321	// Select visibility of elements based on the DCA between daughters.
	322
	323	fMinDaughterDCA = minDaughterDCA;
	324	fMaxDaughterDCA = maxDaughterDCA;
	325
	326	for(List_i i = fChildren.begin(); i != fChildren.end(); ++i)
	327	{
	328	AliEveCascade* cascade = (AliEveCascade) i;
329	Float_t dca = cascade->GetDaughterDCA();
330	Bool_t show = dca >= fMinDaughterDCA && dca <= fMaxDaughterDCA;
331	cascade->SetRnrState(show);
332	}
333	ElementChanged();
334	gEve->Redraw3D();
335	}
336
337	/******************************************************************************/
338
339	//______________________________________________________________________________
340	void AliEveCascadeList::FilterByPt(Float_t minPt, Float_t maxPt)
341	{
342	// Select visibility of elements based on the Cascade pt.
343
344	fMinPt = minPt;
345	fMaxPt = maxPt;
346
347	for(List_i i = fChildren.begin(); i != fChildren.end(); ++i)
348	{
349	AliEveCascade* cascade = (AliEveCascade) i;
350	Float_t pt = cascade->GetPt();
351	Bool_t show = pt >= fMinPt && pt <= fMaxPt;
352	cascade->SetRnrState(show);
353	}
354	ElementChanged();
355	gEve->Redraw3D();
356	}
9851a663	357
	358	/******************************************************************************/
	359
	360	//______________________________________________________________________________
	361	void AliEveCascadeList::FilterByInvariantMass(Float_t minInvariantMass, Float_t maxInvariantMass, Int_t cascadePdgCodeHyp)
	362	{
	363	// Select visibility of elements based on the Cascade invariant mass.
	364
	365	fMinInvariantMass = minInvariantMass;
	366	fMaxInvariantMass = maxInvariantMass;
	367
	368	for(List_i i = fChildren.begin(); i != fChildren.end(); ++i)
	369	{
	370	AliEveCascade* cascade = (AliEveCascade) i;
	371	Float_t invMass = (Float_t)cascade->GetInvMass(cascadePdgCodeHyp);
	372	Bool_t show = invMass >= fMinInvariantMass && invMass <= fMaxInvariantMass;
	373	cascade->SetRnrState(show);
	374	}
	375	ElementChanged();
	376	gEve->Redraw3D();
	377	}
	378