LISummaryAnalyser.cxx

Go to the documentation of this file.

//                                                                     
// Program name LISummaryAnalyser.cxx
//                                                                     
// Package: LISummary           
//  
// Coded by Jeff Hartnell Nov/2002   
//                                                                  
// Purpose: 
//                                                                    
// Contact: jeffrey.hartnell@physics.ox.ac.uk   

#include <iostream>

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

#include "JobControl/JobCModuleRegistry.h"
#include "JobControl/JobCommand.h"
#include "MessageService/MsgService.h"
#include "MinosObjectMap/MomNavigator.h"
#include "RawData/RawLIAdcSummary.h"
#include "RawData/RawDaqHeader.h"
#include "RawData/RawLIAdcSummaryBlock.h"
#include "RawData/RawRecord.h"
#include "RawData/RawLIHeader.h" 

#include "LISummary/LISummaryAnalyser.h"

#include <cassert>

ClassImp(LISummaryAnalyser)

JOBMODULE(LISummaryAnalyser,"LISummaryAnalyser","A module to analyse raw LI data");

CVSID("$Id: LISummaryAnalyser.cxx,v 1.1 2003/07/28 11:05:17 hartnell Exp $");

//......................................................................

LISummaryAnalyser::LISummaryAnalyser() 
{
  MSG("LISummaryAnalyser", Msg::kDebug) 
    << "Running constructor..." << endl;
  
  //initialise data members
  calibPoint=-1;
  calibType=-1;
  led=-1;
  mean=-1;
  numEntries=-1;
  period=-1;
  pulseHeight=-1;
  pulserBox=-1;
  pulses=-1;
  pulseWidth=-1;
  rms=-1;
  timeNanoSec=-1; 
  timeSec=-1;

  plottedOnce=0;
  numSummariesProcessed=0;
  msgLevel=0;

  hMeanPb=new TH1F*[NUMPULSERBOXES];
  hMeanCorrelatedPb=new TH1F*[NUMPULSERBOXES];
  for (Int_t i=0;i<NUMPULSERBOXES;i++){
    string s="Trigger PMT Mean ADC (Pb ";
    string pB=Form("%d",i);
    s=s+pB+")";
    hMeanPb[i]=new TH1F(s.c_str(),s.c_str(),20000,-5,17005);
    hMeanPb[i]->GetXaxis()->SetTitle("Mean");
    hMeanPb[i]->GetXaxis()->CenterTitle();
    hMeanPb[i]->GetYaxis()->SetTitle("Num Entries");
    hMeanPb[i]->GetYaxis()->CenterTitle();
    hMeanPb[i]->SetFillColor(0);
    hMeanPb[i]->SetBit(TH1::kCanRebin);

    s="Trigger PMT Mean ADC (Correlated, Pb ";
    s=s+pB+")";
    hMeanCorrelatedPb[i]=new TH1F(s.c_str(),s.c_str(),150,-5,17005);
    hMeanCorrelatedPb[i]->GetXaxis()->SetTitle("Mean");
    hMeanCorrelatedPb[i]->GetXaxis()->CenterTitle();
    hMeanCorrelatedPb[i]->GetYaxis()->SetTitle("Num Entries");
    hMeanCorrelatedPb[i]->GetYaxis()->CenterTitle();
    hMeanCorrelatedPb[i]->SetFillColor(0);
    hMeanCorrelatedPb[i]->SetLineColor(2);
    hMeanCorrelatedPb[i]->SetBit(TH1::kCanRebin);
  }

  hTimestamp=new TH1F("hTimestamp","Timestamp",
                      6000000,1.0365e9,1.0370e9);
  hTimestamp->GetXaxis()->SetTitle("Timestamp");
  hTimestamp->GetXaxis()->CenterTitle();
  hTimestamp->GetYaxis()->SetTitle("Number of entries");
  hTimestamp->GetYaxis()->CenterTitle();
  hTimestamp->SetFillColor(0);
  hTimestamp->SetBit(TH1::kCanRebin);
  
  hTimestampCorrelated=new TH1F("hTimestampCorrelated",
                                "TimestampCorrelated",
                                3000000,1.035e9,1.038e9);
  hTimestampCorrelated->GetXaxis()->SetTitle("TimestampCorrelated");
  hTimestampCorrelated->GetXaxis()->CenterTitle();
  hTimestampCorrelated->GetYaxis()->SetTitle("Number of entries");
  hTimestampCorrelated->GetYaxis()->CenterTitle();
  hTimestampCorrelated->SetFillColor(0);
  hTimestampCorrelated->SetLineColor(2);
  hTimestampCorrelated->SetBit(TH1::kCanRebin);

  cTimestamp=new TCanvas("cTimestamp","Timestamp",
                                  0,0,1200,600);
  cTimestamp->SetFillColor(0);
  cTimestamp->cd();
  hTimestamp->Draw();
  //hTimestampCorrelated->Draw("same");

  MSG("LISummaryAnalyser", Msg::kDebug) 
    << "Finished constructor" << endl;
}

