[u/mrichter/AliRoot.git] / TRD / TRDbase / AliTRDCalPadStatus.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$ */

///////////////////////////////////////////////////////////////////////////////
//                                                                           //
//  TRD calibration class for the single pad status                          //
//                                                                           //
///////////////////////////////////////////////////////////////////////////////

#include <TStyle.h>
#include <TCanvas.h>
#include <TH1F.h>
#include <TH2F.h>

#include "AliTRDCalPadStatus.h"
#include "AliTRDgeometry.h"
#include "AliTRDpadPlane.h"
#include "AliTRDCalSingleChamberStatus.h"  // test

ClassImp(AliTRDCalPadStatus)

//_____________________________________________________________________________
AliTRDCalPadStatus::AliTRDCalPadStatus()
  :TNamed()
{
  //
  // AliTRDCalPadStatus default constructor
  //

  for (Int_t idet = 0; idet < kNdet; idet++) {
    fROC[idet] = 0;
  }

}

//_____________________________________________________________________________
AliTRDCalPadStatus::AliTRDCalPadStatus(const Text_t *name, const Text_t *title)
  :TNamed(name,title)
{
  //
  // AliTRDCalPadStatus constructor
  //

  for (Int_t isec = 0; isec < kNsect; isec++) {
    for (Int_t ipla = 0; ipla < kNplan; ipla++) {
      for (Int_t icha = 0; icha < kNcham; icha++) {
        Int_t idet = AliTRDgeometry::GetDetector(ipla,icha,isec);
        fROC[idet] = new AliTRDCalSingleChamberStatus(ipla,icha,144);
      }
    }
  }

}

//_____________________________________________________________________________
AliTRDCalPadStatus::AliTRDCalPadStatus(const AliTRDCalPadStatus &c)
  :TNamed(c)
{
  //
  // AliTRDCalPadStatus copy constructor
  //

  ((AliTRDCalPadStatus &) c).Copy(*this);

}

//_____________________________________________________________________________
AliTRDCalPadStatus::~AliTRDCalPadStatus()
{
  //
  // AliTRDCalPadStatus destructor
  //

  for (Int_t idet = 0; idet < kNdet; idet++) {
    if (fROC[idet]) {
      delete fROC[idet];
      fROC[idet] = 0;
    }
  }

}

//_____________________________________________________________________________
AliTRDCalPadStatus &AliTRDCalPadStatus::operator=(const AliTRDCalPadStatus &c)
{
  //
  // Assignment operator
  //

  if (this != &c) ((AliTRDCalPadStatus &) c).Copy(*this);
  return *this;

}

//_____________________________________________________________________________
void AliTRDCalPadStatus::Copy(TObject &c) const
{
  //
  // Copy function
  //

  for (Int_t idet = 0; idet < kNdet; idet++) {
    if (fROC[idet]) {
      fROC[idet]->Copy(*((AliTRDCalPadStatus &) c).fROC[idet]);
    }
  }

  TObject::Copy(c);

}

//_____________________________________________________________________________
Bool_t AliTRDCalPadStatus::CheckStatus(Int_t d, Int_t col, Int_t row, Int_t bitMask) const
{
  //
  // Checks the pad status
  //

  AliTRDCalSingleChamberStatus *roc = GetCalROC(d);
  if (!roc) {
    return kFALSE;
  }
  else {
    return (roc->GetStatus(col, row) & bitMask) ? kTRUE : kFALSE;
  }

}

//_____________________________________________________________________________
AliTRDCalSingleChamberStatus* AliTRDCalPadStatus::GetCalROC(Int_t p, Int_t c, Int_t s) const
{ 
  //
  // Returns the readout chamber of this pad
  //

  return fROC[AliTRDgeometry::GetDetector(p,c,s)];   

}

