InfiniteWindowStrategy.cc

Go to the documentation of this file.

#include "InfiniteWindowStrategy.h"

namespace MdMuonEstimatorAG {

  InfiniteWindowStrategy::InfiniteWindowStrategy(unsigned int nTraceLength) :
    fTraceSize(nTraceLength)
  {
  }
  
  double
  InfiniteWindowStrategy::operator()(std::vector<std::pair<unsigned int, unsigned int>>& vModulePatternMatches, mevt::ModuleRecData& mRecData)
    const
  {
    //std::cout << "In infinite window strategy" << std::endl;
    //std::cout << "Size " << vModulePatternMatches.size() << std::endl;

    double mu = 0;
    double nmu = 0;
    bool isSat = 0;

    if (vModulePatternMatches.size()>0) {
      unsigned int k = 0; // Count the number of channels on (with signal)
      //std::cout << "Channels with signal ";
      //for (std::vector<std::pair<unsigned int, unsigned int>>::const_iterator it = vModulePatternMatches.begin(); it != vModulePatternMatches.end(); ++it) {
        //std::cout << it->first << " ";
      //}
      //std::cout << std::endl;
      for (unsigned int ich = 1 ; ich <= 64 ; ++ich)
      {
        for (std::vector<std::pair<unsigned int, unsigned int>>::const_iterator it = vModulePatternMatches.begin(); it != vModulePatternMatches.end(); ++it) {
          if (it->first == ich) {
            ++k;
            //std::cout << ich << " is on" << std::endl;
            break;
          }
        }
      }
      unsigned int nActiveCh = mRecData.GetSegmentation(); // = number of active channels
      //std::cout << "Total channels on " << k << " out of " << nActiveCh << " active channels " << std::endl;
      if (k >= nActiveCh){  
        isSat = 1;
        //std::cout << "Module is saturated" << std::endl;
      }
      mu = -(double) nActiveCh * log(1. - (double) k/(double) nActiveCh);
      nmu = log(1. - (double) k/(double) nActiveCh) / log(1. - 1./(double) nActiveCh);
    } 
    //std::cout << "Mu " << mu << std::endl;
    //std::cout << "Nmu " << nmu << std::endl;
    //std::cout << "IsSat " << isSat << std::endl;
    
    mRecData.SetMeanMuons(mu);
    mRecData.SetNumberOfEstimatedMuons(nmu);
    mRecData.SetSaturated(isSat);
    
    return nmu;
  } 
}