Skip to content

TRD digitizer signals to ADC O2-455 #2009

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

Merged
sawenzel merged 2 commits into from
May 27, 2019
Merged

TRD digitizer signals to ADC O2-455 #2009

Show file tree
Hide file tree
Changes from all commits
Commits
Show all changes
2 commits
Select commit Hold shift + click to select a range
5125ce5
signal2adc
jolopezl May 23, 2019
21134e2
correct driftv
jolopezl May 24, 2019
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
1 change: 1 addition & 0 deletions Detectors/TRD/simulation/include/TRDSimulation/Digitizer.h
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -29,6 +29,7 @@ class PadResponse;
class Digitizer
{
public:
enum { kTimeBins = 30 };
//
Digitizer();
~Digitizer() = default;
Expand Down
38 changes: 8 additions & 30 deletions Detectors/TRD/simulation/src/Digitizer.cxx
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -150,7 +150,7 @@ bool Digitizer::convertHits(const int det, const std::vector<HitType>& hits, TRD
timeBinTRFend = ((int)(mSimParam->GetTRFhi() * mCommonParam->GetSamplingFrequency())) - 1;
}

const int nTimeTotal = 0; // DigitsManager->GetDigitsParam()->GetNTimeBins(det);
const int nTimeTotal = kTimeBins; // DigitsManager->GetDigitsParam()->GetNTimeBins(det);
const float samplingRate = mCommonParam->GetSamplingFrequency();
const float elAttachProp = mSimParam->GetElAttachProp() / 100;

Expand Down Expand Up @@ -218,7 +218,7 @@ bool Digitizer::convertHits(const int det, const std::vector<HitType>& hits, TRD
absDriftLength /= TMath::Sqrt(1 / (1 + calExBDetValue * calExBDetValue));
}
// double driftVelocity = calVdriftDetValue * calVdriftROC->GetValue(colE, rowE); PLEASE FIX ME when CCDB is ready
double driftVelocity = 2150; // Defaults values from OCDB (AliRoot DrawTrending macro) for 5 TeV pp - 27 runs from LHC15n
double driftVelocity = 2.13; // Defaults values from OCDB (AliRoot DrawTrending macro) for 5 TeV pp - 27 runs from LHC15n

// Loop over all created electrons
const int nElectrons = abs(qTotal);
Expand Down Expand Up @@ -274,7 +274,7 @@ bool Digitizer::convertHits(const int det, const std::vector<HitType>& hits, TRD

// Retrieve drift velocity becuase col and row may have changed
// driftVelocity = calVdriftDetValue* calVdriftROC->GetValue(colE, rowE); PLEASE FIX ME when CCDB is ready
driftVelocity = 2150; // Defaults values from OCDB (AliRoot DrawTrending macro) for 5 TeV pp - 27 runs from LHC15n
driftVelocity = 2.13; // Defaults values from OCDB (AliRoot DrawTrending macro) for 5 TeV pp - 27 runs from LHC15n
// float t0 = calT0DetValue + calT0ROC->getValue(colE, rowE); PLEASE FIX ME when CCDB is ready
const float t0 = -1.38 + 0; // Defaults values from OCDB (AliRoot DrawTrending macro) for 5 TeV pp - 27 runs from LHC15n

Expand Down Expand Up @@ -432,49 +432,27 @@ bool Digitizer::convertSignalsToADC(const int det, int& signals)

int nRowMax = mGeo->getPadPlane(det)->getNrows();
int nColMax = mGeo->getPadPlane(det)->getNcols();
int nTimeTotal = 100; // fDigitsManager->GetDigitsParam()->GetNTimeBins(det);
int nTimeTotal = kTimeBins; // fDigitsManager->GetDigitsParam()->GetNTimeBins(det);
// if (fSDigitsManager->GetDigitsParam()->GetNTimeBins(det)) {
// nTimeTotal = fSDigitsManager->GetDigitsParam()->GetNTimeBins(det);
// } else {
// LOG(FATAL) << "Could not get number of time bins";
// return false;
// }

// The gain factor mCalib objects
// const AliTRDCalDet* calGainFactorDet = mCalib->GetGainFactorDet();
// AliTRDCalROC* calGainFactorROC = 0x0;
float calGainFactorDetValue = 0.0;

// AliTRDarrayADC* digits = 0x0;
int digits = 0;

if (!signals) {
LOG(FATAL) << "Signals array for detector " << det << " does not exis";
return false;
}
// if (signals->HasData()) {
// // Expand the container if neccessary
// signals->Expand();
// } else {
// // Create missing containers
// signals->Allocate(nRowMax, nColMax, nTimeTotal);
// }

// Get the container for the digits of this detector
// if (fDigitsManager->HasSDigits()) {
// AliError("Digits manager has s-digits");
// return false;
// }

// digits = (AliTRDarrayADC*)fDigitsManager->GetDigits(det);
// // Allocate memory space for the digits buffer
// if (!digits->HasData()) {
// digits->Allocate(nRowMax, nColMax, nTimeTotal);
// }

// Get the mCalib objects
// calGainFactorROC = mCalib->GetGainFactorROC(det);
// CalDet* calGainFactorDet = mCalib->GetGainFactorDet();
// CalRoc* calGainFactorROC = mCalib->GetGainFactorROC(det);
// calGainFactorDetValue = calGainFactorDet->GetValue(det);
float calGainFactorDetValue = 0.47; // +/- 0.06 // Defaults value from OCDB (AliRoot DrawTrending macro) for 5 TeV pp - 27 runs from LHC15n

// Create the digits for this chamber
for (row = 0; row < nRowMax; row++) {
Expand All @@ -495,7 +473,7 @@ bool Digitizer::convertSignalsToADC(const int det, int& signals)
// }

// The gain factors
float padgain = calGainFactorDetValue; // * calGainFactorROC->GetValue(col, row);
float padgain = calGainFactorDetValue; // * calGainFactorROC->GetValue(col, row); // PLEASE FIX ME when CCDB is ready
if (padgain <= 0) {
const auto msg = Form("Not a valid gain %f, %d %d %d", padgain, det, col, row);
LOG(FATAL) << msg;
Expand Down