PndEmcCluster Class Reference

a cluster (group of neighboring crystals) of hit emc crystals More...

#include <PndEmcCluster.h>

Inheritance diagram for PndEmcCluster:

Classes
class	LinkScoreBoard

Public Member Functions
	PndEmcCluster ()
virtual void	Print (const Option_t *opt="") const
virtual	~PndEmcCluster ()
const std::vector< Int_t > &	DigiList () const
std::vector< Int_t > &	DigiList ()
const std::map< Int_t, Int_t > &	LocalMaxMap () const
const std::map< Int_t, Int_t > &	MemberDigiMap () const
virtual const PndEmcDigi *	Maxima (const TClonesArray *digiArray) const
virtual PndEmcDigi *	Maxima (const TClonesArray *digiArray)
virtual Double_t	energy () const
Double_t	theta () const
Double_t	phi () const
TVector3	position () const
TVector3	where () const
Double_t	x () const
Double_t	y () const
Double_t	z () const
Double_t	Z20 () const
Double_t	Z53 () const
Double_t	LatMom () const
Short_t	GetModule () const
bool	IsEnergyValid () const
bool	IsPositionValid () const
Int_t	NumberOfDigis () const
Int_t	NBumps () const
bool	isInCluster (PndEmcDigi *theDigi, const TClonesArray *digiArray)
virtual Double_t	DistanceToCentre (const TVector3 &aPoint) const
virtual Double_t	DistanceToCentre (const PndEmcDigi *aDigi) const
virtual void	addDigi (const TClonesArray *digiArray, Int_t iDigi)
virtual void	removeDigi (const TClonesArray *digiArray, Int_t iDigi)
virtual std::vector< Int_t > ::iterator	removeDigi (const TClonesArray *digiArray, std::vector< Int_t >::iterator it)
void	addCluster (PndEmcCluster *cluster, const TClonesArray *digiArray)
virtual void	addLocalMax (const TClonesArray *digiArray, Int_t iDigi)
virtual void	addLocalMax (const PndEmcDigi *digi)
virtual void	SetNBumps (unsigned nbumps)
void	SetEnergy (Double_t en)
void	SetPosition (TVector3 pos)
void	SetZ20 (Double_t z20)
void	SetZ53 (Double_t z53)
void	SetLatMom (Double_t latMom)
void	SetTrackEntering (const FairMultiLinkedData &tracks)
void	SetTrackExiting (const FairMultiLinkedData &tracks)
void	AddTracksEnteringExiting (const FairMultiLinkedData &tracksEntering, const FairMultiLinkedData &tracksExiting)
	Updates the links to entering and exiting tracks. More...
Double_t	GetEnergyCorrected () const
const std::vector< Int_t > &	GetMcList () const
Int_t	GetMcSize ()
Int_t	GetMcIndex (Int_t i=0)
FairMultiLinkedData	GetTrackEntering () const
FairMultiLinkedData	GetTrackExiting () const

Static Public Member Functions
static Double_t	FindPhiDiff (Double_t, Double_t)

Protected Member Functions
void	FillScoreBoard (FairMultiLinkedData tracks, std::map< FairLink, LinkScoreBoard > &scoreBoard, Int_t shift)

Protected Attributes
std::vector< Int_t >	fDigiList
std::vector< Int_t >	fMcList
std::map< Int_t, Int_t >	fMcMap
std::map< Int_t, Int_t >	fMemberDigiMap
std::map< Int_t, Int_t >	fLocalMaxMap
bool	fEnergyValid
Double_t	fEnergy
bool	fWhereValid
TVector3	fWhere
unsigned	fNbumps
Double_t	fZ20
Double_t	fZ53
Double_t	fLatMom
FairMultiLinkedData	fTrackEntering
FairMultiLinkedData	fTrackExiting

Private Member Functions
void	invalidateCache (bool)

Friends
class	PndEmcMakeCluster

Detailed Description

a cluster (group of neighboring crystals) of hit emc crystals

Definition at line 29 of file PndEmcCluster.h.

Constructor & Destructor Documentation

PndEmcCluster::PndEmcCluster

(

)

Definition at line 48 of file PndEmcCluster.cxx.

References fDigiList, fLocalMaxMap, fMcList, and fMemberDigiMap.

                              : 
        fDigiList(),
        fMcList(),
        fMemberDigiMap(),
        fLocalMaxMap(),
        fEnergyValid( false ),
        fEnergy( 0 ),
        fWhereValid( false ),
        fWhere( TVector3(0,0,0) ),
        fNbumps(0),
        fZ20(0),
        fZ53(0),
        fLatMom(0),
        fTrackEntering(),
        fTrackExiting()
{
        fDigiList.clear();
        fMcList.clear();
        fMemberDigiMap.clear();
        fLocalMaxMap.clear();
}

PndEmcCluster::~PndEmcCluster ( )

virtual

Definition at line 75 of file PndEmcCluster.cxx.

{}

Member Function Documentation

void PndEmcCluster::addCluster	(	PndEmcCluster *	cluster,
		const TClonesArray *	digiArray
	)

Definition at line 206 of file PndEmcCluster.cxx.

References addDigi(), AddTracksEnteringExiting(), DigiList(), GetTrackEntering(), and GetTrackExiting().

Referenced by PndEmcMakeClusterOnline::Exec(), and PndEmcMakeCluster::Exec().

{
        const vector<Int_t> tmpList = cluster->DigiList();
        vector<Int_t>::const_iterator digi_iter;
        for (digi_iter=tmpList.begin();digi_iter!=tmpList.end();++digi_iter)
        {
                addDigi(digiArray, *digi_iter);
        }
        AddTracksEnteringExiting(cluster->GetTrackEntering(), cluster->GetTrackExiting());
        AddLinks(cluster->GetLinks()); 
}

