---

TowerArray.h

Go to the documentation of this file.

// $Id: TowerArray.h,v 1.3 2001/01/23 03:41:32 burnett Exp $
//
// This file is part of Gismo 2
//
//
#ifndef __TOWERARRAY_H
#define __TOWERARRAY_H

#include "gismo/DetectorFrame.h"
#include <vector>
class Volume;
class Point;
class MCParticle;

class TowerElement
{
public:
    friend class TowerArray;
    TowerElement(long i=0, float e=0)
        : m_index(i), m_energy(e)
    {}
    
    bool operator<(const TowerElement& rhs)const{return m_index < rhs.m_index;}
    bool operator==(const TowerElement&rhs)const{return m_index == rhs.m_index;}
    float operator+=(float e){return m_energy+= e;}
    ~TowerElement();
    
    float energy()const{return m_energy;}
    long index()const{return m_index;}
    
private:
    long  m_index;
    // unique index reflecting position (TowerArray knows)
    
    float m_energy, simEnergy;
    // accumulated deposited energy, smeared
};

typedef std::vector<TowerElement > TowerElementList;


class TowerArray : public DetectorFrame , public TowerElementList
{
    
public:
    TowerArray( int nx, int ny, float lx, float ly, float z, float eR)
        : xSegmentation(nx), ySegmentation(ny)
        , xWidth(lx) , yWidth(ly)
        , zPlane(z)
        , energyResolution(eR)
    {}
    
    virtual ~TowerArray();
    
    virtual void score(MCParticle* );
    
    float scoreTower(long, float);
    
    virtual void clear();
    
    virtual void draw(gui::DisplayRep&)const;
    
    int hits()const;
    Point xtalPos(unsigned i)const;
    long xtalId(unsigned i)const;
    float energy(unsigned i)const;
    
    virtual void readData(std::istream&);
    virtual void writeData(std::ostream&);
    
    virtual void writeParameters(std::ostream&)const;
    
    virtual void printOn(std::ostream&)const;
    
    const char* nameOf()const {return "TowerArray";}
    
    virtual void generateResponse();
    Point *meanEnergy(float, float *, int *);
    
    float displayScale()const {return s_display_scale;}
    static float s_display_scale;
    
protected:
    
    int xSegmentation,
        ySegmentation;
    // the number of elements in the respective directions
    
    float xWidth,
        yWidth;
    // total size of the array
    
    float zPlane;
    // position of plane where boxes relecting response should be displayed
    
    float energyResolution;
    // root(E) - Gaussian energy resolution
    
    
    virtual long calculateIndex(const Point&)const;
    // return an appropriate index to identify the element
    
    virtual void calculateBin(long,float&,float&)const;
    // determine bin corner from index
    
    TowerElement* towers(int i){return begin()+i;}
    const TowerElement* towers(int i)const{return begin()+i;}
    
    
};




#endif

---