onlinemonitor/framework/T1290Data.cxx (r4864/r3205)

Go to the documentation of this file.

#include "T1290Data.h"

#include <stdio.h>
#include <stdlib.h>
#include <math.h>

#include "TH1.h"
#include "TH2.h"

#include "TGo4EventProcessor.h"

void T1290Data::Clear()
{
   for (int n=0;n<T1290Data::NumChan;n++) {
      lead[n] = -1;
      trail[n] = -1;
      width[n] = -1;

      hit_lead[n] = 0;
      hit_trail[n] = 0;
      for (int m=0;m<T1290Data::MaxHits;m++) {
         lead_multi[n][m] = -1;
         trail_multi[n][m] = -1;
      }
   }
   extended_trigger_time = -1;
}


void T1290Data::CalcWidth()
{
   for (int n=0;n<NumChan;n++)
      width[n] = abs(trail[n]-lead[n]);
}

int T1290Data::FindGeo(int* pdata, unsigned len)
{
   for (unsigned iCha = 0; iCha < len; iCha++) {
      Int_t data32 = pdata[iCha];
      Int_t type_word = (data32 >> 27) & 0x1f; //data type
      if (type_word==8) return (data32 & 0x1f);
   }

   return 0;
}


void T1290Data::Unpack(int* pdata, unsigned len, int expected_geo)
{
   if (len<2) {
      printf("Error unpack V1290 data - too few messages = %u\n", len);
   }

   for (int n=0;n<NumChan;n++) {
      hit_lead[n] = 0;
      hit_trail[n] = 0;
   }

   Int_t geo(0);

   for (unsigned iCha = 0; iCha < len; iCha++) {
      Int_t data32 = pdata[iCha];
      Int_t type_word = (data32 >> 27) & 0x1f; //data type

//      printf("Msg:%2u Type:%d ", iCha, type_word);

      switch (type_word)  {
         case 8: {  //global_header
            // Int_t event_counter = (data32 & 0x07ffffff) >> 5;
            geo = data32 & 0x1f;
//            printf("Found GEO %d", geo);

            // SL: check that geo corresponds to expected number
            if ((expected_geo>0) && (expected_geo!=geo))
               printf("GEO mismatch, expected %d, is %d\n", expected_geo, geo);

            break;
         }
         case 16: { //Global trailer
            // Int_t word_count = (data32 & 0x001fffe0) >> 5;

            CalcWidth();

            static int errcnt = 0;

            if (iCha+1!=len) {
               if (errcnt++<3)
                  printf("Wrong V1290 data format - more messages after globar trailer\n");
               return;
            }

            break;
         }
         case 1: { //TDC Header
            // Int_t tdc = (data32 >> 24) & 0x03;
            // Int_t event_id = (data32 & 0x00fff000) >> 12;
            // Int_t bunch_id = data32 &0x00000fff;
            break;
         }
         case 3: { //TDC Trailer
            // Int_t tdc = (data32 & 0x03000000) >> 24;
            // Int_t event_id = (data32 & 0x00fff000) >> 12;
            // Int_t word_count = data32 &0x00000fff;
            break;
         }
         case 4: {  //TDC Error
            printf("??TDC ERROR  GEO = %d??\n", geo);
            break;
         }
         case 0: {  //TDC Data
            Int_t TDC_input = (data32 >> 21) & 0x1f;
            Int_t ch_data = data32 & 0x1fffff;
            Int_t front = (data32 >> 26) & 0x1;

            if (front==0) {
               if (hit_lead[TDC_input] == 0) {
                  lead[TDC_input] = ch_data;
               }
               if( hit_lead[ TDC_input] <MaxHits )
                  lead_multi[TDC_input][ hit_lead[ TDC_input] ] = ch_data;
               hit_lead[TDC_input] ++;

            } else {
               if (hit_trail[TDC_input] == 0) {
                  trail[TDC_input] = ch_data;
               }
               if(hit_trail[ TDC_input] <MaxHits )
                  trail_multi[TDC_input][hit_trail[ TDC_input] ] = ch_data;
               hit_trail[TDC_input] ++;
            }
            break;
         }
         case 17: { //Global trigger Time tag
            extended_trigger_time = data32 & 0x7ffffff;
            break;
         }
         default: {
            printf("GEO:%d Unknown tag %d\n", geo, type_word);
            break;
         }
      }
//      printf("\n");
   }
}

// =============================
T1290Proc::T1290Proc()
{
   fTDC_chan = 0;
   for (int n=0;n<T1290Data::NumChan;n++) {
      fTDC_lead[n] = 0;
      fTDC_trail[n] = 0;
      fTDC_width[n] = 0;
   }
}

T1290Proc::~T1290Proc()
{
}

