PndEmcMakeCluster Class Reference

Task to cluster PndEmcDigis. More...

#include <PndEmcMakeCluster.h>

Inheritance diagram for PndEmcMakeCluster:

Public Member Functions
	PndEmcMakeCluster (Int_t verbose=0, Bool_t storeclusters=kTRUE)
virtual	~PndEmcMakeCluster ()
virtual InitStatus	Init ()
	Init Task. More...
virtual void	Exec (Option_t *opt)
	Runs the task. More...
virtual void	ExecuteTasks (Option_t *option)
	Finish clusters after subtasks have been executed. More...
void	SetStorageOfData (Bool_t val)
	Method to specify whether clusters are stored or not. More...
void	StoreClusterBaseDigis (Bool_t val=kTRUE)
	Method to specify whether underlying digis are stored or not. More...
void	SetClusterActiveTime (Double_t time)
	Override EmcRecoPar cluster active time parameter ..to be set in ns!!! More...
void	SetPersistency (Bool_t val=kTRUE)
Bool_t	GetPersistency ()

Protected Member Functions
virtual void	SetParContainers ()
virtual void	FinishClusters ()
	Calls FinishCluster() for each cluster. More...

Private Member Functions
void	FinishCluster (PndEmcCluster *tmpcluster)
	Assign final parameters to cluster. More...
bool	HasExpired (PndEmcDigi *latestDigi, PndEmcCluster *theCluster, Int_t clusterIdx)
	Finishes clusters in timebased analysis. More...
void	cleansortmclist (std::vector< Int_t > &newlist, TClonesArray *mcTrackArray)
	Helper function, does not depend on class, identical to the one in PndEmcHitProducer. More...
	PndEmcMakeCluster (const PndEmcMakeCluster &)
PndEmcMakeCluster &	operator= (const PndEmcMakeCluster &)

Private Attributes
TClonesArray *	fDigiArray
TClonesArray *	fHitArray
TClonesArray *	fMCTrackArray
TClonesArray *	fClusterArray
	active clusters More...
TClonesArray *	fWriteOutArray
	expired clusters More...
BinaryFunctor *	fDigiFunctor
std::vector< PndEmcCluster * >	fClusterList
Double_t	fDigiEnergyTresholdBarrel
Double_t	fDigiEnergyTresholdFWD
Double_t	fDigiEnergyTresholdBWD
Double_t	fDigiEnergyTresholdShashlyk
Double_t	fClusterActiveTime
	Defines how long clusters are kept open in timebased reconstruction. More...
std::vector< Double_t >	fClusterPosParam
Int_t	fMapVersion
PndEmcGeoPar *	fGeoPar
PndEmcDigiPar *	fDigiPar
PndEmcRecoPar *	fRecoPar
Int_t	fVerbose
Bool_t	fStoreClusters
Bool_t	fStoreClusterBase
TStopwatch	fTimer

Static Private Attributes
static Int_t	fEventCounter =0

Detailed Description

Task to cluster PndEmcDigis.

Task to cluster PndEmcDigis. In the eventbased version clusters are built within one event, all adjacent digis go in one cluster. The timebased version assigns an active time to every cluster. Once expired, no digis can be added any more to cluster. The spatial cluster condition remains the same as in the eventbased reconstruction

Definition at line 36 of file PndEmcMakeCluster.h.

Constructor & Destructor Documentation

PndEmcMakeCluster::PndEmcMakeCluster	(	Int_t	verbose = 0,
		Bool_t	storeclusters = kTRUE
	)

Definition at line 49 of file PndEmcMakeCluster.cxx.

References fClusterList, fClusterPosParam, and PndPersistencyTask::SetPersistency().

                                                                       :
PndPersistencyTask("EmcClusteringTask", verbose),
fDigiArray(NULL), fHitArray(NULL), fMCTrackArray(NULL), fClusterArray(NULL), fWriteOutArray(NULL), fDigiFunctor(NULL), fClusterList(), fDigiEnergyTresholdBarrel(0), fDigiEnergyTresholdFWD(0), fDigiEnergyTresholdBWD(0), fDigiEnergyTresholdShashlyk(0), fClusterPosParam(), fMapVersion(0), fGeoPar(new PndEmcGeoPar()), fDigiPar(new PndEmcDigiPar()), fRecoPar(new PndEmcRecoPar()), fVerbose(verbose), fStoreClusters(storeclusters)
{
        fClusterList.clear();
        fClusterPosParam.clear();
        SetPersistency(storeclusters);
}

PndEmcMakeCluster::~PndEmcMakeCluster ( )

virtual

Definition at line 61 of file PndEmcMakeCluster.cxx.

{ 
}

PndEmcMakeCluster::PndEmcMakeCluster ( const PndEmcMakeCluster &  )

private

Member Function Documentation

void PndEmcMakeCluster::cleansortmclist ( std::vector< Int_t > &  newlist, TClonesArray *  mcTrackArray  )

private

Helper function, does not depend on class, identical to the one in PndEmcHitProducer.

Currently not used.

Parameters

newlist	...
mcTrackArray	...

Returns

void

Definition at line 432 of file PndEmcMakeCluster.cxx.

References PndMCTrack::GetMotherID(), and pt().

