Skip to content
New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

Sign up for GitHub

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

TPC/simulation: Add PT correction to gas density using configurable p… #13048

Merged
merged 3 commits into from
Apr 23, 2024
Merged

TPC/simulation: Add PT correction to gas density using configurable p… #13048

Show file tree
Hide file tree
Changes from all commits
Commits
Show all changes
3 commits
Select commit Hold shift + click to select a range
1c98c1d
TPC/simulation: Add PT correction to gas density using configurable p…
tubagundem Apr 19, 2024
12d1ba8
Formatting
tubagundem Apr 19, 2024
53c89ab
Implemented feedback
tubagundem Apr 19, 2024
File filter

Filter by extension

Filter by extension
Conversations
Failed to load comments.
Loading
Jump to
Jump to file
Failed to load files.
Loading
Diff view
Diff view
2 changes: 2 additions & 0 deletions Detectors/TPC/base/include/TPCBase/ParameterGas.h
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -39,6 +39,8 @@ struct ParameterGas : public o2::conf::ConfigurableParamHelper<ParameterGas> {
float Nprim = 14.f; ///< Number of primary electrons per MIP and cm [1/cm]
float ScaleFactorG4 = 0.85f; ///< Scale factor to tune WION for GEANT4
float FanoFactorG4 = 0.7f; ///< Parameter for smearing the number of ionizations (nel) using GEANT4
float Pressure = 1013.25f; ///< Pressure [mbar]
float Temperature = 20.0f; ///< Temperature [°C]
float BetheBlochParam[5] = {0.820172e-1f, 9.94795f, 8.97292e-05f, 2.05873f, 1.65272f}; ///< Parametrization of Bethe-Bloch

O2ParamDef(ParameterGas, "TPCGasParam");
Expand Down
29 changes: 24 additions & 5 deletions Detectors/TPC/simulation/src/Detector.cxx
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -315,6 +315,8 @@ void Detector::CreateMaterials()
// Origin: Marek Kowalski IFJ, Krakow, [email protected]
//-----------------------------------------------------------------

const auto& gasParam = ParameterGas::Instance();

Int_t iSXFLD = 2;
Float_t sXMGMX = 10.0;
// init the field tracking params
Expand All @@ -326,6 +328,23 @@ void Detector::CreateMaterials()

Float_t density;

// TODO: load pressure and temperature values from CCDB
const Double_t pressure = gasParam.Pressure; // in mbar
const Double_t temperature = gasParam.Temperature + 273.15; // in K

// densities were taken for these values
const Double_t t1 = 293.15; // 20°C in K
const Double_t p1 = 1013.25; // 1 atm in mbars

// sanity check - temperature between 10 and 30 deg, pressure between 800 and 1200 mbar
Double_t ptCorr = 1.;
if (TMath::Abs(temperature - 293.15) > 10. || TMath::Abs(pressure - 1000.) > 200.) {
ptCorr = 1.;
} else {
ptCorr = (pressure * t1) / (p1 * temperature);
}
LOG(info) << "Setting gas density correction to: " << ptCorr;

//***************** Gases *************************

//--------------------------------------------------------------
Expand All @@ -345,7 +364,7 @@ void Detector::CreateMaterials()

density = 1.842e-3;

o2::base::Detector::Mixture(10, "CO2", amat, zmat, density, 2, wmat);
o2::base::Detector::Mixture(10, "CO2", amat, zmat, density * ptCorr, 2, wmat);
//
// Air
//
Expand All @@ -360,7 +379,7 @@ void Detector::CreateMaterials()
//
density = 0.001205;

o2::base::Detector::Mixture(11, "Air", amat, zmat, density, 2, wmat);
o2::base::Detector::Mixture(11, "Air", amat, zmat, density * ptCorr, 2, wmat);

//----------------------------------------------------------------
// drift gases 5 mixtures, 5 materials
Expand Down Expand Up @@ -466,9 +485,9 @@ void Detector::CreateMaterials()
}

//
o2::base::Detector::Mixture(12, gname1.Data(), amat1, zmat1, density, cnt, wmat1); // nonsensitive
o2::base::Detector::Mixture(13, gname2.Data(), amat1, zmat1, density, cnt, wmat1); // sensitive
o2::base::Detector::Mixture(40, gname3.Data(), amat1, zmat1, density, cnt, wmat1); // sensitive Kr
o2::base::Detector::Mixture(12, gname1.Data(), amat1, zmat1, density * ptCorr, cnt, wmat1); // nonsensitive
o2::base::Detector::Mixture(13, gname2.Data(), amat1, zmat1, density * ptCorr, cnt, wmat1); // sensitive
o2::base::Detector::Mixture(40, gname3.Data(), amat1, zmat1, density * ptCorr, cnt, wmat1); // sensitive Kr

//----------------------------------------------------------------------
// solid materials
Expand Down
Loading