PndEmcXClMoments Class Reference

#include <PndEmcXClMoments.h>

Inheritance diagram for PndEmcXClMoments:

Public Member Functions
	PndEmcXClMoments (const PndEmcCluster &cluster, const TClonesArray *digiArray)
	PndEmcXClMoments (const PndEmcXClMoments &)
virtual	~PndEmcXClMoments ()
virtual Double_t	SecondMoment () const
virtual Double_t	SecondMomentPhi () const
virtual Double_t	SecondMomentTheta () const
virtual Double_t	AbsZernikeMoment (int n, int m, Double_t R0=15) const
virtual Double_t	Lat () const
virtual void	Print (const Option_t *opt="") const

Protected Member Functions
void	Init ()
Double_t	f00 (Double_t r) const
Double_t	f11 (Double_t r) const
Double_t	f20 (Double_t r) const
Double_t	f22 (Double_t r) const
Double_t	f31 (Double_t r) const
Double_t	f33 (Double_t r) const
Double_t	f40 (Double_t r) const
Double_t	f42 (Double_t r) const
Double_t	f44 (Double_t r) const
Double_t	f51 (Double_t r) const
Double_t	f53 (Double_t r) const
Double_t	f55 (Double_t r) const
Double_t	Fast_AbsZernikeMoment (int n, int m, Double_t R0) const
Double_t	Calc_AbsZernikeMoment (int n, int m, Double_t R0) const
Double_t	Fak (int n) const
const PndEmcCluster &	MyCluster () const
const std::map< Int_t, Int_t > &	Members () const
const TClonesArray *	DigiArray () const

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

Private Attributes
std::vector < PndEmcClEnergyDeposition > *	fEnergyDistribution
Int_t	fClusterSize
Double_t(PndEmcXClMoments::*	fFcn [12])(Double_t) const

Detailed Description

Definition at line 39 of file PndEmcXClMoments.h.

Constructor & Destructor Documentation

PndEmcXClMoments::PndEmcXClMoments	(	const PndEmcCluster &	cluster,
		const TClonesArray *	digiArray
	)

Definition at line 38 of file PndEmcXClMoments.cxx.

References Init().

                                                                                             :
        PndEmcAbsClusterProperty( cluster, digiArray ), fEnergyDistribution(), fClusterSize(0)  
{
        Init();
}

PndEmcXClMoments::PndEmcXClMoments

(

const PndEmcXClMoments &

other

)

Definition at line 44 of file PndEmcXClMoments.cxx.

References Init().

                                                                :
        PndEmcAbsClusterProperty( other ), fEnergyDistribution(other.fEnergyDistribution), fClusterSize(other.fClusterSize)
{
        Init();
}

PndEmcXClMoments::~PndEmcXClMoments ( )

virtual

Definition at line 54 of file PndEmcXClMoments.cxx.

References fEnergyDistribution.

{
        fEnergyDistribution->clear();
        delete fEnergyDistribution;
        fEnergyDistribution = 0;
}

PndEmcXClMoments::PndEmcXClMoments ( )

private

Member Function Documentation

Double_t PndEmcXClMoments::AbsZernikeMoment ( int  n, int  m, Double_t  R0 = 15  ) const

virtual

Definition at line 104 of file PndEmcXClMoments.cxx.

References Calc_AbsZernikeMoment(), and Fast_AbsZernikeMoment().

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

{
        // 1. Check if n,m are correctly
        if ((m>n) || ((n-m)%2 != 0) || (n<0) || (m<0)) 
        { 
                return -1;
        }
        
        // 2. Check if n,R0 are within validity Range :
        // n>20 or R0<5cm  just makes so sense !
        if ((n>20) || (R0<=5)) 
        {
                return -1;
        }
        if (n<=5) 
        {
                return Fast_AbsZernikeMoment(n,m,R0);
        } else {
                return Calc_AbsZernikeMoment(n,m,R0);
        }
}

Double_t PndEmcXClMoments::Calc_AbsZernikeMoment ( int  n, int  m, Double_t  R0  ) const

