#include <EdbShowerAlg.h>

Inheritance diagram for EdbShowerAlg:

Collaboration diagram for EdbShowerAlg:

[legend]

Public Member Functions
void	AddRecoShowerArray (TObjArray *RecoShowerArray)

EdbVertex *	CalcVertex (TObjArray *segments)

	ClassDef (EdbShowerAlg, 1)

void	Convert_EdbPVRec_To_InBTArray ()

Double_t	DeltaR_NoPropagation (EdbSegP s, EdbSegP stest)

Double_t	DeltaR_WithoutPropagation (EdbSegP s, EdbSegP stest)

Double_t	DeltaR_WithPropagation (EdbSegP s, EdbSegP stest)

Double_t	DeltaTheta (EdbSegP s1, EdbSegP s2)

Double_t	DeltaThetaComponentwise (EdbSegP s1, EdbSegP s2)

Double_t	DeltaThetaSingleAngles (EdbSegP s1, EdbSegP s2)

	EdbShowerAlg ()

	EdbShowerAlg (TString AlgName, Int_t AlgValue)

virtual void	Execute ()

virtual void	Finalize ()

TString	GetAlgName () const

Int_t	GetAlgValue () const

Double_t	GetMinimumDist (EdbSegP seg1, EdbSegP seg2)

TObjArray *	GetRecoShowerArray () const

Int_t	GetRecoShowerArrayN () const

Double_t	GetSpatialDist (EdbSegP s1, EdbSegP s2)

virtual void	Initialize ()

Bool_t	IsInConeTube (EdbSegP sTest, EdbSegP sStart, Double_t CylinderRadius, Double_t ConeAngle)

void	Print ()

void	PrintAll ()

void	PrintMore ()

void	PrintParameters ()

void	PrintParametersShort ()

void	PrintRecoShowerArray ()

void	SetActualAlgParameterset (Int_t ActualAlgParametersetNr)

void	SetDebug ()

void	SetEdbPVRec (EdbPVRec *Ali)

void	SetEdbPVRecPIDNumbers (Int_t FirstPlate_eAliPID, Int_t LastPlate_eAliPID, Int_t MiddlePlate_eAliPID, Int_t NumberPlate_eAliPID)

void	SetInBTArray (TObjArray *InBTArray)

void	SetParameter (Int_t parNr, Float_t parvalue)

void	SetParameters (Float_t *par)

void	SetRecoShowerArray (TObjArray *RecoShowerArray)

void	SetRecoShowerArrayN (Int_t RecoShowerArrayN)

void	SetUseAliSub (Bool_t UseAliSub)

void	Transform_eAli (EdbSegP *InitiatorBT, Float_t ExtractSize)

virtual	~EdbShowerAlg ()

Protected Member Functions
void	Set0 ()

Protected Attributes
Int_t	eActualAlgParametersetNr

TString	eAlgName

Int_t	eAlgValue

EdbPVRec *	eAli

EdbPVRec *	eAli_Sub

Int_t	eAli_SubNpat

Int_t	eAliNpat

Bool_t	eDebug

Int_t	eFirstPlate_eAliPID

TObjArray *	eInBTArray

Int_t	eInBTArrayN

Int_t	eLastPlate_eAliPID

Int_t	eMiddlePlate_eAliPID

Int_t	eNumberPlate_eAliPID

TString	eParaString [10]

Float_t	eParaValue [10]

EdbTrackP *	eRecoShower

TObjArray *	eRecoShowerArray

Int_t	eRecoShowerArrayN

Int_t	eUseAliSub

Constructor & Destructor Documentation

◆ EdbShowerAlg() [1/2]

EdbShowerAlg::EdbShowerAlg ( )

◆ EdbShowerAlg() [2/2]

EdbShowerAlg::EdbShowerAlg	(	TString	AlgName,
		Int_t	AlgValue
	)

 {
     // Reset all:
     Set0();
  
     eAlgName=AlgName;
     eAlgValue=AlgValue;
     for (int i=0; i<10; ++i) {
         eParaValue[i]=-99999.0;
         eParaString[i]="UNSPECIFIED";
     }
     eAli_Sub=0;
 }

◆ ~EdbShowerAlg()

EdbShowerAlg::~EdbShowerAlg ( )

virtual

 {
     // Default Destructor
     cout << "EdbShowerAlg::~EdbShowerAlg()"<<endl;
 }

Member Function Documentation

◆ AddRecoShowerArray()

void EdbShowerAlg::AddRecoShowerArray ( TObjArray * RecoShowerArray )

 {
     // Add an array to the existing RecoShower array.
     // For retrieving objects and handling them: be carefull that
     // you add objects of the same class.
  
     Log(3,"EdbShowerAlg::AddRecoShowerArray","AddRecoShowerArray...");
  
     if (RecoShowerArray==NULL || RecoShowerArray->GetEntries()==0) return;
  
     // TObject* obj = RecoShowerArray->At(0);
     // if (obj->ClassName()!="EdbTrackP") { cout << "WARNING AddRecoShowerArray ClassNames dont match!"<< endl; cout << obj->ClassName() << endl; return; } // libShower
     // if (obj->ClassName()!="EdbShowerP") { cout << "WARNING AddRecoShowerArray ClassNames dont match!"<< endl; return; } // libShowRec
  
     for (int i=0; i<RecoShowerArray->GetEntries(); ++i) {
         EdbTrackP* obj = (EdbTrackP*)RecoShowerArray->At(i); // libShower
         //EdbShowerP* obj = (EdbShowerP*)RecoShowerArray->At(0); // libShowRec
         eRecoShowerArray->Add(obj);
         ++eRecoShowerArrayN;
     }
  
     return;
 }

