Removing useless flag.
[u/mrichter/AliRoot.git] / TRD / AliTRDqaBuildReference.cxx
CommitLineData
e508dc19 1/**************************************************************************
2 * Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved. *
3 * *
4 * Author: The ALICE Off-line Project. *
5 * Contributors are mentioned in the code where appropriate. *
6 * *
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 **************************************************************************/
15
16/* $Id: AliTRDqaBuildReference.cxx 26344 2008-06-03 10:28:50Z cblume $ */
17
18////////////////////////////////////////////////////////////////////////////
19// //
20// Build the reference histograms //
21// //
22// Author: //
23// Sylwester Radomski (radomski@physi.uni-heidelberg.de) //
24// //
25////////////////////////////////////////////////////////////////////////////
26
27#include "TH1D.h"
28#include "TFile.h"
29
30#include "AliTRDqaBuildReference.h"
31
32//////////////////////////////////////////////////////////////////////////////////////
33
09a939c1 34//_____________________________________________________________________________
35AliTRDqaBuildReference &AliTRDqaBuildReference::operator=(const AliTRDqaBuildReference &qadm)
36{
37 //
38 // Assignment operator
39 //
40
41 if (this != &qadm) {
42 ((AliTRDqaBuildReference &) qadm).Copy(*this);
43 }
44
45 return *this;
46
47}
48
49//_____________________________________________________________________________
50void AliTRDqaBuildReference::Copy(TObject &/*qadm*/) const
51{
52 //
53 // Copy function
54 //
55
56}
57
e508dc19 58void AliTRDqaBuildReference::BuildRefHistos(TFile *file) const
59{
60 //
61 // Build the reference histograms
62 //
63
64 // check
65 if (!file) return;
66 if (file->IsZombie()) return;
67
68 // recpoints
69 Int_t cd = file->cd("TRD/RecPoints");
70 if (cd) {
71
72 // MPV distribution
73
74 TH1D *r = new TH1D("qaTRD_recPoints_ampMPV_ref", "", 150, 0, 150);
75 for(Int_t i=0; i<r->GetNbinsX(); i++) {
76 Double_t x = r->GetBinCenter(i+1);
77 if (x < 10 || x > 70) r->Fill(x, 1);
78 else if (x < 25 || x > 55) r->Fill(x, 0.5);
79 }
80
81 r->Fill(-1e3, 1);
82 r->Fill(1e3, 1);
83
84 // number of clusters
85 r = new TH1D("qaTRD_recPoints_nCls_ref", "", 500, -0.5, 499.5);
86 for(Int_t i=0; i<r->GetNbinsX(); i++) {
87 Double_t x = r->GetBinCenter(i+1);
88 if ( (i+4)%22 > 9 || i < 10) r->Fill(x, 0.5);
89 if (x > 350) r->Fill(x, 1);
90 }
91
92 file->Write();
93 }
94
95 // esds
96
97}
98
99
100//////////////////////////////////////////////////////////////////////////////////////