void T1290Proc::MakeHistos(TGo4EventProcessor* proc, const char* name)
{
   TString obname, obtitle;

   obname.Form("%s/%s_channels", name, name);
   obtitle.Form("Signal counts (L & T in same event) in each channel of %s", name);
   fTDC_chan = proc->MakeTH1('I', obname.Data(), obtitle.Data(), T1290Data::NumChan, 0, T1290Data::NumChan, "channel");
   obname.Form("%s/%s_channelsLeading", name, name);
   obtitle.Form("Leading edge counts in each channel of %s", name);
   fTDC_chanLead = proc->MakeTH1('I', obname.Data(), obtitle.Data(), T1290Data::NumChan, 0, T1290Data::NumChan, "channel; Leading edges");
   obname.Form("%s/%s_channelsTrailing", name, name);
   obtitle.Form("Trailing edge counts in each channel of %s", name);
   fTDC_chanTrail = proc->MakeTH1('I', obname.Data(), obtitle.Data(), T1290Data::NumChan, 0, T1290Data::NumChan, "channel; Trailing edges");

  for (int ch=0;ch<T1290Data::NumChan;ch++) {
     obname.Form("%s/Lead/%s_%02d_Lead", name, name, ch);
     obtitle.Form("Leading edge of %s channel %d", name, ch);
     fTDC_lead[ch] = proc->MakeTH1('I', obname.Data(), obtitle.Data(), 5000, 0, 50000, "time");
     obname.Form("%s/Trail/%s_%02d_Trail", name, name, ch);
     obtitle.Form("Trailing edge of %s channel %d", name, ch);
     fTDC_trail[ch] = proc->MakeTH1('I', obname.Data(), obtitle.Data(), 5000, 0, 50000, "time");
     obname.Form("%s/Width/%s_%02d_Width", name, name, ch);
     obtitle.Form("Width of %s channel %d", name, ch);
     fTDC_width[ch] = proc->MakeTH1('I', obname.Data(), obtitle.Data(), 5000, 0, 50000, "time");
  }

  obname.Form("%s/Multi/%s_channelsLeadingMulti", name, name);
  obtitle.Form("Nb multiple Leading edge counts in each channel of %s", name);
  fTDC_chanLeadMulti = proc->MakeTH2('I', obname.Data(), obtitle.Data(), T1290Data::NumChan, 0, T1290Data::NumChan,
                                      16, 0., 16., "channel; Leading edges");
  obname.Form("%s/Multi/%s_channelsTrailingMulti", name, name);
  obtitle.Form("Nb multiple Trailing edge counts in each channel of %s", name);
  fTDC_chanTrailMulti = proc->MakeTH2('I', obname.Data(), obtitle.Data(), T1290Data::NumChan, 0, T1290Data::NumChan,
                                  16, 0., 16., "channel; Trailing edges");
}
void T1290Proc::MakeDebugHistos(TGo4EventProcessor* proc, const char* name)
{
   TString obname, obtitle;

   for (int ch=0;ch<T1290Data::NumChan;ch++) {

      obname.Form("%s/Multi/LeadNb/%s_%02d_Lead_nb", name, name, ch);
      obtitle.Form("nb multi lead vs Leading edge of %s channel %d", name, ch);
      fTDC_lead_nblead[ch] = proc->MakeTH2('I', obname.Data(), obtitle.Data(), 5000, 0, 50000, 16, 0., 16.,"time; Leading edges");

      obname.Form("%s/Multi/Lead/%s_%02d_Lead_Multi", name, name, ch);
      obtitle.Form("Second and later hits Leading edge of %s channel %d", name, ch);
      fTDC_lead_multi[ch] = proc->MakeTH1('I', obname.Data(), obtitle.Data(), 5000, 0, 50000, "time");
      obname.Form("%s/Multi/Trail/%s_%02d_Trail_Multi", name, name, ch);
      obtitle.Form("Second and later hits Trailing edge of %s channel %d", name, ch);
      fTDC_trail_multi[ch] = proc->MakeTH1('I', obname.Data(), obtitle.Data(), 5000, 0, 50000, "time");
      obname.Form("%s/Multi/Width/%s_%02d_Width_Multi", name, name, ch);
      obtitle.Form("Second and later hits Width of %s channel %d", name, ch);
      fTDC_width_multi[ch] = proc->MakeTH1('I', obname.Data(), obtitle.Data(), 5000, -25000, 25000, "time");
   }
}