◆ CalcVertex()

EdbVertex * EdbShowerAlg::CalcVertex ( TObjArray * segments )

                                                         {
  
     // Exactly same implementation as in EdbEDAUtil.C
     // but I do not want to link in hardcoded all eda libraries.
     // calc vertex point with given segments. just topological calculation.
     // VTA is currently not set.
     // in case of Single-stop. make a vertex 650 micron upstream ob base.
  
  
     double xx,yy,zz,Det,Ax,Ay,Az,Bx,By,Bz,Cx,Cy,Cz,Dx,Dy,Dz;
     Ax=0.;
     Ay=0.;
     Az=0.;
     Bx=0.;
     By=0.;
     Bz=0.;
     Cx=0.;
     Cy=0.;
     Cz=0.;
     Dx=0.;
     Dy=0.;
     Dz=0.;
  
     if (segments->GetEntries()==1) {
         // in case of Single-stop. make a vertex 650 micron upstream ob base.
         EdbSegP *s = new EdbSegP(*((EdbSegP *) segments->At(0)));
         s->PropagateTo(s->Z()-650);
         xx=s->X();
         yy=s->Y();
         zz=s->Z();
         delete s;
     }
     else {
         for (int i=0; i<segments->GetEntries(); i++) {
             EdbSegP *s = (EdbSegP *) segments->At(i);
             double ax = s->TX();
             double ay = s->TY();
             double az = 1.0;
             double x  = s->X();
             double y  = s->Y();
             double z  = s->Z();
             double a = ax*ax+ay*ay+az*az;
             double c = -ax*x-ay*y-az*z;
             double b = (ax*ax+ay*ay);
             //      double w = 1.0/a/a; // weight for small angle tracks.
             double w = 1.0; // no weight
  
             Ax+=2.0*w/b*( a*(ay*ay+az*az) );
             Bx+=2.0*w/b*( -a*ax*ay );
             Cx+=2.0*w/b*( -a*ax*az );
             Dx+=2.0*w/b*( -(a*x+c*ax)*(ax*ax-a)-(a*y+c*ay)*ax*ay-(a*z+c*az)*az*ax );
  
             Ay+=2.0*w/b*( -a*ay*ax );
             By+=2.0*w/b*( a*(az*az+ax*ax) );
             Cy+=2.0*w/b*( -a*ay*az );
             Dy+=2.0*w/b*( -(a*y+c*ay)*(ay*ay-a)-(a*z+c*az)*ay*az-(a*x+c*ax)*ax*ay );
  
             Az+=2.0*w/b*( -a*az*ax );
             Bz+=2.0*w/b*( -a*az*ay );
             Cz+=2.0*w/b*( a*b );
             Dz+=2.0*w/b*( -(a*z+c*az)*(az*az-a)-(a*x+c*ax)*az*ax-(a*y+c*ay)*ay*az );
  
         }
  
         Det=fabs( Ax*(By*Cz-Cy*Bz)-Bx*(Ay*Cz-Cy*Az)+Cx*(Ay*Bz-By*Az) );
         xx=( (By*Cz-Cy*Bz)*Dx-(Bx*Cz-Cx*Bz)*Dy+(Bx*Cy-Cx*By)*Dz)/Det;
         yy=(-(Ay*Cz-Cy*Az)*Dx+(Ax*Cz-Cx*Az)*Dy-(Ax*Cy-Cx*Ay)*Dz)/Det;
         zz=( (Ay*Bz-By*Az)*Dx-(Ax*Bz-Bx*Az)*Dy+(Ax*By-Bx*Ay)*Dz)/Det;
     }
  
     EdbVertex *v = new EdbVertex();
     v->SetXYZ(xx,yy,zz);
  
     return v;
 }

◆ ClassDef()

EdbShowerAlg::ClassDef	(	EdbShowerAlg	,
		1
	)

◆ Convert_EdbPVRec_To_InBTArray()

void EdbShowerAlg::Convert_EdbPVRec_To_InBTArray ( )

 {
     // This function takes all basetracks of the EdbPVRec object,
     // and fills the InBTArray with those.
     // This function is useful, if one wants to start from all
     // possible basetracks.
  
     if (eInBTArray!=NULL) {
         eInBTArray->Clear();
         cout << "EdbShowerAlg::Convert_EdbPVRec_To_InBTArray() eInBTArray->Clear() done." << endl;
     }
     if (eInBTArray==NULL) {
         // (Should never happen.)
         eInBTArray = new TObjArray();
         cout << " eInBTArray = new TObjArray()" << endl;
     }
  
     cout << "EdbShowerAlg::Convert_EdbPVRec_To_InBTArray() eInBTArray at address "<< eInBTArray <<endl;
  
  
     eInBTArray->Clear();
     Int_t npat=eAli->Npatterns();
     EdbSegP* Segment2;
     EdbPattern* pat_one;
     for (Int_t j=0; j<npat; ++j) {
         pat_one=(EdbPattern*)eAli->GetPattern(j);
         for (Int_t i=0; i<  pat_one->N(); i++) {
             Segment2 = (EdbSegP*)pat_one->GetSegment(i);
             eInBTArray->Add(Segment2);
         }
     }
     eInBTArrayN=eInBTArray->GetEntries();
  
     cout << "Convert_EdbPVRec_To_InBTArray   eInBTArray->GetEntries() ";
     cout << eInBTArray->GetEntries()  << endl;
     cout << "EdbShowerAlg::Convert_EdbPVRec_To_InBTArray()...done."<<endl;
     return;
 }

