[u/mrichter/AliRoot.git] / MUON / AliMUONTriggerChamberEfficiency.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$

#include "AliMUONTriggerEfficiencyCells.h"
#include "AliMpConstants.h"
#include "AliMUONConstants.h"

// Classes for display
#include "AliMUONTriggerDisplay.h"
#include "AliCDBManager.h"
#include "AliMpDDLStore.h"

#include "AliLog.h"

#include "TRandom.h"
#include "Riostream.h"
#include "TH1F.h"
#include "TObjArray.h"
#include "TGraphAsymmErrors.h"

#include "TH2F.h"
#include "TCanvas.h"
#include "TROOT.h"

#include "AliMUONTriggerChamberEfficiency.h"

//-----------------------------------------------------------------------------
/// \class AliMUONTriggerChamberEfficiency
/// A class to store and give access to the trigger chamber efficiency.
///
/// Efficiency is stored per cathode on local boards
///
/// The main method of this class is IsTriggered().
///
/// \author Diego Stocco; INFN Torino
//-----------------------------------------------------------------------------

/// \cond CLASSIMP
ClassImp(AliMUONTriggerChamberEfficiency)
/// \endcond

//__________________________________________________________________________
AliMUONTriggerChamberEfficiency::AliMUONTriggerChamberEfficiency(AliMUONTriggerEfficiencyCells* effCells)
:
TObject(),
fIsOwner(kFALSE),
fEfficiencyMap(effCells),
fEfficiencyObjects(0x0),
fDisplayList(0x0)
{
///  Default constructor.
  FillFromList();
}

//__________________________________________________________________________
AliMUONTriggerChamberEfficiency::AliMUONTriggerChamberEfficiency(const Char_t* filename, const Char_t* listname)
:
TObject(),
fIsOwner(kTRUE),
fEfficiencyMap(0x0),
fEfficiencyObjects(0x0),
fDisplayList(0x0)
{
///  Constructor using an ASCII file.
  fEfficiencyMap = new AliMUONTriggerEfficiencyCells(filename, listname);
  FillFromList();
}


//_____________________________________________________________________________
AliMUONTriggerChamberEfficiency::AliMUONTriggerChamberEfficiency(const AliMUONTriggerChamberEfficiency& other)
:
TObject(other),
fIsOwner(other.fIsOwner),
fEfficiencyMap(other.fEfficiencyMap),
fEfficiencyObjects(other.fEfficiencyObjects),
fDisplayList(other.fDisplayList)
{
/// Copy constructor
}

//_____________________________________________________________________________
AliMUONTriggerChamberEfficiency& AliMUONTriggerChamberEfficiency::operator=(const AliMUONTriggerChamberEfficiency& other)
{
  /// Asignment operator
  // check assignement to self
  if (this == &other)
    return *this;

  fIsOwner = other.fIsOwner;
  fEfficiencyMap = other.fEfficiencyMap;
  fEfficiencyObjects = other.fEfficiencyObjects;
  fDisplayList = other.fDisplayList;
    
  return *this;
}

//__________________________________________________________________________
AliMUONTriggerChamberEfficiency::~AliMUONTriggerChamberEfficiency()
{
///  Destructor.
  if ( fIsOwner )
    delete fEfficiencyMap;
  delete fEfficiencyObjects;
  delete fDisplayList;
}


//__________________________________________________________________________
Float_t AliMUONTriggerChamberEfficiency::GetCellEfficiency(Int_t detElemId, Int_t localBoard, Int_t hType) const
{
///  Get the efficiencies of the 2 cathodes at a given local board

  Int_t chamber = FindChamberIndex(detElemId);
  Int_t index = GetIndex(kHboardEff, hType, chamber);
  TGraphAsymmErrors* effGraph = ((TGraphAsymmErrors*)fEfficiencyObjects->At(index));

  // Some graphs are not available in the old implementation
  if ( ! effGraph ) return -1.;

  Double_t xpt, ypt;
  effGraph->GetPoint(localBoard-1, xpt, ypt);
  return ypt;
}


