MEEBGeom.cc

Go to the documentation of this file.

#include <iostream>
#include <stdlib.h>
#include <string>
#include <assert.h>
using namespace std;

#include "CalibCalorimetry/EcalLaserAnalyzer/interface/ME.h"
#include "CalibCalorimetry/EcalLaserAnalyzer/interface/MEEBGeom.h"

#include <TGraph.h>

//GHM ClassImp(MEEBGeom)

int 
MEEBGeom::barrel( EBGlobalCoord ieta, EBGlobalCoord iphi )
{
  int iz=1;
  if( ieta<0 ) iz=-1;
  ieta *= iz;
  assert( ieta>0 && ieta<=85 );
  if( iphi<0 ) iphi+=360;
  assert( iphi>0 && iphi<=360 );
  return iz;
}

int 
MEEBGeom::sm( EBGlobalCoord ieta, EBGlobalCoord iphi )
{
  int iz=1;
  if( ieta<0 ) iz=-1;
  ieta *= iz;
  assert( ieta>0 && ieta<=85 );

  //  int iphi_ = iphi+10;
  int iphi_ = iphi;
  if( iphi_>360 ) iphi_-=360;
  assert( iphi_>0 && iphi_<=360 );

  int ism = (iphi_-1)/20 + 1;
  assert( ism>=1 && ism<=18 );
  //  if( iz==1 ) ism += 18;
  if( iz==-1 ) ism += 18;

  return ism;
}

int 
MEEBGeom::dcc( EBGlobalCoord ieta, EBGlobalCoord iphi )
{
  int ism = sm( ieta, iphi );
  return dccFromSm( ism );
}

int
MEEBGeom::dccFromSm( int ism )
{
  assert( ism>=1 && ism<=36 );
  int iz=1;
  if( ism>18 ) iz=-1;
  if( iz==-1  ) ism-=18;
  assert( ism>=1 && ism<=18 );
  int idcc = 9+ism;
  if( iz==+1  ) idcc+=18;
  return idcc;
}

int
MEEBGeom::smFromDcc( int idcc )
{
  if( idcc>600 ) idcc-=600;  // also works with FEDids
  assert( idcc>=10 && idcc<=45 );
  int ism=idcc-9;
  if( ism>18 ) ism-=18;
  else         ism+=18;
  return ism;
}

TString
MEEBGeom::smName( int ism )
{
  assert( ism>=1 && ism<=36 );
  TString out = "EB+";
  if( ism>18 )
    {
      out = "EB-";
      ism -= 18;
    }
  out += ism;
  return out;
}

int 
MEEBGeom::lmmod( EBGlobalCoord ieta, EBGlobalCoord iphi )
{
  std::pair<EBLocalCoord,EBLocalCoord> ixy = localCoord( ieta, iphi );
  return lm_channel( ixy.first/5, ixy.second/5 );
}

int 
MEEBGeom::tt( EBGlobalCoord ieta, EBGlobalCoord iphi )
{
  std::pair<EBLocalCoord,EBLocalCoord> ixy = localCoord( ieta, iphi );
  return tt_channel( ixy.first/5, ixy.second/5 );
}

int 
MEEBGeom::crystal( EBGlobalCoord ieta, EBGlobalCoord iphi )
{
  std::pair<EBLocalCoord,EBLocalCoord> ixy = localCoord( ieta, iphi );
  return crystal_channel( ixy.first, ixy.second );
}

int 
MEEBGeom::side( EBGlobalCoord ieta, EBGlobalCoord iphi )
{
  int ilmmod = lmmod( ieta, iphi );
  return (ilmmod%2==0)?1:0;
}

int 
MEEBGeom::lmr( EBGlobalCoord ieta, EBGlobalCoord iphi )
{
  int idcc  = dcc( ieta, iphi );
  int ism   = idcc-9;
  int iside = side( ieta, iphi );
  int ilmr = 1 + 2*(ism-1) + iside;
  return ilmr;
}

std::pair< MEEBGeom::EBLocalCoord, MEEBGeom::EBLocalCoord >    
MEEBGeom::localCoord( MEEBGeom::EBGlobalCoord ieta, MEEBGeom::EBGlobalCoord iphi )
{
  int iz=1;
  if( ieta<0 ) iz=-1;
  ieta *= iz;
  assert( ieta>0 && ieta<=85 );

  //  int iphi_ = iphi+10;
  int iphi_ = iphi;
  if( iphi_>360 ) iphi_-=360;
  assert( iphi_>0 && iphi_<=360 );

  int ix = ieta-1;
  
  int iy = (iphi_-1)%20;
  if( iz==-1 ) iy = 19-iy;
  // if( iz==1 ) iy = 19-iy;

  return std::pair< EBLocalCoord, EBLocalCoord >(ix,iy);  
}

