EdbShowAlg Class Reference

#include <EdbShowAlg.h>

Inheritance diagram for EdbShowAlg:

Collaboration diagram for EdbShowAlg:

Public Member Functions
void	AddRecoShowerArray (EdbShowerP *shower)
	ClassDef (EdbShowAlg, 1)
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)
	EdbShowAlg ()
	EdbShowAlg (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
EdbShowerP *	GetShower (Int_t i) const
Double_t	GetSpatialDist (EdbSegP *s1, EdbSegP *s2)
void	Help ()
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	SetActualAlgParameterset (Int_t ActualAlgParametersetNr)
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 par)
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)
void	UpdateShowerIDs ()
void	UpdateShowerMetaData ()
virtual	~EdbShowAlg ()

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
Int_t	eFirstPlate_eAliPID
TObjArray *	eInBTArray
Int_t	eInBTArrayN
Int_t	eLastPlate_eAliPID
Int_t	eMiddlePlate_eAliPID
Int_t	eNumberPlate_eAliPID
Int_t	eParaN
TString	eParaString [10]
Float_t	eParaValue [10]
EdbShowerP *	eRecoShower
TObjArray *	eRecoShowerArray
Int_t	eRecoShowerArrayN
Int_t	eUseAliSub

Constructor & Destructor Documentation

EdbShowAlg() [1/2]

EdbShowAlg::EdbShowAlg

(

)

EdbShowAlg() [2/2]

EdbShowAlg::EdbShowAlg	(	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;
}

~EdbShowAlg()

EdbShowAlg::~EdbShowAlg ( )

virtual

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

Member Function Documentation

AddRecoShowerArray()

void EdbShowAlg::AddRecoShowerArray ( EdbShowerP *  shower )

inline

                                                               {
        eRecoShowerArray->Add(shower);
        ++eRecoShowerArrayN;
    }

ClassDef()

EdbShowAlg::ClassDef	(	EdbShowAlg	,
		1
	)

DeltaR_NoPropagation()

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

{
// SAME function as DeltaR_WithoutPropagation !!!
    return DeltaR_WithoutPropagation(s,stest);
}

DeltaR_WithoutPropagation()

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

{
    return TMath::Sqrt((s->X()-stest->X())*(s->X()-stest->X())+(s->Y()-stest->Y())*(s->Y()-stest->Y()));
}

DeltaR_WithPropagation()

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

{
    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 EdbShowAlg::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 EdbShowAlg::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 EdbShowAlg::DeltaThetaSingleAngles	(	EdbSegP *	s1,
		EdbSegP *	s2
	)

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

Execute()

void EdbShowAlg::Execute ( )

virtual

Reimplemented in EdbShowAlg_SA, EdbShowAlg_CA, EdbShowAlg_OI, EdbShowAlg_RC, EdbShowAlg_BW, EdbShowAlg_GS, EdbShowAlg_NN, and EdbShowAlg_N3.

{
    cout << "EdbShowAlg::Execute()-----------------------------------------------------------------------" << endl;
    return;
}

Finalize()

void EdbShowAlg::Finalize ( )

virtual

Reimplemented in EdbShowAlg_SA, EdbShowAlg_CA, EdbShowAlg_OI, EdbShowAlg_RC, EdbShowAlg_BW, EdbShowAlg_GS, EdbShowAlg_NN, and EdbShowAlg_N3.

{
    return;
}

GetAlgName()

TString EdbShowAlg::GetAlgName ( ) const

inline

                                            {
        return eAlgName;
    }

GetAlgValue()

Int_t EdbShowAlg::GetAlgValue ( ) const

inline

                                             {
        return eAlgValue;
    }

GetMinimumDist()

Double_t EdbShowAlg::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 * EdbShowAlg::GetRecoShowerArray ( ) const

inline

                                                             {
        return eRecoShowerArray;
    }

GetRecoShowerArrayN()

Int_t EdbShowAlg::GetRecoShowerArrayN ( ) const

inline

                                                              {
        return eRecoShowerArrayN;
    }

GetShower()

EdbShowerP * EdbShowAlg::GetShower ( Int_t  i ) const

inline

                                                             {
        EdbShowerP* shower = (EdbShowerP*)eRecoShowerArray->At(i);
        return shower;
    }

GetSpatialDist()

Double_t EdbShowAlg::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;
}

Help()

void EdbShowAlg::Help

(

)

