StRFEmcTrigMaker.cxx

// *-- Author : Renee Fatemi
// 
// $Id: StRFEmcTrigMaker.cxx,v 1.8 2011/04/11 19:35:52 fisyak Exp $

#include "StRFEmcTrigMaker.h"
#include "StChain.h"
#include "StMuDSTMaker/COMMON/StMuEvent.h"
#include "StMuDSTMaker/COMMON/StMuDst.h"
#include "StMuDSTMaker/COMMON/StMuDstMaker.h"
#include "StMuDSTMaker/COMMON/StMuEmcCollection.h"
#include "StMuDSTMaker/COMMON/StMuTrack.h"
#include "StMuDSTMaker/COMMON/StMuEmcPoint.h"
#include "StMuDSTMaker/COMMON/StMuTrack.h"
#include "StEmcUtil/geometry/StEmcGeom.h"
#include "TH1F.h"
#include "TH2F.h"
#include "StEvent.h"
#include "StEventMaker/StEventMaker.h"
#include "StEvent/StEmcDetector.h"
#include "StEvent/StEmcCollection.h"
#include "StEventTypes.h"
#include "StEvent/StBbcTriggerDetector.h"
#include "TClonesArray.h"


/*
http://www.star.bnl.gov/STAR/comp/pkg/dev/StRoot/StUtilities/doc/StMessMgr.html
gMessMgr->Message("This is an error message.","E");

The message type is specified with a single letter:
"E" = "Error"
"I" = "Info"
"W" = "Warning"
"D" = "Debug"
"Q" = "QAInfo"

gMessMgr->Message("","W") << big_num << " seems too big." << endm;
gMessMgr->Info() << big_num << " seems too big." << endm;
char* myText = "hello";
gMessMgr->Message(myText);

*/


ClassImp(StRFEmcTrigMaker)

StRFEmcTrigMaker::StRFEmcTrigMaker(const char *name):StMaker(name){
  activeBBC=false;
  emcGeom =0;
}

StRFEmcTrigMaker::~StRFEmcTrigMaker(){
}

//_____________________________________________________________________________
//_____________________________________________________________________________

Int_t StRFEmcTrigMaker::Init(){
  if (!DataMode) {
    cout << "You are loading StMuDst" << endl;
  } else {
    cout << "You are loading StEvent!" << endl;
  }

  emcGeom = StEmcGeom::getEmcGeom("bemc");
  initHisto();

  return StMaker::Init();
}

//____________________________________________________
//____________________________________________________
void StRFEmcTrigMaker::Sum(int *sum,int *sumadd){

  (*sum)=(*sumadd)+(*sum);
}
//______________________________________________________________________
void StRFEmcTrigMaker::Max(int *max,int *maxcomp){
if ((*max) < (*maxcomp))
   (*max)=(*maxcomp);
}

//______________________________________________________________________

Int_t StRFEmcTrigMaker::getBBCtrig(){
  return bbcTrig;
}
//______________________________________________________________________
Int_t StRFEmcTrigMaker::getBEMCtrigHT(int thres){

  int test=0;
  if (BHTmaxt>thres)  test=1;
  if (BHTmaxt<=thres)  test=0;
  return test;
}
//______________________________________________________________________
Int_t StRFEmcTrigMaker::getBEMCtrigJP(int thres){

  int test=0;
  if (BJPmaxt>thres)  test =  1;
  if (BJPmaxt<=thres)  test =  0;
  return test;
}
//______________________________________________________________________
Int_t StRFEmcTrigMaker::getBEMCtrigTOT(int thres){


  int test=0; 
  if (BJPsumt>thres)  test =  1;
  if (BJPsumt<=thres)  test =  0;
  return test;
}
//______________________________________________________________________
Int_t StRFEmcTrigMaker::getEEMCtrigHT(int thres){

  int test=0;
  if (EHTmaxt>thres) test=1;
  if (EHTmaxt<=thres) test=0;
  return test;

}
//______________________________________________________________________
Int_t StRFEmcTrigMaker::getEEMCtrigJP(int thres ){


  int test=0;  
  if (EJPmaxt>thres)  test =  1;
  if (EJPmaxt<=thres)  test =  0;
  return test;
}
//______________________________________________________________________
Int_t StRFEmcTrigMaker::getEEMCtrigTOT(int thres ){

  int test=0;
  if (EJPsumt>thres)  test =  1;
  if (EJPsumt<=thres)  test =  0;
  return test;
}

