#include <ecalClusterFinder.h>

Inheritance diagram for ecalClusterFinder:

Collaboration diagram for ecalClusterFinder:

[legend]

Public Member Functions
	ecalClusterFinder (const char *name, const Int_t verbose)

	ecalClusterFinder ()

virtual	~ecalClusterFinder ()

virtual void	Finish ()

virtual void	Exec (Option_t *option)

virtual InitStatus	Init ()

TClonesArray *	InitPython (ecalStructure structure, TClonesArray maximums, ecalClusterCalibration *calib)

Double_t	MinMaxE () const

Double_t	MinClusterE () const

void	SetMinMaxE (Double_t minmaxe=0.015)

void	SetMinClusterE (Double_t minmaxe=0.03)

Private Member Functions
void	FormClusters ()

void	FormPreClusters ()

void	ClearPreClusters ()

	ecalClusterFinder (const ecalClusterFinder &)

ecalClusterFinder &	operator= (const ecalClusterFinder &)

Private Attributes
Int_t	fEv

TClonesArray *	fMaximums

TClonesArray *	fClusters

ecalStructure *	fStr

ecalClusterCalibration *	fCalib

ecalInf *	fInf

std::list< ecalPreCluster * >	fPreClusters

Double_t	fMinClusterE

Double_t	fMinMaxE

Detailed Description

Definition at line 18 of file ecalClusterFinder.h.

Constructor & Destructor Documentation

◆ ecalClusterFinder() [1/3]

ecalClusterFinder::ecalClusterFinder	(	const char *	name,
		const Int_t	verbose
	)

Standard constructor

Definition at line 222 of file ecalClusterFinder.cxx.

  : FairTask(name, verbose),
    fEv(0),
    fMaximums(NULL),
    fClusters(NULL),
    fStr(NULL),
    fInf(NULL),
    fPreClusters(),
    fMinClusterE(0.03),
    fMinMaxE(0.015)   
{
  ;
}

◆ ecalClusterFinder() [2/3]

ecalClusterFinder::ecalClusterFinder ( )

Only to comply with frame work.

Definition at line 238 of file ecalClusterFinder.cxx.

  : FairTask(),
    fEv(0),
    fMaximums(NULL),
    fClusters(NULL),
    fStr(NULL),
    fInf(NULL),
    fPreClusters(),
    fMinClusterE(0.03),
    fMinMaxE(0.015)   
{
  ;
}

◆ ~ecalClusterFinder()

ecalClusterFinder::~ecalClusterFinder ( )

virtual

Destructor

Definition at line 81 of file ecalClusterFinder.cxx.

{
  if (fClusters)
  {
    fClusters->Delete();
    delete fClusters;
  }
}

◆ ecalClusterFinder() [3/3]

ecalClusterFinder::ecalClusterFinder ( const ecalClusterFinder & )

private

Member Function Documentation

◆ ClearPreClusters()

void ecalClusterFinder::ClearPreClusters ( )

private

Clear a preclusters list

Definition at line 213 of file ecalClusterFinder.cxx.

{
  list<ecalPreCluster*>::const_iterator p=fPreClusters.begin();
  for(;p!=fPreClusters.end();++p)
    delete (*p);
  fPreClusters.clear();
}

◆ Exec()

void ecalClusterFinder::Exec ( Option_t * option )

virtual

Exec a task

Definition at line 22 of file ecalClusterFinder.cxx.

{
  fEv++;
 
  ClearPreClusters();
  FormPreClusters();
  FormClusters();
}

◆ Finish()

void ecalClusterFinder::Finish ( )

virtual

Finish a task

Definition at line 75 of file ecalClusterFinder.cxx.

{
  ;
}

◆ FormClusters()

void ecalClusterFinder::FormClusters ( )

private

Form clusters from precluster

