rootana/src/main.cpp File Reference

#include <iostream>
#include <string>
#include <map>
#include "SetupNavigator.h"
#include "definitions.h"
#include "TTree.h"
#include "TBranch.h"
#include "TFile.h"
#include "CommandLine.h"
#include "ModulesFactory.h"
#include "ModulesNavigator.h"
#include "BaseModule.h"
#include "TGlobalData.h"
#include "TSetupData.h"
#include "TAnalysedPulse.h"
#include "TDetectorPulse.h"
#include "TMuonEvent.h"
#include "TAnalysedPulseMapWrapper.h"
#include "EventNavigator.h"
#include "debug_tools.h"

Go to the source code of this file.

SourceAnalPulseMap	gAnalysedPulseMap
SourceDetPulseMap	gDetectorPulseMap
MuonEventList	gMuonEvents
Int_t	Main_event_loop (TTree *dataTree, ARGUMENTS &arguments)
void	ClearGlobalData (TGlobalData *)
TTree *	GetTree (TFile *inFile, const char *t_name)
Int_t	PrepareAnalysedPulseMap (TFile *fileOut)
Int_t	PrepareSingletonObjects (const ARGUMENTS &)
void	Finish (TFile *)
int	main (int argc, char **argv)
void	PrintSetupData (TSetupData *s_data)

---

Function Documentation

void ClearGlobalData

(

TGlobalData *

data

)

Definition at line 203 of file main.cpp.

References gAnalysedPulseMap, gDetectorPulseMap, gMuonEvents, and modules::navigator::Instance().

Referenced by LoopSequence::Process().

{
  // TODO: this comment is now gibberish.
  // We could put this into TGlobalData::Clear(), but we need
  // to be sure that's okay at the alcapana level. That is, if
  // Clear() deletes the TPulseIsland objects, but g_event doesn't
  // own the pulses, they would be deleted later. A solution is to
  // be sure that TGlobalData isn't called in alcapana, or ensure
  // that g_event owns the pulse islands at that level.  
//
//---  LoopSequence::Process now takes care of this part by invoking TGlobalData::Clear
//  StringPulseIslandMap::iterator mapIter;
//  StringPulseIslandMap::iterator mapEnd = data->fPulseIslandToChannelMap.end();
//  for(mapIter = data->fPulseIslandToChannelMap.begin(); mapIter != mapEnd; mapIter++) {
//    // The iterator is pointing to a pair<string, vector<TPulseIsland*> >
//    std::vector<TPulseIsland*>& pulse_vector= mapIter->second;
//    for(size_t i=0; i<pulse_vector.size(); i++){
//      delete pulse_vector[i];
//      pulse_vector[i] = NULL;
//    }
//    pulse_vector.clear();
//  }

    // we should not clear TAPs if we're loading them in from a previously made
    // TAP tree (which uses a TClonesArray that is responsible for deleting them )
    static bool should_delete_TAPS=!modules::navigator::Instance()->GetModule("LoadPulses");
    for(SourceAnalPulseMap::iterator mapIter=gAnalysedPulseMap.begin();
            mapIter != gAnalysedPulseMap.end(); ++mapIter) {

        // The iterator is pointing to a pair<string, vector<TPulseIsland*> >
        AnalysedPulseList& pulse_vector= mapIter->second;
        if(should_delete_TAPS){
            for(size_t i=0; i<pulse_vector.size(); i++){
                delete pulse_vector[i];
                pulse_vector[i] = NULL;
            }
        }
    pulse_vector.clear();
  }
  //  gAnalysedPulseMap.clear();

  for(SourceDetPulseMap::iterator mapIter = gDetectorPulseMap.begin();
          mapIter != gDetectorPulseMap.end(); ++mapIter) {
    // The iterator is pointing to a pair<string, vector<TPulseIsland*> >
    DetectorPulseList& pulse_vector= mapIter->second;
    for(size_t i=0; i<pulse_vector.size(); i++){
      delete pulse_vector[i];
      pulse_vector[i] = NULL;
    }
    pulse_vector.clear();
  }
  //gDetectorPulseMap.clear();

  for(MuonEventList::iterator muonEvent = gMuonEvents.begin();
          muonEvent != gMuonEvents.end(); ++muonEvent) {
      delete *muonEvent;
  }
  gMuonEvents.clear();
}

void Finish

(

TFile *

fileOut

)

Definition at line 169 of file main.cpp.

References SetupNavigator::Close(), EventNavigator::Close(), SetupNavigator::Instance(), and EventNavigator::Instance().

Referenced by main().

                           {
  // and finish up
  fileOut->cd();
  //gAnalysedPulseTree->Write();
  fileOut->Write();
  fileOut->Close();
  EventNavigator::Instance().Close();
  SetupNavigator::Instance()->Close();
}