//__________________________________________________________________________
void 
AliMUONTriggerChamberEfficiency::IsTriggered(Int_t detElemId, Int_t localBoard, Bool_t &trigBend, Bool_t &trigNonBend) const
{
///  Whether or not a given local board has a chance to trig, on each cathode.

  // P(B) : probability to fire bending plane
  Float_t effBend = GetCellEfficiency(detElemId, localBoard, AliMUONTriggerEfficiencyCells::kBendingEff);

  // P(BN) : probability to fire bending and non-bending plane
  Float_t effBoth = GetCellEfficiency(detElemId, localBoard, AliMUONTriggerEfficiencyCells::kBothPlanesEff);

  trigBend =  ( gRandom->Rndm() > effBend ) ? kFALSE : kTRUE;

  // P(N) : probability to fire non-bending plane
  Float_t effNonBend = GetCellEfficiency(detElemId, localBoard, AliMUONTriggerEfficiencyCells::kNonBendingEff);

  if ( effBoth > 0 ) {
    effNonBend = ( trigBend ) ? 
      effBoth / effBend :   // P(N|B) = P(BN) / P(B)
      ( effNonBend - effBoth ) / ( 1. - effBend );  // P(N|!B) = ( P(N) - P(BN) ) / ( 1 - P(B) )
  }

  trigNonBend =  ( gRandom->Rndm() > effNonBend ) ? kFALSE : kTRUE;

  AliDebug(2,Form("Ch %i  board %i  resp (%i, %i)  prob (%.2f, %.2f)  effNB %.2f  effBoth %.2f\n", detElemId/100, localBoard, trigBend, trigNonBend, effBend, effNonBend, GetCellEfficiency(detElemId, localBoard, AliMUONTriggerEfficiencyCells::kNonBendingEff), effBoth));


}


//__________________________________________________________________________
Int_t AliMUONTriggerChamberEfficiency::FindChamberIndex(Int_t detElemId) const
{
///  From detElemId to chamber number

  // Int_t iChamber = AliMpDEManager::GetChamberId(detElemId);
  Int_t iChamber = detElemId/100 - 1;
  return iChamber-AliMpConstants::NofTrackingChambers();
}


//__________________________________________________________________________
void
AliMUONTriggerChamberEfficiency::FillFromList(Bool_t useMeanValues)
{
///  Fills internal histos from list.

  if ( fEfficiencyObjects )
    delete fEfficiencyObjects;

  const Int_t kNeffHistos = 
    2 * ( AliMUONTriggerEfficiencyCells::kNcounts - 1 ) * AliMUONConstants::NTriggerCh();

  fEfficiencyObjects = new TObjArray(kNeffHistos);
  fEfficiencyObjects->SetOwner();

  TH1F *histoNum = 0x0, *histoDen=0x0;
  TString histoName = "";
  Int_t deType[2] = {AliMUONTriggerEfficiencyCells::kHboardCount,
		     AliMUONTriggerEfficiencyCells::kHslatCount};
  Int_t deTypeEff[2] = {kHboardEff, kHslatEff};
  Int_t index = -1;

  Bool_t rebuildEfficiency = kTRUE;

  for ( Int_t ide=0; ide<2; ide++){
    Int_t currDe = deType[ide];

    if ( useMeanValues && currDe == AliMUONTriggerEfficiencyCells::kHboardCount ) 
      continue;

    for(Int_t ich=0; ich<AliMUONConstants::NTriggerCh(); ich++){
      histoName = fEfficiencyMap->GetHistoName(currDe, AliMUONTriggerEfficiencyCells::kAllTracks, ich);
      if ( fEfficiencyMap->GetHistoList() ) {
	histoDen = (TH1F*)fEfficiencyMap->GetHistoList()->FindObject(histoName.Data());
	if ( !histoDen ) {
	  AliWarning(Form("Histogram %s not found. Efficiency won't be re-build", histoName.Data()));
	  rebuildEfficiency = kFALSE;
	}
      }
      else {
	AliWarning("Histogram list not present: efficiency won't be re-build");
	rebuildEfficiency = kFALSE;
      }

      Int_t nTypes = ( rebuildEfficiency ) ? AliMUONTriggerEfficiencyCells::kNcounts-1 : 2; 
      for(Int_t hType=0; hType<nTypes; hType++){
	histoName = fEfficiencyMap->GetHistoName(currDe, hType, ich);

	histoNum = ( rebuildEfficiency ) ? 
	  (TH1F*)fEfficiencyMap->GetHistoList()->FindObject(histoName.Data()) :
	  fEfficiencyMap->GetOldEffHisto(currDe, ich, hType);
      
	if ( !histoNum ) {
	  AliWarning(Form("Histogram %s not found. Skip to next", histoName.Data()));
	  continue;
	}

	index = GetIndex(deTypeEff[ide], hType, ich);
	TGraphAsymmErrors* effGraph = GetEfficiencyGraph(histoNum,histoDen);
	histoName.ReplaceAll("Count","Eff");
	effGraph->SetName(histoName.Data());
	fEfficiencyObjects->AddAt(effGraph, index);
	AliDebug(5,Form("Adding object %s (%s/%s) at index %i",effGraph->GetName(),histoNum->GetName(),histoDen->GetName(),index));

	if ( useMeanValues ){
	  Int_t currChamber = ich + AliMpConstants::NofTrackingChambers();
	  histoName = fEfficiencyMap->GetHistoName(AliMUONTriggerEfficiencyCells::kHboardCount, hType, ich);
	  TH1F* auxHistoNum = (TH1F*)fEfficiencyMap->GetHistoList()->FindObject(histoName.Data())->Clone("tempHistoNum");
	  TH1F* auxHistoDen = (TH1F*)fEfficiencyMap->GetHistoList()->FindObject(histoName.Data())->Clone("tempHistoDen");
	  for ( Int_t iBinBoard = 1; iBinBoard<=AliMpConstants::NofLocalBoards(); iBinBoard++){
	    Int_t detElemId = AliMpDDLStore::Instance()->GetDEfromLocalBoard(iBinBoard, currChamber);
	    Int_t iBin = histoNum->FindBin(detElemId%100);

	    auxHistoNum->SetBinContent(iBinBoard, histoNum->GetBinContent(iBin));
	    auxHistoDen->SetBinContent(iBinBoard, histoDen->GetBinContent(iBin));
	  }
	  index = GetIndex(kHboardEff, hType, ich);
	  effGraph = GetEfficiencyGraph(auxHistoNum,auxHistoDen);
	  histoName.ReplaceAll("Count","Eff");
	  effGraph->SetName(histoName.Data());
	  fEfficiencyObjects->AddAt(effGraph, index);
	  AliDebug(5,Form("Adding object %s (%s/%s) at index %i",effGraph->GetName(),histoNum->GetName(),histoDen->GetName(),index));
	  delete auxHistoNum;
	  delete auxHistoDen;
	} // if (useMeanValues)
      } // loop on count type
    } // loop on chamber
  } // loop on detection element histogram
}


