myTrTrack.cxx

Go to the documentation of this file.

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

using namespace std;

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

ClassImp(myTrTrack);
#ifdef _WITHGBATCH_
ClassImp(myTrTrackFiller);
#endif //#ifdef _WITHGBATCH_

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

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

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

void myTrTrack::Clear(Option_t* option){
#ifdef PDEBUG
  printf("In myTrTrack::Clear\n");
#endif
  PRINTDEBUG;
  kDef = -999999;
  BitPatternJ = 0;
  BitPatternXYJ = 0;
  Nhits = -9999;
  NhitsX = -9999;
  NhitsY = -9999;
  NhitsXY = -9999;
  Rigidity=9999999;
  ErrRinv=-999999;
  Chisq=9999999;
  NormChisqY=9999999;
  NormChisqX=9999999;
  Theta=-999999;
  Phi=-999999;
  PRINTDEBUG;
  Q = -999999;
  InnerQ = -999999;
  fill_n(LayerJQ, 9, 0);
  PRINTDEBUG;
  fill_n(TrackXOnLayer, 9, 0);
  fill_n(TrackYOnLayer, 9, 0);
  PRINTDEBUG;
  return;
}

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

Double_t myTrTrack::TkBeta(Double_t MASS) {

  Double_t tkbeta = 0;

  //  printf("Mass is %f\n", MASS);//only for debug
  tkbeta=MASS*MASS;
  tkbeta+=Rigidity*Rigidity;
  tkbeta=Rigidity/tkbeta;
  tkbeta*=Rigidity;//b^2 = p^2/(p^2+m^2)
  tkbeta=sqrt(tkbeta);

  return tkbeta;
}

Short_t myTrTrack::span() {

  Short_t span=-999;
  
  PRINTDEBUG;
  
  bool HasL1 = ((BitPatternJ)&(1<<0))>>0;
  bool HasL2 = ((BitPatternJ)&(1<<1))>>1;
  bool HasL9 = ((BitPatternJ)&(1<<8))>>8;
  // printf("%d %d ", HasL1, HasL2);//only for debug
  // for (int ii=2; ii<8; ii++) printf("%d ", ((BitPatternJ)&(1<<ii))>>ii);//only for debug
  // printf("%d\n", HasL9);//only for debug
  
  if (HasL1&HasL9) {
    //      printf("HasL1&HasL9\n");//only for debug
    if (kDef>0) {
      span=0;
    }
  }
  else if (HasL1) {
    //      printf("HasL1\n");//only for debug
    if (kDef>0) {
      span=1;
    }
  }
  else if (HasL9) {
    //      printf("HasL9\n");//only for debug
    if (kDef>0) {
      if (HasL2) {
        //        printf("HasL2\n");//only for debug
        span=2;
      }
      else span=3;
    }
  }
  else {
    //      printf("Inner\n");//only for debug
    if (kDef>0) {
      if (HasL2) {
        //        printf("HasL2\n");//only for debug
        span=4;
      }
      else span=5;
    }
  }
  
  PRINTDEBUG;

  return span;
}

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

#ifdef _WITHGBATCH_

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

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

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

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

void myTrTrackFiller::Fill(TrTrackR* _track){
#ifdef PDEBUG
  printf("In myTrTrackFiller::Fill\n");
#endif
  PRINTDEBUG;

  track=_track;

  //--------------------TRACKER-----------------------------------  
  if (track) {

    BitPatternJ=track->GetBitPatternJ();
    BitPatternXYJ=track->GetBitPatternXYJ();
    
    NhitsX = track->GetNhitsX();
    NhitsY = track->GetNhitsY();
    NhitsXY = track->GetNhitsXY();
    Nhits = track->GetNhits();
    
    PRINTDEBUG;
    kDef=0;

    Rigidity = track->GetRigidity(kDef);
    ErrRinv = track->GetErrRinv(kDef);
    Chisq = track->GetChisq(kDef);
    NormChisqX = track->GetNormChisqX(kDef);
    NormChisqY = track->GetNormChisqY(kDef);
    Theta = track->GetTheta(kDef);
    // theta = PI - track->GetTheta(kDef);
    // if (theta>PI) {
    //  printf("Folding Theta (%f)\n", theta); //only for debug
    //  theta=PI + track->GetTheta(kDef);
    // }
    // else if (theta<(PI/2.0)) {
    //  printf("Theta not valid (%f)\n", theta); //only for debug
    // }
    Phi = track->GetPhi(kDef);
            
    for (int jlayer=1; jlayer<=9; jlayer++) {
      static AMSPoint global;
      static AMSDir dir;
#ifdef _B524_
      int jlayer=layer-1;
      if (layer==1) jlayer=8;
      else if (layer==9) jlayer=9;
      track->InterpolateLayer(layer-1, global, dir, kDef);
#else
      track->InterpolateLayerJ(jlayer, global, dir, kDef);
#endif
      TrackXOnLayer[jlayer-1]=global.x();
      TrackYOnLayer[jlayer-1]=global.y();
    }
    
    

    Q = track->GetQ();
    InnerQ = track->GetInnerQ();
    for (int ii=0; ii<9; ii++) {
      LayerJQ[ii] = track->GetLayerJQ(ii+1);
    }
    
  }

  PRINTDEBUG;
 
  return;
}

#endif //#ifdef _WITHGBATCH_