TTree* GetTree	(	TFile *	inFile,
		const char *	t_name
	)

int main	(	int	argc,
		char **	argv
	)

Definition at line 76 of file main.cpp.

References analyze_command_line(), SetupNavigator::CacheCalibDB(), EventNavigator::ConnectInput(), EventNavigator::DumpRawTree(), Finish(), EventNavigator::GetRawTree(), ARGUMENTS::infile, EventNavigator::Instance(), modules::navigator::Instance(), SetupNavigator::Instance(), Main_event_loop(), ARGUMENTS::mod_file, ARGUMENTS::outfile, PrepareAnalysedPulseMap(), PrepareSingletonObjects(), and SetupNavigator::SetCommandLineArgs().

{
  //load_config_file("MODULES.txt");

  // Parse the command line
  ARGUMENTS arguments;
  int ret = analyze_command_line (argc, argv,arguments);
  if(ret!=0) return ret;
  std::cout << "Starting event" << std::endl;
  // Give the command line arguments to the navigator
  SetupNavigator::Instance()->SetCommandLineArgs(arguments);
  SetupNavigator::Instance()->CacheCalibDB();

  // Read in the configurations file
  ret= modules::navigator::Instance()->LoadConfigFile(arguments.mod_file.c_str());
  if(ret!=0) {
     std::cout<<"Error: Problem loading MODULES file"<<std::endl;
     return ret;
  }

  // Open the input tree file
  EventNavigator& navi = EventNavigator::Instance();
  if ( ! navi.ConnectInput(arguments.infile.c_str()) ){
    std::cout << "Failed to open input file '"
              << arguments.infile <<"'.  Exiting." << std::endl;
    return 1;
  }
  if( modules::navigator::Instance()->DumpInputFile()){
    navi.DumpRawTree();
  }
    
  // Make an initial call to singleton objects that are very likely to be called at some point.
  ret = PrepareSingletonObjects(arguments);
  if(ret!=0) {
     std::cout<<"Error: Problem creating a singleton object..."<<std::endl;
     return ret;
  }

  // Now that we've loaded the TSetupData for this run check if there are any
  // suggested replacements for the wiremap data
  //CheckSetupData(navi.GetSetupData(), arguments.correction_file);

  //Event Tree
  TTree* eventTree = navi.GetRawTree(); 

  // Let's open the output file for analysis data, histograms and so on.
  TFile *fileOut = new TFile(arguments.outfile.c_str(), "RECREATE");
  if(!fileOut->IsOpen()){
    std::cout << "Could not open ROOT output file " 
              << arguments.outfile << std::endl;
    return 1;
  }
  modules::navigator::Instance()->SetOutFile(fileOut);

  // Now let's setup all the analysis modules we want
  // NOTE: This has to be done after the output file was opened else the
  // modules wont have a directory to store things to.
  ret= modules::navigator::Instance()->MakeModules();
  if(ret!=0) {
    std::cout << "Problem creating analysis modules." << std::endl;
    return ret;
  }
  
  // Setup the analysed pulse map to store / read in the pulses
  ret = PrepareAnalysedPulseMap(fileOut);
  if(ret!=0) {
    std::cout << "Problem preparing analysed pulse map." << std::endl;
    return ret;
  }

  // Finally let's do the main loop
  fileOut->cd();
  try {
    Main_event_loop(eventTree,arguments);
  }
  catch (std::exception& e){
    // Attempt to finish the run but ignore further exceptions
    try{ Finish(fileOut); }catch(...){}
    std::cout << "Terminating due to unexpected exception:\n" 
              << e.what() << std::endl; 
    throw;
  }
  catch (...){
    // Attempt to finish the run but ignore further exceptions
    try{ Finish(fileOut); }catch(...){}
    std::cout << "Terminating due to unknown exception" << std::endl;
    throw;
  }
  
  Finish(fileOut);
  return 0;
}

Int_t PrepareAnalysedPulseMap

(

TFile *

fileOut

)

Definition at line 292 of file main.cpp.

Referenced by main().

{
  // TAnalysedMapWrapper
  //gAnalysedPulseMapWrapper = new TAnalysedPulseMapWrapper(gAnalysedPulseMap);
  
  int split = 1;
  //int bufsize = 64000;
  Int_t branchstyle = 1;
  
  if (split < 0) {branchstyle = 0; split = -1-split;}
  TTree::SetBranchStyle(branchstyle);
  
  //gAnalysedPulseTree = new TTree("AnalysedPulseTree", "AnalysedPulseTree");
  //gAnalysedPulseTree->SetAutoSave(100000000);
  //gAnalysedPulseTree->SetMaxVirtualSize(100000000);
  //gAnalysedPulseBranch = gAnalysedPulseTree->Branch("AnalysedPulseMapWrapper", 
  //                                                  "TAnalysedPulseMapWrapper",
  //                                                  &gAnalysedPulseMapWrapper,
  //                                                  bufsize, split);
  //gAnalysedPulseBranch->SetAutoDelete(kFALSE);
  return 0;
}