std::pair< MEEBGeom::EBLocalCoord, MEEBGeom::EBLocalCoord >    
MEEBGeom::localCoord( int icr )
{
  assert( icr>=1 && icr<=1700 );
  int ix = (icr-1)/20;
  int iy = 19 - (icr-1)%20;
  return std::pair< EBLocalCoord, EBLocalCoord >(ix,iy);  
}

std::pair< MEEBGeom::EBGlobalCoord, MEEBGeom::EBGlobalCoord >    
MEEBGeom::globalCoord( int ism, MEEBGeom::EBLocalCoord ix, MEEBGeom::EBLocalCoord iy )
{
  assert( ism>=1 && ism<=36 );
  assert( ix>=0 && ix<85 );
  assert( iy>=0 && iy<20 );
  //  int iz=-1;
  int iz=1;
  if( ism>18 ) 
    {
      iz=-1;
      ism -= 18;
    }
  if( iz==-1 ) iy = 19-iy;
  // if( iz==1 ) iy = 19-iy;

  int ieta = ix+1;
  ieta *= iz;
  //  int iphi = -9 + iy + 20*(ism-1);
  int iphi = 1 + iy + 20*(ism-1);

  return std::pair< EBGlobalCoord, EBGlobalCoord >(ieta,iphi);  
}

std::pair< float, float >    
MEEBGeom::globalCoord( int ism, float x, float y )
{
  assert( ism>=1 && ism<=36 );
  int iz=1;
  if( ism>18 ) 
    {
      iz=-1;
      ism -= 18;
    }
  if( iz==-1 ) y = 19-y;

  float eta = x+1;
  eta *= iz;
  //  float phi = -9 + y + 20*(ism-1);
  float phi = 1 + y + 20*(ism-1);

  return std::pair< float, float >(eta,phi);  
}

std::pair< MEEBGeom::EBGlobalCoord, MEEBGeom::EBGlobalCoord >    
MEEBGeom::globalCoord( int ism, int icr )
{
  assert( ism>=1 && ism<=36 );
  assert( icr>=1 && icr<=1700 );

  int ix = (icr-1)/20;
  int iy = 19 - (icr-1)%20;

  return globalCoord( ism, ix, iy );  
}

int 
MEEBGeom::lm_channel( EBTTLocalCoord iX, EBTTLocalCoord iY )
{
  static const int 
    idx_[] = 
    {
     // 0  1  2  3  4  5  6  7  8  9 10 11 12 13 14 15 16
        1, 2, 2, 2, 2, 4, 4, 4, 4, 6, 6, 6, 6, 8, 8, 8, 8, // 3
        1, 2, 2, 2, 2, 4, 4, 4, 4, 6, 6, 6, 6, 8, 8, 8, 8, // 2
        1, 3, 3, 3, 3, 5, 5, 5, 5, 7, 7, 7, 7, 9, 9, 9, 9, // 1
        1, 3, 3, 3, 3, 5, 5, 5, 5, 7, 7, 7, 7, 9, 9, 9, 9  // 0
     // 0  1  2  3  4  5  6  7  8  9 10 11 12 13 14 15 16
     };

  int iym, ixm, il, ic, ii;
  iym=4;
  ixm=17;
  int iX_ = iX+1;
  int iY_ = iY+1;
  il=iym-iY_;
  ic=iX_-1;
  ii=il*ixm+ic;
  if(ii < 0 || ii > (int)(sizeof(idx_)/sizeof(int))) { return -1; };
  return idx_[ii];
}

int 
MEEBGeom::tt_type( EBTTLocalCoord iX, EBTTLocalCoord iY )
{
  static const int 
    idx_[] =
    {
     // 0  1  2  3  4  5  6  7  8  9 10 11 12 13 14 15 16
        1, 1, 1, 2, 2, 1, 1, 2, 2, 1, 1, 2, 2, 1, 1, 2, 2, // 3
        1, 1, 1, 2, 2, 1, 1, 2, 2, 1, 1, 2, 2, 1, 1, 2, 2, // 2
        1, 1, 1, 2, 2, 1, 1, 2, 2, 1, 1, 2, 2, 1, 1, 2, 2, // 1
        1, 1, 1, 2, 2, 1, 1, 2, 2, 1, 1, 2, 2, 1, 1, 2, 2  // 0
     // 0  1  2  3  4  5  6  7  8  9 10 11 12 13 14 15 16
     };

  int iym, ixm, il, ic, ii;
  iym=4;
  ixm=17;
  int iX_ = iX+1;
  int iY_ = iY+1;
  il=iym-iY_;
  ic=iX_-1;
  ii=il*ixm+ic;
  if(ii < 0 || ii > (int)(sizeof(idx_)/sizeof(int))) { return -1; };
  return idx_[ii];
}

