#include <RhoSimpleVertexSelector.h>

Inheritance diagram for RhoSimpleVertexSelector:

Public Member Functions
	RhoSimpleVertexSelector (const char *name="SimpleVertexSelector", Double_t d=1.0, Double_t a=3.14159265358979323846, Double_t r1=0.0, Double_t r2=1.E8)

virtual	~RhoSimpleVertexSelector ()

virtual Bool_t	Accept (RhoCandidate &a, RhoCandidate &b)

virtual Bool_t	Accept (RhoCandidate , RhoCandidate )=0

virtual Bool_t	Accept (RhoCandidate , RhoCandidate , RhoCandidate *)

virtual Bool_t	Accept (RhoCandidate , RhoCandidate , RhoCandidate , RhoCandidate )

virtual Bool_t	Accept (RhoCandidate , RhoCandidate , RhoCandidate , RhoCandidate , RhoCandidate *)

void	SetDistanceOfClosestApproach (Double_t d)

void	SetAngle (Double_t d)

void	SetRadius (Double_t r1, Double_t r2=1.E8)

void	SetPrimaryVertex (TVector3 x)

TVector3	GetVertex ()

Double_t	GetDistanceOfClosestApproach ()

Double_t	GetAngle ()

virtual void	PrintOn (std::ostream &o=std::cout) const

const TVector3 &	GetMomentumA () const

const TVector3 &	GetMomentumB () const

TVector3	GetMomentum () const

Double_t	GetRadius () const

void	Combine (RhoCandList &in1, RhoCandList &in2, RhoCandList &out)

Protected Attributes
Double_t	fDoca

Double_t	fVtxip

Double_t	fRmin

Double_t	fRmax

Double_t	fActualR

Double_t	fActualVtxip

Double_t	fActualDoca

TVector3	fVertex

TVector3	fPrimaryVertex

TVector3	fMomA

TVector3	fMomB

Detailed Description

Definition at line 25 of file RhoSimpleVertexSelector.h.

Constructor & Destructor Documentation

RhoSimpleVertexSelector::RhoSimpleVertexSelector	(	const char *	name = `"SimpleVertexSelector"`,
		Double_t	d = `1.0`,
		Double_t	a = `3.14159265358979323846`,
		Double_t	r1 = `0.0`,
		Double_t	r2 = `1.E8`
	)

Definition at line 33 of file RhoSimpleVertexSelector.cxx.

References a, d, RhoVertexSelectorBase::fDoca, RhoVertexSelectorBase::fRmax, RhoVertexSelectorBase::fRmin, RhoVertexSelectorBase::fVtxip, r1, and r2.

                                                                                                                       :
   RhoVertexSelectorBase ( name ) //, fQC(0)
 {
   fDoca = d;
   fVtxip = a;
   fRmin = r1;
   fRmax = r2;
 }

RhoSimpleVertexSelector::~RhoSimpleVertexSelector ( )

virtual

Definition at line 42 of file RhoSimpleVertexSelector.cxx.

42 {};// delete fQC; }

Member Function Documentation

Bool_t RhoSimpleVertexSelector::Accept	(	RhoCandidate &	a,
		RhoCandidate &	b
	)

virtual

Definition at line 44 of file RhoSimpleVertexSelector.cxx.