void PndEmcCluster::addDigi ( const TClonesArray *  digiArray, Int_t  iDigi  )

virtual

Definition at line 157 of file PndEmcCluster.cxx.

References digi, fDigiList, fMcMap, fMemberDigiMap, PndEmcDigi::GetDetectorId(), PndEmcDigi::GetTrackId(), and invalidateCache().

Referenced by addCluster(), PndEmcMakeClusterOnline::Exec(), PndEmcMakeCluster::Exec(), and PndEmcExpClusterSplitter::Exec().

{
        fDigiList.push_back(iDigi);
        PndEmcDigi *digi= (PndEmcDigi *) digiArray->At(iDigi);
        Int_t detectorId =digi->GetDetectorId(); 
        fMemberDigiMap.insert(std::pair<Int_t,Int_t>(detectorId, iDigi));

        fTimeStamp = (digi->GetTimeStamp() > this->GetTimeStamp()) ? digi->GetTimeStamp() : this->GetTimeStamp();
        //std::cout<<"digi belongs to track #"<<digi->GetTrackId()<<std::endl;
        ++fMcMap[digi->GetTrackId()];
        invalidateCache(kFALSE);
        SetInsertHistory(kFALSE);
        if(FairRunAna::Instance()->IsTimeStamp()) {
                AddLink((static_cast<PndEmcDigi*>(digiArray->At(iDigi))->GetEntryNr()));
        } else {
                AddLink(FairLink(-1, FairRootManager::Instance()->GetEntryNr(), "EmcDigi", iDigi));
        }
        SetInsertHistory(kTRUE);
}

void PndEmcCluster::addLocalMax ( const TClonesArray *  digiArray, Int_t  iDigi  )

virtual

Definition at line 218 of file PndEmcCluster.cxx.

References digi, fLocalMaxMap, and PndEmcDigi::GetDetectorId().

Referenced by PndEmc2DLocMaxFinder::Exec().

{
        PndEmcDigi *digi= (PndEmcDigi *) digiArray->At(iDigi);
        Int_t detectorId =digi->GetDetectorId(); 
        fLocalMaxMap.insert(std::pair<Int_t,Int_t>(detectorId, iDigi));
};

void PndEmcCluster::addLocalMax ( const PndEmcDigi *  digi )

virtual

Definition at line 225 of file PndEmcCluster.cxx.

References fLocalMaxMap, fMemberDigiMap, and PndEmcDigi::GetDetectorId().

{
        Int_t detectorId =digi->GetDetectorId();
        Int_t iDigi= fMemberDigiMap.find(detectorId)->second; 
        fLocalMaxMap.insert(std::pair<Int_t,Int_t>(detectorId, iDigi));
};

void PndEmcCluster::AddTracksEnteringExiting	(	const FairMultiLinkedData &	tracksEntering,
		const FairMultiLinkedData &	tracksExiting
	)

Updates the links to entering and exiting tracks.

For each track a bitmask is created which contains the information if the track enters or exits the cluster (bit 1 and 0 for tracks already in the cluster, bit 3 and 2 for the ones added here). Afterwards, the bitmask is checked for each track and the lists of entering and exiting tracks is updated. Tracks which entered before and exit now or vice versa (bitmask=6 or 9) are removed from the list as they only cross between crystals in the same cluster.

Parameters

tracksEntering	Additional tracks entering the cluster (e.g. from adding digi or merging clusters)
tracksExiting	Additional tracks exiting the cluster

Returns

void

Definition at line 535 of file PndEmcCluster.cxx.

References FillScoreBoard(), fTrackEntering, and fTrackExiting.

Referenced by addCluster(), and PndEmcMakeCluster::Exec().

{
        std::map<FairLink, LinkScoreBoard> scoreBoard;
        std::set<FairLink> entering, exiting;

//      std::cout << "tracksEntering " << tracksEntering << std::endl;
//      std::cout << "tracksExiting " << tracksExiting << std::endl;
//      std::cout << "existingEntering " << fTrackEntering << std::endl;
//      std::cout << "existingExiting " << fTrackExiting << std::endl;

        FillScoreBoard(tracksEntering, scoreBoard, 3);
        FillScoreBoard(tracksExiting, scoreBoard, 2);
        FillScoreBoard(fTrackEntering, scoreBoard, 1);
        FillScoreBoard(fTrackExiting, scoreBoard, 0);

        for (std::map<FairLink, LinkScoreBoard>::iterator iter = scoreBoard.begin(); iter != scoreBoard.end(); iter++){
                //std::cout << iter->first << " " << iter->second.score << std::endl;
                switch (iter->second.score){
                        case 15: entering.insert(iter->first); exiting.insert(iter->first); break;
                        case 14: entering.insert(iter->first); break;
                        case 13: exiting.insert(iter->first); break;
                        case 12: entering.insert(iter->first); exiting.insert(iter->first); break;
                        case 11: entering.insert(iter->first); break;
                        case 10: std::cout << "-E- PndEmcCluster::AddTracksEnteringExiting Same particle entering twice!" << std::endl; break;
                        case 9: break;
                        case 8: entering.insert(iter->first); break;
                        case 7: exiting.insert(iter->first); break;
                        case 6: break;
                        case 5: std::cout << "-E- PndEmcCluster::AddTracksEnteringExiting Same particle exiting twice!" << std::endl; break;
                        case 4: exiting.insert(iter->first); break;
                        case 3: entering.insert(iter->first); exiting.insert(iter->first); break;
                        case 2: entering.insert(iter->first);  break;
                        case 1: exiting.insert(iter->first); break;
                        case 0: break;
                        default: std::cout << "-E- PndEmcCluster::AddTracksEnteringExiting wrong score " << iter->second.score << std::endl; break;
                }
        }
        fTrackEntering.SetLinks(entering);
        fTrackExiting.SetLinks(exiting);

//      std::cout << "Entering after merge " << fTrackEntering << std::endl;
//      std::cout << "Exiting after merge " << fTrackExiting << std::endl;
}

