testAnalysis.cxx

Go to the documentation of this file.

#include <iostream>
#include <iomanip>
#include <string>

#include "TFile.h"
#include "TChain.h"
#include "TTree.h"

#include "TGriffin.h"
#include "TChannel.h"
#include "TGRSIOptions.h"
#include "TGRSIMnemonic.h"

int main(int argc, char** argv)
{
   if(argc < 3) {
      std::cerr<<"Usage: "<<argv[0]<<" <entries (-1 = all)> <file(s)>"<<std::endl;
      return 1;
   }

   TChain chain("AnalysisTree");

   for(int i = 2; i < argc; ++i) {
      chain.Add(argv[i]);
   }

   TGriffin* griffin = nullptr;

   if(chain.SetBranchAddress("TGriffin", &griffin) == TTree::kMissingBranch) {
      griffin = new TGriffin;
   }

   int64_t entries = chain.GetEntries();
   int64_t maxEntries = std::stoi(std::string(argv[1]));

   if(maxEntries == -1 || maxEntries > entries) maxEntries = entries;

   static TFile* currentFile = nullptr;
   TChannel::SetMnemonicClass(TGRSIMnemonic::Class());
   for(int64_t e = 0; e < maxEntries; ++e) {
      chain.GetEntry(e);
      if(currentFile != chain.GetCurrentFile()) {
         currentFile = chain.GetCurrentFile();
         TChannel::ReadCalFromFile(currentFile);
         TGRSIOptions::AnalysisOptions()->ReadFromFile(currentFile);
         std::cout<<"Mnemonic class "<<TChannel::GetMnemonicClass().GetClassName()<<std::endl;
      }

      //TChannel::ReadCalFromCurrentFile();

      //if(griffin->GetMultiplicity() == 0) std::cout<<e<<": no GRIFFIN hits"<<std::endl;

      for(int g = 0; g < griffin->GetMultiplicity(); ++g) {
         auto hit = griffin->GetGriffinHit(g);
         std::cout<<"entry "<<e<<", hit "<<g<<std::endl;
         hit->Print("all");
      }
   }

   return 0;
}