//______________________________________________________________________
Int_t StRFEmcTrigMaker::getEEMC_HT_ADC(){

  return EHTmaxt;

}

//______________________________________________________________________
Int_t StRFEmcTrigMaker::getEEMC_JP_ADC(){

  return EJPmaxt;

}

//______________________________________________________________________
Int_t StRFEmcTrigMaker::getEEMC_TOT_ADC(){
  return EJPsumt;  

}

//______________________________________________________________________
Int_t StRFEmcTrigMaker::getBEMC_HT_ADC(){

  return BHTmaxt;

}

//______________________________________________________________________
Int_t StRFEmcTrigMaker::getBEMC_JP_ADC(){

  return BJPmaxt;

}

//______________________________________________________________________
Int_t StRFEmcTrigMaker::getBEMC_TOT_ADC(){

  return BJPsumt;

}

//_____________________________________________________
//_____________________________________________________
void  StRFEmcTrigMaker:: Clear(const char *opt){
  // printf("\n\n Clear is called, JB\n\n");
  
  bbcTrig=0;

   // set array values to zero before each event
   memset(jpBsum,0,sizeof(jpBsum));
   memset(jpBmax,0,sizeof(jpBmax));
   memset(jpB_hit_num,0,sizeof(jpB_hit_num));
   memset(tpBsum,0,sizeof(tpBsum));
   memset(tpBmax,0,sizeof(tpBmax));
   
   
   // set array values to zero before each event
   memset(jpEsum,0,sizeof(jpEsum));
   memset(jpEmax,0,sizeof(jpEmax));
   memset(jpE_hit_num,0,sizeof(jpE_hit_num));
   memset(tpEsum,0,sizeof(tpEsum));
   memset(tpEmax,0,sizeof(tpEmax));
   
   BHTmaxt=0;   //Hold HT for whole Barrel
   BJPmaxt=0;   //Holds max JP sum for whole Barrel
   BJPsumt=0;  //Holds sum of all JP in Barrel
   EHTmaxt=0;   //Hold HT for whole EEMC
   EJPmaxt=0;   //Holds max JP sum for whole EEMC
   EJPsumt=0;  //Holds sum of all JP in EEMC

   // clear pointer to data collections
   bbcCol=0;
   muEmcCol=0;
   stEmcCol=0;
}

//_____________________________________________________
//_____________________________________________________
Int_t StRFEmcTrigMaker::Make(){
  
  // access all existing collections

  switch( DataMode) {
  case 0: { // use muDst
    StMuDstMaker *muDstMaker = (StMuDstMaker*)GetMaker("MuDst");
    if(!muDstMaker) {
      gMessMgr->Warning() <<GetName()<<" no muDstMaker !!! Game Over!"<< endm;
      return kStErr;
    }
    muEvent=muDstMaker->muDst()->event();
    bbcCol=&(muEvent->bbcTriggerDetector());
    muEmcCol=muDstMaker->muDst()->muEmcCollection();
    unpackEmcFromMu();
    } break;
  case 1: {// stEvent    
    stEvent = (StEvent *) GetInputDS("StEvent");
    if(!stEvent) {
      gMessMgr->Warning() <<GetName()<<" no StEvent !!! Game Over!"<< endm;
     return kStErr; 
    }
    StTriggerDetectorCollection *TrigDet=stEvent->triggerDetectorCollection();
    bbcCol=&(TrigDet->bbc());
    stEmcCol =(StEmcCollection*)stEvent->emcCollection(); 
    unpackEmcFromSt();
    } break;
  default: 
    gMessMgr->Error() <<GetName()<<" Logic error 1, kill chain"<< endm;
    return kStFatal;
  }
  
  // common 
  unpackBBC();
  
  if(activeBBC &&  !getBBCtrig()) return kStErr;
  // accept this eve
  
  fillHisto();
  return kStOK;
}