int
MEEBGeom::hv_channel( EBTTLocalCoord iX, EBTTLocalCoord iY )
{
  static const int 
    idx_[] = 
    {
     // 0  1  2  3   4   5   6   7   8   9  10  11  12  13  14  15  16
        1, 3, 5, 7,  9, 11, 13, 15, 17, 19, 21, 23, 25, 27, 29, 31, 33, // 3
        1, 3, 5, 7,  9, 11, 13, 15, 17, 19, 21, 23, 25, 27, 29, 31, 33, // 2
        2, 4, 6, 8, 10, 12, 14, 16, 18, 20, 22, 24, 26, 28, 30, 32, 34, // 1
        2, 4, 6, 8, 10, 12, 14, 16, 18, 20, 22, 24, 26, 28, 30, 32, 34  // 0
     // 0  1  2  3   4   5   6   7   8   9  10  11  12  13  14  15  16
     };

  int iym, ixm, il, ic, ii;
  iym=4;
  ixm=17;
  int iX_ = iX+1;
  int iY_ = iY+1;
  il=iym-iY_;
  ic=iX_-1;
  ii=il*ixm+ic;
  if(ii < 0 || ii > (int)(sizeof(idx_)/sizeof(int))) { return -1; };
  return idx_[ii];
}

int
MEEBGeom::lv_channel( EBTTLocalCoord iX, EBTTLocalCoord iY )
{
  static const int
    idx_[] =
    {
     // 0  1  2  3  4  5  6  7  8   9  10  11  12  13  14  15  16
        1, 2, 2, 4, 4, 6, 6, 8, 8, 10, 10, 12, 12, 14, 14, 16, 16, // 3
        1, 2, 2, 4, 4, 6, 6, 8, 8, 10, 10, 12, 12, 14, 14, 16, 16, // 2
        1, 3, 3, 5, 5, 7, 7, 9, 9, 11, 11, 13, 13, 15, 15, 17, 17, // 1
        1, 3, 3, 5, 5, 7, 7, 9, 9, 11, 11, 13, 13, 15, 15, 17, 17  // 0
     // 0  1  2  3  4  5  6  7  8   9  10  11  12  13  14  15  16
     };

  int iym, ixm, il, ic, ii;
  iym=4;
  ixm=17;
  int iX_ = iX+1;
  int iY_ = iY+1;
  il=iym-iY_;
  ic=iX_-1;
  ii=il*ixm+ic;
  if(ii < 0 || ii > (int)(sizeof(idx_)/sizeof(int))) { return -1; };
  return idx_[ii];
}

int 
MEEBGeom::tt_channel( EBTTLocalCoord iX, EBTTLocalCoord iY )
{

  int itt =  4*iX+4-iY;  // :)

  return itt;

}

int 
MEEBGeom::crystal_channel( EBLocalCoord ix, EBLocalCoord iy )
{
  // "Test beam numbering"
  int icr =  20*ix + 19-iy + 1;  // :)

  return icr;
}

int 
MEEBGeom::electronic_channel( EBLocalCoord ix, EBLocalCoord iy )
{
  int iX = ix/5;
  int iY = iy/5;
  int itt  = tt_channel( iX, iY );
  int type = tt_type( iX, iY );

  int iVFE = ix%5+1;
  int islot = iy%5+1;
  if( iVFE%2==1 ) islot = 6-islot;
  int icr = 5*(iVFE-1)+(islot-1);
  // rotate for type-1 towers
  if( type==1 )
    {
      icr = 24-icr;
    }
  icr += 25*(itt-1);

  return icr;
}