References best, RhoCandidate::Charge(), Double_t, RhoVertexSelectorBase::fActualDoca, RhoVertexSelectorBase::fActualR, RhoVertexSelectorBase::fActualVtxip, RhoVertexSelectorBase::fDoca, RhoVertexSelectorBase::fMomA, RhoVertexSelectorBase::fMomB, RhoVertexSelectorBase::fPrimaryVertex, RhoVertexSelectorBase::fRmax, RhoVertexSelectorBase::fRmin, RhoVertexSelectorBase::fVertex, RhoVertexSelectorBase::fVtxip, RhoCandidate::GetPosition(), p, RhoCandidate::P3(), pnt, r1, r2, sign(), sqrt(), x, and y.

 {
   if ( &a==0 || &b==0 ) { return kFALSE; }
 
   // Position vectors
   TVector3 position1 = a.GetPosition();
   TVector3 position2 = b.GetPosition();
 
   double pnt[3], Bf[3];
   pnt[0]=0.5* ( position1.X() +position2.X() );
   pnt[1]=0.5* ( position1.Y() +position2.Y() );
   pnt[2]=0.5* ( position1.Z() +position2.Z() );
   if(FairRun::Instance()->IsAna()){
     FairRunAna::Instance()->GetField()->GetFieldValue ( pnt, Bf ); //[kGs]
   } else{
     FairRunSim::Instance()->GetField()->GetFieldValue ( pnt, Bf ); //[kGs]
   }
   Float_t bField = Bf[2]; // Retrieve the B-Field
 
   // Momentum vectors
   TVector3 ap3 = a.P3();
   Double_t pPerp1 = ap3.Perp();
   TVector3 d1 = ap3;
   d1.SetZ ( 0 );
   d1*=1.0/pPerp1;
 
   TVector3 bp3 = b.P3();
   Double_t pPerp2 = bp3.Perp();
   TVector3 d2 = bp3;
   d2.SetZ ( 0 );
   d2*=1.0/pPerp2;
 
 
   TVector3 dB ( 0,0,1.0 );
   // Radius and center
   Double_t rho1 = 100. * pPerp1/ ( 0.3*bField ); // Radius in cm //FIXME really?
   TVector3 r1=d1.Cross ( dB );
   r1 *= -a.Charge() *rho1;
   TVector3 center1 = position1 - r1;
   center1.SetZ ( 0 );
 
   Double_t rho2 = 100. * pPerp2/ ( 0.3*bField ); // Radius in cm //FIXME really?
   TVector3 r2=d2.Cross ( dB );
   r2 *= -b.Charge() *rho2;
   TVector3 center2 = position2 - r2;
   center2.SetZ ( 0 );
 
   // distance and angle of the axis between the two centers
   TVector3 ab = center2 - center1;
   Double_t dab = ab.Perp();
   Double_t cosTheAB = ab.X() /dab;
   Double_t sinTheAB = ab.Y() /dab;
 
 
   // x value of intersect at reduced system
   Double_t x = dab/2 + ( rho1*rho1 - rho2*rho2 ) / ( 2*dab );
 
   // y*y value of intersect at reduced system for helix A
   Double_t y2 = ( rho1+x ) * ( rho1-x );
 
   // both circles do not intersect (only one solution)
   Int_t nSolMax=1;
   Double_t y=0;
 
   if ( y2 > 0 ) {
     nSolMax=2;
     y = sqrt ( y2 );
   } else {
     //if( fabs(dab-rho1-rho2) > fActualDoca);
     //return kFALSE;
   }
 
   // now we compute the solution(s)
 
   TVector3 newap3[2];
   TVector3 newbp3[2];
   TVector3 newapos[2];
   TVector3 newbpos[2];
   Int_t best=0;
   fActualDoca=1.E8;
   for ( Int_t ns=0; ns<nSolMax; ns++ ) {   // loop on the solutions
 
     // radius vector of intersection point
     Double_t sign = ns ? 1.0 : -1.0;
     TVector3 rs1 ( cosTheAB*x - sinTheAB*y * sign, sinTheAB*x + cosTheAB*y * sign, 0 );
     TVector3 rs2 ( rs1-ab );
 
     // are we moving forward or backward?
     Double_t adir= ( rs1-r1 ).Dot ( ap3 ) >0 ? 1.0 : -1.0;
     Double_t aangle=adir * r1.Angle ( rs1 );
     // intersection point
     Double_t newaz=position1.Z() + rho1*aangle/pPerp1 * ap3.Z();
     newapos[ns].SetX ( center1.X() + rs1.X() );
     newapos[ns].SetY ( center1.Y() + rs1.Y() );
     newapos[ns].SetZ ( newaz );
 
     // adjust momentum
     newap3[ns]=rs1.Cross ( dB );
     newap3[ns]*=a.Charge() /rho1*pPerp1;
     newap3[ns].SetZ ( ap3.Z() );
 
 
     // same for b
     Double_t bdir= ( rs2-r2 ).Dot ( bp3 ) >0 ? 1.0 : -1.0;
     Double_t bangle=bdir * r2.Angle ( rs2 );
     Double_t newbz=position2.Z() + rho2*bangle/pPerp2 * bp3.Z();
     newbpos[ns].SetX ( center2.X() + rs2.X() ); // ==newapos[ns].X()
     newbpos[ns].SetY ( center2.Y() + rs2.Y() ); // ==newapos[ns].Y()
     newbpos[ns].SetZ ( newbz );
     newbp3[ns]=rs2.Cross ( dB );
     newbp3[ns]*=b.Charge() /rho2*pPerp2;
     newbp3[ns].SetZ ( bp3.Z() );
 
     Double_t delta = ( newapos[ns]-newbpos[ns] ).Mag();
 
     // take the solution of minimal deltaZ
     if ( delta < fActualDoca ) {
       best=ns;
       fActualDoca  = delta;
     }
   }
 
   fVertex=0.5* ( newapos[best]+newbpos[best] );
   fMomA=newap3[best];
   fMomB=newbp3[best];
 
   // Now reconstruct original flight path and
   // calculate the angle to vertex-ip
 
   TVector3 p=fMomA+fMomB;
   TVector3 vtxip = fVertex-fPrimaryVertex;
   fActualVtxip = vtxip.Angle ( p );
   fActualR = vtxip.Mag();
 
   //if (fQC!=0) fQC->Accumulate((Float_t)fActualDoca,(Float_t)fActualVtxip);
 
   if ( fActualDoca < fDoca &&
        fActualR >= fRmin && fActualR < fRmax &&
        fActualVtxip < fVtxip ) {
 
     return kTRUE;
 
   } else {
     fMomA=ap3;
     fMomB=bp3;
     fVertex = fPrimaryVertex;
     return kFALSE;
   }
   return kFALSE; //get rid of warnings
 }

