EdbShowAlg_SA Class Reference

#include <EdbShowAlg.h>

Inheritance diagram for EdbShowAlg_SA:

Collaboration diagram for EdbShowAlg_SA:

Public Member Functions
	ClassDef (EdbShowAlg_SA, 1)
	EdbShowAlg_SA ()
void	Execute ()
void	Finalize ()
void	Init ()
void	Initialize ()
virtual	~EdbShowAlg_SA ()
Public Member Functions inherited from EdbShowAlg
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)
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 ()
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 ()

Additional Inherited Members
Protected Member Functions inherited from EdbShowAlg
void	Set0 ()
Protected Attributes inherited from EdbShowAlg
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_SA()

EdbShowAlg_SA::EdbShowAlg_SA

(

)

{
// Default Constructor
    cout << "EdbShowAlg_SA::EdbShowAlg_SA()   Default Constructor   SA stands for ?????? To be filled..."<<endl;
 
// Reset all:
    Set0();
 
    eAlgName="SA";
    eAlgValue=5; // see default.par_SHOWREC for labeling (labling identical with ShowRec program)
    eParaN =2; // depends on the algorithm, how much parameters it has got.
 
// Set the two Cuts and their names:
    Init();
}

~EdbShowAlg_SA()

EdbShowAlg_SA::~EdbShowAlg_SA ( )

virtual

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

Member Function Documentation

ClassDef()

EdbShowAlg_SA::ClassDef	(	EdbShowAlg_SA	,
		1
	)

Execute()

void EdbShowAlg_SA::Execute ( )

virtual

Reimplemented from EdbShowAlg.