const std::vector<Int_t>& PndEmcCluster::DigiList ( ) const

inline

Definition at line 40 of file PndEmcCluster.h.

References fDigiList.

Referenced by addCluster(), PndEmcClusterEnergySums::E25(), PndEmcClusterEnergySums::E9(), PndEmcClusterProperties::Energy(), PndEmcClusterEnergySums::energy(), PndEmcClusterRemoveDuplCrys::Exec(), PndEmcCorrBump::Exec(), PndEmcAnalysis::Exec(), PndEmcExpClusterSplitter::Exec(), PndEmcPhiBumpSplitter::Exec(), PndEmcMakeClusterOnline::FinishCluster(), for(), PndEmcClusterProperties::GravWhere(), PndEmcXClMoments::Init(), PndEmcClusterProperties::LiloWhere(), PndEmcClusterProperties::LinearWhere(), PndEmcClusterProperties::Major_axis(), PndEmcClusterProperties::Mass(), PndEmcBump::Print(), and PndEmcXClMoments::Print().

{return fDigiList;};

std::vector<Int_t>& PndEmcCluster::DigiList ( )

inline

Definition at line 41 of file PndEmcCluster.h.

References fDigiList.

{return fDigiList;};

Double_t PndEmcCluster::DistanceToCentre ( const TVector3 &  aPoint ) const

virtual

Definition at line 367 of file PndEmcCluster.cxx.

References where().

{
        return  ( where() - aPoint ).Mag();
}

Double_t PndEmcCluster::DistanceToCentre ( const PndEmcDigi *  aDigi ) const

virtual

Definition at line 373 of file PndEmcCluster.cxx.

References where(), and PndEmcDigi::where().

{
        return  ( where() - aDigi->where() ).Mag();
}

Double_t PndEmcCluster::energy ( ) const

virtual

Definition at line 79 of file PndEmcCluster.cxx.

References fEnergy, and fEnergyValid.

Referenced by analysis_cluster_energyCorrection_Ntuple(), PndEmcXClMoments::Calc_AbsZernikeMoment(), PndEmcDigiCalibrator::CalibrationEvtTimeByBump(), PndEmcClusterMoments::Centre1(), PndEmcClusterMoments::Centre2(), PndPidCorrelator::ConstructNeutralCandidate(), emc_correction_hist(), emc_correction_parametrization(), emc_correction_QA(), PndEmcClusterHistCalibrator::Energy(), PndEmcClusterSimpleCalibrator::Energy(), error_matrix_fit(), PndKFParticleFinder::Exec(), PndEmcMakeRecoHit::Exec(), PndEmcHdrFiller::Exec(), PndEmcCorrBump::Exec(), PndEmcMakeCorr::Exec(), PndEmcAnalysis::Exec(), PndEmcXClMoments::Fast_AbsZernikeMoment(), fit_resolution(), for(), PndEmcErrorMatrix::Get4MomentumErrorMatrix(), PndPidCorrelator::GetEmcInfo(), GetEnergyCorrected(), GetEntriesFast(), PndEmcErrorMatrix::GetErrorMatrix(), PndEmcErrorMatrix::GetErrorP7(), PndPidCorrelator::GetFscInfo(), PndEmcDigiCalibrator::GetTimeResolutionOfShower(), PndEmcClusterDraw::GetValue(), PndEmcClusterProperties::GravWhere(), PndEmcClusterMoments::Major1(), PndEmcClusterMoments::Major2(), PndEmcClusterMoments::Minor1(), PndEmcClusterMoments::Minor2(), PndEmcClusterMoments::Phi1(), PndEmcClusterMoments::Phi2(), PndEmcClusterDistances::PndEmcClusterDistances(), Print(), PndEmcBump::Print(), QAmacro_emc_3(), reco_analys(), PndEmcClusterMoments::SecondMoment(), PndEmcXClMoments::SecondMoment(), PndEmcXClMoments::SecondMomentPhi(), PndEmcXClMoments::SecondMomentTheta(), PndEmcClusterMoments::SecondMomentTP(), PndEmcClusterMoments::Theta1(), PndEmcClusterMoments::Theta2(), vis_mc_ids(), and PndEmcClusterHistCalibrator::Where().

{
        if (fEnergyValid)
                return fEnergy;
        else
        {
                std::cout<<"Energy of cluster is not defined"<<std::endl;
                abort();
        }
}

void PndEmcCluster::FillScoreBoard ( FairMultiLinkedData  tracks, std::map< FairLink, LinkScoreBoard > &  scoreBoard, Int_t  shift  )

protected

Definition at line 580 of file PndEmcCluster.cxx.

Referenced by AddTracksEnteringExiting().

{
        std::set<FairLink> links = tracks.GetLinks();
        for (std::set<FairLink>::iterator linkIter = links.begin(); linkIter != links.end(); linkIter++){
                scoreBoard[*linkIter].SetValShift(kTRUE, shift);
        }
}

Double_t PndEmcCluster::FindPhiDiff ( Double_t  phi1, Double_t  phi2  )

static

Definition at line 141 of file PndEmcCluster.cxx.

References Double_t, and Pi.

Referenced by PndEmcClusterMoments::Major1(), PndEmcClusterMoments::Major2(), PndEmcClusterProperties::Major_axis(), PndEmcClusterMoments::Minor1(), PndEmcClusterMoments::Minor2(), PndEmcClusterMoments::Phi1(), and PndEmcClusterMoments::Phi2().

{
        Double_t diff;
        diff = phi1 - phi2;

        while( diff>  TMath::Pi()  ) diff -= 2*TMath::Pi();
        while( diff< -TMath::Pi()  ) diff += 2*TMath::Pi();

        return diff;
}

