StRoot: StPmdUtil/StPmdMapUtil.cxx Source File

00001 /********************************************************
00002  *
00003  * $Id: StPmdMapUtil.cxx,v 1.11 2010/05/02 13:13:46 rashmi Exp $
00004  *
00005  * Author: Subhasis Chattopadhyay
00006  *
00007  *********************************************************
00008  *
00009  * Description: This the class of PMD Utillty which deals with maps of channel to sm,row,col.
00010  * It can be used for reverse mapping as well.
00011  *
00012  *********************************************************
00013  * $Log: StPmdMapUtil.cxx,v $
00014  * Revision 1.11  2010/05/02 13:13:46  rashmi
00015  *  removed some couts
00016  *
00017  * Revision 1.10  2010/04/30 07:05:46  rashmi
00018  * Fixed an array for nSmChainBoards
00019  *
00020  * Revision 1.9  2007/08/31 10:55:47  rashmi
00021  * Changed initialization of channelOR,CR,chainR to -1
00022  *
00023  * Revision 1.8  2007/06/06 04:00:00  perev
00024  * CleanUp
00025  *
00026  * Revision 1.7  2007/05/26 00:40:07  perev
00027  * Initialization added
00028  *
00029  * Revision 1.6  2007/05/23 18:49:55  rashmi
00030  * bug fix for unphysical year  23/05/07
00031  *
00032  * Revision 1.5  2007/05/21 04:38:11  rashmi
00033  * functions for SMChain Combination information 21/05/07
00034  *
00035  * Revision 1.3  2005/12/07 20:35:45  perev
00036  * iostream.h ==> iostream. WarnOff
00037  *
00038  * Revision 1.2  2005/01/27 13:10:08  subhasis
00039  * New map for 2005 data
00040  *
00041  * Revision 1.1  2004/06/24 13:54:39  subhasis
00042  *  Maputils added
00043  *
00044  *
00045  **********************************************************/
00046 #include <iostream> 
00047 #include <strings.h>
00048 #include <stdlib.h>
00049 #include <TROOT.h>
00050 #include<TMatrix.h>
00051 #include<TArrayF.h>
00052 #include <math.h>
00053 #include "StPmdGeom.h"
00054 #include "StPmdMapUtil.h"
00055 #include "StChain.h"
00056 ClassImp(StPmdMapUtil)
00057 
00058   StPmdMapUtil::StPmdMapUtil()         
00059 {
00060   memset(beg,0,end-beg);
00061   mPmdGeom = new StPmdGeom();
00062   StoreMapInfo();
00063 }
00064 
00065 StPmdMapUtil::~StPmdMapUtil(){/*none*/}  
00066 
00067 void StPmdMapUtil::StoreMapInfo(Int_t runno1)
00068 {
00069   /*
00070   int year = runno1/1000000;
00071   int rn   = (runno1/1000)%1000;
00072   */
00073   Int_t year=0;
00074   Int_t rn = 0;
00075   mPmdGeom->GetRunYear(runno1,rn,year);
00076   //  cout<<"runno, rn1 "<<runno1<<" "<<rn<<" "<<year<<endl;
00077   cout<<"PMD runno, rn1 "<<runno1<<" "<<rn<<" "<<year<<endl;
00078   if(year<4){
00079     gMessMgr->Warning("Unphysical year");
00080     return;
00081   }
00082 
00084   Int_t sm=0,row=0,col=0;
00085   //Initialize
00086   //memset(m_ChannelInBoard    [0][0],0,sizeof(m_ChannelInBoard    ));
00087   //   memset(m_TempChannelInBoard[0][0],0,sizeof(m_TempChannelInBoard));
00088   //   memset(m_Chain             [0][0],0,sizeof(m_Chain             ));
00089   
00090   //Initializing channels to -1; channel=0 is a physical channel; RR 2007-08-30
00091   for(Int_t ism=0;ism<PMD_CRAMS_MAX*2;ism++){
00092     for(Int_t irow=0;irow<PMD_ROW_MAX;irow++){
00093       for(Int_t icol=0;icol<PMD_COL_MAX;icol++){
00094         m_ChannelInBoard[ism][irow][icol]=-1;
00095         m_TempChannelInBoard[ism][irow][icol]=-1;
00096         m_Chain[ism][irow][icol]=-1;
00097       }
00098     }
00099   }  
00100   
00101   for(Int_t Chain_No=1;Chain_No<=PMD_CHAIN_MAX;Chain_No++){
00102     //    Int_t AliveBoards = mPmdGeom->GetNBoardsChain(Chain_No);
00103     for(Int_t CHANNEL=0;CHANNEL<PMD_CHAIN_CHANNEL_MAX;CHANNEL++){
00104       Int_t channel=CHANNEL;
00105       Int_t chtemp;
00106       
00107       int fail = 0;
00108       if(year==5) {
00109         fail = mPmdGeom->ChainMapping(Chain_No,channel,sm,col,row,chtemp);
00110       }  else if(year>5){
00111         mPmdGeom->readBoardDetail(runno1);
00112         fail = mPmdGeom->ChainMapping(Chain_No,channel,sm,col,row,chtemp,year);
00113       }
00114       if (sm <=0) fail +=10;
00115       if (col<=0) fail +=20;
00116       if (row<=0) fail +=40;
00117       if (fail) continue;
00118       if(chtemp==-1)continue;
00119       //  if(sm>12)cout<<"chain,ch,sm,col,row,chtemp,year "<<Chain_No<<" "<<channel<<" "<<sm-1<<" "<<col-1<<" "<<row-1<<" "<<chtemp<<" "<<year<<endl;
00120       m_ChannelInBoard[sm-1][row-1][col-1]=channel;
00121       m_TempChannelInBoard[sm-1][row-1][col-1]=chtemp;
00122       m_Chain[sm-1][row-1][col-1]=Chain_No;
00123     }
00124   }
00125 }
00126 
00127 void StPmdMapUtil::StoreMapInfo()
00128 {
00129   Int_t sm,row,col;
00130   //Initialize
00131   memset(m_ChannelInBoard    [0][0],0,sizeof(m_ChannelInBoard    ));
00132   memset(m_TempChannelInBoard[0][0],0,sizeof(m_TempChannelInBoard));
00133   memset(m_Chain             [0][0],0,sizeof(m_Chain             ));
00134   
00135   for(Int_t Chain_No=1;Chain_No<=PMD_CHAIN_MAX;Chain_No++){
00136     for(Int_t CHANNEL=0;CHANNEL<PMD_CHAIN_CHANNEL_MAX;CHANNEL++){
00137       Int_t channel=CHANNEL;
00138       Int_t chtemp;
00139 
00140         mPmdGeom->ChainMapping(Chain_No,channel,sm,col,row,chtemp);
00141         m_ChannelInBoard[sm-1][row-1][col-1]=channel;
00142         m_TempChannelInBoard[sm-1][row-1][col-1]=chtemp;
00143         m_Chain[sm-1][row-1][col-1]=Chain_No;
00144 
00145     }
00146   }
00147 }
00148 
00149 void StPmdMapUtil::SetSMChainCombos(Int_t runno1){
00150   
00151   nSMChain = 0;
00152   mPmdGeom->readBoardDetail(runno1);
00153   Int_t rn,year;
00154   mPmdGeom->GetRunYear(runno1,rn,year);
00155   for(Int_t ichain = 0;ichain<48;ichain++){
00156     for(Int_t ism = 0;ism < 24;ism++){
00157       nBoardSMChain[ism][ichain]  = 0;
00158     }
00159   }
00160   
00161   Int_t mapcheck;
00162   for(Int_t ichain = 0;ichain<48;ichain++){
00163     for(Int_t ibrd = 0;ibrd < 27;ibrd++){    
00164       Int_t channel = 10 + ibrd*64;
00165       Int_t sm=0,col=0,row=0,chmod=0;
00166       Int_t chain = ichain+1;
00167       if (year==5){
00168         mapcheck = mPmdGeom->ChainMapping(chain,channel,sm,col,row,chmod);
00169       }else{
00170         mapcheck = mPmdGeom->ChainMapping(chain,channel,sm,col,row,chmod,year);
00171       }
00172       if(mapcheck==kStOk && sm > 0) {
00173         nBoardSMChain[sm-1][ichain]++;
00174         SMChainExists[sm-1][ichain]=nSMChain+1;
00175         if(nBoardSMChain[sm-1][ichain]==1) {
00176           nSMChain++;
00177           //      cout<<nSMChain-1<<"  "<<ichain+1<<"  "<<ibrd+1<<"  "<<sm<<endl;
00178           FirstBoard[nSMChain-1] = ibrd+1;
00179           SM_Combo[nSMChain-1] = sm;
00180           Chain_Combo[nSMChain-1] = chain;
00181         }
00182       }
00183     }
00184   }
00185   
00186   /*    
00187   for (Int_t ismchain = 0; ismchain<nSMChain;ismchain++){
00188     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;
00189   }
00190   */
00191 }