//_____________________________________________________
//_____________________________________________________
void  StRFEmcTrigMaker::unpackEmcFromMu(){
  if(!muEmcCol) {
    gMessMgr->Warning() <<GetName()<<" no muEmcCollection"<< endm;
    return;
  }  
   
  // According to STAR NOTE#229A all detectors should be numbered first from the +z side (West End) looking toward the interactions region 
  // If a detector needs additional numbering on the -z side then the numbers should be consecutive with the +z elements.
  // Following this, standing on the west side looking at the interaction region, module 58 is at 12 o'clock and in JP0 with JP1 and so on
  // proceeding in a clockwise manner. Standing on the east side looking at the interaction region module 118 is at 12 o'clock 
  // and in JP6 with JP7 and so on proceeding in a clockwise manner.
  
  //JP0 goes from module=53/2 to module=3/1 (TP=0+26-29,30+56-59,60+86-89,90+116-119,120+146-149)
  //JP1 goes from module=3/2 to module=13/1
  //JP2 goes from module=13/2 to module=23/1
  //JP3 goes from module=23/2 to module=33/1
  //JP4 goes from module=33/2 to module=43/1
  //JP5 goes from module=43/2 to module=53/1 (TP=21-25,51-55,81-85,111-115,141-145)
  //JP6 goes from module=113/2 to module=63/1
  //JP7 goes from module=63/2 to module=73/1
  //JP8 goes from module=73/2 to module=83/1
  //JP9 goes from module=83/2 to module=93/1
  //JP10 goes from module=93/2 to module=103/1
  //JP11 goes from module=103/2 to module=113/1
  
  //TP(0-29) for eta bin (1-4) all modules 1-60
  //TP(30-59) for eta bin(5-8) all modules  1-60
  //TP(60-89) for eta bin(9-12) all modules 1-60 
  //TP(90-119) for eta bin(13-16) all modules 1-60 
  //TP(120-149) for eta bin(17-20) all modules 1-60
  //TP(150-179) for eta bin (1-4) all modules 61-120
  //TP(180-209) for eta bin(5-8) all modules 61-120 
  //TP(210-239) for eta bin(9-12) all modules  61-120
  //TP(240-269) for eta bin(13-16) all modules  61-120
  //TP(270-299) for eta bin(17-20) all modules 61-120
  
  det=1;
  for (int n=1; n<=BemcTow; n++){
    emcGeom->getBin(n,Bmod,Beta,Bsub);
    BTowADC =  muEmcCol->getTowerADC(n,det);
    if (BTowADC>0) {
      //printf("n=%d, mod=%d, sub=%d, Beta=%d adc=%d\n",n,Bmod,Bsub,Beta,BTowADC);
      
      int jpBindex=(Bmod+Bsub+5)/10;
      if (((Bmod+Bsub+5)>=60)&&((Bmod+Bsub+5)<=66)) {
        jpBindex=0;
      }
      if ((Bmod == 60)&&(Bsub==2)) jpBindex=0;
      
      int tpBindex=((Bmod+Bsub-3)/2) + 30*((Beta-1)/4);
      if ((Bmod==1)&&(Bsub==1)) {
        tpBindex=(29 + 30*((Beta-1)/4));
      }  
      if (Bmod>60){//need to add 150 to tp# for east side 
        tpBindex=150 + ((Bmod+Bsub-3)/2) + 30*((Beta-1)/4);
        if ((Bmod==61)&&(Bsub==1)) {
          tpBindex=150 + (59 + 30*((Beta+-1)/4));
        }  
      }
      Sum(&jpBsum[jpBindex],&BTowADC);
      Max(&jpBmax[jpBindex],&BTowADC);
      Sum(&tpBsum[tpBindex],&BTowADC);
      Max(&tpBmax[tpBindex],&BTowADC);
      jpB_hit_num[jpBindex]++;
      //printf("jpBindex=%d, jpBsum=%d, jpBmax=%d\n",jpBindex,jpBsum[jpBindex],jpBmax[jpBindex]);
      //printf("tpBindex=%d, tpBsum=%d, tpBmax=%d\n",tpBindex,tpBsum[tpBindex],tpBmax[tpBindex]);
    }
  }
  for (int q=0; q < 6; q++){
    Sum(&BJPsumt,&jpBsum[q]);
    Max(&BHTmaxt,&jpBmax[q]);
    Max(&BJPmaxt,&jpBsum[q]);
  }
  
  
  //Due to changes in EEMC muDST for 2004 sec,sub and eta all start at 1 and go to 12,5,12
  //This code will not work on MuDst before 2004 --actually it will but it will give you GARBAGE!!!
  //Tower 01TA01 has sec=1,sub=1,eta=1
  //Tower 01TA02 has sec=1,sub=1,eta=2
  //Define EMC towers in terms of jet patches (start 0 instead of 1)
  //By definition jp0=towers 11D-1C with Eid #0-119 (m=636-35)
  //              jp1=tower 1D-3C with Eid #120-239 (m=36-155)
  //              jp2=tower 3D-5C with Eid #240-359 (m=156-275)
  //              jp3=tower 3D-5C with Eid #359-479 
  //              jp4=tower 3D-5C with Eid #479-599
  //              jp5=tower 3D-5C with Eid #599-719
  //tp=0 is defined as Eids 0-2 + 12-14
  //tp=1 is defined as Eids 24-26 + 36-38
  //tp=30 is defined as Eid 3-6 + 15-18
  //tp=60 is defined as Eid 7-11 + 19-23        
  
  NumETow=muEmcCol->getNEndcapTowerADC();
  for (int m=0; m<NumETow; m++){
    muEmcCol->getEndcapTowerADC(m,ETowADC,Esec,Esub,Eeta);
    int Eid = m+84;
    if (ETowADC){
      //printf("m=%d,ETowADC=%d,Esec=%d,Esub=%d,Eeta=%d,Eid=%d\n",m,ETowADC,Esec,Esub,Eeta,Eid);
      if (Eid > 719) Eid=Eid-720;
      int jpEindex=Eid/120;
      int tpEindex=Eid/24;
      if (Eeta>7) {
        tpEindex+=60;
      }
      if ((Eeta>3)&&(Eeta<8)) {
        tpEindex+=30;
      }
      
      Sum(&jpEsum[jpEindex],&ETowADC);
      Max(&jpEmax[jpEindex],&ETowADC);
      Sum(&tpEsum[tpEindex],&ETowADC);
      Max(&tpEmax[tpEindex],&ETowADC);
      jpE_hit_num[jpEindex]++;
      
      //printf("Etow=%d, Esec=%d, Esub=%d, Eeta=%d ETowADC=%d\n",m,Esec,Esub,Eeta,ETowADC);
      //printf("Eid=%d, jpEindex=%d, jpEsum=%d, jpEmax=%d\n",Eid,jpEindex,jpEsum[jpEindex],jpEmax[jpEindex]);
      //printf("Eid=%d, tpEindex=%d, tpEsum=%d, tpEmax=%d\n",Eid,tpEindex,tpEsum[tpEindex],tpEmax[tpEindex]);
    }
  }
  
  
  for (int q=0; q < 6; q++){
    
    Sum(&EJPsumt,&jpEsum[q]);
    Max(&EHTmaxt,&jpEmax[q]);
    Max(&EJPmaxt,&jpEsum[q]);
    //printf("q=%d; JPmax=%d ,JPsum=%d\n",q,jpEmax[q],jpEsum[q]);
  }
  
  printf("EJPsum=%d ,EHTmax=%d,EJPmax=%d\n",EJPsumt,EHTmaxt,EJPmaxt);
  printf("BJPsum=%d ,BHTmax=%d,BJPmax=%d\n",BJPsumt,BHTmaxt,BJPmaxt);
}
 

 
//_____________________________________________________
//_____________________________________________________
void  StRFEmcTrigMaker::unpackEmcFromSt(){
  if(!stEmcCol) {
    gMessMgr->Warning() <<GetName()<<" no stEmcCollection"<< endm;
    return;
  }  
  // According to STAR NOTE#229A all detectors should be numbered first from the +z side (West End) looking toward the interactions region 
  // If a detector needs additional numbering on the -z side then the numbers should be consecutive with the +z elements.
  // Following this, standing on the west side looking at the interaction region, module 58 is at 12 o'clock and in JP0 with JP1 and so on
  // proceeding in a clockwise manner. Standing on the east side looking at the interaction region module 118 is at 12 o'clock 
  // and in JP6 with JP7 and so on proceeding in a clockwise manner.
  
  //JP0 goes from module=53/2 to module=3/1 (TP=0+26-29,30+56-59,60+86-89,90+116-119,120+146-149)
  //JP1 goes from module=3/2 to module=13/1
  //JP2 goes from module=13/2 to module=23/1
  //JP3 goes from module=23/2 to module=33/1
  //JP4 goes from module=33/2 to module=43/1
  //JP5 goes from module=43/2 to module=53/1 (TP=21-25,51-55,81-85,111-115,141-145)
  //JP6 goes from module=113/2 to module=63/1
  //JP7 goes from module=63/2 to module=73/1
  //JP8 goes from module=73/2 to module=83/1
  //JP9 goes from module=83/2 to module=93/1
  //JP10 goes from module=93/2 to module=103/1
  //JP11 goes from module=103/2 to module=113/1
  
  //TP(0-29) for eta bin (1-4) all modules 1-60
  //TP(30-59) for eta bin(5-8) all modules  1-60
  //TP(60-89) for eta bin(9-12) all modules 1-60 
  //TP(90-119) for eta bin(13-16) all modules 1-60 
  //TP(120-149) for eta bin(17-20) all modules 1-60
  //TP(150-179) for eta bin (1-4) all modules 61-120
  //TP(180-209) for eta bin(5-8) all modules 61-120 
  //TP(210-239) for eta bin(9-12) all modules  61-120
  //TP(240-269) for eta bin(13-16) all modules  61-120
  //TP(270-299) for eta bin(17-20) all modules 61-120
  
  
  //StDetectorId BemcId=StDetectorId(kBarrelEmcTowerId);
  StDetectorId BemcId=StDetectorId(kBarrelEmcTowerId);
  StEmcDetector *EmcDet = stEmcCol->detector(BemcId); //BEMC tower detector number  
  if(EmcDet) {
    for (UInt_t mod=1;mod<=EmcDet->numberOfModules();mod++){
      StEmcModule* module=EmcDet->module(mod); 
      StSPtrVecEmcRawHit& hit=module->hits();
      for(UInt_t ih=0;ih<hit.size();ih++){
        StEmcRawHit *x=hit[ih];
        Bmod=x->module();
        Bsub=x->sub();
        Beta=x->eta();
        BTowADC=x->adc();
        //printf("ih=%d, mod=%d eta=%d sub=%d adc=%d\n",ih,x->module(),x->eta(),x->sub(),x->adc());       
        if (BTowADC>0) {
          //      printf("mod=%d, sub=%d, Beta=%d adc=%d\n",Bmod,Bsub,Beta,BTowADC);
          int jpBindex=(Bmod+Bsub+5)/10;
          if (((Bmod+Bsub+5)>=60)&&((Bmod+Bsub+5)<=66)) {
            jpBindex=0;
          }
          if ((Bmod == 60)&&(Bsub==2)) jpBindex=0;
          
          int tpBindex=((Bmod+Bsub-3)/2) + 30*((Beta-1)/4);
          if ((Bmod==1)&&(Bsub==1)) {
            tpBindex=(29 + 30*((Beta-1)/4));
          }  
          if (Bmod>60){//need to add 150 to tp# for east side 
            tpBindex=150 + ((Bmod+Bsub-3)/2) + 30*((Beta-1)/4);
            if ((Bmod==61)&&(Bsub==1)) {
              tpBindex=150 + (59 + 30*((Beta+-1)/4));
            }  
          }
          Sum(&jpBsum[jpBindex],&BTowADC);
          Max(&jpBmax[jpBindex],&BTowADC);
          Sum(&tpBsum[tpBindex],&BTowADC);
          Max(&tpBmax[tpBindex],&BTowADC);
          jpB_hit_num[jpBindex]++;
          //printf("jpBindex=%d, jpBsum=%d, jpBmax=%d\n",jpBindex,jpBsum[jpBindex],jpBmax[jpBindex]);
          //printf("tpBindex=%d, tpBsum=%d, tpBmax=%d\n",tpBindex,tpBsum[tpBindex],tpBmax[tpBindex]);
        }
      }
    }
  
    for (int q=0; q < 6; q++){
      Sum(&BJPsumt,&jpBsum[q]);
      Max(&BHTmaxt,&jpBmax[q]);
      Max(&BJPmaxt,&jpBsum[q]);
    }
  }
  
  //Tower 01TA01 has sec=1,sub=1,eta=1
  //Tower 01TA02 has sec=1,sub=1,eta=2
  //Define EMC towers in terms of jet patches (start 0 instead of 1)
  //By definition jp0=towers 11D-1C with Eid #0-119 (m=636-35)
  //              jp1=tower 1D-3C with Eid #120-239 (m=36-155)
  //              jp2=tower 3D-5C with Eid #240-359 (m=156-275)
  //              jp3=tower 3D-5C with Eid #359-479 
  //              jp4=tower 3D-5C with Eid #479-599
  //              jp5=tower 3D-5C with Eid #599-719
  //tp=0 is defined as Eids 0-2 + 12-14
  //tp=1 is defined as Eids 24-26 + 36-38
  //tp=30 is defined as Eid 3-6 + 15-18
  //tp=60 is defined as Eid 7-11 + 19-23        
  
  StDetectorId EemcId=StDetectorId(kEndcapEmcTowerId);
  EmcDet = stEmcCol->detector(EemcId); //EEMC tower detector number
  if(EmcDet) {
    for (UInt_t mod=1;mod<=EmcDet->numberOfModules();mod++){
      StEmcModule* module=EmcDet->module(mod); 
      StSPtrVecEmcRawHit& hit=module->hits();
      for(UInt_t ih=0;ih<hit.size();ih++){
        StEmcRawHit *x=hit[ih];
        Esec=x->module();
        Esub=x->sub();
        Eeta=x->eta();
        ETowADC=x->adc();
        //printf("ih=%d, mod=%d eta=%d sub=%d adc=%d\n",ih,x->module(),x->eta(),x->sub(),x->adc());     
        int m =60*(Esec-1) + 12*(Esub-1) + (Eeta-1);
        int Eid = m+84;
        if (ETowADC){
          //printf("m=%d,ETowADC=%d,Esec=%d,Esub=%d,Eeta=%d,Eid=%d\n",m,ETowADC,Esec,Esub,Eeta,Eid);
          if (Eid > 719) Eid=Eid-720;
          int jpEindex=Eid/120;
          int tpEindex=Eid/24;
          if (Eeta>7) {
            tpEindex+=60;
          }
          if ((Eeta>3)&&(Eeta<8)) {
            tpEindex+=30;
          }
          
          Sum(&jpEsum[jpEindex],&ETowADC);
          Max(&jpEmax[jpEindex],&ETowADC);
          Sum(&tpEsum[tpEindex],&ETowADC);
          Max(&tpEmax[tpEindex],&ETowADC);
          jpE_hit_num[jpEindex]++;
          
          //printf("Etow=%d, Esec=%d, Esub=%d, Eeta=%d ETowADC=%d\n",m,Esec,Esub,Eeta,ETowADC);
          //printf("Eid=%d, jpEindex=%d, jpEsum=%d, jpEmax=%d\n",Eid,jpEindex,jpEsum[jpEindex],jpEmax[jpEindex]);
          //printf("Eid=%d, tpEindex=%d, tpEsum=%d, tpEmax=%d\n",Eid,tpEindex,tpEsum[tpEindex],tpEmax[tpEindex]);
        }      
      }
    }
    for (int q=0; q < 6; q++){
      
      Sum(&EJPsumt,&jpEsum[q]);
      Max(&EHTmaxt,&jpEmax[q]);
      Max(&EJPmaxt,&jpEsum[q]);
      //printf("q=%d; JPmax=%d ,JPsum=%d\n",q,jpEmax[q],jpEsum[q]);
    }   
    printf("EJPsum=%d ,EHTmax=%d,EJPmax=%d\n",EJPsumt,EHTmaxt,EJPmaxt);
    printf("BJPsum=%d ,BHTmax=%d,BJPmax=%d\n",BJPsumt,BHTmaxt,BJPmaxt);
  } // end of EEMC-det
}

 
 