protected

Definition at line 264 of file PndEmcXClMoments.cxx.

References cos(), Double_t, PndEmcCluster::energy(), Fak(), fClusterSize, i, m, PndEmcAbsClusterProperty::MyCluster(), phi, r, s, sin(), and sqrt().

Referenced by AbsZernikeMoment().

{
        Double_t r,ph,e,Re=0,Im=0,f_nm,result;
        Double_t TotalEnergy = MyCluster().energy();
        for (int i=0; i<fClusterSize; i++) 
        {
                r = (*fEnergyDistribution)[i].r / R0;
                if (r<1) 
                { 
                        ph = ((*fEnergyDistribution)[i]).phi;
                        e = (*fEnergyDistribution)[i].deposited_energy;  
                        f_nm=0;
                        for (int s=0; s<=(n-m)/2; s++) {
                                if (s%2==0) {
                                        f_nm = f_nm + Fak(n-s)*pow(r,(Double_t) (n-2*s))/(Fak(s)*Fak((n+m)/2-s)*Fak((n-m)/2-s));
                                }else {
                                        f_nm = f_nm - Fak(n-s)*pow(r,(Double_t) (n-2*s))/(Fak(s)*Fak((n+m)/2-s)*Fak((n-m)/2-s));
                                }
                        }
                        Re = Re + e/TotalEnergy * f_nm * cos( (Double_t) m*ph);
                        Im = Im - e/TotalEnergy * f_nm * sin( (Double_t) m*ph);
                }
        }
        result = sqrt(Re*Re+Im*Im);
        
        return result;
}

const TClonesArray * PndEmcAbsClusterProperty::DigiArray ( ) const

protectedinherited

Definition at line 49 of file PndEmcAbsClusterProperty.cxx.

References PndEmcAbsClusterProperty::fDigiArray.

Referenced by PndEmcClusterMoments::Centre1(), PndEmcClusterMoments::Centre2(), PndEmcClusterEnergySums::E1(), PndEmcClusterEnergySums::E25(), PndEmcClusterEnergySums::E9(), PndEmcClusterProperties::Energy(), PndEmcClusterEnergySums::energy(), PndEmcClusterProperties::GravWhere(), Init(), PndEmcClusterProperties::LiloWhere(), PndEmcClusterProperties::LinearWhere(), PndEmcClusterMoments::Major1(), PndEmcClusterMoments::Major2(), PndEmcClusterProperties::Major_axis(), PndEmcClusterProperties::Mass(), PndEmcClusterMoments::Minor1(), PndEmcClusterMoments::Minor2(), PndEmcClusterMoments::Phi1(), PndEmcClusterMoments::Phi2(), Print(), PndEmcClusterMoments::SecondMoment(), PndEmcClusterMoments::SecondMomentTP(), PndEmcClusterMoments::Theta1(), and PndEmcClusterMoments::Theta2().

{
        return fDigiArray;
}

Double_t PndEmcXClMoments::f00 ( Double_t  r ) const

protected

Definition at line 206 of file PndEmcXClMoments.cxx.

Referenced by Init().

{ return 1; } // r   //[R.K.03/2017] unused variable(s)

Double_t PndEmcXClMoments::f11 ( Double_t  r ) const

protected

Definition at line 209 of file PndEmcXClMoments.cxx.

References r.

Referenced by Init().

{ return r; }

Double_t PndEmcXClMoments::f20 ( Double_t  r ) const

protected

Definition at line 212 of file PndEmcXClMoments.cxx.

Referenced by Init().

{ return 2.0*r*r-1.0; }

Double_t PndEmcXClMoments::f22 ( Double_t  r ) const

protected

Definition at line 215 of file PndEmcXClMoments.cxx.

References r.

Referenced by Init().

{ return r*r; }

Double_t PndEmcXClMoments::f31 ( Double_t  r ) const

protected

Definition at line 218 of file PndEmcXClMoments.cxx.

References r.

Referenced by Init().

{ return 3.0*r*r*r - 2.0*r; }

