---

Material.h

Go to the documentation of this file.

// $Id: Material.h,v 1.3 2001/08/08 05:49:41 chehtman Exp $
// File: Material.h
//
// This file is part of Gismo 2
//
// Contents ----------------------------------------------------------------
//
//      class Material
//
// Description
//

// End --------------------------------------------------------------------

#ifndef __MATERIAL_H
#define __MATERIAL_H


#include <string>
#include <vector>
#include <list>

#include "MatData.h"

class Material
{
public:
    Material()
        : zeff(0), aeff(0), rho(0), rlc(0), lambda(0)
    {};
    

    Material(std::istream&, const char* name);
    
    
    ~Material();
    

    static Material* hatch(const char* name);
    
    int isVacuum()const {return rho==0;}
    
    float dedx(float beta, float eta) const;
    
    struct AtomicElement;
    
    
    AtomicElement* addElementData(char[],float,float,float,float);

    AtomicElement* addElement(const char* bufffer);
    
    const char* name() const {return isVacuum()? "vacuum": materialName;}
    
    const char * title() const {return materialName;}
    
    float Z() const;
    float A() const;
    float density() const;
    
    float radiationLength() const {return rlc;}
    float interactionLength() const { return lambda; }
    
    void setRadiationLength(float r){rlc=r;}
    void setInteractionLength(float r){lambda=r;}
    
    static Material* vacuum;
    
    static void dataPath(const char* path);
    
    static void addPath(const char* path);
    
    static const char *getDataPath();
    


    static unsigned declareMatData(MatData*);
    
    void* getMatData(unsigned index);
    
    typedef std::list<AtomicElement*> Element_list;     
    Element_list elementList;
    
    static void printAll(std::ostream&);
    
    void printOn(std::ostream&);
    
    
    static void deleteAll();
    
    struct AtomicElement 
    {
        AtomicElement(char*,float,float,float,float);
        
        char name[3];
        // Name of the element
        
        float z,a,pz,rhoz;
        // values for Z, A, fraction, and default density
        
        /*  */ operator int()const{return static_cast<int>(z);}
    };
    
    
    protected:
        
        static std::string* s_dataPath;
        // data member pointing to the .mat directory
        
        std::string m_fileName;
        // name of the file
        
        char materialName[16];
        // name of this stuff from data card
        
        float zeff;
        // effective Z of material
        
        float aeff;
        // effective A of material
        
        float rho;
        // density in gm/cm**3
        
        float rlc;
        // radiation length
        
        float lambda;
        // interaction length
        
        typedef std::vector<MatData*> MatData_list; 
        
        MatData_list matDataList;
        // list of pointers to auxialiary data
        
    public:
        Material * nextMaterial;
        static  Material * firstMaterial;
        // implement a linked list of all defined materials
        
    private:
        float energyLoss1,energyLoss2,energyLoss3,IP;
        // constants used by dedx

        float density_correction(float eta) const;


};
inline void*
Material::getMatData(unsigned index)
{
    return matDataList[index];
}
inline float
Material::Z()const{  return zeff;}

inline float
Material::A()const{  return aeff;}

inline float
Material::density()const{  return rho;}

#endif

---