{
    cout << "---------------------------------------------------------" << endl;
    cout << "EdbShowAlg::Help()" << endl << endl;
    cout << "  What is this?" << endl;
    cout << "        Primary Reconstruction Class for the Showers" << endl;
    cout << "        Commonly, it takes as basic data input an EdbPVRec object (gAli)" << endl;
    cout << endl;
    cout << "  What can we do with this?" << endl;
    cout << "        First of all, you can choose different types of algos:" << endl;
    cout << "        ShowAlg_OI (default)" << endl;
    cout << "        ShowAlg_AS" << endl;
    cout << "        ShowAlg_RC" << endl;
    cout << "        ShowAlg_TC" << endl;
    cout << "        ShowAlg_NN" << endl;
    cout << "        ShowAlg_RC (*** to be implemented ***)" << endl;
    cout << "        ShowAlg_BW (*** to be implemented ***)" << endl;
    cout << "        ShowAlg_CL (*** to be skipped ***)" << endl;
    cout << "        ShowAlg_GS (*** to be implemented ***)" << endl;
    cout << "        ...." << endl;
    cout << "        Then you set you input via: " << endl;
    cout << "        ->SetEdbPVRec(gAli) : " << endl;
    cout << "        Then you execute the algorithm via: " << endl;
    cout << "        ->Execute() : " << endl;
    cout << "        Finally you get the TObjArray of reconstructed showers (EdbShowerP*) via: " << endl;
    cout << "        ->GetRecoShowerArray() : " << endl;
    cout << endl;
    cout << "   Technical details, information about the class:" << endl;
    cout << "        ...." << endl;
    cout << "        ...." << endl;
    cout << endl;
    cout << "---------------------------------------------------------" << endl;
    return;
}

Initialize()

void EdbShowAlg::Initialize ( )

virtual

Reimplemented in EdbShowAlg_SA, EdbShowAlg_CA, EdbShowAlg_OI, EdbShowAlg_RC, EdbShowAlg_BW, EdbShowAlg_GS, EdbShowAlg_NN, and EdbShowAlg_N3.

{
    return;
}

IsInConeTube()

Bool_t EdbShowAlg::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 EdbShowAlg::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 EdbShowAlg::Print

(

)

{
    Log(2,"EdbShowAlg::Print()","Print relevant settings for shower reconstruction algorithm");
 
    Log(2,"EdbShowAlg::Print()","Algorithm name abbreviation: %s", eAlgName.Data() );
    Log(2,"EdbShowAlg::Print()","Algorithm number (internal): %2d", eAlgValue );
 
    Log(2,"EdbShowAlg::Print()","Algorithm parameters N:      %2d", eParaN );
    for (int i=0; i<eParaN; i++) Log(2,"EdbShowAlg::Print()","Algorithm parameter: number, name, value:  %d, %40s, %.04f", i, eParaString[i].Data(), eParaValue[i] );
 
//     cout << "Plates: " << endl;
//     cout << eFirstPlate_eAliPID << "  " <<  eLastPlate_eAliPID << "  " << eMiddlePlate_eAliPID << "  " << eNumberPlate_eAliPID << "  " << endl;
//     cout << "UseAliSub=  " <<  eUseAliSub << "  ." << endl;
//     cout << "EdbShowAlg::Print()...done." << endl;
    Log(2,"EdbShowAlg::Print()","Print relevant settings for shower reconstruction algorithm...done.");
    return;
}

PrintAll()

void EdbShowAlg::PrintAll

(

)

{
    return;
}

PrintMore()

void EdbShowAlg::PrintMore

(

)

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

PrintParameters()

void EdbShowAlg::PrintParameters

(

)

{
    Print();
    return;
}

PrintParametersShort()

void EdbShowAlg::PrintParametersShort

(

)

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

Set0()

void EdbShowAlg::Set0 ( )

protected

