EdbShowAlg_RC Class Reference

#include <EdbShowAlg.h>

Inheritance diagram for EdbShowAlg_RC:

Collaboration diagram for EdbShowAlg_RC:

Public Member Functions
	ClassDef (EdbShowAlg_RC, 1)
	EdbShowAlg_RC ()
void	Execute ()
void	Finalize ()
Bool_t	GetConeOrTubeDistanceToBTOfShowerArray (EdbSegP *sa, EdbSegP *InBT, EdbShowerP *shower, Double_t CylinderRadius, Double_t ConeAngle)
void	Initialize ()
virtual	~EdbShowAlg_RC ()
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)
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 ()

Private Member Functions
Bool_t	FindPrecedingBTs (EdbSegP *s, EdbSegP *InBT, EdbPVRec *gAli, EdbShowerP *shower)

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_RC()

EdbShowAlg_RC::EdbShowAlg_RC

(

)

{
// Default Constructor
    cout << "EdbShowAlg_RC::EdbShowAlg_RC()   Default Constructor  Recursive Cone."<<endl;
 
// Reset all:
    Set0();
 
    eAlgName="RC";
    eAlgValue=7; // see default.par_SHOWREC for labeling (labeling identical with ShowRec program)
    eParaN =4; // depends on the algorithm, how much parameters it has got.
 
    eParaValue[0]=700.0;
    eParaString[0]="CylinderRadius";  // micrometer
    eParaValue[1]=0.025;
    eParaString[1]="ConeAngle";       // radiant
    eParaValue[2]=150.0;
    eParaString[2]="DeltaRToPreceedingBT";  // micrometer
    eParaValue[3]=0.13;
    eParaString[3]="DeltaThetaComponentwiseToPreceedingBT";  // radiant
}

~EdbShowAlg_RC()

EdbShowAlg_RC::~EdbShowAlg_RC ( )

virtual

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

Member Function Documentation

ClassDef()

EdbShowAlg_RC::ClassDef	(	EdbShowAlg_RC	,
		1
	)

Execute()

void EdbShowAlg_RC::Execute ( )

virtual

Reimplemented from EdbShowAlg.