Double_t PndEmcCluster::GetEnergyCorrected

(

)

const

Definition at line 379 of file PndEmcCluster.cxx.

References b1, b2, cos(), Double_t, energy(), exp(), log(), p1, p2, Pi, sqrt(), t0, t1, t2, t3, and where().

Referenced by analysis_digi_cluster_7gammaAnalysis_FullEmc(), analysis_digi_cluster_allGammasAnalysis_FullEmc(), analysis_digi_cluster_allGammasAnalysis_fwendcap(), analysis_digi_cluster_fwendcap(), emc_correction_QA(), for(), PndPidBremCorrector::GetSepPhotonE(), and vis_mc_ids().

{
        Double_t e=this->energy();
        TVector3 clusterPosition= this->where();
        Double_t theta_cluster=clusterPosition.Theta();

        Double_t e1=e;
        Double_t theta1=theta_cluster;

        if ( (clusterPosition.Z() < 180.0)&&(theta_cluster<140.*TMath::Pi()/180.))
        {
                if (e<0.03) e1 = 0.03;
                if (e>8.0)  e1 = 8.0 ;
        }

        if ( (clusterPosition.Z() < 180.0)&&(theta_cluster>140.*TMath::Pi()/180.))
        {
                if (e<0.03) e1 = 0.03;
                if (e>2.0)  e1 = 2.0 ;
        }

        if (clusterPosition.Z() > 180.0)
        {
                if (e<0.01) e1 = 0.01;
                if (e>16.0)  e1 = 16.0 ;
        }

        double b0 = 1.45312;
        double b1 = 2.79086e-02;
        double b2 = 3.91932e-04;
        double b3 =-1.23117e-03;
        double b4 = 2.72270e-01;
        double b5 =-1.31540;
        double b6 = 1.44447;
        double b7 =-4.05724e-01;
        double b8 =-2.07396;
        double b9 = 4.80507e-02;

        double p0=0;
        double p1=0;
        double p2=0;
        double p3=0;
        double p4=0;
        double p5=0;
        double p6=0;
        double p7=0;
        double p8=0;
        double p9=0;

        if(e1<1.0)
        {
                p0 =  4.13189e-02;   
                p1 = -2.03834e-02;  
                p2 = -2.58086e-03;   
                p3 = -1.77821e-03;   
                p4 = -1.73738e-02;   
                p5 =  7.40362e-02;   
                p6 = -6.41892e-02;   
                p7 = -9.85564e-02;   
                p8 =  1.50123e-01;   
                p9 = -7.87742e-04;   
        }
        else
        {
                p0 =  5.05003e-02;  
                p1 = -3.47672e-02;   
                p2 =  3.72767e-02;   
                p3 = -1.26492e-02;   
                p4 = -2.16876e-02;   
                p5 =  1.02682e-01;   
                p6 = -9.85242e-02;   
                p7 = -1.39872e-01;   
                p8 =  2.02309e-01;   
                p9 =  1.11696e-03;   
        }


        double t0 = 1.81631;
        double t1 =-1.71202e-02;
        double t2 = 3.59161e-03;
        double t3 =-3.46712e-04;
        double t4 =-3.73691e-01;
        double t5 =-1.56688;
        double t6 =-1.62618;
        double t7 =-4.10972e-01;
        double t8 = 2.2222;
        double t9 = 4.60908e-03;


        double factor1= p0
                +p1*log(e1)
                +p2*log(e1)*log(e1)
                +p3*log(e1)*log(e1)*log(e1)
                +p4*cos(theta1)
                +p5*cos(theta1)*cos(theta1)
                +p6*cos(theta1)*cos(theta1)*cos(theta1)
                +p7*cos(theta1)*cos(theta1)*cos(theta1)*cos(theta1)
                +p8*cos(theta1)*cos(theta1)*cos(theta1)*cos(theta1)*cos(theta1)
                +p9*log(e1)*cos(theta1);

        double factor2= t0
                +t1*log(e1)
                +t2*log(e1)*log(e1)
                +t3*log(e1)*log(e1)*log(e1)
                +t4*cos(theta1)
                +t5*cos(theta1)*cos(theta1)
                +t6*cos(theta1)*cos(theta1)*cos(theta1)
                +t7*cos(theta1)*cos(theta1)*cos(theta1)*cos(theta1)
                +t8*cos(theta1)*cos(theta1)*cos(theta1)*cos(theta1)*cos(theta1)
                +t9*log(e1)*cos(theta1);


        double factor3= b0
                +b1*log(e1)
                +b2*log(e1)*log(e1)
                +b3*log(e1)*log(e1)*log(e1)
                +b4*cos(theta1)
                +b5*cos(theta1)*cos(theta1)
                +b6*cos(theta1)*cos(theta1)*cos(theta1)
                +b7*cos(theta1)*cos(theta1)*cos(theta1)*cos(theta1)
                +b8*cos(theta1)*cos(theta1)*cos(theta1)*cos(theta1)*cos(theta1)
                +b9*log(e1)*cos(theta1);


        double eout1=e* exp(factor1); 
        double eout2=e* exp(factor2); 
        double eout3=e* exp(factor3); 

        double eout4=(3.31694-0.0183379/sqrt(e1)+0.0327113/e1+0.00040156/(e1*e1)-0.00641305/(e1*sqrt(e1)))*e/3.0144;



        if ( clusterPosition.Z() > 500.0)
                return eout4;
        else if ( (clusterPosition.Z() < 180.0)&&(theta_cluster>140.*TMath::Pi()/180.))
                return eout3;
        else  if ( (clusterPosition.Z() < 180.0)&&(theta_cluster<140.*TMath::Pi()/180.))
                return eout1;
        else 
                return eout2;

}

Int_t PndEmcCluster::GetMcIndex ( Int_t  i = 0 )

inline

Definition at line 116 of file PndEmcCluster.h.

