PndEmcClusterMoments Class Reference

#include <PndEmcClusterMoments.h>

Inheritance diagram for PndEmcClusterMoments:

Public Member Functions
	PndEmcClusterMoments (const PndEmcCluster &cluster, const TClonesArray *digiArray)
	PndEmcClusterMoments (const PndEmcClusterMoments &)
virtual	~PndEmcClusterMoments ()
virtual double	SecondMoment () const
virtual double	SecondMomentTP () const
virtual double	Theta1 () const
virtual double	Phi1 () const
virtual double	Theta2 () const
virtual double	Phi2 () const
virtual double	Major1 () const
virtual double	Major2 () const
virtual double	Minor1 () const
virtual double	Minor2 () const
virtual TVector3	Centre1 () const
virtual TVector3	Centre2 () const
virtual void	Print (const Option_t *opt="") const

Protected Member Functions
const PndEmcCluster &	MyCluster () const
const std::map< Int_t, Int_t > &	Members () const
const TClonesArray *	DigiArray () const

Private Member Functions
	PndEmcClusterMoments ()
PndEmcClusterMoments &	operator= (const PndEmcClusterMoments &other)

Detailed Description

Definition at line 32 of file PndEmcClusterMoments.h.

Constructor & Destructor Documentation

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

Definition at line 38 of file PndEmcClusterMoments.cxx.

                                                                                                     :
        PndEmcAbsClusterProperty(cluster, digiArray)            
{
}

PndEmcClusterMoments::PndEmcClusterMoments

(

const PndEmcClusterMoments &

other

)

Definition at line 43 of file PndEmcClusterMoments.cxx.

                                                                           :
        PndEmcAbsClusterProperty( other )
{}

PndEmcClusterMoments::~PndEmcClusterMoments ( )

virtual

Definition at line 51 of file PndEmcClusterMoments.cxx.

{}

PndEmcClusterMoments::PndEmcClusterMoments ( )

private

Member Function Documentation

TVector3 PndEmcClusterMoments::Centre1 ( ) const

virtual

Definition at line 295 of file PndEmcClusterMoments.cxx.

References digi, PndEmcAbsClusterProperty::DigiArray(), PndEmcCluster::energy(), PndEmcDigi::GetEnergy(), PndEmcAbsClusterProperty::Members(), PndEmcAbsClusterProperty::MyCluster(), and PndEmcDigi::where().

{
        TVector3 aVector(0,0,0);
        TVector3 clusterCentre( MyCluster().where() );
        
        std::map<Int_t, Int_t>::const_iterator iter = Members().begin();
        while ( iter != Members().end() ) {
                PndEmcDigi *digi = (PndEmcDigi *) DigiArray()->At(iter->second);
                aVector += ( digi->where() - clusterCentre ) * digi->GetEnergy();
                ++iter;
        }
        aVector *= 1./MyCluster().energy();
        
        return aVector;
}

TVector3 PndEmcClusterMoments::Centre2 ( ) const

virtual

Definition at line 312 of file PndEmcClusterMoments.cxx.

References digi, PndEmcAbsClusterProperty::DigiArray(), PndEmcCluster::energy(), PndEmcDigi::GetEnergy(), PndEmcAbsClusterProperty::Members(), PndEmcAbsClusterProperty::MyCluster(), and PndEmcDigi::where().

{
        TVector3 aVector(0,0,0);
        TVector3 clusterCentre( MyCluster().where() );
        
        std::map<Int_t, Int_t>::const_iterator iter = Members().begin();
        while ( iter != Members().end() ) {
                PndEmcDigi *digi = (PndEmcDigi *) DigiArray()->At(iter->second);
                TVector3 displacement( digi->where() - clusterCentre );
                aVector += TVector3(displacement.X()*displacement.X(), displacement.Y()*displacement.Y(),displacement.Z()*displacement.Z()) * digi->GetEnergy();
                ++iter;
        }
        aVector *= 1./MyCluster().energy();
        
        return aVector;
}

const TClonesArray * PndEmcAbsClusterProperty::DigiArray ( ) const

protectedinherited

Definition at line 49 of file PndEmcAbsClusterProperty.cxx.

References PndEmcAbsClusterProperty::fDigiArray.

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

{
        return fDigiArray;
}

Double_t PndEmcClusterMoments::Major1 ( ) const

virtual

Definition at line 182 of file PndEmcClusterMoments.cxx.

References cos(), digi, PndEmcAbsClusterProperty::DigiArray(), Double_t, PndEmcCluster::energy(), fabs(), PndEmcCluster::FindPhiDiff(), PndEmcDigi::GetEnergy(), PndEmcDigi::GetPhi(), PndEmcDigi::GetTheta(), PndEmcClusterProperties::Major_axis(), PndEmcAbsClusterProperty::Members(), PndEmcAbsClusterProperty::MyCluster(), PndEmcCluster::phi(), sin(), and PndEmcCluster::theta().

