pileup_ana_module.cxx

Go to the documentation of this file.

/********************************************************************\

  Name:         pileup_ana_module.c
  Created by:   Thomas Prokscha
  Date:         05-Oct-2005

  Contents:     Analyzer module for LEM pileup data

\********************************************************************/

/*-- Include files -------------------------------------------------*/

/* standard includes */
#include <stdio.h>
#include <time.h>
#include <string.h>
#include <math.h>

/* midas includes */
#include "midas.h"
#include "rmidas.h"
#include "v1190.h"
#include "experim.h"
#include "nemu_experim.h"

#ifdef HAVE_HBOOK
#include <cfortran.h>
#include <hbook.h>
#else
/* root includes */
#include <TH1F.h>
#include <TH2F.h>
#include <TTree.h>
#endif

/*-- ODB Parameters ------------------------------------------------*/
extern RUNINFO          runinfo;
PILEUPANAMODULE_PARAM   pileup_ana_param;
PILEUPANAMODULE_PARAM_STR(pileup_ana_param_str);

/*-- Module declaration --------------------------------------------*/
INT pileup_init(void);
INT pileup_ana(EVENT_HEADER*,void*);
INT pileup_bor(INT run_number);
INT pileup_eor(INT run_number);


ANA_MODULE pileup_ana_module = {
  "PileUpAnaModule",              /* module name           */
  "Thomas Prokscha",              /* author                */
  pileup_ana,                     /* event routine         */
  pileup_bor,                     /* BOR routine           */
  pileup_eor,                     /* EOR routine           */
  pileup_init,                    /* init routine          */
  NULL,                           /* exit routine          */
  &pileup_ana_param,              /* parameter structure   */
  sizeof(pileup_ana_param),       /* structure size        */
  (const char **)pileup_ana_param_str, /* initial parameters    */
};

/*-- other global variables ----------------------------------------*/
#define N_PUHIST  4
#ifndef HAVE_HBOOK
static TH1F* hPUHist[N_PUHIST];
//the following objects are from mana.c
extern TFolder   *gManaHistosFolder;
extern TObjArray *gHistoFolderStack;
#endif

/*-- INIT routine --------------------------------------------------*/
INT pileup_init(void)
{
  extern char runname[256];
  char title[256], name[256];
  int   nbin;
  float xlow, xhigh;

  for (int i=0; i<N_PUHIST; i++){
    sprintf(name,  "hPileUp%02d", i);
    sprintf(title, "%s Run %s_%04d", pileup_ana_param.histotitles.pileuptitles[i], runname, runinfo.run_number);
    nbin  = pileup_ana_param.histobinning.histonbin[i];
    xlow  = pileup_ana_param.histobinning.histoxlow[i];
    xhigh = pileup_ana_param.histobinning.histoxup[i];
    hPUHist[i] = H1_BOOK(name, title, nbin, xlow, xhigh);
  }

  return SUCCESS;
}

/*-- BOR routine ---------------------------------------------------*/

INT pileup_bor(INT run_number)
{
  char title[256];
  extern char runname[256];

  for (int i=0; i<N_PUHIST; i++){
    sprintf(title, "%s Run %s_%04d", pileup_ana_param.histotitles.pileuptitles[i], runname, run_number);
    hPUHist[i]->SetTitle(title);
  }

  return SUCCESS;
}

/*-- EOR routine ---------------------------------------------------*/

INT pileup_eor(INT run_number)
{
  return SUCCESS;
}

/*-- event routine -------------------------------------------------*/
INT pileup_ana(EVENT_HEADER *pheader, void *pevent)
{
  INT       n, ind, counter[N_TDC_CHANNELS];
  DWORD     *ptdc;
  float     tdc_data[N_TDC_CHANNELS][DATA_N_HITS], pu_time;
  LEM_DATA  *lem_data;

  /* look for TDC0 bank */
  n = bk_locate(pevent, "TDC0", &ptdc);
  if (n == 0 )
    return 1;

  memset(&counter,  0x00, sizeof(counter));
  memset(&tdc_data, 0x00, sizeof(tdc_data));

  for (int i=1; i<n; i++){ // i=0 is a control word, including the event length
    lem_data = (LEM_DATA *) &ptdc[i];
    
    switch ( lem_data->channel ){
      case TD:
        ind = TD;
        tdc_data[ind][counter[ind]] = (float) lem_data->data;
        if ( counter[ind] > 0 && counter[ind] < 3 ){ //use 1st and 2nd pile up hits
          pu_time = tdc_data[ind][counter[ind]] - tdc_data[ind][0];
          if ( pu_time < 0) pu_time = pu_time + TDC_BIT_RANGE;
          hPUHist[counter[ind]-1]->Fill(pu_time,1.);
        }
        counter[ind]++;
        break;

      case MCP2F:
        ind = MCP2F;
        tdc_data[ind][counter[ind]] = (float) lem_data->data;
        if ( counter[ind] > 0 && counter[ind] < 3 ){//use 1st and 2nd pile up hits
          pu_time = tdc_data[ind][counter[ind]] - tdc_data[ind][0];
          if ( pu_time < 0) pu_time = pu_time + TDC_BIT_RANGE;
          hPUHist[counter[ind]+1]->Fill(pu_time,1.);
        }
        counter[ind]++;
        break;

      default:
        break;
    }
  }

  return SUCCESS;
}
/* ----------------------------------------------------------------------------*\
        EOF pileup_ana_module.c
\* ----------------------------------------------------------------------------*/