References fMcList, and i.

{ return fMcList[i]; }

const std::vector< Int_t > & PndEmcCluster::GetMcList

(

)

const

Definition at line 314 of file PndEmcCluster.cxx.

References a, b, fMcList, and i.

Referenced by PndPidCorrelator::ConstructNeutralCandidate(), PndPidCorrelator::GetEmcInfo(), PndPidCorrelator::GetFscInfo(), and vis_mc_ids().

                                                       {
        std::vector< map_ele > sortedVec;
        std::vector<FairLink> mcLinks;
        FairMultiLinkedData mcFairLinks = GetLinksWithType(FairRootManager::Instance()->GetBranchId("MCTrack"));

        for (int i = 0; i < mcFairLinks.GetNLinks(); i++){
                mcLinks.push_back(mcFairLinks.GetLink(i));
        }

        std::sort(mcLinks.begin(), mcLinks.end(), [](const FairLink& a, const FairLink& b) -> bool {
                return a.GetIndex() < b.GetIndex();
        });

        fMcList.clear();
        for (auto link : mcLinks)
                fMcList.push_back(link.GetIndex());
        return fMcList;
}

Int_t PndEmcCluster::GetMcSize ( )

inline

Definition at line 115 of file PndEmcCluster.h.

References fMcList.

{ return fMcList.size(); }

Short_t PndEmcCluster::GetModule

(

)

const

Definition at line 282 of file PndEmcCluster.cxx.

References fLocalMaxMap, and fMemberDigiMap.

Referenced by analysis_cluster_energyCorrection_Ntuple(), analysis_digi_cluster_7gammaAnalysis_FullEmc(), analysis_digi_cluster_allGammasAnalysis_FullEmc(), analysis_digi_cluster_allGammasAnalysis_fwendcap(), analysis_digi_cluster_fwendcap(), PndEmcDigiCalibrator::CalibrationEvtTimeByBump(), PndPidCorrelator::ConstructNeutralCandidate(), emc_correction_hist(), emc_correction_parametrization(), PndEmcClusterHistCalibrator::Energy(), PndEmcAnalysis::Exec(), PndPidCorrelator::Exec(), PndPidCorrelator::GetEmcInfo(), PndEmcErrorMatrix::GetErrorMatrix(), PndPidCorrelator::GetFscInfo(), PndEmcDigiCalibrator::GetTimeResolutionOfShower(), and PndEmcClusterHistCalibrator::Where().

{
        // Get module in which cluster is located
        // If fLocalMaxMap is defined, i.e. after bump splitting procedure
        // the number of module if taken from the first local maxima.
        // Otherwise it is taken from the first digi in cluster.
        Short_t module;
        Int_t detectorId;
        if (fLocalMaxMap.size()>0)
        {
                std::map<Int_t,Int_t>::const_iterator iter=fLocalMaxMap.begin();
                detectorId=iter->first;
        }
        else
        {
                std::map<Int_t,Int_t>::const_iterator iter=fMemberDigiMap.begin();
                detectorId=iter->first;
        }       

        module=detectorId/100000000;

        return module;
};

FairMultiLinkedData PndEmcCluster::GetTrackEntering ( ) const

inline

Definition at line 117 of file PndEmcCluster.h.

References fTrackEntering.

Referenced by addCluster(), PndEmcExpClusterSplitter::Exec(), and PndEmcMakeCluster::FinishCluster().

{ return fTrackEntering;}

FairMultiLinkedData PndEmcCluster::GetTrackExiting ( ) const

inline

Definition at line 118 of file PndEmcCluster.h.

References fTrackExiting.

Referenced by addCluster().

{ return fTrackExiting;}

void PndEmcCluster::invalidateCache ( bool  kState )

private

Definition at line 334 of file PndEmcCluster.cxx.

References fEnergyValid, and fWhereValid.

Referenced by addDigi().

{
        fEnergyValid = kState;
        fWhereValid = kState;
}

bool PndEmcCluster::IsEnergyValid ( ) const

inline

Definition at line 77 of file PndEmcCluster.h.

References fEnergyValid.

Referenced by PndEmcCorrBump::Exec(), and PndEmcExpClusterSplitter::Exec().

{return fEnergyValid;};

bool PndEmcCluster::isInCluster	(	PndEmcDigi *	theDigi,
		const TClonesArray *	digiArray
	)

Definition at line 234 of file PndEmcCluster.cxx.

References fDigiList, and PndEmcDigi::isNeighbour().

Referenced by PndEmcMakeClusterOnline::Exec(), and PndEmcMakeCluster::Exec().

{
        vector<Int_t>::iterator digi_iter;
        for (digi_iter=fDigiList.begin();digi_iter!=fDigiList.end();++digi_iter)
        {
                if(theDigi->isNeighbour((PndEmcDigi *) digiArray->At(*digi_iter))) return true;
        }

        return false;
}

bool PndEmcCluster::IsPositionValid ( ) const

inline

Definition at line 78 of file PndEmcCluster.h.

References fWhereValid.

Referenced by PndEmcCorrBump::Exec(), and PndEmcExpClusterSplitter::Exec().

{return fWhereValid;};

Double_t PndEmcCluster::LatMom ( ) const

inline

Definition at line 71 of file PndEmcCluster.h.

References fLatMom.

Referenced by PndPidCorrelator::ConstructNeutralCandidate(), PndEmcAnalysis::Exec(), PndPidCorrelator::GetEmcInfo(), and PndPidCorrelator::GetFscInfo().

{return fLatMom;}

const std::map<Int_t,Int_t>& PndEmcCluster::LocalMaxMap ( ) const

inline

Definition at line 42 of file PndEmcCluster.h.

References fLocalMaxMap.

Referenced by PndEmc2DLocMaxFinder::Exec(), PndEmcAnalysis::Exec(), and PndEmcExpClusterSplitter::Exec().

