ChamberResponse.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/ChamberResponse.h"
 
#include <cmath>
 
namespace o2
{
namespace mid
{
//______________________________________________________________________________
ChamberResponse::ChamberResponse(const ChamberResponseParams& params, const ChamberHV& hv) : mParams(params), mHV(hv)
{
}
 
//______________________________________________________________________________
double ChamberResponse::getFiredProbability(double distance, int cathode, int deId, double theta) const
{
 
  // Need to convert the distance from cm to mm
  double distMM = distance * 10.;
  double parA = mParams.getParA(mHV.getHV(deId));
  double parB = mParams.getParB(cathode, deId);
  double parC = mParams.getParC(mHV.getHV(deId));
  double costheta = std::cos(theta);
  return (parC + parA / (parA + costheta * std::pow(distMM, parB))) / (1 + parC);
}
 
//______________________________________________________________________________
double ChamberResponse::firedProbabilityFunction(double* var, double* par)
{
 
  int cathode = (int)par[0];
  int deId = (int)par[1];
 
  mParams.setParB(cathode, deId, par[3]);
  mParams.setParA(par[4], par[5]);
  mParams.setParC(par[6], par[7]);
 
  return getFiredProbability(var[0], cathode, deId, par[2]);
}
 
//______________________________________________________________________________
ChamberResponse createDefaultChamberResponse()
{
  ChamberResponseParams params = createDefaultChamberResponseParams();
  ChamberHV hv = createDefaultChamberHV();
  return ChamberResponse(params, hv);
}
 
} // namespace mid
} // namespace o2