TGraph*
MEEBGeom::getGraphBoundary(  int type, int num, bool global )
{
  int ism_=0;
  if( type==iSuperModule )
    {
      ism_ = num; 
    }
  else if( type==iLMRegion )
    {
      ism_ = (num-1)/2+1;
      if( ism_>18 ) 
    ism_-=18;
      else
    ism_+=18;
    }
  else
      abort();

  int ism=1;
  if( global ) ism = ism_;

  //  std::list< std::pair< float, float > > l;
  // getBoundary( l, type, num, global, ism );
  //  int n = l.size();
  //  if( n==0 ) return 0;
  
  float ix[100];
  float iy[100];
  int ixmin =  0;
  int ixmax = 84;
  int iymin =  0;
  int iymax = 19;

  int n(0);
  if( type==iSuperModule )
    {
      n=5;
      ix[0] = ixmin-0.5;
      iy[0] = iymin-0.5;
      ix[1] = ixmin-0.5;
      iy[1] = iymax+0.5;
      ix[2] = ixmax+0.5;
      iy[2] = iymax+0.5;
      ix[3] = ixmax+0.5;
      iy[3] = iymin-0.5;
      ix[4] = ix[0];
      iy[4] = iy[0];
    }
  else if( type==iLMRegion )
    {
      int iside = num;
      if( global )
    {
      iside = (num-1)%2;
    }
      
      if( iside==1 )
    {
      n=5;
      ix[0] = ixmin+5-0.5;
      iy[0] = iymin+10-0.5;
      ix[1] = ixmin+5-0.5;
      iy[1] = iymax+0.5;
      ix[2] = ixmax+0.5;
      iy[2] = iymax+0.5;
      ix[3] = ixmax+0.5;
      iy[3] = iymin+10-0.5;
      ix[4] = ix[0];
      iy[4] = iy[0];
    }
      else
    {
      n=7;
      ix[0] = ixmin-0.5;
      iy[0] = iymin-0.5;
      ix[1] = ixmin-0.5;
      iy[1] = iymax+0.5;
      ix[2] = ixmin+5-0.5;
      iy[2] = iymax+0.5;
      ix[3] = ixmin+5-0.5;
      iy[3] = iymax-10+0.5;
      ix[4] = ixmax+0.5;
      iy[4] = iymax-10+0.5;
      ix[5] = ixmax+0.5;
      iy[5] = iymin-0.5;
      ix[6] = ix[0];
      iy[6] = iy[0];
    }
    }

  if( global )
    {
      for( int ii=0; ii<n; ii ++ )
    {
      std::pair<float,float> xy = globalCoord( ism, ix[ii], iy[ii] ); 
      ix[ii] = xy.first;
      iy[ii] = xy.second;
    }
    }

//   int ii=0;
//   std::list< std::pair< float, float > >::const_iterator l_it;      
//   for( l_it=l.begin(); l_it!=l.end(); l_it++ )
//     {
//       //      std::cout << "[" << l_it->first << "," << l_it->second << "]" << std::endl;
//       ix[ii] = l_it->first;
//       iy[ii] = l_it->second;
//       ii++;
//     }


//  assert( ii==n );
  return new TGraph( n, ix, iy );
}

std::pair< int, int >
MEEBGeom::pn( int ilmmod )
{
  switch( ilmmod )
    {
    case   1: return std::pair<int,int>(  0,  5 );
    case   2: return std::pair<int,int>(  1,  6 );
    case   3: return std::pair<int,int>(  1,  6 );
    case   4: return std::pair<int,int>(  2,  7 );
    case   5: return std::pair<int,int>(  2,  7 );
    case   6: return std::pair<int,int>(  3,  8 );
    case   7: return std::pair<int,int>(  3,  8 );
    case   8: return std::pair<int,int>(  4,  9 );
    case   9: return std::pair<int,int>(  4,  9 );
    default:
      abort();
    }
  return std::pair<int,int>(-1,-1);
}

std::pair<int,int> 
MEEBGeom::memFromLmr( int ilmr )
{
  std::pair< int, int > dccAndSide_ = ME::dccAndSide( ilmr );
  int idcc  = dccAndSide_.first;
  return std::pair<int,int>( idcc, idcc );
}

std::vector<int> 
MEEBGeom::lmmodFromLmr( int ilmr )
{
  std::pair< int, int > dccAndSide_ = ME::dccAndSide( ilmr );
  int iside = dccAndSide_.second;
  std::vector< int > vec;
  for( int ilmmod=1; ilmmod<=9; ilmmod++ )
    {
      if( (ilmmod+iside)%2==1 )  vec.push_back(ilmmod);
    }
  return vec;
}

int
MEEBGeom::apdRefTower( int ilmmod )
{
  switch( ilmmod )
    {
      /* case   1: return 1;
    case   2: return 5; 
    case   3: return 7;
    case   4: return 21; 
    case   5: return 23;
    case   6: return 37;
    case   7: return 39;
    case   8: return 53;
    case   9: return 55;
      */
    case   1: return 2;
    case   2: return 6; 
    case   3: return 8;
    case   4: return 22; 
    case   5: return 24;
    case   6: return 38;
    case   7: return 40;
    case   8: return 54;
    case   9: return 56;
    default:
      abort();
    }
  return 0;
}


std::vector< int >
MEEBGeom::apdRefChannels( int ilmmod )
{
  
  std::vector< int > vec;
  switch( ilmmod )
    {
    case   1:  vec.push_back(0); vec.push_back(24); break;
    case   2:  vec.push_back(0); vec.push_back(24); break;
    case   3:  vec.push_back(0); vec.push_back(24); break;
    case   4:  vec.push_back(0); vec.push_back(24); break;
    case   5:  vec.push_back(0); vec.push_back(24); break;
    case   6:  vec.push_back(0); vec.push_back(24); break;
    case   7:  vec.push_back(0); vec.push_back(24); break;
    case   8:  vec.push_back(0); vec.push_back(24); break;
    case   9:  vec.push_back(0); vec.push_back(24); break;
    default:
      abort();
    }
  return vec;
}