{
        std::vector <Int_t> tmplist;
        // Sort list...
        std::sort( newlist.begin(), newlist.end());
        // and copy every id only once (even though it might be in the list several times)
        std::unique_copy( newlist.begin(), newlist.end(), std::back_inserter( tmplist ) );

        // Now check if mother or (grand)^x-mother are already in the list
        // (which means i am a secondary)... if so, remove myself
        for(Int_t j=tmplist.size()-1; j>=0; j--){
                bool flag;
                PndMCTrack *pt;
                pt=((PndMCTrack*)mcTrackArray->At(tmplist[j]));
                if(pt->GetMotherID()<0) continue;
                flag=false;
                while(!flag){
                        Int_t id;
                        id=pt->GetMotherID();
                        if(id<0) break;
                        pt=(PndMCTrack*)mcTrackArray->At(id);

                        for(Int_t k=j-1; k>=0; k--){
                                if(tmplist[k]==id){
                                        tmplist.erase(tmplist.begin()+j);
                                        flag=true;
                                        break;
                                }
                        }
                }
        }
        newlist=tmplist;
}

void PndEmcMakeCluster::Exec ( Option_t *  opt )

virtual

Runs the task.

The task loops over the digis and adds neighboring digis (PndEmcTwoCoordIndex::IsNeighbour()) to clusters.

Parameters

opt

unused

Returns

void

Definition at line 161 of file PndEmcMakeCluster.cxx.

References PndEmcCluster::addCluster(), PndEmcCluster::addDigi(), PndEmcCluster::AddTracksEnteringExiting(), fClusterActiveTime, fClusterArray, fDigiArray, fDigiEnergyTresholdBarrel, fDigiEnergyTresholdBWD, fDigiEnergyTresholdFWD, fDigiEnergyTresholdShashlyk, fDigiFunctor, fTimer, fVerbose, fWriteOutArray, PndEmcDigi::GetEnergy(), PndEmcDigi::GetModule(), PndEmcHit::GetTrackEntering(), PndEmcHit::GetTrackExiting(), HasExpired(), hit(), i, PndEmcCluster::isInCluster(), PndEmcCluster::SetTrackEntering(), and PndEmcCluster::SetTrackExiting().

{

        if (fVerbose>2){
                fTimer.Start();
        }
        // Reset output array
        if ( ! fClusterArray ) Fatal("Exec", "No Cluster Array");

        fWriteOutArray->Delete();

//      std::cout << "------ Event " << FairRootManager::Instance()->GetEntryNr() << " ------" << std::endl;

        if(FairRunAna::Instance()->IsTimeStamp()) {
                fDigiArray->Delete();
                fDigiArray->AbsorbObjects(FairRootManager::Instance()->GetData("EmcDigiSorted", fDigiFunctor, fClusterActiveTime));
        }       


        Int_t nDigis = fDigiArray->GetEntriesFast();
        if (fVerbose>2){
                cout<<"DigiList length "<<nDigis<<endl;
        }

        //loop to build Cluster
        for (Int_t iDigi=0; iDigi<nDigis; iDigi++)
        {
                PndEmcDigi* theDigi = (PndEmcDigi*) fDigiArray->At(iDigi);
//              std::cout << std::endl << "DigiArray: " << iDigi  << std::endl;
//              theDigi->Print();
//              std::cout << *theDigi->GetPointerToLinks() << std::endl;
//              std::cout << std::endl;
                Int_t module=theDigi->GetModule();

                // In the following lines there is separate threshold for forward endcap
                // and the same for barrel and backward and shashlyk EMC, but for last 2 threshold probably should be different 
                if ((module==1||module==2)&&(theDigi->GetEnergy()< fDigiEnergyTresholdBarrel)) continue;
                if ((module==3)&&(theDigi->GetEnergy()< fDigiEnergyTresholdFWD)) continue;
                if ((module==4)&&(theDigi->GetEnergy()< fDigiEnergyTresholdBWD)) continue;
                if ((module==5)&&(theDigi->GetEnergy()< fDigiEnergyTresholdShashlyk)) continue;

                bool isAdded = false;
                Int_t clustmarker=0;
                Int_t clustLength=fClusterArray->GetEntriesFast();

                for(Int_t i=0;i<clustLength;i++)
                {
                        PndEmcCluster* cluster= (PndEmcCluster*) fClusterArray->At(i);
                        if(HasExpired(theDigi, cluster, i)) {
                                clustLength--;
                                i--;
                                continue;
                        } else {
//                                      std::cout << "Clusterarray: " << i << std::endl;
                        }

                        if(cluster->isInCluster(theDigi, fDigiArray))
                        {
                                if(!isAdded) {
                                        clustmarker=i;
                                        isAdded=true;
                                        cluster->addDigi(fDigiArray, iDigi);
                                        //cluster->AddLink(FairLink("EmcDigi", iDigi));
                                        FairMultiLinkedData hitLinks = theDigi->GetLinksWithType(FairRootManager::Instance()->GetBranchId("EmcHit"));
                                        for (Int_t j = 0; j < hitLinks.GetNLinks(); j++){
                                                PndEmcHit* hit = (PndEmcHit*)FairRootManager::Instance()->GetCloneOfLinkData(hitLinks.GetLink(j));
                                                //std::cout << "Hit: " << hit->GetDetectorID() << std::endl;
                                                if(hit) {
//                                                      std::cout << "Cluster : " << clustmarker << " Hit to add: " << hitLinks.GetLink(j) << std::endl;
                                                        if (cluster->GetLinks().count(hitLinks.GetLink(j)) == 0){
                                                                cluster->AddTracksEnteringExiting(hit->GetTrackEntering(), hit->GetTrackExiting());
                                                                cluster->AddLink(hitLinks.GetLink(j));
//                                                              std::cout << "Links in Cluster: " << *cluster->GetPointerToLinks() << std::endl;
//                                                              std::cout << "EnteringExiting for Hit: " << hitLinks.GetLink(j) << std::endl;
//                                                              std::cout << "TrackEntering: " << hit->GetTrackEntering() << std::endl;
//                                                              std::cout << "TrackExiting: " << hit->GetTrackExiting() << std::endl;
//                                                              std::cout << "Resulting Enter: " << cluster->GetTrackEntering() << std::endl;
//                                                              std::cout << "Resulting Exit: " << cluster->GetTrackExiting() << std::endl;
                                                        } else {
                                                                //std::cout << "Hit Already exists!" << std::endl;
                                                        }
                                                        delete(hit);

                                                } else {
                                                        std::cout << "-E in PndEmcMakeCluster::Exec FairLink " << hitLinks.GetLink(j) << "to EmcHit delivers null" << std::endl;
                                                }
                                        }
                                }
                                else
                                {
                                        PndEmcCluster* clust_clustmarker=(PndEmcCluster*) fClusterArray->At(clustmarker);
                                        PndEmcCluster* clust_i=(PndEmcCluster*) fClusterArray->At(i);
                                        clust_clustmarker->addCluster(clust_i, fDigiArray);
                                        fClusterArray->RemoveAt(i);
                                        fClusterArray->Compress();
                                        clustLength--;
                                        i--;
                                }
                        }
                }

                if (!isAdded)
                {
                        PndEmcCluster* newcluster = new((*fClusterArray)[clustLength]) PndEmcCluster();
                        newcluster->addDigi(fDigiArray, iDigi);

                        FairMultiLinkedData hitLinks = theDigi->GetLinksWithType(FairRootManager::Instance()->GetBranchId("EmcHit"));
                        //std::cout << "HitLinks isNotAdded:  " << hitLinks << std::endl;
                //      theDigi->Print(); std::cout << std::endl;
                        for (Int_t i = 0; i < hitLinks.GetNLinks(); i++){
                                PndEmcHit* hit = (PndEmcHit*)FairRootManager::Instance()->GetCloneOfLinkData(hitLinks.GetLink(i));
                                if(hit) {
//                                      std::cout << "NewCluster : " << clustLength << " Hit Added: " << hitLinks.GetLink(i) << std::endl;
//                                      std::cout << "TrackEntering: " << hit->GetTrackEntering() << std::endl;
//                                      std::cout << "TrackExiting: " << hit->GetTrackExiting() << std::endl;
                                        newcluster->SetTrackEntering(hit->GetTrackEntering());
                                        newcluster->SetTrackExiting(hit->GetTrackExiting());
                                        newcluster->SetInsertHistory(kFALSE);
                                        newcluster->AddLink(hitLinks.GetLink(i));
//                                      std::cout << "Resulting Enter: " << newcluster->GetTrackEntering() << std::endl;
//                                      std::cout << "Resulting Exit: " << newcluster->GetTrackExiting() << std::endl;
//                                      std::cout << "Links in Cluster: " << *newcluster->GetPointerToLinks() << std::endl;
                                        delete(hit);
                                } else {
                                        std::cout << "-E in PndEmcMakeCluster::Exec FairLink " << hitLinks.GetLink(i) << "to EmcHit delivers null" << std::endl;
                                }
                        }
                        //newcluster->SetLink(FairLink("EmcDigi", iDigi));
                }
        }
//      for (int k = 0; k < fClusterArray->GetEntriesFast(); k++){
//              PndEmcCluster* myCluster = (PndEmcCluster*)fClusterArray->At(k);
//              std::cout << k << " : Entering: " << myCluster->GetTrackEntering() << std::endl;
//              std::cout << k << " : Exiting: " << myCluster->GetTrackExiting() << std::endl;
//      }
//      std::cout << std::endl;

        fWriteOutArray->AbsorbObjects(fClusterArray);
}