virtual Bool_t RhoVertexSelectorBase::Accept	(	RhoCandidate *	,
		RhoCandidate *
	)

pure virtualinherited

Referenced by RhoVertexSelectorBase::Combine().

virtual Bool_t RhoVertexSelectorBase::Accept	(	RhoCandidate *	,
		RhoCandidate *	,
		RhoCandidate *
	)

inlinevirtualinherited

Definition at line 35 of file RhoVertexSelectorBase.h.

                                                                         {
       return kTRUE;
     };

virtual Bool_t RhoVertexSelectorBase::Accept	(	RhoCandidate *	,
		RhoCandidate *	,
		RhoCandidate *	,
		RhoCandidate *
	)

inlinevirtualinherited

Definition at line 38 of file RhoVertexSelectorBase.h.

                                                                                       {
       return kTRUE;
     };

virtual Bool_t RhoVertexSelectorBase::Accept	(	RhoCandidate *	,
		RhoCandidate *	,
		RhoCandidate *	,
		RhoCandidate *	,
		RhoCandidate *
	)

inlinevirtualinherited

Definition at line 41 of file RhoVertexSelectorBase.h.

                                                                                                     {
       return kTRUE;
     };

void RhoVertexSelectorBase::Combine	(	RhoCandList &	in1,
		RhoCandList &	in2,
		RhoCandList &	out
	)

inherited

Definition at line 32 of file RhoVertexSelectorBase.cxx.

