PndEmcCorrBump Class Reference

#include <PndEmcCorrBump.h>

Inheritance diagram for PndEmcCorrBump:

Public Member Functions
	PndEmcCorrBump (Int_t verbose=0, Bool_t storeclusters=kTRUE)
virtual	~PndEmcCorrBump ()
virtual InitStatus	Init ()
virtual void	Exec (Option_t *opt)
void	SetStorageOfData (Bool_t val)
void	FinishTask ()
void	RunTimeBased ()

Private Member Functions
virtual void	SetParContainers ()
	PndEmcCorrBump (const PndEmcCorrBump &L)
PndEmcCorrBump &	operator= (const PndEmcCorrBump &)

Private Attributes
TClonesArray *	fDigiArray
TClonesArray *	fSharedDigiArray
TClonesArray *	fBumpArray
Int_t	fMapVersion
PndEmcGeoPar *	fGeoPar
PndEmcDigiPar *	fDigiPar
PndEmcRecoPar *	fRecoPar
Int_t	fVerbose
Bool_t	fStoreClusters
Bool_t	fTimeOrderedDigi
Int_t	HowManyDigi
Int_t	HowManyCluster
PndEmcDigiCalibrator	digiCalibrator
TClonesArray *	fBumpArrayTBD
std::vector< Double_t >	fClusterPosParam

Static Private Attributes
static Int_t	fEventCounter =0

Detailed Description

Definition at line 28 of file PndEmcCorrBump.h.

Constructor & Destructor Documentation

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

Definition at line 73 of file PndEmcCorrBump.cxx.

                                                                 :
fDigiArray(0), fSharedDigiArray(0), fBumpArray(0), fMapVersion(0), fGeoPar(new PndEmcGeoPar()), fDigiPar(new PndEmcDigiPar()), fRecoPar(new PndEmcRecoPar()), fVerbose(verbose), fStoreClusters(storeclusters)
{
}

PndEmcCorrBump::~PndEmcCorrBump ( )

virtual

Definition at line 82 of file PndEmcCorrBump.cxx.

{ 
}

PndEmcCorrBump::PndEmcCorrBump ( const PndEmcCorrBump &  L )

private

Member Function Documentation

void PndEmcCorrBump::Exec ( Option_t *  opt )

virtual

Virtual method Exec

Definition at line 140 of file PndEmcCorrBump.cxx.

References CAMath::Abs(), PndEmcXClMoments::AbsZernikeMoment(), Bool_t, PndEmcDigiCalibrator::CalibrationEvtTimeByDigi(), ctime, digiCalibrator, PndEmcCluster::DigiList(), Double_t, PndEmcClusterProperties::Energy(), PndEmcCluster::energy(), fBumpArray, fBumpArrayTBD, fClusterPosParam, fDigiArray, PndEmcDigi::fDigiArrayTBD, fEventCounter, PndEmcDigi::fEvtNo, fRecoPar, fSharedDigiArray, fVerbose, PndEmcRecoPar::GetEmcClusterPosMethod(), PndEmcDigi::GetEnergy(), PndEmcDigi::GetModule(), PndEmcDigiCalibrator::GetTimeResolutionOfDigi(), HowManyCluster, HowManyDigi, i, PndEmcCluster::IsEnergyValid(), PndEmcCluster::IsPositionValid(), PndEmcXClMoments::Lat(), PndEmcCluster::NumberOfDigis(), PndEmcDigi::Print(), PndEmcCluster::removeDigi(), rtime, PndEmcCluster::SetEnergy(), PndEmcBump::SetEventNo(), PndEmcCluster::SetLatMom(), PndEmcCluster::SetPosition(), PndEmcCluster::SetZ20(), PndEmcCluster::SetZ53(), timer, and PndEmcClusterProperties::Where().

