1 /**************************************************************************
2 * Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved. *
4 * Author: The ALICE Off-line Project. *
5 * Contributors are mentioned in the code where appropriate. *
7 * Permission to use, copy, modify and distribute this software and its *
8 * documentation strictly for non-commercial purposes is hereby granted *
9 * without fee, provided that the above copyright notice appears in all *
10 * copies and that both the copyright notice and this permission notice *
11 * appear in the supporting documentation. The authors make no claims *
12 * about the suitability of this software for any purpose. It is *
13 * provided "as is" without express or implied warranty. *
14 **************************************************************************/
16 /* $Id: AliTRDqaBuildReference.cxx 26344 2008-06-03 10:28:50Z cblume $ */
18 ////////////////////////////////////////////////////////////////////////////
20 // Build the reference histograms //
23 // Sylwester Radomski (radomski@physi.uni-heidelberg.de) //
25 ////////////////////////////////////////////////////////////////////////////
30 #include "AliTRDqaBuildReference.h"
32 //////////////////////////////////////////////////////////////////////////////////////
34 void AliTRDqaBuildReference::BuildRefHistos(TFile *file) const
37 // Build the reference histograms
42 if (file->IsZombie()) return;
45 Int_t cd = file->cd("TRD/RecPoints");
50 TH1D *r = new TH1D("qaTRD_recPoints_ampMPV_ref", "", 150, 0, 150);
51 for(Int_t i=0; i<r->GetNbinsX(); i++) {
52 Double_t x = r->GetBinCenter(i+1);
53 if (x < 10 || x > 70) r->Fill(x, 1);
54 else if (x < 25 || x > 55) r->Fill(x, 0.5);
61 r = new TH1D("qaTRD_recPoints_nCls_ref", "", 500, -0.5, 499.5);
62 for(Int_t i=0; i<r->GetNbinsX(); i++) {
63 Double_t x = r->GetBinCenter(i+1);
64 if ( (i+4)%22 > 9 || i < 10) r->Fill(x, 0.5);
65 if (x > 350) r->Fill(x, 1);
76 //////////////////////////////////////////////////////////////////////////////////////