Double_t PndEmcXClMoments::f33 ( Double_t  r ) const

protected

Definition at line 221 of file PndEmcXClMoments.cxx.

References r.

Referenced by Init().

{ return r*r*r; }

Double_t PndEmcXClMoments::f40 ( Double_t  r ) const

protected

Definition at line 224 of file PndEmcXClMoments.cxx.

Referenced by Init().

{ return 6.0*r*r*r*r-6.0*r*r+1.0; }

Double_t PndEmcXClMoments::f42 ( Double_t  r ) const

protected

Definition at line 227 of file PndEmcXClMoments.cxx.

References r.

Referenced by Init().

{ return 4.0*r*r*r*r-3.0*r*r; }

Double_t PndEmcXClMoments::f44 ( Double_t  r ) const

protected

Definition at line 230 of file PndEmcXClMoments.cxx.

References r.

Referenced by Init().

{ return r*r*r*r; }

Double_t PndEmcXClMoments::f51 ( Double_t  r ) const

protected

Definition at line 233 of file PndEmcXClMoments.cxx.

References r.

Referenced by Init().

{ return 10.0*pow(r,5)-12.0*pow(r,3)+3.0*r; }

Double_t PndEmcXClMoments::f53 ( Double_t  r ) const

protected

Definition at line 236 of file PndEmcXClMoments.cxx.

Referenced by Init().

{ return 5.0*pow(r,5) - 4.0*pow(r,3); }

Double_t PndEmcXClMoments::f55 ( Double_t  r ) const

protected

Definition at line 239 of file PndEmcXClMoments.cxx.

Referenced by Init().

{ return pow(r,5); }

Double_t PndEmcXClMoments::Fak ( int  n ) const

protected

Definition at line 293 of file PndEmcXClMoments.cxx.

References Double_t, i, n, and res.

Referenced by Calc_AbsZernikeMoment().

{
        Double_t res=1.0;
        for (int i=2; i<=n; i++) 
        {
                res=res*(Double_t) i;
        }
        
        return res;
}

Double_t PndEmcXClMoments::Fast_AbsZernikeMoment ( int  n, int  m, Double_t  R0  ) const

protected

Definition at line 242 of file PndEmcXClMoments.cxx.

References cos(), Double_t, PndEmcCluster::energy(), fClusterSize, fFcn, i, PndEmcAbsClusterProperty::MyCluster(), phi, r, sin(), and sqrt().

Referenced by AbsZernikeMoment().

{
        Double_t r,ph,e,Re=0,Im=0,result;
        Double_t TotalEnergy = MyCluster().energy();
        int index = (n/2)*(n/2)+(n/2)+m; 
        for (int i=0; i<fClusterSize; i++)
        {
                r = (*fEnergyDistribution)[i].r / R0;
                if (r<1) 
                { 
                        ph = ((*fEnergyDistribution)[i]).phi;
                        e = (*fEnergyDistribution)[i].deposited_energy;
                        Re = Re + e/TotalEnergy * (this->*fFcn[index])(r) * cos( (Double_t) m * ph);
                        Im = Im - e/TotalEnergy * (this->*fFcn[index])(r) * sin( (Double_t) m * ph);
                }
        }
        result = sqrt(Re*Re+Im*Im);
        
        return result;
}

void PndEmcXClMoments::Init ( )

protected

Definition at line 305 of file PndEmcXClMoments.cxx.

References PndEmcClEnergyDeposition::deposited_energy, digi, PndEmcAbsClusterProperty::DigiArray(), PndEmcCluster::DigiList(), Double_t, f00(), f11(), f20(), f22(), f31(), f33(), f40(), f42(), f44(), f51(), f53(), f55(), fClusterSize, fEnergyDistribution, fFcn, PndEmcDigi::GetEnergy(), PndEmcDigi::getRescaleFactor(), PndEmcClusterProperties::GravWhere(), i, PndEmcCluster::Maxima(), PndEmcAbsClusterProperty::MyCluster(), PndEmcClEnergyDeposition::phi, PndEmcClEnergyDeposition::r, and PndEmcDigi::where().

