myRTI.cxx

Go to the documentation of this file.

// Authors: M.Duranti - INFN di Perugia
#include "myRTI.h"
#include "debug.h"
#include "TClass.h"
#include <iostream>
using namespace std;

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

ClassImp(myRTI);
#ifdef _WITHGBATCH_
ClassImp(myRTIFiller);
#endif //#ifdef _WITHGBATCH_
//--------------------------------------------------------------------

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

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

  run = 999999;
  evno = 999999;
  evnol= 999999;
  lf = -999999;
  memset(cf, 0, 4*2* sizeof(cf[0][0]));
  memset(cfi, 0, 4*2* sizeof(cfi[0][0]));
  memset(nl1l9, 0, 2*2* sizeof(nl1l9[0][0]));
  memset(dl1l9, 0, 2*3* sizeof(dl1l9[0][0]));
  /*
  fill_n(cf, 8, 0);
  fill_n(cfi, 8, 0);
  fill_n( nl1l9, 4, 0);
  fill_n( dl1l9, 6, 0);
  */
  mphe = -999999;
  theta = -999999;
  phi = -999999;
  r = -999999;
  zenith = -999999;
  glat = -999999;
  glong = -999999;
  nev = -999999;
  nerr = -999999;
  ntrig = -999999;
  nhwerr = -999999;
  npart = -999999;
  mtrdh = -999999;
  good = -9999999;
  utime = 999999;
  fill_n(usec, 2, 0);

  PRINTDEBUG;

  return;
}

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

bool myRTI::TPreSelect(){
#ifdef PDEBUG
  printf("In myRTI::TPreSelect\n");
#endif
  PRINTDEBUG;

  bool cut[10]={0};
  cut[0]=(ntrig/nev>0.98);
  cut[1]=(npart/ntrig>0.07/1600*ntrig&&npart/ntrig<0.25);
  cut[2]=(lf>0.5);
  cut[3]=(zenith<40);
  cut[4]=(nerr>=0&&nerr/nev<0.1);
  cut[5]=(npart>0&&nev<1800);
  bool tcut=(cut[0]&&cut[1]&&cut[2]&&cut[3]&&cut[4]&&cut[5]);
  cout<<" 0 "<< cut[0]<<" && "<<cut[1]<<" && "<<cut[2]<<" && "<<cut[3]<<" && "<<cut[4]<<endl;
  return tcut;
}

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

#ifdef _WITHGBATCH_

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

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

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

void myRTIFiller::Fill(AMSSetupR::RTI* ev_RTI){
#ifdef PDEBUG
  printf("In myRTIFiller::Fill\n");
#endif
  PRINTDEBUG;

  if (ev_RTI) {
    
#if !defined _B524_ && !defined _B530_ && !defined _B538_ && !defined _B550_ && !defined _B572_ && !defined _B580_ && !defined _B584_ && !defined _B594_ && !defined _B598_ && !defined _B610_ && !defined _B620_
    run = ev_RTI->run;
    evno = ev_RTI->evno;
    evnol= ev_RTI->evnol;
    //  printf("myRTIFiller::Fill: %d\n", evno);//only for debug
    lf = ev_RTI->lf;
    //  printf("myRTIFiller::Fill: %f\n", lf);//only for debug
    copy_n( ev_RTI->cf[0] ,  4*2,  cf[0]);
    copy_n( ev_RTI->cfi[0],  4*2, cfi[0]);
    mphe = ev_RTI->mphe;
    theta = ev_RTI->theta;
    phi = ev_RTI->phi;
    r = ev_RTI->r;
    zenith = ev_RTI->zenith;
    glat = ev_RTI->glat;
    glong = ev_RTI->glong;
    nev = ev_RTI->nev;
    nerr = ev_RTI->nerr;
    ntrig = ev_RTI->ntrig;
    nhwerr = ev_RTI->nhwerr;
    npart = ev_RTI->npart;
    copy_n(  ev_RTI->nl1l9[0], 2*2, nl1l9[0]);
    copy_n(  ev_RTI->dl1l9[0], 2*3, dl1l9[0]);
    mtrdh = ev_RTI->mtrdh;
    good = ev_RTI->good;
    utime = ev_RTI->utime;
    copy_n(ev_RTI->usec, 2, usec);

#endif

  }
  
  return;
}

#endif