{
    cout << "EdbShowAlg_SA::Execute()" << endl;
    cout << "EdbShowAlg_SA::Execute   DOING MAIN SHOWER RECONSTRUCTION HERE" << endl;
 
    EdbSegP* InBT;
    EdbSegP* Segment;
    EdbShowerP* RecoShower;
 
    Bool_t     StillToLoop=kTRUE;
    Int_t       ActualPID;
    Int_t      newActualPID;
    Int_t      STEP=-1;
    Int_t     NLoopedPattern=0;
    if (eFirstPlate_eAliPID-eLastPlate_eAliPID<0) STEP=1;
    if (gEDBDEBUGLEVEL>3) cout << "EdbShowAlg_SA::Execute--- STEP for patternloop direction =  " << STEP << endl;
 
//--- Loop over InBTs:
 
// Since eInBTArray is filled in ascending ordering by zpositon
// We use the descending loop to begin with BT with lowest z first.
    for (Int_t i=eInBTArrayN-1; i>=0; --i) {
 
// CounterOutPut
        if (gEDBDEBUGLEVEL==2) if ((i%100)==0) cout << eInBTArrayN <<" InBT in total, still to do:"<<Form("%4d",i)<< "\r\r\r\r"<<flush;
 
//-----------------------------------
// 1) Make local_gAli with cut parameters:
//-----------------------------------
 
// Create new EdbShowerP Object for storage;
// See EdbShowerP why I have to call the Constructor as "unique" ideficable value
        RecoShower = new EdbShowerP(i,eAlgValue,eActualAlgParametersetNr);
 
// Get InitiatorBT from eInBTArray
        InBT=(EdbSegP*)eInBTArray->At(i);
        if (gEDBDEBUGLEVEL>2) InBT->PrintNice();
 
// Add InBT to RecoShower:
// This has to be done, since by definition the first BT in the RecoShower is the InBT.
// Otherwise, also the definition of shower axis and transversal profiles is wrong!
        RecoShower -> AddSegment(InBT);
        if (gEDBDEBUGLEVEL>4)  cout << "Segment  (InBT) " << Segment << " was added to RecoShower." <<  endl;
 
// Transform (make size smaller, extract only events having same MC) the  eAli  object:
// See in Execute_CA for description.
        Transform_eAli(InBT,eParaValue[1]); // here with ExtractSize set!
 
// Add InBT to RecoShower // obsolete, since we loop also over the plate containing the InBT
//RecoShower --> AddSegment(InBT);
 
//-----------------------------------
// 2) Loop over (whole) eAli, check BT for Cuts
//-----------------------------------
        ActualPID= InBT->PID() ;
        newActualPID= InBT->PID() ;
 
        while (StillToLoop) {
            if (gEDBDEBUGLEVEL>3) cout << "EdbShowAlg_SA::Execute--- --- Doing patterloop " << ActualPID << endl;
 
            for (Int_t btloop_cnt=0; btloop_cnt<eAli->GetPattern(ActualPID)->GetN(); ++btloop_cnt) {
 
                Segment = (EdbSegP*)eAli->GetPattern(ActualPID)->GetSegment(btloop_cnt);
                if (gEDBDEBUGLEVEL>4) Segment->PrintNice();
 
// Now apply cut conditions: SA Simple Alg  --------------------
//   ___  TO FILL IN HERE ___
                if (Segment->MCEvt()!=InBT->MCEvt()) continue;
                if (Segment->P()<eParaValue[0]) continue;
// end of    cut conditions: SA Simple Alg  --------------------
 
// If we arrive here, Basetrack  Segment  has passed criteria
// and is then added to the RecoShower:
// Check if its not the InBT which is already added:
                if (Segment->X()==InBT->X()&&Segment->Y()==InBT->Y()) {
                    ;    // is InBT, do nothing;
                }
                else {
                    RecoShower -> AddSegment(Segment);
                }
                if (gEDBDEBUGLEVEL>4)  cout << "Segment  " << Segment << " was added at  &Segment : " << &Segment  <<  endl;
 
            } // of btloop_cnt
 
//------------
            newActualPID=ActualPID+STEP;
            ++NLoopedPattern;
 
            if (gEDBDEBUGLEVEL>3) cout << "EdbShowAlg_SA::Execute--- --- newActualPID= " << newActualPID << endl;
            if (gEDBDEBUGLEVEL>3) cout << "EdbShowAlg_SA::Execute--- --- NLoopedPattern= " << NLoopedPattern << endl;
            if (gEDBDEBUGLEVEL>3) cout << "EdbShowAlg_SA::Execute--- --- eNumberPlate_eAliPID= " << eNumberPlate_eAliPID << endl;
            if (gEDBDEBUGLEVEL>3) cout << "EdbShowAlg_SA::Execute--- --- StillToLoop= " << StillToLoop << endl;
 
// This if holds in the case of STEP== +1
            if (STEP==1) {
                if (newActualPID>eLastPlate_eAliPID) StillToLoop=kFALSE;
                if (newActualPID>eLastPlate_eAliPID) cout << "EdbShowAlg_SA::Execute--- ---Stop Loop since: newActualPID>eLastPlate_eAliPID"<<endl;
            }
// This if holds in the case of STEP== -1
            if (STEP==-1) {
                if (newActualPID<eLastPlate_eAliPID) StillToLoop=kFALSE;
                if (newActualPID<eLastPlate_eAliPID && gEDBDEBUGLEVEL>3) cout << "EdbShowAlg_SA::Execute--- ---Stop Loop since: newActualPID<eLastPlate_eAliPID"<<endl;
            }
// This if holds  general, since eNumberPlate_eAliPID is not dependent of the structure of the gAli subject:
            if (NLoopedPattern>eNumberPlate_eAliPID) StillToLoop=kFALSE;
            if (NLoopedPattern>eNumberPlate_eAliPID && gEDBDEBUGLEVEL>3) cout << "EdbShowAlg_SA::Execute--- ---Stop Loop since: NLoopedPattern>eNumberPlate_eAliPID"<<endl;
 
            ActualPID=newActualPID;
        } // of // while (StillToLoop)
 
// Obligatory when Shower Reconstruction is finished!
        RecoShower ->Update();
//RecoShower ->PrintBasics();
 
 
//     if (gEDBDEBUGLEVEL>3) cout << "EdbShowAlg_SA::Execute--- Before adding to array delete the histograms...by finalize() of shower."<<endl;
//     RecoShower ->Finalize();
//     if (gEDBDEBUGLEVEL>3) cout << "EdbShowAlg_SA::Execute--- Finalize done..."<<endl;
 
 
// Add Shower Object to Shower Reco Array.
// Not, if its empty:
// Not, if its containing only one BT:
    if (RecoShower->N()>1) {
      AddRecoShowerArray(RecoShower);
    }
    else {
      // Delete the created shower object. Saves _lots_ of memory and speeds up reconstruction for large #InBTs!
      // TODO Also Include this in the other RecoAlgs ...
      Log(3,"EdbShowAlg_SA::Execute()","InBT # %d RecoShower->N()<=1: Dont add shower to RecoShowerArray. Delete this shower for memory safing.");
      delete RecoShower;
      RecoShower=NULL;
    }
 
// Set back loop values:
        StillToLoop=kTRUE;
        NLoopedPattern=0;
    } // of  //   for (Int_t i=eInBTArrayN-1; i>=0; --i) {
 
 
// Set new value for  eRecoShowerArrayN  (may now be < eInBTArrayN).
    SetRecoShowerArrayN(eRecoShowerArray->GetEntries());
 
    cout << "EdbShowAlg_SA::eRecoShowerArray() Entries: " << eRecoShowerArray->GetEntries() << endl;
    cout << "EdbShowAlg_SA::Execute()...done." << endl;
    return;
}

Finalize()

void EdbShowAlg_SA::Finalize ( )

virtual

Reimplemented from EdbShowAlg.

{
    cout << "EdbShowAlg_SA::Finalize()" << endl;
    return;
}

Init()

void EdbShowAlg_SA::Init

(

void

)

{
    Log(2,"EdbShowAlg_SA::EdbShowAlg_SA","Init()");
 
//  Init with values according to ?? Alg:
 
// threshold cut on s->P()
    eParaValue[0]=0;
    eParaString[0]="MINIMUM_P";
 
// size cut on the eAliSub to choose
    eParaValue[1]=1000;
    eParaString[1]="EXTRACTIONSIZE";
 
    return;
}

Initialize()

void EdbShowAlg_SA::Initialize ( )

virtual

Reimplemented from EdbShowAlg.

{
    cout << "EdbShowAlg_SA::Initialize()" << endl;
    return;
}

---

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

• src/libShowRec/EdbShowAlg.h
• src/libShowRec/EdbShowAlg.cxx