Referenced by PndEmcXClMoments().

{
        // need to get one digi to ask what kind of digi position method is used
        // or to be more precise what scaling factor has do be used
        Double_t rescaleFactor = 1.;
        PndEmcClusterProperties *properties = new PndEmcClusterProperties(MyCluster(), DigiArray());

        const PndEmcDigi* pDigi = MyCluster().Maxima(DigiArray());
        if ( pDigi != 0 ) rescaleFactor = pDigi->getRescaleFactor();
        
        // make sure we always use gravWhere for the cluster moments
        //  TVector3 cl( MyCluster().where() );
        TVector3 cl( properties->GravWhere() );
        TVector3 ClusDirection(cl.x(),cl.y(),cl.z());
        ClusDirection *= 1.0/ClusDirection.Mag();
        TVector3 theta_axis(ClusDirection.y(),-ClusDirection.x(),0.0);
        theta_axis *= 1.0/theta_axis.Mag();
        TVector3 phi_axis = theta_axis.Cross(ClusDirection);
        
        std::vector<Int_t> digiList = MyCluster().DigiList();
        std::vector<Int_t>::iterator current;

        fClusterSize = digiList.size();
        fEnergyDistribution = new std::vector<PndEmcClEnergyDeposition> (fClusterSize);
        
        int i=0;
        PndEmcClEnergyDeposition clEdep;
 
        for (current=digiList.begin();current<digiList.end();++current)
        {
                PndEmcDigi *digi = (PndEmcDigi *) DigiArray()->At(*current);
                clEdep.deposited_energy = digi->GetEnergy();
                
                TVector3 diff = digi->where() - cl;
                TVector3 DigiVect = diff - diff.Dot(ClusDirection) * ClusDirection;
                clEdep.r = DigiVect.Mag() / rescaleFactor;
                clEdep.phi = DigiVect.Angle(theta_axis);
                if (DigiVect.Dot(phi_axis)<0) 
                { 
                        clEdep.phi = 2*M_PI - clEdep.phi; 
                }
        
                (*fEnergyDistribution)[i] = clEdep;
                i++;
        }
        
        fFcn[0]  = &PndEmcXClMoments::f00 ;
        fFcn[1]  = &PndEmcXClMoments::f11 ;
        fFcn[2]  = &PndEmcXClMoments::f20 ;
        fFcn[3]  = &PndEmcXClMoments::f31 ;
        fFcn[4]  = &PndEmcXClMoments::f22 ;
        fFcn[5]  = &PndEmcXClMoments::f33 ;
        fFcn[6]  = &PndEmcXClMoments::f40 ;
        fFcn[7]  = &PndEmcXClMoments::f51 ;
        fFcn[8]  = &PndEmcXClMoments::f42 ;
        fFcn[9]  = &PndEmcXClMoments::f53 ;
        fFcn[10] = &PndEmcXClMoments::f44 ;
        fFcn[11] = &PndEmcXClMoments::f55 ;
}

Double_t PndEmcXClMoments::Lat ( ) const

virtual

Definition at line 127 of file PndEmcXClMoments.cxx.

References Double_t, fClusterSize, fEnergyDistribution, i, n, and r.

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

{
        Double_t r,redmoment=0;
        int n,n1,n2,tmp;
        if (fClusterSize<3) 
        { 
                return 0; 
        }
        
        n1=0; n2=1;
        if ((*fEnergyDistribution)[1].deposited_energy > (*fEnergyDistribution)[0].deposited_energy) 
        {
                tmp=n2; n2=n1; n1=tmp;
        }
        for (int i=2; i<fClusterSize; i++) 
        {
                n=i;
                if ((*fEnergyDistribution)[i].deposited_energy > (*fEnergyDistribution)[n1].deposited_energy) 
                {
                        tmp = n2;
                        n2 = n1; n1 = i; n=tmp;
                } else {
                        if ((*fEnergyDistribution)[i].deposited_energy > (*fEnergyDistribution)[n2].deposited_energy) 
                        {
                                tmp=n2; n2=i; n=tmp;
                        }
                }
                r = (*fEnergyDistribution)[n].r;
                redmoment+= r*r* (*fEnergyDistribution)[n].deposited_energy;
        }
        Double_t e1 = (*fEnergyDistribution)[n1].deposited_energy;
        Double_t e2 = (*fEnergyDistribution)[n2].deposited_energy;
        Double_t r0 = 2.5;//cm, the average distance between crystals
        Double_t lat = redmoment/(redmoment+r0*r0*(e1+e2));
        
        return lat;
}

