Commented unnecessary include AliLog
[u/mrichter/AliRoot.git] / TRD / AliTRDqaBuildReference.cxx
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
34 //_____________________________________________________________________________
35 AliTRDqaBuildReference &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 //_____________________________________________________________________________
50 void AliTRDqaBuildReference::Copy(TObject &/*qadm*/) const
51 {
52   //
53   // Copy function
54   //
55
56 }
57
58 void  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 //////////////////////////////////////////////////////////////////////////////////////