{
        TStopwatch timer;
        if (fVerbose>2){
                timer.Start();
        }
        // Reset output array
        if ( ! fBumpArray ) Fatal("Exec", "No Bump Array");

        Double_t fEventTime = FairRootManager::Instance()->GetEventTime();
        if (fVerbose>0){
                Int_t nDigis = fDigiArray->GetEntriesFast();
                Int_t nSharedDigis = fSharedDigiArray->GetEntriesFast();
                cout<<"*************************PndEmcCorrBump*****************************"<<endl;
                cout<<"Event NO. #"<<fEventCounter<<", EventTime#"<<fEventTime
                        <<", Digis#"<<nDigis<<", SharedDigis#"<<nSharedDigis<<endl;
                cout<<"********************************************************************"<<endl;
                cout<<"Before correction, bumps#"<<fBumpArray->GetEntriesFast()
                        <<", BumptoBeDetermined #"<<fBumpArrayTBD->GetEntriesFast()
                        <<", DigitoBeDetermined #"<<PndEmcDigi::fDigiArrayTBD->GetEntriesFast()<<endl;
                cout<<"********************************************************************"<<endl;
                //cout<<"------------> fDigiArray#"<<fDigiArray->GetEntriesFast()<<endl;
                //for(Int_t iDigi=0; iDigi < fDigiArray->GetEntriesFast(); ++iDigi){
                //      PndEmcDigi* theDigi = (PndEmcDigi*)fDigiArray->At(iDigi);
                //      theDigi->Print();
                //      cout<<endl;
                //}
                //cout<<"------------> fSharedDigiArray#"<<fSharedDigiArray->GetEntriesFast()<<endl;
                //for(Int_t iDigi=0; iDigi < fSharedDigiArray->GetEntriesFast(); ++iDigi){
                //      PndEmcSharedDigi* theDigi = (PndEmcSharedDigi*)fSharedDigiArray->At(iDigi);
                //      theDigi->Print();
                //      cout<<endl;
                //}
                //if(fSharedDigiArray->GetEntriesFast() != fDigiArray->GetEntriesFast())
                //      cout<<"fSharedDigiArray.size() - fDigiArray.size() ##="<<(fDigiArray->GetEntriesFast()-fSharedDigiArray->GetEntriesFast())<<endl;
        }
        //merge buffer bumps into current fBumpArray
        //Int_t nBumpTBD = fBumpArrayTBD->GetEntriesFast();
        //for(Int_t ibump = 0; ibump < nBumpTBD; ++ibump ){
        //      PndEmcBump* copy = (PndEmcBump*)fBumpArrayTBD->At(ibump);
        //      new ((*fBumpArray)[nBump++]) PndEmcBump(*copy);
        //}
        //fBumpArrayTBD->Delete();

        //variabes definition
        std::vector<PndEmcBump*> taggedBumpofToBeDeleted;
        PndEmcDigi* seedDigi(0);
        Double_t    EnergyOfSeedDigi(-1.), fTimeError(0.), CalibTimeOfaDigi(0.), CalibTimeOfSeedDigi(0.);
        Int_t       TheIndexOfModule;//TheIndexOfEnergy,  //[R.K. 01/2017] unused variable?
        Int_t       iDigi = PndEmcDigi::fDigiArrayTBD->GetEntriesFast();
        Int_t iBump = 0;
        Bool_t isBumpOK = kTRUE;
        for(iBump=0; iBump< fBumpArray->GetEntriesFast(); ++iBump)
        {
                EnergyOfSeedDigi = -1.;
                //isBumpOK = kTRUE;
                PndEmcBump* theBump = (PndEmcBump*)fBumpArray->At(iBump);
                //find the seed digi with maximum energy
                const std::vector<Int_t>& listOfDigi = theBump->DigiList();
                for(size_t id=0;id <listOfDigi.size();++id){
                        PndEmcDigi* theDigi = (PndEmcDigi*)fSharedDigiArray->At(listOfDigi[id]);
                        if(theDigi->GetEnergy() > EnergyOfSeedDigi){
                                seedDigi = theDigi;
                                EnergyOfSeedDigi = theDigi->GetEnergy();
                        }
                        if(fVerbose>1){
                                theDigi->Print();
                                CalibTimeOfaDigi = digiCalibrator.CalibrationEvtTimeByDigi(theDigi, kFALSE);
                                cout<<", CalibrationTimeOfEvent#"<<CalibTimeOfaDigi<<endl;
                        }
                }
                theBump->SetEventNo(seedDigi->fEvtNo);
                //check current bump belongs to current event by event time
                fTimeError          = digiCalibrator.GetTimeResolutionOfDigi(seedDigi);
                CalibTimeOfSeedDigi = digiCalibrator.CalibrationEvtTimeByDigi(seedDigi, kFALSE);
                if( CalibTimeOfSeedDigi - fEventTime < 5.*fTimeError){//out of region
                        isBumpOK = kTRUE;
                }else{
                        taggedBumpofToBeDeleted.push_back(theBump);
                        isBumpOK = kFALSE;
                }
                if(fVerbose>0){
                        std::cout<<"Bump["<<iBump<<"], (Nd, E, Es, CalT, EvtT, TimeR, EvtNo)=("
                                <<listOfDigi.size()
                                <<", "<<theBump->energy()
                                <<", "<<EnergyOfSeedDigi
                                <<", "<<CalibTimeOfSeedDigi
                                <<", "<<fEventTime
                                <<", "<<fTimeError
                                <<", "<<seedDigi->fEvtNo//theBump->where()
                                <<") ... ";
                        cout<<(isBumpOK ? "[ OK ]" : "[ Veto ]")<<endl;
                }


                //remove digis those time out of region compared to time of seed digi.
                if(isBumpOK)
                {
                        Double_t WeightedFactor = 0.;
                        Bool_t DigiStatus = kTRUE;
                        for(size_t id=0;id <listOfDigi.size();++id){
                                PndEmcDigi* theDigi = (PndEmcDigi*)fSharedDigiArray->At(listOfDigi[id]);
                                CalibTimeOfaDigi = digiCalibrator.CalibrationEvtTimeByDigi(theDigi, kFALSE);
                                fTimeError = digiCalibrator.GetTimeResolutionOfDigi(theDigi);
                                if(CalibTimeOfaDigi - CalibTimeOfSeedDigi  < 5.*fTimeError )
                                {//time  ok
                                        if(TMath::Abs(CalibTimeOfaDigi - CalibTimeOfSeedDigi) < 3.*fTimeError){
                                                WeightedFactor += 1./fTimeError/fTimeError;
                                        }
                                        DigiStatus = kTRUE;
                                }else{
                                        theBump->removeDigi(fSharedDigiArray, listOfDigi[id]);
                                        //save digis to a buffer for analysis in next event
                                        //if(CalibTimeOfaDigi - CalibTimeOfSeedDigi > 5.*fTimeError)
                                        new((*PndEmcDigi::fDigiArrayTBD)[iDigi++]) PndEmcDigi(*theDigi);
                                        DigiStatus = kFALSE;
                                }
                                if(fVerbose>0){
                                        TheIndexOfModule = theDigi->GetModule();
                                        cout<<"Digi["<<id<<"] of bump["<<iBump<<"], (E, Mod, CalT, EvtT, SigmaT, EvtNO) = ("
                                        <<theDigi->GetEnergy()
                                        <<", "<<TheIndexOfModule
                                        <<", "<<CalibTimeOfaDigi
                                        <<", "<<fEventTime
                                        <<", "<<fTimeError
                                        <<", "<<theDigi->fEvtNo
                                        <<(DigiStatus == kTRUE ? ") ... [ OK ]": ") ... [ Veto ]")<<endl;
                                //theDigi->Print();
                                //cout<<endl;
                                }
                        }
                }
        }
        //remove tagged bumps
        for(size_t ix=0;ix <taggedBumpofToBeDeleted.size();++ix){
                const std::vector<Int_t>& listOfDigi = taggedBumpofToBeDeleted[ix]->DigiList();
                //copy all digis of this bump to the buffer PndEmcDigi::fDigiArrayTBD;
                for(size_t id=0;id <listOfDigi.size();++id){
                        PndEmcDigi* theDigi = (PndEmcDigi*)fSharedDigiArray->At(listOfDigi[id]);
                        new((*PndEmcDigi::fDigiArrayTBD)[iDigi++]) PndEmcDigi(*theDigi);
                }
                fBumpArray->Remove(taggedBumpofToBeDeleted[ix]);
                fBumpArray->Compress();
        }

        if(fVerbose > 0){
                cout<<"******************************************************"<<endl;
                cout<<"After correction, bumps#"<<fBumpArray->GetEntriesFast()
                        <<", BumptoBeDetermined #"<<fBumpArrayTBD->GetEntriesFast()
                        <<", DigitoBeDetermined #"<<PndEmcDigi::fDigiArrayTBD->GetEntriesFast()<<endl;
                cout<<"**********************PndEmcCorrBump******************"<<endl;
        }
        //remake bump 
        for (Int_t i=0; i<fBumpArray->GetEntriesFast(); i++){
                PndEmcBump *tmpbump = (PndEmcBump*) fBumpArray->At(i);
                PndEmcClusterProperties clusterProperties(*tmpbump, fSharedDigiArray);
                HowManyDigi += tmpbump->NumberOfDigis();

                if (!tmpbump->IsEnergyValid())
                        tmpbump->SetEnergy(clusterProperties.Energy());
                if (!tmpbump->IsPositionValid())
                        tmpbump->SetPosition(clusterProperties.Where(fRecoPar->GetEmcClusterPosMethod(), fClusterPosParam));
                PndEmcXClMoments xClMoments(*tmpbump, fSharedDigiArray);
                tmpbump->SetZ20(xClMoments.AbsZernikeMoment(2, 0, 15));
                tmpbump->SetZ53(xClMoments.AbsZernikeMoment(5, 3, 15));
                tmpbump->SetLatMom(xClMoments.Lat());
        }
        HowManyCluster += fBumpArray->GetEntriesFast();

        fEventCounter++;
        if (fVerbose>2){
                timer.Stop();
                Double_t rtime = timer.RealTime();
                Double_t ctime = timer.CpuTime();
                cout << "PndEmcCorrBump, Real time " << rtime << " s, CPU time " << ctime << " s" << endl;
        }       
}