◆ DeltaR_NoPropagation()

Double_t EdbShowerAlg::DeltaR_NoPropagation	(	EdbSegP *	s,
		EdbSegP *	stest
	)

 {
     // Calculate the distance between two basetracks by
     // sqrt((deltaX)^2+(deltaY)^2)
     // Do NOT propagate tracks onto same Z, in case they are different.
     // SAME function as DeltaR_WithoutPropagation !!!
     return DeltaR_WithoutPropagation(s,stest);
 }

◆ DeltaR_WithoutPropagation()

Double_t EdbShowerAlg::DeltaR_WithoutPropagation	(	EdbSegP *	s,
		EdbSegP *	stest
	)

 {
     // Calculate the distance between two basetracks by
     // sqrt((deltaX)^2+(deltaY)^2)
     // Do NOT propagate tracks onto same Z, in case they are different.
     return TMath::Sqrt((s->X()-stest->X())*(s->X()-stest->X())+(s->Y()-stest->Y())*(s->Y()-stest->Y()));
 }

◆ DeltaR_WithPropagation()

Double_t EdbShowerAlg::DeltaR_WithPropagation	(	EdbSegP *	s,
		EdbSegP *	stest
	)

 {
     // Calculate the distance between two basetracks by
     // sqrt((deltaX)^2+(deltaY)^2)
     // DO propagate tracks onto same Z, in case they are different.
     // (propagation is done from stest onto Z-position of s)
     if (s->Z()==stest->Z()) return TMath::Sqrt((s->X()-stest->X())*(s->X()-stest->X())+(s->Y()-stest->Y())*(s->Y()-stest->Y()));
     Double_t zorig;
     Double_t dR;
     zorig=s->Z();
     s->PropagateTo(stest->Z());
     dR=TMath::Sqrt( (s->X()-stest->X())*(s->X()-stest->X())+(s->Y()-stest->Y())*(s->Y()-stest->Y()) );
     s->PropagateTo(zorig);
     return dR;
 }

◆ DeltaTheta()

Double_t EdbShowerAlg::DeltaTheta	(	EdbSegP *	s1,
		EdbSegP *	s2
	)

 {
     // Be aware that this DeltaTheta function returns the abs() difference between the
     // ABSOLUTE values of dTheta!!! (not componentwise!
     Double_t tx1,tx2,ty1,ty2;
     tx1=s1->TX();
     tx2=s2->TX();
     ty1=s1->TY();
     ty2=s2->TY();
     Double_t dt= TMath::Abs(TMath::Sqrt(tx1*tx1+ty1*ty1) - TMath::Sqrt(tx2*tx2+ty2*ty2));
     return dt;
 }

◆ DeltaThetaComponentwise()

Double_t EdbShowerAlg::DeltaThetaComponentwise	(	EdbSegP *	s1,
		EdbSegP *	s2
	)

 {
     // Be aware that this DeltaTheta function returns the difference between the
     // component values of dTheta!!!
     // Acutally this function should be the normal way to calculate dTheta correctly...
     Double_t tx1,tx2,ty1,ty2;
     tx1=s1->TX();
     tx2=s2->TX();
     ty1=s1->TY();
     ty2=s2->TY();
     Double_t dt= TMath::Sqrt( (tx1-tx2)*(tx1-tx2) + (ty1-ty2)*(ty1-ty2) );
     return dt;
 }

◆ DeltaThetaSingleAngles()

Double_t EdbShowerAlg::DeltaThetaSingleAngles	(	EdbSegP *	s1,
		EdbSegP *	s2
	)

 {
     // SAME function as DeltaThetaComponentwise !!!
     return DeltaThetaComponentwise(s1,s2);
 }

◆ Execute()

void EdbShowerAlg::Execute ( )

virtual

Reimplemented in EdbShowerAlg_GS.

 {
     cout << "EdbShowerAlg::Execute()    -----------------------------------------------------------" << endl;
  
     cout << "EdbShowerAlg::Execute()   This function should be implemented in the inheriting classes!!" << endl;
     return;
 }

◆ Finalize()

void EdbShowerAlg::Finalize ( )

virtual

Reimplemented in EdbShowerAlg_GS.

 {
  
     return;
 }

◆ GetAlgName()

TString EdbShowerAlg::GetAlgName ( ) const

inline

                                             {
         return eAlgName;
     }

◆ GetAlgValue()

Int_t EdbShowerAlg::GetAlgValue ( ) const

inline

                                              {
         return eAlgValue;
     }

◆ GetMinimumDist()