virtual void PndEmcMakeCluster::ExecuteTasks ( Option_t *  option )

inlinevirtual

Finish clusters after subtasks have been executed.

Definition at line 48 of file PndEmcMakeCluster.h.

References FinishClusters().

                                              {
                PndPersistencyTask::ExecuteTasks(option);
                FinishClusters();
  }

void PndEmcMakeCluster::FinishCluster ( PndEmcCluster *  tmpclust )

private

Assign final parameters to cluster.

Calculates properties of the cluster and sets them.

Assign final parameters to cluster. Subtasks might introduce changes in cluster composition. To keep track of them, function is called after subtasks have been executed

The calculated properties include the energy, position, timestamp and Zernike moments.

Parameters

tmpclust

The cluster to work on

Returns

void

Definition at line 310 of file PndEmcMakeCluster.cxx.

References PndEmcXClMoments::AbsZernikeMoment(), PndEmcClusterProperties::Energy(), fClusterPosParam, fDigiArray, fRecoPar, PndEmcRecoPar::GetEmcClusterPosMethod(), PndEmcCluster::GetTrackEntering(), PndEmcXClMoments::Lat(), PndEmcCluster::Maxima(), PndEmcCluster::SetEnergy(), PndEmcCluster::SetLatMom(), PndEmcCluster::SetPosition(), PndEmcCluster::SetZ20(), PndEmcCluster::SetZ53(), and PndEmcClusterProperties::Where().

Referenced by FinishClusters().

