+++ /dev/null
-#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;
-}
-