// $Id$
//
//    File: DBCALGeometry.h
// Created: Thu Nov 17 15:10:51 CST 2005
// Creator: gluexuser (on Linux hydra.phys.uregina.ca 2.4.20-8smp i686)
//

#ifndef _DBCALGeometry_
#define _DBCALGeometry_

#include <JANA/JObject.h>
#include <JANA/JFactory.h>
using namespace jana;

// create a single number channel id which is useful in algorithms
// if M L S are module layer sector the bit map looks like:
// MMMM MMMM LLLL SSSS

#define MODULE_SHIFT 8
#define MODULE_MASK 0xFF00
#define LAYER_SHIFT 4
#define LAYER_MASK 0x00F0
#define SECTOR_SHIFT 0
#define SECTOR_MASK 0X000F

class DBCALGeometry : public JObject{
    
public:
    
    JOBJECT_PUBLIC(DBCALGeometry);
    
    DBCALGeometry();
    
    enum End { kUpstream, kDownstream };
        
    static int NBCALMODS;         ///> number of modules
    static int NBCALLAYS1;        ///> number of layers in first 10 ccm 
    static int NBCALLAYS2;        ///> number of layers in last  15 cm 
    static int NBCALSECS1;        ///> number of sectors in first 10cm of Mod 
    static int NBCALSECS2;        ///> number of sectors in last 15cm of Mod 
    static float BCALINNERRAD;    ///> innner radius of BCAL in cm
    static int BCALMID;         ///> first outer layer (default 7)
    static float m_radius[11];
    float BCALMIDRAD;      ///> mid radius of BCAL in cm
    static float BCALOUTERRAD;    ///> outer radius of BCAL in cm
    static float BCALFIBERLENGTH; ///> BCAL Scintilator fiber lenth in cm
    static float GLOBAL_CENTER;  ///> center of BCAL in gloobal coordinate system
    
    static float ATTEN_LENGTH;   ///> attenuation length
    static float C_EFFECTIVE;    ///> speed of light in fibers 
    
    static inline int module( int cellId ) { return ( cellId & MODULE_MASK ) >> MODULE_SHIFT; };
    static inline int layer( int cellId ) { return ( cellId & LAYER_MASK ) >> LAYER_SHIFT; };
    static inline int sector( int cellId ) { return ( cellId & SECTOR_MASK ) >> SECTOR_SHIFT; };                
    
    static inline int cellId( int module, int layer, int sector ) {
        return ( ( module << MODULE_SHIFT ) | ( layer << LAYER_SHIFT ) | 
                 ( sector << SECTOR_SHIFT ) ); }

};

#endif // _DBCALGeometry_