Double_t EdbShowerAlg::GetMinimumDist	(	EdbSegP *	seg1,
		EdbSegP *	seg2
	)

 {
     // Calculate minimum distance of 2 lines.
     // Use the data of (the selected object)->X(), Y(), Z(), TX(), TY().
     // means, if the selected object == segment, use the data of the segment. or it == track, the use the fitted data.
     // Original code from Tomoko Ariga
     // (double calc_dmin(EdbSegP *seg1, EdbSegP *seg2, double *dminz = NULL))
  
     double x1,y1,z1,ax1,ay1;
     double x2,y2,z2,ax2,ay2;
     double s1,s2,s1bunsi,s1bunbo,s2bunsi,s2bunbo;
     double p1x,p1y,p1z,p2x,p2y,p2z,p1p2;
  
     if (seg1->ID()==seg2->ID()&&seg1->PID()==seg2->PID()) return 0.0;
  
     x1 = seg1->X();
     y1 = seg1->Y();
     z1 = seg1->Z();
     ax1= seg1->TX();
     ay1= seg1->TY();
  
     x2 = seg2->X();
     y2 = seg2->Y();
     z2 = seg2->Z();
     ax2= seg2->TX();
     ay2= seg2->TY();
  
     s1bunsi=(ax2*ax2+ay2*ay2+1)*(ax1*(x2-x1)+ay1*(y2-y1)+z2-z1) - (ax1*ax2+ay1*ay2+1)*(ax2*(x2-x1)+ay2*(y2-y1)+z2-z1);
     s1bunbo=(ax1*ax1+ay1*ay1+1)*(ax2*ax2+ay2*ay2+1) - (ax1*ax2+ay1*ay2+1)*(ax1*ax2+ay1*ay2+1);
     s2bunsi=(ax1*ax2+ay1*ay2+1)*(ax1*(x2-x1)+ay1*(y2-y1)+z2-z1) - (ax1*ax1+ay1*ay1+1)*(ax2*(x2-x1)+ay2*(y2-y1)+z2-z1);
     s2bunbo=(ax1*ax1+ay1*ay1+1)*(ax2*ax2+ay2*ay2+1) - (ax1*ax2+ay1*ay2+1)*(ax1*ax2+ay1*ay2+1);
     s1=s1bunsi/s1bunbo;
     s2=s2bunsi/s2bunbo;
     p1x=x1+s1*ax1;
     p1y=y1+s1*ay1;
     p1z=z1+s1*1;
     p2x=x2+s2*ax2;
     p2y=y2+s2*ay2;
     p2z=z2+s2*1;
     p1p2=sqrt( (p1x-p2x)*(p1x-p2x)+(p1y-p2y)*(p1y-p2y)+(p1z-p2z)*(p1z-p2z) );
  
     return p1p2;
 }

◆ GetRecoShowerArray()

TObjArray* EdbShowerAlg::GetRecoShowerArray ( ) const

inline

                                                              {
         return eRecoShowerArray;
     }

◆ GetRecoShowerArrayN()

Int_t EdbShowerAlg::GetRecoShowerArrayN ( ) const

inline

                                                               {
         return eRecoShowerArrayN;
     }

◆ GetSpatialDist()

Double_t EdbShowerAlg::GetSpatialDist	(	EdbSegP *	s1,
		EdbSegP *	s2
	)

 {
     // Mainly Z values should dominate...
     // since the are at the order of 10k microns and x,y of 1k microns
     Double_t x1,x2,y1,y2,z1,z2;
     x1=s1->X();
     x2=s2->X();
     y1=s1->Y();
     y2=s2->Y();
     z1=s1->Z();
     z2=s2->Z();
     Double_t dist= TMath::Sqrt( (x1-x2)*(x1-x2) + (y1-y2)*(y1-y2) + (z1-z2)*(z1-z2)  );
     //cout << "dist = "  <<  dist << endl;
     return dist;
 }

◆ Initialize()

void EdbShowerAlg::Initialize ( )

virtual

Reimplemented in EdbShowerAlg_GS.

 {
  
     return;
 }

◆ IsInConeTube()

Bool_t EdbShowerAlg::IsInConeTube	(	EdbSegP *	sTest,
		EdbSegP *	sStart,
		Double_t	CylinderRadius,
		Double_t	ConeAngle
	)

Outside if angle greater than ConeAngle (to be fulfilled for Cone and Tube in both cases)

