In trigger:
authorivana <ivana@f7af4fe6-9843-0410-8265-dc069ae4e863>
Tue, 3 Nov 2009 10:42:26 +0000 (10:42 +0000)
committerivana <ivana@f7af4fe6-9843-0410-8265-dc069ae4e863>
Tue, 3 Nov 2009 10:42:26 +0000 (10:42 +0000)
- Add possibility to fill the trigger board efficiency with the average
  value of the corresponding RPC.
- Update MUONTriggerChamberEfficiency.C
- Update READMEtrigger.txt
(Diego)

MUON/AliMUONTriggerEfficiencyCells.cxx
MUON/AliMUONTriggerEfficiencyCells.h
MUON/MUONTriggerChamberEfficiency.C
MUON/READMEtrigger.txt

index d1a5f99..652449d 100755 (executable)
@@ -21,6 +21,7 @@
 // Classes for display
 #include "AliMUONTriggerDisplay.h"
 #include "AliCDBManager.h"
+#include "AliMpCDB.h"
 #include "AliMpDDLStore.h"
 #include "AliMpLocalBoard.h"
 #include "AliMpPad.h"
@@ -106,9 +107,12 @@ fFiredDisplayHistoList(0x0)
 ///  Constructor using an ASCII file.
   CheckConstants();
   Reset();
+  InitHistos();
   ReadFile(filename);
 }
 
+
+//__________________________________________________________________________
 AliMUONTriggerEfficiencyCells::AliMUONTriggerEfficiencyCells(TList *countHistoList,
                                                             TList *noCountHistoList)
 :
@@ -121,7 +125,7 @@ fBoardLabelList(0x0),
 fFiredFitHistoList(0x0),
 fFiredDisplayHistoList(0x0)
 {
-///  Constructor using an ASCII file.
+///  Constructor using a list of histograms with counts.
   CheckConstants();
   Reset();
   InitHistos();
@@ -177,6 +181,10 @@ AliMUONTriggerEfficiencyCells& AliMUONTriggerEfficiencyCells::operator=(const Al
 AliMUONTriggerEfficiencyCells::~AliMUONTriggerEfficiencyCells()
 {
 ///  Destructor.
+
+  delete [] fBoardEfficiency;
+  delete [] fSlatEfficiency;
+  Reset();
 }
 
 
@@ -219,7 +227,6 @@ void AliMUONTriggerEfficiencyCells::ReadFile(const Char_t* filename)
       return;
   }
 
-  InitHistos();
   ifstream file(fileName.Data());
   Char_t dat[50];
   if (file.good()){
@@ -228,7 +235,7 @@ void AliMUONTriggerEfficiencyCells::ReadFile(const Char_t* filename)
       else AliWarning("File .dat in wrong format");
       file.close();
   } else {
-      AliWarning(Form("Can't read file %s",fileName.Data()));
+    AliError(Form("Can't read file %s",fileName.Data()));
   }
 }
 