{
        using namespace std;

        //cout << "\tfinalizing cluster: ";
        //cout << "\tmaxTime: " << tmpclust->GetTimeStamp() << endl;
        //cout << "\tcontains the following digis: " << endl;
        //for(unsigned int i=0; i<tmpclust->DigiList().size(); ++i) {
        //cout << "\t#:" << i << "\ttime: " << ((PndEmcDigi*) (fDigiArray->UncheckedAt(tmpclust->DigiList()[i])))->GetTimeStamp() << "\tenergy: " << ((PndEmcDigi*) (fDigiArray->UncheckedAt(tmpclust->DigiList()[i])))->GetEnergy() << "\tidx: " << ((PndEmcDigi*) (fDigiArray->UncheckedAt(tmpclust->DigiList()[i])))->GetDetectorId() << endl;
        //}

        PndEmcClusterProperties clustProperties(*tmpclust, fDigiArray);

        tmpclust->SetEnergy(clustProperties.Energy());
        tmpclust->SetTimeStamp(tmpclust->Maxima(fDigiArray)->GetTimeStamp());
        TVector3 tmpbumppos = clustProperties.Where(fRecoPar->GetEmcClusterPosMethod(), fClusterPosParam);
        tmpclust->SetPosition(tmpbumppos);
        PndEmcXClMoments xClMoments(*tmpclust, fDigiArray);
        tmpclust->SetZ20(xClMoments.AbsZernikeMoment(2, 0, 15));
        tmpclust->SetZ53(xClMoments.AbsZernikeMoment(5, 3, 15));
        tmpclust->SetLatMom(xClMoments.Lat());
        tmpclust->SetLinks(tmpclust->GetTrackEntering());
        //      std::cout << "ClusterOutput: " << std::endl;
        //      tmpclust->Print();
        //const std::vector<Int_t>& MCTruth = tmpclust->GetMcList();
        //std::cout<<"The cluster #"<<i<<" produced by MC Track #";
        //for(Int_t j=0;j<MCTruth.size();++j)
        //      cout<<MCTruth[j]<<' ';
        //cout<<endl;

        //tmpclust->fMcList.clear();
        //if(fHitArray && fMCTrackArray && !FairRunAna::Instance()->IsTimeStamp()){
        //      // BS: this is a first order approximation only !!!!

        //      std::vector <Int_t> newlist;
        //      newlist.clear();
        //      for(Int_t j=0; j<tmpclust->fDigiList.size(); j++){
        //              PndEmcDigi*m;
        //              m=(PndEmcDigi*)fDigiArray->At(tmpclust->fDigiList[j]);

        //              Int_t inx;
        //              inx=m->GetHitIndex();

        //              if(inx>=0){
                //                      if(fHitArray->At(inx) != 0 ){//add by hujf
        //                      const std::vector <Int_t> &tmplist=((PndEmcHit*) fHitArray->At(inx))->GetMcList();

        //                      // I copy the complete list instead of only the highest energy particle
        //                      // highest energy might be a problem for Hits which belong to two real clusters?
        //                      for(Int_t k=0; k<tmplist.size(); k++){
        //                              newlist.push_back(tmplist[k]);
        //                      }
        //              }
        //      }
        //      // if( fMCTrackArray) checked above already
        //      cleansortmclist(newlist,fMCTrackArray);
        //      tmpclust->fMcList=newlist;
        //}
}

void PndEmcMakeCluster::FinishClusters ( )

protectedvirtual

Calls FinishCluster() for each cluster.

Returns

void

Definition at line 376 of file PndEmcMakeCluster.cxx.

References ctime, Double_t, fEventCounter, FinishCluster(), fTimer, fVerbose, fWriteOutArray, i, and rtime.

Referenced by ExecuteTasks().

{
        Int_t nCluster = fWriteOutArray->GetEntriesFast();
        for (Int_t i=0; i<nCluster; i++) {
                FinishCluster((PndEmcCluster*) (fWriteOutArray->At(i)));
        }

        fEventCounter++;
        if (fVerbose>0)
                cout<<"PndEmcMakeCluster, event: "<<fEventCounter<<endl;

        if (fVerbose>2){
                fTimer.Stop();
                Double_t rtime = fTimer.RealTime();
                Double_t ctime = fTimer.CpuTime();
                fTimer.Reset();
                cout << "PndEmcMakeCluster, Real time " << rtime << " s, CPU time " << ctime << " s" << endl;
        }       
}

Bool_t PndPersistencyTask::GetPersistency ( )

inlineinherited

Definition at line 32 of file PndPersistencyTask.h.

References PndPersistencyTask::fPersistency.

