StPmdMapUtil.cxx

/********************************************************
 *
 * $Id: StPmdMapUtil.cxx,v 1.11 2010/05/02 13:13:46 rashmi Exp $
 *
 * Author: Subhasis Chattopadhyay
 *
 *********************************************************
 *
 * Description: This the class of PMD Utillty which deals with maps of channel to sm,row,col.
 * It can be used for reverse mapping as well.
 *
 *********************************************************
 * $Log: StPmdMapUtil.cxx,v $
 * Revision 1.11  2010/05/02 13:13:46  rashmi
 * removed some couts
 *
 * Revision 1.10  2010/04/30 07:05:46  rashmi
 * Fixed an array for nSmChainBoards
 *
 * Revision 1.9  2007/08/31 10:55:47  rashmi
 * Changed initialization of channelOR,CR,chainR to -1
 *
 * Revision 1.8  2007/06/06 04:00:00  perev
 * CleanUp
 *
 * Revision 1.7  2007/05/26 00:40:07  perev
 * Initialization added
 *
 * Revision 1.6  2007/05/23 18:49:55  rashmi
 * bug fix for unphysical year  23/05/07
 *
 * Revision 1.5  2007/05/21 04:38:11  rashmi
 * functions for SMChain Combination information 21/05/07
 *
 * Revision 1.3  2005/12/07 20:35:45  perev
 * iostream.h ==> iostream. WarnOff
 *
 * Revision 1.2  2005/01/27 13:10:08  subhasis
 * New map for 2005 data
 *
 * Revision 1.1  2004/06/24 13:54:39  subhasis
 *  Maputils added
 *
 *
 **********************************************************/
#include <iostream> 
#include <strings.h>
#include <stdlib.h>
#include <TROOT.h>
#include<TMatrix.h>
#include<TArrayF.h>
#include <math.h>
#include "StPmdGeom.h"
#include "StPmdMapUtil.h"
#include "StChain.h"
ClassImp(StPmdMapUtil)

  StPmdMapUtil::StPmdMapUtil()         
{
  memset(beg,0,end-beg);
  mPmdGeom = new StPmdGeom();
  StoreMapInfo();
}

StPmdMapUtil::~StPmdMapUtil(){/*none*/}  

void StPmdMapUtil::StoreMapInfo(Int_t runno1)
{
  /*
  int year = runno1/1000000;
  int rn   = (runno1/1000)%1000;
  */
  Int_t year=0;
  Int_t rn = 0;
  mPmdGeom->GetRunYear(runno1,rn,year);
  //  cout<<"runno, rn1 "<<runno1<<" "<<rn<<" "<<year<<endl;
  cout<<"PMD runno, rn1 "<<runno1<<" "<<rn<<" "<<year<<endl;
  if(year<4){
    gMessMgr->Warning("Unphysical year");
    return;
  }

  Int_t sm=0,row=0,col=0;
  //Initialize
  //memset(m_ChannelInBoard    [0][0],0,sizeof(m_ChannelInBoard    ));
  //   memset(m_TempChannelInBoard[0][0],0,sizeof(m_TempChannelInBoard));
  //   memset(m_Chain             [0][0],0,sizeof(m_Chain             ));
  
  //Initializing channels to -1; channel=0 is a physical channel; RR 2007-08-30
  for(Int_t ism=0;ism<PMD_CRAMS_MAX*2;ism++){
    for(Int_t irow=0;irow<PMD_ROW_MAX;irow++){
      for(Int_t icol=0;icol<PMD_COL_MAX;icol++){
        m_ChannelInBoard[ism][irow][icol]=-1;
        m_TempChannelInBoard[ism][irow][icol]=-1;
        m_Chain[ism][irow][icol]=-1;
      }
    }
  }  
  
  for(Int_t Chain_No=1;Chain_No<=PMD_CHAIN_MAX;Chain_No++){
    //    Int_t AliveBoards = mPmdGeom->GetNBoardsChain(Chain_No);
    for(Int_t CHANNEL=0;CHANNEL<PMD_CHAIN_CHANNEL_MAX;CHANNEL++){
      Int_t channel=CHANNEL;
      Int_t chtemp;
      
      int fail = 0;
      if(year==5) {
        fail = mPmdGeom->ChainMapping(Chain_No,channel,sm,col,row,chtemp);
      }  else if(year>5){
        mPmdGeom->readBoardDetail(runno1);
        fail = mPmdGeom->ChainMapping(Chain_No,channel,sm,col,row,chtemp,year);
      }
      if (sm <=0) fail +=10;
      if (col<=0) fail +=20;
      if (row<=0) fail +=40;
      if (fail) continue;
      if(chtemp==-1)continue;
      //  if(sm>12)cout<<"chain,ch,sm,col,row,chtemp,year "<<Chain_No<<" "<<channel<<" "<<sm-1<<" "<<col-1<<" "<<row-1<<" "<<chtemp<<" "<<year<<endl;
      m_ChannelInBoard[sm-1][row-1][col-1]=channel;
      m_TempChannelInBoard[sm-1][row-1][col-1]=chtemp;
      m_Chain[sm-1][row-1][col-1]=Chain_No;
    }
  }
}