if angle smaller than ConeAngle, then you can differ between Tuberadius and CylinderRadius

 {
     // General Function which returns Bool if the Testing BaeTrack is in a cone defined
     // by the StartingBaseTrack. In case of starting same Z position, a distance cut of
     // 20microns is assumed....
     // In case of  TestingSegment==StartingSegment this function should correctly return kTRUE also...
     if (gEDBDEBUGLEVEL>3) cout << "Bool_t EdbShowerAlg::IsInConeTube() Test Segment " << TestingSegment << " vs. Starting Segment " << StartingSegment << endl;
  
     // We reject any TestingSegment segments which have lower Z than the StartingSegment .
     if (StartingSegment->Z()>TestingSegment->Z() ) return kFALSE;
  
     TVector3 x1(StartingSegment->X(),StartingSegment->Y(),StartingSegment->Z());
     TVector3 x2(TestingSegment->X(),TestingSegment->Y(),TestingSegment->Z());
     TVector3 direction_x1(StartingSegment->TX()*1300,StartingSegment->TY()*1300,1300);
  
     // u1 is the difference vector of the position!
     TVector3 u1=x2-x1;
  
     Double_t direction_x1_norm= direction_x1.Mag();
     Double_t cosangle=  (direction_x1*u1)/(u1.Mag()*direction_x1_norm);
     Double_t angle = TMath::ACos(cosangle);
  
     // This is the old version of angle calculation. It does not give the same results as in ROOT
     // when use TVector3.Angle(&TVector3). // For this IsInConeTube() we use therefore the ROOT calculation.
     angle=u1.Angle(direction_x1);
  
     // For the case where the two basetracks have same z position
     // the angle is about 90 degree so it makes no sense to calculate it...
     // therefore we set it artificially to zero:
     if (StartingSegment->Z()==TestingSegment->Z() ) {
         if (gEDBDEBUGLEVEL>3) cout << "same Z position of TestingSegment and StartingSegment, Set angle artificially to zero" << endl;
         angle=0.0;
         // Check here for dR manually:
         //cout << DeltaR_WithoutPropagation(StartingSegment,TestingSegment) << endl;
         //StartingSegment->PrintNice();
         //TestingSegment->PrintNice();
 //   cout << StartingSegment->Flag() << " " << TestingSegment->Flag() << endl;
 //   cout << StartingSegment->P() << " " << TestingSegment->P() << endl;
  
         // Check for position distance for 20microns if
         // Testing Segment is in same Z as StartingSegment
         if (gEDBDEBUGLEVEL>3) cout << "Check for position distance for 20microns if Testing Segment is in same Z as StartingSegment" << endl;
         if (gEDBDEBUGLEVEL>3) cout << "DeltaR_WithoutPropagation(StartingSegment,TestingSegment) = "<< DeltaR_WithoutPropagation(StartingSegment,TestingSegment) << endl;
         if (DeltaR_WithoutPropagation(StartingSegment,TestingSegment)<20) return kTRUE;
         if (DeltaR_WithoutPropagation(StartingSegment,TestingSegment)>=20) return kFALSE;
     }
  
     if (gEDBDEBUGLEVEL>3) cout << "Check if AngleVector now within the ConeAngleVector (<"<< ConeAngle<<"): " <<   angle << endl;
     if (angle>ConeAngle) {
         return kFALSE;
     }
  
     Double_t TubeDistance = 1.0/direction_x1_norm  *  ( (x2-x1).Cross(direction_x1) ).Mag();
  
     if (gEDBDEBUGLEVEL>3) cout << "Check if TestingSegment is now within the Tube (<"<< CylinderRadius<<"): " <<   TubeDistance << endl;
  
     if (TubeDistance>CylinderRadius) {
         return kFALSE;
     }
  
     return kTRUE;
 }

◆ Print()

void EdbShowerAlg::Print ( )

 {
     cout << "EdbShowerAlg::Print()" << endl;
     cout << eAlgName << "  ;  AlgValue=  " <<  eAlgName << "  ." << endl;
     for (int i=0; i<10; i++) cout << eParaString[i] << "=" << eParaValue[i] << ";  ";
     cout << endl;
     cout << "PID numbers: ..." << eFirstPlate_eAliPID << "  " <<  eLastPlate_eAliPID << "  " << eMiddlePlate_eAliPID << "  " << eNumberPlate_eAliPID << "  " << endl;
     cout << "UseAliSub=  " <<  eUseAliSub << "  ." << endl;
     cout << "EdbShowerAlg::Print()...done." << endl;
     return;
 }

◆ PrintAll()

void EdbShowerAlg::PrintAll ( )

 {
     cout << "------------------------------------------------------------" << endl;
     cout << "EdbShowerAlg::PrintAll()" << endl;
     cout << "------------------------" << endl;
     Print();
     PrintMore();
     cout << "------------------------------------------------------------" << endl;
     return;
 }

◆ PrintMore()

void EdbShowerAlg::PrintMore ( )

 {
     cout << "EdbShowerAlg::PrintMore()" << endl;
     cout << "eInBTArray->GetEntries()       =  " <<  eInBTArray->GetEntries() << "  ." << endl;
     cout << "eRecoShowerArray->GetEntries() =  " <<  eRecoShowerArray->GetEntries() << "  ." << endl;
     cout << "EdbShowerAlg::PrintMore()...done." << endl;
     return;
 }

◆ PrintParameters()

void EdbShowerAlg::PrintParameters ( )

 {
     cout << "EdbShowerAlg::PrintParameters()" << endl;
     cout << eAlgName<< " :" << endl;
     for (int i=0; i<5; i++) cout << setw(6) << eParaString[i];
     cout << endl;
     for (int i=0; i<5; i++) cout << setw(6) << eParaValue[i];
     cout <<  " ."<<endl;
     return;
 }

◆ PrintParametersShort()

void EdbShowerAlg::PrintParametersShort ( )

 {
     cout << eAlgName<< " :";
     for (int i=0; i<5; i++) cout << setw(6) << eParaValue[i];
     cout <<  " ."<<endl;
     return;
 }

◆ PrintRecoShowerArray()

void EdbShowerAlg::PrintRecoShowerArray ( )

 {
     cout << "EdbShowerAlg::PrintRecoShowerArray()..." << endl;
     cout << "EdbShowerAlg::PrintRecoShowerArray()   eRecoShowerArray->GetEntries() =  " <<  eRecoShowerArray->GetEntries() << "  ." << endl;
  
     if (eRecoShowerArrayN<1) return;
     printf("i     X   Y   Z    NBT    \n");
  
     Bool_t caseMC=kFALSE;
  
 //      case a) MC event info there:
     EdbTrackP* tr = (EdbTrackP*) eRecoShowerArray->At(0);
     EdbSegP* seg= tr->GetSegment(0);
  
     if (seg->MCEvt()>0) caseMC =kTRUE;
  
     for (Int_t i=0; i<eRecoShowerArray->GetEntries(); ++i) {
         if (caseMC==kTRUE) {
  
             EdbTrackP* tr = (EdbTrackP*) eRecoShowerArray->At(i);
             printf("%06d %06f %06f %06f  %06d  \n", i, tr->X(), tr->Y(), tr->Z(), tr->N());
         }
         else {
             EdbTrackP* tr = (EdbTrackP*) eRecoShowerArray->At(i);
             printf("%06d %06f %06f %06f  %06d  \n", i, tr->X(), tr->Y(), tr->Z(), tr->N());
         }
     }
  
 //      case a) MC event info there:
     cout << "EdbShowerAlg::PrintRecoShowerArray()...done." << endl;
     return;
 }

