SimpleHistograms Class Reference

#include <SimpleHistograms.h>

Inheritance diagram for SimpleHistograms:

List of all members.

Public Member Functions
	SimpleHistograms (const char *HistogramDirectoryName)
	~SimpleHistograms ()
int	ProcessGenericEntry (TGlobalData *gData, const TSetupData *gSetup)
int	Preprocess (TGlobalData *gData, const TSetupData *gSetup)
	Method called by the main pre-process loop.
int	Postprocess (TGlobalData *gData, const TSetupData *gSetup)
	Method called by the main pre-process loop.
virtual int	BeforeFirstEntry (TGlobalData *gData, const TSetupData *setup)=0
virtual int	AfterLastEntry (TGlobalData *gData, const TSetupData *setup)=0
void	SetAlias (const std::string &alias)
std::string	GetAlias () const
const char *	GetName () const
	Get the name of this module as given to the constructor of the base class.
Protected Member Functions
bool	Debug () const
TDirectory *	GetDirectory () const
	Get the TDirectory for this module.
TDirectory *	GetDirectory (const std::string &name="")
Protected Attributes
TDirectory *	dir
Private Member Functions
virtual int	ProcessEntry (TGlobalData *gData)

---

Detailed Description

Definition at line 7 of file SimpleHistograms.h.

---

Constructor & Destructor Documentation

SimpleHistograms::SimpleHistograms

(

const char *

HistogramDirectoryName

)

Definition at line 22 of file SimpleHistograms.cpp.

References banks, hBankSize, and hNumberOfNoCoincidences.

                                                                     :
  BaseModule(HistogramDirectoryName){

  // The following are the banks we want to compare in this example
  banks.push_back("Nec0");
  banks.push_back("Nfc0");
  banks.push_back("Ngc0");
  banks.push_back("Nhc0");
  banks.push_back("Nae0");
  banks.push_back("Nbe0");
  banks.push_back("Nce0");
  banks.push_back("Nee0");
//  banks.push_back("Nfe0");
//  banks.push_back("Nge0");
//  banks.push_back("Nhe0");
  
  hBankSize = new TH2F("hBankSize", "Bank size",
                       1500, -1.5, 1498.5, banks.size(), 0.5, banks.size()+0.5);

  for(unsigned int i = 1; i<=banks.size(); ++i) hBankSize->GetYaxis()->SetBinLabel(i, banks.at(i-1).data());

  // Let's compare the number of comparisons. If N is the size of the bank vector, then we need 
  // (N-1) + (N-2) +... + 1 comparisons
  int NrComparisons = banks.size() * (banks.size()-1) / 2; 
  
  hNumberOfNoCoincidences = new TH2F("hNumberOfNoCoincidences", "Number of missing coincidences per tree entry",
                                     1000, -1.5, 998.5, NrComparisons, 0.5, NrComparisons+0.5);
}

SimpleHistograms::~SimpleHistograms

(

)

Definition at line 51 of file SimpleHistograms.cpp.

                                   {
}

---

Member Function Documentation

virtual int BaseModule::AfterLastEntry	(	TGlobalData *	gData,
		const TSetupData *	setup
	)			[pure virtual, inherited]

Optional method which is called once after the main event loop Can be used to dump a summary or finalise a histogram

Returns:

0 on sucess and non-zero if a problem occurred

Implemented in TME_Al50_EvdE, TME_SiArmTDiff, TMEViewer, IslandAmplitude, IslandLength, FastSlowCompare, GeSpectrum, LoadPulses, PlotIntegralRatios, CalcCoarseTimeOffset, CheckTMEs, ExportPulse, PlotAmpVsMuScTDiff, PlotAmpVsTDiff, PlotTAP_Amplitude, PlotTAP_Energy, PlotTAP_Integral, PlotTAP_selfCorrelation, PlotTAP_TDiff, PlotTAP_Time, PlotTDP_TDiff, PlotTDPs, PlotTME_activeSiRmuStops, PlotTPI_PedestalAndNoise, PulseViewer, SavePulses, MakeAnalysedPulses, PulseCandidateFinder_InvestigateParameters, TemplateCreator, MakeDetectorPulses, MakeMuonEvents, and MakeSiliconEvents.

Referenced by BaseModule::Postprocess().

virtual int BaseModule::BeforeFirstEntry	(	TGlobalData *	gData,
		const TSetupData *	setup
	)			[pure virtual, inherited]

Optional method which is called once before the main event loop Can be used to parse options and setup histograms.

Returns:

0 on sucess and non-zero if a problem occurred

Implemented in TME_Al50_EvdE, TME_SiArmTDiff, TMEViewer, IslandAmplitude, IslandLength, FastSlowCompare, GeSpectrum, LoadPulses, PlotIntegralRatios, CalcCoarseTimeOffset, CheckTMEs, ExportPulse, PlotAmpVsMuScTDiff, PlotAmpVsTDiff, PlotTAP_Amplitude, PlotTAP_Energy, PlotTAP_Integral, PlotTAP_selfCorrelation, PlotTAP_TDiff, PlotTAP_Time, PlotTDP_TDiff, PlotTDPs, PlotTME_activeSiRmuStops, PlotTPI_PedestalAndNoise, PulseViewer, SavePulses, MakeAnalysedPulses, PulseCandidateFinder_InvestigateParameters, TemplateCreator, MakeDetectorPulses, MakeMuonEvents, and MakeSiliconEvents.

