mEmcMaskDeadTowers.C

Go to the documentation of this file.

#include "mEmcMaskDeadTowers.h"
#include "PHCompositeNode.h"

#include "emcBadModules.h"
#include "emcTowerContainer.h"
#include "emcTowerContent.h"
#include "emcNodeHelper.h"
#include "EmcIndexer.h"
#include "recoConsts.h"
#include "getClass.h"

#include <fstream>
#include <iostream>

using namespace std;

//_____________________________________________________________________________
mEmcMaskDeadTowers::mEmcMaskDeadTowers()
{
  name = "mEmcMaskDeadTowers";
  recoConsts *rc = recoConsts::instance();

  unsigned int err, warn;
  unsigned int towerid;

  fErrorRaw.resize(NCHANNELS, 0);
  fWarnRaw.resize(NCHANNELS, 0);

  fErrorRaw.resize(NCHANNELS, 0);
  fWarnRaw.resize(NCHANNELS, 0);

  fErrorMap.resize(fErrorRaw.size() , 0 );

  ifstream fin(rc->get_CharFlag("EMCDEADRECALDATASOURCE"));
  while (fin >> towerid >> err >> warn)
    {
      cout << towerid << " " << err << endl;
      fErrorRaw[towerid] = err;
      fWarnRaw[towerid] = warn;
    }
  fin.close();
}

//_____________________________________________________________________________
mEmcMaskDeadTowers::~mEmcMaskDeadTowers()
{
}

//_____________________________________________________________________________
PHBoolean
mEmcMaskDeadTowers::event(PHCompositeNode* topNode)
{
  emcTowerContainer* towers =
    findNode::getClass<emcTowerContainer>(topNode, "emcTowerContainer");
  
  if (!towers)
    return 0;
  
  for ( size_t i = 0; i < towers->size(); ++i )
    {
      emcTowerContent* twr = towers->getTower(i);
      if (twr == 0)
        {
          cout << PHWHERE << " SEVERE ERROR: Got NULL pointer for tower, contact the emc group with this"
               << endl;
          exit(1);
        }
      
      int towerid = twr->TowerID();
      unsigned int cdead = ErrorFast(towerid);
      if (cdead & 0x400)
        {
          //          twr->SetCalibrated(0.0,0.0);
          //          cout << towerid << "  " << cdead << endl;
          twr->Zero();
          twr->SetNeighbours(cdead,0);
        }
    }
}

//_____________________________________________________________________________
unsigned int
mEmcMaskDeadTowers::ErrorFast(int towerID)
{
  return fErrorRaw[towerID];
}

//_____________________________________________________________________________
unsigned int
mEmcMaskDeadTowers::DeadmapFast(int towerID)
{
  return fErrorMap[towerID];
}

//_____________________________________________________________________________
unsigned int
mEmcMaskDeadTowers::WarningFast(int towerID)
{
  return fWarnRaw[towerID];
}