{
// Reset all internal variables.
    eParaN =10; // depends on the algorithm, how much parameters it has got. to be set in the constructor of the specific class...
    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 EdbShowAlg::SetActualAlgParameterset ( Int_t  ActualAlgParametersetNr )

inline

                                                                            {
        eActualAlgParametersetNr=ActualAlgParametersetNr;
    }

SetEdbPVRec()

void EdbShowAlg::SetEdbPVRec ( EdbPVRec *  Ali )

inline

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

SetEdbPVRecPIDNumbers()

void EdbShowAlg::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 EdbShowAlg::SetInBTArray ( TObjArray *  InBTArray )

inline

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

SetParameter()

void EdbShowAlg::SetParameter	(	Int_t	parNr,
		Float_t	par
	)

{
// SetParameter from the given value.
    if (parNr>9) {
        cout << "EdbShowAlg::SetParameters() parNr>9 not supported. Dont set par value!"<<endl;
    }
    eParaValue[parNr]=par;
    cout << "EdbShowAlg::SetParameter()...done"<<endl;
}

SetParameters()

void EdbShowAlg::SetParameters

(

Float_t *

par

)

{
// SetParameters from the given array.
    for (int i=0; i<10; ++i) {
        eParaValue[i]=par[i];
    }
    cout << "EdbShowAlg::SetParameters()...done"<<endl;
}

SetRecoShowerArray()

void EdbShowAlg::SetRecoShowerArray ( TObjArray *  RecoShowerArray )

inline

                                                                       {
        eRecoShowerArray=RecoShowerArray;
    }

SetRecoShowerArrayN()

void EdbShowAlg::SetRecoShowerArrayN ( Int_t  RecoShowerArrayN )

inline

                                                                    {
        eRecoShowerArrayN=RecoShowerArrayN;
    }

SetUseAliSub()

void EdbShowAlg::SetUseAliSub ( Bool_t  UseAliSub )

inline

                                                       {
        eUseAliSub=UseAliSub;
    }

Transform_eAli()

void EdbShowAlg::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
// ---                   the InBT. on the todo list....
// ---
// --------------------------------------------------------------------------------------
 
 
    Log(3,"EdbShowAlg::Transform_eAli(EdbSegP* InitiatorBT, Float_t ExtractSize=1500)","EdbShowAlg::Transform_eAli");
 
// 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 << "EdbShowAlg::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 << "EdbShowAlg::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 << "EdbShowAlg::Transform_eAli   eAli_Sub  created at : "<< eAli_Sub << endl;
        cout << "EdbShowAlg::Transform_eAli   eAli_Sub->Npatterns(): "<< eAli_Sub->Npatterns() << endl;
    }
 
    Log(3,"EdbShowAlg::Transform_eAli(EdbSegP* InitiatorBT, Float_t ExtractSize=1500)...done.","EdbShowAlg::Transform_eAli...done.");
 
    return;
}

UpdateShowerIDs()

void EdbShowAlg::UpdateShowerIDs

(

)

{
    // Shower IDs are numbered in the ordering in the RecoShowerArray
    if ( NULL == eRecoShowerArray ) return;
    Int_t idNumber=0;
    for (int i=0; i<eRecoShowerArrayN; ++i) {
        EdbShowerP* shower = GetShower(i);
        shower->SetID(i);
    }
    return;
}

UpdateShowerMetaData()

void EdbShowAlg::UpdateShowerMetaData

(

)

{
    // To know which Algorithm and Parameters were used 
    // for this shower to be reconstructed.
    if ( NULL == eRecoShowerArray ) return;
    Int_t idNumber=0;
    for (int i=0; i<eRecoShowerArrayN; ++i) {
        EdbShowerP* shower = GetShower(i);
        shower->SetAlgName(eAlgName);
    shower->SetAlgValue(eAlgValue);
    for (Int_t j=0; j<10; ++j) {
      shower->SetAlgParaValue(j, eParaValue[j]);
      shower->SetAlgParaString(j, eParaString[j]);
    }
    }
    return;
}

Member Data Documentation

eActualAlgParametersetNr

Int_t EdbShowAlg::eActualAlgParametersetNr

protected

eAlgName

TString EdbShowAlg::eAlgName

protected

eAlgValue

Int_t EdbShowAlg::eAlgValue

protected

eAli

EdbPVRec* EdbShowAlg::eAli

protected

eAli_Sub

EdbPVRec* EdbShowAlg::eAli_Sub

protected

eAli_SubNpat

Int_t EdbShowAlg::eAli_SubNpat

protected

eAliNpat

Int_t EdbShowAlg::eAliNpat

protected

eFirstPlate_eAliPID

Int_t EdbShowAlg::eFirstPlate_eAliPID

protected

eInBTArray

TObjArray* EdbShowAlg::eInBTArray

protected

eInBTArrayN

Int_t EdbShowAlg::eInBTArrayN

protected

eLastPlate_eAliPID

Int_t EdbShowAlg::eLastPlate_eAliPID

protected

eMiddlePlate_eAliPID

Int_t EdbShowAlg::eMiddlePlate_eAliPID

protected

eNumberPlate_eAliPID

Int_t EdbShowAlg::eNumberPlate_eAliPID

protected

eParaN

Int_t EdbShowAlg::eParaN

protected

eParaString

TString EdbShowAlg::eParaString[10]

protected

eParaValue

Float_t EdbShowAlg::eParaValue[10]

protected

eRecoShower

EdbShowerP* EdbShowAlg::eRecoShower

protected

eRecoShowerArray

TObjArray* EdbShowAlg::eRecoShowerArray

protected

eRecoShowerArrayN

Int_t EdbShowAlg::eRecoShowerArrayN

protected

eUseAliSub

Int_t EdbShowAlg::eUseAliSub

protected

---

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

• /home/antonio/fedra_doxygen/src/libShowRec/EdbShowAlg.h
• /home/antonio/fedra_doxygen/src/libShowRec/EdbShowAlg.cxx