Referenced by PndLmdPixelHitProducerFast::GetPersistance(), PndMdtDigitization::Init(), PndMdtHitProducerIdeal::Init(), PndMdtClusterTask::Init(), PndFtsHitProducerRealFast::Init(), PndDiscTaskReconstruction::Init(), PndRichHitProducer::Init(), PndSttHitProducerRealFast::Init(), PndSttHelixHitProducer::Init(), PndDiscTaskPID::Init(), PndIdealTrackFinder::Init(), PndSttMvdGemTracking::Init(), PndMdtTrkProducer::Init(), PndFtsHitProducerRealFull::Init(), PndLmdPixelClusterTask::Init(), PndSttHitProducerRealFull::Init(), PndLmdStripClusterTask::Init(), PndEmcApdHitProducer::Init(), PndMissingPzCleanerTask::Init(), PndEmcMakeRecoHit::Init(), PndEmcMakeClusterOnline::Init(), PndTrackSmearTask::Init(), PndEmcFWEndcapTimebasedWaveforms::Init(), PndSttHitProducerIdeal::Init(), PndEmcFWEndcapDigi::Init(), PndFtsHitProducerIdeal::Init(), Init(), PndMdtPointsToWaveform::Init(), PndDiscTaskDigitization::Init(), PndEmcMakeDigi::Init(), PndSdsTimeWalkCorrTask::Init(), PndLmdPixelHitProducerFast::Init(), PndDrcHitFinder::Init(), PndRichHitFinder::Init(), PndEmcMakeCorr::Init(), PndFtofHitProducerIdeal::Init(), PndEmcHitsToWaveform::Init(), PndSciTDigiTask::Init(), PndDrcHitProducerIdeal::Init(), PndSdsHitProducerIdeal::Init(), PndSciTHitProducerIdeal::Init(), PndRecoMultiKalmanTask2::Init(), PndEmcHitProducer::Init(), PndDrcHitProducerReal::Init(), PndDskFLGHitProducerIdeal::Init(), PndEmcTmpWaveformToDigi::Init(), PndDrcDigiTask::Init(), PndEmcWaveformToDigi::Init(), PndSttMatchTracks::Init(), PndEmcWaveformToCalibratedDigi::Init(), PndTrkTracking2::Init(), PndSttFindTracks::Init(), PndEmcMultiWaveformToCalibratedDigi::Init(), PndRecoKalmanTask2::Init(), PndDrcTimeDigiTask::Init(), PndEmcExpClusterSplitter::Init(), PndFtsHoughTrackerTask::Init(), PndSdsNoiseProducer::Init(), PndEmcPhiBumpSplitter::Init(), PndSdsIdealRecoTask::Init(), PndSdsHybridHitProducer::Init(), PndRecoMultiKalmanTask::Init(), PndSdsIdealClusterTask::Init(), PndRecoKalmanTask::Init(), PndSdsStripHitProducerDif::Init(), PndGemDigitize::Init(), PndSdsStripHitProducer::Init(), PndGemFindHits::Init(), PndSdsPixelClusterTask::Init(), PndSdsStripClusterTask::Init(), PndMvdGemTrackFinderOnHits::Init(), PndBarrelTrackFinder::Init(), PndEmcFullDigiTask::PndEmcFullDigiTask(), PndEmcMakeBump::PndEmcMakeBump(), PndUnassignedHitsTask::RegisterBranches(), PndMvdClusterTask::SetPersistance(), PndMvdDigiTask::SetPersistance(), PndEmcMakeBump::SetStorageOfData(), and PndEmcFullDigiTask::StoreDigi().

{ return fPersistency; }

bool PndEmcMakeCluster::HasExpired ( PndEmcDigi *  latestDigi, PndEmcCluster *  theCluster, Int_t  clusterIdx  )

private

Finishes clusters in timebased analysis.

Checks if the timestamp of latestDigi is later than the active time of theCluster. If yes, the cluster is written and true is returned.

Parameters

latestDigi	Current digi to check timestamp of.
theCluster	Current cluster for which the expiration is checked.
clusterIdx	Index of theCluster in the cluster TClonesArray.

Returns

bool

Return values

true	Cluster has expired and was written.
false	Cluster not expired.

Definition at line 409 of file PndEmcMakeCluster.cxx.

References fClusterActiveTime, fClusterArray, and fWriteOutArray.

Referenced by Exec().

{
        if(!FairRunAna::Instance()->IsTimeStamp()) return false;

        if(latestDigi->GetTimeStamp() > theCluster->GetTimeStamp() + fClusterActiveTime) {
                fWriteOutArray->AbsorbObjects(fClusterArray, clusterIdx, clusterIdx);
                //cout << "cluster " << clusterIdx << "has been expired" << std::endl;
                return true;
        } else {
                return false;
        }       
}

InitStatus PndEmcMakeCluster::Init ( )

virtual

Init Task.

Prepares the TClonesArray of PndEmcDigi for reading and of PndEmcCluster for writing.

Returns

InitStatus

Return values

kSUCCESS

success

Definition at line 74 of file PndEmcMakeCluster.cxx.

References fClusterActiveTime, fClusterArray, fClusterPosParam, fDigiArray, fDigiEnergyTresholdBarrel, fDigiEnergyTresholdBWD, fDigiEnergyTresholdFWD, fDigiEnergyTresholdShashlyk, fDigiFunctor, fGeoPar, fHitArray, fMCTrackArray, fRecoPar, fStoreClusterBase, fWriteOutArray, PndEmcRecoPar::GetClusterActiveTime(), PndEmcRecoPar::GetEmcClusterPosMethod(), PndEmcRecoPar::GetEnergyThresholdBarrel(), PndEmcRecoPar::GetEnergyThresholdBWD(), PndEmcRecoPar::GetEnergyThresholdFWD(), PndEmcRecoPar::GetEnergyThresholdShashlyk(), PndEmcRecoPar::GetOffsetParmA(), PndEmcRecoPar::GetOffsetParmB(), PndEmcRecoPar::GetOffsetParmC(), PndPersistencyTask::GetPersistency(), PndEmcGeoPar::InitEmcMapper(), and PndEmcStructure::Instance().

