Atlfast::CalSection Class Reference

Provides rudimentary geometry information for a calorimeter section. More...

#include <CalSection.h>

Collaboration diagram for Atlfast::CalSection:

[legend]List of all members.

Public Member Functions
	CalSection (MsgStream &log, const double minEta, const double maxEta, const double granEta, const double granPhi, const double minPhi=-M_PI, const double maxPhi=M_PI)
	CalSection (const CalSection &c)
	~CalSection ()
void	reset ()
void	deposit (ITransportedParticleCollectionIter &, ITransportedParticleCollectionIter &)
void	deposit (std::vector< Atlfast::EPileupDeposit * >::iterator &, std::vector< Atlfast::EPileupDeposit * >::iterator &)
void	depositEcal (std::vector< const GridletElement * >::iterator &, std::vector< const GridletElement * >::iterator &)
void	depositHcal (std::vector< const GridletElement * >::iterator &, std::vector< const GridletElement * >::iterator &)
void	depositEgen (std::vector< const Gridlet * >::iterator &, std::vector< const Gridlet * >::iterator &)
void	giveHits (const ICellSelector *, ITwoCptCellCollection *) const
void	smearCells (ISmearer *)
Private Member Functions
void	newHit (const ITransportedParticle *)
void	newHit (const EPileupDeposit *)
void	newEHit (const GridletElement *)
void	newHHit (const GridletElement *)
void	setEgen (const Gridlet *)
int	index (const Gridlet *) const
int	index (const GridletElement *) const
int	index (const ITransportedParticle *) const
int	index (const Atlfast::EPileupDeposit *) const
int	iindex (double phi, double eta) const
Private Attributes
MsgStream &	m_log
const double	m_minEta
const double	m_maxEta
const double	m_minPhi
const double	m_maxPhi
double	m_granEta
double	m_granPhi
int	m_nPhi
CalSectionReject	m_calSectionReject
std::vector< ITwoCptCell * >	m_cells
std::map< int, ITwoCptCell *, std::less< int > >	m_hitCells

---

Detailed Description

Provides rudimentary geometry information for a calorimeter section.

It is used by the Calorimeter object to deposit HepMc particles

Definition at line 43 of file CalSection.h.

---

Constructor & Destructor Documentation

Atlfast::CalSection::CalSection	(	MsgStream &	log,
		const double	minEta,
		const double	maxEta,
		const double	granEta,
		const double	granPhi,
		const double	minPhi = -M_PI,
		const double	maxPhi = M_PI
	)

Constructor takes sensitivity parameters

Definition at line 34 of file CalSection.cxx.