{
        Double_t moment=0;
        PndEmcClusterProperties *properties = new PndEmcClusterProperties(MyCluster(), DigiArray());
        Double_t axis=properties->Major_axis();
        Double_t clusT=MyCluster().theta(),clusP=MyCluster().phi();
        
        if ( axis==-999. ) return ( 0. );
        
        TVector2 maj(cos( axis ), sin( axis ));
        
        std::map<Int_t, Int_t>::const_iterator iter = Members().begin();
        
        while ( iter != Members().end() ) {
                PndEmcDigi *digi = (PndEmcDigi *) DigiArray()->At(iter->second);
                TVector2 coord(PndEmcCluster::FindPhiDiff(digi->GetPhi(),clusP), digi->GetTheta()-clusT);

                moment += fabs( coord*maj ) * digi->GetEnergy();
                ++iter;
        }

        moment/=MyCluster().energy();
        
        return( moment );
}

Double_t PndEmcClusterMoments::Major2 ( ) const

virtual

Definition at line 210 of file PndEmcClusterMoments.cxx.

References cos(), digi, PndEmcAbsClusterProperty::DigiArray(), Double_t, PndEmcCluster::energy(), PndEmcCluster::FindPhiDiff(), PndEmcDigi::GetEnergy(), PndEmcDigi::GetPhi(), PndEmcDigi::GetTheta(), PndEmcClusterProperties::Major_axis(), PndEmcAbsClusterProperty::Members(), PndEmcAbsClusterProperty::MyCluster(), PndEmcCluster::phi(), sin(), and PndEmcCluster::theta().

{
        Double_t moment=0;
        PndEmcClusterProperties *properties = new PndEmcClusterProperties(MyCluster(), DigiArray());
        Double_t axis=properties->Major_axis();
        Double_t clusT=MyCluster().theta(),clusP=MyCluster().phi();
        
        if ( axis==-999. ) return ( 0. );
        
        TVector2 maj(cos( axis ), sin( axis ));
        
        std::map<Int_t, Int_t>::const_iterator iter = Members().begin();
        
        while ( iter != Members().end() ) {
                PndEmcDigi *digi = (PndEmcDigi *) DigiArray()->At(iter->second);
                TVector2 coord(PndEmcCluster::FindPhiDiff(digi->GetPhi(),clusP), digi->GetTheta()-clusT);

                moment+=(coord*maj)*(coord*maj)*digi->GetEnergy();
                ++iter;
        }

        moment/=MyCluster().energy();
        
        return( moment );
}

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

protectedinherited

Definition at line 44 of file PndEmcAbsClusterProperty.cxx.

References PndEmcAbsClusterProperty::fMembers.

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

{
  return fMembers;
}

Double_t PndEmcClusterMoments::Minor1 ( ) const

virtual

Definition at line 237 of file PndEmcClusterMoments.cxx.

References cos(), digi, PndEmcAbsClusterProperty::DigiArray(), Double_t, PndEmcCluster::energy(), fabs(), PndEmcCluster::FindPhiDiff(), PndEmcDigi::GetEnergy(), PndEmcDigi::GetPhi(), PndEmcDigi::GetTheta(), PndEmcClusterProperties::Major_axis(), PndEmcAbsClusterProperty::Members(), min(), PndEmcAbsClusterProperty::MyCluster(), PndEmcCluster::phi(), Pi, sin(), and PndEmcCluster::theta().

{
        Double_t moment=0;
        PndEmcClusterProperties *properties = new PndEmcClusterProperties(MyCluster(), DigiArray());
        Double_t axis=properties->Major_axis();
        Double_t clusT=MyCluster().theta(),clusP=MyCluster().phi();
        
        
        if ( axis==-999. ) return ( 0. );
        
        axis+=TMath::Pi()/2.0;
        TVector2 min(cos( axis ), sin( axis ));
        
        std::map<Int_t, Int_t>::const_iterator iter = Members().begin();
        
        while ( iter != Members().end() ) {
                PndEmcDigi *digi = (PndEmcDigi *) DigiArray()->At(iter->second);
                TVector2 coOrd(PndEmcCluster::FindPhiDiff(digi->GetPhi(),clusP), digi->GetTheta()-clusT);

                moment+=fabs( coOrd*min )*digi->GetEnergy();
                ++iter;
        }
        
        moment/=MyCluster().energy();
        
        return( moment );
}

Double_t PndEmcClusterMoments::Minor2 ( ) const

