ecalMatch Class Reference

#include <ecalMatch.h>

Inheritance diagram for ecalMatch:

Collaboration diagram for ecalMatch:

Public Member Functions
	ecalMatch (const char *name, const Int_t verbose)
	ecalMatch ()
Int_t	GetUse3x3 () const
void	SetUse3x3 (Int_t use3x3=0)
virtual void	Finish ()
virtual void	Exec (Option_t *option, TClonesArray *reconstructed, TClonesArray *mctracks)
virtual InitStatus	Init ()
void	InitPython (ecalStructure *str, TClonesArray *reconstructed, TClonesArray *mctracks)
	~ecalMatch ()

Private Attributes
Int_t	fEv
Int_t	fN
Int_t	fRejected
Int_t	fUse3x3
TClonesArray *	fReconstucted
TClonesArray *	fMCTracks
ecalStructure *	fStr

Detailed Description

Definition at line 11 of file ecalMatch.h.

Constructor & Destructor Documentation

ecalMatch() [1/2]

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

Standard constructor

Definition at line 120 of file ecalMatch.cxx.

  : FairTask(name, verbose), fEv(0), fN(0), fRejected(0), fUse3x3(0), 
    fReconstucted(NULL), fMCTracks(NULL), fStr(NULL) 
{
  ;
}

ecalMatch() [2/2]

ecalMatch::ecalMatch

(

)

Only to comply with frame work.

Definition at line 128 of file ecalMatch.cxx.

  : FairTask(), fEv(0), fN(0), fRejected(0), fUse3x3(0), 
    fReconstucted(NULL), fMCTracks(NULL), fStr(NULL)
{
  ;
}

~ecalMatch()

ecalMatch::~ecalMatch

(

)

Destructor

Definition at line 142 of file ecalMatch.cxx.

{
  ;
}

Member Function Documentation

Exec()

void ecalMatch::Exec ( Option_t *  option, TClonesArray *  reconstructed, TClonesArray *  mctracks  )

virtual

Exec a task

Definition at line 17 of file ecalMatch.cxx.

{
  fReconstucted=reconstructed;
  fMCTracks=mctracks;
  fEv++; fN=0; fRejected=0;
 
  Int_t n=fReconstucted->GetEntries();
  Int_t i;
  ecalReconstructed* rc;
  ecalCell* cell;
  ecalCellMC* mccell;
  list<ecalCell*> cells;
  list<ecalCell*>::const_iterator p;
  map<Int_t, Float_t> e;
  map<Int_t, Float_t> e2;
  map<Int_t, Float_t>::const_iterator ep;
  map<Int_t, Float_t>::reverse_iterator rp;
  ShipMCTrack* tr;
  Int_t trn;
  Float_t max;
//  if (fVerbose>0) Info("Exec", "Event %d.", fEv);
  for(i=0;i<n;i++)
  {
    rc=(ecalReconstructed*)fReconstucted->At(i);
    cell=fStr->GetHitCell(rc->CellNum());
    cells.clear();
 
    if (fUse3x3)
    {
      cell->GetNeighborsList(cells); //3x3 without maximum
      cells.push_back(cell);
    }
    else
      cell->Get5x5Cluster(cells);
 
    e.clear(); e2.clear();
    //Counting energy depositions for all particles
    for(p=cells.begin();p!=cells.end();++p)
    {
      mccell=(ecalCellMC*)(*p);
      for(ep=mccell->GetTrackEnergyBegin();ep!=mccell->GetTrackEnergyEnd();++ep)
      {
    if (e.find(ep->first)==e.end())
      e[ep->first]=ep->second;
    else
      e[ep->first]+=ep->second;
      }
    }
 
    //...and parent photons and electrons/positrons
    for(ep=e.begin();ep!=e.end();++ep)
    {
      if (e2.find(ep->first)==e.end())
    e2[ep->first]=ep->second;
      else
    e2[ep->first]+=ep->second;
      if (ep->first<0&&fVerbose==0) continue;
      tr=(ShipMCTrack*)fMCTracks->At(ep->first);
      if (tr==NULL)
      {
    Info("Exec", "Event %d. Can't find MCTrack %d.", fEv, ep->first);
    continue;
      }
      for(;;)
      {
    if (tr->GetPdgCode()!=22&&TMath::Abs(tr->GetPdgCode())!=11) break;
    trn=tr->GetMotherId();
    if (trn<0) break;
    tr=(ShipMCTrack*)fMCTracks->At(trn);
    if (tr==NULL)
    {
      Info("Exec", "Event %d. Can't find MCTrack %d.", fEv, ep->first);
      break;
    }
    if (tr->GetPdgCode()!=22&&TMath::Abs(tr->GetPdgCode())!=11) break;
    if (e2.find(trn)==e2.end())
      e2[trn]=ep->second;
    else
      e2[trn]+=ep->second;
      }
    }
 
    //Maximum location
    max=-1e11; trn=-1111;
    for(rp=e2.rbegin();rp!=e2.rend();++rp)
    {
      if (rp->second>max)
    { max=rp->second; trn=rp->first;}
    }
 
    if (trn>=0)
    {
      rc->SetMCTrack(trn);
      fN++;
    }
    else
      fRejected++;
  }
 
  if (fVerbose>0) Info("Exec", "Event %d. Matched %d. Skipped %d maxs.", fEv, fN, fRejected);
}

