[u/mrichter/AliRoot.git] / TRD / AliTRDqaBuildReference.cxx

/**************************************************************************
 * Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved. *
 *                                                                        *
 * Author: The ALICE Off-line Project.                                    *
 * Contributors are mentioned in the code where appropriate.              *
 *                                                                        *
 * Permission to use, copy, modify and distribute this software and its   *
 * documentation strictly for non-commercial purposes is hereby granted   *
 * without fee, provided that the above copyright notice appears in all   *
 * copies and that both the copyright notice and this permission notice   *
 * appear in the supporting documentation. The authors make no claims     *
 * about the suitability of this software for any purpose. It is          *
 * provided "as is" without express or implied warranty.                  *
 **************************************************************************/

/* $Id: AliTRDqaBuildReference.cxx 26344 2008-06-03 10:28:50Z cblume $ */

////////////////////////////////////////////////////////////////////////////
//                                                                        //
//  Build the reference histograms                                        //
//                                                                        //
//  Author:                                                               //
//    Sylwester Radomski (radomski@physi.uni-heidelberg.de)               //
//                                                                        //
////////////////////////////////////////////////////////////////////////////

#include "TH1D.h"
#include "TFile.h"

#include "AliTRDqaBuildReference.h"

//////////////////////////////////////////////////////////////////////////////////////

void  AliTRDqaBuildReference::BuildRefHistos(TFile *file) const
{
  //
  // Build the reference histograms
  //

  // check
  if (!file) return;
  if (file->IsZombie()) return;
  
  // recpoints
  Int_t cd = file->cd("TRD/RecPoints");
  if (cd) {
    
    // MPV distribution

    TH1D *r = new TH1D("qaTRD_recPoints_ampMPV_ref", "", 150, 0, 150);
    for(Int_t i=0; i<r->GetNbinsX(); i++) {
      Double_t x = r->GetBinCenter(i+1);
      if (x < 10 || x > 70) r->Fill(x, 1);
      else if (x < 25 || x > 55) r->Fill(x, 0.5);
    }
    
    r->Fill(-1e3, 1);
    r->Fill(1e3, 1);
    
    // number of clusters
    r = new TH1D("qaTRD_recPoints_nCls_ref", "", 500, -0.5, 499.5);
    for(Int_t i=0; i<r->GetNbinsX(); i++) {
      Double_t x = r->GetBinCenter(i+1);
      if ( (i+4)%22 > 9 || i < 10) r->Fill(x, 0.5);
      if (x > 350) r->Fill(x, 1);
    }

    file->Write();
  }

  // esds

}


//////////////////////////////////////////////////////////////////////////////////////
Commit	Line	Data
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
	34	void AliTRDqaBuildReference::BuildRefHistos(TFile *file) const
	35	{
	36	//
	37	// Build the reference histograms
	38	//
	39
	40	// check
	41	if (!file) return;
	42	if (file->IsZombie()) return;
	43
	44	// recpoints
	45	Int_t cd = file->cd("TRD/RecPoints");
	46	if (cd) {
	47
	48	// MPV distribution
	49
	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);
	55	}
	56
	57	r->Fill(-1e3, 1);
	58	r->Fill(1e3, 1);
	59
	60	// number of clusters
	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);
66	}
67
68	file->Write();
69	}
70
71	// esds
72
73	}
74
75
76	//////////////////////////////////////////////////////////////////////////////////////