{return fLocalMaxMap;};

const PndEmcDigi * PndEmcCluster::Maxima ( const TClonesArray *  digiArray ) const

virtual

Definition at line 247 of file PndEmcCluster.cxx.

References digi, Double_t, fDigiList, PndEmcDigi::GetEnergy(), and max().

Referenced by PndEmcClusterEnergySums::E1(), PndEmcClusterEnergySums::E25(), PndEmcClusterEnergySums::E9(), PndEmc2DLocMaxFinder::Exec(), PndEmcMakeCluster::FinishCluster(), and PndEmcXClMoments::Init().

{
        Double_t max=0;
        PndEmcDigi *biggest=0;

        std::vector<Int_t>::const_iterator digipos;
        for (digipos=fDigiList.begin();digipos!=fDigiList.end();++digipos){
                PndEmcDigi *digi = (PndEmcDigi *) digiArray->At(*digipos);
                if ( max < digi->GetEnergy() ) {
                        max=digi->GetEnergy();
                        biggest=digi;
                }
        }

        return( biggest );
}

PndEmcDigi * PndEmcCluster::Maxima ( const TClonesArray *  digiArray )

virtual

Definition at line 265 of file PndEmcCluster.cxx.

References digi, Double_t, fDigiList, PndEmcDigi::GetEnergy(), and max().

{
        Double_t max=0;
        PndEmcDigi *biggest=0;

        std::vector<Int_t>::iterator digipos;
        for (digipos=fDigiList.begin();digipos!=fDigiList.end();++digipos){
                PndEmcDigi *digi = (PndEmcDigi *) digiArray->At(*digipos);
                if ( max < digi->GetEnergy() ) {
                        max=digi->GetEnergy();
                        biggest=digi;
                }
        }

        return( biggest );
}

const std::map<Int_t,Int_t>& PndEmcCluster::MemberDigiMap ( ) const

inline

Definition at line 43 of file PndEmcCluster.h.

References fMemberDigiMap.

Referenced by PndEmcClusterRemoveDuplCrys::Exec(), PndEmc2DLocMaxFinder::Exec(), PndEmcMakeCorr::Exec(), PndEmcAnalysis::Exec(), PndEmcExpClusterSplitter::Exec(), and PndEmc2DLocMaxFinder::isALocalMax().

{return fMemberDigiMap;};

Int_t PndEmcCluster::NBumps

(

)

const

Definition at line 346 of file PndEmcCluster.cxx.

References fNbumps.

Referenced by PndPidCorrelator::ConstructNeutralCandidate(), PndEmcAnalysis::Exec(), and PndPidCorrelator::GetEmcInfo().

{
        return fNbumps;
}

Int_t PndEmcCluster::NumberOfDigis

(

)

const

Definition at line 341 of file PndEmcCluster.cxx.

References fDigiList.

Referenced by analysis_digi_cluster_7gammaAnalysis_FullEmc(), analysis_digi_cluster_allGammasAnalysis_FullEmc(), analysis_digi_cluster_allGammasAnalysis_fwendcap(), analysis_digi_cluster_fwendcap(), PndEmcDigiCalibrator::CalibrationEvtTimeByBump(), PndPidCorrelator::ConstructNeutralCandidate(), PndEmcCorrBump::Exec(), PndEmcAnalysis::Exec(), PndEmcExpClusterSplitter::Exec(), fit_resolution(), for(), PndPidCorrelator::GetEmcInfo(), PndPidCorrelator::GetFscInfo(), and QAmacro_emc_3().

{
        return fDigiList.size();
}

Double_t PndEmcCluster::phi

(

)

const

Definition at line 98 of file PndEmcCluster.cxx.

References where().

Referenced by analysis_cluster_energyCorrection_Ntuple(), PndEmcHdrFiller::Exec(), PndEmcErrorMatrix::GetErrorMatrix(), PndEmcClusterMoments::Major1(), PndEmcClusterMoments::Major2(), PndEmcClusterProperties::Major_axis(), PndEmcClusterMoments::Minor1(), PndEmcClusterMoments::Minor2(), PndEmcClusterMoments::Phi1(), PndEmcClusterMoments::Phi2(), and PndEmcClusterDistances::PndEmcClusterDistances().

{
        return where().Phi();
}

TVector3 PndEmcCluster::position

(

)

const

Definition at line 104 of file PndEmcCluster.cxx.

References where().

Referenced by PndEmcMakeRecoHit::Exec(), PndEmcPhiBumpSplitter::Exec(), and PndPidBremCorrector::GetSepPhotonE().

{
        return this->where();
}

void PndEmcCluster::Print ( const Option_t *  opt = "" ) const

virtual

Definition at line 357 of file PndEmcCluster.cxx.

References energy(), fTrackEntering, and fTrackExiting.

{
        std::cout<<"*********************************"<< endl;
        std::cout<<"total energy of cluster: "<< energy() << endl;
        std::cout <<"TrackEntering: " << fTrackEntering << std::endl;
        std::cout <<"TrackExiting: " << fTrackExiting << std::endl;
}

void PndEmcCluster::removeDigi ( const TClonesArray *  digiArray, Int_t  iDigi  )

virtual

Definition at line 177 of file PndEmcCluster.cxx.

References digi, fDigiList, fMcMap, fMemberDigiMap, PndEmcDigi::GetDetectorId(), PndEmcDigi::GetTrackId(), and i.

Referenced by PndEmcClusterRemoveDuplCrys::Exec(), and PndEmcCorrBump::Exec().

{
        for(size_t i=0;i<fDigiList.size();++i){
                if(fDigiList[i] == iDigi){
                        fDigiList.erase(fDigiList.begin()+i);
                        break;
                }
        }
        PndEmcDigi *digi= (PndEmcDigi *) digiArray->At(iDigi);
        Int_t detectorId =digi->GetDetectorId(); 
        fMemberDigiMap.erase(detectorId);
        --fMcMap[digi->GetTrackId()];
}

