PndEmcClusterRemoveDuplCrys Class Reference

Task to remove duplicate crystals. More...

#include <PndEmcClusterRemoveDuplCrys.h>

Inheritance diagram for PndEmcClusterRemoveDuplCrys:

Public Member Functions
	PndEmcClusterRemoveDuplCrys (Int_t verbose=0)
virtual	~PndEmcClusterRemoveDuplCrys ()
virtual InitStatus	Init ()
virtual void	Exec (Option_t *opt)

Private Member Functions
	PndEmcClusterRemoveDuplCrys (const PndEmcClusterRemoveDuplCrys &)
PndEmcClusterRemoveDuplCrys &	operator= (const PndEmcClusterRemoveDuplCrys &)

Private Attributes
TClonesArray *	fClusterArray
TClonesArray *	fDigiArray

Detailed Description

Task to remove duplicate crystals.

In timebased cluster reconstruction, same detectorId might appear multiple times in the same cluster, which could cause undefined behavior in classical (eventbased) routines. The class provides a workaround, adding up the energies of all digis with the same origin (=crystal, detectorId). (One could also think off dumping the second hit, etc...)

Author

Ph. Mahlberg mahlb.nosp@m.erg@.nosp@m.hiskp.nosp@m..uni.nosp@m.-bonn.nosp@m..de

Definition at line 26 of file PndEmcClusterRemoveDuplCrys.h.

Constructor & Destructor Documentation

PndEmcClusterRemoveDuplCrys::PndEmcClusterRemoveDuplCrys

(

Int_t

verbose = 0

)

Definition at line 14 of file PndEmcClusterRemoveDuplCrys.cxx.

                                                                      :
        FairTask("PndEmcClusterRemoveDuplCrys", verbose) {
}

virtual PndEmcClusterRemoveDuplCrys::~PndEmcClusterRemoveDuplCrys ( )

inlinevirtual

Definition at line 36 of file PndEmcClusterRemoveDuplCrys.h.

{};

PndEmcClusterRemoveDuplCrys::PndEmcClusterRemoveDuplCrys ( const PndEmcClusterRemoveDuplCrys &  )

private

Member Function Documentation

void PndEmcClusterRemoveDuplCrys::Exec ( Option_t *  opt )

virtual

Virtual method Exec

Definition at line 49 of file PndEmcClusterRemoveDuplCrys.cxx.

References PndEmcDigi::AbsorbEnergy(), PndEmcCluster::DigiList(), fClusterArray, fDigiArray, PndEmcDigi::GetDetectorId(), PndEmcCluster::MemberDigiMap(), and PndEmcCluster::removeDigi().

                                                {

        for(Int_t iClu=0; iClu<fClusterArray->GetEntriesFast(); ++iClu) {
                PndEmcCluster* theCluster = static_cast<PndEmcCluster*>(fClusterArray->UncheckedAt(iClu));
                
                for(std::vector<Int_t>::reverse_iterator rit = theCluster->DigiList().rbegin(); rit!=theCluster->DigiList().rend(); ) {

                        PndEmcDigi* theDigi = static_cast<PndEmcDigi*>(fDigiArray->UncheckedAt(*rit));
                        Int_t detId = theDigi->GetDetectorId();

                        //check wether other digi with same detector Id was already included in cluster
                        //..we do a reverse iteration here, because only the first (i.e. with the smallest idx) digi out of all digis with same detId in fDigiArry will be seen in the MemberDigiMap
                        //(cf PndEmcMakeCluster::addDigi...)
                        
                        Int_t tmp_pos = theCluster->MemberDigiMap().find(detId)->second;
                        if(tmp_pos != *rit) { //in Member Map ..other digi is related to this detId...add digis up
                                PndEmcDigi* oldDigi = static_cast<PndEmcDigi*>(fDigiArray->UncheckedAt(tmp_pos));
                                oldDigi->AbsorbEnergy(*theDigi);

                                //std::cerr << "Absorbing energy (ptr|detId): (" << oldDigi << " | "<< oldDigi->GetDetectorId() << ")" << std::flush;
                                theCluster->removeDigi(fDigiArray, (++rit).base());     //remove theDigi from cluster
                        } else {
                                ++rit;
                        }
                }
        }
}

InitStatus PndEmcClusterRemoveDuplCrys::Init ( )

virtual

Virtual method Init

Definition at line 19 of file PndEmcClusterRemoveDuplCrys.cxx.

References fClusterArray, and fDigiArray.

                                             {

        FairRootManager* ioman = FairRootManager::Instance();
        if (!ioman) {
                cout << "-E- PndEmcClusterRemoveDuplCrys::Init: "
                        << "RootManager not instantiated!" << endl;
                return kFATAL;
        }

        // Get inputarray
        fClusterArray = (TClonesArray*) ioman->GetObject("EmcCluster");
        if (!fClusterArray) {
                cout << "-W- PndEmcClusterRemoveDuplCrys::Init: "
                        << "No PndEmcHit array! Needed for MC Truth" << endl;
        }

        // Get input/output array
        fDigiArray = FairRunAna::Instance()->IsTimeStamp() ? (TClonesArray*) ioman->GetObject("EmcDigiClusterBase") : (TClonesArray*) ioman->GetObject("EmcDigi");
        if (!fDigiArray) {
                cout << "-W- PndEmcClusterRemoveDuplCrys::Init: "
                        << "No PndEmcDigi array!" << endl;
                return kERROR;
        }

        cout <<"-I- PndEmcClusterRemoveDuplCrys: "
                << "Got both arrays" << endl;

        return kSUCCESS;
}

PndEmcClusterRemoveDuplCrys& PndEmcClusterRemoveDuplCrys::operator= ( const PndEmcClusterRemoveDuplCrys &  )

private

Member Data Documentation

TClonesArray* PndEmcClusterRemoveDuplCrys::fClusterArray

private

Get parameter containers Verbosity level

Definition at line 61 of file PndEmcClusterRemoveDuplCrys.h.

Referenced by Exec(), and Init().

TClonesArray* PndEmcClusterRemoveDuplCrys::fDigiArray

private

Definition at line 62 of file PndEmcClusterRemoveDuplCrys.h.

Referenced by Exec(), and Init().

---

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

• PndEmcClusterRemoveDuplCrys.h
• PndEmcClusterRemoveDuplCrys.cxx