References RhoVertexSelectorBase::Accept(), Bool_t, c, RhoCandList::Cleanup(), Double_t, RhoCandidate::E(), RhoCandList::GetLength(), RhoVertexSelectorBase::GetMomentum(), RhoVertexSelectorBase::GetVertex(), RhoCandList::Put(), RhoCandidate::SetE(), RhoCandidate::SetMarker(), RhoCandidate::SetPosition(), and RhoCandidate::SetVect().

 {
   TLorentzVector vl;
   Double_t charge;
   Bool_t nearby = kTRUE;
   //fEntries = 0;
   //end of unlock
   Int_t endpos1 = l1.GetLength();
   Int_t endpos2 = l2.GetLength();
   out.Cleanup();
 
   if ( &l1==&l2 ) {
     //combination of a list with itself
     for ( Int_t comb_i=0; comb_i<endpos1; ++comb_i ) {
       for ( Int_t comb_k=comb_i+1; comb_k<endpos2; ++comb_k ) {
         if ( !l1[comb_i]->Overlaps ( l2[comb_k] ) ) {
           vl=l1[comb_i]->P4() +l2[comb_k]->P4();
           charge=l1[comb_i]->Charge() +l2[comb_k]->Charge();
           nearby = Accept ( l1[comb_i],l2[comb_k] );
           if ( nearby ) {
             //fill list with new candidate
             RhoCandidate c ( vl,charge );
             c.SetMarker ( l1[comb_i]->GetMarker ( 0 ) |l2[comb_k]->GetMarker ( 0 ),0 );
             c.SetMarker ( l1[comb_i]->GetMarker ( 1 ) |l2[comb_k]->GetMarker ( 1 ),1 );
             c.SetMarker ( l1[comb_i]->GetMarker ( 2 ) |l2[comb_k]->GetMarker ( 2 ),2 );
             c.SetMarker ( l1[comb_i]->GetMarker ( 3 ) |l2[comb_k]->GetMarker ( 3 ),3 );
             c.SetPosition ( GetVertex() );
             c.SetVect ( GetMomentum() );
             c.SetE ( c.E() );
             out.Put ( &c );
           }
         }
       }
     }
   } else {
     //combination of 2 different lists
     for ( Int_t comb_i=0; comb_i<endpos1; ++comb_i ) {
       for ( Int_t comb_k=0; comb_k<endpos2; ++comb_k ) {
         if ( !l1[comb_i]->Overlaps ( l2[comb_k] ) ) {
           vl=l1[comb_i]->P4() +l2[comb_k]->P4();
           charge=l1[comb_i]->Charge() +l2[comb_k]->Charge();
           nearby = Accept ( l1[comb_i],l2[comb_k] );
           if ( nearby ) {
             //fill list with new candidate
             RhoCandidate c ( vl,charge );
             c.SetMarker ( l1[comb_i]->GetMarker ( 0 ) |l2[comb_k]->GetMarker ( 0 ),0 );
             c.SetMarker ( l1[comb_i]->GetMarker ( 1 ) |l2[comb_k]->GetMarker ( 1 ),1 );
             c.SetMarker ( l1[comb_i]->GetMarker ( 2 ) |l2[comb_k]->GetMarker ( 2 ),2 );
             c.SetMarker ( l1[comb_i]->GetMarker ( 3 ) |l2[comb_k]->GetMarker ( 3 ),3 );
             c.SetPosition ( GetVertex() );
             c.SetVect ( GetMomentum() );
             c.SetE ( c.E() );
             out.Put ( &c );
           }
         }
       }
     }
   }
 }

Double_t RhoVertexSelectorBase::GetAngle ( )

inlineinherited

Definition at line 63 of file RhoVertexSelectorBase.h.

References RhoVertexSelectorBase::fActualVtxip.

                         {
       return fActualVtxip;
     }

Double_t RhoVertexSelectorBase::GetDistanceOfClosestApproach ( )

inlineinherited

Definition at line 60 of file RhoVertexSelectorBase.h.

References RhoVertexSelectorBase::fActualDoca.

                                             {
       return fActualDoca;
     }

TVector3 RhoVertexSelectorBase::GetMomentum ( ) const

inlineinherited

Definition at line 74 of file RhoVertexSelectorBase.h.

References RhoVertexSelectorBase::fMomA, and RhoVertexSelectorBase::fMomB.

Referenced by RhoVertexSelectorBase::Combine().

                                  {
       return fMomA+fMomB;
     }

const TVector3& RhoVertexSelectorBase::GetMomentumA ( ) const

inlineinherited

Definition at line 68 of file RhoVertexSelectorBase.h.

References RhoVertexSelectorBase::fMomA.

                                          {
       return fMomA;
     }

const TVector3& RhoVertexSelectorBase::GetMomentumB ( ) const

inlineinherited

Definition at line 71 of file RhoVertexSelectorBase.h.

References RhoVertexSelectorBase::fMomB.

                                          {
       return fMomB;
     }

Double_t RhoVertexSelectorBase::GetRadius ( ) const

inlineinherited

Definition at line 77 of file RhoVertexSelectorBase.h.

References RhoVertexSelectorBase::fActualR.

                                {
       return fActualR;
     }

TVector3 RhoVertexSelectorBase::GetVertex ( )

inlineinherited

Definition at line 57 of file RhoVertexSelectorBase.h.

References RhoVertexSelectorBase::fVertex.

Referenced by RhoVertexSelectorBase::Combine().

                          {
       return fVertex;
     }

void RhoVertexSelectorBase::PrintOn ( std::ostream & o = std::cout ) const

virtualinherited

Definition at line 92 of file RhoVertexSelectorBase.cxx.

References RhoVertexSelectorBase::fDoca, RhoVertexSelectorBase::fRmax, RhoVertexSelectorBase::fRmin, and RhoVertexSelectorBase::fVtxip.

Referenced by operator<<().

 {
   o << GetName() << " settings:" << endl;
   o << "Distance of closest approach < " << fDoca << " cm" << endl;
   o << "Angle vertex-IP->momentum    < " << fVtxip << " rad" << endl;
   o << "Radius of acceptance         = " << fRmin << " ... " << fRmax << " cm" << endl;
   o << endl;
 }

