raster and PSTH are inconsistent

[agbondy]

When plotting PSTHs/rasters using a trial file, the timing of the raster and PSTH does not seem to match up. This can be seen in the image below, where the raster shows a time-locked response (to an auditory pulse) but no corresponding bump in the PSTH. Digging into the code, it seems that the PSTH is plotting a narrower time slice of the data around the events than the raster is, and so the broad bump in the PSTH actually corresponds to the time-locked response seen in the raster. I spent a half hour in the code and I couldn't identify the source of the problem.

Output of ver

---

MATLAB Version: 9.4.0.813654 (R2018a)
MATLAB License Number: 296492
Operating System: Microsoft Windows 10 Pro for Workstations Version 10.0 (Build 17763)
Java Version: Java 1.8.0_144-b01 with Oracle Corporation Java HotSpot(TM) 64-Bit Server VM mixed mode

MATLAB Version 9.4 (R2018a)
Simulink Version 9.1 (R2018a)
Control System Toolbox Version 10.4 (R2018a)
Curve Fitting Toolbox Version 3.5.7 (R2018a)
DSP System Toolbox Version 9.6 (R2018a)
Data Acquisition Toolbox Version 3.13 (R2018a)
GUI Layout Toolbox Version 2.3.4 (R2018b)
Image Processing Toolbox Version 10.2 (R2018a)
Instrument Control Toolbox Version 3.13 (R2018a)
Optimization Toolbox Version 8.1 (R2018a)
Parallel Computing Toolbox Version 6.12 (R2018a)
Predictive Maintenance Toolbox Version 1.0 (R2018a)
Signal Processing Toolbox Version 8.0 (R2018a)
Simulink Control Design Version 5.1 (R2018a)
Statistics and Machine Learning Toolbox Version 11.3 (R2018a)
Symbolic Math Toolbox Version 8.1 (R2018a)
System Identification Toolbox Version 9.8 (R2018a)

Output of gpuDevice

                  Name: 'Quadro P5000'
                 Index: 1
     ComputeCapability: '6.1'
        SupportsDouble: 1
         DriverVersion: 10.2000
        ToolkitVersion: 9
    MaxThreadsPerBlock: 1024
      MaxShmemPerBlock: 49152
    MaxThreadBlockSize: [1024 1024 64]
           MaxGridSize: [2.1475e+09 65535 65535]
             SIMDWidth: 32
           TotalMemory: 1.7180e+10
       AvailableMemory: 1.4096e+10
   MultiprocessorCount: 20
          ClockRateKHz: 1733500
           ComputeMode: 'Default'
  GPUOverlapsTransfers: 1
KernelExecutionTimeout: 1
      CanMapHostMemory: 1
       DeviceSupported: 1
        DeviceSelected: 1