PhotonNdAfterpulseModel.cxx

Go to the documentation of this file.

#include "PhotonNdAfterpulseModel.h"
#include "PhotonTransportMaker.h"
#include "ScintPhoton.h"
#include "MessageService/MsgService.h"
#include "Plex/PlexHandle.h"
#include "Digitization/DigiSignal.h"

CVSID( " $Id: PhotonNdAfterpulseModel.cxx,v 1.3 2006/12/01 20:12:11 rhatcher Exp $" );

ClassImp(PhotonNdAfterpulseModel)

PhotonTransportMakerProxy<PhotonNdAfterpulseModel>
        gPhotonNdAfterpulseProxy("photonNdAfterpulseModel"); 


// Constructor
PhotonNdAfterpulseModel::PhotonNdAfterpulseModel() 
{
  Configure(PhotonConfiguration()); // Set up default values.
}

// Reset
void 
PhotonNdAfterpulseModel::Reset( const VldContext& /* cx */)
{
  // Called once per event. Use this to get or update
  // any database tables used by the model.
}


// Configure the model.
void 
PhotonNdAfterpulseModel::Configure( const Registry& r )
{
  // The other classses in the model need configuring:
  ScintPhoton::Configure(r);

  MSG("Photon",Msg::kDebug) << "Configuring NdAfterpulseModel." << std::endl;
  
  // Load any configuration options from the provided registry.
  // Called only when configuration changes.
  // Put default values into the file PhotonConfiguration.cxx

  fAfterpulseFraction = r.GetDouble("AfterpulseFraction");
  fAfterpulseStart    = r.GetDouble("AfterpulseStart");
  fAfterpulseEnd      = r.GetDouble("AfterpulseEnd");
}


void PhotonNdAfterpulseModel::MakeAfterpulses( const std::vector<DigiPE*>& inPeList,
                                               std::vector<DigiPE*>& outPeList )
{

  // Simple model: generate a pe at random on the same tube,
  // randomly on pixels 1-8 or 57-64, between 1 and 4us after the
  // intial pulse, 4% of the time.

  for(UInt_t i=0;i<inPeList.size(); i++) {
    const DigiPE* inPE = inPeList[i];
    PlexPixelSpotId inPix = inPE->GetPixelSpotId();
    if(inPix.GetElecType()==ElecType::kQIE) {
      // This is an M64.
      
      // Roll to see if an afterpulse is made:
      if(fRandom->Rndm() < fAfterpulseFraction) {
        
        // Make the new one.

        // Chose the pixel:
        int pixel = fRandom->Integer(8)+1; // 1-8
        if(fRandom->Rndm()>0.5) pixel+= 56; // 57-64

        // Set the PixelSpotId:
        PlexPixelSpotId outPix = inPix;
        outPix.SetSpot(1);
        outPix.SetPixel(pixel);

        // Chose the time of the PE:
        double time = inPE->GetTime();
        time += fRandom->Uniform(fAfterpulseStart,fAfterpulseEnd);
        
        // Create it and put it into the output list.
        DigiPE* outPE = new DigiPE(time, outPix, DigiSignal::kAfterpulse );
        outPeList.push_back(outPE);
      }

    }
  }

  return;
}




void PhotonNdAfterpulseModel::Print(Option_t* /* option */) const
{
}

---