//......................................................................

LISummaryAnalyser::~LISummaryAnalyser()
{
  MSG("LISummaryAnalyser", Msg::kDebug) 
    << "Running destructor..." << endl;

  MSG("LISummaryAnalyser", Msg::kDebug) 
    << "Finished destructor" << endl;
}

//......................................................................

JobCResult LISummaryAnalyser::Ana(const MomNavigator *mom)
{
  JobCResult result(JobCResult::kPassed);

  //check that mom exists.
  assert(mom);
  
  //find raw record in mom
  RawRecord *rawrec = dynamic_cast<RawRecord *>
    (mom->GetFragment("RawRecord"));
  
  if (!rawrec) {
    result.SetError().SetFailed();
    return result;
  }

  //get LI Summary block from rawdatablock
  const RawLIAdcSummaryBlock *rawliadcsummaryblock = 
    dynamic_cast<const RawLIAdcSummaryBlock*>
    (rawrec->FindRawBlock("RawLIAdcSummaryBlock"));

  //check li block exists
  if(!rawliadcsummaryblock){  
    result.SetError().SetFailed();
    return result;
  }

  // Work through rawliadcsummaryblock..  
  timeStamp=rawliadcsummaryblock->GetTimeStamp();
  timeSec=static_cast<Int_t>(timeStamp.GetSec());
  timeNanoSec=static_cast<Int_t>(timeStamp.GetNanoSec());
  pulserBox=rawliadcsummaryblock->GetPulserBox();
  led=rawliadcsummaryblock->GetLed();
  pulseHeight=rawliadcsummaryblock->GetPulseHeight();
  pulses=rawliadcsummaryblock->GetPulses();
     
  if (numSummariesProcessed>msgLevel){
    MSG("LIAnalysis",Msg::kInfo) 
      <<"Number of summaries processed = "<<numSummariesProcessed
      <<endl;
    msgLevel+=500000;
  }
  numSummariesProcessed+=rawliadcsummaryblock->GetNumberOfSummaries();

  //Work through all summaries
  for(Int_t n=0;n<rawliadcsummaryblock->GetNumberOfSummaries();n++){
    const RawLIAdcSummary *rawliadcsummary=
      rawliadcsummaryblock->At(n);
      
    mean=rawliadcsummary->GetMean();
    rms=rawliadcsummary->GetRms();
    numEntries=rawliadcsummary->GetEntries();
    
    Int_t crate=rawliadcsummary->GetChannel().GetCrate();
    Int_t varc=rawliadcsummary->GetChannel().GetVarcId();
    Int_t vfb=rawliadcsummary->GetChannel().GetVaAdcSel();
    Int_t vmm=rawliadcsummary->GetChannel().GetVmm();
    Int_t chip=rawliadcsummary->GetChannel().GetVaChip();

    //only fill if trigger pmt
    //lilookup.Pb2Channel(pulserBox)
    //only plot correct channel
    if (crate==1 && varc==2 && vmm==5 && vfb==0 && chip==1){
      hTimestamp->Fill(timeSec+(1e-9)*timeNanoSec);
      cout<<"seconds="<<timeSec<<", nano="<<timeNanoSec
          <<", sum="<<timeSec+(1e-9)*timeNanoSec
          <<endl;
    }
  }

  //if (crate>1 && plottedOnce==0){

  //plottedOnce=1;
  //}
  cTimestamp->Update();

  return result;
}

//......................................................................

void LISummaryAnalyser::EndJob() 
{

}

//......................................................................

---