//_____________________________________________________________________________
TH1F *AliTRDCalPadStatus::MakeHisto1D()
{
  //
  // Make 1D histo
  //

  char  name[1000];
  snprintf(name,1000,"%s Pad 1D",GetTitle());
  TH1F * his = new TH1F(name,name,6, -0.5,5.5);
  his->GetXaxis()->SetBinLabel(1,"Good");
  his->GetXaxis()->SetBinLabel(2,"Masked");
  his->GetXaxis()->SetBinLabel(3,"PadBridgedLeft");
  his->GetXaxis()->SetBinLabel(4,"PadBridgedRight");
  his->GetXaxis()->SetBinLabel(5,"ReadSecond");
  his->GetXaxis()->SetBinLabel(6,"NotConnected");

  for (Int_t idet = 0; idet < kNdet; idet++) 
    {
      if (fROC[idet])
	{
	  for (Int_t ichannel=0; ichannel<fROC[idet]->GetNchannels(); ichannel++)
	    {
	      Int_t status = (Int_t) fROC[idet]->GetStatus(ichannel);
	      if(status==2)  status= 1;
	      if(status==4)  status= 2;
	      if(status==8)  status= 3;
	      if(status==16) status= 4;
	      if(status==32) status= 5;
	      his->Fill(status);
	    }
	}
    }

  return his;

}

//_____________________________________________________________________________
TH2F *AliTRDCalPadStatus::MakeHisto2DSmPl(Int_t sm, Int_t pl)
{
  //
  // Make 2D graph
  //

  gStyle->SetPalette(1);
  AliTRDgeometry *trdGeo = new AliTRDgeometry();
  AliTRDpadPlane *padPlane0 = trdGeo->GetPadPlane(pl,0);
  Double_t row0    = padPlane0->GetRow0();
  Double_t col0    = padPlane0->GetCol0();

  char  name[1000];
  snprintf(name,1000,"%s Pad 2D sm %d pl %d",GetTitle(),sm,pl);
  TH2F * his = new TH2F( name, name, 88,-TMath::Abs(row0),TMath::Abs(row0)
                                   ,148,-TMath::Abs(col0),TMath::Abs(col0));

  // Where we begin
  Int_t offsetsmpl = 30*sm+pl;

  for (Int_t k = 0; k < kNcham; k++){
    Int_t det = offsetsmpl+k*6;
    if (fROC[det]){
      AliTRDCalSingleChamberStatus * calRoc = fROC[det];
      for (Int_t icol=0; icol<calRoc->GetNcols(); icol++){
	for (Int_t irow=0; irow<calRoc->GetNrows(); irow++){
	  Int_t binz     = 0;
	  Int_t kb       = kNcham-1-k;
	  Int_t krow     = calRoc->GetNrows()-1-irow;
	  Int_t kcol     = calRoc->GetNcols()-1-icol;
	  if(kb > 2) binz = 16*(kb-1)+12+krow+1+2*(kb+1);
	  else binz = 16*kb+krow+1+2*(kb+1); 
	  Int_t biny = kcol+1+2;
	  Float_t value = calRoc->GetStatus(icol,irow);
	  his->SetBinContent(binz,biny,value);
	}
      }
      for(Int_t icol = 1; icol < 147; icol++){
	for(Int_t l = 0; l < 2; l++){
	  Int_t binz     = 0;
	  Int_t kb       = kNcham-1-k;
	  if(kb > 2) binz = 16*(kb-1)+12+1+2*(kb+1)-(l+1);
	  else binz = 16*kb+1+2*(kb+1)-(l+1); 
	  his->SetBinContent(binz,icol,50.0);
	}
      }
    }
  }
  for(Int_t icol = 1; icol < 147; icol++){
    his->SetBinContent(88,icol,50.0);
    his->SetBinContent(87,icol,50.0);
  }
  for(Int_t irow = 1; irow < 89; irow++){
    his->SetBinContent(irow,1,50.0);
    his->SetBinContent(irow,2,50.0);
    his->SetBinContent(irow,147,50.0);
    his->SetBinContent(irow,148,50.0);
  }

  his->SetXTitle("z (cm)");
  his->SetYTitle("y (cm)");
  his->SetMaximum(50);
  his->SetMinimum(0.0);
  his->SetStats(0);

  return his;

}

