Double_t KVHarpeeIC::GetCalibE() { // Calculate energy in MeV from coder values. Only the 7 segments B // are used to calculate calibrated energy. // Returns 0 if calibration not ready for fired acquisition parameters // or if no acquisition parameter is fired. // (we require that at least one acquisition parameter has a value // greater than the current pedestal value). Bool_t ok = kTRUE; Double_t E = 0; TIter next(GetListOfCalibrators()); KVFunctionCal* cal = NULL; while (ok && (cal = (KVFunctionCal*)next())) { //Energy only coded in ACQParam with label 'A' if (cal->GetLabel()[0] != 'A') continue; // Int_t num = cal->GetUniqueID(); // if( ((num%100)/10) != 1 ) continue; if (cal->GetACQParam()->Fired("P")) { if (cal->GetStatus()) { E += cal->Compute(); } else ok = kFALSE; } } return (ok ? E : 0.); }
Double_t KVHarpeeIC::GetCalibE(Int_t seg_num) { // Calculate energy in MeV from the coder value of the acquisition parameter // with number 'seg_num' (i.e. segment number) and its corresponding calibrator. // Returns 0 if calibration not ready for the acquisition parameter. KVFunctionCal* cal = (KVFunctionCal*)GetListOfCalibrators()->FindObjectWithMethod(Form("%d", seg_num), "GetNumber"); if (cal && cal->GetStatus() && cal->GetACQParam()->Fired("P")) return cal->Compute(); return 0.; }
Double_t KVVAMOSDetector::GetCalibE() { // Calculate energy in MeV from coder values. // Returns 0 if calibration not ready or acquisition parameter not fired // (we require that the acquisition parameter has a value // greater than the current pedestal value) KVFunctionCal* cal = GetECalibrator(); if (cal && cal->GetStatus() && cal->GetACQParam()->Fired("P")) return cal->Compute(); return 0; }
Double_t KVVAMOSDetector::GetCalibT(const Char_t* type) { // Calculate calibrated time in ns of type 'type' (SED_HF, SI_HF, // SI_SED1, ...) for coder values. // Returns 0 if calibration not ready or time ACQ parameter not fired. // (we require that the acquisition parameter has a value // greater than the current pedestal value) KVFunctionCal* cal = GetTCalibrator(type); if (cal && cal->GetStatus() && cal->GetACQParam()->Fired("P")) return cal->Compute() + GetT0(type); return 0; }
Double_t KVHarpeeSi::GetCalibT(const Char_t *type){ // Calculate calibrated time in ns of type 'type' (SI_HF, SI_SED1, // SI_INDRA, SI_MCP, ...) for coder values. // Returns 0 if calibration not ready or time ACQ parameter not fired. // (we require that the acquisition parameter has a value // greater than the current pedestal value). // The reference time T0 used for the calibration is the one of the // fired silicon detector of Harpee which stopped the time. KVHarpeeSi *firedSi = GetFiredHarpeeSi(); if( !firedSi ) return 0; KVFunctionCal *cal = GetTCalibrator( type ); if( cal && cal->GetStatus() && cal->GetACQParam()->Fired("P")) return cal->Compute() + firedSi->GetT0( type ); return 0; }