myTrdK.cxx

Go to the documentation of this file.

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

using namespace std;

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

ClassImp(myTrdK)
#ifdef _WITHGBATCH_
//ClassImp(template<class SeedType> myTrdKFiller)
#endif

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

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

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

void myTrdK::Clear(Option_t* option){
#ifdef PDEBUG
  printf("In myTrdK::Clear\n");
#endif
  PRINTDEBUG;
  IsReadAlignmentOK = -10;
  IsReadCalibOK = -10;
  LHRValid = false;
  LHRValidCorr = false;
  fill_n(LHR, 3, (int)(-10));
  fill_n(LHRCorr, 3, (int)(-10));
  fill_n(LH, 3, (int)(-10));
  fill_n(LHCorr, 3, (int)(-10));
  NHitsUsedInLH = -10;
  NHitsUsedInLHCorr = -10;
  OffTrackNHits = -10;
  OffTrackAmps = -10;
  LHRrefitValid = LHRrefitCorrValid = false;
  NHitsUsedInLHrefit = NHitsUsedInLHrefitCorr = -1;
  fill_n(LHRrefit, 3, -2);
  fill_n(LHrefit, 3, -2);
  fill_n(LHRrefitCorr, 3, -2);
  fill_n(LHrefitCorr, 3, -2);
  PRINTDEBUG;
  TRDCharge = -10;
  TRDChargeError = -10;
  QTRDHitCollectionNucleiSize = -10;
  PRINTDEBUG;
  return;
}

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

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

#ifdef _WITHGBATCH_

template<class SeedType>
myTrdKFiller<SeedType>::myTrdKFiller(){
#ifdef PDEBUG
  printf("In template<class SeedType> myTrdKFiller::myTrdKFromTrTrackFiller\n");
#endif
  PRINTDEBUG;
  init();
  PRINTDEBUG;
}

template<class SeedType>
myTrdKFiller<SeedType>::~myTrdKFiller(){
#ifdef PDEBUG
  printf("In template<class SeedType> myTrdKFiller::~myTrdKFromTrTrackFiller\n");
#endif
  PRINTDEBUG;
  Clear();
}

template<class SeedType>
void myTrdKFiller<SeedType>::Clear(Option_t* option){
#ifdef PDEBUG
  printf("In template<class SeedType> myTrdKFiller::Clear\n");
#endif
  PRINTDEBUG;
  return;
}

template<class SeedType>
void myTrdKFiller<SeedType>::init(){
#ifdef PDEBUG
  printf("In template<class SeedType> myTrdKFiller::init\n");
#endif
  PRINTDEBUG;
  Clear();
  PRINTDEBUG;
  return;
}

template<class SeedType>
void myTrdKFiller<SeedType>::Fill(short int i_track, SeedType* track, float Energy, float kDef){
#ifdef PDEBUG
  printf("In template<class SeedType> myTrdKFiller::Fill\n");
#endif
  PRINTDEBUG;

  AMSEventR* pev = AMSEventR::Head();

  int NTrdRawHits = pev->nTrdRawHit();

#if !defined _B524_ && !defined _B550_ && !defined _B572_ 
  
  if (kDef>=0 && track && NTrdRawHits>0) {  
    
    TrdKCluster trdkclu(pev, track, kDef); 
    
    IsReadAlignmentOK= trdkclu.IsReadAlignmentOK;  // 0: Alignment not performed,  1: Static Alignment of Layer level,  2: Dynamic Alignment for entire TRD 
    IsReadCalibOK= trdkclu.IsReadCalibOK;  // 0: Gain Calibration not performed,  1: Gain Calibration Succeeded
    
    if (IsReadAlignmentOK && IsReadCalibOK) {
      LHRValid = trdkclu.GetLikelihoodRatio_TrTrack(trdkADCth, LHR, NHitsUsedInLH, 0, LH);
      LHRValidCorr = trdkclu.GetLikelihoodRatio_TrTrack(trdkADCth, LHRCorr, NHitsUsedInLHCorr, Energy, LHCorr);
      //Get Number of surrounding fired tubes, excluding the ones crossed by the current TrTrack prediction
      trdkclu.GetOffTrackHit_TrTrack(OffTrackNHits, OffTrackAmps);
      //      printf("Likelihood = {%f, %f %f}\n", LHR[0],  LHR[1],  LHR[2]);//only for debug
      LHRrefitValid = trdkclu.GetLikelihoodRatio_TRDRefit(trdkADCth, LHRrefit, NHitsUsedInLHrefit, 0, LHrefit);
      LHRrefitCorrValid = trdkclu.GetLikelihoodRatio_TRDRefit(trdkADCth, LHRrefitCorr, NHitsUsedInLHrefitCorr, Energy, LHrefitCorr);
      
#if !defined _B524_ && !defined _B530_ && !defined _B538_ && !defined _B550_ && !defined _B572_ && !defined _B580_ && !defined _B584_ && !defined _B594_ && !defined _B598_ && !defined _B610_ && !defined _B620_
      //Calculate TRD charge, this automatically call a new Trd Refit according to  "charge nuclei" hypothesis
      trdkclu.CalculateTRDCharge();
      //Get value of TRD charge, suggest discarding events with QTRD=0
      TRDCharge = trdkclu.GetTRDCharge();
      //Get one sigma error of TRD charge
      TRDChargeError = trdkclu.GetTRDChargeError();
      //Number of hits used in likelihood calculation, suggest discarding events with NHits<15
      QTRDHitCollectionNucleiSize = trdkclu.GetQTRDHitCollectionNuclei().size();
    }
#endif //after B630

    PRINTDEBUG;
        
  }

#endif //#if !defined _B524_ && !defined _B550_ && !defined _B572_ //both for TRDK and TRDQT

  PRINTDEBUG;

  return;
}

/* explicit instantiation */
template class myTrdKFiller<TrTrackR>;
template class myTrdKFiller<TrdTrackR>;

#endif