std::vector< Int_t >::iterator PndEmcCluster::removeDigi ( const TClonesArray *  digiArray, std::vector< Int_t >::iterator  it  )

virtual

Definition at line 191 of file PndEmcCluster.cxx.

References digi, fDigiList, fMcMap, fMemberDigiMap, PndEmcDigi::GetDetectorId(), and PndEmcDigi::GetTrackId().

{
        PndEmcDigi *digi= (PndEmcDigi *) digiArray->At(*iDigi);
        //std::cout << "removing digi (ptr|detId): (" << digi << " | " << digi->GetDetectorId() << ")" << std::endl;

        std::map<Int_t, Int_t>::iterator it = fMemberDigiMap.find(digi->GetDetectorId());
        if(it != fMemberDigiMap.end() && it->second == *iDigi) {
                fMemberDigiMap.erase(it);
        }

        --fMcMap[digi->GetTrackId()];
        return fDigiList.erase(iDigi);
}

void PndEmcCluster::SetEnergy ( Double_t  en )

inline

Definition at line 100 of file PndEmcCluster.h.

References fEnergy, and fEnergyValid.

Referenced by PndEmcCorrBump::Exec(), PndEmcExpClusterSplitter::Exec(), PndEmcPhiBumpSplitter::Exec(), PndEmcMakeClusterOnline::FinishCluster(), and PndEmcMakeCluster::FinishCluster().

{fEnergy=en; fEnergyValid=true;}

void PndEmcCluster::SetLatMom ( Double_t  latMom )

inline

Definition at line 104 of file PndEmcCluster.h.

References fLatMom.

Referenced by PndEmcCorrBump::Exec(), PndEmcExpClusterSplitter::Exec(), and PndEmcMakeCluster::FinishCluster().

{fLatMom=latMom;}

void PndEmcCluster::SetNBumps ( unsigned  nbumps )

virtual

Definition at line 352 of file PndEmcCluster.cxx.

References fNbumps.

Referenced by PndEmcExpClusterSplitter::Exec(), and PndEmcBump::PndEmcBump().

                                        {
        fNbumps = nbumps;
}

void PndEmcCluster::SetPosition ( TVector3  pos )

inline

Definition at line 101 of file PndEmcCluster.h.

References fWhere, fWhereValid, and pos.

Referenced by PndEmcCorrBump::Exec(), PndEmcExpClusterSplitter::Exec(), PndEmcPhiBumpSplitter::Exec(), and PndEmcMakeCluster::FinishCluster().

{fWhere=pos; fWhereValid=true;}

void PndEmcCluster::SetTrackEntering ( const FairMultiLinkedData &  tracks )

inline

Definition at line 105 of file PndEmcCluster.h.

References fTrackEntering.

Referenced by PndEmcMakeCluster::Exec().

{fTrackEntering = tracks;}

void PndEmcCluster::SetTrackExiting ( const FairMultiLinkedData &  tracks )

inline

Definition at line 106 of file PndEmcCluster.h.

References fTrackExiting.

Referenced by PndEmcMakeCluster::Exec().

{fTrackExiting = tracks;}

void PndEmcCluster::SetZ20 ( Double_t  z20 )

inline

Definition at line 102 of file PndEmcCluster.h.

References fZ20.

Referenced by PndEmcCorrBump::Exec(), PndEmcExpClusterSplitter::Exec(), and PndEmcMakeCluster::FinishCluster().

{fZ20=z20;}

void PndEmcCluster::SetZ53 ( Double_t  z53 )

inline

Definition at line 103 of file PndEmcCluster.h.

References fZ53.

Referenced by PndEmcCorrBump::Exec(), PndEmcExpClusterSplitter::Exec(), and PndEmcMakeCluster::FinishCluster().

{fZ53=z53;}

Double_t PndEmcCluster::theta

(

)

const

Definition at line 92 of file PndEmcCluster.cxx.

References where().

Referenced by analysis_cluster_energyCorrection_Ntuple(), PndEmcHdrFiller::Exec(), PndEmcErrorMatrix::GetErrorMatrix(), PndEmcClusterMoments::Major1(), PndEmcClusterMoments::Major2(), PndEmcClusterProperties::Major_axis(), PndEmcClusterMoments::Minor1(), PndEmcClusterMoments::Minor2(), and PndEmcClusterDistances::PndEmcClusterDistances().

{
        return where().Theta();
}

TVector3 PndEmcCluster::where

(

)

const

Definition at line 111 of file PndEmcCluster.cxx.

References fWhere, and fWhereValid.

Referenced by analysis_digi_cluster_7gammaAnalysis_FullEmc(), analysis_digi_cluster_allGammasAnalysis_FullEmc(), analysis_digi_cluster_allGammasAnalysis_fwendcap(), analysis_digi_cluster_fwendcap(), PndEmcDigiCalibrator::CalibrationEvtTimeByBump(), PndPidCorrelator::ConstructNeutralCandidate(), DistanceToCentre(), emc_correction_hist(), emc_correction_parametrization(), emc_correction_QA(), PndEmcClusterHistCalibrator::Energy(), PndEmcClusterSimpleCalibrator::Energy(), error_matrix_fit(), PndEmcMakeCorr::Exec(), PndEmcAnalysis::Exec(), for(), PndEmcErrorMatrix::Get4MomentumErrorMatrix(), PndPidCorrelator::GetEmcInfo(), GetEnergyCorrected(), GetEntriesFast(), PndEmcErrorMatrix::GetErrorP7(), PndPidCorrelator::GetFscInfo(), PndEmcClusterDraw::GetVector(), phi(), PndEmcClusterDistances::PndEmcClusterDistances(), position(), QAmacro_emc_3(), reco_analys(), PndEmcClusterMoments::SecondMomentTP(), theta(), PndEmcClusterHistCalibrator::Where(), PndEmcClusterSimpleCalibrator::Where(), x(), y(), and z().