const std::map< Int_t, Int_t > & PndEmcAbsClusterProperty::Members ( ) const

protectedinherited

Definition at line 44 of file PndEmcAbsClusterProperty.cxx.

References PndEmcAbsClusterProperty::fMembers.

Referenced by PndEmcClusterMoments::Centre1(), PndEmcClusterMoments::Centre2(), PndEmcClusterMoments::Major1(), PndEmcClusterMoments::Major2(), PndEmcClusterMoments::Minor1(), PndEmcClusterMoments::Minor2(), PndEmcClusterMoments::Phi1(), PndEmcClusterMoments::Phi2(), PndEmcClusterMoments::SecondMoment(), PndEmcClusterMoments::SecondMomentTP(), PndEmcClusterMoments::Theta1(), and PndEmcClusterMoments::Theta2().

{
  return fMembers;
}

const PndEmcCluster & PndEmcAbsClusterProperty::MyCluster ( ) const

protectedinherited

Definition at line 38 of file PndEmcAbsClusterProperty.cxx.

References PndEmcAbsClusterProperty::fMyCluster.

Referenced by Calc_AbsZernikeMoment(), PndEmcClusterMoments::Centre1(), PndEmcClusterMoments::Centre2(), PndEmcClusterEnergySums::E1(), PndEmcClusterEnergySums::E25(), PndEmcClusterEnergySums::E9(), PndEmcClusterProperties::Energy(), PndEmcClusterEnergySums::energy(), Fast_AbsZernikeMoment(), PndEmcClusterProperties::GravWhere(), Init(), PndEmcClusterProperties::LiloWhere(), PndEmcClusterProperties::LinearWhere(), PndEmcClusterMoments::Major1(), PndEmcClusterMoments::Major2(), PndEmcClusterProperties::Major_axis(), PndEmcClusterProperties::Mass(), PndEmcClusterMoments::Minor1(), PndEmcClusterMoments::Minor2(), PndEmcClusterMoments::Phi1(), PndEmcClusterMoments::Phi2(), Print(), PndEmcClusterMoments::SecondMoment(), SecondMoment(), SecondMomentPhi(), SecondMomentTheta(), PndEmcClusterMoments::SecondMomentTP(), PndEmcClusterMoments::Theta1(), and PndEmcClusterMoments::Theta2().

{
  return fMyCluster;
}

PndEmcXClMoments& PndEmcXClMoments::operator= ( const PndEmcXClMoments &  )

private

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

virtual

Definition at line 166 of file PndEmcXClMoments.cxx.

References digi, PndEmcAbsClusterProperty::DigiArray(), PndEmcCluster::DigiList(), fFcn, i, PndEmcAbsClusterProperty::MyCluster(), PndEmcClEnergyDeposition::phi, PndEmcDigi::Print(), and PndEmcClEnergyDeposition::r.