//_____________________________________________________________________________
void AliTRDCalPadStatus::PlotHistos2DSm(Int_t sm, const Char_t *name)
{
  //
  // Make 2D graph
  //

  gStyle->SetPalette(1);
  TCanvas *c1 = new TCanvas(name,name,50,50,600,800);
  c1->Divide(3,2);
  c1->cd(1);
  MakeHisto2DSmPl(sm,0)->Draw("colz");
  c1->cd(2);
  MakeHisto2DSmPl(sm,1)->Draw("colz");
  c1->cd(3);
  MakeHisto2DSmPl(sm,2)->Draw("colz");
  c1->cd(4);
  MakeHisto2DSmPl(sm,3)->Draw("colz");
  c1->cd(5);
  MakeHisto2DSmPl(sm,4)->Draw("colz");
  c1->cd(6);
  MakeHisto2DSmPl(sm,5)->Draw("colz");

}
Commit	Line	Data
7754cd1f	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$ */
	17
	18	///////////////////////////////////////////////////////////////////////////////
	19	// //
2745a409	20	// TRD calibration class for the single pad status //
7754cd1f	21	// //
	22	///////////////////////////////////////////////////////////////////////////////
	23
3a0f6479	24	#include <TStyle.h>
	25	#include <TCanvas.h>
	26	#include <TH1F.h>
	27	#include <TH2F.h>
	28
7754cd1f	29	#include "AliTRDCalPadStatus.h"
2745a409	30	#include "AliTRDgeometry.h"
3a0f6479	31	#include "AliTRDpadPlane.h"
6eddce12	32	#include "AliTRDCalSingleChamberStatus.h" // test
7754cd1f	33
	34	ClassImp(AliTRDCalPadStatus)
	35
	36	//_____________________________________________________________________________
2745a409	37	AliTRDCalPadStatus::AliTRDCalPadStatus()
2745a409	38	:TNamed()
7754cd1f	39	{
	40	//
	41	// AliTRDCalPadStatus default constructor
	42	//
	43
	44	for (Int_t idet = 0; idet < kNdet; idet++) {
	45	fROC[idet] = 0;
	46	}
	47
	48	}
	49
	50	//_____________________________________________________________________________
	51	AliTRDCalPadStatus::AliTRDCalPadStatus(const Text_t name, const Text_t title)
2745a409	52	:TNamed(name,title)
7754cd1f	53	{
	54	//
	55	// AliTRDCalPadStatus constructor
	56	//
	57
	58	for (Int_t isec = 0; isec < kNsect; isec++) {
	59	for (Int_t ipla = 0; ipla < kNplan; ipla++) {
	60	for (Int_t icha = 0; icha < kNcham; icha++) {
	61	Int_t idet = AliTRDgeometry::GetDetector(ipla,icha,isec);
	62	fROC[idet] = new AliTRDCalSingleChamberStatus(ipla,icha,144);
	63	}
	64	}
	65	}
	66
	67	}
	68
7754cd1f	69	//_____________________________________________________________________________
2745a409	70	AliTRDCalPadStatus::AliTRDCalPadStatus(const AliTRDCalPadStatus &c)
2745a409	71	:TNamed(c)
7754cd1f	72	{
	73	//
	74	// AliTRDCalPadStatus copy constructor
	75	//
	76
	77	((AliTRDCalPadStatus &) c).Copy(*this);
	78
	79	}
	80
