myBetaHPlus.cxx

Go to the documentation of this file.

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

using namespace std;

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

ClassImp(myBetaHPlus);
#ifdef _WITHGBATCH_
ClassImp(myBetaHPlusFiller);
#endif

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

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

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

void myBetaHPlus::Clear(Option_t* option){
#ifdef PDEBUG
  printf("In myBetaHPlus::Clear\n");
#endif
  PRINTDEBUG;
  Q = -999999;
  QRMS = -999999;
  QNPoints = -999999;
  PRINTDEBUG;
  Z = -999999;
  ZNPoints = -999999;
  ZProb = -999999;
  PRINTDEBUG;
  ZLHNPoints = -999999;
  fill_n(QL,       4, 0);
  fill_n(ZLH,     10, 0);
  fill_n(ZLHProb, 10, 0);
  PRINTDEBUG;
  NormChi2T=-999999;
  NormChi2C=-999999;
  PRINTDEBUG;

  return;
}

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

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

#ifdef _WITHGBATCH_

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

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

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

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

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

  betah = _betah;

  if (betah) {

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

    for (int ii=0; ii<4; ii++) {
      QL[ii] = betah->GetQL(ii);
    }

    Z=betah->GetZ(ZNPoints, ZProb);

    TofChargeHR *tofhc = betah->pTofChargeH();
    for(int iZ=1; iZ<10; iZ++){
      ZLH[iZ-1] = tofhc->GetLkh(iZ, ZLHNPoints);
      ZLHProb[iZ-1] = tofhc->GetProbZ(iZ, ZLHNPoints);
    }
  
    NormChi2T = betah->GetNormChi2T();
    NormChi2C = betah->GetNormChi2C();

  }
  
  PRINTDEBUG;

  return;
}

#endif