{
    cout << "EdbShowAlg_RC::Execute()" << endl;
    cout << "EdbShowAlg_RC::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:
    if (gEDBDEBUGLEVEL>2) cout << "EdbShowAlg_RC::Execute    Loop over InBTs:" << endl;
 
 
 
// 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();
        Float_t X0=InBT->X();
        Float_t Y0=InBT->Y();
 
// 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:
        Transform_eAli(InBT);
        if (gEDBDEBUGLEVEL>3) eAli_Sub->Print();
 
//-----------------------------------
// 2) Loop over (whole) eAli, check BT for Cuts
// eAli_Sub
//-----------------------------------
        ActualPID= InBT->PID() ;
        newActualPID= InBT->PID() ;
 
        while (StillToLoop) {
            if (gEDBDEBUGLEVEL>3) cout << "EdbShowAlg_RC::Execute--- --- Doing patterloop " << ActualPID << " for patterns Z position=" << eAli_Sub->GetPattern(ActualPID)->Z() << endl;
 
            for (Int_t btloop_cnt=0; btloop_cnt<eAli_Sub->GetPattern(ActualPID)->GetN(); ++btloop_cnt) {
 
                Segment = (EdbSegP*)eAli_Sub->GetPattern(ActualPID)->GetSegment(btloop_cnt);
                if (gEDBDEBUGLEVEL>4) Segment->PrintNice();
 
 
// Now apply cut conditions: RC RECURSIVE CONE Alg  --------------------
                if (Segment->MCEvt() > 0 ) if ( Segment->MCEvt()!=InBT->MCEvt() ) continue; // MCEvtNr (>0) or BgMCNr (-999)
                if ( Abs(Segment->X()-X0) > 7000 ) continue;
                if ( Abs(Segment->Y()-Y0) > 7000 ) continue;
                if (!GetConeOrTubeDistanceToBTOfShowerArray(Segment, InBT, RecoShower, eParaValue[0], eParaValue[1])) continue;
                if (!FindPrecedingBTs(Segment, InBT, eAli_Sub, RecoShower)) continue;
// end of    cut conditions: RC RECURSIVE CONE 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
 
            if (gEDBDEBUGLEVEL>3) cout << "EdbShowAlg_RC::Execute--- --- ActualPID= " << newActualPID << "  done. Reconstructed shower has up to now: " << RecoShower->N()  << " Segments." << endl;
 
//------------
            newActualPID=ActualPID+STEP;
            ++NLoopedPattern;
 
            if (gEDBDEBUGLEVEL>3) cout << "EdbShowAlg_RC::Execute--- --- newActualPID= " << newActualPID << endl;
            if (gEDBDEBUGLEVEL>3) cout << "EdbShowAlg_RC::Execute--- --- NLoopedPattern= " << NLoopedPattern << endl;
            if (gEDBDEBUGLEVEL>3) cout << "EdbShowAlg_RC::Execute--- --- eNumberPlate_eAliPID= " << eNumberPlate_eAliPID << endl;
            if (gEDBDEBUGLEVEL>3) cout << "EdbShowAlg_RC::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_RC::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_RC::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_RC::Execute--- ---Stop Loop since: NLoopedPattern>eNumberPlate_eAliPID"<<endl;
 
            ActualPID=newActualPID;
        } // of // while (StillToLoop)
 
// Obligatory when Shower Reconstruction is finished:
        RecoShower ->Update();
 
        if (gEDBDEBUGLEVEL>3) RecoShower ->PrintBasics();
        if (gEDBDEBUGLEVEL>3) RecoShower ->PrintNice();
        if (gEDBDEBUGLEVEL>3) RecoShower ->PrintSegments();
 
//     if (gEDBDEBUGLEVEL>3) cout << "EdbShowAlg_RC::Execute--- Before adding to array delete the histograms...by finalize() of shower."<<endl;
//     RecoShower ->Finalize();
//     if (gEDBDEBUGLEVEL>3) cout << "EdbShowAlg_RC::Execute--- Finalize done..."<<endl;
 
 
// Add Shower Object to Shower Reco Array.
// Not, if its empty or containing only one BT:
        if (RecoShower->N()>1) eRecoShowerArray->Add(RecoShower);
    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_RC::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_RC::eRecoShowerArray() Entries: " << eRecoShowerArray->GetEntries() << endl;
    cout << "EdbShowAlg_RC::Execute()...done." << endl;
    return;
}

Finalize()

void EdbShowAlg_RC::Finalize ( )

virtual

Reimplemented from EdbShowAlg.

{
    return;
}

FindPrecedingBTs()

Bool_t EdbShowAlg_RC::FindPrecedingBTs ( EdbSegP *  s, EdbSegP *  InBT, EdbPVRec *  gAli, EdbShowerP *  shower  )

private

{
    if (gEDBDEBUGLEVEL>3) cout << "EdbShowAlg_RC::FindPrecedingBTs   Find BTs to be connected with the Test BT:" << endl;
 
    EdbSegP* s_TestBT;
    Int_t nentries=shower->N();
    Double_t dZ;
 
// In case the shower is empty we do not have to search for a preceeding BT:
    if (nentries==0) return kTRUE;
 
// We do not test BaseTracks which are before the initiator BaseTrack!
    if (s->Z()<InBT->Z()) {
        if (gEDBDEBUGLEVEL>3) cout << "--- --- EdbShowAlg_RC::FindPrecedingBTs(): s->Z()<InBT->Z()..: We do not test BaseTracks which are before the initiator BaseTrack!  return kFALSE;" << endl;
        return kFALSE;
    }
 
// For the very first Z position we do not test
// if testBT has Preceeders, only if it it has a BT around (case for e+e- coming from gammma):
// Take 50microns and 80mrad in (dR/dT) around.
// This does not affect the normal results, but helps for
// events which may have a second BT close to InBT (like in e+e-)
    if (s->Z()==InBT->Z()) {
//cout << DeltaThetaComponentwise(s, InBT) << "  " << DeltaR_WithPropagation(s, InBT) << endl;
        if (DeltaThetaComponentwise(s, InBT) < 0.08 && DeltaR_WithoutPropagation(s, InBT) < 50.0 ) return kTRUE;
    }
 
    if (gEDBDEBUGLEVEL>3) cout << "--- --- EdbShowAlg_RC::FindPrecedingBTs(): Testing " << nentries << " entries of the shower." << endl;
 
    for (Int_t i=nentries-1; i>=0; --i) {
        s_TestBT = (EdbSegP*)( shower->GetSegment(i) );
 
        if (gEDBDEBUGLEVEL>3) cout << "--- --- EdbShowAlg_RC::FindPrecedingBTs(): Do  s_TestBT->ID()  s->ID()  s_TestBT->MCEvt() s_TestBT->Z()  s->Z() "<< s_TestBT->ID() << " " << s->ID() << " " << s_TestBT->MCEvt() <<"  " << s_TestBT->Z() << "  " <<  s->Z()  <<endl;
        if (gEDBDEBUGLEVEL>3) s_TestBT->PrintNice();
        if (gEDBDEBUGLEVEL>3) s->PrintNice();
 
 
        dZ=TMath::Abs(s_TestBT->Z()-s->Z());
        if (dZ<30) continue;                  // Exclude the case of same Zpositions...
        if (dZ>(3*1300.0)+30.0) continue;     // Exclude the case of more than 3 plates before...
 
        if (gEDBDEBUGLEVEL>3) cout << "--- --- EdbShowAlg_RC::FindPrecedingBTs(): Checking dT,dR and dZ for i:  " << i << "  " << DeltaThetaComponentwise(s, s_TestBT)  << "  " << DeltaR_WithPropagation(s, s_TestBT) << "  "<<dZ << endl;
 
        if (DeltaThetaComponentwise(s, s_TestBT) >  eParaValue[3] ) continue;
        if (DeltaR_WithPropagation(s, s_TestBT) > eParaValue[2]) continue;
 
        if (gEDBDEBUGLEVEL>3) cout << "--- --- EdbShowAlg_RC::FindPrecedingBTs(): Checking dT,dR and dZ for i:  " << i << "  " << DeltaThetaComponentwise(s, s_TestBT)  << "  " << DeltaR_WithPropagation(s, s_TestBT) << "  "<<dZ << "   ok!"<<endl;
        return kTRUE;
    }
//---------------------------------------------
    return kFALSE;
}

GetConeOrTubeDistanceToBTOfShowerArray()

Bool_t EdbShowAlg_RC::GetConeOrTubeDistanceToBTOfShowerArray	(	EdbSegP *	sa,
		EdbSegP *	InBT,
		EdbShowerP *	shower,
		Double_t	CylinderRadius,
		Double_t	ConeAngle
	)

{
    Bool_t isTrueForBT=kFALSE;
    Int_t lastI=0;
    Double_t factor=1.0;
 
    EdbSegP* s_TestBT;
    Int_t nentries=shower->GetNBT();
 
// Now call IsInConeTube for every BT which was reconstructed up to now in the shower:
    for (Int_t i=0; i<nentries; i++) {
        s_TestBT = (EdbSegP*)( shower->GetSegment(i) );
// Dont check the BT which is before or same Z as the BaseTrack(i) position:
// But ---not--- for the first Basetrack. There we allow it!
//if (s_TestBT->Z()>=sa->Z() && i>0) continue;
        if (s_TestBT->Z()>sa->Z() ) continue;
        lastI=i;
        if (i==0) {
            factor=1.0;
        }
        else {
            factor=3.0;
        }
 
//     Bool_t EdbShowAlg::IsInConeTube(EdbSegP* TestingSegment, EdbSegP* StartingSegment, Double_t CylinderRadius, Double_t ConeAngle)
//     if (GetConeOrTubeDistanceToInBT(sa, s_TestBT, factor*CylinderRadius, factor*ConeAngle)==kTRUE) {
        if (IsInConeTube(sa, s_TestBT, factor*CylinderRadius, factor*ConeAngle)==kTRUE) {
            isTrueForBT=kTRUE;
            break;
        }
    }
 
    if (isTrueForBT) {
//if (lastI>0) cout <<"for i= " << lastI<< " return true with factor "<< factor << endl;
        return kTRUE;
    }
 
    return kFALSE;
}

Initialize()

void EdbShowAlg_RC::Initialize ( )

virtual

Reimplemented from EdbShowAlg.

{
    return;
}

---

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