ExampleRawAnalysis.cpp

Example of an analysis class for raw (i.e. unreconstructed) data

This is the analysis class generated by default by KaliVedaGUI for raw data analysis.

#pragma once
 
#include "KVRawDataAnalyser.h"
#include <vector>
#include <string>
 
class ExampleRawAnalysis : public KVRawDataAnalyser
{
   int n_sig;
   Int_t Mult = 0;
 
public:
   ExampleRawAnalysis() = default;
 
   void InitAnalysis() override;
   void InitRun() override;
   Bool_t Analysis() override;
   void EndRun() override {}
   void EndAnalysis() override {}
 
   ClassDefOverride(ExampleRawAnalysis,1)//Analysis of raw data
};
 

#include "ExampleRawAnalysis.h"
#include "KVMultiDetArray.h"
 
ClassImp(ExampleRawAnalysis)
 
// This class is derived from the KaliVeda class KVRawDataAnalyser.
// It is to be used for analysis of raw (i.e. unreconstructed) data.
//
// The following member functions are called in turn:
//
//    InitAnalysis():  called at the very beginning of the analysis
//    InitRun():       called everytime a run starts
//    Analysis():      called for each event
//    EndRun():        called everytime a run ends
//    EndAnalysis():   called after all requested data has been read
//
// Modify these methods as you wish in order to create your analysis class.
 
void ExampleRawAnalysis::InitAnalysis()
{
   // Declaration of histograms, trees, etc.
   // Called at the beginning of the analysis
   // The examples given are compatible with interactive and batch analyses.
 
   /*** DECLARING SOME HISTOGRAMS ***/
   AddHisto<TH1F>("Mult", "Number of fired detectors in each event", 1000, -.5, 999.5);
 
   /*** USING A TREE ***/
   // Here we fill a TTree every event with
   //  - the number of fired detector signals
   //  - for each signal:
   //    + detector index
   //    + detector type
   //    + signal type
   //    + signal value
   auto tt = AddTree("raw_data", "");
   tt->Branch("n_sig", &n_sig);
   AddVectorBranch<int>(tt, "det_index");
   AddVectorBranch<std::string>(tt, "det_type");
   AddVectorBranch<std::string>(tt, "sig_type");
   AddVectorBranch<float>(tt, "sig_val");
 
   /*** DEFINE WHERE TO SAVE THE RESULTS ***/
   // This filename will be used for interactive jobs.
   // When running in batch mode, this will automatically use the job name.
   SetJobOutputFileName("ExampleRawAnalysis_results.root");
}
 
//____________________________________________________________________________//
 
void ExampleRawAnalysis::InitRun()
{
   // Initialisation performed at beginning of each runfile. The experimental setup (KVMultiDetArray
   // or KVExpSetUp object) is initialised with the parameters of the run just before calling this method.
   //
   // Some useful global pointers/methods:
   //
   //  - gMultiDetArray : global pointer to experiment setup description;
   //  - gExpDB : global pointer to database (KVExpDB) for current dataset;
   //  - gDataSet : global pointer to current dataset (KVDataSet)
   //  - GetRunNumber() : returns current run number;
   //  - GetRunIndexNumber() : returns full run/index of current runfile;
   //  - GetCurrentRun() : returns pointer to KVDBRun object describing current run in database;
   //  - GetCurrentRunFile() : return pointer to KVDBRunFile object describing current runfile in database;
   //  - GetSystem() : return pointer to KVDBSystem object describing current reaction in database.
   //
   // See parent classes KVRawDataAnalyser, KVDataSetAnalyser and KVDataAnalyser for more.
 
   Info("InitRun", "Beginning analysis of runfile %s containing %llu events", GetRunIndexNumber().as_string().Data(), GetCurrentRunFile()->GetEvents());
   GetCurrentRunFile()->ls();
 
   // Set tree title to name of system being analysed
   GetTree("raw_data")->SetTitle(GetSystem()->GetName());
}
 
//____________________________________________________________________________//
 
Bool_t ExampleRawAnalysis::Analysis()
{
   //Analysis method called for each event
   //
   //  - GetEventNumber() : returns current event number;
   //  - GetRunFileReader() : returns object used to read data (KVRawDataReader child class);
   //  - gMultiDetArray->HandledRawData() : returns kTRUE if interesting data was read from event;
   //
   //  Processing will stop if this method returns kFALSE
 
   if (gMultiDetArray->HandledRawData()) {
 
      ClearVectorBranches();
 
      // get number of fired detectors
      Mult = gMultiDetArray->GetFiredDetectors()->GetEntries();
      // Note - if experimental set up is a combination of two or more detector arrays,
      //        and you are interested in only one of them, you could use
      //           gMultiDetArray->GetArray("[name]")->GetFiredDetectors()
      FillHisto("Mult", Mult);
 
      // get list of fired signals
      auto sigs = gMultiDetArray->GetFiredSignals();
 
      KVDetectorSignal* det_sig;
      TIter it(sigs);
      n_sig=0;
      while ((det_sig = (KVDetectorSignal*)it())) {
         ++n_sig;
         FillVectorBranch("det_index",det_sig->GetDetector()->GetIndex());
         FillVectorBranch("det_type",(std::string)det_sig->GetDetector()->GetLabel());
         FillVectorBranch("sig_type",(std::string)det_sig->GetType());
         FillVectorBranch("sig_val",det_sig->GetValue());
      }
      FillTree();
   }
 
   return kTRUE;
}