void PndEmcCorrBump::FinishTask

(

)

Definition at line 344 of file PndEmcCorrBump.cxx.

References HowManyCluster, and HowManyDigi.

{
        std::cout<<"*********************************************************"<<std::endl;
        std::cout<<"PndEmcCorrBump::FinishTask"<<std::endl;
        std::cout<<"*********************************************************"<<std::endl;
        std::cout<<"Read digis# "<<HowManyDigi<<std::endl;
        std::cout<<"Produce cluster# "<<HowManyCluster<<std::endl;
        std::cout<<"*********************************************************"<<std::endl;
}

InitStatus PndEmcCorrBump::Init ( )

virtual

Virtual method Init

Definition at line 87 of file PndEmcCorrBump.cxx.

References fBumpArray, fBumpArrayTBD, fClusterPosParam, fDigiArray, fRecoPar, fSharedDigiArray, PndEmcRecoPar::GetEmcClusterPosMethod(), PndEmcRecoPar::GetOffsetParmA(), PndEmcRecoPar::GetOffsetParmB(), PndEmcRecoPar::GetOffsetParmC(), HowManyCluster, HowManyDigi, and PndEmcDigi::InitDigiArrayTBD().

{
        PndEmcDigi::InitDigiArrayTBD();
        // Get RootManager
        FairRootManager* ioman = FairRootManager::Instance();
        if ( ! ioman )
        {
                cout << "-E- PndEmcCorrBump::Init: "
                << "RootManager not instantiated!" << endl;
                return kFATAL;
        }
        
        // Get input array
        fDigiArray = dynamic_cast<TClonesArray *>(ioman->GetObject("EmcDigi"));
        if ( ! fDigiArray ) {
                cout << "-W- PndEmcCorrBump::Init: "
                << "No PndEmcDigi array!" << endl;
                return kERROR;
        }
        // Get input array
        fSharedDigiArray =dynamic_cast<TClonesArray *> ( ioman->GetObject("EmcSharedDigi"));
        if ( ! fSharedDigiArray ) {
                cout << "-W- PndEmcCorrBump::Init: "
                << "No PndEmcSharedDigi array!" << endl;
                return kERROR;
        }

        // Get input array
        fBumpArray = dynamic_cast<TClonesArray *> (ioman->GetObject("EmcBump"));
        if ( ! fBumpArray ) {
                cout << "-W- PndEmcCorrBump::Init: "
                << "No EmcBump array!" << endl;
                return kERROR;
        }
        fBumpArrayTBD = new TClonesArray("PndEmcBump");

        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());
        }

        HowManyDigi = 0;
        HowManyCluster = 0;

        cout << "-I- PndEmcCorrBump: Intialization successfull" << endl;

        return kSUCCESS;
}