◆ Set0()

void EdbShowerAlg::Set0 ( )

protected

 {
  
     // Reset internal variable to default values.
  
     cout << "EdbShowerAlg::Set0()"<<endl;
  
     eAlgName="UNSPECIFIED";
     eAlgValue=-999999;
     for (int i=0; i<10; ++i) {
         eParaValue[i]=-99999.0;
         eParaString[i]="UNSPECIFIED";
     }
     eAli_Sub=0;
  
     // do not use use small eAli Object by default:
     // (this solution is memory safe....)
     eUseAliSub=0;
  
     eActualAlgParametersetNr=0;
  
     eRecoShowerArrayN=0;
     eRecoShowerArray=NULL;
     eInBTArrayN=0;
     eInBTArray=NULL;
 }

◆ SetActualAlgParameterset()

void EdbShowerAlg::SetActualAlgParameterset ( Int_t ActualAlgParametersetNr )

inline

                                                                                         {
         eActualAlgParametersetNr=ActualAlgParametersetNr;
     }

◆ SetDebug()

void EdbShowerAlg::SetDebug ( )

inline

                            {
         eDebug=kTRUE;
     }

◆ SetEdbPVRec()

void EdbShowerAlg::SetEdbPVRec ( EdbPVRec * Ali )

inline

                                                               {
         eAli = Ali;
         eAliNpat=eAli->Npatterns();
     }

◆ SetEdbPVRecPIDNumbers()

void EdbShowerAlg::SetEdbPVRecPIDNumbers	(	Int_t	FirstPlate_eAliPID,
		Int_t	LastPlate_eAliPID,
		Int_t	MiddlePlate_eAliPID,
		Int_t	NumberPlate_eAliPID
	)

inline

                                                                                                                                                        {
         eFirstPlate_eAliPID=FirstPlate_eAliPID;
         eLastPlate_eAliPID=LastPlate_eAliPID;
         eMiddlePlate_eAliPID=MiddlePlate_eAliPID;
         eNumberPlate_eAliPID=NumberPlate_eAliPID;
     }

◆ SetInBTArray()

void EdbShowerAlg::SetInBTArray ( TObjArray * InBTArray )

inline

                                                              {
         eInBTArray = InBTArray;
         eInBTArrayN=eInBTArray->GetEntries();
     }

◆ SetParameter()

void EdbShowerAlg::SetParameter	(	Int_t	parNr,
		Float_t	parvalue
	)

 {
     if (parNr>9) {
         cout << "EdbShowerAlg::SetParameter()   WARNING  parNr>9 .  Do nothing."<<endl;
         return;
     }
     eParaValue[parNr]=parvalue;
     cout << "EdbShowerAlg::SetParameter()...done."<<endl;
 }

◆ SetParameters()

void EdbShowerAlg::SetParameters ( Float_t * par )

 {
     cout << "EdbShowerAlg::SetParameters()..."<<endl;
     // SetParameters
     for (int i=0; i<10; ++i) {
         eParaValue[i]=par[i];
         cout << "EdbShowerAlg::SetParameters()...Parameter " << i << " set to: " << eParaValue[i]  <<endl;
     }
     cout << "EdbShowerAlg::SetParameters()...done."<<endl;
 }

◆ SetRecoShowerArray()

void EdbShowerAlg::SetRecoShowerArray ( TObjArray * RecoShowerArray )

inline

                                                                        {
         eRecoShowerArray=RecoShowerArray;
         eRecoShowerArrayN=eRecoShowerArray->GetEntries();
     }

◆ SetRecoShowerArrayN()

void EdbShowerAlg::SetRecoShowerArrayN ( Int_t RecoShowerArrayN )

inline

                                                                     {
         eRecoShowerArrayN=RecoShowerArrayN;
     }

◆ SetUseAliSub()

void EdbShowerAlg::SetUseAliSub ( Bool_t UseAliSub )

inline

                                                        {
         eUseAliSub=UseAliSub;
     }

◆ Transform_eAli()

void EdbShowerAlg::Transform_eAli	(	EdbSegP *	InitiatorBT,
		Float_t	ExtractSize = `1500`
	)

DEBUG if (eAli_Sub) { delete eAli_Sub;eAli_Sub=0;} // original, but keeps not adresses of segment in eAli.

do nothing now... let it live... delete eAli_Sub;eAli_Sub=0;}