3a0f6479	81	//_____________________________________________________________________________
7754cd1f	82	AliTRDCalPadStatus::~AliTRDCalPadStatus()
	83	{
	84	//
	85	// AliTRDCalPadStatus destructor
	86	//
	87
	88	for (Int_t idet = 0; idet < kNdet; idet++) {
	89	if (fROC[idet]) {
	90	delete fROC[idet];
	91	fROC[idet] = 0;
	92	}
	93	}
	94
	95	}
	96
	97	//_____________________________________________________________________________
	98	AliTRDCalPadStatus &AliTRDCalPadStatus::operator=(const AliTRDCalPadStatus &c)
	99	{
	100	//
	101	// Assignment operator
	102	//
	103
	104	if (this != &c) ((AliTRDCalPadStatus &) c).Copy(*this);
	105	return *this;
	106
	107	}
	108
	109	//_____________________________________________________________________________
	110	void AliTRDCalPadStatus::Copy(TObject &c) const
	111	{
	112	//
	113	// Copy function
	114	//
	115
	116	for (Int_t idet = 0; idet < kNdet; idet++) {
	117	if (fROC[idet]) {
	118	fROC[idet]->Copy(*((AliTRDCalPadStatus &) c).fROC[idet]);
	119	}
	120	}
	121
	122	TObject::Copy(c);
2745a409	123
7754cd1f	124	}
7754cd1f	125
2745a409	126	//_____________________________________________________________________________
	127	Bool_t AliTRDCalPadStatus::CheckStatus(Int_t d, Int_t col, Int_t row, Int_t bitMask) const
	128	{
	129	//
	130	// Checks the pad status
	131	//
	132
	133	AliTRDCalSingleChamberStatus *roc = GetCalROC(d);
	134	if (!roc) {
	135	return kFALSE;
	136	}
	137	else {
	138	return (roc->GetStatus(col, row) & bitMask) ? kTRUE : kFALSE;
	139	}
	140
	141	}
	142
	143	//_____________________________________________________________________________
	144	AliTRDCalSingleChamberStatus* AliTRDCalPadStatus::GetCalROC(Int_t p, Int_t c, Int_t s) const
	145	{
	146	//
	147	// Returns the readout chamber of this pad
	148	//
	149
	150	return fROC[AliTRDgeometry::GetDetector(p,c,s)];
	151
	152	}
