myBetaH.cxx

Go to the documentation of this file.

// Authors: M.Duranti - INFN di Perugia
#include "myBetaH.h"
#include "debug.h"
#include "TClass.h"

using namespace std;

//--------------------------------------------------------------------

ClassImp(myBetaH);
#ifdef _WITHGBATCH_
ClassImp(myBetaHFiller);
#endif //#ifdef _WITHGBATCH_

//--------------------------------------------------------------------

myBetaH::myBetaH(){
#ifdef PDEBUG
  printf("In myBetaH::myBetaH\n");
#endif
  PRINTDEBUG;
  init();
  PRINTDEBUG;
}

myBetaH::~myBetaH(){
#ifdef PDEBUG
  printf("In myBetaH::~myBetaH\n");
#endif
  PRINTDEBUG;
}

void myBetaH::Clear(Option_t* option){
#ifdef PDEBUG
  printf("In myBetaH::Clear\n");
#endif
  PRINTDEBUG;
  Q = -999999;
  QRMS = -999999;
  Beta = -999999;
  Chi2T = 9999;
  Chi2C = 9999;
  EBetaV = -999999;
  BetaPattern = 999999;
  Status=-999999;
  PRINTDEBUG;
  SumHit = -999;
  UseHit = -999;
  PRINTDEBUG;
  _iTrTrack = -9999;
  _iEcalShower = -9999;
  _iTrdTrack = -9999;
  PRINTDEBUG;
  fill_n(TofTimeOnLayer,  9, 0);
  fill_n(TofXOnLayer,     9, 0);
  fill_n(TofYOnLayer,     9, 0);
  PRINTDEBUG;
  return;
}

void myBetaH::init(){
#ifdef PDEBUG
  printf("In myBetaH::init\n");
#endif
  PRINTDEBUG;
  Clear();
  PRINTDEBUG;
  return;
}

//----------------------------------------------------------------

#ifdef _WITHGBATCH_

myBetaHFiller::myBetaHFiller(){
#ifdef PDEBUG
  printf("In myBetaHFiller::myBetaHFiller\n");
#endif
  PRINTDEBUG;
  init();
  PRINTDEBUG;
}

myBetaHFiller::~myBetaHFiller(){
#ifdef PDEBUG
  printf("In myBetaHFiller::~myBetaHFiller\n");
#endif
  PRINTDEBUG;
}

void myBetaHFiller::Clear(Option_t* option){
#ifdef PDEBUG
  printf("In myBetaHFiller::Clear\n");
#endif
  PRINTDEBUG;
  betah=0;
  PRINTDEBUG;
  return;
}

void myBetaHFiller::init(){
#ifdef PDEBUG
  printf("In myBetaHFiller::init\n");
#endif
  PRINTDEBUG;
  Clear();
  PRINTDEBUG;
  return;
}

void myBetaHFiller::Fill(BetaHR* _betah){
#ifdef PDEBUG
  printf("In myBetaHFiller::Fill\n");
#endif
  PRINTDEBUG;

  betah = _betah;

  //--------------------------TOF-------------------------------
  if (betah) {

    int QNPoints = -999999;
    Q = betah->GetQ(QNPoints, QRMS);

    Beta = betah->GetBeta();
    EBetaV = betah->GetEBetaV();
    BetaPattern = betah->GetBetaPattern();
    Status  =  betah->GetStatus();
    Chi2T = betah->GetChi2T();
    Chi2C = betah->GetChi2C();

    SumHit = betah->GetSumHit();
    UseHit = betah->GetUseHit();

    
    _iTrTrack = betah->iTrTrack();
    PRINTDEBUG;
    
    _iEcalShower = betah->iEcalShower();
    PRINTDEBUG;
    
    _iTrdTrack = betah->iTrdTrack();
    PRINTDEBUG;
    double toftime=0.;
    for (int jlayer=1; jlayer<=9; jlayer++) {
      static AMSPoint global;
      static AMSDir dir;

      double zA=TkDBc::Head->GetZlayerAJ (jlayer);
      double z= TkDBc::Head->GetZlayerJ (jlayer);
      
      betah->TInterpolate (zA, global, dir, toftime, false);
      TofTimeOnLayer[jlayer-1]= (Double32_t )toftime;
      TofXOnLayer[jlayer-1]=global.x();
      TofYOnLayer[jlayer-1]=global.y();
    }
  
    

    
  }
  
  PRINTDEBUG;

  return;
}

Double_t myBetaHFiller::GetTofCharge(int pattern){
  //pattern =0 ALL Tof Layer
  //pattern =1 Upper Tof
  //pattern =2 Lower Tof
#ifdef PDEBUG
  printf("In myBetaHFiller::GetTofCharge\n");
#endif
  PRINTDEBUG;
  if(!betah)  return 0;
  int first_layer=0;
  int last_layer=4;
  if (pattern==1){
    first_layer=0;
    last_layer=2;
  }else if (pattern==2){
    first_layer=2;
    last_layer=4;
  }
  
  Double_t tofh_q_upp = 0;
  Double_t tofh_n_upp = 0; 
  Double_t tofh_r_upp = 0;
  AMSPoint pnt;AMSDir dir; double time;
  for (int ilay=first_layer; ilay<last_layer; ilay++) {
    if(!betah->GetClusterHL(ilay)) continue;
    betah->TInterpolate(betah->GetClusterHL(ilay)->Coo[2],pnt,dir,time);
    if (!TOFGeom::IsGoodQGeom(ilay,betah->GetClusterHL(ilay)->Bar,pnt)) continue;
    Double_t q = betah->GetQL(ilay,2);
    if (q<=1e-6) continue;
    tofh_q_upp += q;
    tofh_r_upp += q*q;
    tofh_n_upp++;
  }
  
  if (tofh_n_upp>0) { 
    tofh_q_upp /= tofh_n_upp;
    tofh_r_upp /= tofh_n_upp;
    tofh_r_upp = sqrt(tofh_r_upp - tofh_q_upp*tofh_q_upp);
  }   
 
  return tofh_q_upp;
}

#endif //#ifdef _WITHGBATCH_