In This Package:

FastTrigReadSimTool.cc

Go to the documentation of this file.

#include "FastTrigReadSimTool.h"
#include "Context/TimeStamp.h"
#include "Conventions/Trigger.h"
#include "MuonProphet/MpTriggerStat.h"
#include "MuonProphet/MpMuonFate.h"
#include "Event/IHeader.h"
#include "Event/GenHeader.h"

using namespace DayaBay;

FastTrigReadSimTool::FastTrigReadSimTool(const std::string& type,
                                         const std::string& name, 
                                         const IInterface* parent)
    : GaudiTool(type,name,parent)
{
  declareInterface< IFastTrigReadSimTool >(this) ;
}

FastTrigReadSimTool::~FastTrigReadSimTool()
{}

StatusCode FastTrigReadSimTool::initialize()
{
  return StatusCode::SUCCESS;
}

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


StatusCode FastTrigReadSimTool::mutate(SimReadoutHeader* sroHeader,
                                       SimTrigHeader* stHeader,
                                       const ElecHeader& elecHeader)
{
    SimTrigCommandHeader* trigCommandHdr = new SimTrigCommandHeader;
    SimTrigCommandCollection* tcc = new SimTrigCommandCollection;
    trigCommandHdr->addCollection( tcc );

    std::vector<SimReadout*> readouts;
    readouts.clear();

    const IHeader* iSheader = (elecHeader.inputHeaders())[0];
    const IHeader* iGheader = (iSheader->inputHeaders()) [0];
    const GenHeader* genHeader = dynamic_cast<const GenHeader*>(iGheader);
    const HepMC::GenEvent* genEvent = genHeader->event();

    HepMC::GenEvent::particle_const_iterator pci, pci_end = genEvent->particles_end();

    for( pci = genEvent->particles_begin(); pci != pci_end; ++pci )  {
      if( (*pci)->pdg_id() == 13 || (*pci)->pdg_id() == -13 )  {   // muon

        MpMuonFate fate( (*pci)->status() );    
        if( fate.getRpcTriggerStat().getCode() == MpTriggerStat::kTriggered ||
            fate.getPoolTriggerStat().getCode() == MpTriggerStat::kTriggered )  {

          /* According to current RootIO for readout, 
             for unknown detector, a Readout won't be converted to a data file.
          */
          Site::Site_t siteId = elecHeader.context().GetSite();
          DetectorId::DetectorId_t detId = DetectorId::kUnknown;
          Detector det(siteId,detId);

          Trigger::TriggerType_t type = Trigger::kUnknown;

          SimTrigCommand *tc = new SimTrigCommand(det,type,0);
          std::vector<SimTrigCommand *> vtc;
          vtc.push_back(tc);
          SimTrigCommandCollection* tcc = new SimTrigCommandCollection( trigCommandHdr, det, vtc );
          trigCommandHdr -> addCollection( tcc );

          SimReadout* pSimRO = new SimReadout;
          readouts.push_back( pSimRO );
          
          Readout* pRO = new Readout;
          pRO -> setTriggerTime( elecHeader.earliest() );
          pRO -> setDetector( det );
          pSimRO -> setReadout( pRO );
          pSimRO -> setHeader( sroHeader );
        }
      }
    }

    stHeader -> setCommandHeader( trigCommandHdr );
    trigCommandHdr -> setHeader( stHeader );
  
    sroHeader->setReadouts( readouts );

    return StatusCode::SUCCESS;
}

---