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 //////////////////////////////////////////////////////////////////////////////////////