#include "AliMUONTriggerEfficiencyCells.h"
#include "AliMpConstants.h"
-#include "AliMpDEManager.h"
+
+// Classes for display
+#include "AliMUONGeometryTransformer.h"
+#include "AliMpCDB.h"
+#include "AliCDBManager.h"
+#include "AliMpDDLStore.h"
+#include "AliMpDDL.h"
+#include "AliMpTriggerCrate.h"
+#include "AliMpLocalBoard.h"
+#include "AliMpPad.h"
+#include "AliMpVSegmentation.h"
+#include "AliMpSegmentation.h"
#include "AliLog.h"
#include "TMath.h"
#include "TH2F.h"
+#include "TH3F.h"
+#include "TF1.h"
#include "TStyle.h"
#include "TPaveLabel.h"
#include "TCanvas.h"
/// \class AliMUONTriggerEfficiencyCells
/// A class to store and give access to the trigger chamber efficiency.
///
-/// Efficiency is stored per cathode on local boards, or, alternatively,
-/// on "cells" of a given size.
+/// Efficiency is stored per cathode on local boards
///
/// The main method of this class is IsTriggered().
///
-/// $ALICE_ROOT/MUON/data/TriggerChamberefficiencyCells.dat contains efficiency
+/// $ALICE_ROOT/MUON/data/efficiencyCells.dat contains efficiency
/// for each chamber (i.e. DetElement).
///
/// In the case of local boards, efficiency is stored from left to right
/// per increasing board number (from 1 to 234)
///
-/// Otherwise, he efficiency cells goes from right to left and
-/// from bottom to top of the chamber, namely, the efficiencies tabulated in the
-/// file refers to the following reference frame:
-///
-/// <pre>
-/// x
-/// <----------------------------------|
-/// |
-/// --------------------------- |
-/// | 0.97 | 0.97 | 0.97 | 0.97 | |
-/// --------------------------- |
-/// | 0.97 | 0.97 | 0.97 | 0.97 | |
-/// --------------------------- |
-/// | 0.97 | 0.97 | 0.97 | 0.97 | |
-/// --------------------------- |
-/// |
-/// \/ y
-/// </pre>
-///
-/// In both cases, the file can be edited in order to change efficiency
-/// in a chosen local board/region of the chamber.
+/// The file can be edited in order to change efficiency
+/// in a chosen local board/region of the chamber.
///
///
/// But please note that this object is also available from the CDB
-/// (generated using the MUONCDB.C macro)
///
/// \author Diego Stocco; INFN Torino
//-----------------------------------------------------------------------------
:
TObject(),
fCountHistoList(0x0),
-fNoCountHistoList(0x0)
+fNoCountHistoList(0x0),
+fFiredStrips(0x0),
+fDisplayHistoList(0x0),
+fBoardLabelList(0x0),
+fFiredFitHistoList(0x0),
+fFiredDisplayHistoList(0x0)
{
/// Default constructor.
CheckConstants();
}
//__________________________________________________________________________
-AliMUONTriggerEfficiencyCells::AliMUONTriggerEfficiencyCells(const char* filename)
+AliMUONTriggerEfficiencyCells::AliMUONTriggerEfficiencyCells(const Char_t* filename)
:
TObject(),
fCountHistoList(0x0),
-fNoCountHistoList(0x0)
+fNoCountHistoList(0x0),
+fFiredStrips(0x0),
+fDisplayHistoList(0x0),
+fBoardLabelList(0x0),
+fFiredFitHistoList(0x0),
+fFiredDisplayHistoList(0x0)
{
/// Constructor using an ASCII file.
CheckConstants();
:
TObject(),
fCountHistoList(countHistoList),
-fNoCountHistoList(noCountHistoList)
+fNoCountHistoList(noCountHistoList),
+fFiredStrips(0x0),
+fDisplayHistoList(0x0),
+fBoardLabelList(0x0),
+fFiredFitHistoList(0x0),
+fFiredDisplayHistoList(0x0)
{
/// Constructor using an ASCII file.
CheckConstants();
:
TObject(other),
fCountHistoList(other.fCountHistoList),
-fNoCountHistoList(other.fNoCountHistoList)
+fNoCountHistoList(other.fNoCountHistoList),
+fFiredStrips(other.fFiredStrips),
+fDisplayHistoList(other.fDisplayHistoList),
+fBoardLabelList(other.fBoardLabelList),
+fFiredFitHistoList(other.fFiredFitHistoList),
+fFiredDisplayHistoList(other.fFiredDisplayHistoList)
{
/// Copy constructor
for(Int_t chCath=0; chCath<fgkNplanes; chCath++){
- for(Int_t slat=0; slat<fgkNslats; slat++){
- fCellContent[chCath][slat] = other.fCellContent[chCath][slat];
- }
- fCellSize[chCath] = other.fCellSize[chCath];
- fCellNumber[chCath] = other.fCellNumber[chCath];
fBoardEfficiency[chCath] = other.fBoardEfficiency[chCath];
fSlatEfficiency[chCath] = other.fSlatEfficiency[chCath];
}
return *this;
for(Int_t chCath=0; chCath<fgkNplanes; chCath++){
- for(Int_t slat=0; slat<fgkNslats; slat++){
- fCellContent[chCath][slat] = other.fCellContent[chCath][slat];
- }
- fCellSize[chCath] = other.fCellSize[chCath];
- fCellNumber[chCath] = other.fCellNumber[chCath];
fBoardEfficiency[chCath] = other.fBoardEfficiency[chCath];
fSlatEfficiency[chCath] = other.fSlatEfficiency[chCath];
}
fCountHistoList = other.fCountHistoList;
fNoCountHistoList = other.fNoCountHistoList;
+ fFiredStrips = other.fFiredStrips;
+ fDisplayHistoList = other.fDisplayHistoList;
+ fBoardLabelList = other.fBoardLabelList;
+ fFiredFitHistoList = other.fFiredFitHistoList;
+ fFiredDisplayHistoList = other.fFiredDisplayHistoList;
+
return *this;
}
}
-//__________________________________________________________________________
-void AliMUONTriggerEfficiencyCells::GetCellEfficiency(Int_t detElemId, Float_t x, Float_t y, Float_t &eff1, Float_t &eff2) const
-{
-/// Get the efficiencies of the 2 cathodes at a given location (x,y)
-
- Int_t chamber = FindChamberIndex(detElemId);
- Int_t slat = FindSlatIndex(detElemId);
- TArrayI cell = CellByCoord(detElemId,x,y);
- eff1 = 0.0;
- eff2 = 0.0;
- if(cell.At(0)>=0 && cell.At(1)>=0)
- {
- eff1 = fCellContent[chamber][slat][cell.At(0)][cell.At(1)];
- eff2 = fCellContent[fgkNchambers+chamber][slat][cell.At(0)][cell.At(1)];
- }
-}
-
-
//__________________________________________________________________________
void AliMUONTriggerEfficiencyCells::GetCellEfficiency(Int_t detElemId, Int_t localBoard, Float_t &eff1, Float_t &eff2) const
{
}
-//__________________________________________________________________________
-void
-AliMUONTriggerEfficiencyCells::IsTriggered(Int_t detElemId, Float_t x, Float_t y, Bool_t &trig1, Bool_t &trig2) const
-{
-/// Whether or not a given location (x,y) has a chance to trig, on each cathode.
-
- Float_t eff1 = 0.0;
- Float_t eff2 = 0.0;
- GetCellEfficiency(detElemId, x, y, eff1, eff2);
- trig1 = kTRUE;
- trig2 = kTRUE;
- if(gRandom->Rndm()>eff1)trig1 = kFALSE;
- if(gRandom->Rndm()>eff2)trig2 = kFALSE;
-}
-
-
//__________________________________________________________________________
void
AliMUONTriggerEfficiencyCells::IsTriggered(Int_t detElemId, Int_t localBoard, Bool_t &trig1, Bool_t &trig2) const
if(gRandom->Rndm()>eff2)trig2 = kFALSE;
}
-//__________________________________________________________________________
-TArrayI AliMUONTriggerEfficiencyCells::CellByCoord(Int_t detElemId, Float_t x, Float_t y) const
-{
-/// Get the efficiencies at a given location.
-
- Int_t chamber = FindChamberIndex(detElemId);
- Int_t slat = FindSlatIndex(detElemId);
- Int_t cell[fgkNcathodes]={-1,-1};
- Float_t maxX = fCellSize[chamber][slat]*((Float_t)fCellNumber[chamber][slat]);
- Float_t maxY = fCellSize[fgkNchambers+chamber][slat]*((Float_t)fCellNumber[fgkNchambers+chamber][slat]);
- if(x>=0 & x<maxX & y>=0 & y<maxY)
- {
- cell[0] = (Int_t)(x/fCellSize[chamber][slat]);
- cell[1] = (Int_t)(y/fCellSize[fgkNchambers+chamber][slat]);
- }
- return TArrayI(fgkNcathodes,cell);
-}
//__________________________________________________________________________
-void AliMUONTriggerEfficiencyCells::ReadFile(const char* filename)
+void AliMUONTriggerEfficiencyCells::ReadFile(const Char_t* filename)
{
/// Reads a file containing the efficiency map.
InitHistos();
ifstream file(fileName.Data());
- char dat[50];
+ Char_t dat[50];
if (file.good()){
file >> dat;
if(!strcmp(dat,"localBoards"))ReadFileBoards(file);
- else ReadFileXY(file);
+ else AliWarning("File .dat in wrong format");
file.close();
} else {
AliWarning(Form("Can't read file %s",fileName.Data()));
}
-//__________________________________________________________________________
-void AliMUONTriggerEfficiencyCells::ReadFileXY(ifstream &file)
-{
-/// Structure of file (.dat) containing geometrical efficency
- Int_t datInt=0, detEl=0, chamber=0, rpc=0, chCath=0;
- Float_t datFloat=0.0;
- char dat[50];
-
- while (file >> dat) {
- file >> detEl;
- chamber = FindChamberIndex(detEl);
- rpc = FindSlatIndex(detEl);
- file >> dat;
- for(Int_t i=0; i<fgkNcathodes; i++){
- chCath = fgkNchambers*i + chamber;
- file >> datInt;
- fCellNumber[chCath][rpc] = datInt;
- file >> dat;
- }
- for(Int_t i=0; i<fgkNcathodes; i++){
- chCath = fgkNchambers*i + chamber;
- file >> datFloat;
- fCellSize[chCath][rpc] = datFloat;
- if(i==0)file >> dat;
- }
- for(Int_t cath=0; cath<fgkNcathodes; cath++){
- chCath = fgkNchambers*cath + chamber;
- file >> dat;
- file >> datInt;
- for(Int_t iy=0; iy<fCellNumber[fgkNchambers+chamber][rpc]; iy++){
- for(Int_t ix=0; ix<fCellNumber[chamber][rpc]; ix++){
- file >> datFloat;
- fCellContent[chCath][rpc][ix][iy] = datFloat;
- }
- }
- }
- }
-}
-
//__________________________________________________________________________
void AliMUONTriggerEfficiencyCells::ReadFileBoards(ifstream &file)
{
/// Structure of file (.dat) containing local board efficency
Int_t datInt=0, detEl=0, chamber=0, chCath=0, bin=0;
Float_t datFloat=0.0;
- char dat[50];
+ Char_t dat[50];
while (file >> dat) {
file >> detEl;
chamber = FindChamberIndex(detEl);
- //rpc = FindSlatIndex(detEl);
for(Int_t cath=0; cath<fgkNcathodes; cath++){
chCath = fgkNchambers*cath + chamber;
file >> dat;
//__________________________________________________________________________
-void AliMUONTriggerEfficiencyCells::ReadHistoBoards(const char *filename)
+void AliMUONTriggerEfficiencyCells::ReadHistoBoards(const Char_t *filename)
{
/// Structure of file (.root) containing local board efficency
TFile *file = new TFile(filename, "read");
AliWarning(Form("Can't read file %s",filename));
return;
}
- char histoName[30];
- char *cathCode[fgkNcathodes] = {"bendPlane", "nonBendPlane"};
+ Char_t histoName[40];
+ Char_t *cathCode[fgkNcathodes] = {"bendPlane", "nonBendPlane"};
+ enum {kAllChEff, kChNonEff, kNumOfHistoTypes};
+ Char_t *histoTypeName[2] = {"CountInCh", "NonCountInCh"};
- for(Int_t ch=0; ch<fgkNchambers; ch++){
- for(Int_t cath=0; cath<fgkNcathodes; cath++){
- sprintf(histoName, "%sBoardEffChamber%i", cathCode[cath], 11+ch);
- if(!(TH1F *)file->Get(histoName)) {
- AliWarning(Form("Can't find histo %s in file %s",histoName, filename));
- continue;
- }
- Int_t chCath = fgkNchambers*cath + ch;
- fBoardEfficiency[chCath] = (TH1F *)file->Get(histoName);
+ if(!fCountHistoList) fCountHistoList = new TList();
+ else fCountHistoList->Delete();
+ if(!fNoCountHistoList) fNoCountHistoList = new TList();
+ else fNoCountHistoList->Delete();
+
+ TList *currList[2] = {fCountHistoList, fNoCountHistoList};
+
+ TH1F *histo = 0x0;
+
+ for(Int_t cath=0; cath<fgkNcathodes; cath++){
+ for(Int_t hType=0; hType<kNumOfHistoTypes; hType++){
+ sprintf(histoName, "%sChamber%s", cathCode[cath], histoTypeName[hType]);
+ histo = (TH1F*)file->Get(histoName);
+ currList[hType]->Add(histo);
+ }
+ }
+
+ for(Int_t cath=0; cath<fgkNcathodes; cath++){
+ for(Int_t ch=0; ch<fgkNchambers; ch++){
+ for(Int_t hType=0; hType<kNumOfHistoTypes; hType++){
+ sprintf(histoName, "%sSlat%s%i", cathCode[cath], histoTypeName[hType], 11+ch);
+ histo = (TH1F*)file->Get(histoName);
+ currList[hType]->Add(histo);
}
+ }
}
+
+ for(Int_t cath=0; cath<fgkNcathodes; cath++){
+ for(Int_t ch=0; ch<fgkNchambers; ch++){
+ for(Int_t hType=0; hType<kNumOfHistoTypes; hType++){
+ sprintf(histoName, "%sBoard%s%i", cathCode[cath], histoTypeName[hType], 11+ch);
+ histo = (TH1F*)file->Get(histoName);
+ currList[hType]->Add(histo);
+ }
+ }
+ }
+
+ InitHistos();
+ FillHistosFromList();
}
Int_t AliMUONTriggerEfficiencyCells::FindChamberIndex(Int_t detElemId) const
{
/// From detElemId to chamber number
- Int_t iChamber = AliMpDEManager::GetChamberId(detElemId);
+
+ // Int_t iChamber = AliMpDEManager::GetChamberId(detElemId);
+ Int_t iChamber = detElemId/100 - 1;
Int_t chamber = iChamber-AliMpConstants::NofTrackingChambers();
return chamber;
}
-//__________________________________________________________________________
-Int_t AliMUONTriggerEfficiencyCells::FindSlatIndex(Int_t detElemId) const
-{
-/// From detElemId to slat index.
- Int_t slat = detElemId%100;
- return slat;
-}
-
-
-//__________________________________________________________________________
-TVector2 AliMUONTriggerEfficiencyCells::ChangeReferenceFrame(Float_t x, Float_t y, Float_t x0, Float_t y0)
-{
-/// (x0,y0) position of the local reference frame (center of the chamber)
-
- Float_t x1 = x0-x;//reflection of axis
- Float_t y1 = y+y0;
- return TVector2(x1,y1);
-}
-
//__________________________________________________________________________
void
AliMUONTriggerEfficiencyCells::Reset()
{
/// Sets our internal array contents to zero.
- for(Int_t chCath=0; chCath<fgkNplanes; chCath++){
- fCellSize[chCath].Set(fgkNslats);
- fCellNumber[chCath].Set(fgkNslats);
- for(Int_t slat=0; slat<fgkNslats; slat++){
- fCellContent[chCath][slat].ResizeTo(fgkNcells,fgkNcells);
- }
- }
-
- for(Int_t chCath=0; chCath<fgkNplanes; chCath++){
- fCellSize[chCath].Reset();
- fCellNumber[chCath].Reset();
- for(Int_t slat=0; slat<fgkNslats; slat++){
- fCellContent[chCath][slat].Zero();
- }
- }
-
- for(Int_t chCath=0; chCath<fgkNplanes; chCath++){
- fBoardEfficiency[chCath] = 0x0;
- fSlatEfficiency[chCath] = 0x0;
- }
-
+ for(Int_t chCath=0; chCath<fgkNplanes; chCath++){
+ fBoardEfficiency[chCath] = 0x0;
+ fSlatEfficiency[chCath] = 0x0;
+ }
}
/// Sets our internal array contents to zero.
const Int_t kNumOfBoards = AliMpConstants::NofLocalBoards();
+ const Int_t kNslats = 18;
Int_t chCath=0;
Char_t histoName[40];
sprintf(histoName, "%sBoardEffChamber%i", cathCode[cath], 11+ch);
fBoardEfficiency[chCath] = new TH1F(histoName, histoName, kNumOfBoards, 1-0.5, kNumOfBoards+1.-0.5);
sprintf(histoName, "%sSlatEffChamber%i", cathCode[cath], 11+ch);
- fSlatEfficiency[chCath] = new TH1F(histoName, histoName, fgkNslats, 0-0.5, fgkNslats-0.5);
+ fSlatEfficiency[chCath] = new TH1F(histoName, histoName, kNslats, 0-0.5, kNslats-0.5);
}
}
}
void
AliMUONTriggerEfficiencyCells::FillHistosFromList()
{
-/// Sets our internal array contents to zero.
+/// Fills internal histos from list.
Int_t nHistoBins=0;
TH1F *histoNum = 0x0, *histoDen=0x0, *currHisto = 0x0;
}
-
//_____________________________________________________________________________
void AliMUONTriggerEfficiencyCells::CalculateEfficiency(Int_t trigger44, Int_t trigger34,
Float_t &efficiency, Float_t &error,
else error = TMath::Sqrt(efficiency*q/((Double_t)trigger34));
}
+
//_____________________________________________________________________________
-void AliMUONTriggerEfficiencyCells::DisplayEfficiency(Bool_t perSlat)
+void AliMUONTriggerEfficiencyCells::CheckFiredStrips(const Char_t* geoFilename,
+ const Char_t* cdbStorage,
+ Int_t runNumber)
{
//
- /// Display calculated efficiency.
+ /// Check for fired strips participating to efficiency
+ /// calculation (when available).
+ /// Strips inside a local board should be quite homogeneously hit
+ /// If not, this could be a problem of electronics (i.e. ADULT board off).
//
- const Int_t kNumOfBoards = AliMpConstants::NofLocalBoards();
+ if(!fFiredStrips) {
+ AliWarning("List of fired pads not present. Check not performable.");
+ return;
+ }
+
+ GetListsForCheck(geoFilename, cdbStorage, runNumber);
+
+ Char_t histoName[40], histoTitle[90];
+
+ // Check fired pads (when available)
+ if(fFiredFitHistoList){
+ TH1F *histo1D = 0x0;
+ TF1 *fitFunc = 0x0;
+ TCanvas *histoFiredCan[20];
+ Int_t nEntries = fFiredFitHistoList->GetEntries();
+ for(Int_t iEntry=0; iEntry<nEntries; iEntry++){
+ histo1D = (TH1F*)fFiredFitHistoList->At(iEntry);
+ printf("Problems found in %s\n", histo1D->GetTitle());
+ }
+ Int_t nPrintCan = nEntries;
+ if(nPrintCan>20) {
+ AliWarning("Too many boards with problems: only 20 will be shown");
+ nPrintCan = 20;
+ }
+ for(Int_t iCan=0; iCan<nPrintCan; iCan++){
+ histo1D = (TH1F*)fFiredFitHistoList->At(iCan);
+ histoFiredCan[iCan] = new TCanvas(histoName, histoTitle, 100+10*iCan, 10*iCan, 700, 700);
+ histoFiredCan[iCan]->SetRightMargin(0.14);
+ histoFiredCan[iCan]->SetLeftMargin(0.12);
+ histo1D->Draw("E");
+ fitFunc = histo1D->GetFunction("pol0");
+ fitFunc->SetLineColor(2);
+ fitFunc->Draw("same");
+ }
+ if(nEntries==0){
+ printf("\nAll local boards seem ok!!\n\n");
+ }
+ }
+}
+
+
+//_____________________________________________________________________________
+void AliMUONTriggerEfficiencyCells::DisplayEfficiency(Bool_t perSlat,
+ const Char_t* geoFilename,
+ const Char_t* cdbStorage,
+ Int_t runNumber)
+{
+ //
+ /// Display calculated efficiency.
+ //
Bool_t isInitSlat = kFALSE, isInitBoard = kFALSE;
for(Int_t chCath=0; chCath<fgkNplanes; chCath++){
printf("Trigger efficiency not initialized for slat.\nPlease try option kFALSE.\n");
return;
}
+
+ GetListsForCheck(geoFilename, cdbStorage, runNumber);
+
+ const Int_t kNumOfBoards = AliMpConstants::NofLocalBoards();
+
+ TH2F *histo = 0x0;
+ Char_t histoName[40], histoTitle[90];
+ TPaveLabel *boardLabel = 0x0;
+
+ // Plot fired strips (when available)
+ if(fFiredDisplayHistoList){
+ TCanvas *displayFiredCan[fgkNplanes];
+ for(Int_t chCath=0; chCath<fgkNplanes; chCath++){
+ histo = (TH2F*)fFiredDisplayHistoList->At(chCath);
+ sprintf(histoName, "%sCan", histo->GetName());
+ sprintf(histoTitle, "%s", histo->GetTitle());
+ displayFiredCan[chCath] = new TCanvas(histoName, histoTitle, 100+10*chCath, 10*chCath, 700, 700);
+ displayFiredCan[chCath]->SetRightMargin(0.14);
+ displayFiredCan[chCath]->SetLeftMargin(0.12);
+ histo->GetYaxis()->SetTitleOffset(1.4);
+ histo->SetStats(kFALSE);
+ histo->Draw("COLZ");
+ for (Int_t board = 0; board < kNumOfBoards; board++) {
+ Int_t currLabel = chCath * kNumOfBoards + board;
+ boardLabel = (TPaveLabel*)fBoardLabelList->At(currLabel);
+ boardLabel->Draw("same");
+ }
+ }
+ }
+
+ // Plot efficiency
+ if(fDisplayHistoList){
+ TCanvas *can[fgkNplanes];
+ for(Int_t chCath=0; chCath<fgkNplanes; chCath++){
+ Int_t currChCath = chCath;
+ if(perSlat==kTRUE) currChCath += fgkNplanes;
+ histo = (TH2F*)fDisplayHistoList->At(currChCath);
+ sprintf(histoName, "%sCan", histo->GetName());
+ sprintf(histoTitle, "%s", histo->GetTitle());
+ can[chCath] = new TCanvas(histoName, histoTitle, 100+10*chCath, 10*chCath, 700, 700);
+ can[chCath]->SetRightMargin(0.14);
+ can[chCath]->SetLeftMargin(0.12);
+ histo->GetZaxis()->SetRangeUser(0.,1.);
+ histo->GetYaxis()->SetTitleOffset(1.4);
+ histo->SetStats(kFALSE);
+ histo->Draw("COLZ");
+ for (Int_t board = 0; board < kNumOfBoards; board++) {
+ Int_t currLabel = chCath * kNumOfBoards + board;
+ if(perSlat==kTRUE) currLabel += kNumOfBoards * fgkNplanes;
+ boardLabel = (TPaveLabel*)fBoardLabelList->At(currLabel);
+ boardLabel->Draw("same");
+ }
+ }
+ }
+}
+
+
+//__________________________________________________________________________
+Bool_t AliMUONTriggerEfficiencyCells::GetListsForCheck(const Char_t* geoFilename,
+ const Char_t* cdbStorage,
+ Int_t runNumber)
+{
+ //
+ /// Getting histograms for efficiency,
+ /// map of fired strips entering efficiency calculations,
+ /// fits for checking switched-off elements in chambers.
+ //
+
+ const Int_t kNumOfBoards = AliMpConstants::NofLocalBoards();
+ const Float_t kChi2RedMax = 1.5;
+ const Float_t kDummyFired = 1e-5;
+
+ if(fDisplayHistoList || fBoardLabelList || fFiredFitHistoList || fFiredDisplayHistoList) return kTRUE;
+
+ if(!fDisplayHistoList) fDisplayHistoList = new TList();
+ if(!fBoardLabelList) fBoardLabelList = new TList();
+ if(!fFiredFitHistoList && fFiredStrips) fFiredFitHistoList = new TList();
+ if(!fFiredDisplayHistoList && fFiredStrips) fFiredDisplayHistoList = new TList();
+
+ AliCDBManager::Instance()->SetDefaultStorage(cdbStorage);
+ AliCDBManager::Instance()->SetRun(runNumber);
+
+ AliMUONGeometryTransformer *transform = new AliMUONGeometryTransformer();
+ transform->LoadGeometryData(geoFilename);
+
+ if(!AliMpCDB::LoadDDLStore()) // Load DDL store from OCDB
+ AliFatal("Could not access mapping OCDB");
- Int_t side, col, line, nbx, slat;
- Float_t xCenter, yCenter, zCenter, xWidth, yWidth;
- Float_t x1Label, x2Label, y1Label, y2Label;
- Int_t x1, y1, x2, y2, board=0;
- Char_t name[8], text[200];
+ AliMpDDLStore *ddlStore = AliMpDDLStore::Instance();
+ Int_t line, slat;
+ Float_t xLocal1=0., yLocal1=0., xLocal2=0., yLocal2=0.;
+ Float_t xg1, yg1, zg1, xg2, yg2, zg2;
+ Float_t xWidth=0., yWidth=0.;
+ Float_t x1Label=0., x2Label=0., y1Label=0., y2Label=0.;
+ Float_t x1LabelSlat=0., x2LabelSlat=0., y1LabelSlat=0., y2LabelSlat=0.;
+ Int_t x1=0, y1=0, x2=0, y2=0, localId;
gStyle->SetPalette(1);
Char_t *cathCode[fgkNcathodes] = {"bendPlane", "nonBendPlane"};
- Float_t boardsX = 257.00; // cm
- Float_t boardsY = 307.00; // cm
+ Float_t boardsX = 280.00; // cm
+ Float_t boardsY = 335.00; // cm
+
+ // Check fired pads (when available)
+ Int_t maxY[fgkNcathodes] = {64,1};
+ TH3F *padFired[fgkNplanes];
+ TH1F *histoFired[fgkNplanes][234];
+ TH2F *histoFiredDisplay[fgkNplanes];
+ TF1 *fitFunc = 0x0;
+ Bool_t isStripOffInBoard[fgkNplanes][234];
TH2F *histo[fgkNplanes];
+ TH2F *histoSlat[fgkNplanes];
TPaveLabel *boardLabel[fgkNplanes][234];
- assert(kNumOfBoards==234);
- TArrayI boardsPerColumn[9];
- for(Int_t iLine=0; iLine<9; iLine++){
- boardsPerColumn[iLine].Set(7);
- boardsPerColumn[iLine].Reset();
- }
+ TPaveLabel *boardLabelSlat[fgkNplanes][234];
+ assert(kNumOfBoards==234);
- Char_t histoName[40], histoTitle[90], labelTxt[5];
+ Char_t histoName[40], histoTitle[90], labelTxt[5], labelSlatTxt[5];
- Float_t efficiency, efficiencyError;
+ Float_t efficiency, efficiencyError, efficiencySlat, efficiencySlatError;
+ // Book histos
for(Int_t cath=0; cath<fgkNcathodes; cath++){
+ padFired[cath] = 0x0;
+ if(fFiredStrips) padFired[cath] = (TH3F*)fFiredStrips->At(cath);
for(Int_t ch=0; ch<fgkNchambers; ch++){
Int_t chCath = fgkNchambers*cath + ch;
sprintf(histoName, "%sChamber%i", cathCode[cath], 11+ch);
histo[chCath] = new TH2F(histoName, histoTitle, (Int_t)boardsX, -boardsX, boardsX, (Int_t)boardsY, -boardsY, boardsY);
histo[chCath]->SetXTitle("X (cm)");
histo[chCath]->SetYTitle("Y (cm)");
- }
- }
-
- TString mapspath = gSystem->Getenv("ALICE_ROOT");
- mapspath.Append("/MUON/data");
+ fDisplayHistoList->Add(histo[chCath]);
+
+ if(!fFiredStrips) continue;
+ sprintf(histoName, "firedPads%sChamber%i", cathCode[cath], 11+ch);
+ sprintf(histoTitle, "Chamber %i: Fired pads %s", 11+ch, cathCode[cath]);
+ histoFiredDisplay[chCath] = new TH2F(histoName, histoTitle, (Int_t)boardsX, -boardsX, boardsX, (Int_t)boardsY, -boardsY, boardsY);
+ histoFiredDisplay[chCath]->SetXTitle("X (cm)");
+ histoFiredDisplay[chCath]->SetYTitle("Y (cm)");
+ fFiredDisplayHistoList->Add(histoFiredDisplay[chCath]);
+
+ for(Int_t ib=0; ib<kNumOfBoards; ib++){
+ sprintf(histoName, "%sChamber%iBoard%i", cathCode[cath], 11+ch, ib+1);
+ sprintf(histoTitle, "Chamber %i: fired pads %s board = %i", 11+ch, cathCode[cath], ib+1);
+ histoFired[chCath][ib] = new TH1F(histoName, histoTitle, 16, -0.5, 15.5);
+ histoFired[chCath][ib]->SetXTitle("board");
+ isStripOffInBoard[chCath][ib] = kFALSE;
+ } // loop on board
+ } // loop on chamber
+ } // loop on cathode
- sprintf(text,"%s/guimapp11.txt",mapspath.Data());
-
- FileStat_t fs;
- if(gSystem->GetPathInfo(text,fs)){
- AliWarning(Form("Map file %s not found. Nothing done",text));
- return;
- }
-
- FILE *fmap = fopen(text,"r");
-
- for (Int_t ib = 0; ib < kNumOfBoards; ib++) {
- fscanf(fmap,"%d %d %d %d %f %f %f %f %f %s \n",&side,&col,&line,&nbx,&xCenter,&yCenter,&xWidth,&yWidth,&zCenter,&name[0]);
- if(side==0) continue;
- boardsPerColumn[line-1][col-1]++;
+ for(Int_t cath=0; cath<fgkNcathodes; cath++){
+ for(Int_t ch=0; ch<fgkNchambers; ch++){
+ Int_t chCath = fgkNchambers*cath + ch;
+ sprintf(histoName, "%sChamber%iSlatEff", cathCode[cath], 11+ch);
+ sprintf(histoTitle, "Chamber %i: efficiency %s per slat", 11+ch, cathCode[cath]);
+ histoSlat[chCath] = new TH2F(histoName, histoTitle, (Int_t)boardsX, -boardsX, boardsX, (Int_t)boardsY, -boardsY, boardsY);
+ histoSlat[chCath]->SetXTitle("X (cm)");
+ histoSlat[chCath]->SetYTitle("Y (cm)");
+ fDisplayHistoList->Add(histoSlat[chCath]);
+ }
}
- rewind(fmap);
+ // loop over the trigger DDL (Right: 20, Left: 21)
+ for (Int_t iDDL = 20; iDDL <= 21; ++iDDL) {
+ AliMpDDL* ddl = ddlStore->GetDDL(iDDL);
+ Int_t nCrate = ddl->GetNofTriggerCrates();
+ // loop over the number of crates in DDL
+ for (Int_t index = 0; index < nCrate; ++index) {
+ // get crate object
+ AliMpTriggerCrate* crate = ddlStore->GetTriggerCrate(iDDL, index);
+ Int_t nLocal = crate->GetNofLocalBoards();
+ for (Int_t iLocal = 0; iLocal < nLocal; ++iLocal) {
+ // get local board Id from crate object
+ localId = crate->GetLocalBoardId(iLocal);
+ // get local board object
+ AliMpLocalBoard* localBoard = ddlStore->GetLocalBoard(localId);
+
+ if (!localBoard->IsNotified()) continue;
+
+ // get detection element connected to this board
+ for (Int_t ch = 0; ch < fgkNchambers; ++ch) {
+ Int_t iCh = ch + AliMpConstants::NofTrackingChambers();
+ Int_t detElemId = ddlStore->GetDEfromLocalBoard(localId, iCh);
+
+ if (!detElemId) continue;
+
+ // get segmentation
+ for (Int_t cath = 0; cath < fgkNcathodes; ++cath) {
+ const AliMpVSegmentation* seg = AliMpSegmentation::Instance()
+ ->GetMpSegmentation(detElemId,
+ AliMp::GetCathodType(cath));
- for (Int_t ib = 0; ib < kNumOfBoards; ib++) {
- fscanf(fmap,"%d %d %d %d %f %f %f %f %f %s \n",&side,&col,&line,&nbx,&xCenter,&yCenter,&xWidth,&yWidth,&zCenter,&name[0]);
-
- board=0;
- for(Int_t iCol=1; iCol<=col; iCol++){
- Int_t lastLine = 9;
- if(iCol==col) lastLine = line-1;
- for(Int_t iLine=1; iLine<=lastLine; iLine++){
- board += boardsPerColumn[iLine-1][iCol-1];
- }
- }
- if(side==0) board += kNumOfBoards/2;
- board += nbx - 1;
-
- slat = (line+13)%fgkNslats;
- if(side==0) slat = 14-line;
+ Int_t chCath = fgkNchambers*cath + ch;
+ slat = detElemId%100;
+ Int_t nStrips=0;
+
+ // loop over strips
+ for (Int_t ibitxy = 0; ibitxy < 16; ++ibitxy) {
+ // get pad from electronics
+ AliMpPad pad = seg->PadByLocation(AliMpIntPair(localId,
+ ibitxy),kFALSE);
+
+ if (!pad.IsValid()) continue;
+ if(cath==0){ // Geometry info from bending plane only
+ if(ibitxy==0) {
+ xLocal1 = pad.Position().X();
+ yLocal1 = pad.Position().Y();
+ xWidth = pad.Dimensions().X();
+ yWidth = pad.Dimensions().Y();
+ }
+ xLocal2 = pad.Position().X();
+ yLocal2 = pad.Position().Y();
+ }
+
+ // Check fired pads (when available)
+ if(padFired[cath]) {
+ Int_t padX = pad.GetIndices().GetFirst();
+ Int_t padY = pad.GetIndices().GetSecond();
+ Int_t currPair = padX*maxY[cath] + padY;
+ nStrips++;
+ //printf("cath %i board = %i (%2i, %2i) -> %i\n",
+ //cath, localId, padX, padY, currPair);
+ Int_t chBin = padFired[cath]->GetXaxis()->FindBin(ch);
+ Int_t slatBin = padFired[cath]->GetYaxis()->FindBin(slat);
+ Int_t pairBin = padFired[cath]->GetZaxis()->FindBin(currPair);
+ Float_t nFired = padFired[cath]->GetBinContent(chBin, slatBin, pairBin);;
+
+ Float_t dimX = pad.Dimensions().X();
+ Float_t dimY = pad.Dimensions().Y();
+
+ Float_t stripX1 = pad.Position().X();
+ Float_t stripY1 = pad.Position().Y();
+ Float_t stripX2 = pad.Position().X();
+ Float_t stripY2 = pad.Position().Y();
+
+ transform->Local2Global(detElemId, stripX1, stripY1, 0, xg1, yg1, zg1);
+ transform->Local2Global(detElemId, stripX2, stripY2, 0, xg2, yg2, zg2);
+
+ Float_t x1Float = TMath::Min(xg1,xg2) - dimX;
+ Float_t y1Float = TMath::Min(yg1,yg2) - dimY;
+ Float_t x2Float = TMath::Max(xg1,xg2) + dimX;
+ Float_t y2Float = TMath::Max(yg1,yg2) + dimY;
+
+ Int_t x1Int = histoFiredDisplay[chCath]->GetXaxis()->FindBin(x1Float)+1;
+ Int_t y1Int = histoFiredDisplay[chCath]->GetYaxis()->FindBin(y1Float)+1;
+ Int_t x2Int = histoFiredDisplay[chCath]->GetXaxis()->FindBin(x2Float)-1;
+ Int_t y2Int = histoFiredDisplay[chCath]->GetYaxis()->FindBin(y2Float)-1;
+
+ for(Int_t binX=x1Int; binX<=x2Int; binX++){
+ for(Int_t binY=y1Int; binY<=y2Int; binY++){
+ histoFiredDisplay[chCath]->SetBinContent(binX, binY, nFired);
+ }
+ }
- for(Int_t chCath=0; chCath<fgkNplanes; chCath++){
- x1 = histo[chCath]->GetXaxis()->FindBin(xCenter-xWidth/2.)+1;
- y1 = histo[chCath]->GetYaxis()->FindBin(yCenter-yWidth/2.)+1;
- x2 = histo[chCath]->GetXaxis()->FindBin(xCenter+xWidth/2.)-1;
- y2 = histo[chCath]->GetYaxis()->FindBin(yCenter+yWidth/2.)-1;
-
- x1Label = xCenter-xWidth/2.;
- y1Label = yCenter-yWidth/2.;
- x2Label = xCenter+xWidth/2.;
- y2Label = yCenter+yWidth/2.;
- sprintf(labelTxt,"%3d", board+1);
- if(perSlat){
- x1Label = 140.;
- x2Label = x1Label + 40.;
- y1Label = -285. + ((Float_t)(line - 1)) * 68;
- y2Label = y1Label + 34.;
- if(side==0){
- x1Label = -x2Label;
- x2Label = x1Label + 40.;
- }
- sprintf(labelTxt,"%2d", slat);
- }
-
- boardLabel[chCath][board] = new TPaveLabel(x1Label, y1Label, x2Label, y2Label, labelTxt);
- boardLabel[chCath][board]->SetFillStyle(0);
- boardLabel[chCath][board]->SetBorderSize(0);
+ if(nFired==0.) nFired = kDummyFired;
+ histoFired[chCath][localId-1]->Fill(ibitxy, nFired);
+ }
+ } // loop on strips
+
+ if(cath==0){ // Geometry info from bending plane only
+ transform->Local2Global(detElemId, xLocal1, yLocal1, 0, xg1, yg1, zg1);
+ transform->Local2Global(detElemId, xLocal2, yLocal2, 0, xg2, yg2, zg2);
+
+ // Per board
+ x1Label = TMath::Min(xg1,xg2) - xWidth;
+ y1Label = TMath::Min(yg1,yg2) - yWidth;
+ x2Label = TMath::Max(xg1,xg2) + xWidth;
+ y2Label = TMath::Max(yg1,yg2) + yWidth;
+
+ x1 = histo[ch]->GetXaxis()->FindBin(x1Label)+1;
+ y1 = histo[ch]->GetYaxis()->FindBin(y1Label)+1;
+ x2 = histo[ch]->GetXaxis()->FindBin(x2Label)-1;
+ y2 = histo[ch]->GetYaxis()->FindBin(y2Label)-1;
+
+ sprintf(labelTxt,"%3d", localId);
+
+ // Per slat
+ line = localBoard->GetPosition().GetFirst();
+ x1LabelSlat = 140.;
+ x2LabelSlat = x1LabelSlat + 40.;
+ y1LabelSlat = -285. + ((Float_t)(line - 1)) * 68;
+ y2LabelSlat = y1LabelSlat + 34.;
+ if(localId>kNumOfBoards/2){
+ x1LabelSlat = -x2LabelSlat;
+ x2LabelSlat = x1LabelSlat + 40.;
+ }
+ sprintf(labelSlatTxt,"%2d", slat);
+ }
- Int_t histoBin = board+1;
- if(perSlat) histoBin = slat+1;
+ boardLabel[chCath][localId-1] = new TPaveLabel(x1Label, y1Label, x2Label, y2Label, labelTxt);
+ boardLabel[chCath][localId-1]->SetFillStyle(0);
+ boardLabel[chCath][localId-1]->SetBorderSize(0);
+
+ boardLabelSlat[chCath][localId-1] = new TPaveLabel(x1LabelSlat, y1LabelSlat, x2LabelSlat, y2LabelSlat, labelSlatTxt);
+ boardLabelSlat[chCath][localId-1]->SetFillStyle(0);
+ boardLabelSlat[chCath][localId-1]->SetBorderSize(0);
+
+ Int_t histoBin = localId;
+ efficiency = fBoardEfficiency[chCath]->GetBinContent(histoBin);
+ efficiencyError = fBoardEfficiency[chCath]->GetBinError(histoBin);
+ if(efficiency==0.) efficiency = kDummyFired; // It allows to graphically distinguish
+ // efficiency = 0 from no hit on board
+
+ histoBin = slat+1;
+ efficiencySlat = fSlatEfficiency[chCath]->GetBinContent(histoBin);
+ efficiencySlatError = fSlatEfficiency[chCath]->GetBinError(histoBin);
+ if(efficiencySlat==0.) efficiencySlat = kDummyFired; // It allows to graphically distinguish
+ // efficiency = 0 from no hit on slat
+
+ for(Int_t binX=x1; binX<=x2; binX++){
+ for(Int_t binY=y1; binY<=y2; binY++){
+ histo[chCath]->SetBinContent(binX, binY, efficiency);
+ histo[chCath]->SetBinError(binX, binY, efficiencyError);
+ histoSlat[chCath]->SetBinContent(binX, binY, efficiencySlat);
+ histoSlat[chCath]->SetBinError(binX, binY, efficiencySlatError);
+ }
+ }
- if(!perSlat){
- efficiency = fBoardEfficiency[chCath]->GetBinContent(histoBin);
- efficiencyError = fBoardEfficiency[chCath]->GetBinError(histoBin);
- }
- else {
- efficiency = fSlatEfficiency[chCath]->GetBinContent(histoBin);
- efficiencyError = fSlatEfficiency[chCath]->GetBinError(histoBin);
- }
+ // Check fired pads (when available)
+ if(padFired[cath]) {
+ histoFired[chCath][localId-1]->Fit("pol0","Q0R","",0., (Float_t)nStrips-1.);
+ fitFunc = histoFired[chCath][localId-1]->GetFunction("pol0");
+ Float_t chi2 = fitFunc->GetChisquare();
+ Float_t ndf = (Float_t)fitFunc->GetNDF();
+ Float_t reducedChi2 = chi2/ndf;
+ if(reducedChi2>kChi2RedMax) {
+ isStripOffInBoard[chCath][localId-1] = kTRUE;
+ //printf("Chamber = %i Cath = %i Board %i: chi2/NDF = %f\tparam = %f\n", ch, cath, localId, reducedChi2, fitFunc->GetParameter(0));
+ fFiredFitHistoList->Add(histoFired[chCath][localId-1]);
+ }
+ }
+ } // loop on cathodes
+ } // loop on chambers
+ } // loop on boards
+ } // loop on crates
+ } // loop on DDL
- for(Int_t binX=x1; binX<=x2; binX++){
- for(Int_t binY=y1; binY<=y2; binY++){
- histo[chCath]->SetBinContent(binX, binY, efficiency);
- histo[chCath]->SetBinError(binX, binY, efficiencyError);
- }
- }
+ for(Int_t chCath=0; chCath<fgkNplanes; chCath++){
+ for(Int_t ib=0; ib<kNumOfBoards; ib++){
+ fBoardLabelList->Add(boardLabel[chCath][ib]);
}
}
-
- fclose(fmap);
-
- TCanvas *can[8];
for(Int_t chCath=0; chCath<fgkNplanes; chCath++){
- sprintf(histoName, "%sCan", histo[chCath]->GetName());
- sprintf(histoTitle, "%s", histo[chCath]->GetTitle());
- can[chCath] = new TCanvas(histoName, histoTitle, 100+10*chCath, 10*chCath, 700, 700);
- can[chCath]->SetRightMargin(0.14);
- can[chCath]->SetLeftMargin(0.12);
- histo[chCath]->GetZaxis()->SetRangeUser(0.,1.);
- histo[chCath]->GetYaxis()->SetTitleOffset(1.4);
- histo[chCath]->SetStats(kFALSE);
- histo[chCath]->Draw("COLZ");
- for (Int_t board = 0; board < kNumOfBoards; board++) {
- boardLabel[chCath][board]->Draw("same");
+ for(Int_t ib=0; ib<kNumOfBoards; ib++){
+ fBoardLabelList->Add(boardLabelSlat[chCath][ib]);
}
}
+
+ return kTRUE;
}
}
FillHistosFromList();
+
+ if(!fFiredStrips) {
+ AliWarning("Histograms for fired region check not implemented in object");
+ return kFALSE;
+ }
+ if(!other.fFiredStrips) {
+ AliWarning("Histograms for fired region check not implemented in object passed as argument");
+ return kFALSE;
+ }
+
+ TH3F *currFired = 0x0, *otherFired = 0x0;
+ nentries = fFiredStrips->GetEntries();
+ for(Int_t iEntry=0; iEntry<nentries; iEntry++){
+ currFired = (TH3F*)fFiredStrips->At(iEntry);
+ otherFired = (TH3F*)fFiredStrips->At(iEntry);
+ currFired->Add(otherFired);
+ }
+
return kTRUE;
}