SimPmtM64CrosstalkTable.cxx

Go to the documentation of this file.

#include "SimPmtM64CrosstalkTable.h"
#include "MessageService/MsgService.h"
#include "Conventions/Munits.h"

CVSID("$Id: SimPmtM64CrosstalkTable.cxx,v 1.9 2008/05/06 12:49:25 tagg Exp $");

ClassImp(SimPmtM64CrosstalkTable)

SimPmtM64CrosstalkTable::SimPmtM64CrosstalkTable( const VldContext& vc )
    : fContext(vc),
      fCache(0),
      fResPtr(vc)
{
  // Zero arrays.
  for(int i=0;i<65;i++) {
    fOptTotProb[i]=0;
    for(int j=0;j<65;j++) {
      fOptDist[i][j] = 0;
      fElecFrac[i][j] = 0;
      fElecK[i][j] = 0;
    }
  }

  RebuildCache();
  MSG("DetSim",Msg::kDebug) << "Creating SimPmtM64CrosstalkTable.\n";
}

void SimPmtM64CrosstalkTable::Reset(const VldContext& vc )
{
  if(fContext==vc) return;
  fResPtr.NewQuery(vc);

  VldRange newrange = fResPtr.GetValidityRec()->GetVldRange();

  if( (newrange.GetTimeStart()!=fLastRange.GetTimeStart()) ||
      (newrange.GetTimeEnd()!=fLastRange.GetTimeEnd()) ||
      ( newrange.GetDetectorMask()!=fLastRange.GetDetectorMask()) ||
      ( newrange.GetSimMask()!=fLastRange.GetSimMask()) ) {
    
    // The vld range has changed, indicating the table has been reloaded.
    fLastRange = newrange;
    RebuildCache();
  }
    
  fContext=vc;
}

const SimPmtM64Crosstalk* 
SimPmtM64CrosstalkTable::GetRow( Int_t injPix, Int_t xPix )
{
  Int_t index = injPix*100 + xPix;
  //  return fResPtr.GetRowByIndex(index);

  if(!fCache) {
    MSG("DetSim",Msg::kError) << "Crosstalk cache not yet built!" << endl;
    return 0;
  }

  return fCache[index];
}

void SimPmtM64CrosstalkTable::RebuildCache()
{
  if(fResPtr.GetNumRows()==0) return; // No data to cache.

  const int kCacheSize = 6500;
  if(fCache) delete [] fCache;
  fCache = new const SimPmtM64Crosstalk*[kCacheSize];
  for(int i=0; i<kCacheSize; i++) {
    fCache[i] = fResPtr.GetRowByIndex(i);
  }

  for(int inj=1;inj<=64;inj++) {
    fOptTotProb[inj] = 0;

    for(int xpix=1; xpix<=64; xpix++) {
      fOptDist[inj][xpix] = 0;
      fElecFrac[inj][xpix] = 0;
      fElecK[inj][xpix] = 0;    
      
      const SimPmtM64Crosstalk* row = GetRow(inj,xpix);
      if((row)&&(inj!=xpix)) {
        float optprob = row->GetOptProb();
        if(optprob < 0) optprob = 0; // No negative optical crosstalk.
        fOptTotProb[inj] += optprob;

        fElecFrac[inj][xpix] = row->GetElecFrac();
        fElecK[inj][xpix] = row->GetElecK()*Munits::fC; //convert from fC.      
      }
      fOptDist[inj][xpix] = fOptTotProb[inj]; // Cumulative probability.        
    }

    // Normalize the opt distribution.
    // .. protect against floating point error. When the FD is simualted,
    // the M64 model gets initialized too, even though there's no Xtalk table.
    float norm = 1.0;
    if(fOptTotProb[inj]<=0) 
      MSG("DetSim",Msg::kError) << "Got zero total opt xtalk on pixel " << inj << ".. that shouldn't happen."<< endl;
    else 
      norm = 1.0/fOptTotProb[inj];
    
    for(int xpix=1; xpix<=64; xpix++) {
      fOptDist[inj][xpix] *= norm;
    }
  }
}

---