virtual

Definition at line 267 of file PndEmcClusterMoments.cxx.

References cos(), digi, PndEmcAbsClusterProperty::DigiArray(), Double_t, PndEmcCluster::energy(), PndEmcCluster::FindPhiDiff(), PndEmcDigi::GetEnergy(), PndEmcDigi::GetPhi(), PndEmcDigi::GetTheta(), PndEmcClusterProperties::Major_axis(), PndEmcAbsClusterProperty::Members(), min(), PndEmcAbsClusterProperty::MyCluster(), PndEmcCluster::phi(), Pi, sin(), and PndEmcCluster::theta().

{
        Double_t moment=0;
        PndEmcClusterProperties *properties = new PndEmcClusterProperties(MyCluster(), DigiArray());
        Double_t axis=properties->Major_axis();
        Double_t clusT=MyCluster().theta(),clusP=MyCluster().phi();
                
        if ( axis==-999. ) return ( 0. );
        
        axis+=TMath::Pi()/2.0;
        TVector2 min(cos( axis ), sin( axis));
        
        std::map<Int_t, Int_t>::const_iterator iter = Members().begin();
        
        while ( iter != Members().end() ) {
                PndEmcDigi *digi = (PndEmcDigi *) DigiArray()->At(iter->second);
                TVector2 coOrd(PndEmcCluster::FindPhiDiff(digi->GetPhi(),clusP), digi->GetTheta()-clusT);
        
                moment+=(coOrd*min)*(coOrd*min)*digi->GetEnergy();
                iter++;
        }
        
        moment/=MyCluster().energy();
        
        return( moment );
}

const PndEmcCluster & PndEmcAbsClusterProperty::MyCluster ( ) const

protectedinherited

Definition at line 38 of file PndEmcAbsClusterProperty.cxx.

References PndEmcAbsClusterProperty::fMyCluster.

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

{
  return fMyCluster;
}

PndEmcClusterMoments& PndEmcClusterMoments::operator= ( const PndEmcClusterMoments &  other )

private

Double_t PndEmcClusterMoments::Phi1 ( ) const

virtual

Definition at line 121 of file PndEmcClusterMoments.cxx.

References digi, PndEmcAbsClusterProperty::DigiArray(), Double_t, PndEmcCluster::energy(), PndEmcCluster::FindPhiDiff(), PndEmcDigi::GetEnergy(), PndEmcDigi::GetPhi(), PndEmcAbsClusterProperty::Members(), PndEmcAbsClusterProperty::MyCluster(), p, and PndEmcCluster::phi().

Referenced by PndEmcClusterProperties::Major_axis().

{
        Double_t p=0;
        Double_t clus_phi=MyCluster().phi();
        
        std::map<Int_t, Int_t>::const_iterator iter = Members().begin();
        
        while ( iter != Members().end() ) {
                PndEmcDigi *digi = (PndEmcDigi *) DigiArray()->At(iter->second);
                p+=digi->GetEnergy()*PndEmcCluster::FindPhiDiff(digi->GetPhi(),clus_phi);
                ++iter;
        }
        
        p/=MyCluster().energy();
        
        return(p);
}

Double_t PndEmcClusterMoments::Phi2 ( ) const

virtual

Definition at line 161 of file PndEmcClusterMoments.cxx.

References digi, PndEmcAbsClusterProperty::DigiArray(), Double_t, PndEmcCluster::energy(), PndEmcCluster::FindPhiDiff(), PndEmcDigi::GetEnergy(), PndEmcDigi::GetPhi(), PndEmcAbsClusterProperty::Members(), PndEmcAbsClusterProperty::MyCluster(), p, PndEmcCluster::phi(), and phi_diff.

{
        Double_t p=0;
        Double_t clus_phi=MyCluster().phi();
        
        std::map<Int_t, Int_t>::const_iterator iter = Members().begin();
        
        while ( iter != Members().end() ) {
                PndEmcDigi *digi = (PndEmcDigi *) DigiArray()->At(iter->second);
                Double_t phi_diff=PndEmcCluster::FindPhiDiff(digi->GetPhi(),clus_phi);
                p+=digi->GetEnergy()*phi_diff*phi_diff;
                ++iter;
        }
        
        p/=MyCluster().energy();
        
        return(p);
 
}

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

virtual

Definition at line 330 of file PndEmcClusterMoments.cxx.