{

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

        // Get input array
        fDigiArray = (TClonesArray*) ioman->GetObject("EmcDigi");
        if ( ! fDigiArray ) {
                fDigiArray = (TClonesArray*) ioman->GetObject("EmcDigiSorted");
                if ( ! fDigiArray ) {
                        cout << "-W- PndEmcMakeCluster::Init: "
                                << "No PndEmcDigi array!" << endl;
                        return kERROR;
                }
        }

        // Get input array
        fHitArray = (TClonesArray*) ioman->GetObject("EmcHit");
        if ( ! fHitArray ) {
                cout << "-W- PndEmcMakeCluster::Init: "
                        << "No PndEmcHit array! Needed for MC Truth" << endl;
        }

        // Get input array
        fMCTrackArray = (TClonesArray*) ioman->GetObject("MCTrack");
        if ( ! fMCTrackArray ) {
                cout << "-W- PndEmcMakeCluster::Init: "
                        << "No MCTrack array! Needed for MC Truth" << endl;
        }

        // Create and register output array
        fClusterArray = new TClonesArray("PndEmcCluster");
        fWriteOutArray = new TClonesArray(fClusterArray->GetClass());

        ioman->Register("EmcCluster","Emc", fWriteOutArray, GetPersistency());

        if(FairRunAna::Instance()->IsTimeStamp()) {     
                // for subsequent methods we need the "event grouping" as given by the TS buffer --> fDigiArray becomes an output array. 
                fDigiArray = new TClonesArray(fDigiArray->GetClass());          //still needed to activate array 
                ioman->Register("EmcDigiClusterBase", "Emc", fDigiArray, fStoreClusterBase);
        }
        fDigiFunctor = new TimeGap();

        fGeoPar->InitEmcMapper();  
        PndEmcStructure::Instance();

        fDigiEnergyTresholdBarrel=fRecoPar->GetEnergyThresholdBarrel();
        fDigiEnergyTresholdFWD=fRecoPar->GetEnergyThresholdFWD();
        fDigiEnergyTresholdBWD=fRecoPar->GetEnergyThresholdBWD();
        fDigiEnergyTresholdShashlyk=fRecoPar->GetEnergyThresholdShashlyk();
        fClusterActiveTime=fRecoPar->GetClusterActiveTime() * 1.0e9; //convert from seconds to nanoseconds

        cout<<"PndEmcMakeCluster::fDigiEnergyTresholdBarrel: "<<fDigiEnergyTresholdBarrel<<endl;
        cout<<"PndEmcMakeCluster::fDigiEnergyTresholdFWD: "<<fDigiEnergyTresholdFWD<<endl;
        cout<<"PndEmcMakeCluster::fDigiEnergyTresholdBWD: "<<fDigiEnergyTresholdBWD<<endl;
        cout<<"PndEmcMakeCluster::fDigiEnergyTresholdShashlyk: "<<fDigiEnergyTresholdShashlyk<<endl;
        cout<<"PndEmcMakeCluster::fClusterActiveTime in ns(!): "<<fClusterActiveTime<<endl;
        cout<<"PndEmcMakeCluster::ClusterPosMethod(): " << fRecoPar->GetEmcClusterPosMethod() << std::endl;

        if (!strcmp(fRecoPar->GetEmcClusterPosMethod(),"lilo"))
        {
                cout<<"Lilo cluster position method"<<endl;
                fClusterPosParam.push_back(fRecoPar->GetOffsetParmA());
                fClusterPosParam.push_back(fRecoPar->GetOffsetParmB());
                fClusterPosParam.push_back(fRecoPar->GetOffsetParmC());
        }       

        cout << "-I- PndEmcMakeCluster: Intialization successfull" << endl;
        return kSUCCESS;
}

PndEmcMakeCluster& PndEmcMakeCluster::operator= ( const PndEmcMakeCluster &  )

private

void PndEmcMakeCluster::SetClusterActiveTime ( Double_t  time )

inline

Override EmcRecoPar cluster active time parameter ..to be set in ns!!!

Definition at line 62 of file PndEmcMakeCluster.h.

References fClusterActiveTime.

void PndEmcMakeCluster::SetParContainers ( )

protectedvirtual

Get parameter containers

Definition at line 466 of file PndEmcMakeCluster.cxx.

References fDigiPar, fGeoPar, fRecoPar, and run.

                                         {

        // Get run and runtime database
        FairRun* run = FairRun::Instance();
        if ( ! run ) Fatal("SetParContainers", "No analysis run");

        FairRuntimeDb* db = run->GetRuntimeDb();
        if ( ! db ) Fatal("SetParContainers", "No runtime database");

        // Get Emc geometry parameter container
        fGeoPar = (PndEmcGeoPar*) db->getContainer("PndEmcGeoPar");

        // Get Emc digitisation parameter container
        fDigiPar = (PndEmcDigiPar*) db->getContainer("PndEmcDigiPar");

        // Get Emc reconstruction parameter container
        fRecoPar = (PndEmcRecoPar*) db->getContainer("PndEmcRecoPar");
}

void PndPersistencyTask::SetPersistency ( Bool_t  val = kTRUE )

inlineinherited

Definition at line 31 of file PndPersistencyTask.h.

References PndPersistencyTask::fPersistency, and val.