//_____________________________________________________________________________
void AliMUONTriggerChamberEfficiency::DisplayEfficiency(Bool_t perSlat, Bool_t show2Dhisto)
{
  //
  /// Display calculated efficiency.
  //

  if ( !AliCDBManager::Instance()->GetDefaultStorage() ){
    AliWarning("Please set default CDB storage (needed for mapping).");
    return;
  }
  if ( AliCDBManager::Instance()->GetRun() < 0 ){
    AliWarning("Please set CDB run number (needed for mapping).");
    return;
  }

  TString baseCanName = "MTRtrigChEffCan";
  TString histoName;

  // Remove previously created canvases
  TCanvas* can = 0x0;
  TIter next(gROOT->GetListOfCanvases());
  while ((can = (TCanvas *)next())) {
    histoName = can->GetName();
    if ( histoName.Contains(baseCanName.Data()))
      delete can;
  }

  delete fDisplayList;
  fDisplayList = new TList();
  fDisplayList->SetOwner();

  TH2F* displayHisto = 0x0;

  AliMUONTriggerDisplay triggerDisplay;

  Int_t deType = ( perSlat ) ? kHslatEff : kHboardEff;
  AliMUONTriggerDisplay::EDisplayType displayType = ( perSlat ) ? 
    AliMUONTriggerDisplay::kDisplaySlats : AliMUONTriggerDisplay::kDisplayBoards;
  Int_t index = -1;

  TGraph* graph = 0x0;

  // Book histos
  for(Int_t ich=0; ich<AliMUONConstants::NTriggerCh(); ich++){
    Int_t currCh = 11 + ich;
    for(Int_t hType=0; hType<AliMUONTriggerEfficiencyCells::kNcounts - 1; hType++){
      index = GetIndex(deType, hType, ich);
      graph = (TGraph*)fEfficiencyObjects->At(index);
      if ( ! graph ) continue;
      histoName = graph->GetName();
      histoName += baseCanName;
      Int_t shift = 10*(index%((AliMUONTriggerEfficiencyCells::kNcounts - 1)*
			       AliMUONConstants::NTriggerCh()));
      can = new TCanvas(histoName.Data(), histoName.Data(), 100+shift, shift, 700, 700);
      can->SetRightMargin(0.14);
      can->SetLeftMargin(0.12);
      histoName.ReplaceAll(baseCanName.Data(), "Display");
      if ( show2Dhisto ) {
	displayHisto = 
	  (TH2F*)triggerDisplay.GetDisplayHistogram(graph, histoName,
						    displayType,
						    hType,currCh,histoName,
						    AliMUONTriggerDisplay::kShowZeroes);
	displayHisto->SetDirectory(0);
      }

      if ( show2Dhisto ){
	displayHisto->GetZaxis()->SetRangeUser(0.,1.);
	displayHisto->GetYaxis()->SetTitleOffset(1.4);
	displayHisto->SetStats(kFALSE);
	displayHisto->DrawCopy("COLZ");
	delete displayHisto;

	if ( deType == kHboardEff ){
	  histoName = Form("labels%iChamber%i", hType, currCh);
	  displayHisto = 
	    (TH2F*)triggerDisplay.GetBoardNumberHisto(histoName,currCh);
	  displayHisto->SetDirectory(0);
	  displayHisto->DrawCopy("textsame");
	  delete displayHisto;
	}
      }
      else {
	TGraphAsymmErrors* drawGraph = (TGraphAsymmErrors*)graph->Clone(histoName.Data());
	drawGraph->SetMarkerStyle(20);
	drawGraph->SetMarkerSize(0.7);
	drawGraph->SetMarkerColor(kRed);
	fDisplayList->Add(drawGraph);
	drawGraph->Draw("ap");
      } // loop on chamber
    } // loop on count type
  } // loop on chamber
}