Int_t PrepareSingletonObjects

(

const ARGUMENTS &

)

Definition at line 316 of file main.cpp.

References EventNavigator::Instance(), and modules::navigator::Instance().

Referenced by main().

                                               {
  // Set up the modules navigator
  if(!modules::navigator::Instance()) return 1;
  
  // Make sure the instance of TSetupData get's loaded 
  // (though in practice we can do this lazily)
  if( !EventNavigator::Instance().GetSetupData() ) return 2;
  
  return 0;
}

void PrintSetupData

(

TSetupData *

s_data

)

Definition at line 265 of file main.cpp.

References TSetupData::fBankToDetectorMap, TSetupData::GetClockTick(), TSetupData::GetDetectorName(), TSetupData::GetNBits(), TSetupData::GetPedestal(), TSetupData::GetTimeShift(), and TSetupData::GetTriggerPolarity().

{
  if(!s_data) return;
  // print things out
  std::map<std::string, std::string>::iterator it_info;
  std::cout << "### TSetupData ###\n"
            << "Bank  Detector Name  Clock Tick (ns)  Pedestal (ADC)"
            << "  Trigger Polarity  Time Shift (ns)  No. of Bits\n"
            << "----  -------------  ---------------  --------------"
            << "  ----------------  ---------------  -----------" 
            << std::endl;

  for (it_info=s_data->fBankToDetectorMap.begin(); 
       it_info!=s_data->fBankToDetectorMap.end(); ++it_info){
    const std::string& bankname = it_info->first;
    std::cout << bankname << "  " 
              << s_data->GetDetectorName(bankname) << "\t\t" 
              << s_data->GetClockTick(bankname) << "\t" 
              << s_data->GetPedestal(bankname) << "\t\t" 
              << s_data->GetTriggerPolarity(bankname) << "\t\t" 
              << s_data->GetTimeShift(bankname) << "\t\t" 
              << s_data->GetNBits(bankname)
              << std::endl;
  }
}

---

Variable Documentation

SourceAnalPulseMap gAnalysedPulseMap

Definition at line 72 of file main.cpp.

Referenced by MakeDetectorPulses::BeforeFirstEntry(), MakeAnalysedPulses::BeforeFirstEntry(), SavePulses::BeforeFirstEntry(), PlotTAP_TDiff::BeforeFirstEntry(), PlotTAP_selfCorrelation::BeforeFirstEntry(), PlotTAP_Integral::BeforeFirstEntry(), PlotAmpVsMuScTDiff::BeforeFirstEntry(), CalcCoarseTimeOffset::BeforeFirstEntry(), PlotIntegralRatios::BeforeFirstEntry(), LoadPulses::BeforeFirstEntry(), ClearGlobalData(), MakeDetectorPulses::ProcessEntry(), MakeAnalysedPulses::ProcessEntry(), SavePulses::ProcessEntry(), PulseViewer::ProcessEntry(), PlotTAP_Time::ProcessEntry(), PlotTAP_TDiff::ProcessEntry(), PlotTAP_selfCorrelation::ProcessEntry(), PlotTAP_Integral::ProcessEntry(), PlotTAP_Energy::ProcessEntry(), PlotTAP_Amplitude::ProcessEntry(), PlotAmpVsTDiff::ProcessEntry(), PlotAmpVsMuScTDiff::ProcessEntry(), CalcCoarseTimeOffset::ProcessEntry(), PlotIntegralRatios::ProcessEntry(), and LoadPulses::ProcessEntry().

SourceDetPulseMap gDetectorPulseMap

Definition at line 73 of file main.cpp.

Referenced by MakeDetectorPulses::BeforeFirstEntry(), PlotTDPs::BeforeFirstEntry(), ClearGlobalData(), MakeSiliconEvents::GetTDPSource(), MakeMuonEvents::ProcessEntry(), MakeDetectorPulses::ProcessEntry(), PlotTDPs::ProcessEntry(), and PlotTDP_TDiff::ProcessEntry().

MuonEventList gMuonEvents

Definition at line 74 of file main.cpp.

Referenced by ClearGlobalData(), MakeSiliconEvents::ProcessEntry(), MakeMuonEvents::ProcessEntry(), PlotTME_activeSiRmuStops::ProcessEntry(), CheckTMEs::ProcessEntry(), TMEViewer::ProcessEntry(), TME_SiArmTDiff::ProcessEntry(), and TME_Al50_EvdE::ProcessEntry().