PndEmcCorrBump& PndEmcCorrBump::operator= ( const PndEmcCorrBump &  )

inlineprivate

Definition at line 88 of file PndEmcCorrBump.h.

{return *this;};

void PndEmcCorrBump::RunTimeBased ( )

inline

Definition at line 52 of file PndEmcCorrBump.h.

References fTimeOrderedDigi.

{fTimeOrderedDigi = kTRUE;}

void PndEmcCorrBump::SetParContainers ( )

privatevirtual

Reconstruction parameter container Get parameter containers

Definition at line 320 of file PndEmcCorrBump.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 PndEmcCorrBump::SetStorageOfData

(

Bool_t

val

)

Definition at line 339 of file PndEmcCorrBump.cxx.

References fStoreClusters, and val.

{
        fStoreClusters=val;
        return;
}

Member Data Documentation

PndEmcDigiCalibrator PndEmcCorrBump::digiCalibrator

private

Definition at line 102 of file PndEmcCorrBump.h.

Referenced by Exec().

TClonesArray* PndEmcCorrBump::fBumpArray

private

Definition at line 61 of file PndEmcCorrBump.h.

Referenced by Exec(), and Init().

TClonesArray* PndEmcCorrBump::fBumpArrayTBD

private

Definition at line 105 of file PndEmcCorrBump.h.

