LMD/TFndHLmd.cxx

// @(#)fROOT/LMD:$Name:  $:$Id: TFndHLmd.cxx,v 1.17 2007/09/24 07:32:41 Diego_Faso Exp $
// Revision Author: Diego Faso <mailto:faso@to.infn.it>, 2005/06/24
// Oroginal Author: P. Cerello <mailto:cerello@to.infn.it>, 1999; D. Faso <mailto:faso@to.infn.it>, 2006

//                                          //
//                 TFndHLmd                 //
//                                          //
//  Description of the hit on Lmd           //
//                                          //

#include <Riostream.h>

#include "TMath.h"

#include "TFndRun.h"
#include "TFndHLmd.h"

ClassImp(TFndHLmd)

//__________________________________________________
TFndHLmd::TFndHLmd(const Int_t &chamber,const Int_t &wire):
  TFndHit(),
  fLayer(),fChamber(),fWire(),
  fAdc(),fTdc(),fBothAdc(),fBothTdc(),
  fCalibAdc(),
  fCalibTdc_T0(),fCalibTdc_T0C(),fCalibMeanTimer_T0(),fCalibMeanTimer_T0C()
{
  SetDchIdentifiers(chamber,wire);
  
  Init();
  
  
}

//__________________________________________________
TFndHLmd::TFndHLmd(const Int_t &channel):
  TFndHit(channel),
  fLayer(),fChamber(),fWire(),
  fAdc(),fTdc(),fBothAdc(),fBothTdc(),
  fCalibAdc(),
  fCalibTdc_T0(),fCalibTdc_T0C(),fCalibMeanTimer_T0(),fCalibMeanTimer_T0C()
{

  Int_t w;
  TFndDB::Lmd_ChNum_To_DetChan(fChannel,fLayer,fChamber,fWire,w);
  
  Init();
}

//__________________________________________________
TFndHLmd::~TFndHLmd() {
  // ---
}

//__________________________________________________
void TFndHLmd::Init(){
  // initialize all variables to defaults
  
  // do not set fLayer,fChamber,fWire here
  for(Int_t i=0;i<2;i++){
    fAdc[i] = -1;
    fTdc[i] = -1;
    //
    fCalibTdc_T0[i]  = -1;
    fCalibTdc_T0C[i] = -1;
    fCalibAdc[i] = -1;
  }
  fCalibMeanTimer_T0  = -1;
  fCalibMeanTimer_T0C = -1;
  fBothAdc = -1;
  fBothTdc = -1;
  
}

//_______________________________
void TFndHLmd::EvalHit(){
  // evaluate all hit attributes
  // (get info from DB (TFndDB) when needed)
  
  // --- "both-ADC/TDC" flags
  fBothAdc = E_HIT_NO_CH;
  if(fAdc[E_FIN_SIDE_P] >0  && fAdc[E_FIN_SIDE_E] >0)  fBothAdc = E_HIT_BOTH_CH;
  if(fAdc[E_FIN_SIDE_P] <=0 && fAdc[E_FIN_SIDE_E] >0)  fBothAdc = E_HIT_E_CH;
  if(fAdc[E_FIN_SIDE_P] >0  && fAdc[E_FIN_SIDE_E] <=0) fBothAdc = E_HIT_P_CH;
  
  fBothTdc = E_HIT_NO_CH;
  if(fTdc[E_FIN_SIDE_P] >0  && fTdc[E_FIN_SIDE_E] >0)  fBothTdc = E_HIT_BOTH_CH;
  if(fTdc[E_FIN_SIDE_P] <=0 && fTdc[E_FIN_SIDE_E] >0)  fBothTdc = E_HIT_E_CH;
  if(fTdc[E_FIN_SIDE_P] >0  && fTdc[E_FIN_SIDE_E] <=0) fBothTdc = E_HIT_P_CH;
  
  // --- evaluation of the zeta coordinate
  //     (x and y set to zero for every wire)
  // //   Double_t fPosLRS[3]; // Position in the Local Reference System (cm) [x,y,z] 
  // //   Double_t fPosMRS[3]; // Position in the Main Reference System (cm) [x,y,z] 

  fPosLRS[0] = 0;
  fPosLRS[1] = 0;
  fPosLRS[2] = 0;

  TFndFeeMap *fDB = fndrun->GetFeeMap();

  for(Int_t sid=0;sid<2;sid++){ // loop on the two sides
    if(fAdc[sid] == -1) continue;
    fCalibAdc[sid] = ( (Float_t) (fAdc[sid]) ) - ( (Float_t) fDB->GetLmdCal_Q0(fLayer,fChamber,fWire,sid) );
    fCalibAdc[sid] = fCalibAdc[sid] * ( (Float_t) fDB->GetLmdCal_Q1(fLayer,fChamber,fWire,sid) );
  }
  
  if( fBothAdc == E_HIT_BOTH_CH) {
    fPosLRS[2] = 
      LmdZetaGeoAdjust((E_Fnd_Layer)fLayer) *
      (Double_t)(
                 (
                  (fCalibAdc[E_FIN_SIDE_E] - fCalibAdc[E_FIN_SIDE_P]) /
                  (fCalibAdc[E_FIN_SIDE_E] + fCalibAdc[E_FIN_SIDE_P]) )
                 );
    // position in MRS:
    gFndGeom->RefSys_LmdWireToMRS(fChannel,fPosLRS,fPosMRS);
    
  }

  // TDC calibrations (both COSM and central triggers)
  for(Int_t sid=0;sid<2;sid++){ // loop on the two sides
    if(fTdc[sid] == -1) continue;
    fCalibTdc_T0[sid]  = ((Float_t) (fDB->GetLmdCal_T0(fLayer,fChamber,fWire,sid))  ) - ((Float_t)fTdc[sid]) ;
    fCalibTdc_T0C[sid] = ((Float_t) (fDB->GetLmdCal_T0C(fLayer,fChamber,fWire,sid)) ) - ((Float_t)fTdc[sid]) ;
  }
  
  if(fBothTdc != E_HIT_NO_CH){
    Float_t meantim = (Float_t)( (fTdc[0] + fTdc[1]) ) / 2;
    fCalibMeanTimer_T0  = ((Float_t) (fDB->GetLmdCal_T0(fLayer,fChamber,fWire,2))  ) - meantim ;
    fCalibMeanTimer_T0C = ((Float_t) (fDB->GetLmdCal_T0C(fLayer,fChamber,fWire,2)) ) - meantim ;
  }
  
  // ---
  fDB = 0; // owned by TFndRun
  //  PrintHit(3); // this line can be useful for debug
}