Finish()

void ecalMatch::Finish ( )

virtual

Finish a task

Definition at line 136 of file ecalMatch.cxx.

{
  ;
}

GetUse3x3()

Int_t ecalMatch::GetUse3x3 ( ) const

inline

Should we use 3x3 cluster instead of 5x5

Definition at line 20 of file ecalMatch.h.

{return fUse3x3;}

Init()

InitStatus ecalMatch::Init ( )

virtual

Initialization

Definition at line 147 of file ecalMatch.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;
  }
  fReconstucted=(TClonesArray*)io->GetObject("EcalReco");
  if (!fReconstucted)
  {
    Fatal("Init", "Can't find array of reconstructed calorimeter objects in the system.");
    return kFATAL;
  }
  fMCTracks=(TClonesArray*)io->GetObject("MCTrack");
  if (!fMCTracks)
  {
    Fatal("Init", "Can't find array of MC tracks in the system.");
    return kFATAL;
  }
 
  fEv=0;
  return kSUCCESS;
}

InitPython()

void ecalMatch::InitPython	(	ecalStructure *	str,
		TClonesArray *	reconstructed,
		TClonesArray *	mctracks
	)

Definition at line 178 of file ecalMatch.cxx.

{
  fStr=str;
  fReconstucted=reconstructed;
  fMCTracks=mctracks;
  
  fEv=0;
}

SetUse3x3()

void ecalMatch::SetUse3x3 ( Int_t  use3x3 = 0 )

inline

Definition at line 21 of file ecalMatch.h.

{fUse3x3=use3x3;}

Member Data Documentation

fEv

Int_t ecalMatch::fEv

private

Current event

Definition at line 34 of file ecalMatch.h.

fMCTracks

TClonesArray* ecalMatch::fMCTracks

private

Array of MC tracks

Definition at line 45 of file ecalMatch.h.

fN

Int_t ecalMatch::fN

private

Current reconstructed particle

Definition at line 36 of file ecalMatch.h.

fReconstucted

TClonesArray* ecalMatch::fReconstucted

private

Array of reconstructed objects

Definition at line 43 of file ecalMatch.h.

fRejected

Int_t ecalMatch::fRejected

private

Number of rejected clusters

Definition at line 38 of file ecalMatch.h.

fStr

ecalStructure* ecalMatch::fStr

private

A calorimeter structure

Definition at line 47 of file ecalMatch.h.

fUse3x3

Int_t ecalMatch::fUse3x3

private

Use 3x3 cluster instead of 5x5

Definition at line 40 of file ecalMatch.h.

---

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

• ecal/ecalMatch.h
• ecal/ecalMatch.cxx