Referenced by barrelTrackFinder(), digi_complete(), digi_complete_newSTT(), digiOnly_complete(), PndBarrelTrackFinder::PndBarrelTrackFinder(), PndCATracking::PndCATracking(), PndDrcHitFinder::PndDrcHitFinder(), PndEmc2DLocMaxFinder::PndEmc2DLocMaxFinder(), PndEmcExpClusterSplitter::PndEmcExpClusterSplitter(), PndEmcFullDigiTask::PndEmcFullDigiTask(), PndEmcFWEndcapDigi::PndEmcFWEndcapDigi(), PndEmcFWEndcapTimebasedWaveforms::PndEmcFWEndcapTimebasedWaveforms(), PndEmcHitProducer::PndEmcHitProducer(), PndEmcHitsToWaveform::PndEmcHitsToWaveform(), PndEmcMakeBump::PndEmcMakeBump(), PndEmcMakeCluster(), PndEmcMakeClusterOnline::PndEmcMakeClusterOnline(), PndEmcMakeDigi::PndEmcMakeDigi(), PndEmcMakeRecoHit::PndEmcMakeRecoHit(), PndEmcMultiWaveformToCalibratedDigi::PndEmcMultiWaveformToCalibratedDigi(), PndEmcPhiBumpSplitter::PndEmcPhiBumpSplitter(), PndEmcTmpWaveformToDigi::PndEmcTmpWaveformToDigi(), PndEmcWaveformToCalibratedDigi::PndEmcWaveformToCalibratedDigi(), PndEmcWaveformToDigi::PndEmcWaveformToDigi(), PndFtofHitProducerIdeal::PndFtofHitProducerIdeal(), PndFtsCATracking::PndFtsCATracking(), PndFtsHitProducerIdeal::PndFtsHitProducerIdeal(), PndFtsHitProducerRealFast::PndFtsHitProducerRealFast(), PndFtsHitProducerRealFull::PndFtsHitProducerRealFull(), PndFtsHoughTrackerTask::PndFtsHoughTrackerTask(), PndGemDigitize::PndGemDigitize(), PndGemFindHits::PndGemFindHits(), PndIdealTrackFinder::PndIdealTrackFinder(), PndLmdPixelClusterTask::PndLmdPixelClusterTask(), PndLmdPixelHitProducerFast::PndLmdPixelHitProducerFast(), PndMdtClusterTask::PndMdtClusterTask(), PndMdtDigitization::PndMdtDigitization(), PndMdtHitProducerIdeal::PndMdtHitProducerIdeal(), PndMdtPointsToWaveform::PndMdtPointsToWaveform(), PndMdtTrkProducer::PndMdtTrkProducer(), PndMissingPzCleanerTask::PndMissingPzCleanerTask(), PndMvdGemTrackFinderOnHits::PndMvdGemTrackFinderOnHits(), PndMvdHitProducerIdeal::PndMvdHitProducerIdeal(), PndMvdPixelClusterTask::PndMvdPixelClusterTask(), PndMvdTimeWalkCorrTask::PndMvdTimeWalkCorrTask(), PndMvdToPix4ClusterTask::PndMvdToPix4ClusterTask(), PndRecoKalmanTask::PndRecoKalmanTask(), PndRecoKalmanTask2::PndRecoKalmanTask2(), PndRecoMultiKalmanTask::PndRecoMultiKalmanTask(), PndRecoMultiKalmanTask2::PndRecoMultiKalmanTask2(), PndRichHitFinder::PndRichHitFinder(), PndRichHitProducer::PndRichHitProducer(), PndSciTDigiTask::PndSciTDigiTask(), PndSciTHitProducerIdeal::PndSciTHitProducerIdeal(), PndSdsHitProducerIdeal::PndSdsHitProducerIdeal(), PndSdsHybridHitProducer::PndSdsHybridHitProducer(), PndSdsIdealClusterTask::PndSdsIdealClusterTask(), PndSdsIdealRecoTask::PndSdsIdealRecoTask(), PndSdsNoiseProducer::PndSdsNoiseProducer(), PndSdsPixelClusterTask::PndSdsPixelClusterTask(), PndSdsStripClusterTask::PndSdsStripClusterTask(), PndSdsStripHitProducer::PndSdsStripHitProducer(), PndSdsTimeWalkCorrTask::PndSdsTimeWalkCorrTask(), PndSttFindTracks::PndSttFindTracks(), PndSttHelixHitProducer::PndSttHelixHitProducer(), PndSttHitProducerIdeal::PndSttHitProducerIdeal(), PndSttHitProducerRealFast::PndSttHitProducerRealFast(), PndSttHitProducerRealFull::PndSttHitProducerRealFull(), PndSttMatchTracks::PndSttMatchTracks(), PndSttMvdGemTracking::PndSttMvdGemTracking(), PndTrackSmearTask::PndTrackSmearTask(), PndTrkTracking2::PndTrkTracking2(), reco(), reco_complete(), reco_complete_gf2(), reco_complete_newSTT(), reco_complete_sec(), recoideal_complete(), PndMvdClusterTask::SetPersistance(), PndMvdDigiTask::SetPersistance(), PndLmdPixelHitProducerFast::SetPersistance(), PndSdsHitProducerIdeal::SetPersistance(), PndSttMvdGemTracking::SetPersistenc(), PndMdtClusterTask::SetPersistence(), PndSttHelixHitProducer::SetPersistence(), PndMissingPzCleanerTask::SetPersistence(), PndFtsHitProducerRealFast::SetPersistence(), PndFtsHitProducerRealFull::SetPersistence(), PndSttHitProducerRealFull::SetPersistence(), PndSttHitProducerIdeal::SetPersistence(), PndSttHitProducerRealFast::SetPersistence(), PndFtsHitProducerIdeal::SetPersistence(), PndTrackSmearTask::SetPersistence(), PndSciTHitProducerIdeal::SetPersistence(), PndIdealTrackFinder::SetPersistence(), PndSttMatchTracks::SetPersistence(), PndSttFindTracks::SetPersistence(), PndFtsHoughTrackerTask::SetPersistence(), PndTrkTracking2::SetPersistence(), PndEmcMakeRecoHit::SetStorageOfData(), PndEmcFWEndcapDigi::SetStorageOfData(), PndEmcMakeClusterOnline::SetStorageOfData(), PndEmcFWEndcapTimebasedWaveforms::SetStorageOfData(), PndEmcMakeDigi::SetStorageOfData(), PndMdtPointsToWaveform::SetStorageOfData(), PndEmc2DLocMaxFinder::SetStorageOfData(), SetStorageOfData(), PndEmcHitsToWaveform::SetStorageOfData(), PndEmcMakeBump::SetStorageOfData(), PndEmcTmpWaveformToDigi::SetStorageOfData(), PndEmcWaveformToDigi::SetStorageOfData(), PndEmcWaveformToCalibratedDigi::SetStorageOfData(), PndEmcMultiWaveformToCalibratedDigi::SetStorageOfData(), PndEmcExpClusterSplitter::SetStorageOfData(), PndEmcPhiBumpSplitter::SetStorageOfData(), standard_tracking(), and PndEmcFullDigiTask::StoreDigi().