//_______________________________
void TFndHLmd::PrintHit(Int_t depth)
{
  // to be completed

  if( depth <=0 ) return;

  TString msg = "";
  TString linemsg = " --- --- --- --- ---\n";
  msg+=linemsg;

  //  if(depth>0) Printf("Lmd Channel: %d",fChannel);
  
  linemsg.Form("Lmd Channel: %d\n",fChannel);
  msg+=linemsg;
  
  if(depth>1){
    linemsg.Form(" ---> Adc E: %d\n",fAdc[E_FIN_SIDE_E]); msg+=linemsg;
    linemsg.Form(" ---> Adc P: %d\n",fAdc[E_FIN_SIDE_P]); msg+=linemsg;
    linemsg.Form(" ---> Tdc E: %d\n",fTdc[E_FIN_SIDE_E]); msg+=linemsg;
    linemsg.Form(" ---> Tdc P: %d\n",fTdc[E_FIN_SIDE_P]); msg+=linemsg;
    linemsg.Form(" ---> zeta : %.2f\n",fPosLRS[2]); msg+=linemsg;
    if(depth>2){
      linemsg.Form(" ------> Calib. Adc E     : %.2f\n",fCalibAdc[E_FIN_SIDE_E]); msg+=linemsg;
      linemsg.Form(" ------> Calib. Adc P     : %.2f\n",fCalibAdc[E_FIN_SIDE_P]); msg+=linemsg;
      linemsg.Form(" ------> Calib. Tdc E(T0) : %.2f\n",fCalibTdc_T0[E_FIN_SIDE_E]); msg+=linemsg;
      linemsg.Form(" ------> Calib. Tdc P(T0) : %.2f\n",fCalibTdc_T0[E_FIN_SIDE_P]); msg+=linemsg;
      linemsg.Form(" ------> Calib. Tdc E(T0C): %.2f\n",fCalibTdc_T0C[E_FIN_SIDE_E]); msg+=linemsg;
      linemsg.Form(" ------> Calib. Tdc P(T0C): %.2f\n",fCalibTdc_T0C[E_FIN_SIDE_P]); msg+=linemsg;
      linemsg.Form(" ------> Calib. MT(T0)    : %.2f\n",fCalibMeanTimer_T0); msg+=linemsg;
      linemsg.Form(" ------> Calib. MT(T0C)   : %.2f\n",fCalibMeanTimer_T0C); msg+=linemsg;
    }
  }
  // ---
  msg.Resize(msg.Sizeof()-2);
  Printf("%s",msg.Data());

  if(depth>2) PrintPosition(0);
  // ---
}

---