void StPmdMapUtil::StoreMapInfo()
{
  Int_t sm,row,col;
  //Initialize
  memset(m_ChannelInBoard    [0][0],0,sizeof(m_ChannelInBoard    ));
  memset(m_TempChannelInBoard[0][0],0,sizeof(m_TempChannelInBoard));
  memset(m_Chain             [0][0],0,sizeof(m_Chain             ));
  
  for(Int_t Chain_No=1;Chain_No<=PMD_CHAIN_MAX;Chain_No++){
    for(Int_t CHANNEL=0;CHANNEL<PMD_CHAIN_CHANNEL_MAX;CHANNEL++){
      Int_t channel=CHANNEL;
      Int_t chtemp;

        mPmdGeom->ChainMapping(Chain_No,channel,sm,col,row,chtemp);
        m_ChannelInBoard[sm-1][row-1][col-1]=channel;
        m_TempChannelInBoard[sm-1][row-1][col-1]=chtemp;
        m_Chain[sm-1][row-1][col-1]=Chain_No;

    }
  }
}

void StPmdMapUtil::SetSMChainCombos(Int_t runno1){
  
  nSMChain = 0;
  mPmdGeom->readBoardDetail(runno1);
  Int_t rn,year;
  mPmdGeom->GetRunYear(runno1,rn,year);
  for(Int_t ichain = 0;ichain<48;ichain++){
    for(Int_t ism = 0;ism < 24;ism++){
      nBoardSMChain[ism][ichain]  = 0;
    }
  }
  
  Int_t mapcheck;
  for(Int_t ichain = 0;ichain<48;ichain++){
    for(Int_t ibrd = 0;ibrd < 27;ibrd++){    
      Int_t channel = 10 + ibrd*64;
      Int_t sm=0,col=0,row=0,chmod=0;
      Int_t chain = ichain+1;
      if (year==5){
        mapcheck = mPmdGeom->ChainMapping(chain,channel,sm,col,row,chmod);
      }else{
        mapcheck = mPmdGeom->ChainMapping(chain,channel,sm,col,row,chmod,year);
      }
      if(mapcheck==kStOk && sm > 0) {
        nBoardSMChain[sm-1][ichain]++;
        SMChainExists[sm-1][ichain]=nSMChain+1;
        if(nBoardSMChain[sm-1][ichain]==1) {
          nSMChain++;
          //      cout<<nSMChain-1<<"  "<<ichain+1<<"  "<<ibrd+1<<"  "<<sm<<endl;
          FirstBoard[nSMChain-1] = ibrd+1;
          SM_Combo[nSMChain-1] = sm;
          Chain_Combo[nSMChain-1] = chain;
        }
      }
    }
  }
  
  /*    
  for (Int_t ismchain = 0; ismchain<nSMChain;ismchain++){
    cout<<"ismchain="<<ismchain<<" chain="<<Chain_Combo[ismchain]<<" sm="<<SM_Combo[ismchain]<<" has number of boards = "<<nBoardSMChain[SM_Combo[ismchain]-1][Chain_Combo[ismchain]-1]<<" first board="<<FirstBoard[ismchain]<<endl;
  }
  */
}