{ fPersistency = val; }

void PndEmcMakeCluster::SetStorageOfData

(

Bool_t

val

)

Method to specify whether clusters are stored or not.

Definition at line 485 of file PndEmcMakeCluster.cxx.

References PndPersistencyTask::SetPersistency().

Referenced by digi_sttcombi(), emc(), emc_complete(), emc_complete_corr(), emc_correction_data_production(), emc_correction_QA_data_production(), and RecoComplete().

{
        SetPersistency(val);
        return;
}

void PndEmcMakeCluster::StoreClusterBaseDigis ( Bool_t  val = kTRUE )

inline

Method to specify whether underlying digis are stored or not.

Restoring digis makes only sense if reconstruction is done timebased

Definition at line 60 of file PndEmcMakeCluster.h.

References fStoreClusterBase, and val.

Referenced by emc_reco(), run_reco(), and run_reco_timebased().

{ fStoreClusterBase = val; }

Member Data Documentation

Double_t PndEmcMakeCluster::fClusterActiveTime

private

Defines how long clusters are kept open in timebased reconstruction.

Definition at line 101 of file PndEmcMakeCluster.h.

Referenced by Exec(), HasExpired(), Init(), and SetClusterActiveTime().

TClonesArray* PndEmcMakeCluster::fClusterArray

private

active clusters

Output array of PndEmcClusters

Definition at line 91 of file PndEmcMakeCluster.h.

Referenced by Exec(), HasExpired(), and Init().

std::vector<PndEmcCluster*> PndEmcMakeCluster::fClusterList

private

Definition at line 95 of file PndEmcMakeCluster.h.

Referenced by PndEmcMakeCluster().

std::vector<Double_t> PndEmcMakeCluster::fClusterPosParam

private

Definition at line 103 of file PndEmcMakeCluster.h.

Referenced by FinishCluster(), Init(), and PndEmcMakeCluster().

TClonesArray* PndEmcMakeCluster::fDigiArray

private

Input array of CbmDigis

Definition at line 84 of file PndEmcMakeCluster.h.

Referenced by Exec(), FinishCluster(), and Init().

Double_t PndEmcMakeCluster::fDigiEnergyTresholdBarrel

private

Definition at line 97 of file PndEmcMakeCluster.h.

Referenced by Exec(), and Init().

Double_t PndEmcMakeCluster::fDigiEnergyTresholdBWD

private

Definition at line 99 of file PndEmcMakeCluster.h.

Referenced by Exec(), and Init().

Double_t PndEmcMakeCluster::fDigiEnergyTresholdFWD

private

Definition at line 98 of file PndEmcMakeCluster.h.

Referenced by Exec(), and Init().

Double_t PndEmcMakeCluster::fDigiEnergyTresholdShashlyk

private

Definition at line 100 of file PndEmcMakeCluster.h.

Referenced by Exec(), and Init().

BinaryFunctor* PndEmcMakeCluster::fDigiFunctor

private

Definition at line 94 of file PndEmcMakeCluster.h.

Referenced by Exec(), and Init().

PndEmcDigiPar* PndEmcMakeCluster::fDigiPar

private

Geometry parameter container

Definition at line 109 of file PndEmcMakeCluster.h.

Referenced by SetParContainers().

Int_t PndEmcMakeCluster::fEventCounter =0

staticprivate

Definition at line 106 of file PndEmcMakeCluster.h.

Referenced by FinishClusters().

PndEmcGeoPar* PndEmcMakeCluster::fGeoPar

private

Definition at line 108 of file PndEmcMakeCluster.h.

Referenced by Init(), and SetParContainers().

TClonesArray* PndEmcMakeCluster::fHitArray

private

Input array of Hits and MC Tracks ... needed for MC

Definition at line 87 of file PndEmcMakeCluster.h.

Referenced by Init().

Int_t PndEmcMakeCluster::fMapVersion

private

Definition at line 104 of file PndEmcMakeCluster.h.

TClonesArray* PndEmcMakeCluster::fMCTrackArray

private

Definition at line 88 of file PndEmcMakeCluster.h.

Referenced by Init().

PndEmcRecoPar* PndEmcMakeCluster::fRecoPar

private

Digitisation parameter container

Definition at line 110 of file PndEmcMakeCluster.h.

Referenced by FinishCluster(), Init(), and SetParContainers().

Bool_t PndEmcMakeCluster::fStoreClusterBase

private

Definition at line 116 of file PndEmcMakeCluster.h.

Referenced by Init(), and StoreClusterBaseDigis().

Bool_t PndEmcMakeCluster::fStoreClusters

private

Definition at line 115 of file PndEmcMakeCluster.h.

TStopwatch PndEmcMakeCluster::fTimer

private

Definition at line 118 of file PndEmcMakeCluster.h.

Referenced by Exec(), and FinishClusters().

Int_t PndEmcMakeCluster::fVerbose

private

Reconstruction parameter container Verbosity level

Definition at line 113 of file PndEmcMakeCluster.h.

Referenced by Exec(), and FinishClusters().

TClonesArray* PndEmcMakeCluster::fWriteOutArray

private

expired clusters

Definition at line 92 of file PndEmcMakeCluster.h.

Referenced by Exec(), FinishClusters(), HasExpired(), and Init().

---

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

• PndEmcMakeCluster.h
• PndEmcMakeCluster.cxx