ecalCluster needs a destructor call :-(

Definition at line 144 of file ecalClusterFinder.cxx.

{
  fClusters->Delete();
  Int_t fN=0;
  list<ecalPreCluster*>::const_iterator p1=fPreClusters.begin();
  list<ecalPreCluster*>::const_iterator p2;
  list<ecalCell*> cluster;
  list<ecalMaximum*> maxs;
  list<ecalCell*>::const_iterator pc;
  list<ecalCell*>::const_iterator pc1;
  UInt_t oldsize;
  Int_t MaxSize=0;
  Int_t Maximums=0;
  Int_t max;
  Int_t type;
  
  if (fVerbose>9)
  {
    Info("FormClusters", "Total %d preclusters found.", (Int_t)fPreClusters.size());
  }
  for(;p1!=fPreClusters.end();++p1)
  if ((*p1)->fMark==0)
  {
    cluster.clear(); oldsize=0; maxs.clear();
    cluster=(*p1)->fCells; maxs.push_back((*p1)->fMax); type=(*p1)->fMaximum->GetType();
    max=1;
    while(cluster.size()!=oldsize)
    {
      oldsize=cluster.size();
      p2=p1;
      for(++p2;p2!=fPreClusters.end();++p2)
      if ((*p2)->fMark==0)
      {
        pc=cluster.begin();
    for(;pc!=cluster.end();++pc)
    {
      pc1=find((*p2)->fCells.begin(), (*p2)->fCells.end(), (*pc));
      if (pc1==(*p2)->fCells.end()) continue;
      break;
    }
    if (pc!=cluster.end())
    {
      (*p2)->fMark=1;
      pc=(*p2)->fCells.begin();
      for(;pc!=(*p2)->fCells.end();++pc)
        if (find(cluster.begin(), cluster.end(), (*pc))==cluster.end())
          cluster.push_back(*pc);
      maxs.push_back((*p2)->fMax);
      max++;
    }
      }
    }
    (*p1)->fMark=1;
    if ((Int_t)cluster.size()>MaxSize)
      MaxSize=cluster.size();
    if (max>Maximums) Maximums=max;
    ecalCluster* cls=new ((*fClusters)[fN]) ecalCluster(fN, cluster, maxs); fN++;
    cls->fPreCalibrated=fCalib->Calibrate(type, cls->fEnergy);
  }
  if (fVerbose>0)
  {
    Info("FormClusters", "Total %d clusters formed.", fN);
    Info("FormClusters", "Maximum size of cluster is %d cells.",  MaxSize);
    Info("FormClusters", "Maximum number of photons per cluster is %d.",  Maximums);
  }
}

◆ FormPreClusters()

void ecalClusterFinder::FormPreClusters ( )

private

Form a preclusters

Form a preclusters. A precluster — a group of cells neighbor to maximum cell. A cluster is a group of preclusters with common cells.

Remove maximums matched with charged tracks

Remove low energy maximums

Remove low energy clusters

Definition at line 93 of file ecalClusterFinder.cxx.

{
  Int_t nm=fMaximums->GetEntriesFast();
  Int_t i=0;
  ecalMaximum* max;
  list<ecalCell*> all;
  list<ecalCell*>::const_iterator p;
  list<ecalCell*>::const_iterator p2;
  list<ecalCell*> cls;
  list<ecalCell*> cls2;
  ecalCell* cell;
  ecalCell* min;
  Double_t e;
  Double_t ecls;
  ecalPreCluster* precluster;
 
 
  for(;i<nm;i++)
  {
    max=(ecalMaximum*)fMaximums->At(i);
    if (max==NULL) continue;
    if (max->Mark()!=0) continue;
    cell=max->Cell();
    ecls=cell->GetEnergy();
//    cout << ecls << endl;
    if (ecls<fMinMaxE) continue;
/*
    cell->GetNeighborsList(all);
    cls.clear();
    for(p=all.begin();p!=all.end();++p)
    {
      (*p)->GetNeighborsList(cls2);
      for(p2=cls2.begin();p2!=cls2.end();++p2)
        if (find(cls.begin(), cls.end(), *p2)==cls.end()) cls.push_back(*p2);
    }
*/
    cell->Get5x5Cluster(cls);
    ecls=0.0;
    for(p=cls.begin();p!=cls.end();++p)
      ecls+=(*p)->GetEnergy();
//    cout << ":" << ecls << endl;
    if (ecls<fMinClusterE) continue;
    precluster=new ecalPreCluster(cls, max);
    fPreClusters.push_back(precluster);
  }
}

◆ Init()

InitStatus ecalClusterFinder::Init ( )

virtual

Initialization

Definition at line 31 of file ecalClusterFinder.cxx.

{
  FairRootManager* io=FairRootManager::Instance();
  if (!io)
  {
    Fatal("Init", "Can't find IOManager.");
    return kFATAL;
  }
  fStr=(ecalStructure*)io->GetObject("EcalStructure");
  if (!fStr) 
  {
    Fatal("Init()", "Can't find calorimeter structure in the system.");
    return kFATAL;
  }
  fInf=fStr->GetEcalInf();
  fMaximums=(TClonesArray*)io->GetObject("EcalMaximums");
  if (!fMaximums)
  {
    Fatal("Init", "Can't find array of calorimeter maximums in the system.");
    return kFATAL;
  }
  fCalib=(ecalClusterCalibration*)io->GetObject("ecalClusterCalibration");
  if (!fCalib)
  {
    Fatal("Init", "Can't find ecalClusterCalibration in the system.");
    return kFATAL;
  }
 
  fClusters=new TClonesArray("ecalCluster", 2000);
  io->Register("EcalClusters", "ECAL", fClusters, kTRUE);
  fEv=0;
  return kSUCCESS;
}

◆ InitPython()

TClonesArray * ecalClusterFinder::InitPython	(	ecalStructure *	structure,
		TClonesArray *	maximums,
		ecalClusterCalibration *	calib
	)

Definition at line 65 of file ecalClusterFinder.cxx.

{
  fMaximums=maximums;
  fStr=structure;
  fCalib=calib;
  fClusters=new TClonesArray("ecalCluster", 2000);
  return fClusters;
}

◆ MinClusterE()

Double_t ecalClusterFinder::MinClusterE ( ) const

inline

Definition at line 40 of file ecalClusterFinder.h.

40{return fMinClusterE;}

◆ MinMaxE()

Double_t ecalClusterFinder::MinMaxE ( ) const

inline

Definition at line 39 of file ecalClusterFinder.h.

39{return fMinMaxE;}

◆ operator=()

ecalClusterFinder & ecalClusterFinder::operator= ( const ecalClusterFinder & )

private

◆ SetMinClusterE()

void ecalClusterFinder::SetMinClusterE ( Double_t minmaxe = 0.03 )

inline

Minimum uncalibrated energy of precluster maximum for consideration

Definition at line 45 of file ecalClusterFinder.h.

45{fMinClusterE=minmaxe;}

◆ SetMinMaxE()

void ecalClusterFinder::SetMinMaxE ( Double_t minmaxe = 0.015 )

inline

Minimum precluster uncalibrated energy

Definition at line 43 of file ecalClusterFinder.h.

43{fMinMaxE=minmaxe;}

Member Data Documentation

◆ fCalib

ecalClusterCalibration* ecalClusterFinder::fCalib

private

Cluster calibration object for photons

Definition at line 64 of file ecalClusterFinder.h.

◆ fClusters

TClonesArray* ecalClusterFinder::fClusters

private

An array of clusters

Definition at line 60 of file ecalClusterFinder.h.

◆ fEv

Int_t ecalClusterFinder::fEv

private

Current event

Definition at line 54 of file ecalClusterFinder.h.

◆ fInf

ecalInf* ecalClusterFinder::fInf

private

An information about calorimeter

Definition at line 66 of file ecalClusterFinder.h.

◆ fMaximums

TClonesArray* ecalClusterFinder::fMaximums

private

Array of maximums in calorimeter. Maximums belong to charged tracks excluded.

Definition at line 58 of file ecalClusterFinder.h.

◆ fMinClusterE

Double_t ecalClusterFinder::fMinClusterE

private

Minimum precluster uncalibrated energy

Definition at line 73 of file ecalClusterFinder.h.

◆ fMinMaxE

Double_t ecalClusterFinder::fMinMaxE

private

Minimum uncalibrated energy of precluster maximum for consideration

Definition at line 75 of file ecalClusterFinder.h.

◆ fPreClusters

std::list<ecalPreCluster*> ecalClusterFinder::fPreClusters

private

A list of preclusters May be better use TClonesArray?

Definition at line 70 of file ecalClusterFinder.h.

◆ fStr

ecalStructure* ecalClusterFinder::fStr

private

A calorimeter structure

Definition at line 62 of file ecalClusterFinder.h.

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

ecal/ecalClusterFinder.h
ecal/ecalClusterFinder.cxx

Public Member Functions

Private Member Functions

Private Attributes

Detailed Description

Constructor & Destructor Documentation

◆ ecalClusterFinder() [1/3]

◆ ecalClusterFinder() [2/3]

◆ ~ecalClusterFinder()

◆ ecalClusterFinder() [3/3]

Member Function Documentation

◆ ClearPreClusters()

◆ Exec()

◆ Finish()

◆ FormClusters()

◆ FormPreClusters()

◆ Init()

◆ InitPython()

◆ MinClusterE()

◆ MinMaxE()

◆ operator=()

◆ SetMinClusterE()

◆ SetMinMaxE()

Member Data Documentation

◆ fCalib

◆ fClusters

◆ fEv

◆ fInf

◆ fMaximums

◆ fMinClusterE

◆ fMinMaxE

◆ fPreClusters

◆ fStr