if (eAli_Sub) { delete eAli_Sub;eAli_Sub=0;} Try not to delete it maybe then it works.....

 {
     // --------------------------------------------------------------------------------------
     //      Transform eAli to eAli_Sub:
     //      the lenght of eAli_Sub is not changed.
     //      Only XY-size (and MC) cuts are applied.
     // --------------------------------------------------------------------------------------
     // ---
     // --- Whereas in ShowRec.cpp the treebranch file is written directly after each
     // --- BT reconstruction, it was not a problem when the eAliSub was deleted each time.
     // --- But now situation is different, since the BTs of the showers have the adresses from
     // --- the eAli_Sub !each! so if one eAli_Sub is deleted the EdbShowerP object has lost its
     // --- BT adresses.
     // --- So if we do not delete the eAli_Sub, reconstrtuction is fast (compared to eAli) but
     // --- memory increases VERY fast.
     // --- If we do use eAli, then memory consumption will not increase fast, but
     // --- reconstruction is slow.
     // ---
     // --- Possible workarounds:
     // ---  For few InBTs (EDA use for data...):        use eAliSub and dont delete it.
     // ---  For many InBTs (parameterstudies):          use eAli---.
     // ---  Use always eAliSub, search BT correspond in eAli, add BT from eAli....May take also long time...
     // ---
     // ---
     // --- SEVERE warning:  IF gAli is in wrong order it can be that no showers
     // --- SEVERE warning:  are reconstructed since most implemented algorithms
     // --- SEVERE warning:  rely on the order that plate 2 comes directly behind
     // --- SEVERE warning:  the InBT.  THIS IS STILL ON BUGFIXING LIST!!!
     // --------------------------------------------------------------------------------------
  
  
     /*
     cout << "void EdbShowerAlg::Transform_eAli()   SEVERE WARNING:  IF gAli is in wrong order it can be that no showers " << endl;
     cout << "void EdbShowerAlg::Transform_eAli()   SEVERE WARNING:  are reconstructed since most implemented algorithms " << endl;
     cout << "void EdbShowerAlg::Transform_eAli()   SEVERE WARNING:  rely on the order that plate 2 comes directly behind " << endl;
     cout << "void EdbShowerAlg::Transform_eAli()   SEVERE WARNING:  the InBT.   /// DEBUG  TODO  // DEBUG TODO " << endl;
     */
  
  
 //  gEDBDEBUGLEVEL=3;
 //  cout << "void EdbShowerAlg::Transform_eAli()     ===== eUseAliSub "  << eUseAliSub << endl;
  
     // IF TO RECREATE THE gALI_SUB in RECONSTRUCTION  or if to use gAli global (slowlier but maybe memory leak safe).
     if (!eUseAliSub)
     {
         if (gEDBDEBUGLEVEL>2) cout << "EdbShowerAlg::Transform_eAli   UseAliSub==kFALSE No new eAli_Sub created. Use eAli instead. "<<endl;
         eAli_Sub=eAli;
         eAli_SubNpat=eAli_Sub->Npatterns();  //number of plates
         if (gEDBDEBUGLEVEL>3) eAli_Sub->Print();
         return;
     }
  
     Int_t npat;
     npat = eAli->Npatterns();  //number of plates
  
     // has to be deleted in some part of the script outside this function...
     // Dont forget , otherwise memory heap overflow!
     if (eAli_Sub) {
         ;
     } 
 //  eAli_Sub = new EdbPVRec();
  
  
     if (eUseAliSub) {
         if (gEDBDEBUGLEVEL>2) cout << "EdbShowerAlg::Transform_eAli   UseAliSub==kTRUE Will now create new eAli_Sub. "<<endl;
  
         eAli_Sub = new EdbPVRec();
 //      cout << "Adress of eAli_Sub " << eAli_Sub << endl;
 //      cout << "eAli_Sub->Npatterns();  " << eAli_Sub->Npatterns() << endl;
 //      eAli_Sub->Print();
     }
  
  
  
 //  cout << "PROBLE THRANFOMING, CAUSE ADRESES ARE DELETED ALSO...." << endl;
 //      cout << "TEMPORAY SOLUTION:     comment the delete  eAli_Sub  ( will lead to large memory consumption when run for long time" << endl;
  
  
     // Create SubPattern objects
     EdbSegP* ExtrapolateInitiatorBT=0;
     ExtrapolateInitiatorBT = (EdbSegP*)InitiatorBT->Clone();
  
     Int_t InitiatorBTMCEvt=InitiatorBT->MCEvt();
  
     // Create Variables For ExtractSubpattern boundaries
     Float_t mini[5];
     Float_t maxi[5];
     //Float_t ExtractSize=1000;// now in fucntion header
     mini[0]=ExtrapolateInitiatorBT->X()-ExtractSize;
     mini[1]=ExtrapolateInitiatorBT->Y()-ExtractSize;
     maxi[0]=ExtrapolateInitiatorBT->X()+ExtractSize;
     maxi[1]=ExtrapolateInitiatorBT->Y()+ExtractSize;
     mini[2]=-0.7;
     mini[3]=-0.7;
     mini[4]=0.0;
     maxi[2]=0.7;
     maxi[3]=0.7;
     maxi[4]=100.0;
  
     EdbPattern* singlePattern;
     Float_t ExtrapolateInitiatorBT_zpos_orig=ExtrapolateInitiatorBT->Z();
  
     // Add the subpatterns in a loop for the plates:
     // in reverse ordering.due to donwstream behaviour (!):
     // (Only downstream is supported now...)
     for (Int_t ii=0; ii<npat; ++ii) {
  
         Float_t zpos=eAli->GetPattern(ii)->Z();
  
         ExtrapolateInitiatorBT->PropagateTo(zpos);
  
         mini[0]=ExtrapolateInitiatorBT->X()-ExtractSize;
         mini[1]=ExtrapolateInitiatorBT->Y()-ExtractSize;
         maxi[0]=ExtrapolateInitiatorBT->X()+ExtractSize;
         maxi[1]=ExtrapolateInitiatorBT->Y()+ExtractSize;
  
         singlePattern=(EdbPattern*)eAli->GetPattern(ii)->ExtractSubPattern(mini,maxi,InitiatorBTMCEvt);
         // This sets PID() analogue to (upstream), nut not PID of the BTs !
         singlePattern-> SetID(eAli->GetPattern(ii)->ID());
         // This sets PID() analogue to (upstream), nut not PID of the BTs !
         singlePattern-> SetPID(eAli->GetPattern(ii)->PID());
  
         eAli_Sub->AddPattern(singlePattern);
  
         // Propagate back...! (in order not to change the original BT)
         ExtrapolateInitiatorBT->PropagateTo(ExtrapolateInitiatorBT_zpos_orig);
     }
  
     delete ExtrapolateInitiatorBT;
  
     eAli_SubNpat=eAli_Sub->Npatterns();  //number of plates
  
  
     if (gEDBDEBUGLEVEL>2) {
         cout << "EdbShowerAlg::Transform_eAli   eAli_Sub  created."<<endl;
         cout << "Adress of eAli_Sub " << eAli_Sub << endl;
         cout << "eAli_Sub->Npatterns();  " << eAli_Sub->Npatterns() << endl;
     }
     return;
 }

