AlphanumericHelper.cxx

Go to the documentation of this file.

#include "AlphanumericHelper.hh"
 
void AlphanumericHelper::CreateHistograms(unsigned int slot)
{
   // get the channel map and loop over it, building a map of array number to channel address
   // the NOLINT is there to keep clang-tidy from wrongly complaining that these variables need to be initialized
   std::unordered_map<unsigned int, TChannel*>* channelMap = TChannel::GetChannelMap();   // NOLINT(cppcoreguidelines-init-variables)
   std::cout << "Got " << channelMap->size() << " channels" << std::endl;
   std::map<int, unsigned int> address;   // NOLINT(cppcoreguidelines-init-variables)
 
   for(auto channel : *channelMap) {
      // use only GRIFFIN channels
      if(channel.second->GetClassType() == TGriffin::Class()) {
         // use only the "A" channel, not the "B" channel
         if(channel.second->GetMnemonic()->OutputSensor() == TMnemonic::EMnemonic::kA) {
            address[4 * (channel.second->GetDetectorNumber() - 1) + channel.second->GetCrystalNumber() + 1] = channel.second->GetAddress();
         }
      }
   }
 
   std::cout << "x-axis: " << address.rbegin()->first - address.begin()->first + 1 << " bins from " << address.begin()->first << " to " << address.rbegin()->first + 1 << std::endl;
 
   // Define Histograms
   fH2[slot]["EnergyVsChannel"] = new TH2D("EnergyVsChannel", "#gamma singles energy vs. channel", address.rbegin()->first - address.begin()->first + 1, address.begin()->first, address.rbegin()->first + 1, 4000, 0, 2000);
   fH2[slot]["ChargeVsChannel"] = new TH2D("ChargeVsChannel", "#gamma singles charge vs. channel", address.rbegin()->first - address.begin()->first + 1, address.begin()->first, address.rbegin()->first + 1, 4000, 0, 4000);
 
   for(int bin = 1; bin <= fH2[slot]["EnergyVsChannel"]->GetNbinsX(); ++bin) {
      if(address.find(bin) == address.end()) {
         std::cout << "Couldn't find array number " << bin << " in address map" << std::endl;
         continue;
      }
      fH2[slot]["EnergyVsChannel"]->GetXaxis()->SetBinLabel(bin, Form("0x%x", address.at(bin)));
      fH2[slot]["ChargeVsChannel"]->GetXaxis()->SetBinLabel(bin, Form("0x%x", address.at(bin)));
   }
}
 
void AlphanumericHelper::Exec(unsigned int slot, TGriffin& grif)
{
   // Loop over all Griffin Hits
   for(auto i = 0; i < grif.GetMultiplicity(); ++i) {
      auto grif1 = grif.GetGriffinHit(i);
      fH2[slot].at("EnergyVsChannel")->Fill(grif1->GetArrayNumber(), grif1->GetEnergy());
      fH2[slot].at("ChargeVsChannel")->Fill(grif1->GetArrayNumber(), grif1->GetCharge());
   }
}