{
        if (fWhereValid)
                return fWhere;
        else
        {
                std::cout<<"Position of cluster is not defined"<<std::endl;
                abort();
        }
}

Double_t PndEmcCluster::x

(

)

const

Definition at line 123 of file PndEmcCluster.cxx.

References where().

Referenced by analysis_digi_cluster_7gammaAnalysis_FullEmc(), analysis_digi_cluster_allGammasAnalysis_FullEmc(), analysis_digi_cluster_allGammasAnalysis_fwendcap(), analysis_digi_cluster_fwendcap(), PndKFParticleFinderQA::Exec(), and vis_mc_ids().

{
        return where().x();
}

Double_t PndEmcCluster::y

(

)

const

Definition at line 129 of file PndEmcCluster.cxx.

References where().

Referenced by analysis_digi_cluster_7gammaAnalysis_FullEmc(), analysis_digi_cluster_allGammasAnalysis_FullEmc(), analysis_digi_cluster_allGammasAnalysis_fwendcap(), analysis_digi_cluster_fwendcap(), PndKFParticleFinderQA::Exec(), and vis_mc_ids().

{
        return where().y();
}

Double_t PndEmcCluster::z

(

)

const

Definition at line 135 of file PndEmcCluster.cxx.

References where().

Referenced by analysis_digi_cluster_7gammaAnalysis_FullEmc(), analysis_digi_cluster_allGammasAnalysis_FullEmc(), analysis_digi_cluster_allGammasAnalysis_fwendcap(), analysis_digi_cluster_fwendcap(), PndKFParticleFinderQA::Exec(), and for().

{
        return where().z();
}

Double_t PndEmcCluster::Z20 ( ) const

inline

Definition at line 67 of file PndEmcCluster.h.

References fZ20.

Referenced by PndPidCorrelator::ConstructNeutralCandidate(), PndEmcAnalysis::Exec(), PndPidCorrelator::GetEmcInfo(), and PndPidCorrelator::GetFscInfo().

{return fZ20;}

Double_t PndEmcCluster::Z53 ( ) const

inline

Definition at line 69 of file PndEmcCluster.h.

References fZ53.

Referenced by PndPidCorrelator::ConstructNeutralCandidate(), PndEmcAnalysis::Exec(), PndPidCorrelator::GetEmcInfo(), and PndPidCorrelator::GetFscInfo().

{return fZ53;}

Friends And Related Function Documentation

friend class PndEmcMakeCluster

friend

Definition at line 122 of file PndEmcCluster.h.

Member Data Documentation

std::vector<Int_t> PndEmcCluster::fDigiList

protected

Definition at line 140 of file PndEmcCluster.h.

Referenced by addDigi(), DigiList(), isInCluster(), Maxima(), NumberOfDigis(), PndEmcCluster(), and removeDigi().

Double_t PndEmcCluster::fEnergy

protected

Definition at line 147 of file PndEmcCluster.h.

Referenced by energy(), and SetEnergy().

bool PndEmcCluster::fEnergyValid

protected

Definition at line 146 of file PndEmcCluster.h.

Referenced by energy(), invalidateCache(), IsEnergyValid(), and SetEnergy().

Double_t PndEmcCluster::fLatMom

protected

Definition at line 153 of file PndEmcCluster.h.

Referenced by LatMom(), and SetLatMom().

std::map<Int_t,Int_t> PndEmcCluster::fLocalMaxMap

protected

Definition at line 144 of file PndEmcCluster.h.

Referenced by addLocalMax(), GetModule(), LocalMaxMap(), and PndEmcCluster().

std::vector<Int_t> PndEmcCluster::fMcList

mutableprotected

Definition at line 141 of file PndEmcCluster.h.

Referenced by GetMcIndex(), GetMcList(), GetMcSize(), and PndEmcCluster().

std::map<Int_t, Int_t> PndEmcCluster::fMcMap

mutableprotected

Definition at line 142 of file PndEmcCluster.h.

Referenced by addDigi(), and removeDigi().

std::map<Int_t,Int_t> PndEmcCluster::fMemberDigiMap

protected

Definition at line 143 of file PndEmcCluster.h.

Referenced by addDigi(), addLocalMax(), GetModule(), MemberDigiMap(), PndEmcCluster(), and removeDigi().

unsigned PndEmcCluster::fNbumps

protected

Definition at line 150 of file PndEmcCluster.h.

Referenced by NBumps(), and SetNBumps().

FairMultiLinkedData PndEmcCluster::fTrackEntering

protected

Definition at line 155 of file PndEmcCluster.h.

Referenced by AddTracksEnteringExiting(), GetTrackEntering(), Print(), and SetTrackEntering().

FairMultiLinkedData PndEmcCluster::fTrackExiting

protected

Definition at line 156 of file PndEmcCluster.h.

Referenced by AddTracksEnteringExiting(), GetTrackExiting(), Print(), and SetTrackExiting().

TVector3 PndEmcCluster::fWhere

protected

Definition at line 149 of file PndEmcCluster.h.

Referenced by SetPosition(), and where().

bool PndEmcCluster::fWhereValid

protected

Definition at line 148 of file PndEmcCluster.h.

Referenced by invalidateCache(), IsPositionValid(), SetPosition(), and where().

Double_t PndEmcCluster::fZ20

protected

Definition at line 151 of file PndEmcCluster.h.

Referenced by SetZ20(), and Z20().

Double_t PndEmcCluster::fZ53

protected

Definition at line 152 of file PndEmcCluster.h.

Referenced by SetZ53(), and Z53().

---

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

• PndEmcCluster.h
• PndEmcCluster.cxx