1 /**************************************************************************
2 * Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved. *
4 * Author: The ALICE Off-line Project. *
5 * Contributors are mentioned in the code where appropriate. *
7 * Permission to use, copy, modify and distribute this software and its *
8 * documentation strictly for non-commercial purposes is hereby granted *
9 * without fee, provided that the above copyright notice appears in all *
10 * copies and that both the copyright notice and this permission notice *
11 * appear in the supporting documentation. The authors make no claims *
12 * about the suitability of this software for any purpose. It is *
13 * provided "as is" without express or implied warranty. *
14 **************************************************************************/
16 /** @file AliFMDParameters.cxx
17 @author Christian Holm Christensen <cholm@nbi.dk>
18 @date Mon Mar 27 12:44:26 2006
19 @brief Manager of FMD parameters
21 //____________________________________________________________________
23 // Forward Multiplicity Detector based on Silicon wafers.
25 // This class is a singleton that handles various parameters of
27 // The manager normally serves the parameters from the Conditions
28 // Database (CDB). These are retrivied by the member function
29 // `Init'. Optionally, the class can serve hard-coded constants, if
30 // no CDB is available.
32 #include "AliFMDDebug.h" // ALILOG_H
33 #include "AliFMDParameters.h" // ALIFMDPARAMETERS_H
34 #include "AliFMDGeometry.h" // ALIFMDGEOMETRY_H
35 #include "AliFMDRing.h" // ALIFMDRING_H
36 #include "AliFMDCalibGain.h" // ALIFMDCALIBGAIN_H
37 #include "AliFMDCalibPedestal.h" // ALIFMDCALIBPEDESTAL_H
38 #include "AliFMDCalibSampleRate.h" // ALIFMDCALIBSAMPLERATE_H
39 #include "AliFMDCalibStripRange.h" // ALIFMDCALIBSTRIPRANGE_H
40 #include "AliFMDAltroMapping.h" // ALIFMDALTROMAPPING_H
41 #include <AliCDBManager.h> // ALICDBMANAGER_H
42 #include <AliCDBEntry.h> // ALICDBMANAGER_H
43 #include <AliFMDPreprocessor.h>
45 #include <Riostream.h>
51 //====================================================================
52 ClassImp(AliFMDParameters)
54 ; // This is here to keep Emacs for indenting the next line
57 //____________________________________________________________________
58 AliFMDParameters* AliFMDParameters::fgInstance = 0;
60 //____________________________________________________________________
61 const char* AliFMDParameters::fgkPulseGain = "FMD/Calib/PulseGain";
62 const char* AliFMDParameters::fgkPedestal = "FMD/Calib/Pedestal";
63 const char* AliFMDParameters::fgkDead = "FMD/Calib/Dead";
64 const char* AliFMDParameters::fgkSampleRate = "FMD/Calib/SampleRate";
65 const char* AliFMDParameters::fgkAltroMap = "FMD/Calib/AltroMap";
66 const char* AliFMDParameters::fgkZeroSuppression = "FMD/Calib/ZeroSuppression";
67 const char* AliFMDParameters::fgkStripRange = "FMD/Calib/StripRange";
68 const char* AliFMDParameters::fkPedestalShuttleID = "pedestals";
69 const char* AliFMDParameters::fkGainShuttleID = "gains";
70 const char* AliFMDParameters::fkConditionsShuttleID = "conditions";
72 //____________________________________________________________________
74 AliFMDParameters::Instance()
76 // Get static instance
77 if (!fgInstance) fgInstance = new AliFMDParameters;
81 //____________________________________________________________________
82 AliFMDParameters::AliFMDParameters()
90 fFixedPedestalWidth(0),
91 fFixedZeroSuppression(0),
98 fHasRcuTrailer(kTRUE),
99 fHasCompleteHeader(kTRUE),
108 // Default constructor
110 SetAltroChannelSize();
111 SetChannelsPerAltro();
112 SetZeroSuppression();
119 fAltroMap = new AliFMDAltroMapping;
122 //__________________________________________________________________
124 AliFMDParameters::Init(Bool_t forceReInit, UInt_t what)
126 // Initialize the parameters manager. We need to get stuff from the
128 if (forceReInit) fIsInit = kFALSE;
130 if (what & kPulseGain) InitPulseGain();
131 if (what & kPedestal) InitPedestal();
132 if (what & kDeadMap) InitDeadMap();
133 if (what & kSampleRate) InitSampleRate();
134 if (what & kZeroSuppression) InitZeroSuppression();
135 if (what & kAltroMap) InitAltroMap();
138 //__________________________________________________________________
140 AliFMDParameters::Init(AliFMDPreprocessor* pp, Bool_t forceReInit, UInt_t what)
142 // Initialize the parameters manager. We need to get stuff from the
144 if (forceReInit) fIsInit = kFALSE;
146 if (what & kPulseGain) InitPulseGain(pp);
147 if (what & kPedestal) InitPedestal(pp);
148 if (what & kDeadMap) InitDeadMap(pp);
149 if (what & kSampleRate) InitSampleRate(pp);
150 if (what & kZeroSuppression) InitZeroSuppression(pp);
151 if (what & kAltroMap) InitAltroMap(pp);
155 //__________________________________________________________________
156 #define DET2IDX(det,ring,sec,str) \
157 (det * 1000 + (ring == 'I' ? 0 : 512) + str)
159 //__________________________________________________________________
161 AliFMDParameters::Draw(Option_t* option)
165 kPulseGain, // Path to PulseGain calib object
166 kThreshold, // Path to PulseGain calib object
167 kPedestal, // Path to Pedestal calib object
168 kPedestalWidth, // Path to Pedestal calib object
169 kDead, // Path to Dead calib object
170 kSampleRate, // Path to SampleRate calib object
171 kAltroMap, // Path to AltroMap calib object
172 kZeroSuppression, // Path to ZeroSuppression cal object
173 kMinStripRange, // Path to strip range cal object
174 kMaxStripRange // Path to strip range cal object
178 if (opt.Contains("dead", TString::kIgnoreCase))
180 else if (opt.Contains("threshold",TString::kIgnoreCase))
182 else if (opt.Contains("gain",TString::kIgnoreCase))
184 else if (opt.Contains("pedestal",TString::kIgnoreCase))
186 else if (opt.Contains("noise",TString::kIgnoreCase))
187 what = kPedestalWidth;
188 else if (opt.Contains("zero",TString::kIgnoreCase))
189 what = kZeroSuppression;
190 else if (opt.Contains("rate",TString::kIgnoreCase))
192 else if (opt.Contains("min",TString::kIgnoreCase))
193 what = kMinStripRange;
194 else if (opt.Contains("max",TString::kIgnoreCase))
195 what = kMaxStripRange;
196 else if (opt.Contains("map",TString::kIgnoreCase))
199 Warning("Draw", "unknown parameter: %s\n\tShould be one of\n\t"
200 "dead, threshold, gain, pedestal, noise, zero, rate, "
206 TArrayD xbins(3 * 512 + 2 * 256 + 5);
209 for (UShort_t det = 1; det <= 3; det++) {
210 UShort_t nRings = (det == 1 ? 1 : 2);
211 for (UShort_t iring = 0; iring < nRings; iring++) {
212 UShort_t nStrip = (iring == 0 ? 512 : 256);
213 Char_t ring = (iring == 0 ? 'I' : 'O');
214 for (UShort_t str = 0; str < nStrip; str++) {
215 // UShort_t nSec = (iring == 0 ? 20 : 40);
216 // Char_t ring = (iring == 0 ? 'I' : 'O');
217 // for (UShort_t sec = 0; sec < nSec; sec++) {
218 Int_t idx = DET2IDX(det, ring, 0, str);
219 // Int_t idx = DET2IDX(det, ring, sec, 0);
221 xbins[i-1] = idx - .5;
232 for (Int_t i = 0; i < ybins.fN; i++) ybins[i] = Float_t(i - .5);
233 TH2D* hist = new TH2D("calib", Form("Calibration %s", option),
234 xbins.fN-1, xbins.fArray,
235 ybins.fN-1, ybins.fArray);
236 hist->GetXaxis()->SetTitle("1000 #times detector + 512 #times ring + strip");
237 hist->GetYaxis()->SetTitle("sector");
239 // hist->Draw("Lego");
242 for (UShort_t det = 1; det <= 3; det++) {
243 UShort_t nRings = (det == 1 ? 1 : 2);
244 for (UShort_t iring = 0; iring < nRings; iring++) {
245 UShort_t nSector = (iring == 0 ? 20 : 40);
246 UShort_t nStrip = (iring == 0 ? 512 : 256);
247 Char_t ring = (iring == 0 ? 'I' : 'O');
248 for (UShort_t sec = 0; sec < nSector; sec++) {
249 for (UShort_t str = 0; str < nStrip; str++) {
250 Int_t idx = DET2IDX(det, ring, sec, str);
254 case kPulseGain: // Path to PulseGain calib object
255 val = GetPulseGain(det,ring,sec,str); break;
256 case kThreshold: // Path to PulseGain calib object
257 val = GetThreshold(); break;
258 case kPedestal: // Path to Pedestal calib object
259 val = GetPedestal(det,ring,sec,str); break;
260 case kPedestalWidth: // Path to Pedestal calib object
261 val = GetPedestalWidth(det,ring,sec,str); break;
262 case kDead: // Path to Dead calib object
263 val = IsDead(det,ring,sec,str); break;
264 case kSampleRate: // Path to SampleRate calib object
265 val = GetSampleRate(det,ring,sec,str); break;
266 case kAltroMap: // Path to AltroMap calib object
267 Detector2Hardware(det,ring,sec,str, ddl, addr);
269 case kZeroSuppression: // Path to ZeroSuppression cal object
270 val = GetZeroSuppression(det,ring,sec,str); break;
271 case kMinStripRange: // Path to strip range cal object
272 val = GetMinStrip(det,ring,sec,str); break;
273 case kMaxStripRange: // Path to strip range cal object
274 val = GetMaxStrip(det,ring,sec,str); break;
276 hist->Fill(idx,sec,val);
277 // hist->Fill(idx,str,val);
285 //__________________________________________________________________
287 AliFMDParameters::Print(Option_t* option) const
289 // Print information.
290 // If option contains an 'A' then everything is printed.
291 // If the option contains the string "FMD" the function will search
292 // for detector, ring, sector, and strip numbers to print, in the
295 // FMD<detector><ring>[<sector>,<string>]
297 // The wild card '*' means all of <detector>, <ring>, <sector>, or
300 Bool_t showStrips = opt.Contains("a", TString::kIgnoreCase);
301 UShort_t ds[] = { 1, 2, 3, 0 };
302 Char_t rs[] = { 'I', 'O', '\0' };
303 UShort_t minStrip = 0;
304 UShort_t maxStrip = 512;
305 UShort_t minSector = 0;
306 UShort_t maxSector = 40;
309 if (opt.Contains("fmd",TString::kIgnoreCase)) {
311 size_t i = opt.Index("fmd",TString::kIgnoreCase);
312 size_t j = opt.Index("]",TString::kIgnoreCase);
323 std::stringstream s(opt(i+4, j-i-3).Data());
324 while (state != kEnd) {
325 Char_t tmp = s.peek();
326 if (tmp == ' ' || tmp == '\t') {
331 case kReadDet: { // First, try to kRead the detector
332 if (tmp == '*') s.get();
341 state = (s.bad() ? kEnd : kReadRing);
343 case kReadRing: { // Then try to read the ring;
346 if (ring != '*' && !s.bad()) {
350 state = (s.bad() ? kEnd : kReadLbrack);
352 case kReadLbrack: { // Try to read a left bracket
355 state = (s.bad() ? kEnd : kReadSector);
357 case kReadSector: { // Try to read a sector
358 if (tmp == '*') s.get();
367 state = (s.bad() ? kEnd : kReadComma);
369 case kReadComma: { // Try to read a left bracket
372 state = (s.bad() ? kEnd : kReadStrip);
374 case kReadStrip: { // Try to read a strip
375 if (tmp == '*') s.get();
384 state = (s.bad() ? kEnd : kReadRbrack);
386 case kReadRbrack: { // Try to read a left bracket
398 while ((det = *(dp++))) {
402 while ((ring = *(rp++))) {
403 if (det == 1 && ring == 'O') continue;
404 UShort_t min = GetMinStrip(det, ring, 0, 0);
405 UShort_t max = GetMaxStrip(det, ring, 0, 0);
406 UShort_t rate = GetSampleRate(det, ring, 0, 0);
407 std::cout << "FMD" << det << ring
409 << std::setw(3) << min << ","
410 << std::setw(3) << max << " Rate: "
411 << std::setw(2) << rate << std::endl;
413 if (!showStrips) continue;
414 UShort_t nSec = ( ring == 'I' ? 20 : 40 );
415 UShort_t nStr = ( ring == 'I' ? 512 : 256 );
416 for (UShort_t sec = minSector; sec < maxSector && sec < nSec; sec++) {
418 << " Strip | Pedestal | Gain | ZS thr. | Address\n"
419 << "--------+-------------------+------------+---------+---------"
421 for (UShort_t str = minStrip; str < nStr && str < maxStrip; str++) {
422 if (str == minStrip) std::cout << std::setw(3) << sec << ",";
423 else std::cout << " ";
424 std::cout << std::setw(3) << str << " | ";
425 if (IsDead(det, ring, sec, str)) {
426 std::cout << "dead" << std::endl;
430 Detector2Hardware(det, ring, sec, str, ddl, addr);
431 std::cout << std::setw(7) << GetPedestal(det, ring, sec, str)
432 << "+/-" << std::setw(7)
433 << GetPedestalWidth(det, ring, sec, str)
434 << " | " << std::setw(10)
435 << GetPulseGain(det, ring, sec, str)
436 << " | " << std::setw(7)
437 << GetZeroSuppression(det, ring, sec, str)
438 << " | 0x" << std::hex << std::setw(4)
439 << std::setfill('0') << ddl << ",0x" << std::setw(3)
440 << addr << std::dec << std::setfill(' ') << std::endl;
444 << "============================================================="
451 //__________________________________________________________________
453 AliFMDParameters::SetStripRange(UShort_t min, UShort_t max)
455 // Set fixed strip range
456 fFixedMinStrip = min;
457 fFixedMaxStrip = max;
460 //__________________________________________________________________
462 AliFMDParameters::GetEntry(const char* path, AliFMDPreprocessor* pp,
465 // Get an entry from the CDB or via preprocessor
466 AliCDBEntry* entry = 0;
468 AliCDBManager* cdb = AliCDBManager::Instance();
469 entry = cdb->Get(path);
472 const char* third = gSystem->BaseName(path);
473 const char* second = gSystem->BaseName(gSystem->DirName(path));
474 entry = pp->GetFromCDB(second, third);
477 TString msg(Form("No %s found in CDB, perhaps you need to "
478 "use AliFMDCalibFaker?", path));
479 if (fatal) { AliFatal(msg.Data()); }
480 else AliLog::Message(AliLog::kWarning, msg.Data(), "FMD",
481 "AliFMDParameters", "GetEntry", __FILE__,
489 //__________________________________________________________________
491 AliFMDParameters::InitPulseGain(AliFMDPreprocessor* pp)
493 // Get pulse gain from CDB or used fixed
494 AliCDBEntry* gain = GetEntry(fgkPulseGain, pp);
497 AliFMDDebug(1, ("Got gain from CDB"));
498 fPulseGain = dynamic_cast<AliFMDCalibGain*>(gain->GetObject());
499 if (!fPulseGain) AliFatal("Invalid pulser gain object from CDB");
501 //__________________________________________________________________
503 AliFMDParameters::InitPedestal(AliFMDPreprocessor* pp)
505 // Initialize the pedestals from CDB
506 AliCDBEntry* pedestal = GetEntry(fgkPedestal, pp);
507 if (!pedestal) return;
509 AliFMDDebug(1, ("Got pedestal from CDB"));
510 fPedestal = dynamic_cast<AliFMDCalibPedestal*>(pedestal->GetObject());
511 if (!fPedestal) AliFatal("Invalid pedestal object from CDB");
514 //__________________________________________________________________
516 AliFMDParameters::InitDeadMap(AliFMDPreprocessor* pp)
518 // Get Dead-channel-map from CDB
519 AliCDBEntry* deadMap = GetEntry(fgkDead, pp);
520 if (!deadMap) return;
522 AliFMDDebug(1, ("Got dead map from CDB"));
523 fDeadMap = dynamic_cast<AliFMDCalibDeadMap*>(deadMap->GetObject());
524 if (!fDeadMap) AliFatal("Invalid dead map object from CDB");
527 //__________________________________________________________________
529 AliFMDParameters::InitZeroSuppression(AliFMDPreprocessor* pp)
531 // Get 0-suppression from CDB
532 AliCDBEntry* zeroSup = GetEntry(fgkZeroSuppression, pp);
533 if (!zeroSup) return;
534 AliFMDDebug(1, ("Got zero suppression from CDB"));
536 dynamic_cast<AliFMDCalibZeroSuppression*>(zeroSup->GetObject());
537 if (!fZeroSuppression)AliFatal("Invalid zero suppression object from CDB");
540 //__________________________________________________________________
542 AliFMDParameters::InitSampleRate(AliFMDPreprocessor* pp)
544 // get Sample rate from CDB
545 AliCDBEntry* sampRat = GetEntry(fgkSampleRate, pp);
546 if (!sampRat) return;
547 AliFMDDebug(1, ("Got zero suppression from CDB"));
548 fSampleRate = dynamic_cast<AliFMDCalibSampleRate*>(sampRat->GetObject());
549 if (!fSampleRate) AliFatal("Invalid zero suppression object from CDB");
552 //__________________________________________________________________
554 AliFMDParameters::InitAltroMap(AliFMDPreprocessor* pp)
556 // Get hardware mapping from CDB
561 AliCDBEntry* hwMap = GetEntry(fgkAltroMap, pp);
564 AliFMDDebug(1, ("Got ALTRO map from CDB"));
565 fAltroMap = dynamic_cast<AliFMDAltroMapping*>(hwMap->GetObject());
567 AliFatal("Invalid ALTRO map object from CDB");
568 fAltroMap = new AliFMDAltroMapping;
572 //__________________________________________________________________
574 AliFMDParameters::InitStripRange(AliFMDPreprocessor* pp)
576 // Get strips read-out from CDB
577 AliCDBEntry* range = GetEntry(fgkStripRange, pp);
579 AliFMDDebug(1, ("Got strip range from CDB"));
580 fStripRange = dynamic_cast<AliFMDCalibStripRange*>(range->GetObject());
581 if (!fStripRange) AliFatal("Invalid strip range object from CDB");
585 //__________________________________________________________________
587 AliFMDParameters::GetThreshold() const
589 // Get threshold from CDB
590 if (!fPulseGain) return fFixedThreshold;
591 return fPulseGain->Threshold();
594 //__________________________________________________________________
596 AliFMDParameters::GetPulseGain(UShort_t detector, Char_t ring,
597 UShort_t sector, UShort_t strip) const
599 // Returns the pulser calibrated gain for strip # strip in sector #
600 // sector or ring id ring of detector # detector.
602 // For simulation, this is normally set to
605 // ------------------ * MIP_Energy_Loss
606 // ALTRO_channel_size
609 if (fFixedPulseGain <= 0)
610 fFixedPulseGain = fVA1MipRange * GetEdepMip() / fAltroChannelSize;
611 return fFixedPulseGain;
613 AliFMDDebug(50, ("pulse gain for FMD%d%c[%2d,%3d]=%f",
614 detector, ring, sector, strip,
615 fPulseGain->Value(detector, ring, sector, strip)));
616 return fPulseGain->Value(detector, ring, sector, strip);
619 //__________________________________________________________________
621 AliFMDParameters::IsDead(UShort_t detector, Char_t ring,
622 UShort_t sector, UShort_t strip) const
624 // Check if the channel is dead
625 if (!fDeadMap) return kFALSE;
626 AliFMDDebug(50, ("Dead for FMD%d%c[%2d,%3d]=%s",
627 detector, ring, sector, strip,
628 fDeadMap->operator()(detector, ring, sector, strip) ?
630 return fDeadMap->operator()(detector, ring, sector, strip);
633 //__________________________________________________________________
635 AliFMDParameters::GetZeroSuppression(UShort_t detector, Char_t ring,
636 UShort_t sector, UShort_t strip) const
638 // Get zero suppression threshold
639 if (!fZeroSuppression) return fFixedZeroSuppression;
640 // Need to map strip to ALTRO chip.
641 AliFMDDebug(50, ("zero sup. for FMD%d%c[%2d,%3d]=%f",
642 detector, ring, sector, strip,
643 fZeroSuppression->operator()(detector, ring,
645 return fZeroSuppression->operator()(detector, ring, sector, strip/128);
648 //__________________________________________________________________
650 AliFMDParameters::GetSampleRate(UShort_t det, Char_t ring, UShort_t sector,
654 if (!fSampleRate) return fFixedSampleRate;
655 // Need to map sector to digitizier card.
656 UInt_t ret = fSampleRate->Rate(det, ring, sector, str);
657 AliFMDDebug(50, ("Sample rate for FMD%d%c[%2d,%3d]=%d",
658 det, ring, sector, str, ret));
662 //__________________________________________________________________
664 AliFMDParameters::GetMinStrip(UShort_t det, Char_t ring, UShort_t sector,
667 // Get strip range read out
668 if (!fStripRange) return fFixedMinStrip;
669 // Need to map sector to digitizier card.
670 UInt_t ret = fStripRange->Min(det, ring, sector, str);
671 AliFMDDebug(50, ("Min strip # for FMD%d%c[%2d,%3d]=%d",
672 det, ring, sector, str, ret));
676 //__________________________________________________________________
678 AliFMDParameters::GetMaxStrip(UShort_t det, Char_t ring, UShort_t sector,
681 // Get strip range read out
682 if (!fStripRange) return fFixedMaxStrip;
683 // Need to map sector to digitizier card.
684 UInt_t ret = fStripRange->Max(det, ring, sector, str);
685 AliFMDDebug(50, ("Max strip # for FMD%d%c[%2d,%3d]=%d",
686 det, ring, sector, str, ret));
690 //__________________________________________________________________
692 AliFMDParameters::GetPedestal(UShort_t detector, Char_t ring,
693 UShort_t sector, UShort_t strip) const
696 if (!fPedestal) return fFixedPedestal;
697 AliFMDDebug(50, ("pedestal for FMD%d%c[%2d,%3d]=%f",
698 detector, ring, sector, strip,
699 fPedestal->Value(detector, ring, sector, strip)));
700 return fPedestal->Value(detector, ring, sector, strip);
703 //__________________________________________________________________
705 AliFMDParameters::GetPedestalWidth(UShort_t detector, Char_t ring,
706 UShort_t sector, UShort_t strip) const
709 if (!fPedestal) return fFixedPedestalWidth;
710 AliFMDDebug(50, ("pedetal width for FMD%d%c[%2d,%3d]=%f",
711 detector, ring, sector, strip,
712 fPedestal->Width(detector, ring, sector, strip)));
713 return fPedestal->Width(detector, ring, sector, strip);
716 //__________________________________________________________________
718 AliFMDParameters::GetAltroMap() const
720 // Get the hardware address to detector index map
725 //__________________________________________________________________
727 AliFMDParameters::Hardware2Detector(UInt_t ddl, UInt_t board,
728 UInt_t chip, UInt_t chan,
729 UShort_t& det, Char_t& ring,
730 UShort_t& sec, UShort_t& str) const
732 // Map hardware address to detector index
733 if (!fAltroMap) return kFALSE;
734 return fAltroMap->Hardware2Detector(ddl,board,chip,chan, det,ring,sec,str);
736 //__________________________________________________________________
738 AliFMDParameters::Hardware2Detector(UInt_t ddl, UInt_t addr,
739 UShort_t& det, Char_t& ring,
740 UShort_t& sec, UShort_t& str) const
742 // Map hardware address to detector index
743 if (!fAltroMap) return kFALSE;
744 return fAltroMap->Hardware2Detector(ddl, addr, det, ring, sec, str);
747 //__________________________________________________________________
749 AliFMDParameters::Detector2Hardware(UShort_t det, Char_t ring,
750 UShort_t sec, UShort_t str,
751 UInt_t& ddl, UInt_t& board,
752 UInt_t& chip, UInt_t& chan) const
754 // Map detector index to hardware address
755 if (!fAltroMap) return kFALSE;
756 return fAltroMap->Detector2Hardware(det,ring,sec,str, ddl,board,chip,chan);
759 //__________________________________________________________________
761 AliFMDParameters::Detector2Hardware(UShort_t det, Char_t ring,
762 UShort_t sec, UShort_t str,
763 UInt_t& ddl, UInt_t& addr) const
765 // Map detector index to hardware address
766 if (!fAltroMap) return kFALSE;
767 return fAltroMap->Detector2Hardware(det, ring, sec, str, ddl, addr);
771 //__________________________________________________________________
773 AliFMDParameters::GetEdepMip() const
775 // Get energy deposited by a MIP in the silicon sensors
777 AliFMDGeometry* fmd = AliFMDGeometry::Instance();
778 fEdepMip = (fkSiDeDxMip
779 * fmd->GetRing('I')->GetSiThickness()
780 * fmd->GetSiDensity());
789 //____________________________________________________________________