3a0f6479	153
	154	//_____________________________________________________________________________
	155	TH1F *AliTRDCalPadStatus::MakeHisto1D()
	156	{
	157	//
6eddce12	158	// Make 1D histo
3a0f6479	159	//
	160
	161	char name[1000];
02f3bfcc	162	snprintf(name,1000,"%s Pad 1D",GetTitle());
6eddce12	163	TH1F * his = new TH1F(name,name,6, -0.5,5.5);
	164	his->GetXaxis()->SetBinLabel(1,"Good");
	165	his->GetXaxis()->SetBinLabel(2,"Masked");
	166	his->GetXaxis()->SetBinLabel(3,"PadBridgedLeft");
	167	his->GetXaxis()->SetBinLabel(4,"PadBridgedRight");
	168	his->GetXaxis()->SetBinLabel(5,"ReadSecond");
	169	his->GetXaxis()->SetBinLabel(6,"NotConnected");
	170
	171	for (Int_t idet = 0; idet < kNdet; idet++)
	172	{
	173	if (fROC[idet])
	174	{
	175	for (Int_t ichannel=0; ichannel<fROC[idet]->GetNchannels(); ichannel++)
	176	{
	177	Int_t status = (Int_t) fROC[idet]->GetStatus(ichannel);
	178	if(status==2) status= 1;
	179	if(status==4) status= 2;
	180	if(status==8) status= 3;
	181	if(status==16) status= 4;
	182	if(status==32) status= 5;
	183	his->Fill(status);
	184	}
	185	}
3a0f6479	186	}
3a0f6479	187
	188	return his;
	189
	190	}
	191
	192	//_____________________________________________________________________________
	193	TH2F *AliTRDCalPadStatus::MakeHisto2DSmPl(Int_t sm, Int_t pl)
	194	{
	195	//
	196	// Make 2D graph
	197	//
	198
	199	gStyle->SetPalette(1);
	200	AliTRDgeometry *trdGeo = new AliTRDgeometry();
	201	AliTRDpadPlane *padPlane0 = trdGeo->GetPadPlane(pl,0);
	202	Double_t row0 = padPlane0->GetRow0();
	203	Double_t col0 = padPlane0->GetCol0();
	204
	205	char name[1000];
02f3bfcc	206	snprintf(name,1000,"%s Pad 2D sm %d pl %d",GetTitle(),sm,pl);
3a0f6479	207	TH2F * his = new TH2F( name, name, 88,-TMath::Abs(row0),TMath::Abs(row0)
	208	,148,-TMath::Abs(col0),TMath::Abs(col0));
	209
	210	// Where we begin
	211	Int_t offsetsmpl = 30*sm+pl;
	212
	213	for (Int_t k = 0; k < kNcham; k++){
	214	Int_t det = offsetsmpl+k*6;
	215	if (fROC[det]){
	216	AliTRDCalSingleChamberStatus * calRoc = fROC[det];
	217	for (Int_t icol=0; icol<calRoc->GetNcols(); icol++){
	218	for (Int_t irow=0; irow<calRoc->GetNrows(); irow++){
	219	Int_t binz = 0;
	220	Int_t kb = kNcham-1-k;
7dcf6933	221	Int_t krow = calRoc->GetNrows()-1-irow;
	222	Int_t kcol = calRoc->GetNcols()-1-icol;
	223	if(kb > 2) binz = 16(kb-1)+12+krow+1+2(kb+1);
	224	else binz = 16kb+krow+1+2(kb+1);
	225	Int_t biny = kcol+1+2;
3a0f6479	226	Float_t value = calRoc->GetStatus(icol,irow);
	227	his->SetBinContent(binz,biny,value);
	228	}
	229	}
	230	for(Int_t icol = 1; icol < 147; icol++){
	231	for(Int_t l = 0; l < 2; l++){
	232	Int_t binz = 0;
	233	Int_t kb = kNcham-1-k;
	234	if(kb > 2) binz = 16(kb-1)+12+1+2(kb+1)-(l+1);
	235	else binz = 16kb+1+2(kb+1)-(l+1);
7dcf6933	236	his->SetBinContent(binz,icol,50.0);
3a0f6479	237	}
	238	}
	239	}
	240	}
	241	for(Int_t icol = 1; icol < 147; icol++){
7dcf6933	242	his->SetBinContent(88,icol,50.0);
7dcf6933	243	his->SetBinContent(87,icol,50.0);
3a0f6479	244	}
3a0f6479	245	for(Int_t irow = 1; irow < 89; irow++){
7dcf6933	246	his->SetBinContent(irow,1,50.0);
	247	his->SetBinContent(irow,2,50.0);
	248	his->SetBinContent(irow,147,50.0);
	249	his->SetBinContent(irow,148,50.0);
3a0f6479	250	}
	251
	252	his->SetXTitle("z (cm)");
	253	his->SetYTitle("y (cm)");
7dcf6933	254	his->SetMaximum(50);
3a0f6479	255	his->SetMinimum(0.0);
	256	his->SetStats(0);
	257
	258	return his;
	259
	260	}
	261
	262	//_____________________________________________________________________________
	263	void AliTRDCalPadStatus::PlotHistos2DSm(Int_t sm, const Char_t *name)
	264	{
	265	//
	266	// Make 2D graph
	267	//
	268
	269	gStyle->SetPalette(1);
	270	TCanvas *c1 = new TCanvas(name,name,50,50,600,800);
	271	c1->Divide(3,2);
	272	c1->cd(1);
	273	MakeHisto2DSmPl(sm,0)->Draw("colz");
	274	c1->cd(2);
	275	MakeHisto2DSmPl(sm,1)->Draw("colz");
	276	c1->cd(3);
	277	MakeHisto2DSmPl(sm,2)->Draw("colz");
	278	c1->cd(4);
	279	MakeHisto2DSmPl(sm,3)->Draw("colz");
	280	c1->cd(5);
	281	MakeHisto2DSmPl(sm,4)->Draw("colz");
	282	c1->cd(6);
	283	MakeHisto2DSmPl(sm,5)->Draw("colz");
	284
	285	}