00042                           :
00043     m_log(log),
00044     m_minEta(minEta), m_maxEta(maxEta),
00045     m_minPhi(minPhi), m_maxPhi(maxPhi),
00046     m_calSectionReject(minEta, maxEta, minPhi, maxPhi){
00047     
00048     log <<MSG::DEBUG<< "CalSection: construction. Eta range: "
00049         <<minEta<<" - "<<maxEta<<endreq; 
00050     log <<MSG::DEBUG<< "                          Phi range: "
00051         <<minPhi<<" - "<<maxPhi<<endreq; 
00052     log <<MSG::DEBUG<< " starting granularity Eta: "<<granEta
00053         << " phi: "<<granPhi
00054         <<endreq;
00055 
00056     log <<MSG::DEBUG<<"CalSectionReject parameters:"<<endreq;
00057     log <<MSG::DEBUG<<m_calSectionReject<<endreq;
00058 
00059     int nEta    = int(((maxEta-minEta)/granEta)+0.5);
00060     m_granEta   = (maxEta-minEta)/int(nEta);
00061     m_nPhi      = 1+int(((m_maxPhi-m_minPhi)/granPhi)+0.5);
00062     m_granPhi   = (m_maxPhi-m_minPhi)/double(m_nPhi); 
00063     log <<MSG::DEBUG<< " final granularity Eta: "<<m_granEta
00064         << " phi: "<<m_granPhi
00065         <<endreq; 
00066     log <<MSG::DEBUG<< " final Bins in Eta: "<<nEta
00067         << " phi: "<<m_nPhi
00068         <<endreq; 
00069 
00070     for (int ieta=0; ieta<nEta; ieta++){
00071       for (int iphi=0; iphi<m_nPhi; iphi++){
00072         
00073         double eta        = minEta+((double(ieta)+0.5)*double(m_granEta));
00074         double phi        = m_minPhi + ((double(iphi)+0.5)*double(m_granPhi));
00075         CellDescriptor id = CellDescriptor(eta,phi);
00076         ITwoCptCell* cell = new TwoCptCell(id);
00077         
00078         m_cells.push_back(cell);
00079       }
00080     }
00081   }
  CalSection::~CalSection(){

Atlfast::CalSection::CalSection

(

const CalSection &

c

)

Definition at line 87 of file CalSection.cxx.

00087                                            :
00088     m_log(c.m_log),
00089     m_minEta(c.m_minEta),
00090     m_maxEta(c.m_maxEta),
00091     m_minPhi(c.m_minPhi),
00092     m_maxPhi(c.m_maxPhi),
00093     m_granEta(c.m_granEta),
00094     m_granPhi(c.m_granPhi),
00095     m_nPhi(c.m_nPhi),
00096     m_calSectionReject(c.m_calSectionReject){
00097     std::vector<ITwoCptCell*>::const_iterator i=c.m_cells.begin();
00098     for(;i<c.m_cells.end();++i) m_cells.push_back((*i)->cloneITCC());
00099   }
  //----------------------------------------------------------------

Atlfast::CalSection::~CalSection

(

)

Definition at line 82 of file CalSection.cxx.

                         {
    m_log<<MSG::DEBUG<<"CalSection destructor starts "<<endreq;
    std::for_each(m_cells.begin(), m_cells.end(), DeleteObject());
    m_log<<MSG::DEBUG<<"CalSection destructor ends "<<endreq;
  }

---

Member Function Documentation

void Atlfast::CalSection::reset

(

)

Resets the contents of the calorimeter section.

Definition at line 312 of file CalSection.cxx.

                        {
    //reset the cells of this CalSection

    std::map<int, ITwoCptCell*>::const_iterator i = m_hitCells.begin();
    for(;i!=m_hitCells.end();++i) (*i).second->resetCell();
    m_hitCells.erase(m_hitCells.begin(), m_hitCells.end());
  }

void Atlfast::CalSection::deposit	(	ITransportedParticleCollectionIter &	,
		ITransportedParticleCollectionIter &
	)

Takes the ITransportedParticle (transported through mag field) and deposits it within the calorimeter section.

Definition at line 101 of file CalSection.cxx.

                                                                 {
    // remove particles not in acceptance
    ITransportedParticleCollectionIter divider;
    divider = std::partition(f, l, m_calSectionReject);
    
    // fill cells 
    for(ITransportedParticleCollectionCIter i=divider; i<l; ++i){newHit(*i);}

    // assume sections do not overlap: do not reprocess.
    l=divider;
  }

void Atlfast::CalSection::deposit	(	std::vector< Atlfast::EPileupDeposit * >::iterator &	,
		std::vector< Atlfast::EPileupDeposit * >::iterator &
	)

Takes the Energy pileup map and deposits it within the calorimeter section.

Definition at line 255 of file CalSection.cxx.

                                                                        {
    // remove pileup deposits not in acceptance
    std::vector<Atlfast::EPileupDeposit*>::iterator divider;
    divider = std::partition(f, l, m_calSectionReject);

    // fill cells 
    std::vector<Atlfast::EPileupDeposit*>::iterator  i=divider;
    for(; i<l; ++i) newHit(*i);

    // assume sections do not overlap: do not reprocess.
    // hits processed here: set last to end of unprocessed hits
    l=divider;
  }

void Atlfast::CalSection::depositEcal	(	std::vector< const GridletElement * >::iterator &	,
		std::vector< const GridletElement * >::iterator &
	)

a deposit method into Ecal comaprtment (used with FastShower)

Definition at line 115 of file CalSection.cxx.

                           {
    //    cerr<<"Calsection depositEcal begin"<<endl;
  
    //    cerr<<"  no of elements on entry "<<l-f<<endl;
    std::vector<const GridletElement*>::iterator divider;
    divider = std::partition(f, l, m_calSectionReject);
    //    cerr<<"Calsection depositEcal after partition"<<endl;
  
    std::vector<const GridletElement*>::iterator i;
    for( i=divider; i<l; ++i){newEHit(*i);}
    //    cerr<<"Calsection depositEcal after newhit call"<<endl;
  

    // assume sections do not overlap: do not reprocess.
    l=divider;
    //    cerr<<"  no of elements on exit "<<l-f<<endl;
    //    cerr<<"Calsection depositEcal end"<<endl;
  }

void Atlfast::CalSection::depositHcal	(	std::vector< const GridletElement * >::iterator &	,
		std::vector< const GridletElement * >::iterator &
	)

a deposit method into Hcal comaprtment (used with FastShower)

Definition at line 136 of file CalSection.cxx.

                                {
    
    std::vector<const GridletElement*>::iterator divider;
    divider = std::partition(f, l, m_calSectionReject);
  
    std::vector<const GridletElement*>::iterator i;
    for( i=divider; i<l; ++i){newHHit(*i);}

    // assume sections do not overlap: do not reprocess.
    l=divider;
  }

void Atlfast::CalSection::depositEgen	(	std::vector< const Gridlet * >::iterator &	,
		std::vector< const Gridlet * >::iterator &
	)

set Generated Energy (used with FastShower)

Definition at line 150 of file CalSection.cxx.

                                {
    std::vector<const Gridlet*>::iterator divider;
    divider = std::partition(f, l, m_calSectionReject);
    
    
    std::vector<const Gridlet*>::iterator i;
    for( i=divider; i<l; ++i){setEgen(*i);}
    
    // assume sections do not overlap: do not reprocess.
    l=divider;
  }

void Atlfast::CalSection::giveHits	(	const ICellSelector *	,
		ITwoCptCellCollection *
	)			const

Returns a list of hit cells within the calorimeter section.

Definition at line 285 of file CalSection.cxx.

                                                                           {
    std::map<int, ITwoCptCell*, std::less<int> >::const_iterator i;
    for(i=m_hitCells.begin(); i!=m_hitCells.end(); ++i){
      //test cell before newing and filling vector
      if( (*p_cellselect) (((*i).second))){ 
        ITwoCptCell* cell = ((*i).second)->cloneITCC();
      //Need to copy cells because Athena deletes container AND POINTERS
      // AND POINTEES
        cells->push_back(cell);
      }
    }

  }

void Atlfast::CalSection::smearCells

(

ISmearer *

)

Smear the cells resolution

Definition at line 300 of file CalSection.cxx.

                          {
    std::map<int, ITwoCptCell*, std::less<int> >::const_iterator i;
    for(i=m_hitCells.begin(); i!=m_hitCells.end(); ++i){
      if((*((*i).second)).momentum().e() > 0){
        HepLorentzVector temp = p_smearer->smear((*((*i).second)).momentum());
          (*((*i).second)).setPt(temp);
      }
    }
  }

void Atlfast::CalSection::newHit

(

const ITransportedParticle *

)

[private]

Registers a new hit in the calorimeter section.

Definition at line 165 of file CalSection.cxx.

                                                        {
    //check cell is reasonable
    const HepMC::GenParticle* particle = tp->particle();
    //const double eta = particle->momentum().pseudoRapidity();
    const double eta = tp->eta();
    const double phi = tp->phi();
    int ind = iindex(phi, eta);
    ITwoCptCell* cell = m_cells[ind];
    assert(std::abs( eta-(cell->eta()) ) <m_granEta);
    assert(std::abs( phi-(cell->phi()) ) <m_granPhi);
    cell->newHit(particle);
    //use a map to store pointers to hit cells to avoid
    //duplicate entries if the cell is hit more than once.
    m_hitCells[ind] = cell;

    return;
  }

void Atlfast::CalSection::newHit

(

const EPileupDeposit *

)

[private]

Registers a new energy pileup in the calorimeter section.

Definition at line 271 of file CalSection.cxx.

                                                         {
    //check cell is reasonable
    int ind = index(ep);
    ITwoCptCell* cell = m_cells[ind];
    assert(std::abs( (ep->eta()) - (cell->eta()) )<m_granEta);
    assert(std::abs( (ep->phi()) - (cell->phi()) )<m_granPhi);
    cell->newHit(ep);
     //use a map to store pointers to hit cells to avoid
    //duplicate entries if the cell is hit more than once.
    m_hitCells[ind] = cell;
    return;
  }

void Atlfast::CalSection::newEHit

(

const GridletElement *

)

[private]

Registers a new deposit in the Ecal in the calorimeter section.

Definition at line 183 of file CalSection.cxx.

                                                   {
    int ind = index(ge);
    ITwoCptCell* cell = m_cells[ind];
    //    cerr<<"gridletEl, phi, eta, eT:     "
    //    <<ge->phi()<<" "<<ge->eta()<<" "<<ge->et()<<endl;
    //    cerr<<"cell, phi, eta, granularity: "<<cell->phi()<<" "<<cell->eta()
    //  <<" "<<m_granPhi<<" "<<m_granEta<<endl;
    assert(std::abs( (ge->eta()-cell->eta()) )<m_granEta);
    assert(std::abs( (ge->phi()-cell->phi()) )<m_granPhi);
    cell->depositEcal(ge->et());
    //use a map to store pointers to hit cells to avoid
    //duplicate entries if the cell is hit more than once.
    m_hitCells[ind] = cell;

    return;
  }

void Atlfast::CalSection::newHHit

(

const GridletElement *

)

[private]

Registers a new deposit in the Ecal in the calorimeter section.

Definition at line 200 of file CalSection.cxx.

                                                   {
    int ind = index(ge);
    ITwoCptCell* cell = m_cells[ind];

    assert(std::abs( (ge->eta())-(cell->eta()) )<m_granEta);
    assert(std::abs( (ge->phi())-(cell->phi()) )<m_granPhi);

    cell->depositHcal(ge->et());
    //use a map to store pointers to hit cells to avoid
    //duplicate entries if the cell is hit more than once.
    m_hitCells[ind] = cell;

    return;
  }

void Atlfast::CalSection::setEgen

(

const Gridlet *

)

[private]

find the cell index for a variety of input arguments

Definition at line 215 of file CalSection.cxx.

                                            {
    int ind = index(gr);
    ITwoCptCell* cell = m_cells[ind];
    cell->addEgen(gr->eGen());
    m_hitCells[ind] = cell;
    return;
  }

int Atlfast::CalSection::index

(

const Gridlet *

)

const [private]

find the cell index for a variety of input arguments

Definition at line 223 of file CalSection.cxx.

                                               {
    int ind = this->iindex( gr->phi0(), gr->eta0() );
    return ind;
  }

int Atlfast::CalSection::index

(

const GridletElement *

)

const [private]

Definition at line 228 of file CalSection.cxx.

                                                      {
    int ind = this->iindex( ge->phi(), ge->eta() );
    return ind;
  }

int Atlfast::CalSection::index

(

const ITransportedParticle *

)

const [private]

Definition at line 233 of file CalSection.cxx.

                                                             {
    const HepMC::GenParticle* particle = tp->particle();
    const double eta = particle->momentum().pseudoRapidity();
    const double phi = tp->phi();
    int ind = iindex(phi, eta);
    return ind;
  }

int Atlfast::CalSection::index

(

const Atlfast::EPileupDeposit *

)

const [private]

Definition at line 241 of file CalSection.cxx.

                                                             {
    int ind = iindex(ep->phi(),ep->eta());
    return ind;
  }

int Atlfast::CalSection::iindex	(	double	phi,
		double	eta
	)			const [private]

Definition at line 246 of file CalSection.cxx.

                                                    {
    int iphi = int( (phi-m_minPhi)/m_granPhi );
    int ieta = int( (eta-m_minEta)/m_granEta );
    
    return  (m_nPhi*ieta  + iphi);
  }

---

Member Data Documentation

MsgStream& Atlfast::CalSection::m_log [private]

Definition at line 94 of file CalSection.h.

const double Atlfast::CalSection::m_minEta [private]

minimum Eta sensitivity

Definition at line 112 of file CalSection.h.

const double Atlfast::CalSection::m_maxEta [private]

maximum Eta sensitivity

Definition at line 114 of file CalSection.h.

const double Atlfast::CalSection::m_minPhi [private]

minimum phi sensitivity

Definition at line 116 of file CalSection.h.

const double Atlfast::CalSection::m_maxPhi [private]

maximum phi sensitivity

Definition at line 118 of file CalSection.h.

double Atlfast::CalSection::m_granEta [private]

granularity in eta

Definition at line 120 of file CalSection.h.

double Atlfast::CalSection::m_granPhi [private]

granularity in phi

Definition at line 122 of file CalSection.h.

int Atlfast::CalSection::m_nPhi [private]

number of cells in phi space

Definition at line 124 of file CalSection.h.

CalSectionReject Atlfast::CalSection::m_calSectionReject [private]

Calorimeter selection object.

Definition at line 128 of file CalSection.h.

std::vector<ITwoCptCell*> Atlfast::CalSection::m_cells [private]

Cells owned by calorimeter section.

Definition at line 132 of file CalSection.h.

std::map<int, ITwoCptCell*, std::less<int> > Atlfast::CalSection::m_hitCells [private]

Definition at line 133 of file CalSection.h.

---

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

• /afs/cern.ch/user/s/sdean/project_builds/AtlasOffline-rel_1/Simulation/Atlfast/AtlfastAlgs/AtlfastAlgs/CalSection.h
• /afs/cern.ch/user/s/sdean/project_builds/AtlasOffline-rel_1/Simulation/Atlfast/AtlfastAlgs/src/CalSection.cxx

---