CCaloOrganization Class Reference

#include <CCaloOrganization.hh>

Public Member Functions
	CCaloOrganization ()
virtual	~CCaloOrganization ()
unsigned int	packindex (int det, int z, int eta, int phi) const
unsigned int	packindex (int det, int depth, int z, int eta, int phi) const
void	unpackindex (const unsigned int &idx, int &det, int &z, int &eta, int &phi) const
void	unpackindex (const unsigned int &idx, int &det, int &depth, int &z, int &eta, int &phi) const
int	getUnitWithMaxEnergy (std::map< int, float, std::less< int > > &themap)
float	energyInMatrix (int nCellInEta, int nCellInPhi, int crystalWithMaxEnergy, std::map< int, float, std::less< int > > &themap)

Detailed Description

Definition at line 38 of file CCaloOrganization.hh.

Constructor & Destructor Documentation

CCaloOrganization::CCaloOrganization ( )

inline

Definition at line 41 of file CCaloOrganization.hh.

{};

virtual CCaloOrganization::~CCaloOrganization ( )

inlinevirtual

Definition at line 42 of file CCaloOrganization.hh.

{};

Member Function Documentation

float CCaloOrganization::energyInMatrix	(	int	nCellInEta,
		int	nCellInPhi,
		int	crystalWithMaxEnergy,
		std::map< int, float, std::less< int > > &	themap
	)

Definition at line 121 of file CCaloOrganization.cc.

References G4cout, G4endl, packindex(), unpackindex(), and z.

                                                            {

  int det,z,eta,phi;
  this->unpackindex(crystalWithMaxEnergy, det, z, eta, phi);
  int ncristals=0;
    
  int goBackInEta = nCellInEta/2;
  int goBackInPhi = nCellInPhi/2;
  int startEta = eta-goBackInEta;
  int startPhi = phi-goBackInPhi;

  float totalEnergy = 0.;
  
  for(int ieta=startEta; ieta<startEta+nCellInEta; ieta++){
    for(int iphi=startPhi; iphi<startPhi+nCellInPhi; iphi++){
      
      int index = this->packindex(det,z,ieta,iphi);
      totalEnergy += themap[index];
      ncristals+=1;
#ifdef debug
      G4cout << "ieta - iphi - E = " << ieta << "  " << iphi << " " 
         << themap[index] << G4endl;
#endif
    }
  }
    
#ifdef debug
  G4cout << "Energy in " << nCellInEta << " cells in eta times "
     << nCellInPhi << " cells in phi matrix = " << totalEnergy
     << " for " << ncristals << " crystals" << G4endl;
#endif
  return totalEnergy;

}   

int CCaloOrganization::getUnitWithMaxEnergy

(

std::map< int, float, std::less< int > > &

themap

)

Definition at line 94 of file CCaloOrganization.cc.

References G4cout, G4endl, unpackindex(), and z.

                                                                                  {

  //look for max
  int UnitWithMaxEnergy = 0;
  float maxEnergy = 0.;
    
  for(std::map<int,float,std::less<int> >::iterator iter = themap.begin();
      iter != themap.end(); iter++){
        
    if( maxEnergy < (*iter).second) {
      maxEnergy = (*iter).second;   
      UnitWithMaxEnergy = (*iter).first;
    }               
  } 
#ifdef debug
  G4cout << " *** max energy of " << maxEnergy << " MeV was found in Unit id "
     << UnitWithMaxEnergy;
  int det,z,eta,phi;
  unpackindex(UnitWithMaxEnergy, det, z, eta, phi);
  G4cout << " corresponding to z= " << z << " eta= " << eta << " phi = " << phi
     << G4endl;
#endif
  return UnitWithMaxEnergy;

}

unsigned int CCaloOrganization::packindex	(	int	det,
		int	z,
		int	eta,
		int	phi
	)		const

Definition at line 35 of file CCaloOrganization.cc.

References G4cout, and G4endl.

Referenced by energyInMatrix().

                                     {
  //So this is the actual encoding of the index:
  //top 4 bits encode Detector type
  //
  // Should work for all calorimeter with no depth information

  unsigned int idx=(det&15)<<28;   //bits 28-31   (21-23 are free for now)
  idx+=(((z+1)/2)&1)<<20;          //bits 20
  idx+=(eta&1023)<<10;             //bits 10-19
  idx+=(phi&1023);                 //bits  0-9
#ifdef debug
  G4cout << " ECAL packing " << det << " " << z << " " << eta << " " << phi 
     << "  into " << idx << G4endl;
#endif
  return idx;
}

unsigned int CCaloOrganization::packindex	(	int	det,
		int	depth,
		int	z,
		int	eta,
		int	phi
	)		const

Definition at line 53 of file CCaloOrganization.cc.

References G4cout, and G4endl.

                                     {
  //So this is the actual encoding of the index:
  //top 4 bits encode Detector type
  //next 4 bits encode depth information
  // Should work for all calorimeter with no depth information

  unsigned int idx=(det&15)<<28;  //bits 28-31   (21-23 are free for now)
  idx+=(depth&15)<<24;            //bits 24-27  
  idx+=(z&1)<<20;                 //bits 20
  idx+=(eta&1023)<<10;            //bits 10-19
  idx+=(phi&1023);                //bits  0-9
#ifdef debug
  G4cout << " HCAL packing " << det << " " << depth << " " << z << " " << eta 
     << " " << phi  << "  into " << idx << G4endl;
#endif
  return idx;
}

void CCaloOrganization::unpackindex	(	const unsigned int &	idx,
		int &	det,
		int &	z,
		int &	eta,
		int &	phi
	)		const

Definition at line 73 of file CCaloOrganization.cc.

Referenced by energyInMatrix(), and getUnitWithMaxEnergy().

                                              {
  det = (idx>>28)&15;
  z   = (idx>>20)&1;
  z   = 2*z-1;
  eta = (idx>>10)&1023;
  phi = (idx&1023);
}

void CCaloOrganization::unpackindex	(	const unsigned int &	idx,
		int &	det,
		int &	depth,
		int &	z,
		int &	eta,
		int &	phi
	)		const

Definition at line 83 of file CCaloOrganization.cc.

                                    {
  det = (idx>>28)&15;
  depth=(idx>>24)&15;
  z   = (idx>>20)&1;
  eta = (idx>>10)&1023;
  phi = (idx&1023);
}

---

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

• CCaloOrganization.hh
• CCaloOrganization.cc