Referenced by BaseModule::Preprocess().

bool BaseModule::Debug

(

)

const [inline, protected, inherited]

Check whether this module was asked to print extra debug information

Returns:

true if 'debug' was given in the modules file, false if not

Definition at line 71 of file BaseModule.h.

References BaseModule::fDebug.

Referenced by MakeAnalysedPulses::AddGenerator(), ExportPulse::AddToExportList(), PulseCandidateFinder_InvestigateParameters::AfterLastEntry(), PlotTPI_PedestalAndNoise::AfterLastEntry(), PlotTDPs::AfterLastEntry(), PlotTDP_TDiff::AfterLastEntry(), IslandLength::AfterLastEntry(), IslandAmplitude::AfterLastEntry(), MakeDetectorPulses::BeforeFirstEntry(), TemplateCreator::BeforeFirstEntry(), PulseCandidateFinder_InvestigateParameters::BeforeFirstEntry(), PlotTPI_PedestalAndNoise::BeforeFirstEntry(), PlotTDPs::BeforeFirstEntry(), PlotTDP_TDiff::BeforeFirstEntry(), IslandLength::BeforeFirstEntry(), IslandAmplitude::BeforeFirstEntry(), PulseViewer::ConsiderDrawing(), MakeAnalysedPulses::MakeAnalysedPulses(), MakeDetectorPulses::MakeGenerator(), ExportPulse::PlotTPI(), MakeDetectorPulses::ProcessEntry(), TemplateCreator::ProcessEntry(), PulseCandidateFinder_InvestigateParameters::ProcessEntry(), MakeAnalysedPulses::ProcessEntry(), PulseViewer::ProcessEntry(), PlotTDP_TDiff::ProcessEntry(), and TemplateCreator::StartTemplate().

{return fDebug;};

std::string BaseModule::GetAlias

(

)

const [inline, inherited]

Returns a string for the alias of this module. May be empty if no alias was given in the modules file.

Definition at line 63 of file BaseModule.h.

References BaseModule::fAlias.

Referenced by BaseModule::BaseModule().

{return fAlias;};

TDirectory * BaseModule::GetDirectory

(

const std::string &

name = ""

)

[protected, inherited]

Definition at line 77 of file BaseModule.cpp.

References BaseModule::dir, and BaseModule::fDirectory.

                                                         {
  if(name.empty()) return fDirectory;
  TDirectory* dir=fDirectory->GetDirectory(name.c_str());
  if(dir) return dir;
  return fDirectory->mkdir(name.c_str());
}

TDirectory* BaseModule::GetDirectory

(

)

const [inline, protected, inherited]

Get the TDirectory for this module.

Definition at line 74 of file BaseModule.h.

References BaseModule::fDirectory.

Referenced by TemplateCreator::AfterLastEntry(), SavePulses::AfterLastEntry(), TemplateCreator::BeforeFirstEntry(), PlotTDPs::BeforeFirstEntry(), ExportPulse::ExportPulse(), and TemplateCreator::ProcessEntry().

{return fDirectory;}

const char* BaseModule::GetName

(

)

const [inline, inherited]

Get the name of this module as given to the constructor of the base class.

Definition at line 66 of file BaseModule.h.

References BaseModule::fName.

Referenced by SavePulses::BeforeFirstEntry(), PlotTDPs::BeforeFirstEntry(), PlotTDP_TDiff::BeforeFirstEntry(), PlotTAP_Time::ProcessEntry(), PlotTAP_Energy::ProcessEntry(), PlotTAP_Amplitude::ProcessEntry(), and LoopSequence::Run().

{return fName.c_str();};

int BaseModule::Postprocess	(	TGlobalData *	gData,
		const TSetupData *	gSetup
	)			[inherited]

Method called by the main pre-process loop.

Does some simple work, then hooks into the derived class through AfterLastEntry.

Returns:

0 on sucess and non-zero if a problem occurred

Definition at line 66 of file BaseModule.cpp.

References BaseModule::AfterLastEntry(), and BaseModule::fDirectory.

Referenced by LoopSequence::Postprocess().

                                                                       {
  // This is called by our main routine and would allow later to split into different 
  // process routines if we have more than one Tree and hence different tpyes of data input.

  if(fDirectory) fDirectory->cd();
  int ret = AfterLastEntry(gData, gSetup);
  gDirectory->cd("/");

  return ret;
}

int BaseModule::Preprocess	(	TGlobalData *	gData,
		const TSetupData *	gSetup
	)			[inherited]

Method called by the main pre-process loop.

Does some simple work, then hooks into the derived class through BeforeFirstEntry.

Returns:

0 on sucess and non-zero if a problem occurred

Definition at line 55 of file BaseModule.cpp.

References BaseModule::BeforeFirstEntry(), and BaseModule::fDirectory.