Referenced by Exec(), and Init().

std::vector<Double_t> PndEmcCorrBump::fClusterPosParam

private

Definition at line 107 of file PndEmcCorrBump.h.

Referenced by Exec(), and Init().

TClonesArray* PndEmcCorrBump::fDigiArray

private

Input array of CbmDigis

Definition at line 57 of file PndEmcCorrBump.h.

Referenced by Exec(), and Init().

PndEmcDigiPar* PndEmcCorrBump::fDigiPar

private

Geometry parameter container

Definition at line 73 of file PndEmcCorrBump.h.

Referenced by SetParContainers().

Int_t PndEmcCorrBump::fEventCounter =0

staticprivate

Definition at line 70 of file PndEmcCorrBump.h.

Referenced by Exec().

PndEmcGeoPar* PndEmcCorrBump::fGeoPar

private

Definition at line 72 of file PndEmcCorrBump.h.

Referenced by SetParContainers().

Int_t PndEmcCorrBump::fMapVersion

private

Output array of PndEmcClusters

Definition at line 68 of file PndEmcCorrBump.h.

PndEmcRecoPar* PndEmcCorrBump::fRecoPar

private

Digitisation parameter container

Definition at line 74 of file PndEmcCorrBump.h.

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

TClonesArray* PndEmcCorrBump::fSharedDigiArray

private

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

Definition at line 60 of file PndEmcCorrBump.h.

Referenced by Exec(), and Init().

Bool_t PndEmcCorrBump::fStoreClusters

private

Definition at line 81 of file PndEmcCorrBump.h.

Referenced by SetStorageOfData().

Bool_t PndEmcCorrBump::fTimeOrderedDigi

private

Definition at line 82 of file PndEmcCorrBump.h.

Referenced by RunTimeBased().

Int_t PndEmcCorrBump::fVerbose

private

Verbosity level

Definition at line 79 of file PndEmcCorrBump.h.

Referenced by Exec().

Int_t PndEmcCorrBump::HowManyCluster

private

Definition at line 92 of file PndEmcCorrBump.h.

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

Int_t PndEmcCorrBump::HowManyDigi

private

Definition at line 88 of file PndEmcCorrBump.h.

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

---

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

• PndEmcCorrBump.h
• PndEmcCorrBump.cxx