void T1290Proc::FillHistos(T1290Data& data)
{
   if (fTDC_chan==0) return;

   for (int ch=0; ch<T1290Data::NumChan; ch++) {
      if (data.lead[ch]>0) {
         fTDC_lead[ch]->Fill(data.lead[ch]);
         fTDC_chanLead->Fill(ch);
      }
      if (data.trail[ch]>0) {
         fTDC_trail[ch]->Fill(data.trail[ch]);
         fTDC_chanTrail->Fill(ch);
      }
      if( data.lead[ch]>0 && data.trail[ch]>0)
         fTDC_chan->Fill(ch);
      if (data.width[ch]>0) {
         fTDC_width[ch]->Fill(data.width[ch]);
      }
      fTDC_chanLeadMulti->Fill(ch, data.hit_lead[ch]);
      fTDC_chanTrailMulti->Fill(ch, data.hit_trail[ch]);
   }
}
void T1290Proc::FillDebugHistos(T1290Data& data)
{
   if (fTDC_chan==0) return;

   for (int ch=0; ch<T1290Data::NumChan; ch++) {
      if (data.lead[ch]>0)
         fTDC_lead_nblead[ch]->Fill(data.lead[ch], data.hit_lead[ch]);
      for(int iMultiIndex = 0; iMultiIndex < data.hit_lead[ch]-1 && iMultiIndex<T1290Data::MaxHits; iMultiIndex++ )
      {
         fTDC_lead_multi[ch]->Fill(data.lead_multi[ch][iMultiIndex]);
         if( iMultiIndex < data.hit_trail[ch] -1)
            fTDC_width_multi[ch]->Fill(data.trail_multi[ch][iMultiIndex] - data.lead_multi[ch][iMultiIndex] );
      }
      for(int iMultiIndex = 0; iMultiIndex < data.hit_trail[ch]-1 && iMultiIndex<T1290Data::MaxHits; iMultiIndex++ )
         fTDC_trail_multi[ch]->Fill(data.trail_multi[ch][iMultiIndex]);
   }
}

