EmcRecTofDigitCalib.cxx

Go to the documentation of this file.

//
//  Tof digit to hit
//
//  He Miao 2006, 9,26 
//
#include "EmcRec/EmcRecTofDigitCalib.h"
#include "Identifier/TofID.h"
#include "RawEvent/RawDataUtil.h"
 
EmcRecTofDigitCalib::EmcRecTofDigitCalib()
{}
 
EmcRecTofDigitCalib::~EmcRecTofDigitCalib()
{}
 
void EmcRecTofDigitCalib::Convert(TofDigiCol* tofDigiCol,
    RecEmcTofHitMap& tofHitMap)
{
  TofDigiCol::iterator iter2;
  for(iter2=tofDigiCol->begin();
      iter2!=tofDigiCol->end();
      iter2++,iter2++) {
    
    Identifier id=(*iter2)->identify();
    bool   is_barrel    = TofID::is_barrel(id);
    int    layer        = TofID::layer(id);
    int    im           = TofID::phi_module(id);
    
    if(!is_barrel) return;
    
    double adc0 = RawDataUtil::TofCharge((*iter2)->getChargeChannel())*2;
    double tdc0 = RawDataUtil::TofTime((*iter2)->getTimeChannel());
    double adc1 = RawDataUtil::TofCharge((*(iter2+1))->getChargeChannel())*2;
    double tdc1 = RawDataUtil::TofTime((*(iter2+1))->getTimeChannel());
    
    double rand0 = (double)(rand()%100)/1000.-0.05;   //-0.05 --> 0.05
    double rand1 = (double)(rand()%100)/1000.-0.05;   //-0.05 --> 0.05
 
    //smear 5%
    adc0 *= (1+rand0);
    adc1 *= (1+rand1);
 
    double zpos, energy, e0, e1;
    zpos = 0.08331*(tdc1-tdc0);
    if(fabs(zpos)>1.15) continue;
 
    //ADC->energy
    if((adc0>46.&&adc0<920.) && (adc1<=46.||adc1>=920.)) {
      energy = adc0/exp(2.783+0.3133*zpos);
    } else if((adc1>46.&&adc1<920.) && (adc0<=46.||adc0>=920.)) {
      energy = adc1/exp(2.784-0.3079*zpos);
    } else if((adc0>46.&&adc0<920.) && (adc1>46.&&adc1<920.)) {
      e0 = adc0/exp(2.783+0.3133*zpos);
      e1 = adc1/exp(2.784-0.3079*zpos);
      energy = (e0+e1)/2;
    } else if(adc0>=920.&&adc1>=920.) {
      e0 = 920./exp(2.783+0.3133*zpos);
      e1 = 920./exp(2.784-0.3079*zpos);
      energy = e0>e1?e0:e1;
    } else {
      continue;
    }
 
    double phi,r;
    if(layer==0) {
      phi=(im+0.5)*CLHEP::twopi/88.;
      r=838.5;
    } else {
      phi=im*CLHEP::twopi/88.;
      r=895.5;
    }
 
    zpos*=1000; //zpos: m->mm
    HepPoint3D pos(r*cos(phi), r*sin(phi), zpos);
     
    tofHitMap[im+layer*88] = RecEmcTofHit(id, energy, pos);
  }
}