//__________________________________________________________________________
Bool_t AliMUONTriggerChamberEfficiency::LowStatisticsSettings(Bool_t useMeanValues)
{
  //
  /// In case of low statistics, fill the local board efficiency with
  /// the average value of the RPC
  //

  if ( useMeanValues )
    AliInfo("Boards filled with the average efficiency of the RPC");

  FillFromList(useMeanValues);

  return kTRUE;
}


//__________________________________________________________________________
Int_t
AliMUONTriggerChamberEfficiency::GetIndex(Int_t histoType, Int_t countType, 
					  Int_t chamber) const
{
  //
  /// Return the index of the object in the array
  //

  const Int_t kNtypes = AliMUONTriggerEfficiencyCells::kNcounts - 1;
  const Int_t kNchambers = AliMUONConstants::NTriggerCh();
  return 
    histoType * kNtypes * kNchambers + 
    chamber * kNtypes +
    countType;
  //countType * kNchambers +
  //chamber;
}

//_____________________________________________________________________________
TGraphAsymmErrors* AliMUONTriggerChamberEfficiency::GetEfficiencyGraph(TH1* histoNum, TH1* histoDen)
{
  //
  /// Create the graph of efficiency from the numerator and denominator
  /// histogram in such a way to have a point set also for
  /// detection elements with efficiency = 0 or non calculated
  //

  TGraphAsymmErrors* auxGraph = 0x0;
  if ( histoDen ) auxGraph = new TGraphAsymmErrors(histoNum,histoDen,"cp");
  else auxGraph = new TGraphAsymmErrors(histoNum);

  Int_t npoints = histoNum->GetNbinsX();
  TGraphAsymmErrors* effGraph = new TGraphAsymmErrors(npoints);
  Double_t oldX, oldY;
  for ( Int_t ibin=0; ibin<npoints; ibin++ ) {
    Int_t foundPoint = -1;
    for (Int_t ipt=0; ipt<auxGraph->GetN(); ipt++) {
      auxGraph->GetPoint(ipt, oldX, oldY);
      if ( oldX > histoNum->GetBinLowEdge(ibin+1) &&
	   oldX < histoNum->GetBinLowEdge(ibin+2) ) {
	foundPoint = ipt;
	break;
      }
    }
    Double_t currX = ( foundPoint < 0 ) ? histoNum->GetBinCenter(ibin+1) : oldX; 
    Double_t currY = ( foundPoint < 0 ) ? 0. : oldY;
    Double_t eyl   = ( foundPoint < 0 ) ? 0. : auxGraph->GetErrorYlow(foundPoint);
    Double_t eyh   = ( foundPoint < 0 ) ? 0. : auxGraph->GetErrorYhigh(foundPoint);
    effGraph->SetPoint(ibin, currX, currY);
    effGraph->SetPointError(ibin, 0., 0., eyl, eyh);
  }

  delete auxGraph;

  return effGraph;
}
Commit	Line	Data
f3d288ab	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	#include "AliMUONTriggerEfficiencyCells.h"
	19	#include "AliMpConstants.h"
	20	#include "AliMUONConstants.h"
	21
	22	// Classes for display
	23	#include "AliMUONTriggerDisplay.h"
	24	#include "AliCDBManager.h"
	25	#include "AliMpDDLStore.h"
	26
	27	#include "AliLog.h"
	28
	29	#include "TRandom.h"
	30	#include "Riostream.h"
	31	#include "TH1F.h"
	32	#include "TObjArray.h"
	33	#include "TGraphAsymmErrors.h"
	34
	35	#include "TH2F.h"
	36	#include "TCanvas.h"
	37	#include "TROOT.h"
	38
	39	#include "AliMUONTriggerChamberEfficiency.h"
	40
	41	//-----------------------------------------------------------------------------
	42	/// \class AliMUONTriggerChamberEfficiency
	43	/// A class to store and give access to the trigger chamber efficiency.
	44	///
	45	/// Efficiency is stored per cathode on local boards
	46	///
	47	/// The main method of this class is IsTriggered().
	48	///
	49	/// \author Diego Stocco; INFN Torino
	50	//-----------------------------------------------------------------------------
	51
	52	/// \cond CLASSIMP
	53	ClassImp(AliMUONTriggerChamberEfficiency)
	54	/// \endcond
	55
	56	//__________________________________________________________________________
	57	AliMUONTriggerChamberEfficiency::AliMUONTriggerChamberEfficiency(AliMUONTriggerEfficiencyCells* effCells)
	58	:
	59	TObject(),
	60	fIsOwner(kFALSE),
	61	fEfficiencyMap(effCells),
	62	fEfficiencyObjects(0x0),
	63	fDisplayList(0x0)
	64	{
65	/// Default constructor.
66	FillFromList();
67	}
68
69	//__________________________________________________________________________
70	AliMUONTriggerChamberEfficiency::AliMUONTriggerChamberEfficiency(const Char_t* filename, const Char_t* listname)
71	:
72	TObject(),
73	fIsOwner(kTRUE),
74	fEfficiencyMap(0x0),
75	fEfficiencyObjects(0x0),
76	fDisplayList(0x0)
77	{
78	/// Constructor using an ASCII file.
79	fEfficiencyMap = new AliMUONTriggerEfficiencyCells(filename, listname);
80	FillFromList();
81	}
82
83
84	//_____________________________________________________________________________
85	AliMUONTriggerChamberEfficiency::AliMUONTriggerChamberEfficiency(const AliMUONTriggerChamberEfficiency& other)
86	:
87	TObject(other),
88	fIsOwner(other.fIsOwner),
89	fEfficiencyMap(other.fEfficiencyMap),
90	fEfficiencyObjects(other.fEfficiencyObjects),
91	fDisplayList(other.fDisplayList)
92	{
93	/// Copy constructor
94	}
95
96	//_____________________________________________________________________________
97	AliMUONTriggerChamberEfficiency& AliMUONTriggerChamberEfficiency::operator=(const AliMUONTriggerChamberEfficiency& other)
98	{
99	/// Asignment operator
100	// check assignement to self
101	if (this == &other)
102	return *this;
103
104	fIsOwner = other.fIsOwner;
105	fEfficiencyMap = other.fEfficiencyMap;
106	fEfficiencyObjects = other.fEfficiencyObjects;
107	fDisplayList = other.fDisplayList;
108
109	return *this;
110	}
111
112	//__________________________________________________________________________
113	AliMUONTriggerChamberEfficiency::~AliMUONTriggerChamberEfficiency()
114	{
115	/// Destructor.
116	if ( fIsOwner )
117	delete fEfficiencyMap;
118	delete fEfficiencyObjects;
119	delete fDisplayList;
120	}
121
122
123	//__________________________________________________________________________
124	Float_t AliMUONTriggerChamberEfficiency::GetCellEfficiency(Int_t detElemId, Int_t localBoard, Int_t hType) const
125	{
126	/// Get the efficiencies of the 2 cathodes at a given local board
127
128	Int_t chamber = FindChamberIndex(detElemId);
129	Int_t index = GetIndex(kHboardEff, hType, chamber);
130	TGraphAsymmErrors* effGraph = ((TGraphAsymmErrors*)fEfficiencyObjects->At(index));
131
132	// Some graphs are not available in the old implementation
133	if ( ! effGraph ) return -1.;
134
135	Double_t xpt, ypt;
136	effGraph->GetPoint(localBoard-1, xpt, ypt);
137	return ypt;
138	}
139
140
141	//__________________________________________________________________________
142	void
143	AliMUONTriggerChamberEfficiency::IsTriggered(Int_t detElemId, Int_t localBoard, Bool_t &trigBend, Bool_t &trigNonBend) const
144	{
145	/// Whether or not a given local board has a chance to trig, on each cathode.
146
147	// P(B) : probability to fire bending plane
148	Float_t effBend = GetCellEfficiency(detElemId, localBoard, AliMUONTriggerEfficiencyCells::kBendingEff);
149
150	// P(BN) : probability to fire bending and non-bending plane
151	Float_t effBoth = GetCellEfficiency(detElemId, localBoard, AliMUONTriggerEfficiencyCells::kBothPlanesEff);
152
153	trigBend = ( gRandom->Rndm() > effBend ) ? kFALSE : kTRUE;
154
155	// P(N) : probability to fire non-bending plane
156	Float_t effNonBend = GetCellEfficiency(detElemId, localBoard, AliMUONTriggerEfficiencyCells::kNonBendingEff);
157
158	if ( effBoth > 0 ) {
159	effNonBend = ( trigBend ) ?
160	effBoth / effBend : // P(N\|B) = P(BN) / P(B)
161	( effNonBend - effBoth ) / ( 1. - effBend ); // P(N\|!B) = ( P(N) - P(BN) ) / ( 1 - P(B) )
162	}
163
164	trigNonBend = ( gRandom->Rndm() > effNonBend ) ? kFALSE : kTRUE;
165
166	AliDebug(2,Form("Ch %i board %i resp (%i, %i) prob (%.2f, %.2f) effNB %.2f effBoth %.2f\n", detElemId/100, localBoard, trigBend, trigNonBend, effBend, effNonBend, GetCellEfficiency(detElemId, localBoard, AliMUONTriggerEfficiencyCells::kNonBendingEff), effBoth));
167
168
169	}
170
171
172	//__________________________________________________________________________
173	Int_t AliMUONTriggerChamberEfficiency::FindChamberIndex(Int_t detElemId) const
174	{
175	/// From detElemId to chamber number
176
177	// Int_t iChamber = AliMpDEManager::GetChamberId(detElemId);
178	Int_t iChamber = detElemId/100 - 1;
179	return iChamber-AliMpConstants::NofTrackingChambers();
180	}
181
182
183	//__________________________________________________________________________
184	void
185	AliMUONTriggerChamberEfficiency::FillFromList(Bool_t useMeanValues)
186	{
187	/// Fills internal histos from list.
188
189	if ( fEfficiencyObjects )
190	delete fEfficiencyObjects;
191
192	const Int_t kNeffHistos =
193	2 * ( AliMUONTriggerEfficiencyCells::kNcounts - 1 ) * AliMUONConstants::NTriggerCh();
194
195	fEfficiencyObjects = new TObjArray(kNeffHistos);
196	fEfficiencyObjects->SetOwner();
197
198	TH1F histoNum = 0x0, histoDen=0x0;
199	TString histoName = "";
200	Int_t deType[2] = {AliMUONTriggerEfficiencyCells::kHboardCount,
201	AliMUONTriggerEfficiencyCells::kHslatCount};
202	Int_t deTypeEff[2] = {kHboardEff, kHslatEff};
203	Int_t index = -1;
204
205	Bool_t rebuildEfficiency = kTRUE;
206
207	for ( Int_t ide=0; ide<2; ide++){
208	Int_t currDe = deType[ide];
209
210	if ( useMeanValues && currDe == AliMUONTriggerEfficiencyCells::kHboardCount )
211	continue;
212
213	for(Int_t ich=0; ich<AliMUONConstants::NTriggerCh(); ich++){
214	histoName = fEfficiencyMap->GetHistoName(currDe, AliMUONTriggerEfficiencyCells::kAllTracks, ich);
215	if ( fEfficiencyMap->GetHistoList() ) {
216	histoDen = (TH1F*)fEfficiencyMap->GetHistoList()->FindObject(histoName.Data());
217	if ( !histoDen ) {
218	AliWarning(Form("Histogram %s not found. Efficiency won't be re-build", histoName.Data()));
219	rebuildEfficiency = kFALSE;
220	}
221	}
222	else {
223	AliWarning("Histogram list not present: efficiency won't be re-build");
224	rebuildEfficiency = kFALSE;
225	}
226
227	Int_t nTypes = ( rebuildEfficiency ) ? AliMUONTriggerEfficiencyCells::kNcounts-1 : 2;
228	for(Int_t hType=0; hType<nTypes; hType++){
229	histoName = fEfficiencyMap->GetHistoName(currDe, hType, ich);
230
231	histoNum = ( rebuildEfficiency ) ?
232	(TH1F*)fEfficiencyMap->GetHistoList()->FindObject(histoName.Data()) :
233	fEfficiencyMap->GetOldEffHisto(currDe, ich, hType);
234
235	if ( !histoNum ) {
236	AliWarning(Form("Histogram %s not found. Skip to next", histoName.Data()));
237	continue;
238	}
239
240	index = GetIndex(deTypeEff[ide], hType, ich);
241	TGraphAsymmErrors* effGraph = GetEfficiencyGraph(histoNum,histoDen);
242	histoName.ReplaceAll("Count","Eff");
243	effGraph->SetName(histoName.Data());
244	fEfficiencyObjects->AddAt(effGraph, index);
245	AliDebug(5,Form("Adding object %s (%s/%s) at index %i",effGraph->GetName(),histoNum->GetName(),histoDen->GetName(),index));
246
247	if ( useMeanValues ){
248	Int_t currChamber = ich + AliMpConstants::NofTrackingChambers();
249	histoName = fEfficiencyMap->GetHistoName(AliMUONTriggerEfficiencyCells::kHboardCount, hType, ich);
250	TH1F* auxHistoNum = (TH1F*)fEfficiencyMap->GetHistoList()->FindObject(histoName.Data())->Clone("tempHistoNum");
251	TH1F* auxHistoDen = (TH1F*)fEfficiencyMap->GetHistoList()->FindObject(histoName.Data())->Clone("tempHistoDen");
252	for ( Int_t iBinBoard = 1; iBinBoard<=AliMpConstants::NofLocalBoards(); iBinBoard++){
253	Int_t detElemId = AliMpDDLStore::Instance()->GetDEfromLocalBoard(iBinBoard, currChamber);
254	Int_t iBin = histoNum->FindBin(detElemId%100);
255
256	auxHistoNum->SetBinContent(iBinBoard, histoNum->GetBinContent(iBin));
257	auxHistoDen->SetBinContent(iBinBoard, histoDen->GetBinContent(iBin));
258	}
259	index = GetIndex(kHboardEff, hType, ich);
260	effGraph = GetEfficiencyGraph(auxHistoNum,auxHistoDen);
261	histoName.ReplaceAll("Count","Eff");
262	effGraph->SetName(histoName.Data());
263	fEfficiencyObjects->AddAt(effGraph, index);
264	AliDebug(5,Form("Adding object %s (%s/%s) at index %i",effGraph->GetName(),histoNum->GetName(),histoDen->GetName(),index));
265	delete auxHistoNum;
266	delete auxHistoDen;
267	} // if (useMeanValues)
268	} // loop on count type
269	} // loop on chamber
270	} // loop on detection element histogram
271	}
272
273
274	//_____________________________________________________________________________
275	void AliMUONTriggerChamberEfficiency::DisplayEfficiency(Bool_t perSlat, Bool_t show2Dhisto)
276	{
277	//
278	/// Display calculated efficiency.
279	//
280
281	if ( !AliCDBManager::Instance()->GetDefaultStorage() ){
282	AliWarning("Please set default CDB storage (needed for mapping).");
283	return;
284	}
285	if ( AliCDBManager::Instance()->GetRun() < 0 ){
286	AliWarning("Please set CDB run number (needed for mapping).");
287	return;
288	}
289
290	TString baseCanName = "MTRtrigChEffCan";
291	TString histoName;
292
293	// Remove previously created canvases
294	TCanvas* can = 0x0;
295	TIter next(gROOT->GetListOfCanvases());
296	while ((can = (TCanvas *)next())) {
297	histoName = can->GetName();
298	if ( histoName.Contains(baseCanName.Data()))
299	delete can;
300	}
301
302	delete fDisplayList;
303	fDisplayList = new TList();
304	fDisplayList->SetOwner();
305
306	TH2F* displayHisto = 0x0;
307
308	AliMUONTriggerDisplay triggerDisplay;
309
310	Int_t deType = ( perSlat ) ? kHslatEff : kHboardEff;
311	AliMUONTriggerDisplay::EDisplayType displayType = ( perSlat ) ?
312	AliMUONTriggerDisplay::kDisplaySlats : AliMUONTriggerDisplay::kDisplayBoards;
313	Int_t index = -1;
314
315	TGraph* graph = 0x0;
316
317	// Book histos
318	for(Int_t ich=0; ich<AliMUONConstants::NTriggerCh(); ich++){
319	Int_t currCh = 11 + ich;
320	for(Int_t hType=0; hType<AliMUONTriggerEfficiencyCells::kNcounts - 1; hType++){
321	index = GetIndex(deType, hType, ich);
322	graph = (TGraph*)fEfficiencyObjects->At(index);
323	if ( ! graph ) continue;
324	histoName = graph->GetName();
325	histoName += baseCanName;
326	Int_t shift = 10(index%((AliMUONTriggerEfficiencyCells::kNcounts - 1)
327	AliMUONConstants::NTriggerCh()));
328	can = new TCanvas(histoName.Data(), histoName.Data(), 100+shift, shift, 700, 700);
329	can->SetRightMargin(0.14);
330	can->SetLeftMargin(0.12);
331	histoName.ReplaceAll(baseCanName.Data(), "Display");
332	if ( show2Dhisto ) {
333	displayHisto =
334	(TH2F*)triggerDisplay.GetDisplayHistogram(graph, histoName,
335	displayType,
336	hType,currCh,histoName,
337	AliMUONTriggerDisplay::kShowZeroes);
338	displayHisto->SetDirectory(0);
339	}
340
341	if ( show2Dhisto ){
342	displayHisto->GetZaxis()->SetRangeUser(0.,1.);
343	displayHisto->GetYaxis()->SetTitleOffset(1.4);
344	displayHisto->SetStats(kFALSE);
345	displayHisto->DrawCopy("COLZ");
346	delete displayHisto;
347
348	if ( deType == kHboardEff ){
349	histoName = Form("labels%iChamber%i", hType, currCh);
350	displayHisto =
351	(TH2F*)triggerDisplay.GetBoardNumberHisto(histoName,currCh);
352	displayHisto->SetDirectory(0);
353	displayHisto->DrawCopy("textsame");
354	delete displayHisto;
355	}
356	}
357	else {
358	TGraphAsymmErrors* drawGraph = (TGraphAsymmErrors*)graph->Clone(histoName.Data());
359	drawGraph->SetMarkerStyle(20);
360	drawGraph->SetMarkerSize(0.7);
361	drawGraph->SetMarkerColor(kRed);
362	fDisplayList->Add(drawGraph);
363	drawGraph->Draw("ap");
364	} // loop on chamber
365	} // loop on count type
366	} // loop on chamber
367	}
368
369
370	//__________________________________________________________________________
371	Bool_t AliMUONTriggerChamberEfficiency::LowStatisticsSettings(Bool_t useMeanValues)
372	{
373	//
374	/// In case of low statistics, fill the local board efficiency with
375	/// the average value of the RPC
376	//
377
378	if ( useMeanValues )
379	AliInfo("Boards filled with the average efficiency of the RPC");
380
381	FillFromList(useMeanValues);
382
383	return kTRUE;
384	}
385
386
387	//__________________________________________________________________________
388	Int_t
389	AliMUONTriggerChamberEfficiency::GetIndex(Int_t histoType, Int_t countType,
390	Int_t chamber) const
391	{
392	//
393	/// Return the index of the object in the array
394	//
395
396	const Int_t kNtypes = AliMUONTriggerEfficiencyCells::kNcounts - 1;
397	const Int_t kNchambers = AliMUONConstants::NTriggerCh();
398	return
399	histoType * kNtypes * kNchambers +
400	chamber * kNtypes +
401	countType;
402	//countType * kNchambers +
403	//chamber;
404	}
405
406	//_____________________________________________________________________________
407	TGraphAsymmErrors* AliMUONTriggerChamberEfficiency::GetEfficiencyGraph(TH1* histoNum, TH1* histoDen)
408	{
409	//
410	/// Create the graph of efficiency from the numerator and denominator
411	/// histogram in such a way to have a point set also for
412	/// detection elements with efficiency = 0 or non calculated
413	//
414
415	TGraphAsymmErrors* auxGraph = 0x0;
4784d180	416	if ( histoDen ) auxGraph = new TGraphAsymmErrors(histoNum,histoDen,"cp");
f3d288ab	417	else auxGraph = new TGraphAsymmErrors(histoNum);
	418
	419	Int_t npoints = histoNum->GetNbinsX();
	420	TGraphAsymmErrors* effGraph = new TGraphAsymmErrors(npoints);
	421	Double_t oldX, oldY;
	422	for ( Int_t ibin=0; ibin<npoints; ibin++ ) {
	423	Int_t foundPoint = -1;
	424	for (Int_t ipt=0; ipt<auxGraph->GetN(); ipt++) {
	425	auxGraph->GetPoint(ipt, oldX, oldY);
	426	if ( oldX > histoNum->GetBinLowEdge(ibin+1) &&
	427	oldX < histoNum->GetBinLowEdge(ibin+2) ) {
	428	foundPoint = ipt;
	429	break;
	430	}
	431	}
	432	Double_t currX = ( foundPoint < 0 ) ? histoNum->GetBinCenter(ibin+1) : oldX;
	433	Double_t currY = ( foundPoint < 0 ) ? 0. : oldY;
	434	Double_t eyl = ( foundPoint < 0 ) ? 0. : auxGraph->GetErrorYlow(foundPoint);
	435	Double_t eyh = ( foundPoint < 0 ) ? 0. : auxGraph->GetErrorYhigh(foundPoint);
	436	effGraph->SetPoint(ibin, currX, currY);
	437	effGraph->SetPointError(ibin, 0., 0., eyl, eyh);
	438	}
	439
	440	delete auxGraph;
	441
	442	return effGraph;
	443	}