mcp1_ana_module.cxx

Go to the documentation of this file.

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

  Name:         mcp1_ana_module.c
  Created by:   Thomas Prokscha
  Date:         07-Oct-2005

  Contents:     Analyzer module for LEM experiment BC-MCP1-(e+) TOF

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

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

/* standard includes */
#include <stdio.h>
#include <stdlib.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;
extern TRIGGER_SETTINGS trigger_settings;
MCP1ANAMODULE_PARAM     mcp1_ana_param;
MCP1ANAMODULE_PARAM_STR(mcp1_ana_param_str);

/*-- Module declaration --------------------------------------------*/
INT mcp1_init(void);
INT mcp1_ana(EVENT_HEADER*,void*);
INT mcp1_bor(INT run_number);
INT mcp1_eor(INT run_number);

ANA_MODULE mcp1_ana_module = {
  "MCP1AnaModule",                /* module name           */
  "Thomas Prokscha",              /* author                */
  mcp1_ana,                       /* event routine         */
  mcp1_bor,                       /* BOR routine           */
  mcp1_eor,                       /* EOR routine           */
  mcp1_init,                      /* init routine          */
  NULL,                           /* exit routine          */
  &mcp1_ana_param,                /* parameter structure   */
  sizeof(mcp1_ana_param),         /* structure size        */
  (const char **)mcp1_ana_param_str, /* initial parameters    */
};

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

static DWORD dataWindowMuonDecay;  //needed for post-pileup rejection

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

  for (int i=0; i<N_MCP1HIST; i++){
    sprintf(name,  "hMCP1%02d", i);
    sprintf(title, "%s Run %s_%04d", mcp1_ana_param.histotitles.toftitles[i], runname, runinfo.run_number);
    nbin  = mcp1_ana_param.histobinning.histonbin[i];
    xlow  = mcp1_ana_param.histobinning.histoxlow[i];
    xhigh = mcp1_ana_param.histobinning.histoxup[i];
    hMCP1Hist[i] = H1_BOOK(name, title, nbin, xlow, xhigh);
  }

  dataWindowMuonDecay = 0;
  
  return SUCCESS;
}

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

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

  for (int i=0; i<N_MCP1HIST; i++){
    sprintf(title, "%s Run %s_%04d", mcp1_ana_param.histotitles.toftitles[i], runname, run_number);
    hMCP1Hist[i]->SetTitle(title);
  }

  dataWindowMuonDecay = trigger_settings.event_0_settings.data_window; 

  return SUCCESS;
}

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

INT mcp1_eor(INT run_number)
{
  return SUCCESS;
}

/*-- event routine -------------------------------------------------*/
INT mcp1_ana(EVENT_HEADER *pheader, void *pevent)
{
  INT        n, counter[N_TDC_CHANNELS];
  WORD       channel;
  DWORD      *ptdc;
  DWORD      tdc_data[N_TDC_CHANNELS][DATA_N_HITS];
  INT        difftime_last, difftime_now;
  BOOL       foundBC, foundPOS;
  LEM_DATA   *lem_data;

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

  // return if we don't have the NEW_EVENT_FLAG present
  if ( !(ptdc[0] & NEW_EVENT_MASK)) return 1;

  // return if not EVENT_0_TYPE (BC-MCP1)
  if ( !(ptdc[0] & (EVENT_0_TYPE<<16)) ) return 1;

  memset(&counter,  0x00, sizeof(counter));
  memset(&tdc_data, 0x00, sizeof(tdc_data));
  foundBC = foundPOS = FALSE;

  for (int i=1; i<n; i++){ // i=0 is a control word, including the event length
    lem_data = (LEM_DATA *) &ptdc[i];
    channel    = lem_data->channel;
    if ( channel > N_TDC_CHANNELS-1 ) continue; //not understood how this can happen, but it can happen
                                  //after a HV break at the moderator (all channel get "many" hits, like
                                  //a "flash" in the detectors
    tdc_data[channel][counter[channel]] = lem_data->data;
    counter[channel]++;

    switch ( channel ){
      case BCL1: case BCL2: case BCL3: case BCL4: case BCR1: case BCR2: case BCR3: case BCR4:
        foundBC = TRUE;
        break;

      case POSMCP1I: case POSMCP1O:
        foundPOS = TRUE;
        break;

      default:
        break;
    }
  }

  if ( !foundBC && !foundPOS ) return 1; // no BeamCounter or e+ data found


  if ( foundBC ){
    // find the nearest BC hit with respect to MCP1 and fill it in the TOF histogram
    difftime_now  = difftime_last = -2*trigger_settings.tof_bc_window;
    for (int i=BCL1; i<=BCR4; i++){
      if ( counter[i] == 0 ) continue;
      for (int j=0; j<counter[i]; j++){
          difftime_now = tdc_data[MCP1][0] - tdc_data[i][j];
          if ( (difftime_now < 0) && (fabs(difftime_now) < trigger_settings.tof_bc_window) &&
               (difftime_now > difftime_last)) //we have negative differences
             difftime_last = difftime_now;
      }
    }
    difftime_last += trigger_settings.tof_bc_window; //to have "forward" timing

    hMCP1Hist[0]->Fill( (float) difftime_last, 1.);
  }
  
  if ( tdc_data[POSMCP1I][0] > 0 &&  tdc_data[POSMCP1O][0] > 0) {
    difftime_last = tdc_data[POSMCP1I][0] - tdc_data[MCP1][0];
    hMCP1Hist[1]->Fill( (float) difftime_last, 1.);
    
    // check for post-pileup in MCP1: if post-pileup, return
    if ( tdc_data[MCP1][1] != 0 && 
      (tdc_data[MCP1][1] - tdc_data[MCP1][0]) < dataWindowMuonDecay) return SUCCESS;
  
    // check, if there was no 2nd e+ during data window   
    if ( tdc_data[POSMCP1I][1] != 0 ){
      if ( (tdc_data[POSMCP1I][1] - tdc_data[MCP1][0]) < dataWindowMuonDecay &&
          abs(tdc_data[POSMCP1I][1]-tdc_data[POSMCP1O][1]) < 50 ) return SUCCESS; //2nd e+
    }
    hMCP1Hist[2]->Fill( (float) tdc_data[POSMCP1I][0] - tdc_data[MCP1][0], 1.);  
  }
  return SUCCESS;
}
/* ----------------------------------------------------------------------------*\
        EOF mcp1_ana_module.c
\* ----------------------------------------------------------------------------*/