myTrdQtFromTrTrack.cxx

Go to the documentation of this file.

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

using namespace std;

#define NEP 2.71828182845904

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

ClassImp(myTrdQtFromTrTrack);
#ifdef _WITHGBATCH_
ClassImp(myTrdQtFromTrTrackFiller);
#endif

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

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

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

void myTrdQtFromTrTrack::Clear(Option_t* option){
#ifdef PDEBUG
  printf("In myTrdQtFromTrTrack::Clear\n");
#endif
  PRINTDEBUG;
  IfProcessed = false;
  NumberOfActiveLayers = -999999;
  NHitsUsedInLH = -999999;
  fill_n(LHR, 3, (int)(-999999));
  IsInsideTrdGeometricalAcceptance = false;
  PassesTrdPreselectionCuts = false;
  UsefulForTrdParticleId = false;
  PRINTDEBUG;
  return;
}

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

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

#ifdef _WITHGBATCH_

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

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

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

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

void myTrdQtFromTrTrackFiller::Fill(short int i_track, TrTrackR* track, int kDef, float Energy){
#ifdef PDEBUG
  printf("In myTrdQtFromTrTrackFiller::Fill\n");
#endif
  PRINTDEBUG;

  AMSEventR* pev = AMSEventR::Head();
  
  int NTrdRawHits = pev->nTrdRawHit();
  
#ifdef AMS_ACQT_INTERFACE
  
#if !defined _B524_ && !defined _B550_ && !defined _B572_
  
  if (kDef>0 && track && NTrdRawHits>0) {  
  
#if !defined _B584_ && !defined _B610_ && !defined _B620_
  
    static TrdQt* trdQt = new TrdQt(0);
    
    // Query TRD information by passing the trackIndex 0 to operate on the first tracker track, and the 1,3,3 tracker track fit.
    IfProcessed = trdQt->ProcessEvent(pev, i_track, 1, 0, 1);
    NumberOfActiveLayers = trdQt->GetNumberOfActiveLayers();
    NHitsUsedInLH = trdQt->GetNumberOfActiveStraws();
    // test if chosen track is inside TRD acceptance
    IsInsideTrdGeometricalAcceptance = trdQt->IsInsideTrdGeometricalAcceptance();
    // TRD preselection
    PassesTrdPreselectionCuts = trdQt->PassesTrdPreselectionCuts();
    // enough information in TRD for reliable particle ID?
    UsefulForTrdParticleId = trdQt->UsefulForTrdParticleId();
    if (IfProcessed) {
      // get the e/p log-likelihood ratio
      LHR[0]=trdQt->LogLikelihoodRatioElectronProton();
      LHR[1]=-log(1.0-pow(NEP,-trdQt->LogLikelihoodRatioHeliumElectron()));
      LHR[2]=-log(1.0-pow(NEP,-trdQt->LogLikelihoodRatioHeliumProton()));
    }

#endif //#if !defined _B584_ && !defined _B610_ && !defined _B620_

  }

#endif //#if !defined _B524_ && !defined _B550_ && !defined _B572_

#endif //#ifdef AMS_ACQT_INTERFACE

  PRINTDEBUG;

  return;
}

#endif