//_____________________________________________________
//_____________________________________________________
void  StRFEmcTrigMaker::unpackBBC(){
  if(!bbcCol) {
    gMessMgr->Warning() <<GetName()<<" no BBC Collection"<< endm;
    return;
  }  
   int Npmt=bbcCol->numberOfPMTs();
   //  bbc->dump();
   
   int Wbbc=0;
   int Ebbc=0;
   for (int pmt=0;pmt<Npmt;pmt++){    
     BBCadc[pmt]=bbcCol->adc(pmt);
     int bbcadc=bbcCol->adc(pmt);
     if (bbcadc>5) {
       if (pmt<16) {
         //cout << "BBC EAST = true" << endl;
         Ebbc=1;
       }
       if (23<pmt && pmt<40) {
         //cout << "BBC WEST = true" << endl;
         Wbbc=1;
       }
     }
  }
   
   if ((Ebbc==1)&&(Wbbc==1)){
     bbcTrig=1;
   }

   gMessMgr->Info() <<GetName()<<" BBCtrig="<<getBBCtrig()<<endm;
   
}

//____________________________________________________
//____________________________________________________
void StRFEmcTrigMaker::initHisto() {
  int mxAdc=300;
  int nAdc=mxAdc/4;
  ha[0]=new TH1F("bHT","BEMC HT ADC",nAdc,-.5, mxAdc-.5);
  ha[1]=new TH1F("bJP","BEMC JP ADC",nAdc,-.5, mxAdc-.5);
  ha[2]=new TH1F("bTot","BEMC Tot ADC",nAdc,-.5, mxAdc-.5);
  ha[3]=new TH1F("eHT","EEMC HT ADC",nAdc,-.5, mxAdc-.5);
  ha[4]=new TH1F("eJP","EEMC JP ADC",nAdc,-.5, mxAdc-.5);
  ha[5]=new TH1F("eTot","EEMC Tot ADC",nAdc,-.5, mxAdc-.5);
}

//____________________________________________________
//____________________________________________________
void StRFEmcTrigMaker::fillHisto() {

  ha[0]->Fill(getBEMC_HT_ADC());
  ha[1]->Fill(getBEMC_JP_ADC());
  ha[2]->Fill(getBEMC_TOT_ADC());
  ha[3]->Fill(getEEMC_HT_ADC());
  ha[4]->Fill(getEEMC_JP_ADC());
  ha[5]->Fill(getEEMC_TOT_ADC());
}


// $Log: StRFEmcTrigMaker.cxx,v $
// Revision 1.8  2011/04/11 19:35:52  fisyak
// Replace uint by UInt_t, use TMath
//
// Revision 1.7  2004/10/21 13:31:34  balewski
// to match new name of emcCollection in muDst
//
// Revision 1.6  2004/08/18 19:52:49  balewski
// works for BBC
//
// Revision 1.5  2004/08/18 14:56:29  balewski
// trying to get BBC working
//