Following Gene's suggestions,  I divided the sample into 4 parts  based on their dcaGlobal variable.  

1) for BHT1 matched tracks. 

2) For "minimum-biased" sample, ! BHT0,1,2->didFire()

PDF:  fit e/p for p slice by slice