{
//  std::cout << " centre1=" << Centre1() << ", centre2=" << Centre2();
//  std::cout << endl;
//  std::cout << " theta1=" << Theta1() << ", phi1=" << Phi1();
//  std::cout << ", theta2=" << Theta2() << ", phi2=" << Phi2();
//  std::cout << endl;
//  std::cout << " major1=" << Major1() << ", minor1=" << Minor1();
//  std::cout << ", major2=" << Major2() << ", minor2=" << Minor2();
//  std::cout << endl;
//  std::cout << " secondMoment=" << SecondMoment() << ", secondMomentTP=" << SecondMomentTP();
//  std::cout << endl;
}

Double_t PndEmcClusterMoments::SecondMoment ( ) const

virtual

Definition at line 59 of file PndEmcClusterMoments.cxx.

References digi, PndEmcAbsClusterProperty::DigiArray(), Double_t, dx, dy, dz, PndEmcCluster::energy(), PndEmcDigi::GetEnergy(), PndEmcAbsClusterProperty::Members(), PndEmcAbsClusterProperty::MyCluster(), and PndEmcDigi::where().

{
        Double_t sum=0;
        
        std::map<Int_t, Int_t>::const_iterator iter = Members().begin();
        
        TVector3 clusterLocation( MyCluster().where() );
        while ( iter != Members().end() ) {
                PndEmcDigi *digi = (PndEmcDigi *) DigiArray()->At(iter->second);
                TVector3 digiLocation( digi->where() );
                Double_t dx = digiLocation.x() - clusterLocation.x();
                Double_t dy = digiLocation.y() - clusterLocation.y();
                Double_t dz = digiLocation.z() - clusterLocation.z();    
                sum+=digi->GetEnergy()*(dx*dx + dy*dy + dz*dz);
                ++iter;
        }

        sum/=MyCluster().energy();
        return sum;
}

Double_t PndEmcClusterMoments::SecondMomentTP ( ) const

virtual

Definition at line 81 of file PndEmcClusterMoments.cxx.

References digi, PndEmcAbsClusterProperty::DigiArray(), Double_t, PndEmcCluster::energy(), PndEmcDigi::GetEnergy(), PndEmcAbsClusterProperty::Members(), PndEmcAbsClusterProperty::MyCluster(), PndEmcCluster::where(), and PndEmcDigi::where().

{
        TVector3 cluster = MyCluster().where();
        
        Double_t sum =0;
        std::map<Int_t, Int_t>::const_iterator iter = Members().begin();
        while ( iter != Members().end() ) {
                PndEmcDigi *digi = (PndEmcDigi *) DigiArray()->At(iter->second);
                Double_t da = cluster.Angle(digi->where());
                sum+=digi->GetEnergy()*da*da;
                ++iter;
        }

        sum/=MyCluster().energy();
        return sum;
}

Double_t PndEmcClusterMoments::Theta1 ( ) const

virtual

Definition at line 99 of file PndEmcClusterMoments.cxx.

References digi, PndEmcAbsClusterProperty::DigiArray(), Double_t, PndEmcCluster::energy(), PndEmcDigi::GetEnergy(), PndEmcDigi::GetTheta(), PndEmcAbsClusterProperty::Members(), PndEmcAbsClusterProperty::MyCluster(), t, and theta.

Referenced by PndEmcClusterProperties::Major_axis().

{
        Double_t t=0;
        
        std::map<Int_t, Int_t>::const_iterator iter = Members().begin();
        
        Double_t clusterTheta(MyCluster().theta());
        
        while ( iter != Members().end() ) {
                PndEmcDigi *digi = (PndEmcDigi *) DigiArray()->At(iter->second);
                t+=digi->GetEnergy()*(digi->GetTheta()- clusterTheta);
                ++iter;
        }
        
        t/=MyCluster().energy();
        
        return t;
 
}

Double_t PndEmcClusterMoments::Theta2 ( ) const

virtual

Definition at line 140 of file PndEmcClusterMoments.cxx.

References digi, PndEmcAbsClusterProperty::DigiArray(), Double_t, PndEmcCluster::energy(), PndEmcDigi::GetEnergy(), PndEmcDigi::GetTheta(), PndEmcAbsClusterProperty::Members(), PndEmcAbsClusterProperty::MyCluster(), t, theta, and theta_diff.

{
        Double_t t=0;
        
        std::map<Int_t, Int_t>::const_iterator iter = Members().begin();
        Double_t clusterTheta(MyCluster().theta());
        
        while ( iter != Members().end() ) {
                PndEmcDigi *digi = (PndEmcDigi *) DigiArray()->At(iter->second);
                Double_t theta_diff=digi->GetTheta()-clusterTheta;
                t+=digi->GetEnergy()*theta_diff*theta_diff;
                ++iter;
        }

        t/=MyCluster().energy();
        
        return(t);
 
}

---

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

• PndEmcClusterMoments.h
• PndEmcClusterMoments.cxx