//    File: DBCALGeometry.h
// Created: Fri Nov 26 15:10:51 CST 2010
// Creator: dwbennet

#include "DBCALGeometry.h"

int DBCALGeometry::NBCALMODS = 48;
int DBCALGeometry::NBCALLAYS1 =  2;
int DBCALGeometry::NBCALLAYS2 =  1; 
int DBCALGeometry::NBCALSECS1 =  4; 
int DBCALGeometry::NBCALSECS2 =  4;
int DBCALGeometry::BCALMID = 7; //Enter the first layer in "outer cells" (default 7)
float DBCALGeometry::BCALINNERRAD = 64.3;   
float DBCALGeometry::BCALOUTERRAD = 86.17;
float DBCALGeometry::BCALFIBERLENGTH = 390.0;
float DBCALGeometry::GLOBAL_CENTER = 212;
  
float DBCALGeometry::ATTEN_LENGTH = 300.;
   
float DBCALGeometry::C_EFFECTIVE  = 16.75;

float DBCALGeometry::m_radius[] = { 64.3, 
				  66.3,
				  68.3,
				  70.3,
				  72.3,
				  74.3,
				  76.3,
				  78.77,
				  81.24,
				  83.70,
				  86.17};

DBCALGeometry::DBCALGeometry()  
{
  BCALMIDRAD = m_radius[BCALMID-1];
  /// End if groupings do not evenly divide SiPM cells
  if(  (BCALMID-1) % NBCALLAYS1 != 0 
       || (11-BCALMID) % NBCALLAYS2 != 0
       ||  4 % NBCALSECS1 != 0
       ||  4 % NBCALSECS2 != 0)
    {
      std::cout<<"ERROR: Bad BCAL fADC groupings, do not evenly divide cells";
      assert (false);
    }
}

int
DBCALGeometry::module( int cellId ) {
  
  return ( cellId & MODULE_MASK ) >> MODULE_SHIFT;
}

int
DBCALGeometry::layer( int cellId ) {
  
  return ( cellId & LAYER_MASK ) >> LAYER_SHIFT;
}

int
DBCALGeometry::sector( int cellId ) {
  
  return ( cellId & SECTOR_MASK ) >> SECTOR_SHIFT;
}

int
DBCALGeometry::cellId( int module, int layer, int sector ) {
  
  return ( ( module << MODULE_SHIFT ) | 
           ( layer << LAYER_SHIFT   ) | 
           ( sector << SECTOR_SHIFT ) );
}

float
DBCALGeometry::r( int cell ) {
 
  int lay = layer( cell );
  
  float innerRad, outerRad;
 
  if( lay <= NBCALLAYS1 ) { // we're in the inner region
   
    // the size of the layer in SiPM units
    int laySize = ( BCALMID - 1 ) / NBCALLAYS1;
    
    innerRad = m_radius[ laySize * ( lay - 1 ) ];
    
    // index below is same as:  laySize * ( lay - 1 ) + laySize
    outerRad = m_radius[ laySize * lay ];
  }
  else{
    
    // find the size of the outer layer in SiPM units
    int laySize = ( 11 - BCALMID ) / NBCALLAYS2;
    
    // subtract the number of inner layers to get
    // layer number within the outer section
    lay -= NBCALLAYS1;
    
    innerRad = m_radius[ BCALMID - 1 + laySize * ( lay - 1 ) ];
    outerRad = m_radius[ BCALMID - 1 + laySize * lay ];
  }
  
  return 0.5 * ( innerRad + outerRad );
}

float
DBCALGeometry::rSize( int cell ) {
  
  int lay = layer( cell );
  
  float innerRad, outerRad;
  
  if( lay <= NBCALLAYS1 ) { // in the inner region
    
    // the size of the layer in SiPM units
    int laySize = ( BCALMID - 1 ) / NBCALLAYS1;
    
    innerRad = m_radius[ laySize * ( lay - 1 ) ];
    
    // index below is same as:  laySize * ( lay - 1 ) + laySize
    outerRad = m_radius[ laySize * lay ];
  }
  else{ // in the outer region
    
    // find the size of the outer layer in SiPM units
    int laySize = ( 11 - BCALMID ) / NBCALLAYS2;
    
    // subtract the number of inner layers to get
    // layer number within the outer section
    lay -= NBCALLAYS1;
    
    innerRad = m_radius[ BCALMID - 1 + laySize * ( lay - 1 ) ];
    outerRad = m_radius[ BCALMID - 1 + laySize * lay ];
  }
  
  return ( outerRad - innerRad );
}

float
DBCALGeometry::phi( int cell ) {
 
  int lay = layer( cell );
  int sect = sector( cell );
  
  float sectSize;
  
  if( lay <= NBCALLAYS1 ) { 
    
    sect += NBCALSECS1 * ( module( cell ) - 1 );
    sectSize = 2 * PI / ( NBCALMODS * NBCALSECS1 );
  }
  else {
    
    sect += NBCALSECS2 * ( module( cell ) - 1 );
    sectSize = 2 * PI / ( NBCALMODS * NBCALSECS2 );
  }
  
  return sectSize * ( sect - 0.5 );
}

float
DBCALGeometry::phiSize( int cell ) {
  
  int lay = layer( cell );
  
  float sectSize;
  
  if( lay <= NBCALLAYS1 ) { 
    
    sectSize = 2 * 3.1416 / ( NBCALMODS * NBCALSECS1 );
  }
  else {
    
    sectSize = 2 * 3.1416 / ( NBCALMODS * NBCALSECS2 );
  }
  
  return sectSize;
}