In This Package:

EsIdealPulseTool.cc

Go to the documentation of this file.

#include "EsIdealPulseTool.h"

#include "Event/SimHit.h"
#include "Event/ElecPulse.h"
#include "Event/ElecPmtPulse.h"
#include "Event/ElecRpcPulse.h"
#include "Event/SimTrackReference.h"
#include "Event/SimTrack.h"

#include "Conventions/PulseType.h"

EsIdealPulseTool::EsIdealPulseTool(const std::string& type,
                                   const std::string& name, 
                                   const IInterface* parent)
    : GaudiTool(type,name,parent)
    , m_cableSvc(0)
{
    declareInterface< IEsPulseTool >(this) ;   
    declareProperty("CableSvcName",m_cableSvcName="StaticCableSvc",
                    "Name of service to map between detector, hardware, and electronic IDs");
}

EsIdealPulseTool::~EsIdealPulseTool(){}

StatusCode EsIdealPulseTool::initialize()
{
  // Get Cable Service
  m_cableSvc = svc<ICableSvc>(m_cableSvcName,true);

  return StatusCode::SUCCESS;
}

StatusCode EsIdealPulseTool::finalize()
{
  return StatusCode::SUCCESS;
}

StatusCode EsIdealPulseTool::generatePulses(DayaBay::SimHitCollection* hits,
                                           DayaBay::ElecPulseCollection* pulses)
{    
  debug() << "running generatePulses()" << endreq; 
   
  const DayaBay::SimHitCollection::hit_container& htcon = hits->collection();
  DayaBay::SimHitCollection::hit_container::const_iterator hcIter, done = htcon.end();
  DayaBay::ElecPulseCollection::PulseContainer pulsecon;
  pulsecon.clear();

  // Context, ServiceMode for this data
  Context context(pulses->detector().site(), SimFlag::kMC, 
                  pulses->header()->header()->timeStamp(),
                  pulses->detector().detectorId());
  int task = 0;
  ServiceMode svcMode(context, task);

  for (hcIter=htcon.begin(); hcIter != done; ++hcIter) {
    DayaBay::DetectorSensor sensDetId((*hcIter)->sensDetId());
    if (sensDetId.isAD()) {
        DayaBay::AdPmtSensor pmtid((*hcIter)->sensDetId());
        if (pmtid.bogus()) {
            warning() << "generatePulses() got bogus PMT ID: " << pmtid << endreq;
        }
    }
    const DayaBay::ElecChannelId& elecChannelId = m_cableSvc->elecChannelId( 
                                                                     sensDetId,
                                                                     svcMode);
    DayaBay::ElecPulse* simpulse = 0;
    if(sensDetId.detectorId() != DetectorId::kRPC){
      simpulse = new DayaBay::ElecPmtPulse();
    }else{
      simpulse = new DayaBay::ElecRpcPulse();
    }

    // Primary vertex time (absolute, in seconds)
    TimeStamp vertexTime = (*hcIter)->hc()->header()->header()->timeStamp();
    // Need to subtract ElecHeader time (absolute, in seconds) to get pulse time wrt ElecHeader
    TimeStamp headerTime = pulses->header()->header()->timeStamp();
    // Set pulse time (nanoseconds)
    simpulse->setTime((vertexTime-headerTime)*1.e9 + (*hcIter)->hitTime());

    simpulse->setAmplitude((*hcIter)->weight());
    simpulse->setAncestor((*hcIter));
    simpulse->setType(PulseType::kPmtHit);
    simpulse->setChannelId(elecChannelId);

    pulsecon.push_back(simpulse);
  }

  pulses->setPulses(pulsecon);

  return StatusCode::SUCCESS;
}

---