Referenced by LoopSequence::Preprocess().

                                                                      {
  // This is called by our main routine and would allow later to split into different 
  // process routines if we have more than one Tree and hence different tpyes of data input.

  if(fDirectory) fDirectory->cd();
  int ret = BeforeFirstEntry(gData, gSetup);
  gDirectory->cd("/");

  return ret;
}

int SimpleHistograms::ProcessEntry

(

TGlobalData *

gData

)

[private, virtual]

Definition at line 54 of file SimpleHistograms.cpp.

References banks, TGlobalData::fPulseIslandToChannelMap, hBankSize, hNumberOfNoCoincidences, last, and verbose.

                                                    {
  typedef StringPulseIslandMap::iterator map_iterator;

  vector<TPulseIsland*> islands[banks.size()];
  for(unsigned int b1 = 0; b1 < banks.size(); ++b1){
    string bank = banks.at(b1);
    for(map_iterator iter = gData->fPulseIslandToChannelMap.begin();
        iter != gData->fPulseIslandToChannelMap.end(); iter++){
      if(strcmp((iter->first).data(), bank.data()) == 0){
        islands[b1] = iter->second;
      }  
    }
    hBankSize->Fill(islands[b1].size(), b1+1);
  }

  PulseIslandList islands1, islands2;

  int count  = 1;

  for(unsigned int b1 = 0; b1 < banks.size()-1; ++b1){
    for(unsigned int b2 = b1+1; b2 < banks.size(); ++b2){
      char title[100];
      sprintf(title, "%s - %s", banks.at(b1).data(), banks.at(b2).data());
      hNumberOfNoCoincidences->GetYaxis()->SetBinLabel(count, title);

      int last = 0;
      int NumberOfNoCoincidences = 0;
      
      if(verbose){
        for(unsigned int i = 0; i < islands[b1].size(); ++i){
          std::cout << "i = " << i << ": " <<  islands[b1][i]->GetTimeStamp() * islands[b1][i]->GetClockTickInNs() << "\n";
        }
        for(unsigned int i = 0; i < islands[b2].size(); ++i){
          std::cout << "i = " << i << ": " <<  islands[b2][i]->GetTimeStamp() * islands[b2][i]->GetClockTickInNs() << "\n";
        }
      }
      
      for(unsigned int i = 0; i < islands[b1].size(); ++i){
        double t1 = islands[b1][i]->GetTimeStamp(); // * islands1[i]->GetClockTickInNs();
        bool foundCoincidence = false;
        for(unsigned int j = last; j < islands[b2].size(); ++j){
          double t2 = islands[b2][j]->GetTimeStamp(); // * islands2[j]->GetClockTickInNs();
          double tdiff = t2 - t1;
          
          if(tdiff < -200.){
            last = j;
            continue;
          }
          if(tdiff > 200.){
            break;
          }
          // Found a coincidence
          foundCoincidence = true;      
        }
        // Check if we found a coincidence
        if(!foundCoincidence){
          ++NumberOfNoCoincidences;
        }
      }
      if(verbose) 
        printf("Found %d tree entries with no coincidence for banks %s and %s.\n", NumberOfNoCoincidences, 
               banks.at(b1).data(), banks.at(b2).data());
      hNumberOfNoCoincidences->Fill(NumberOfNoCoincidences, count++);
    }
  }

  return 0;
}

int BaseModule::ProcessGenericEntry	(	TGlobalData *	gData,
		const TSetupData *	gSetup
	)			[inherited]

Method called by the main event loop for each entry in the input root tree. Does some simple work, then hooks into the derived class through ProcessEntry.

Returns:

0 on sucess and non-zero if a problem occurred

Definition at line 44 of file BaseModule.cpp.

References BaseModule::fDirectory, and BaseModule::ProcessEntry().

Referenced by LoopSequence::Process().

                                                                               {
  // This is called by our main routine and would allow later to split into different 
  // process routines if we have more than one Tree and hence different tpyes of data input.

  if(fDirectory) fDirectory->cd();
  int ret = ProcessEntry(gData, gSetup);
  gDirectory->cd("/");

  return ret;
}

void BaseModule::SetAlias

(

const std::string &

alias

)

[inline, inherited]

Sets the alias for this module, which should be provided in the modules file

Definition at line 59 of file BaseModule.h.

References BaseModule::fAlias.

Referenced by BaseModule::BaseModule().

{fAlias=alias;};

---

Member Data Documentation

TDirectory* BaseModule::dir [protected, inherited]

Many modules use 'dir' still which was the old protected pointer to the modules directory. To prevent things being broken so soon, we keep this pointer available, but be warned that it will be removed shortly...

Definition at line 98 of file BaseModule.h.

Referenced by BaseModule::BaseModule(), FastSlowCompare::FastSlowCompare(), GeSpectrum::GeSpectrum(), BaseModule::GetDirectory(), MakeMuonEvents::MakeMuonEvents(), and PlotAmpVsTDiff::PlotAmpVsTDiff().

---

The documentation for this class was generated from the following files:

• rootana/src/examples/SimpleHistograms.h
• rootana/src/examples/SimpleHistograms.cpp