{
        TVector3 cl( MyCluster().where() );
        TVector3 ClusDirection(cl.x(),cl.y(),cl.z());
        ClusDirection *= 1.0/ClusDirection.Mag();
        TVector3 theta_axis(ClusDirection.y(),-ClusDirection.x(),0.0);
        theta_axis *= 1.0/theta_axis.Mag();
        TVector3 phi_axis = theta_axis.Cross(ClusDirection);
        
        std::vector<Int_t> digiList = MyCluster().DigiList();
        std::vector<Int_t>::iterator current;
        
        PndEmcClEnergyDeposition clEdep;
        std::cout << "Dump of PndEmcCluster with CoG at (" << cl.x() << "," << cl.y()
                << "," << cl.z() << ")" << endl;
        std::cout << "-------------------------------------------------------------"<< endl;
        std::cout << " theta - axis : (" << theta_axis.x() << "," << theta_axis.y() << ","
                << theta_axis.z() << ")" << endl;
        std::cout << " phi - axis : (" << phi_axis.x() << "," << phi_axis.y() << ","
                << phi_axis.z() << ")" << endl;
        int i=0;
        
        for (current=digiList.begin();current<digiList.end();++current){
                clEdep = (*fEnergyDistribution)[i];
                PndEmcDigi *digi = (PndEmcDigi *) DigiArray()->At(*current);
                digi->Print();
                std::cout << "Polar coords : (" << clEdep.r << "," << clEdep.phi << ")" << endl;
                i++;
        }  
        
        for (int j=0; j<=8; j++)
        {
                std::cout << "f[" << j << "](5) = " << (this->*fFcn[j])(5.0) << endl;
        }
        
        std::cout << "--------------------------------------------------------------"<< endl;
}

Double_t PndEmcXClMoments::SecondMoment ( ) const

virtual

Definition at line 62 of file PndEmcXClMoments.cxx.

References Double_t, PndEmcCluster::energy(), fClusterSize, i, PndEmcAbsClusterProperty::MyCluster(), and r.

{
        Double_t r,sum=0;
        for (int i=0; i<fClusterSize; i++) 
        {
                r=(*fEnergyDistribution)[i].r;
                sum+=(*fEnergyDistribution)[i].deposited_energy*r*r;
        }
        sum/=MyCluster().energy();
        
        return sum;
}

Double_t PndEmcXClMoments::SecondMomentPhi ( ) const

virtual

Definition at line 76 of file PndEmcXClMoments.cxx.

References Double_t, PndEmcCluster::energy(), fClusterSize, fEnergyDistribution, i, PndEmcAbsClusterProperty::MyCluster(), phi, r, and sin().

{
        Double_t r,sum=0;
        for (int i=0; i<fClusterSize; i++) 
        {
                r=(*fEnergyDistribution)[i].r*sin((*fEnergyDistribution)[i].phi);
                sum+=(*fEnergyDistribution)[i].deposited_energy*r*r;
        }
        sum/=MyCluster().energy();
        
        return sum;
}

Double_t PndEmcXClMoments::SecondMomentTheta ( ) const

virtual

Definition at line 90 of file PndEmcXClMoments.cxx.

References cos(), Double_t, PndEmcCluster::energy(), fClusterSize, fEnergyDistribution, i, PndEmcAbsClusterProperty::MyCluster(), phi, and r.

{
        Double_t r,sum=0;
        for (int i=0; i<fClusterSize; i++) 
        {
                r=(*fEnergyDistribution)[i].r*cos((*fEnergyDistribution)[i].phi);
                sum+=(*fEnergyDistribution)[i].deposited_energy*r*r;
        }
        sum/=MyCluster().energy();
        
        return sum;
}  

Member Data Documentation

Int_t PndEmcXClMoments::fClusterSize

private

Definition at line 115 of file PndEmcXClMoments.h.

Referenced by Calc_AbsZernikeMoment(), Fast_AbsZernikeMoment(), Init(), Lat(), SecondMoment(), SecondMomentPhi(), and SecondMomentTheta().

std::vector<PndEmcClEnergyDeposition>* PndEmcXClMoments::fEnergyDistribution

private

Definition at line 114 of file PndEmcXClMoments.h.

Referenced by Init(), Lat(), SecondMomentPhi(), SecondMomentTheta(), and ~PndEmcXClMoments().

Double_t(PndEmcXClMoments::* PndEmcXClMoments::fFcn[12])(Double_t) const

private

Definition at line 116 of file PndEmcXClMoments.h.

Referenced by Fast_AbsZernikeMoment(), Init(), and Print().

---

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

• PndEmcXClMoments.h
• PndEmcXClMoments.cxx