ChamberEfficiencyResponse.cxx

Go to the documentation of this file.

// Copyright 2019-2020 CERN and copyright holders of ALICE O2.
// See https://alice-o2.web.cern.ch/copyright for details of the copyright holders.
// All rights not expressly granted are reserved.
//
// This software is distributed under the terms of the GNU General Public
// License v3 (GPL Version 3), copied verbatim in the file "COPYING".
//
// In applying this license CERN does not waive the privileges and immunities
// granted to it by virtue of its status as an Intergovernmental Organization
// or submit itself to any jurisdiction.
 
 
#include "MIDSimulation/ChamberEfficiencyResponse.h"
 
namespace o2
{
namespace mid
{
//______________________________________________________________________________
ChamberEfficiencyResponse::ChamberEfficiencyResponse(const ChamberEfficiency& efficiencyMap) : mEfficiency(efficiencyMap)
{
}
 
//______________________________________________________________________________
bool ChamberEfficiencyResponse::isEfficient(int deId, int columnId, int line, bool& isEfficientBP, bool& isEfficientNBP)
{
 
  // P(B) : probability to fire bending plane
  double effBP = mEfficiency.getEfficiency(deId, columnId, line, ChamberEfficiency::EffType::BendPlane);
 
  // P(BN) : probability to fire bending and non-bending plane
  double effBoth = mEfficiency.getEfficiency(deId, columnId, line, ChamberEfficiency::EffType::BothPlanes);
 
  isEfficientBP = (mRandom(mGenerator) <= effBP);
 
  // P(N) : probability to fire non-bending plane
  double effNBP = 0.;
  if (isEfficientBP) {
    // P(N|B) = P(BN) / P(B)
    effNBP = effBoth / effBP;
  } else {
    // P(N|!B) = ( P(N) - P(BN) ) / ( 1 - P(B) )
    effNBP = (mEfficiency.getEfficiency(deId, columnId, line, ChamberEfficiency::EffType::NonBendPlane) - effBoth) / (1. - effBP);
  }
 
  isEfficientNBP = (mRandom(mGenerator) <= effNBP);
 
  return (isEfficientBP || isEfficientNBP);
}
 
ChamberEfficiencyResponse createDefaultChamberEfficiencyResponse()
{
  return ChamberEfficiencyResponse(createDefaultChamberEfficiency());
}
 
} // namespace mid
} // namespace o2