@@ -245,7 +252,7 @@ void AliMUONTriggerEfficiencyCells::ReadFileBoards(ifstream &file)
            file >> detEl;
            chamber = FindChamberIndex(detEl);
            for(Int_t cath=0; cath<fgkNcathodes; cath++){
-               chCath = fgkNchambers*cath + chamber;
+               chCath = GetPlane(chamber, cath);
                file >> dat;
                file >> datInt;
                for(Int_t board=1; board<=AliMpConstants::NofLocalBoards(); board++){
@@ -263,10 +270,11 @@ void AliMUONTriggerEfficiencyCells::ReadHistoBoards(const Char_t *filename)
 {
 ///  Structure of file (.root) containing local board efficency
     TFile *file = new TFile(filename, "read");
-    if(!file) {
-       AliWarning(Form("Can't read file %s",filename));
-       return;
+    if(!file || !file->IsOpen()) {
+      AliError(Form("Can't read file %s",filename));
+      return;
     }
+
     TString histoName;
     TString cathCode[fgkNcathodes] = {"bendPlane", "nonBendPlane"};
     enum {kAllChEff, kChNonEff, kNumOfHistoTypes};
@@ -309,7 +317,6 @@ void AliMUONTriggerEfficiencyCells::ReadHistoBoards(const Char_t *filename)
       }
     }
 
-    InitHistos();
     FillHistosFromList();
 }
 
@@ -339,14 +346,36 @@ Int_t AliMUONTriggerEfficiencyCells::FindChamberIndex(Int_t detElemId) const
 
 //__________________________________________________________________________
 void
-AliMUONTriggerEfficiencyCells::Reset()
+AliMUONTriggerEfficiencyCells::Reset(Bool_t resetAll)
 {
 ///  Sets our internal array contents to zero.
 
-  for(Int_t chCath=0; chCath<fgkNplanes; chCath++){
-    fBoardEfficiency[chCath] = 0x0;
-    fSlatEfficiency[chCath] = 0x0;
+  if ( resetAll ){
+    for(Int_t chCath=0; chCath<fgkNplanes; chCath++){
+      if ( fBoardEfficiency[chCath] ) {
+       fBoardEfficiency[chCath] = 0x0;
+      }
+      if ( fSlatEfficiency[chCath] ) {
+       fSlatEfficiency[chCath] = 0x0;
+      }
+    }
+
+    if ( fCountHistoList )
+      delete fCountHistoList;
+    if ( fNoCountHistoList )
+      delete fNoCountHistoList;
   }
+
+  if ( fFiredStrips )
+    delete fFiredStrips;
+  if ( fDisplayHistoList )
+    delete fDisplayHistoList;
+  if ( fBoardLabelList )
+    delete fBoardLabelList;
+  if ( fFiredFitHistoList )
+    delete fFiredFitHistoList;
+  if ( fFiredDisplayHistoList )
+    delete fFiredDisplayHistoList;
 }
 
 
@@ -365,7 +394,7 @@ AliMUONTriggerEfficiencyCells::InitHistos()
 
   for(Int_t ch=0; ch<fgkNchambers; ch++){
     for(Int_t cath=0; cath<fgkNcathodes; cath++){
-      chCath = fgkNchambers*cath + ch;
+      chCath = GetPlane(ch, cath);
       histoName = Form("%sBoardEffChamber%i", cathCode[cath].Data(), 11+ch);
       fBoardEfficiency[chCath] = new TH1F(histoName.Data(), histoName.Data(), kNumOfBoards, 1-0.5, kNumOfBoards+1.-0.5);
       histoName = Form("%sSlatEffChamber%i", cathCode[cath].Data(), 11+ch);
@@ -377,18 +406,21 @@ AliMUONTriggerEfficiencyCells::InitHistos()
 
 //__________________________________________________________________________
 void
-AliMUONTriggerEfficiencyCells::FillHistosFromList()
+AliMUONTriggerEfficiencyCells::FillHistosFromList(Bool_t useMeanValues)
 {
 ///  Fills internal histos from list.
 
   Int_t nHistoBins=0;
   TH1F *histoNum = 0x0, *histoDen=0x0, *currHisto = 0x0;
-  TString slatName = "Slat", boardName = "Board", histoName;
-  Int_t iHistoBoard = -1, iHistoSlat = -1;
+  TString slatName = "SLAT", boardName = "BOARD", histoName;
   Float_t efficiency, efficiencyError;
+  TString chName;
 
   Int_t nentries = fCountHistoList->GetEntries();
 
+  if ( useMeanValues )
+    AliInfo("Boards filled with the average efficiency of the RPC");
+
   for(Int_t iEntry=0; iEntry<nentries; iEntry++){
     histoNum = (TH1F*)fCountHistoList->At(iEntry);
     histoDen = (TH1F*)fNoCountHistoList->At(iEntry);
@@ -403,17 +435,31 @@ AliMUONTriggerEfficiencyCells::FillHistosFromList()
     }
 
     histoName = histoNum->GetName();
+    histoName.ToUpper();
     nHistoBins = histoNum->GetNbinsX();
 
+
+    Int_t currCh = 0;
+    for(Int_t ich = 11; ich<=14; ich++){
+      chName = Form("%i", ich);
+      if ( histoName.Contains(chName.Data())){
+       currCh = ich-11;
+       break;
+      }
+    }
+    Int_t currCath = ( histoName.Contains("NONBEND") ) ? 1 : 0;
+
+    Int_t chCath = GetPlane(currCh, currCath);
+
     if(histoName.Contains(boardName)){
-      iHistoBoard++;
-      currHisto = fBoardEfficiency[iHistoBoard];
+      if ( useMeanValues ) continue;
+      currHisto = fBoardEfficiency[chCath];
     }
     else if(histoName.Contains(slatName)){
-      iHistoSlat++;
-      currHisto = fSlatEfficiency[iHistoSlat];
+      currHisto = fSlatEfficiency[chCath];
     }
-    else continue;
+    else
+      continue;
 
     for(Int_t iBin=1; iBin<=nHistoBins; iBin++){
       CalculateEfficiency((Int_t)histoNum->GetBinContent(iBin), (Int_t)histoNum->GetBinContent(iBin) + (Int_t)histoDen->GetBinContent(iBin), efficiency, efficiencyError, kFALSE);
@@ -421,6 +467,21 @@ AliMUONTriggerEfficiencyCells::FillHistosFromList()
       currHisto->SetBinContent(iBin, efficiency);
       currHisto->SetBinError(iBin, efficiencyError);
     }
+
+    if ( useMeanValues ){
+      currHisto = fBoardEfficiency[chCath];
+      Int_t nBinsBoards = currHisto->GetNbinsX();
+      Int_t currChamber = currCh + AliMpConstants::NofTrackingChambers();
+      for ( Int_t iBinBoard = 1; iBinBoard<= nBinsBoards; iBinBoard++){
+       Int_t detElemId = AliMpDDLStore::Instance()->GetDEfromLocalBoard(iBinBoard, currChamber);
+       Int_t iBin = histoNum->FindBin(detElemId%100);
+
+       CalculateEfficiency((Int_t)histoNum->GetBinContent(iBin), (Int_t)histoNum->GetBinContent(iBin) + (Int_t)histoDen->GetBinContent(iBin), efficiency, efficiencyError, kFALSE);
+
+       currHisto->SetBinContent(iBinBoard, efficiency);
+       currHisto->SetBinError(iBinBoard, efficiencyError);
+      }
+    }
   }
 }
 
@@ -447,8 +508,7 @@ void AliMUONTriggerEfficiencyCells::CalculateEfficiency(Int_t trigger44, Int_t t
 
 
 //_____________________________________________________________________________
-void AliMUONTriggerEfficiencyCells::CheckFiredStrips(const Char_t* cdbStorage,
-                                                    Int_t runNumber)
+void AliMUONTriggerEfficiencyCells::CheckFiredStrips()
 {
   //
   /// Check for fired strips participating to efficiency
@@ -462,7 +522,7 @@ void AliMUONTriggerEfficiencyCells::CheckFiredStrips(const Char_t* cdbStorage,
     return;
   }
 
-  GetListsForCheck(cdbStorage, runNumber);
+  GetListsForCheck();
 
   TString histoName;
 
@@ -501,9 +561,7 @@ void AliMUONTriggerEfficiencyCells::CheckFiredStrips(const Char_t* cdbStorage,
 
 
 //_____________________________________________________________________________
-void AliMUONTriggerEfficiencyCells::DisplayEfficiency(Bool_t perSlat,
-                                                     const Char_t* cdbStorage,
-                                                     Int_t runNumber)
+void AliMUONTriggerEfficiencyCells::DisplayEfficiency(Bool_t perSlat)
 {
   //
   /// Display calculated efficiency.
@@ -516,15 +574,23 @@ void AliMUONTriggerEfficiencyCells::DisplayEfficiency(Bool_t perSlat,
   }
 
   if(!isInitBoard){
-    printf("Trigger efficiency not initialized per board.\nDisplay not yet implemented.\n");
+    AliWarning("Trigger efficiency not initialized per board.\nDisplay not yet implemented.");
     return;
   }
   if(!isInitSlat && perSlat){
-    printf("Trigger efficiency not initialized for slat.\nPlease try option kFALSE.\n");
+    AliWarning("Trigger efficiency not initialized for slat.\nPlease try option kFALSE.");
+    return;
+  }
+  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;
   }
   
-  GetListsForCheck(cdbStorage, runNumber);
+  GetListsForCheck();
 
   //const Int_t kNumOfBoards = AliMpConstants::NofLocalBoards();
 
@@ -572,8 +638,7 @@ void AliMUONTriggerEfficiencyCells::DisplayEfficiency(Bool_t perSlat,
 
 
 //__________________________________________________________________________
-Bool_t AliMUONTriggerEfficiencyCells::GetListsForCheck(const Char_t* cdbStorage,
-                                                      Int_t runNumber)
+Bool_t AliMUONTriggerEfficiencyCells::GetListsForCheck()
 {
   //
   /// Getting histograms for efficiency, 
@@ -583,15 +648,18 @@ Bool_t AliMUONTriggerEfficiencyCells::GetListsForCheck(const Char_t* cdbStorage,
   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);
+  Reset(kFALSE);
+
+  fDisplayHistoList = new TList();
+  fDisplayHistoList->SetOwner();
+  fBoardLabelList = new TList();
+  fBoardLabelList->SetOwner();
+  if ( fFiredStrips ){
+    fFiredFitHistoList = new TList();
+    fFiredFitHistoList->SetOwner();
+    fFiredDisplayHistoList = new TList();
+    fFiredDisplayHistoList->SetOwner();
+  }
 
   TH3F* padFired = 0x0;
   TH2F* displayHisto = 0x0;
@@ -607,7 +675,7 @@ Bool_t AliMUONTriggerEfficiencyCells::GetListsForCheck(const Char_t* cdbStorage,
   for(Int_t iCath=0; iCath<fgkNcathodes; iCath++){
     if(fFiredStrips) padFired = (TH3F*)fFiredStrips->At(iCath);
     for(Int_t iCh=0; iCh<fgkNchambers; iCh++){
-      Int_t chCath = fgkNchambers*iCath + iCh;
+      Int_t chCath = GetPlane(iCh, iCath);
       Int_t currCh = 11 + iCh;
       histoName = Form("%sChamber%i", cathName[iCath].Data(), currCh);
       histoTitle = Form("Chamber %i: efficiency %s", currCh, cathName[iCath].Data());
@@ -646,7 +714,7 @@ Bool_t AliMUONTriggerEfficiencyCells::GetListsForCheck(const Char_t* cdbStorage,
 
   for(Int_t iCath=0; iCath<fgkNcathodes; iCath++){
     for(Int_t iCh=0; iCh<fgkNchambers; iCh++){
-      Int_t chCath = fgkNchambers*iCath + iCh;
+      Int_t chCath = GetPlane(iCh, iCath);
       Int_t currCh = 11+iCh;
       histoName = Form("%sChamber%iSlatEff", cathName[iCath].Data(), currCh);
       histoTitle = Form("Chamber %i: efficiency %s per slat", currCh, cathName[iCath].Data());
@@ -660,6 +728,8 @@ Bool_t AliMUONTriggerEfficiencyCells::GetListsForCheck(const Char_t* cdbStorage,
 
   if(!fFiredStrips) return kTRUE;
 
+  AliMpCDB::LoadDDLStore();
+
   // Check fired pads (when available)
   for(Int_t iLoc = 0; iLoc < AliMpConstants::NofLocalBoards(); iLoc++) {  
     Int_t iBoard = iLoc+1;
@@ -677,7 +747,7 @@ Bool_t AliMUONTriggerEfficiencyCells::GetListsForCheck(const Char_t* cdbStorage,
        continue;
       
       for(Int_t iCath=0; iCath<AliMpConstants::NofCathodes(); iCath++){
-       Int_t chCath = fgkNchambers*iCath + iCh;
+       Int_t chCath = GetPlane(iCh, iCath);
        // loop over strips
        const AliMpVSegmentation* seg = 
          AliMpSegmentation::Instance()->GetMpSegmentation(detElemId, AliMp::GetCathodType(iCath));
@@ -754,3 +824,37 @@ Bool_t AliMUONTriggerEfficiencyCells::SumRunEfficiency(const AliMUONTriggerEffic
     
   return kTRUE;
 }
+
+
+//__________________________________________________________________________
+Bool_t AliMUONTriggerEfficiencyCells::LowStatisticsSettings(Bool_t useMeanValues)
+{
+  //
+  /// In case of low statistics, fill the local board efficiency with
+  /// the average value of the RPC
+  //
+
+  if ( ! fCountHistoList ){
+    AliWarning("Cannot find histograms for RPC efficiency calculation. Nothing done!");
+    return kFALSE;
+  }
+  if ( !AliCDBManager::Instance()->GetDefaultStorage() ){
+    AliWarning("Please set default CDB storage (needed for mapping).");
+    return kFALSE;
+  }
+  if ( AliCDBManager::Instance()->GetRun() < 0 ){
+    AliWarning("Please set CDB run number (needed for mapping).");
+    return kFALSE;
+  }
+
+  AliMpCDB::LoadDDLStore();
+  
+  for(Int_t chCath=0; chCath<fgkNplanes; chCath++){
+    fBoardEfficiency[chCath]->Reset();
+    fSlatEfficiency[chCath]->Reset();
+  }
+
+  FillHistosFromList(useMeanValues);
+
+  return kTRUE;
+}
index d5c929e..bb48381 100755 (executable)
@@ -35,21 +35,18 @@ public:
 
   Bool_t SumRunEfficiency(const AliMUONTriggerEfficiencyCells &other);
 
+  Bool_t LowStatisticsSettings(Bool_t useMeanValues=kTRUE);
 
   // Methods for display
-  void DisplayEfficiency(Bool_t perSlat=kFALSE,
-                        const Char_t* cdbStorage = "local://$ALICE_ROOT/OCDB",
-                        Int_t runNumber=0);
+  void DisplayEfficiency(Bool_t perSlat=kFALSE);
 
   // Methods for efficiency check
   /// Set the list of fired strips
   void SetFiredStrips(TList *firedStrips){fFiredStrips = firedStrips;}
-  void CheckFiredStrips(const Char_t* cdbStorage = "local://$ALICE_ROOT/OCDB",
-                       Int_t runNumber=0);
-                                // Check for strips with lower counts than others:
-                                // syntomatic of possible read-out problems in boards
+  void CheckFiredStrips(); // Check for strips with lower counts than others:
+                           // syntomatic of possible read-out problems in boards
 protected:
-    void Reset();
+    void Reset(Bool_t resetAll = kTRUE);
     void ReadFile(const Char_t* filename="$ALICE_ROOT/MUON/data/efficiencyCells.dat");
     void CalculateEfficiency(Int_t trigger44, Int_t trigger34,
                             Float_t &efficiency, Float_t &error,
@@ -62,13 +59,15 @@ private:
     void ReadFileBoards(ifstream &file);
     void ReadHistoBoards(const Char_t* filename="MUON.TriggerEfficiencyMap.root");
     void InitHistos();
-    void FillHistosFromList();
-    Bool_t GetListsForCheck(const Char_t* cdbStorage, Int_t runNumber);
-    
+    void FillHistosFromList(Bool_t useMeanValues = kFALSE);
+    Bool_t GetListsForCheck();
+
     static const Int_t fgkNcathodes=2; ///<Number of cathodes
     static const Int_t fgkNchambers=4; ///<Number of chambers
     static const Int_t fgkNplanes=8;   ///<Number of planes
 
+    /// Return the current plane
+    Int_t GetPlane(Int_t chamber, Int_t cathode){ return fgkNchambers * cathode + chamber; }
     
     TH1F *fBoardEfficiency[fgkNplanes];///< the boards content
     TH1F *fSlatEfficiency[fgkNplanes];///< the slats content
index cf0086a..d5b65f3 100644 (file)
 
 #if !defined(__CINT__) || defined(__MAKECINT__)
 // ROOT includes
-#include "TFile.h"
-#include "TTree.h"
+#include "TGrid.h"
+#include "TString.h"
 
 // MUON includes
 #include "AliMUONCDB.h"
+#include "AliMUONCalibrationData.h"
 #include "AliMUONTriggerEfficiencyCells.h"
 #include "AliCDBManager.h"
 #include "AliCDBRunRange.h"
 ///
 /// \author Diego Stocco, INFN Torino
 
-void MUONTriggerChamberEfficiency(Bool_t addMapInSimulation=kFALSE,
-                                 const char *inputDir=".")
+void MUONTriggerChamberEfficiency(Char_t* inputFile="./MUON.TriggerEfficiencyMap.root",
+                                 Bool_t addMapInLocalOCDB = kFALSE,
+                                 Int_t firstRun=0, Int_t lastRun = AliCDBRunRange::Infinity(),
+                                 Bool_t useMeanValues = kFALSE)
 {
-/// \param addMapInSimulation (default kFALSE);
-///     kTRUE: creates file MUON/Calib/TriggerEfficiency/Run0_99999999_v0_s?.root
-///            with calculated chamber efficiency which can be used in the next simulation
-/// \param inputDir 
-///     path to AliESDs.root (default ".")
+/// \param inputFile (default "./MUON.TriggerEfficiencyMaps.root")
+///     File with the numerator and denominator histos for efficiency calculation
+///     (It is the output of the PWG3/muon/AliAnalysisTaskTrigChEff analysis
+/// \param addMapInLocalOCDB (default kFALSE)
+///     if the map is OK, add it in local OCDB
+/// \param firstRun (default 0)
+///     first run of validity for CDB object
+/// \param lastRun (default AliCDBRunRange::Infinity())
+///     last run of validity for CDB Object
+/// \param useMeanValues (default kFALSE)
+///     Fill local board maps with the average RPC efficiency
+///    (useful when the statistics is low (i.e. for cosmic runs))
 
-    Char_t filename[150], *className = "AliMUONTriggerEfficiencyCells";
-    sprintf(filename,"%s/AliESDs.root",inputDir);
+  AliCDBManager::Instance()->SetDefaultStorage("local://$ALICE_ROOT/OCDB");
 
-    TFile *file = new TFile(filename,"read");
-    if(!file){
-       cerr << "Cannot find " << filename << "\nExit!" << endl;
-       return;
-    }
-    
-    TTree *esdTree = (TTree*)file->Get("esdTree");
-    if(!esdTree){
-       cerr << "Cannot find esdTree in " << filename << "\nExit!" << endl;
-       return;
-    }
-    
-    AliMUONTriggerEfficiencyCells *effMap = 
-      (AliMUONTriggerEfficiencyCells*)esdTree->GetUserInfo()->FindObject(className);
-    if(!effMap){
-       cerr << "Cannot find " << className << " in esdTree.\nExit!" << endl;
-       return;
-    }
+  AliMUONTriggerEfficiencyCells* effMap = new AliMUONTriggerEfficiencyCells(inputFile);
+  if ( useMeanValues ){
+    AliCDBManager::Instance()->SetRun(firstRun);
+    effMap->LowStatisticsSettings();
+  }
 
+  if ( !addMapInLocalOCDB ){
+    AliCDBManager::Instance()->SetRun(firstRun);
     effMap->DisplayEfficiency();
+  }
+  else
+    AliMUONCDB::WriteToCDB("MUON/Calib/TriggerEfficiency", effMap, firstRun, lastRun,true);
+}
+
+//____________________________________________________________
+void ShowOCDBmap(Int_t runNumber = 0, TString ocdbPath="local://$ALICE_ROOT/OCDB")
+{
+/// \param runNumber (default 0)
+///     run number
+/// \param ocdbPath(default "local://$ALICE_ROOT/OCDB")
+///     path to OCDB
+  if ( ocdbPath.BeginsWith("alien://") || ocdbPath.BeginsWith("raw://"))
+    TGrid::Connect("alien://");
 
-    if(!addMapInSimulation) return;
+  AliCDBManager::Instance()->SetDefaultStorage(ocdbPath.Data());
+  AliCDBManager::Instance()->SetRun(runNumber);
+  AliMUONCalibrationData calib(runNumber);
 
-    AliCDBManager::Instance()->SetDefaultStorage("local://$ALICE_ROOT/OCDB");
-    AliMUONCDB::WriteToCDB("MUON/Calib/TriggerEfficiency",effMap,0,AliCDBRunRange::Infinity(),true);
+  calib.TriggerEfficiency()->DisplayEfficiency();
 }
+
index 8d6c6e6..4de1b68 100644 (file)
@@ -180,7 +180,8 @@ root [1] MUONTriggerEfficiencyPt()
 
 \section trigger_s5 How to get trigger chamber efficiency from data
 
-Trigger chamber efficiency map is calculated during reconstruction and saved in AliESDs.root
+Trigger chamber efficiency map is calculated from AliESDs.root with the analysis task: PWG3/muon/AliMUONTriggerChamberEff.
+The output of the task is a file filename="MUON.TriggerEfficiencyMap.root" containing the numerator and denominator histograms for the efficiency calculation.
 In order to view and save the map, use macro MUONTriggerChamberEfficiency.C
 
 To compile MUONTriggerChamberEfficiency.C
@@ -191,12 +192,24 @@ To compile MUONTriggerChamberEfficiency.C
 
 To run MUONTriggerChamberEfficiency.C
 <pre>
-MUONTriggerChamberEfficiency();
+MUONTriggerChamberEfficiency("filename");
 </pre>
 
-If you want to make the calculated map available for next simulation use option kTRUE, i.e.
+If you want to make the calculated map available for next simulation use option kTRUE, and specify the run validity:
 <pre>
-MUONTriggerChamberEfficiency(kTRUE);
+MUONTriggerChamberEfficiency("filename",kTRUE,firstRun,lastRun);
+</pre>
+
+The last option allows to set the efficiency of each board of an RPC with the mean value of the chamber: this is useful when the efficiency was calculated with low statistics (like in cosmic runs).
+<pre>
+MUONTriggerChamberEfficiency("filename",kFALSE,firstRun,lastRun,kTRUE);
+</pre>
+
+Finally, to have a look to the map stored in the OCDB for the run runNumber, do:
+<pre>
+.includepath $ALICE_ROOT/MUON
+.L $ALICE_ROOT/MUON/MUONTriggerChamberEfficiency.C+
+ShowOCDBmap(runNumber)
 </pre>
 
 When running next simulation, please remember to activate trigger efficiency