// Functions to allow for finer control on histograms
void T1290Proc::MakeHistos(TGo4EventProcessor* proc, Bool_t bWithDebugHistos, Bool_t bWithTriggerTime,
                             const char* name)
{
   TString obname, obtitle;

   obname.Form("CAEN/%s_channels", name);
   obtitle.Form("Signal counts (L & T in same event) in each channel of %s", name);
   fTDC_chan = proc->MakeTH1('I', obname.Data(), obtitle.Data(), T1290Data::NumChan, 0, T1290Data::NumChan, "channel");
   obname.Form("CAEN/%s_channelsLeading", name);
   obtitle.Form("Leading edge counts in each channel of %s", name);
   fTDC_chanLead = proc->MakeTH1('I', obname.Data(), obtitle.Data(), T1290Data::NumChan, 0, T1290Data::NumChan, "channel; Leading edges");
   obname.Form("CAEN/%s_channelsTrailing", name);
   obtitle.Form("Trailing edge counts in each channel of %s", name);
   fTDC_chanTrail = proc->MakeTH1('I', obname.Data(), obtitle.Data(), T1290Data::NumChan, 0, T1290Data::NumChan, "channel; Trailing edges");

  for (int ch=0;ch<T1290Data::NumChan;ch++) {
     obname.Form("CAEN/%s/Lead/%s_%02d_Lead", name, name, ch);
     obtitle.Form("Leading edge of %s channel %d", name, ch);
     fTDC_lead[ch] = proc->MakeTH1('I', obname.Data(), obtitle.Data(), 5000, 0, 50000, "time");
     obname.Form("CAEN/%s/Trail/%s_%02d_Trail", name, name, ch);
     obtitle.Form("Trailing edge of %s channel %d", name, ch);
     fTDC_trail[ch] = proc->MakeTH1('I', obname.Data(), obtitle.Data(), 5000, 0, 50000, "time");
     obname.Form("CAEN/%s/Width/%s_%02d_Width", name, name, ch);
     obtitle.Form("Width of %s channel %d", name, ch);
     fTDC_width[ch] = proc->MakeTH1('I', obname.Data(), obtitle.Data(), 5000, 0, 50000, "time");

     if( kTRUE == bWithDebugHistos ) {
        obname.Form("CAEN/%s/Multi/LeadNb/%s_%02d_Lead_nb", name, name, ch);
        obtitle.Form("nb multi lead vs Leading edge of %s channel %d", name, ch);
        fTDC_lead_nblead[ch] = proc->MakeTH2('I', obname.Data(), obtitle.Data(), 5000, 0, 50000, 16, 0., 16.,"time; Leading edges");

        obname.Form("CAEN/%s/Multi/Lead/%s_%02d_Lead_Multi", name, name, ch);
        obtitle.Form("Second and later hits Leading edge of %s channel %d", name, ch);
        fTDC_lead_multi[ch] = proc->MakeTH1('I', obname.Data(), obtitle.Data(), 5000, 0, 50000, "time");
        obname.Form("CAEN/%s/Multi/Trail/%s_%02d_Trail_Multi", name, name, ch);
        obtitle.Form("Second and later hits Trailing edge of %s channel %d", name, ch);
        fTDC_trail_multi[ch] = proc->MakeTH1('I', obname.Data(), obtitle.Data(), 5000, 0, 50000, "time");
        obname.Form("CAEN/%s/Multi/Width/%s_%02d_Width_Multi", name, name, ch);
        obtitle.Form("Second and later hits Width of %s channel %d", name, ch);
        fTDC_width_multi[ch] = proc->MakeTH1('I', obname.Data(), obtitle.Data(), 5000, -25000, 25000, "time");
     } // if( kTRUE == bWithDebugHistos )
  } // for (int ch=0;ch<T1290Data::NumChan;ch++)

  obname.Form("CAEN/%s/Multi/%s_channelsLeadingMulti", name, name);
  obtitle.Form("Nb multiple Leading edge counts in each channel of %s", name);
  fTDC_chanLeadMulti = proc->MakeTH2('I', obname.Data(), obtitle.Data(), T1290Data::NumChan, 0, T1290Data::NumChan,
                                      16, 0., 16., "channel; Leading edges");
  obname.Form("CAEN/%s/Multi/%s_channelsTrailingMulti", name, name);
  obtitle.Form("Nb multiple Trailing edge counts in each channel of %s", name);
  fTDC_chanTrailMulti = proc->MakeTH2('I', obname.Data(), obtitle.Data(), T1290Data::NumChan, 0, T1290Data::NumChan,
                                  16, 0., 16., "channel; Trailing edges");

  if( kTRUE == bWithTriggerTime )
  {
     obname.Form("CAEN/%s/%s_TriggerTimeHigh", name, name);
     obtitle.Form("Higher bits for Trigger time from extended trigger time tag for CAEN %s", name);
     fTDC_CAEN_triggertimeHigh = proc->MakeTH1('I', obname.Data(), obtitle.Data(), 2047, 0, 134152192, "Bin []");
     obname.Form("CAEN/%s/%s_TriggerTimeLow", name, name);
     obtitle.Form("Lower bits for Trigger time from extended trigger time tag for CAEN %s", name);
     fTDC_CAEN_triggertimeLow = proc->MakeTH1('I', obname.Data(), obtitle.Data(), 0xFFFF, 0, 0xFFFF, "Bin []");
  }

}
void T1290Proc::FillHistos(T1290Data& data, Bool_t bWithDebugHistos, Bool_t bWithTriggerTime)
{
   if (fTDC_chan==0) return;

   for (int ch=0; ch<T1290Data::NumChan; ch++) {
      if (data.lead[ch]>0) {
         fTDC_lead[ch]->Fill(data.lead[ch]);
         fTDC_chanLead->Fill(ch);
      }
      if (data.trail[ch]>0) {
         fTDC_trail[ch]->Fill(data.trail[ch]);
         fTDC_chanTrail->Fill(ch);
      }
      if( data.lead[ch]>0 && data.trail[ch]>0)
         fTDC_chan->Fill(ch);
      if (data.width[ch]>0) {
         fTDC_width[ch]->Fill(data.width[ch]);
      }
      fTDC_chanLeadMulti->Fill(ch, data.hit_lead[ch]);
      fTDC_chanTrailMulti->Fill(ch, data.hit_trail[ch]);

      if( kTRUE == bWithDebugHistos ) {
         if (data.lead[ch]>0)
            fTDC_lead_nblead[ch]->Fill(data.lead[ch], data.hit_lead[ch]);
         for(int iMultiIndex = 1; iMultiIndex < data.hit_lead[ch] && iMultiIndex<T1290Data::MaxHits; iMultiIndex++ )
         {
            fTDC_lead_multi[ch]->Fill(data.lead_multi[ch][iMultiIndex]);
            if( iMultiIndex < data.hit_trail[ch] -1)
               fTDC_width_multi[ch]->Fill(data.trail_multi[ch][iMultiIndex] - data.lead_multi[ch][iMultiIndex] );
         }
         for(int iMultiIndex = 1; iMultiIndex < data.hit_trail[ch] && iMultiIndex<T1290Data::MaxHits; iMultiIndex++ )
            fTDC_trail_multi[ch]->Fill(data.trail_multi[ch][iMultiIndex]);
      } // if( kTRUE == bWithDebugHistos ) {
   } // for (int ch=0; ch<T1290Data::NumChan; ch++)

   if( kTRUE == bWithTriggerTime && -1 < data.extended_trigger_time )
   {
      fTDC_CAEN_triggertimeHigh->Fill( data.extended_trigger_time ) ;
      fTDC_CAEN_triggertimeLow->Fill( data.extended_trigger_time % 0x10000  ) ;
   }
}