void RhoVertexSelectorBase::SetAngle ( Double_t d )

inlineinherited

Definition at line 47 of file RhoVertexSelectorBase.h.

References d, and RhoVertexSelectorBase::fVtxip.

                                  {
       fVtxip = d;
     }

void RhoVertexSelectorBase::SetDistanceOfClosestApproach ( Double_t d )

inlineinherited

Definition at line 44 of file RhoVertexSelectorBase.h.

References d, and RhoVertexSelectorBase::fDoca.

                                                      {
       fDoca = d;
     }

void RhoVertexSelectorBase::SetPrimaryVertex ( TVector3 x )

inlineinherited

Definition at line 54 of file RhoVertexSelectorBase.h.

References RhoVertexSelectorBase::fPrimaryVertex, and x.

                                          {
       fPrimaryVertex = x;
     }

void RhoVertexSelectorBase::SetRadius	(	Double_t	r1,
		Double_t	r2 = `1.E8`
	)

inlineinherited

Definition at line 50 of file RhoVertexSelectorBase.h.

References RhoVertexSelectorBase::fRmax, RhoVertexSelectorBase::fRmin, r1, and r2.

                                                        {
       fRmin = r1;
       fRmax = r2;
     }

Member Data Documentation

Double_t RhoVertexSelectorBase::fActualDoca

protectedinherited

Definition at line 89 of file RhoVertexSelectorBase.h.

Referenced by Accept(), and RhoVertexSelectorBase::GetDistanceOfClosestApproach().

Double_t RhoVertexSelectorBase::fActualR

protectedinherited

Definition at line 87 of file RhoVertexSelectorBase.h.

Referenced by Accept(), and RhoVertexSelectorBase::GetRadius().

Double_t RhoVertexSelectorBase::fActualVtxip

protectedinherited

Definition at line 88 of file RhoVertexSelectorBase.h.

Referenced by Accept(), and RhoVertexSelectorBase::GetAngle().

Double_t RhoVertexSelectorBase::fDoca

protectedinherited

Definition at line 83 of file RhoVertexSelectorBase.h.

Referenced by Accept(), RhoVertexSelectorBase::PrintOn(), RhoSimpleVertexSelector(), and RhoVertexSelectorBase::SetDistanceOfClosestApproach().

TVector3 RhoVertexSelectorBase::fMomA

protectedinherited

Definition at line 92 of file RhoVertexSelectorBase.h.

Referenced by Accept(), RhoVertexSelectorBase::GetMomentum(), and RhoVertexSelectorBase::GetMomentumA().

TVector3 RhoVertexSelectorBase::fMomB

protectedinherited

Definition at line 92 of file RhoVertexSelectorBase.h.

Referenced by Accept(), RhoVertexSelectorBase::GetMomentum(), and RhoVertexSelectorBase::GetMomentumB().

TVector3 RhoVertexSelectorBase::fPrimaryVertex

protectedinherited

Definition at line 91 of file RhoVertexSelectorBase.h.

Referenced by Accept(), and RhoVertexSelectorBase::SetPrimaryVertex().

Double_t RhoVertexSelectorBase::fRmax

protectedinherited

Definition at line 86 of file RhoVertexSelectorBase.h.

Referenced by Accept(), RhoVertexSelectorBase::PrintOn(), RhoSimpleVertexSelector(), and RhoVertexSelectorBase::SetRadius().

Double_t RhoVertexSelectorBase::fRmin

protectedinherited

Definition at line 85 of file RhoVertexSelectorBase.h.

Referenced by Accept(), RhoVertexSelectorBase::PrintOn(), RhoSimpleVertexSelector(), and RhoVertexSelectorBase::SetRadius().

TVector3 RhoVertexSelectorBase::fVertex

protectedinherited

Definition at line 90 of file RhoVertexSelectorBase.h.

Referenced by Accept(), and RhoVertexSelectorBase::GetVertex().

Double_t RhoVertexSelectorBase::fVtxip

protectedinherited

Definition at line 84 of file RhoVertexSelectorBase.h.

Referenced by Accept(), RhoVertexSelectorBase::PrintOn(), RhoSimpleVertexSelector(), and RhoVertexSelectorBase::SetAngle().

The documentation for this class was generated from the following files:

Public Member Functions

Protected Attributes

Detailed Description

Constructor & Destructor Documentation

Member Function Documentation

Member Data Documentation