Member Data Documentation

◆ eActualAlgParametersetNr

Int_t EdbShowerAlg::eActualAlgParametersetNr

protected

◆ eAlgName

TString EdbShowerAlg::eAlgName

protected

◆ eAlgValue

Int_t EdbShowerAlg::eAlgValue

protected

◆ eAli

EdbPVRec* EdbShowerAlg::eAli

protected

◆ eAli_Sub

EdbPVRec* EdbShowerAlg::eAli_Sub

protected

◆ eAli_SubNpat

Int_t EdbShowerAlg::eAli_SubNpat

protected

◆ eAliNpat

Int_t EdbShowerAlg::eAliNpat

protected

◆ eDebug

Bool_t EdbShowerAlg::eDebug

protected

◆ eFirstPlate_eAliPID

Int_t EdbShowerAlg::eFirstPlate_eAliPID

protected

◆ eInBTArray

TObjArray* EdbShowerAlg::eInBTArray

protected

◆ eInBTArrayN

Int_t EdbShowerAlg::eInBTArrayN

protected

◆ eLastPlate_eAliPID

Int_t EdbShowerAlg::eLastPlate_eAliPID

protected

◆ eMiddlePlate_eAliPID

Int_t EdbShowerAlg::eMiddlePlate_eAliPID

protected

◆ eNumberPlate_eAliPID

Int_t EdbShowerAlg::eNumberPlate_eAliPID

protected

◆ eParaString

TString EdbShowerAlg::eParaString[10]

protected

◆ eParaValue

Float_t EdbShowerAlg::eParaValue[10]

protected

◆ eRecoShower

EdbTrackP* EdbShowerAlg::eRecoShower

protected

◆ eRecoShowerArray

TObjArray* EdbShowerAlg::eRecoShowerArray

protected

◆ eRecoShowerArrayN

Int_t EdbShowerAlg::eRecoShowerArrayN

protected

◆ eUseAliSub

Int_t EdbShowerAlg::eUseAliSub

protected

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

src/libShower/EdbShowerAlg.h
src/libShower/EdbShowerAlg.cxx

Public Member Functions

Protected Member Functions

Protected Attributes

Constructor & Destructor Documentation

◆ EdbShowerAlg() [1/2]

◆ EdbShowerAlg() [2/2]

◆ ~EdbShowerAlg()

Member Function Documentation

◆ AddRecoShowerArray()

◆ CalcVertex()

◆ ClassDef()

◆ Convert_EdbPVRec_To_InBTArray()

◆ DeltaR_NoPropagation()

◆ DeltaR_WithoutPropagation()

◆ DeltaR_WithPropagation()

◆ DeltaTheta()

◆ DeltaThetaComponentwise()

◆ DeltaThetaSingleAngles()

◆ Execute()

◆ Finalize()

◆ GetAlgName()

◆ GetAlgValue()

◆ GetMinimumDist()

◆ GetRecoShowerArray()

◆ GetRecoShowerArrayN()

◆ GetSpatialDist()

◆ Initialize()

◆ IsInConeTube()

◆ Print()

◆ PrintAll()

◆ PrintMore()

◆ PrintParameters()

◆ PrintParametersShort()

◆ PrintRecoShowerArray()

◆ Set0()

◆ SetActualAlgParameterset()

◆ SetDebug()

◆ SetEdbPVRec()

◆ SetEdbPVRecPIDNumbers()

◆ SetInBTArray()

◆ SetParameter()

◆ SetParameters()

◆ SetRecoShowerArray()

◆ SetRecoShowerArrayN()

◆ SetUseAliSub()

◆ Transform_eAli()

Member Data Documentation

◆ eActualAlgParametersetNr

◆ eAlgName

◆ eAlgValue

◆ eAli

◆ eAli_Sub

◆ eAli_SubNpat

◆ eAliNpat

◆ eDebug

◆ eFirstPlate_eAliPID

◆ eInBTArray

◆ eInBTArrayN

◆ eLastPlate_eAliPID

◆ eMiddlePlate_eAliPID

◆ eNumberPlate_eAliPID

◆ eParaString

◆ eParaValue

◆ eRecoShower

◆ eRecoShowerArray

◆ eRecoShowerArrayN

◆ eUseAliSub