#if !defined(__CINT__) || defined(__MAKECINT__)
#include <TRD/AliTRDgeometry.h>
#include <TRD/AliTRDcluster.h>
#include <TRD/AliTRDseedV1.h>
#include <TRD/AliTRDtrackV1.h>
#include <TProfile.h>
#endif

TH1* PH(const AliTRDtrackV1* track)
{
  if (!track){
    Error("PH()", "Missing track.");
    return NULL;
  }
  
  AliTRDcluster* cls(NULL);
  AliTRDseedV1 *tracklet(NULL);

  const Int_t ntb(27);
  TProfile* ph = new TProfile("PH", "Average PH", ntb, -.5, ntb-.5);
  ph->GetXaxis()->SetTitle("drift time [1/100ns]");
  ph->GetYaxis()->SetTitle("<PH> [a.u.]");

  for (Int_t ily = 0; ily < AliTRDgeometry::kNlayer; ily++) {
    if(!(tracklet = track->GetTracklet(ily))) continue;
    if(!tracklet->IsOK()) continue;
    
    for (Int_t icl = 0; icl < AliTRDseedV1::kNclusters; icl++) {
      if(!(cls = tracklet->GetClusters(icl))) continue;
      ph->Fill(cls->GetLocalTimeBin(), cls->GetQ());
    }
  }

  return ph;
}