QwBeamLine.cc

Go to the documentation of this file.

/*!
 * \file   QwBeamLine.cc
 * \brief  Beam line subsystem implementation for full beam instrumentation
 *
 * Subsystem managing the full beam instrumentation line: creation from map
 * files, geometry input, mock-data parameters, event processing, event cuts,
 * error propagation, publication of variables, and data encoding/decoding.
 * Documentation-only edits; runtime behavior unchanged.
 */
 
#include "QwBeamLine.h"
 
// System headers
#include <iostream>
#include <stdexcept>
 
#ifdef HAS_RNTUPLE_SUPPORT
// ROOT headers for RNTuple support
#include <ROOT/RNTupleModel.hxx>
#include <ROOT/RNTupleWriter.hxx>
#endif
 
// Qweak headers
#include "QwLog.h"
#ifdef __USE_DATABASE__
#include "QwParitySchema.h"
#include "QwParityDB.h"
#endif // __USE_DATABASE__
 
#include "QwPromptSummary.h"
 
// Forward declarations
#ifdef __USE_DATABASE__
class QwParityDB;
//class QwDBInterface;
#endif // __USE_DATABASE__
 
//*****************************************************************//
/** Parse and handle beamline-specific command-line options. */
void QwBeamLine::ProcessOptions(QwOptions &options){
      //Handle command line options
}
 
//*****************************************************************//
template <>
Int_t QwBeamLine::AddToElementList<VQwClock_ptr>(std::vector<VQwClock_ptr> &elementlist,
                   QwBeamDetectorID &detector_id)
{
  if (detector_id.fTypeID == kQwClock){
    VQwClock_ptr element(
               VQwClock::Create(GetName(),
                      detector_id.fdetectorname,
                      detector_id.fmoduletype) );
    elementlist.push_back(element);
  }
  detector_id.fIndex = elementlist.size()-1;
  return detector_id.fIndex;
}
//*****************************************************************//
template <>
Int_t QwBeamLine::AddToElementList<VQwBCM_ptr>(std::vector<VQwBCM_ptr> &elementlist,
                   QwBeamDetectorID &detector_id)
{
  if (detector_id.fTypeID == kQwBCM){
    VQwBCM_ptr element(
               VQwBCM::Create(GetName(),
                      detector_id.fdetectorname,
                      detector_id.fmoduletype) );
    elementlist.push_back(element);
  }
  detector_id.fIndex = elementlist.size()-1;
  return detector_id.fIndex;
}
//*****************************************************************//
template <>
Int_t QwBeamLine::AddToElementList<VQwBPM_ptr>(std::vector<VQwBPM_ptr> &elementlist,
                   QwBeamDetectorID &detector_id)
{
  if (detector_id.fTypeID == kQwBPMStripline){
    VQwBPM_ptr element(
    VQwBPM::CreateStripline(GetName(),
                detector_id.fdetectorname,
                detector_id.fmoduletype) );
    elementlist.push_back(element);
  }
  detector_id.fIndex = elementlist.size()-1;
  return detector_id.fIndex;
}
 
//*****************************************************************//
template <typename TT>
Int_t QwBeamLine::AddToElementList(std::vector<TT> &elementlist,
                   QwBeamDetectorID &detector_id)
{
  TT element(GetName(), detector_id.fdetectorname);
  elementlist.push_back(element);
  detector_id.fIndex = elementlist.size()-1;
  return detector_id.fIndex;
}
 
 
 
//*****************************************************************//
/**
 * Load and decode the beamline channel map, creating physical and combined
 * devices, recording detector IDs, and wiring publish lists and clocks.
 *
 * @param mapfile Path to the map file.
 * @return 0 on success.
 */
Int_t QwBeamLine::LoadChannelMap(TString mapfile)
{
  Bool_t ldebug=kFALSE;
 
  TString varname, varvalue;
 
  TString combotype, comboname, dev_name;
 
  Int_t   index = 0;
  Bool_t  combolistdecoded;
  Bool_t  deviceok;
 
  std::vector<TString>  fDeviceName;
  std::vector<TString>  fProperty;
  std::vector<TString>  fType;
  std::vector<Double_t> fQWeight;
  std::vector<Double_t> fXWeight;
  std::vector<Double_t> fYWeight;
  Double_t sumQweights = 0.0;
 
  std::vector<QwBeamDetectorID> clock_needed_list;
 
  QwParameterFile mapstr(mapfile.Data());  //Open the file
  fDetectorMaps.insert(mapstr.GetParamFileNameContents());
  mapstr.EnableGreediness();
  mapstr.SetCommentChars("!");
  mapstr.AddBreakpointKeyword("begin");
  mapstr.AddBreakpointKeyword("end");
 
  Int_t buffer_offset = 0; /*  Allow some extra words at the start of a bank.
                *  The buffer_offset value will be reset at the
                *  start of each ROC or bank declaration, so should
                *  be relisted for each bank.
                */
  while (mapstr.ReadNextLine() && mapstr.SkipSection("PUBLISH")) {
    RegisterRocBankMarker(mapstr);
    //  Remove the "vqwk_buffer_offset" and "scaler_buffer_offset"
    //  keywords from the parameter file's listing
    if (mapstr.ReturnValue("vqwk_buffer_offset",buffer_offset)) {
      QwDebug << "QwBeamLine::LoadChannelMap: "
          << "ROC " << fCurrentROC_ID
          << ", Bank " << fCurrentBank_ID
          << "; vqwk_buffer_offset:  " << buffer_offset << QwLog::endl;
    }
 
    if (mapstr.HasVariablePair("=",varname,varvalue)){ //  This is a declaration line.  Decode it.
      varname.ToLower();
 
      if (varname=="begin"){
 
          // Start to decode derived beamline devices combined+energy
          deviceok = kTRUE;
          combotype = varvalue;
          combolistdecoded = kFALSE;
 
          TString dettype;
 
          while(mapstr.ReadNextLine()&&!combolistdecoded) {
            if (mapstr.HasVariablePair("=",varname,varvalue)) {
              varname.ToLower();
              if (varname=="end"){
              // calculate the total weights of the charge
                sumQweights = 0.0;
                for(size_t i=0;i<fDeviceName.size();i++)
                    sumQweights+=fabs(fQWeight[i]);
                combolistdecoded = kTRUE;
                break;
              }
            }
 
            if (mapstr.PopValue("name",varvalue)) {
              comboname = varvalue;
            }
 
            dev_name = mapstr.GetTypedNextToken<TString>();
            dev_name.ToLower();
            dettype  = mapstr.GetTypedNextToken<TString>();
            dettype.ToLower();
 
            // Check to see if the device being read is a valid physical device.
            // If not, discard the combination.
            index=GetDetectorIndex(GetQwBeamInstrumentType(dettype),dev_name);
 
            if (index == -1) {
              QwError << "QwBeamLine::LoadChannelMap:  Unknown device: "
                  <<  dev_name <<" used in "<< comboname
                  <<". This combination  will not be decoded!"
                  <<  QwLog::endl;
              deviceok = kFALSE;
              combolistdecoded = kTRUE;
            }else{
              // Found the device
              // Add to the array of names
              fDeviceName.push_back(dev_name);
 
              // Read in the weights.
              // For combined bpms and combined bcms these are charge weights.
              // For the energy calculator these are the ratios of the transport matrix elements.
              fQWeight.push_back(mapstr.GetTypedNextToken<Double_t>());
 
              // For combined BPMs,in addition, there are weights for the X & Y positions.
              if(combotype == "combinedbpm"){
                fXWeight.push_back(mapstr.GetTypedNextToken<Double_t>());
                fYWeight.push_back(mapstr.GetTypedNextToken<Double_t>());
              }
 
              // For the enrgy calculator there are device type and the specified beam parameters.
              if(combotype == "energycalculator"){
                fType.push_back(dettype);
                fProperty.push_back(mapstr.GetTypedNextToken<TString>());
              }
            }
          }
 
    // Now create the combined device
    QwBeamDetectorID localComboID(-1, -1, comboname, combotype,
                      fBeamDetectorID.at(index).fmoduletype );
 
    localComboID.fdetectorname=comboname(0,comboname.Sizeof()-1);
    localComboID.fIndex = GetDetectorIndex(localComboID.fTypeID,localComboID.fdetectorname);
 
    if(localComboID.fTypeID==kQwUnknownDeviceType){
      QwError << "QwBeamLine::LoadChannelMap:  Unknown detector type: "
          << combotype <<", the detector "<<comboname<<" will not be decoded "
          << QwLog::endl;
      deviceok = kFALSE;
      continue;
    }
 
    if((localComboID.fIndex==-1) && deviceok) {
 
      // Decoding combined BCM array
      if (localComboID.fTypeID == kQwCombinedBCM){
 
        VQwBCM_ptr localbcmcombo(
        VQwBCM::CreateCombo(GetName(),
                    localComboID.fdetectorname,localComboID.fmoduletype));
        fBCMCombo.push_back(localbcmcombo);
 
        for(size_t i=0;i<fDeviceName.size();i++){
          index=GetDetectorIndex(GetQwBeamInstrumentType(dettype),fDeviceName[i]);
          fBCMCombo[fBCMCombo.size()-1].get()->SetBCMForCombo(fBCM.at(index).get(),
            fQWeight[i],sumQweights );
        }
        fDeviceName.clear();
        fQWeight.clear();
        localComboID.fIndex=fBCMCombo.size()-1;
      }
 
      // Decoding combined BPM array.
      if(localComboID.fTypeID== kQwCombinedBPM){
        VQwBPM_ptr localbpmcombo(
                VQwBPM::CreateCombo(GetName(),
                    localComboID.fdetectorname,localComboID.fmoduletype));
        fBPMCombo.push_back(localbpmcombo);
 
        for(size_t i=0;i<fDeviceName.size();i++){
          index=GetDetectorIndex(GetQwBeamInstrumentType(dettype),fDeviceName[i]);
          fBPMCombo[fBPMCombo.size()-1].get()->SetBPMForCombo(
              fStripline.at(index).get(), fQWeight[i],fXWeight[i],
              fYWeight[i],sumQweights  );
 
        }
        fDeviceName.clear();
        fQWeight.clear();
        fXWeight.clear();
        fYWeight.clear();
        localComboID.fIndex=fBPMCombo.size()-1;
      }
 
      // Decoding energy calculator.
      if(localComboID.fTypeID== kQwEnergyCalculator){
 
        QwEnergyCalculator localecalculator(GetName(), localComboID.fdetectorname);
        fECalculator.push_back(localecalculator);
 
        for(size_t i=0;i<fDeviceName.size();i++){
          index=GetDetectorIndex(GetQwBeamInstrumentType(fType[i]),fDeviceName[i]);
 
          if ( GetQwBeamInstrumentType(fType[i]) == kQwBPMStripline)
        fECalculator[fECalculator.size()-1].Set(fStripline.at(index).get(),fType[i],fProperty[i],fQWeight[i]);
 
          if ( GetQwBeamInstrumentType(fType[i]) == kQwCombinedBPM)
        fECalculator[fECalculator.size()-1].Set(fBPMCombo.at(index).get(),fType[i],fProperty[i],fQWeight[i]);
 
        }
 
        fDeviceName.clear();
        fQWeight.clear();
 
        fProperty.clear();
        fType.clear();
        localComboID.fIndex=fECalculator.size()-1;
      }
    }
    // Use only the combinations that are of known type and has known physical devices.
    if(deviceok)
      fBeamDetectorID.push_back(localComboID);
      }
 
      QwDebug << "At end of processing the combined device " << QwLog::endl;
 
    } else{
      // Start to decode the physical beamline devices
      QwBeamDetectorID localBeamDetectorID(GetSubbankIndex(), mapstr);
      Bool_t lineok = localBeamDetectorID.ReportInitErrors();
      if (! lineok) continue;
 
      localBeamDetectorID.fIndex=
    GetDetectorIndex(localBeamDetectorID.fTypeID,
             localBeamDetectorID.fdetectorname);
 
      if(localBeamDetectorID.fIndex==-1){
    Int_t index;
    VQwDataElement* local_element = NULL;
 
    if(localBeamDetectorID.fTypeID == kQwHaloMonitor){
      index = AddToElementList(fHaloMonitor, localBeamDetectorID);
      local_element = &(fHaloMonitor.at(index));
    }
    if(localBeamDetectorID.fTypeID==kQwBPMCavity){
      index = AddToElementList(fCavity, localBeamDetectorID);
      local_element = &(fCavity.at(index));
    }
    if(localBeamDetectorID.fTypeID== kQwBPMStripline){
      index = AddToElementList(fStripline, localBeamDetectorID);
      local_element = fStripline.at(index).get();
    }
    if(localBeamDetectorID.fTypeID== kQwBCM){
      index = AddToElementList(fBCM, localBeamDetectorID);
      local_element = fBCM.at(index).get();
    }
    if(localBeamDetectorID.fTypeID== kQwClock ) {
      index = AddToElementList(fClock, localBeamDetectorID);
      local_element = fClock.at(index).get();
    }
    if(localBeamDetectorID.fTypeID== kQwQPD){
      index = AddToElementList(fQPD, localBeamDetectorID);
      local_element = &(fQPD.at(index));
    }
    if(localBeamDetectorID.fTypeID== kQwLinearArray){
      index = AddToElementList(fLinearArray, localBeamDetectorID);
      local_element = &(fLinearArray.at(index));
    }
 
    local_element->LoadChannelParameters(mapstr);
    if (local_element->NeedsExternalClock()){
      QwDebug << "Try to push device "
          << local_element->GetElementName()
          << " (address=" << std::hex << local_element << std::dec
          << ") onto the clock_needed_list"
          << QwLog::endl;
      clock_needed_list.push_back(localBeamDetectorID);
    }
      }
 
 
 
      fBeamDetectorID.push_back(localBeamDetectorID);
    }
  }
 
  // Now load the variables to publish
  mapstr.RewindToFileStart();
  std::unique_ptr<QwParameterFile> section;
  std::vector<TString> publishinfo;
  while ((section=mapstr.ReadNextSection(varvalue))) {
    if (varvalue == "PUBLISH") {
      fPublishList.clear();
      while (section->ReadNextLine()) {
        section->TrimComment(); // Remove everything after a comment character
        section->TrimWhitespace(); // Get rid of leading and trailing spaces
        for (int ii = 0; ii < 4; ii++) {
          varvalue = section->GetTypedNextToken<TString>();
          if (varvalue.Length()) {
            publishinfo.push_back(varvalue);
          }
        }
        if (publishinfo.size() == 4)
          fPublishList.push_back(publishinfo);
        publishinfo.clear();
      }
    }
  }
  // Print list of variables to publish
  if (fPublishList.size()>0){
    QwMessage << "Variables to publish:" << QwLog::endl;
    for (size_t jj = 0; jj < fPublishList.size(); jj++)
      QwMessage << fPublishList.at(jj).at(0) << " " << fPublishList.at(jj).at(1) << " "
        << fPublishList.at(jj).at(2) << " " << fPublishList.at(jj).at(3) << QwLog::endl;
  }
 
  if(ldebug){
    std::cout<<"QwBeamLine::Done with Load map channel \n";
    for(size_t i=0;i<fBeamDetectorID.size();i++)
      fBeamDetectorID[i].Print();
  }
 
  // Now propagate clock pointers to those channels that need it
  index = 0;
  VQwDataElement* local_element;
  std::string clockname;
  for (size_t i=0; i<clock_needed_list.size();i++ ) {
    local_element = GetElement(clock_needed_list[i]);
    clockname = local_element->GetExternalClockName();
    if (clockname.empty()){
    QwWarning << "QwBeamLine::LoadChannelMap  "
          << "Device, " << local_element->GetElementName()
          << " needs a reference clock, but the reference clock name is empty"
          << QwLog::endl;
    } else {
      index = GetDetectorIndex(kQwClock,clockname);
      if( index >= 0 ){
    QwMessage << "QwBeamLine::LoadChannelMap  "
          << "Setting " << fClock.at(index).get()->GetElementName()
          << " as the reference clock for channel "
          << local_element->GetElementName()
          << QwLog::endl;
    local_element->SetExternalClockPtr(fClock.at(index).get()->GetTime());
      } else {
    QwWarning << "QwBeamLine::LoadChannelMap  "
          << "Cannot find clock, " << local_element->GetExternalClockName()
          << ", needed by device, " << local_element->GetElementName()
          << QwLog::endl;
      }
    }
  }
  ldebug=kFALSE;
 
  mapstr.Close(); // Close the file (ifstream)
  return 0;
}
 
 
//*****************************************************************//
/**
 * Decode a single event-cuts line for devices and subchannels, routing the
 * resulting thresholds and error flags to the proper element.
 */
void QwBeamLine::LoadEventCuts_Line(QwParameterFile &mapstr, TString &varvalue, Int_t &eventcut_flag) {
  TString device_type = mapstr.GetTypedNextToken<TString>();
  device_type.ToLower();
  TString device_name = mapstr.GetTypedNextToken<TString>();
  device_name.ToLower();
 
  Int_t det_index = GetDetectorIndex(GetQwBeamInstrumentType(device_type),device_name);
  if (det_index == -1) {
    QwWarning << " Device not found " << device_name << " of type " << device_type << QwLog::endl;
    //continue;
  }
 
  TString channel_name;
 
  if(device_type == GetQwBeamInstrumentTypeName(kQwBPMStripline) ||
     device_type == GetQwBeamInstrumentTypeName(kQwQPD) ||
     device_type == GetQwBeamInstrumentTypeName(kQwLinearArray)||
     device_type ==  GetQwBeamInstrumentTypeName(kQwBPMCavity)||
     device_type == GetQwBeamInstrumentTypeName(kQwCombinedBPM)){
 
    channel_name =  mapstr.GetTypedNextToken<TString>();
      channel_name.ToLower();
  }
 
    Double_t LLX = mapstr.GetTypedNextToken<Double_t>();   //lower limit for BCM value
    Double_t ULX = mapstr.GetTypedNextToken<Double_t>();   //upper limit for BCM value
    varvalue = mapstr.GetTypedNextToken<TString>();//global/local
  varvalue.ToLower();
 
  Double_t stabilitycut = mapstr.GetTypedNextToken<Double_t>();
  Double_t burplevel = mapstr.GetTypedNextToken<Double_t>();
 
 
  if (device_type == GetQwBeamInstrumentTypeName(kQwBCM)){
      QwMessage<<"QwBeamLine Error Code passing to QwBCM "<<GetGlobalErrorFlag(varvalue,eventcut_flag,stabilitycut)<<", burp = "<<burplevel<<QwLog::endl;
    fBCM[det_index].get()->SetSingleEventCuts(GetGlobalErrorFlag(varvalue,eventcut_flag,stabilitycut),LLX,ULX,stabilitycut,burplevel);//(fBCMEventCuts);
 
  } else if (device_type == GetQwBeamInstrumentTypeName(kQwHaloMonitor)){
      QwMessage<<"QwBeamLine Error Code passing to QwHaloMonitor "<<GetGlobalErrorFlag(varvalue,eventcut_flag,stabilitycut)<<QwLog::endl;
    fHaloMonitor[det_index].SetSingleEventCuts(GetGlobalErrorFlag(varvalue,eventcut_flag,stabilitycut),LLX,ULX,stabilitycut,burplevel);//(fBCMEventCuts);
 
  } else if (device_type == GetQwBeamInstrumentTypeName(kQwEnergyCalculator)){
      QwMessage<<"QwBeamLine Error Code passing to QwEnergyCalculator "<<GetGlobalErrorFlag(varvalue,eventcut_flag,stabilitycut)<<QwLog::endl;
    fECalculator[det_index].SetSingleEventCuts(GetGlobalErrorFlag(varvalue,eventcut_flag,stabilitycut),LLX,ULX,stabilitycut,burplevel);//(fEnergyEventCuts);
 
  } else if (device_type == GetQwBeamInstrumentTypeName(kQwBPMStripline)){
      QwMessage<<"QwBeamLine:QwBPMStripline "<<channel_name<<" "<<varvalue<<" "<<stabilitycut<<QwLog::endl;
      //QwMessage<<"QwBeamLine Error Code passing to QwBPMStripline "<<GetGlobalErrorFlag(varvalue,eventcut_flag,stabilitycut)<<" stability  "<<stabilitycut <<QwLog::endl;
      fStripline[det_index].get()->SetSingleEventCuts(channel_name, GetGlobalErrorFlag(varvalue,eventcut_flag,stabilitycut), LLX, ULX, stabilitycut, burplevel);
 
  } else if (device_type == GetQwBeamInstrumentTypeName(kQwQPD)){
      QwMessage<<"QwBeamLine Error Code passing to QwQPD "<<GetGlobalErrorFlag(varvalue,eventcut_flag,stabilitycut)<<QwLog::endl;
      fQPD[det_index].SetSingleEventCuts(channel_name, GetGlobalErrorFlag(varvalue,eventcut_flag,stabilitycut), LLX, ULX, stabilitycut, burplevel);
 
  } else if (device_type == GetQwBeamInstrumentTypeName(kQwLinearArray)){
      QwMessage<<"QwBeamLine Error Code passing to QwLinearArray "<<GetGlobalErrorFlag(varvalue,eventcut_flag,stabilitycut)<<QwLog::endl;
      fLinearArray[det_index].SetSingleEventCuts(channel_name, GetGlobalErrorFlag(varvalue,eventcut_flag,stabilitycut), LLX, ULX, stabilitycut, burplevel);
 
  } else if (device_type ==  GetQwBeamInstrumentTypeName(kQwBPMCavity)){
    varvalue.ToLower();
    QwMessage<<"QwBeamLine Error Code passing to QwBPMCavity "<<GetGlobalErrorFlag(varvalue,eventcut_flag,stabilitycut)<<" "<<det_index<<QwLog::endl;
    fCavity[det_index].SetSingleEventCuts(channel_name, GetGlobalErrorFlag(varvalue,eventcut_flag,stabilitycut), LLX, ULX, stabilitycut, burplevel);
 
  } else if (device_type == GetQwBeamInstrumentTypeName(kQwCombinedBCM)){
    varvalue.ToLower();
    QwMessage<<"QwBeamLine Error Code passing to QwCombinedBCM "<<GetGlobalErrorFlag(varvalue,eventcut_flag,stabilitycut)<<QwLog::endl;
    fBCMCombo[det_index].get()->PrintInfo();
    fBCMCombo[det_index].get()->SetSingleEventCuts(GetGlobalErrorFlag(varvalue,eventcut_flag,stabilitycut), LLX, ULX, stabilitycut, burplevel);
 
  } else if (device_type == GetQwBeamInstrumentTypeName(kQwCombinedBPM)){
    varvalue.ToLower();
    QwMessage<<"QwBeamLine Error Code passing to QwCombinedBPM "<<GetGlobalErrorFlag(varvalue,eventcut_flag,stabilitycut)<<QwLog::endl;
    fBPMCombo[det_index].get()->SetSingleEventCuts(channel_name, GetGlobalErrorFlag(varvalue,eventcut_flag,stabilitycut), LLX, ULX, stabilitycut, burplevel);
  }
}
 
/** Apply the event-cut mode bitmask to all managed devices. */
void QwBeamLine::LoadEventCuts_Fin(Int_t &eventcut_flag) {
  for (size_t i=0;i<fStripline.size();i++)
    fStripline[i].get()->SetEventCutMode(eventcut_flag);
 
  for (size_t i=0;i<fQPD.size();i++)
    fQPD[i].SetEventCutMode(eventcut_flag);
 
  for (size_t i=0;i<fLinearArray.size();i++)
    fLinearArray[i].SetEventCutMode(eventcut_flag);
 
  for (size_t i=0;i<fCavity.size();i++)
    fCavity[i].SetEventCutMode(eventcut_flag);
 
  for (size_t i=0;i<fBCM.size();i++)
    fBCM[i].get()->SetEventCutMode(eventcut_flag);
 
  for (size_t i=0;i<fClock.size();i++)
    fClock[i].get()->SetEventCutMode(eventcut_flag);
 
  for (size_t i=0;i<fHaloMonitor.size();i++)
    fHaloMonitor[i].SetEventCutMode(eventcut_flag);
 
  for (size_t i=0;i<fBCMCombo.size();i++)
    fBCMCombo[i].get()->SetEventCutMode(eventcut_flag);
 
  for (size_t i=0;i<fBPMCombo.size();i++)
    fBPMCombo[i].get()->SetEventCutMode(eventcut_flag);
 
  for (size_t i=0;i<fECalculator.size();i++)
    fECalculator[i].SetEventCutMode(eventcut_flag);
 
  fQwBeamLineErrorCount=0; //set the error counter to zero
}
 
//*****************************************************************//
/**
 * Load geometry and electronics calibration parameters from a map file and
 * apply them to the proper devices (offsets, rotation, gains).
 */
Int_t QwBeamLine::LoadGeometryDefinition(TString mapfile){
  Bool_t ldebug=kFALSE;
  TString varname, varvalue;
  Int_t lineread=1;
  Int_t index = 0;
  TString  devname,devtype;
  TString  melement;
  Double_t devOffsetX = 0,devOffsetY = 0, devOffsetZ = 0;
  Double_t devSENfactor = 0, devAlphaX = 0, devAlphaY = 0;
  TString   rotation_stat;
  VQwBPM * bpm = NULL;
 
  if(ldebug)std::cout<<"QwBeamLine::LoadGeometryParameters("<< mapfile<<")\n";
 
  QwParameterFile mapstr(mapfile.Data());  //Open the file
  fDetectorMaps.insert(mapstr.GetParamFileNameContents());
  while (mapstr.ReadNextLine()){
    lineread+=1;
    if(ldebug)std::cout<<" line read so far ="<<lineread<<"\n";
    mapstr.TrimComment('!');
    mapstr.TrimWhitespace();
 
    if (mapstr.LineIsEmpty())  continue;
 
    Bool_t notfound=kTRUE;
 
    devtype = mapstr.GetTypedNextToken<TString>();
    devtype.ToLower();
    devtype.Remove(TString::kBoth,' ');
    devname = mapstr.GetTypedNextToken<TString>();
    devname.ToLower();
    devname.Remove(TString::kBoth,' ');
 
    if (GetQwBeamInstrumentType(devtype)==kQwUnknownDeviceType){
      QwError << "Error! Unknown detector type '"<<devtype
          << "' in Geometry file!"<<QwLog::endl;
      /*If the device type is unknown there is no point in going through the rest of the specs for that device*/
      continue;
    }
 
    index=GetDetectorIndex(GetQwBeamInstrumentType(devtype),devname);
    if( index<0 ) {
      /*Detector name isn't recognized. Ignore it!*/
      QwWarning << "Unrecognized detector name '" << devname
        << "' in Geometry file.  This may not be a problem, "
        << "if we're using a reduced channel map."
        << QwLog::endl;
    }
    else {
      devOffsetX   = mapstr.GetTypedNextToken<Double_t>(); // X offset
      devOffsetY   = mapstr.GetTypedNextToken<Double_t>(); // Y offset
      devOffsetZ   = mapstr.GetTypedNextToken<Double_t>(); // Z offset
      devSENfactor = mapstr.GetTypedNextToken<Double_t>(); // sensitivity scaling factor
      devAlphaX    = mapstr.GetTypedNextToken<Double_t>(); // alpha X
      devAlphaY    = mapstr.GetTypedNextToken<Double_t>(); // alpha Y
 
 
      /*If the device is a bpm stripline, assign the rotations and gains*/
      if(GetQwBeamInstrumentType(devtype)==kQwBPMStripline){
    bpm = fStripline.at(index).get();
    AssignGeometry(&mapstr,bpm);
      }
 
 
      if(ldebug==1){
    std::cout<<"####################\n";
    std::cout<<"! device type, device_name, Xoffset, Yoffset, Zoffset, BSEN scaling factor, AlpaX, AlpaY\n"<<std::endl;
    std::cout<<GetQwBeamInstrumentType(devtype)<<" / "
         <<devname    <<" / "
         <<devOffsetX <<" / "
         <<devOffsetY <<" / "
         <<devOffsetZ <<" / "
         <<devSENfactor <<" / "
         <<devAlphaX <<" / "
         <<devAlphaY <<" / "
         <<std::endl;
      }
 
 
      while(notfound){
    if(GetQwBeamInstrumentType(devtype)==kQwBPMStripline){
      //Load bpm offsets
      if(index == -1){
        QwWarning << "QwBeamLine::LoadGeometryDefinition:  Unknown bpm in qweak_beamline_geometry.map: "
              <<devname
              <<QwLog::endl;
        notfound=kFALSE;
        continue;
      }
 
      TString localname = fStripline.at(index).get()->GetElementName();
      localname.ToLower();
      if(ldebug)  std::cout<<"element name =="<<localname
                   <<"== to be compared to =="<<devname<<"== \n";
 
      if(localname==devname){
        if(ldebug) std::cout<<" I found the bpm !\n";
        bpm->GetSurveyOffsets(devOffsetX,devOffsetY,devOffsetZ);
        bpm->GetElectronicFactors(devSENfactor,devAlphaX, devAlphaY);
 
        // If nothing is specified, a default rotation of 45 degrees is implied.
        notfound=kFALSE;
      }
    }
    else if (GetQwBeamInstrumentType(devtype)==kQwCombinedBPM){
      //Load combined bpm offsets which are, ofcourse, target position in the beamline
      if(index == -1){
        QwError << "QwBeamLine::LoadGeometryDefinition:  Unknown combined bpm in qweak_beamline_geometry.map: "
            <<devname<<" Check the combined bpm names!\n "
            << QwLog::endl;
        notfound=kFALSE;
        continue;
      }
 
      TString localname = fBPMCombo.at(index).get()->GetElementName();
      localname.ToLower();
      if(ldebug)
        std::cout<<"element name =="<<localname<<"== to be compared to =="<<devname<<"== \n";
 
      if(localname==devname){
        if(ldebug) std::cout<<" I found the combinedbpm !\n";
        fBPMCombo.at(index).get()->GetSurveyOffsets(devOffsetX,devOffsetY,devOffsetZ);
        notfound=kFALSE;
      }
    }
    else if(GetQwBeamInstrumentType(devtype)==kQwBPMCavity){
      //Load cavity bpm offsets
      if(index == -1){
        QwError << "QwBeamLine::LoadGeometryDefinition:  Unknown bpm : "
            <<devname<<" will not be assigned with geometry parameters. \n"
            <<QwLog::endl;
        notfound=kFALSE;
        continue;
      }
      TString localname = fCavity.at(index).GetElementName();
      localname.ToLower();
      if(ldebug)  std::cout<<"element name =="<<localname
                   <<"== to be compared to =="<<devname<<"== \n";
 
      if(localname==devname){
        if(ldebug) std::cout<<" I found the cavity bpm !\n";
        fCavity.at(index).GetSurveyOffsets(devOffsetX,devOffsetY,devOffsetZ);
        notfound=kFALSE;
      }
    }
    else if(GetQwBeamInstrumentType(devtype)==kQwQPD){
      //Load QPD calibration factors
      if(index == -1){
        QwError << "QwBeamLine::LoadGeometryDefinition:  Unknown QPD : "
            <<devname<<" will not be assigned with calibration factors. \n"
            <<QwLog::endl;
        notfound=kFALSE;
        continue;
      }
      TString localname = fQPD.at(index).GetElementName();
      localname.ToLower();
      if(ldebug)  std::cout<<"element name =="<<localname
                   <<"== to be compared to =="<<devname<<"== \n";
 
      if(localname==devname){
        if(ldebug) std::cout<<"I found the QPD !\n";
        fQPD.at(index).GetCalibrationFactors(devAlphaX, devAlphaY);
        notfound=kFALSE;
      }
    }
    else QwError << "QwBeamLine::LoadGeometryDefinition: Unknown device type : "<< devtype <<
                  ". Are you sure we have this in the beamline? I am skipping this."<< QwLog::endl;
      }
    }
  }
 
  if(ldebug) std::cout << "line read in the geometry file = " << lineread<< " \n";
 
  ldebug=kFALSE;
  mapstr.Close(); // Close the file (ifstream)
  return 0;
 
}
 
//--------------------------------------------------------------------------------------------------------
/** Load per-device mock-data generation parameters from a file. */
void QwBeamLine::LoadMockDataParameters(TString mapfile) {
 
  Bool_t   ldebug=kFALSE;
  TString  varname, varvalue;
  Int_t    lineread=1;
  Int_t    index = 0;
  TString  devname, devtype;
  TString  melement;
 
  if(ldebug) std::cout << "QwBeamLine::LoadMockDataParameters(" << mapfile << ") \n" << std::endl;
 
  QwParameterFile mapstr(mapfile.Data());  //Open the file
  fDetectorMaps.insert(mapstr.GetParamFileNameContents());
 
  while (mapstr.ReadNextLine()) {
    lineread+=1;
    //    std::cerr << "Line:  " << mapstr.GetLine() << std::endl;
    if(ldebug) std::cout << "Line read so far = " << lineread << "\n" << std::endl;
    mapstr.TrimComment('!');
    mapstr.TrimWhitespace();
 
    if(mapstr.LineIsEmpty()) continue;
 
    devtype = mapstr.GetTypedNextToken<TString>();
    devtype.ToLower();
    devtype.Remove(TString::kBoth,' ');
    devname = mapstr.GetTypedNextToken<TString>();
    devname.ToLower();
    devname.Remove(TString::kBoth,' ');
 
    if (GetQwBeamInstrumentType(devtype)==kQwUnknownDeviceType) {
      /*If the device type is unknown there is no point in going through the rest of the specs for that device*/
      QwError << "Error! Unknown detector type '" << devtype << "' in MockDataParameters file!" << QwLog::endl;
      continue;
    }
    index = GetDetectorIndex(GetQwBeamInstrumentType(devtype),devname);
    if (index<0) {
      /*Detector name isn't recognized. Ignore it!*/
      QwWarning << "Unrecognized detector name '" << devname << "' in MockDataParameters file." << QwLog::endl;
      continue;
    }
    //  The device should process the reminder of this line.
    GetElement(GetQwBeamInstrumentType(devtype),index)->LoadMockDataParameters(mapstr);
  }
}
 
//--------------------------------------------------------------------------------------------------------
 
/** Parse rotation and gain tokens for a BPM and apply to the instance. */
void QwBeamLine::AssignGeometry(QwParameterFile* mapstr, VQwBPM * bpm)
{
 
  Bool_t ldebug = kFALSE;
 
  TString token = "0";
  TString   angle,xgain,ygain;
  Double_t  rotation_angle = 0;
 
  while(token!=""){
    token= mapstr->GetTypedNextToken<TString>();
    token.Remove(TString::kBoth,'\0');
 
    if(token.Contains("unrotated")){
      if(ldebug) std::cout<<" unrotated "<<std::endl;
      bpm->SetRotationOff();
    }
    else if(token.Contains("rotation")){
      // If the status is 'rotated'
 
      // If a specific rotation angle is given read that
      if(token.Contains("=")){
    angle = token.Remove(0,9);
    rotation_angle = atof(angle);
    if(ldebug) std::cout<<"Rotation angle = "<<rotation_angle<<std::endl;
    bpm->SetRotation(rotation_angle);
      }
    }
    // If nothing is specified for rotation, a default rotation of 45 degrees is implied.
 
    if(token.Contains("xgain")){
      xgain = token.Remove(0,6);
      if(ldebug) std::cout<<" xgain ="<<xgain<<std::endl;
      bpm->SetGains("X", atof(xgain));
    }
 
    if(token.Contains("ygain")){
      ygain = token.Remove(0,6);
      if(ldebug) std::cout<<" ygain ="<<ygain<<std::endl;
      bpm->SetGains("Y", atof(ygain));
    }
  }
}
 
//*****************************************************************//
/**
 * Load pedestals and calibration factors per subelement/device.
 *
 * @param pedestalfile Path to the input parameter file.
 * @return 0 on success.
 */
Int_t QwBeamLine::LoadInputParameters(TString pedestalfile)
{
  Bool_t ldebug=kFALSE;
 
  Int_t lineread=1;
 
  if(ldebug)std::cout<<"QwBeamLine::LoadInputParameters("<< pedestalfile<<")\n";
 
  QwParameterFile mapstr(pedestalfile.Data());  //Open the file
  fDetectorMaps.insert(mapstr.GetParamFileNameContents());
 
   while (mapstr.ReadNextLine())
    {
      lineread+=1;
      if(ldebug)std::cout<<" line read so far ="<<lineread<<"\n";
      mapstr.TrimComment('!');   // Remove everything after a '!' character.
      mapstr.TrimWhitespace();   // Get rid of leading and trailing spaces.
      if (mapstr.LineIsEmpty())  continue;
      else
    {
      TString varname = mapstr.GetTypedNextToken<TString>();   //name of the channel
      varname.ToLower();
      varname.Remove(TString::kBoth,' ');
      Double_t varped = mapstr.GetTypedNextToken<Double_t>(); // value of the pedestal
      Double_t varcal = mapstr.GetTypedNextToken<Double_t>(); // value of the calibration factor
      /*Double_t varweight = */ mapstr.GetTypedNextToken<Double_t>(); // value of the statistical weight
 
      //if(ldebug) std::cout<<"inputs for channel "<<varname
      //          <<": ped="<<varped<<": cal="<<varcal<<": weight="<<varweight<<"\n";
      Bool_t notfound=kTRUE;
 
      if(notfound) {
        for(size_t i=0;i<fStripline.size();i++)
          {
        for(int j=0;j<4;j++)
          {
            TString localname = fStripline[i].get()->GetSubElementName(j);
            localname.ToLower();
            if(ldebug)  std::cout<<"Stripline element name =="<<localname
                     <<"== to be compared to =="<<varname<<"== \n";
            if(localname==varname)
              {
            if(ldebug) std::cout<<" I found it !\n";
            fStripline[i].get()->SetSubElementPedestal(j,varped);
            fStripline[i].get()->SetSubElementCalibrationFactor(j,varcal);
            notfound=kFALSE;
            j=5;
            i=fStripline.size()+1;
              }
          }
          }
        for(size_t i=0;i<fQPD.size();i++)
          {
        for(int j=0;j<4;j++)
          {
            TString localname = fQPD[i].GetSubElementName(j);
            localname.ToLower();
            if(ldebug)  std::cout<<"QPD element name =="<<localname
                     <<"== to be compared to =="<<varname<<"== \n";
            if(localname==varname)
              {
            if(ldebug) std::cout<<" I found it !\n";
            fQPD[i].SetSubElementPedestal(j,varped);
            fQPD[i].SetSubElementCalibrationFactor(j,varcal);
            notfound=kFALSE;
            j=5;
            i=fQPD.size()+1;
              }
          }
          }
        for(size_t i=0;i<fLinearArray.size();i++)
          {
        for(int j=0;j<8;j++)
          {
            TString localname = fLinearArray[i].GetSubElementName(j);
            localname.ToLower();
            if(ldebug)  std::cout<<"LinearArray element name =="<<localname
                     <<"== to be compared to =="<<varname<<"== \n";
            if(localname==varname)
              {
            if(ldebug) std::cout<<" I found it !\n";
            fLinearArray[i].SetSubElementPedestal(j,varped);
            fLinearArray[i].SetSubElementCalibrationFactor(j,varcal);
            notfound=kFALSE;
            j=9;
            i=fLinearArray.size()+1;
              }
          }
          }
        for(size_t i=0;i<fCavity.size();i++)
          {
        for(size_t j=0;j<QwBPMCavity::kNumElements;j++)
          {
            TString localname = fCavity[i].GetSubElementName(j);
            localname.ToLower();
            if(ldebug)  std::cout<<"Cavity element name =="<<localname
                     <<"== to be compared to =="<<varname<<"== \n";
            if(localname==varname)
              {
            if(ldebug) std::cout<<" I found it !\n";
            fCavity[i].SetSubElementPedestal(j,varped);
            fCavity[i].SetSubElementCalibrationFactor(j,varcal);
            notfound=kFALSE;
            j=3;
            i=fCavity.size()+1;
              }
          }
          }
        for(size_t i=0;i<fBCM.size();i++) {
          if(fBCM[i].get()->GetElementName()==varname)
        {
          fBCM[i].get()->SetPedestal(varped);
          fBCM[i].get()->SetCalibrationFactor(varcal);
          i=fBCM.size()+1;
          notfound=kFALSE;
          i=fBCM.size()+1;
        }
            }
        for(size_t i=0;i<fClock.size();i++) {
          if(fClock[i].get()->GetElementName()==varname)
        {
          fClock[i].get()->SetPedestal(varped);
          fClock[i].get()->SetCalibrationFactor(varcal);
          i=fClock.size()+1;
          notfound=kFALSE;
          i=fClock.size()+1;
        }
            }
        for(size_t i=0;i<fHaloMonitor.size();i++) {
          if(fHaloMonitor[i].GetElementName()==varname)
        {
          std::cout<<varname<<" I found it ! "<<varcal<<" ped. "<<varped<<"\n";
          fHaloMonitor[i].SetPedestal(varped);
          fHaloMonitor[i].SetCalibrationFactor(varcal);
          i=fHaloMonitor.size()+1;
          notfound=kFALSE;
          i=fHaloMonitor.size()+1;
        }
            }
          }
    }
    }
  if(ldebug) std::cout<<" line read in the pedestal + cal file ="<<lineread<<" \n";
 
  ldebug=kFALSE;
  mapstr.Close(); // Close the file (ifstream)
  return 0;
}
 
 
//*****************************************************************//
/** Randomize event data for all managed devices for mock runs. */
void QwBeamLine::RandomizeEventData(int helicity, double time)
{
  // Randomize all QwBPMStripline buffers
  for (size_t i = 0; i < fStripline.size(); i++)
  {
    fStripline[i].get()->RandomizeEventData(helicity, time);
     // fStripline[i].get()->PrintInfo();
  }
 
  for (size_t i = 0; i < fCavity.size(); i++)
    fCavity[i].RandomizeEventData(helicity, time);
 
  // Randomize all QwBCM buffers
  for (size_t i = 0; i < fBCM.size(); i++)
  {
    fBCM[i].get()->RandomizeEventData(helicity, time);
    //    fBCM[i].get()->PrintInfo();
  }
 
  // Randomize all QwHaloMonitor buffers
  //for (size_t i = 0; i < fHaloMonitor.size(); i++)
    //fHaloMonitor[i].RandomizeEventData(helicity, time);
 
//-------------------------------------------------------------------------------------------
  for(size_t i=0;i<fBCMCombo.size();i++){
    fBCMCombo[i].get()->RandomizeEventData(helicity, time);
    for (size_t j=0; j<fBCMCombo[i].get()->GetNumberOfElements(); j++){
      VQwBCM * bcm = GetBCM(fBCMCombo[i].get()->GetSubElementName(j));
      if (bcm){
        fBCMCombo[i].get()->GetProjectedCharge(bcm);
      }
    }
  }
//-------------------------------------------------------------------------------------------
  for (size_t i=0; i<fBPMCombo.size(); i++){
    fBPMCombo[i].get()->RandomizeEventData(helicity, time);
    //fBPMCombo[i].get()->PrintValue();
    for (size_t j=0; j<fBPMCombo[i].get()->GetNumberOfElements(); j++){
      VQwBPM * bpm = GetBPMStripline(fBPMCombo[i].get()->GetSubElementName(j));
      if (bpm){
        fBPMCombo[i].get()->GetProjectedPosition(bpm);
//  print the bpm device name, and get the x and y and z? values and print them
//      std::cout << "bpm " << bpm->GetElementName() << std::endl;
//      std::cout << "xpos= " << bpm->GetPosition(VQwBPM::kXAxis) << std::endl;
//      std::cout << "ypos= " << bpm->GetPosition(VQwBPM::kYAxis) << std::endl;
//      std::cout << "zpos= " << bpm->GetPosition(VQwBPM::kZAxsi) << std::endl;
//        bpm->PrintInfo();
        //  Call the new function in stripline class to fill all internal variables from the fAbsPos for the bpm object
      }
    }
  }
//-------------------------------------------------------------------------------------------
  for (size_t i=0;i<fECalculator.size();i++)
  {
   fECalculator[i].RandomizeEventData(helicity, time);
   //fECalculator[i].PrintValue();
   for (size_t j=0; j<fECalculator[i].GetNumberOfElements(); j++)
    {
      // std::cout << "In loop over ec subelements; device name = " << fECalculator[i].GetSubElementName(j) << std::endl;
      VQwBPM * bpm = GetBPMStripline(fECalculator[i].GetSubElementName(j));
      if (bpm)
      {
        fECalculator[i].GetProjectedPosition(bpm);
      }
    }
  }
}
 
 
//*****************************************************************//
/**
 * Serialize active element buffers into a CODA-style bank/subbank layout and
 * append to the provided buffer.
 */
void QwBeamLine::EncodeEventData(std::vector<UInt_t> &buffer)
{
  std::vector<UInt_t> elements;
  elements.clear();
 
  // Get all buffers in the order they are defined in the map file
  for (size_t i = 0; i < fBeamDetectorID.size(); i++) {
    // This is a QwBCM
    if (fBeamDetectorID.at(i).fTypeID == kQwBCM){
      fBCM[fBeamDetectorID.at(i).fIndex].get()->EncodeEventData(elements);
      //std::cout << "" << fBCM[fBeamDetectorID.at(i).fIndex].get()->GetElementName() << std::endl;
    }
    // This is a QwBPMStripline (which has 4 entries, only process the first one)
    if (fBeamDetectorID.at(i).fTypeID == kQwBPMStripline
     && fBeamDetectorID.at(i).fSubelement == 0){
      fStripline[fBeamDetectorID.at(i).fIndex].get()->EncodeEventData(elements);
      //  Print the HWsum absolute position values for the BPM
      //fStripline[fBeamDetectorID.at(i).fIndex].get()->PrintValue();
    }
 
    //  If this is a combined BPM, let's try to print the position and angle HWsum values
  }
 
  // If there is element data, generate the subbank header
  std::vector<UInt_t> subbankheader;
  std::vector<UInt_t> rocheader;
  if (elements.size() > 0) {
 
    // Form CODA subbank header
    subbankheader.clear();
    subbankheader.push_back(elements.size() + 1);   // subbank size
    subbankheader.push_back((fCurrentBank_ID << 16) | (0x01 << 8) | (1 & 0xff));
        // subbank tag | subbank type | event number
 
    // Form CODA bank/roc header
    rocheader.clear();
    rocheader.push_back(subbankheader.size() + elements.size() + 1);    // bank/roc size
    rocheader.push_back((fCurrentROC_ID << 16) | (0x10 << 8) | (1 & 0xff));
        // bank tag == ROC | bank type | event number
 
    // Add bank header, subbank header and element data to output buffer
    buffer.insert(buffer.end(), rocheader.begin(), rocheader.end());
    buffer.insert(buffer.end(), subbankheader.begin(), subbankheader.end());
    buffer.insert(buffer.end(), elements.begin(), elements.end());
  }
}
 
//*****************************************************************//
/**
 * Route a raw ROC/bank buffer to the correct devices and subelements.
 *
 * @return 0 on success.
 */
/*!
 * \brief Route raw ROC/bank buffer data to the correct devices and subelements.
 * \param roc_id ReadOut Controller identifier.
 * \param bank_id Data bank identifier within the ROC.
 * \param buffer Pointer to raw data buffer from DAQ.
 * \param num_words Number of 32-bit words in the buffer.
 * \return 0 on success.
 *
 * This is a critical event processing function (~130 lines) that routes raw data
 * from the data acquisition system to the appropriate beam line devices. It:
 *
 * - Maps ROC/bank combinations to internal subbank indices using GetSubbankIndex()
 * - Handles data alignment by skipping padding words (0xf0f0f0f0) at buffer start
 * - Iterates through all registered beam detector IDs to find matching devices
 * - Routes data to device-specific ProcessEvBuffer() methods based on device type:
 *   - kQwBPMStripline: 4-channel stripline beam position monitors
 *   - kQwQPD: Quad photodiode detectors
 *   - kQwLinearArray: Linear diode arrays
 *   - kQwBPMCavity: Cavity-type beam position monitors
 *   - kQwBCM: Beam current monitors (various ADC/scaler types)
 *   - kQwClock: Timing reference channels
 *
 * Each device receives a pointer to its specific data section within the buffer,
 * calculated using the device's fWordInSubbank offset. The remaining word count
 * is passed to prevent buffer overruns during device-specific decoding.
 *
 * Buffer Layout Handling:
 * - Accounts for variable-length device data (e.g., 6 words for VQWK, 1 for scalers)
 * - Manages subelement indexing for multi-channel devices
 * - Handles both integrating ADCs and scaler-based measurements
 *
 * Error Conditions:
 * - Invalid ROC/bank combinations are silently ignored (index < 0)
 * - Empty buffers (num_words == 0) are skipped
 * - Buffer alignment issues are detected and corrected automatically
 *
 * \warning This function must be called after LoadChannelMap() to ensure
 * proper device registration and buffer offset calculations.
 */
Int_t QwBeamLine::ProcessEvBuffer(const ROCID_t roc_id, const BankID_t bank_id, UInt_t* buffer, UInt_t num_words)
{
  Bool_t lkDEBUG=kFALSE;
 
  Int_t index = GetSubbankIndex(roc_id,bank_id);
  if (index>=0 && num_words>0){
    //  We want to process this ROC.  Begin looping through the data.
    if (lkDEBUG)
      std::cout << "QwBeamLine::ProcessEvBuffer:  "
        << "Begin processing ROC" << roc_id
        << " and subbank "<<bank_id
        << " number of words="<<num_words<<std::endl;
    if (buffer[0]==0xf0f0f0f0 && num_words%2==1){
      buffer++;
      if (lkDEBUG)
    std::cout << "QwBeamLine::ProcessEvBuffer:  "
          << "Skipped padding word 0xf0f0f0f0 at beginning of buffer."
          << std::endl;
    }
 
    for(size_t i=0;i<fBeamDetectorID.size();i++)
      {
    if(fBeamDetectorID[i].fSubbankIndex==index)
      {
 
        if(fBeamDetectorID[i].fTypeID==kQwBPMStripline)
          {
        if (lkDEBUG)
          {
            std::cout<<"found stripline data for "<<fBeamDetectorID[i].fdetectorname<<std::endl;
            std::cout<<"word left to read in this buffer:"<<num_words-fBeamDetectorID[i].fWordInSubbank<<std::endl;
          }
        fStripline[fBeamDetectorID[i].fIndex].get()->
          ProcessEvBuffer(&(buffer[fBeamDetectorID[i].fWordInSubbank]),
                  num_words-fBeamDetectorID[i].fWordInSubbank,
                  fBeamDetectorID[i].fSubelement);
          }
 
        if(fBeamDetectorID[i].fTypeID==kQwQPD)
          {
        if (lkDEBUG)
          {
            std::cout<<"found qpd data for "<<fBeamDetectorID[i].fdetectorname<<std::endl;
            std::cout<<"word left to read in this buffer:"<<num_words-fBeamDetectorID[i].fWordInSubbank<<std::endl;
          }
        fQPD[fBeamDetectorID[i].fIndex].
          ProcessEvBuffer(&(buffer[fBeamDetectorID[i].fWordInSubbank]),
                  num_words-fBeamDetectorID[i].fWordInSubbank,
                  fBeamDetectorID[i].fSubelement);
          }
 
        if(fBeamDetectorID[i].fTypeID==kQwLinearArray)
          {
        if (lkDEBUG)
          {
            std::cout<<"found linear array data for "<<fBeamDetectorID[i].fdetectorname<<fBeamDetectorID[i].fIndex<<std::endl;
            std::cout<<"word left to read in this buffer:"<<num_words-fBeamDetectorID[i].fWordInSubbank<<std::endl;
          }
        fLinearArray[fBeamDetectorID[i].fIndex].
          ProcessEvBuffer(&(buffer[fBeamDetectorID[i].fWordInSubbank]),
                  num_words-fBeamDetectorID[i].fWordInSubbank,
                  fBeamDetectorID[i].fSubelement);
 
          }
 
        if(fBeamDetectorID[i].fTypeID==kQwBPMCavity)
          {
        if (lkDEBUG)
          {
            std::cout<<"found stripline data for "<<fBeamDetectorID[i].fdetectorname<<std::endl;
            std::cout<<"word left to read in this buffer:"<<num_words-fBeamDetectorID[i].fWordInSubbank<<std::endl;
          }
        fCavity[fBeamDetectorID[i].fIndex].
          ProcessEvBuffer(&(buffer[fBeamDetectorID[i].fWordInSubbank]),
                  num_words-fBeamDetectorID[i].fWordInSubbank,
                  fBeamDetectorID[i].fSubelement);
          }
 
        if(fBeamDetectorID[i].fTypeID==kQwBCM)
          {
        if (lkDEBUG)
          {
            std::cout<<"found bcm data for "<<fBeamDetectorID[i].fdetectorname<<std::endl;
            std::cout<<"word left to read in this buffer:"<<num_words-fBeamDetectorID[i].fWordInSubbank<<std::endl;
          }
        fBCM[fBeamDetectorID[i].fIndex].get()->
          ProcessEvBuffer(&(buffer[fBeamDetectorID[i].fWordInSubbank]),
                  num_words-fBeamDetectorID[i].fWordInSubbank);
          }
 
        if(fBeamDetectorID[i].fTypeID==kQwClock)
          {
        if (lkDEBUG)
          {
            std::cout<<"found clock data for "<<fBeamDetectorID[i].fdetectorname<<std::endl;
            std::cout<<"word left to read in this buffer:"<<num_words-fBeamDetectorID[i].fWordInSubbank<<std::endl;
          }
        fClock[fBeamDetectorID[i].fIndex].get()->
          ProcessEvBuffer(&(buffer[fBeamDetectorID[i].fWordInSubbank]),
                  num_words-fBeamDetectorID[i].fWordInSubbank);
          }
 
        if(fBeamDetectorID[i].fTypeID==kQwHaloMonitor)
          {
        if (lkDEBUG)
          {
            std::cout<<"found halo monitor data for "<<fBeamDetectorID[i].fdetectorname<<std::endl;
            std::cout<<"word left to read in this buffer:"<<num_words-fBeamDetectorID[i].fWordInSubbank<<std::endl;
          }
        fHaloMonitor[fBeamDetectorID[i].fIndex].
          ProcessEvBuffer(&(buffer[fBeamDetectorID[i].fWordInSubbank]),
                  num_words-fBeamDetectorID[i].fWordInSubbank);
          }
 
      }
      }
  }
 
  return 0;
}
 
 
 
//*****************************************************************//
/** Apply single-event cuts to all devices and count failures. */
Bool_t QwBeamLine::ApplySingleEventCuts(){
 
  Bool_t status=kTRUE;
 
  for(size_t i=0;i<fBCM.size();i++){
    status &= fBCM[i].get()->ApplySingleEventCuts();
    if(!status && bDEBUG) std::cout<<"******* QwBeamLine::SingleEventCuts()->BCM[ "<<i
                   <<" , "<<fBCM[i].get()->GetElementName()<<" ] ******\n";
  }
 
  for(size_t i=0;i<fClock.size();i++){
    status &= fClock[i].get()->ApplySingleEventCuts();
    if(!status && bDEBUG) std::cout<<"******* QwBeamLine::SingleEventCuts()->Clock[ "<<i
                   <<" , "<<fClock[i].get()->GetElementName()<<" ] ******\n";
  }
 
  for(size_t i=0;i<fHaloMonitor.size();i++){
    status &= fHaloMonitor[i].ApplySingleEventCuts();
    if(!status && bDEBUG) std::cout<<"******* QwBeamLine::SingleEventCuts()->HaloMonitor[ "<<i
                   <<" , "<<fHaloMonitor[i].GetElementName()<<" ] ******\n";
  }
 
  for(size_t i=0;i<fStripline.size();i++){
    status &= fStripline[i].get()->ApplySingleEventCuts();
    if(!status && bDEBUG) std::cout<<"******** QwBeamLine::SingleEventCuts()->BPMStripline[ "<<i
                   <<" , "<<fStripline[i].get()->GetElementName()<<" ] *****\n";
 
    }
 
  for(size_t i=0;i<fQPD.size();i++){
    status &= fQPD[i].ApplySingleEventCuts();
    if(!status && bDEBUG) std::cout<<"******** QwBeamLine::SingleEventCuts()->QPD[ "<<i
                   <<" , "<<fQPD[i].GetElementName()<<" ] *****\n";
  }
  for(size_t i=0;i<fLinearArray.size();i++){
    status &= fLinearArray[i].ApplySingleEventCuts();
    if(!status && bDEBUG) std::cout<<"******** QwBeamLine::SingleEventCuts()->LinearArray[ "<<i
                   <<" , "<<fLinearArray[i].GetElementName()<<" ] *****\n";
  }
 
  for(size_t i=0;i<fCavity.size();i++){
    status &= fCavity[i].ApplySingleEventCuts();
    if(!status && bDEBUG) std::cout<<"******** QwBeamLine::SingleEventCuts()->BPMCavity[ "<<i
                   <<" , "<<fCavity[i].GetElementName()<<" ] *****\n";
  }
 
  for(size_t i=0;i<fBCMCombo.size();i++){
    status &= fBCMCombo[i].get()->ApplySingleEventCuts();
    if(!status && bDEBUG) std::cout<<"******* QwBeamLine::SingleEventCuts()->CombinedBCM[ "<<i
                   <<" , "<<fBCMCombo[i].get()->GetElementName()<<" ] ******\n";
  }
 
  for(size_t i=0;i<fBPMCombo.size();i++){
    status &= fBPMCombo[i].get()->ApplySingleEventCuts();
    if(!status && bDEBUG) std::cout<<"******* QwBeamLine::SingleEventCuts()->CombinedBPM[ "<<i
                   <<" , "<<fBPMCombo[i].get()->GetElementName()<<" ] ******\n";
 
  }
  for(size_t i=0;i<fECalculator.size();i++){
    status &= fECalculator[i].ApplySingleEventCuts();
    if(!status && bDEBUG) std::cout<<"******* QwBeamLine::SingleEventCuts()->EnergyCalculator[ "<<i
                   <<" , "<<fECalculator[i].GetElementName()<<" ] ******\n";
 
  }
 
 
  //If at least one of the devices failed  event cuts, increment error counter for QwBeamLine
  if (!status)
    fQwBeamLineErrorCount++;
 
 
  return status;
 
}
 
/** Polymorphic per-device burp/burst check vs. a reference subsystem. */
Bool_t QwBeamLine::CheckForBurpFail(const VQwSubsystem *subsys){
  Bool_t burpstatus = kFALSE;
  VQwSubsystem* tmp = const_cast<VQwSubsystem *>(subsys);
  if(Compare(tmp)) {
    const QwBeamLine* input = dynamic_cast<const QwBeamLine*>(subsys);
    for(size_t i=0;i<input->fClock.size();i++){
      //QwError << "************* test Clock *****************" << QwLog::endl;
      burpstatus |= (this->fClock[i].get())->CheckForBurpFail(input->fClock[i].get());
    }
    for(size_t i=0;i<input->fStripline.size();i++){
      //QwError << "************* test stripline *****************" << QwLog::endl;
      burpstatus |= (this->fStripline[i].get())->CheckForBurpFail(input->fStripline[i].get());
    }
    for(size_t i=0;i<input->fQPD.size();i++){
      //QwError << "************* test QPD *****************" << QwLog::endl;
      burpstatus |= (this->fQPD[i]).CheckForBurpFail(&(input->fQPD[i]));
    }
    for(size_t i=0;i<input->fLinearArray.size();i++){
      //QwError << "************* test Lin *****************" << QwLog::endl;
      burpstatus |= (this->fLinearArray[i]).CheckForBurpFail(&(input->fLinearArray[i]));
    }
    for(size_t i=0;i<input->fCavity.size();i++){
      //QwError << "************* test Cavity *****************" << QwLog::endl;
      burpstatus |= (this->fCavity[i]).CheckForBurpFail(&(input->fCavity[i]));
    }
    for(size_t i=0;i<input->fBCM.size();i++){
      //QwError << "************* test BCM *****************" << QwLog::endl;
      burpstatus |= (this->fBCM[i].get())->CheckForBurpFail(input->fBCM[i].get());
    }
    for(size_t i=0;i<input->fBCMCombo.size();i++){
      //QwError << "************* test BCMC *****************" << QwLog::endl;
      burpstatus |= (this->fBCMCombo[i].get())->CheckForBurpFail(input->fBCMCombo[i].get());
    }
    for(size_t i=0;i<input->fBPMCombo.size();i++){
      //QwError << "************* test BPMC *****************" << QwLog::endl;
      burpstatus |= (this->fBPMCombo[i].get())->CheckForBurpFail(input->fBPMCombo[i].get());
    }
    for(size_t i=0;i<input->fECalculator.size();i++){
      //QwError << "************* test ECalc *****************" << QwLog::endl;
      burpstatus |= (this->fECalculator[i]).CheckForBurpFail(&(input->fECalculator[i]));
    }
    for(size_t i=0;i<input->fHaloMonitor.size();i++){
      //QwError << "************* test Halo *****************" << QwLog::endl;
      burpstatus |= (this->fHaloMonitor[i]).CheckForBurpFail(&(input->fHaloMonitor[i]));
    }
  }
  return burpstatus;
}
 
 
 
//*****************************************************************//
/** Print a beamline-wide summary of persistent error counters. */
void QwBeamLine::PrintErrorCounters() const{//inherited from the VQwSubsystemParity; this will display the error summary
 
  QwMessage<<"*********QwBeamLine Error Summary****************"<<QwLog::endl;
  QwVQWK_Channel::PrintErrorCounterHead();
 
  for(size_t i=0;i<fClock.size();i++){
    fClock[i].get()->PrintErrorCounters();
  }
  printf("\n");
  for(size_t i=0;i<fBCM.size();i++){
    fBCM[i].get()->PrintErrorCounters();
  }
  printf("\n");
  for(size_t i=0;i<fHaloMonitor.size();i++){
    fHaloMonitor[i].PrintErrorCounters();
  }
  printf("\n");
  for(size_t i=0;i<fStripline.size();i++){
    fStripline[i].get()->PrintErrorCounters();
    printf("\n");
  }
  for(size_t i=0;i<fQPD.size();i++){
    fQPD[i].PrintErrorCounters();
  }
  printf("\n");
  for(size_t i=0;i<fLinearArray.size();i++){
    fLinearArray[i].PrintErrorCounters();
  }
  printf("\n");
  for(size_t i=0;i<fCavity.size();i++){
    fCavity[i].PrintErrorCounters();
  }
  printf("\n");
  for(size_t i=0;i<fBCMCombo.size();i++){
    fBCMCombo[i].get()->PrintErrorCounters();
  }
  printf("\n");
  for(size_t i=0;i<fBPMCombo.size();i++){
    fBPMCombo[i].get()->PrintErrorCounters();
  }
  printf("\n");
  for(size_t i=0;i<fECalculator.size();i++){
    fECalculator[i].PrintErrorCounters();
  }
  QwVQWK_Channel::PrintErrorCounterTail();
}
 
//*****************************************************************//
/** Increment error counters across all managed devices. */
void QwBeamLine::IncrementErrorCounters()
{
  for(size_t i=0;i<fClock.size();i++){
    fClock[i].get()->IncrementErrorCounters();
  }
  for(size_t i=0;i<fBCM.size();i++){
    fBCM[i].get()->IncrementErrorCounters();
  }
  for(size_t i=0;i<fHaloMonitor.size();i++){
    fHaloMonitor[i].IncrementErrorCounters();
  }
  for(size_t i=0;i<fStripline.size();i++){
    fStripline[i].get()->IncrementErrorCounters();
    }
  for(size_t i=0;i<fQPD.size();i++){
    fQPD[i].IncrementErrorCounters();
  }
  for(size_t i=0;i<fLinearArray.size();i++){
    fLinearArray[i].IncrementErrorCounters();
  }
  for(size_t i=0;i<fCavity.size();i++){
    fCavity[i].IncrementErrorCounters();
  }
  for(size_t i=0;i<fBCMCombo.size();i++){
    fBCMCombo[i].get()->IncrementErrorCounters();
  }
  for(size_t i=0;i<fBPMCombo.size();i++){
    fBPMCombo[i].get()->IncrementErrorCounters();
  }
  for(size_t i=0;i<fECalculator.size();i++){
    fECalculator[i].IncrementErrorCounters();
  }
}
 
//*****************************************************************//
/** Return the OR of per-device event-cut error flags. */
UInt_t QwBeamLine::GetEventcutErrorFlag(){//return the error flag
  UInt_t ErrorFlag;
  UInt_t ErrorFlagtmp;
  ErrorFlag=0;
  for(size_t i=0;i<fBCM.size();i++){
    ErrorFlagtmp = fBCM[i].get()->GetEventcutErrorFlag();
    ErrorFlag |=ErrorFlagtmp;
  }
  for(size_t i=0;i<fStripline.size();i++){
    ErrorFlag |= fStripline[i].get()->GetEventcutErrorFlag();
  }
  for(size_t i=0;i<fQPD.size();i++){
    ErrorFlag |= fQPD[i].GetEventcutErrorFlag();
  }
  for(size_t i=0;i<fLinearArray.size();i++){
    ErrorFlag |= fLinearArray[i].GetEventcutErrorFlag();
  }
  for(size_t i=0;i<fCavity.size();i++){
    ErrorFlag |= fCavity[i].GetEventcutErrorFlag();
  }
  for(size_t i=0;i<fBCMCombo.size();i++){
    ErrorFlag |= fBCMCombo[i].get()->GetEventcutErrorFlag();
  }
  for(size_t i=0;i<fBPMCombo.size();i++){
    ErrorFlag |= fBPMCombo[i].get()->GetEventcutErrorFlag();
  }
  for(size_t i=0;i<fECalculator.size();i++){
    ErrorFlag |= fECalculator[i].GetEventcutErrorFlag();
  }
 
  return ErrorFlag;
 
}
 
//*****************************************************************//
/** Update and return the OR of per-device event-cut error flags. */
UInt_t QwBeamLine::UpdateErrorFlag(){//return the error flag
  UInt_t ErrorFlag;
  UInt_t ErrorFlagtmp;
  ErrorFlag=0;
  for(size_t i=0;i<fBCM.size();i++){
    ErrorFlagtmp = fBCM[i].get()->UpdateErrorFlag();
    ErrorFlag |=ErrorFlagtmp;
  }
  for(size_t i=0;i<fStripline.size();i++){
    ErrorFlag |= fStripline[i].get()->UpdateErrorFlag();
  }
  for(size_t i=0;i<fQPD.size();i++){
    ErrorFlag |= fQPD[i].UpdateErrorFlag();
  }
  for(size_t i=0;i<fLinearArray.size();i++){
    ErrorFlag |= fLinearArray[i].UpdateErrorFlag();
  }
  for(size_t i=0;i<fCavity.size();i++){
    ErrorFlag |= fCavity[i].UpdateErrorFlag();
  }
  for(size_t i=0;i<fBCMCombo.size();i++){
    ErrorFlag |= fBCMCombo[i].get()->UpdateErrorFlag();
  }
  for(size_t i=0;i<fBPMCombo.size();i++){
    ErrorFlag |= fBPMCombo[i].get()->UpdateErrorFlag();
  }
  for(size_t i=0;i<fECalculator.size();i++){
    ErrorFlag |= fECalculator[i].UpdateErrorFlag();
  }
 
  return ErrorFlag;
 
}
 
//*****************************************************************//
/** Copy error flags from a reference beamline into this instance. */
void QwBeamLine::UpdateErrorFlag(const VQwSubsystem *ev_error){
  VQwSubsystem* tmp = const_cast<VQwSubsystem*>(ev_error);
  if(Compare(tmp))
    {
 
      const QwBeamLine* input = dynamic_cast<const QwBeamLine*>(ev_error);
 
       /*
      for(size_t i=0;i<input->fClock.size();i++)
    *(this->fClock[i].get())=*(input->fClock[i].get());
    */
      for(size_t i=0;i<input->fStripline.size();i++)
    (this->fStripline[i].get())->UpdateErrorFlag(input->fStripline[i].get());
      for(size_t i=0;i<input->fQPD.size();i++)
    (this->fQPD[i]).UpdateErrorFlag(&(input->fQPD[i]));
      for(size_t i=0;i<input->fLinearArray.size();i++)
    (this->fLinearArray[i]).UpdateErrorFlag(&(input->fLinearArray[i]));
      for(size_t i=0;i<input->fCavity.size();i++)
    (this->fCavity[i]).UpdateErrorFlag(&(input->fCavity[i]));
      for(size_t i=0;i<input->fBCM.size();i++){
    (this->fBCM[i].get())->UpdateErrorFlag(input->fBCM[i].get());
      }
      for(size_t i=0;i<input->fBCMCombo.size();i++)
    (this->fBCMCombo[i].get())->UpdateErrorFlag(input->fBCMCombo[i].get()); //*(this->fBCMCombo[i].get())=*(input->fBCMCombo[i].get());
      for(size_t i=0;i<input->fBPMCombo.size();i++)
    (this->fBPMCombo[i].get())->UpdateErrorFlag(input->fBPMCombo[i].get()); //=*(input->fBPMCombo[i].get());
      for(size_t i=0;i<input->fECalculator.size();i++)
    (this->fECalculator[i]).UpdateErrorFlag(&(input->fECalculator[i]));
      for(size_t i=0;i<input->fHaloMonitor.size();i++)
    (this->fHaloMonitor[i]).UpdateErrorFlag(&(input->fHaloMonitor[i]));
 
 
 
    }
}
 
 
//*****************************************************************//
/** Process one physics event: clocks first, then all devices. */
void  QwBeamLine::ProcessEvent()
{
  // Make sure this one comes first! The clocks are needed by
  // other elements.
  for(size_t i=0;i<fClock.size();i++)
    fClock[i].get()->ProcessEvent();
 
  for(size_t i=0;i<fStripline.size();i++){
     fStripline[i].get()->ProcessEvent();
    // fStripline[i].get()->PrintInfo();
  }
  for(size_t i=0;i<fCavity.size();i++)
    fCavity[i].ProcessEvent();
 
  for(size_t i=0;i<fBCM.size();i++){
     fBCM[i].get()->ProcessEvent();
    // fBCM[i].get()->PrintInfo();
  }
  for(size_t i=0;i<fQPD.size();i++)
    fQPD[i].ProcessEvent();
 
  for(size_t i=0;i<fLinearArray.size();i++)
    fLinearArray[i].ProcessEvent();
 
  for(size_t i=0;i<fHaloMonitor.size();i++){
    fHaloMonitor[i].ProcessEvent();
  }
 
  for(size_t i=0;i<fBCMCombo.size();i++){
    fBCMCombo[i].get()->ProcessEvent();
     // fBCMCombo[i].get()->PrintInfo();
  }
  for(size_t i=0;i<fBPMCombo.size();i++){
     fBPMCombo[i].get()->ProcessEvent();
    // fBPMCombo[i].get()->PrintInfo();
  }
  for(size_t i=0;i<fECalculator.size();i++){
    fECalculator[i].ProcessEvent();
  }
  return;
}
 
 
//*****************************************************************//
/**
 * Handle configuration banks if present. Placeholder for future use.
 */
Int_t QwBeamLine::ProcessConfigurationBuffer(const ROCID_t roc_id, const BankID_t bank_id, UInt_t* buffer, UInt_t num_words)
{
 
  return 0;
}
 
//*****************************************************************//
/** Publish selected device channels according to the map file. */
Bool_t QwBeamLine::PublishInternalValues() const
{
  // Publish variables
  Bool_t status = kTRUE;
 
  // Publish variables through map file
  // This should work with bcm, bpmstripline, bpmcavity, combo bpm and combo bcm
  for (size_t pp = 0; pp < fPublishList.size(); pp++) {
    TString publish_name = fPublishList.at(pp).at(0);
    TString device_type = fPublishList.at(pp).at(1);
    TString device_name = fPublishList.at(pp).at(2);
    TString device_prop = fPublishList.at(pp).at(3);
    device_type.ToLower();
    device_prop.ToLower();
 
    const VQwHardwareChannel* tmp_channel = 0;
 
    EQwBeamInstrumentType type_id;
    if (device_type == "combobpm")
      type_id = kQwCombinedBPM;
    else if (device_type == "combobcm")
      type_id = kQwCombinedBCM;
    else if (device_type == "comboenergy")
      type_id = kQwEnergyCalculator;
    else if (device_type == "scaler")
      type_id = kQwHaloMonitor;
    else
      type_id = GetQwBeamInstrumentType(device_type);
 
    Int_t index = GetDetectorIndex(type_id,device_name);
 
    if (type_id != kQwUnknownDeviceType
    && index != -1){
      tmp_channel = GetChannel(type_id,index,device_prop);
    } else
      QwError << "QwBeamLine::PublishInternalValues() error "<< QwLog::endl;
 
     if (tmp_channel == NULL) {
       QwError << "QwBeamLine::PublishInternalValues(): " << publish_name << " not found" << QwLog::endl;
       status = kFALSE;
     } else {
       QwDebug << "QwBeamLine::PublishInternalValues(): " << publish_name << " found" << QwLog::endl;
 
       status = PublishInternalValue(publish_name, "published-value", tmp_channel);
     }
  }
 
  return status;
}
 
//*****************************************************************//
/** Publish a device channel on-demand by inferred property suffix. */
Bool_t QwBeamLine::PublishByRequest(TString device_name)
{
  Bool_t status = kFALSE;
  const VQwHardwareChannel* tmp_channel = 0;
 
  std::vector<TString> publishinfo(4,TString(""));
  publishinfo.at(0) = device_name;
 
  EQwBeamInstrumentType type_id;
  Int_t index=-1;
 
  TString name = device_name;
  TString device_prop = "value";
  if (device_name.EndsWith("WS")){
    name = device_name(0,device_name.Length()-2);
    device_prop = "ef";
  } else if (device_name.EndsWith("Q")){
    name = device_name(0,device_name.Length()-1);
    device_prop = "ef";
  } else if (device_name.EndsWith("XSlope")){
    name = device_name(0,device_name.Length()-6);
    device_prop = "xp";
  } else if (device_name.EndsWith("YSlope")){
    name = device_name(0,device_name.Length()-6);
    device_prop = "yp";
  } else if (device_name.EndsWith("X")){
    name = device_name(0,device_name.Length()-1);
    device_prop = "x";
  } else if (device_name.EndsWith("Y")){
    name = device_name(0,device_name.Length()-1);
    device_prop = "y";
  }
 
  for(size_t i=0;i<fBeamDetectorID.size();i++) {
    if(fBeamDetectorID[i].fdetectorname==name
       || fBeamDetectorID[i].fdetectorname==device_name){
      index   = fBeamDetectorID[i].fIndex;
      type_id = fBeamDetectorID[i].fTypeID;
 
      publishinfo.at(1) = GetQwBeamInstrumentTypeName(type_id);
      publishinfo.at(2) = fBeamDetectorID[i].fdetectorname;
      publishinfo.at(3) = device_prop;
      break;
    }
  }
 
  if (index != -1){
    tmp_channel = GetChannel(type_id,index,publishinfo.at(3));
    fPublishList.push_back(publishinfo);
    status = PublishInternalValue(publishinfo.at(0), "published-by-request",
                  tmp_channel);
  }
 
  return status;
}
 
 
//*****************************************************************//
/** Clear current-event state for all managed devices. */
void QwBeamLine::ClearEventData()
{
  for(size_t i=0;i<fClock.size();i++)
    fClock[i].get()->ClearEventData();
  for(size_t i=0;i<fStripline.size();i++)
    fStripline[i].get()->ClearEventData();
  for(size_t i=0;i<fCavity.size();i++)
    fCavity[i].ClearEventData();
  for(size_t i=0;i<fBCM.size();i++)
    fBCM[i]->ClearEventData();
  for(size_t i=0;i<fQPD.size();i++)
    fQPD[i].ClearEventData();
  for(size_t i=0;i<fLinearArray.size();i++)
    fLinearArray[i].ClearEventData();
  for(size_t i=0;i<fHaloMonitor.size();i++)
    fHaloMonitor[i].ClearEventData();
 
  for(size_t i=0;i<fBCMCombo.size();i++)
    fBCMCombo[i].get()->ClearEventData();
  for(size_t i=0;i<fBPMCombo.size();i++)
    fBPMCombo[i].get()->ClearEventData();
  for(size_t i=0;i<fECalculator.size();i++)
    fECalculator[i].ClearEventData();
  return;
}
 
//*****************************************************************//
/** Lookup the element index for a given device type and name. */
Int_t QwBeamLine::GetDetectorIndex( EQwBeamInstrumentType type_id, TString name) const
{
  Bool_t ldebug=kFALSE;
  Int_t result=-1;
  if(ldebug) {
    std::cout<<"QwBeamLine::GetDetectorIndex\n";
    std::cout<<"type_id=="<<type_id<<" name="<<name<<"\n";
    std::cout<<fBeamDetectorID.size()<<" already registered detector\n";
  }
  for(size_t i=0;i<fBeamDetectorID.size();i++) {
    if(ldebug){
      std::cout<<"testing against ("<<fBeamDetectorID[i].fTypeID
           <<","<<fBeamDetectorID[i].fdetectorname<<")=>"<<result<<"\n";
    }
    if(fBeamDetectorID[i].fTypeID==type_id
       && fBeamDetectorID[i].fdetectorname==name){
      result=fBeamDetectorID[i].fIndex;
      break;
    }
  }
  return result;
}
 
/** Overload: retrieve a data element by detector ID. */
VQwDataElement* QwBeamLine::GetElement(QwBeamDetectorID det_id)
{
  return GetElement(det_id.fTypeID, det_id.fIndex);
};
 
/** Overload: retrieve a data element by type and name. */
VQwDataElement* QwBeamLine::GetElement(EQwBeamInstrumentType TypeID, TString name)
{
  Int_t index = GetDetectorIndex(TypeID,name);
  return GetElement(TypeID,index);
};
 
/** Retrieve a data element pointer by type and index with dispatch. */
VQwDataElement* QwBeamLine::GetElement(EQwBeamInstrumentType TypeID, Int_t index)
{
  VQwDataElement* tmp_ptr;
  switch (TypeID){
  case kQwBPMStripline:
    tmp_ptr = fStripline.at(index).get();
    break;
  case kQwQPD:
    tmp_ptr = &(fQPD.at(index));
    break;
  case kQwLinearArray:
    tmp_ptr = &(fLinearArray.at(index));
    break;
  case kQwBCM:
    tmp_ptr = fBCM.at(index).get();
    break;
  case kQwCombinedBCM:
    tmp_ptr = fBCMCombo.at(index).get();
    break;
  case kQwCombinedBPM:
    tmp_ptr = fBPMCombo.at(index).get();
    break;
  case kQwEnergyCalculator:
    tmp_ptr = &(fECalculator.at(index));
    break;
  case kQwHaloMonitor:
    tmp_ptr = &(fHaloMonitor.at(index));
    break;
  case kQwBPMCavity:
    tmp_ptr = &(fCavity.at(index));
    break;
  case kQwClock:
    tmp_ptr = fClock.at(index).get();
    break;
  default:
    TString loc="QwBeamLine::GetElement called by "
      +this->GetName()+" with invalid arguments: "
      +GetQwBeamInstrumentTypeName(TypeID)+" "
      +Form("%d",index);
    throw std::invalid_argument(loc.Data());
  }
  return tmp_ptr;
};
 
/** Const overload: retrieve a data element by type and index. */
const VQwDataElement* QwBeamLine::GetElement(EQwBeamInstrumentType TypeID, Int_t index) const
{
  return const_cast<QwBeamLine*>(this)->GetElement(TypeID,index);
}
 
/**
 * Retrieve a hardware channel (value/x/y/ef/xp/yp) from a specific device.
 */
const VQwHardwareChannel* QwBeamLine::GetChannel(EQwBeamInstrumentType TypeID, Int_t index, TString device_prop) const
{
  const VQwHardwareChannel* tmp_channel = 0;
 
  if (TypeID==kQwBPMStripline || TypeID==kQwBPMCavity || TypeID==kQwQPD){
    const VQwBPM* tmp_ptr = dynamic_cast<const VQwBPM*>(GetElement(TypeID, index));
    if (device_prop == "x")
      tmp_channel = tmp_ptr->GetPosition(VQwBPM::kXAxis);
    else if (device_prop == "y")
      tmp_channel = tmp_ptr->GetPosition(VQwBPM::kYAxis);
    else if (device_prop == "ef")
      tmp_channel = tmp_ptr->GetEffectiveCharge();
  } else if (TypeID==kQwCombinedBPM){
    const VQwBPM* tmp_ptr = dynamic_cast<const VQwBPM*>(GetElement(TypeID, index));
    if (device_prop == "x")
      tmp_channel = tmp_ptr->GetPosition(VQwBPM::kXAxis);
    else if (device_prop == "y")
      tmp_channel = tmp_ptr->GetPosition(VQwBPM::kYAxis);
    else if (device_prop == "ef")
      tmp_channel = fBPMCombo.at(index)->GetEffectiveCharge();
    else if (device_prop == "xp")
      tmp_channel = fBPMCombo.at(index)->GetAngleX();
    else if (device_prop == "yp")
      tmp_channel = fBPMCombo.at(index)->GetAngleY();
  } else if (TypeID==kQwLinearArray){
    /// TODO: QwBeamLine::GetChannel
    /// How do we access linear array channel outputs?
  } else if (TypeID==kQwBCM || TypeID==kQwCombinedBCM){
    tmp_channel = dynamic_cast<const VQwBCM*>(GetElement(TypeID,index))->GetCharge();
  } else if (TypeID==kQwEnergyCalculator){
    tmp_channel = fECalculator.at(index).GetEnergy();
  } else if (TypeID==kQwHaloMonitor){
    tmp_channel = fHaloMonitor.at(index).GetScaler();
  } else if (TypeID==kQwClock){
    tmp_channel = fClock.at(index)->GetTime();
  } else {
    TString loc="QwBeamLine::GetChannel called by "
      +this->GetName()+" with invalid arguments: "
      +GetQwBeamInstrumentTypeName(TypeID)+" "
      +Form("%d",index);
    throw std::invalid_argument(loc.Data());
  }
  return tmp_channel;
}
 
//*****************************************************************//
/** Find a BPMStripline by name. */
VQwBPM* QwBeamLine::GetBPMStripline(const TString name)
{
  if (! fStripline.empty()) {
    for (std::vector<VQwBPM_ptr >::iterator stripline = fStripline.begin(); stripline != fStripline.end(); ++stripline) {
      if ((*stripline).get()->GetElementName() == name) {
    return (*stripline).get();
      }
    }
  }
  return 0;
}
 
//*****************************************************************//
 
/** Find a BPM Cavity by name. */
QwBPMCavity* QwBeamLine::GetBPMCavity(const TString name)
{
  if (! fCavity.empty()) {
    for (std::vector<QwBPMCavity>::iterator cavity = fCavity.begin(); cavity != fCavity.end(); ++cavity) {
      if (cavity->GetElementName() == name) {
    return &(*cavity);
      }
    }
  }
  return 0;
}
 
 
//*****************************************************************//
/** Find a BCM by name. */
VQwBCM* QwBeamLine::GetBCM(const TString name)
{
  //QwWarning << "QwBeamLine::GetBCM" << QwLog::endl;
  if (! fBCM.empty()) {
    for (std::vector<VQwBCM_ptr >::iterator bcm = fBCM.begin(); bcm != fBCM.end(); ++bcm) {
      if ((*bcm).get()->GetElementName() == name) {
    return (*bcm).get();
      }
    }
 
    //QwWarning << "BCM Found" << QwLog::endl;
    return 0;
  }
  return 0;
}
 
 
//*****************************************************************//
/** Find a Clock by name. */
VQwClock* QwBeamLine::GetClock(const TString name)
{
  //QwWarning << "QwBeamLine::GetClock" << QwLog::endl;
  if (! fClock.empty()) {
    for (std::vector<VQwClock_ptr >::iterator clock = fClock.begin(); clock != fClock.end(); ++clock) {
      if ((*clock).get()->GetElementName() == name) {
    return (*clock).get();
      }
    }
 
    //QwWarning << "Clock Found" << QwLog::endl;
    return 0;
  }
  return 0;
}
 
//*****************************************************************//
/** Find a CombinedBCM by name. */
VQwBCM* QwBeamLine::GetCombinedBCM(const TString name)
{
  //QwWarning << "QwBeamLine::GetCombinedBCM" << QwLog::endl;
  if (! fBCMCombo.empty()) {
 
    for (std::vector<VQwBCM_ptr>::iterator cbcm = fBCMCombo.begin(); cbcm != fBCMCombo.end(); ++cbcm) {
      if ((*cbcm).get()->GetElementName() == name) {
    return (*cbcm).get();
      }
    }
 
 
  }
  return 0;
}
 
//*****************************************************************//
/** Find a CombinedBPM by name. */
VQwBPM* QwBeamLine::GetCombinedBPM(const TString name)
{
  //QwWarning << "QwBeamLine::GetCombinedBPM" << QwLog::endl;
  if (! fBPMCombo.empty()) {
 
    for (std::vector<VQwBPM_ptr>::iterator cbpm = fBPMCombo.begin(); cbpm != fBPMCombo.end(); ++cbpm) {
      if ((*cbpm).get()->GetElementName() == name) {
    return (*cbpm).get();
      }
    }
 
 
  }
  return 0;
}
 
//*****************************************************************//
/** Find an EnergyCalculator by name. */
QwEnergyCalculator* QwBeamLine::GetEnergyCalculator(const TString name){
   if (! fECalculator.empty()) {
 
    for (std::vector<QwEnergyCalculator>::iterator ecal = fECalculator.begin(); ecal != fECalculator.end(); ++ecal) {
      if (ecal->GetElementName() == name) {
    return &(*ecal);
      }
    }
 
 
  }
  return 0;
}
 
//*****************************************************************//
/** Const overload: find a BPMStripline by name. */
const VQwBPM* QwBeamLine::GetBPMStripline(const TString name) const
{
  return const_cast<QwBeamLine*>(this)->GetBPMStripline(name);
}
 
QwHaloMonitor* QwBeamLine::GetScalerChannel(const TString name){
   if (! fHaloMonitor.empty()) {
 
    for (std::vector<QwHaloMonitor>::iterator halo = fHaloMonitor.begin(); halo != fHaloMonitor.end(); ++halo) {
      if (halo->GetElementName() == name) {
       return &(*halo);
      }
    }
 
 
  }
  return 0;
};
 
//*****************************************************************//
const QwBPMCavity* QwBeamLine::GetBPMCavity(const TString name) const
{
  return const_cast<QwBeamLine*>(this)->GetBPMCavity(name);
}
 
//*****************************************************************//
const VQwBCM* QwBeamLine::GetBCM(const TString name) const
{
  return const_cast<QwBeamLine*>(this)->GetBCM(name);
}
 
//*****************************************************************//
const VQwClock* QwBeamLine::GetClock(const TString name) const
{
  return const_cast<QwBeamLine*>(this)->GetClock(name);
}
 
//*****************************************************************//
const VQwBCM* QwBeamLine::GetCombinedBCM(const TString name) const{
  return const_cast<QwBeamLine*>(this)->GetCombinedBCM(name);
}
 
//*****************************************************************//
const VQwBPM* QwBeamLine::GetCombinedBPM(const TString name) const{
  return const_cast<QwBeamLine*>(this)->GetCombinedBPM(name);
}
 
//*****************************************************************//
const QwEnergyCalculator* QwBeamLine::GetEnergyCalculator(const TString name) const{
  return const_cast<QwBeamLine*>(this)->GetEnergyCalculator(name);
}
 
 //*****************************************************************//
const QwHaloMonitor* QwBeamLine::GetScalerChannel(const TString name)const {
  return const_cast<QwBeamLine*>(this)->GetScalerChannel(name);
};
 
 
//*****************************************************************//
VQwSubsystem&  QwBeamLine::operator=  (VQwSubsystem *value)
{
  //  std::cout<<" here in QwBeamLine::operator= \n";
  if(Compare(value))
    {
 
      VQwSubsystem::operator=(value);
      QwBeamLine* input = dynamic_cast<QwBeamLine*>(value);
 
      for(size_t i=0;i<input->fClock.size();i++)
    *(this->fClock[i].get())=*(input->fClock[i].get());
      for(size_t i=0;i<input->fStripline.size();i++)
    *(this->fStripline[i].get())=*(input->fStripline[i].get());
      for(size_t i=0;i<input->fQPD.size();i++)
    this->fQPD[i]=input->fQPD[i];
      for(size_t i=0;i<input->fLinearArray.size();i++)
    this->fLinearArray[i]=input->fLinearArray[i];
      for(size_t i=0;i<input->fCavity.size();i++)
    this->fCavity[i]=input->fCavity[i];
      for(size_t i=0;i<input->fBCM.size();i++)
    *(this->fBCM[i].get())=*(input->fBCM[i].get());
      for(size_t i=0;i<input->fHaloMonitor.size();i++)
    this->fHaloMonitor[i]=input->fHaloMonitor[i];
      for(size_t i=0;i<input->fBCMCombo.size();i++)
    *(this->fBCMCombo[i].get())=*(input->fBCMCombo[i].get());
      for(size_t i=0;i<input->fBPMCombo.size();i++)
    *(this->fBPMCombo[i].get())=*(input->fBPMCombo[i].get());
      for(size_t i=0;i<input->fECalculator.size();i++)
    this->fECalculator[i]=input->fECalculator[i];
 
      if (input->fPublishList.size()>0){
    this->fPublishList.resize(input->fPublishList.size());
    for(size_t i=0;i<input->fPublishList.size();i++){
      this->fPublishList.at(i).resize(input->fPublishList.at(i).size());
      for(size_t j=0;j<input->fPublishList.at(i).size();j++){
        this->fPublishList.at(i).at(j)=input->fPublishList.at(i).at(j);
      }
    }
      }
 
    }
  return *this;
}
 
 
//*****************************************************************//
VQwSubsystem&  QwBeamLine::operator+=  (VQwSubsystem *value)
{
  if(Compare(value))
    {
      //QwBeamLine* input= (QwBeamLine*)value ;
      QwBeamLine* input = dynamic_cast<QwBeamLine*>(value);
 
      for(size_t i=0;i<input->fClock.size();i++)
    *(this->fClock[i].get())+=*(input->fClock[i].get());
      for(size_t i=0;i<input->fStripline.size();i++)
    *(this->fStripline[i].get())+=*(input->fStripline[i].get());
      for(size_t i=0;i<input->fCavity.size();i++)
    this->fCavity[i]+=input->fCavity[i];
      for(size_t i=0;i<input->fQPD.size();i++)
    this->fQPD[i]+=input->fQPD[i];
      for(size_t i=0;i<input->fLinearArray.size();i++)
    this->fLinearArray[i]+=input->fLinearArray[i];
      for(size_t i=0;i<input->fBCM.size();i++)
    *(this->fBCM[i].get())+=*(input->fBCM[i].get());
      for(size_t i=0;i<input->fHaloMonitor.size();i++)
    this->fHaloMonitor[i]+=input->fHaloMonitor[i];
      for(size_t i=0;i<input->fBCMCombo.size();i++)
    *(this->fBCMCombo[i].get())+=*(input->fBCMCombo[i].get());
      for(size_t i=0;i<input->fBPMCombo.size();i++)
    *(this->fBPMCombo[i].get())+=*(input->fBPMCombo[i].get());
      for(size_t i=0;i<input->fECalculator.size();i++)
    this->fECalculator[i]+=input->fECalculator[i];
 
      if (input->fPublishList.size()>0){
    this->fPublishList.resize(input->fPublishList.size());
    for(size_t i=0;i<input->fPublishList.size();i++){
      this->fPublishList.at(i).resize(input->fPublishList.at(i).size());
      for(size_t j=0;j<input->fPublishList.at(i).size();j++){
        this->fPublishList.at(i).at(j)=input->fPublishList.at(i).at(j);
      }
    }
      }
 
    }
  return *this;
}
 
//*****************************************************************//
VQwSubsystem&  QwBeamLine::operator-=  (VQwSubsystem *value)
{
 
  if(Compare(value))
    {
      QwBeamLine* input = dynamic_cast<QwBeamLine*>(value);
 
      for(size_t i=0;i<input->fClock.size();i++)
    *(this->fClock[i].get())-=*(input->fClock[i].get());
      for(size_t i=0;i<input->fStripline.size();i++)
    *(this->fStripline[i].get())-=*(input->fStripline[i].get());
      for(size_t i=0;i<input->fCavity.size();i++)
    this->fCavity[i]-=input->fCavity[i];
      for(size_t i=0;i<input->fQPD.size();i++)
    this->fQPD[i]-=input->fQPD[i];
      for(size_t i=0;i<input->fLinearArray.size();i++)
    this->fLinearArray[i]-=input->fLinearArray[i];
      for(size_t i=0;i<input->fBCM.size();i++)
    *(this->fBCM[i].get())-=*(input->fBCM[i].get());
      for(size_t i=0;i<input->fHaloMonitor.size();i++)
    this->fHaloMonitor[i]-=input->fHaloMonitor[i];
      for(size_t i=0;i<input->fBCMCombo.size();i++)
    *(this->fBCMCombo[i].get())-=*(input->fBCMCombo[i].get());
      for(size_t i=0;i<input->fBPMCombo.size();i++)
    *(this->fBPMCombo[i].get())-=*(input->fBPMCombo[i].get());
      for(size_t i=0;i<input->fECalculator.size();i++)
    this->fECalculator[i]-=input->fECalculator[i];
 
      if (input->fPublishList.size()>0){
    this->fPublishList.resize(input->fPublishList.size());
    for(size_t i=0;i<input->fPublishList.size();i++){
      this->fPublishList.at(i).resize(input->fPublishList.at(i).size());
      for(size_t j=0;j<input->fPublishList.at(i).size();j++){
        this->fPublishList.at(i).at(j)=input->fPublishList.at(i).at(j);
      }
    }
      }
 
    }
  return *this;
}
 
//*****************************************************************//
void QwBeamLine::Ratio(VQwSubsystem  *numer, VQwSubsystem  *denom)
{
  if(Compare(numer)&&Compare(denom))
    {
      QwBeamLine* innumer = dynamic_cast<QwBeamLine*>(numer);
      QwBeamLine* indenom = dynamic_cast<QwBeamLine*>(denom);
 
      for(size_t i=0;i<innumer->fClock.size();i++)
    this->fClock[i].get()->Ratio(*(innumer->fClock[i].get()),
      *(indenom->fClock[i].get()));
      for(size_t i=0;i<innumer->fStripline.size();i++)
    this->fStripline[i].get()->Ratio(*(innumer->fStripline[i].get()),
      *(indenom->fStripline[i].get()));
      for(size_t i=0;i<innumer->fCavity.size();i++)
    this->fCavity[i].Ratio(innumer->fCavity[i],indenom->fCavity[i]);
      for(size_t i=0;i<innumer->fQPD.size();i++)
    this->fQPD[i].Ratio(innumer->fQPD[i],indenom->fQPD[i]);
      for(size_t i=0;i<innumer->fLinearArray.size();i++)
    this->fLinearArray[i].Ratio(innumer->fLinearArray[i],indenom->fLinearArray[i]);
      for(size_t i=0;i<innumer->fBCM.size();i++)
    this->fBCM[i].get()->Ratio(*(innumer->fBCM[i].get()),
      *(indenom->fBCM[i].get()));
     for(size_t i=0;i<innumer->fHaloMonitor.size();i++)
    this->fHaloMonitor[i].Ratio(innumer->fHaloMonitor[i],indenom->fHaloMonitor[i]);
      for(size_t i=0;i<innumer->fBCMCombo.size();i++)
    this->fBCMCombo[i].get()->Ratio(*(innumer->fBCMCombo[i].get()),
      *(indenom->fBCMCombo[i]));
      for(size_t i=0;i<innumer->fBPMCombo.size();i++)
    this->fBPMCombo[i].get()->Ratio(*(innumer->fBPMCombo[i].get()),
      *(indenom->fBPMCombo[i].get()));
     for(size_t i=0;i<innumer->fECalculator.size();i++)
       this->fECalculator[i].Ratio(innumer->fECalculator[i],indenom->fECalculator[i]);
 
      // For the combined bcm, maybe we might want to think about getting
      // the asymmetry using the asymmetries of the individual bcms with a
      // weight. But right now it is unclear if really need to have that
      // option.
    }
  return;
}
 
//*****************************************************************//
void QwBeamLine::Scale(Double_t factor)
{
  for(size_t i=0;i<fClock.size();i++)         fClock[i].get()->Scale(factor);
  for(size_t i=0;i<fStripline.size();i++)   fStripline[i].get()->Scale(factor);
  for(size_t i=0;i<fCavity.size();i++)      fCavity[i].Scale(factor);
  for(size_t i=0;i<fQPD.size();i++)         fQPD[i].Scale(factor);
  for(size_t i=0;i<fLinearArray.size();i++) fLinearArray[i].Scale(factor);
  for(size_t i=0;i<fBCM.size();i++)         fBCM[i].get()->Scale(factor);
  for(size_t i=0;i<fHaloMonitor.size();i++) fHaloMonitor[i].Scale(factor);
  for(size_t i=0;i<fBCMCombo.size();i++)    fBCMCombo[i].get()->Scale(factor);
  for(size_t i=0;i<fBPMCombo.size();i++)    fBPMCombo[i].get()->Scale(factor);
  for(size_t i=0;i<fECalculator.size();i++) fECalculator[i].Scale(factor);
  return;
}
 
//*****************************************************************//
void QwBeamLine::CalculateRunningAverage()
{
  for (size_t i = 0; i < fClock.size();          i++) fClock[i].get()->CalculateRunningAverage();
  for (size_t i = 0; i < fStripline.size();    i++) fStripline[i].get()->CalculateRunningAverage();
  for (size_t i = 0; i < fCavity.size();       i++) fCavity[i].CalculateRunningAverage();
  for (size_t i = 0; i < fQPD.size();          i++) fQPD[i].CalculateRunningAverage();
  for (size_t i = 0; i < fLinearArray.size();  i++) fLinearArray[i].CalculateRunningAverage();
  for (size_t i = 0; i < fBCM.size();          i++) fBCM[i].get()->CalculateRunningAverage();
  for (size_t i = 0; i < fHaloMonitor.size();  i++) fHaloMonitor[i].CalculateRunningAverage();
  for (size_t i = 0; i < fBCMCombo.size();     i++) fBCMCombo[i].get()->CalculateRunningAverage();
  for (size_t i = 0; i < fBPMCombo.size();     i++) fBPMCombo[i].get()->CalculateRunningAverage();
  for (size_t i = 0; i < fECalculator.size();  i++) fECalculator[i].CalculateRunningAverage();
}
 
//*****************************************************************//
void QwBeamLine::PrintValue() const
{
  QwMessage << "=== QwBeamLine: " << GetName() << " ===" << QwLog::endl;
  QwMessage << "Clock" << QwLog::endl;
  for (size_t i = 0; i < fClock.size();       i++) fClock[i].get()->PrintValue();
  QwMessage << "BPM stripline" << QwLog::endl;
  for (size_t i = 0; i < fStripline.size(); i++) fStripline[i].get()->PrintValue();
  QwMessage << "QPD" << QwLog::endl;
  for (size_t i = 0; i < fQPD.size(); i++) fQPD[i].PrintValue();
  QwMessage << "LinearArray" << QwLog::endl;
  for (size_t i = 0; i < fLinearArray.size(); i++) fLinearArray[i].PrintValue();
  QwMessage << "BPM cavity" << QwLog::endl;
  for (size_t i = 0; i < fCavity.size(); i++) fCavity[i].PrintValue();
  QwMessage << "BCM" << QwLog::endl;
  for (size_t i = 0; i < fBCM.size();       i++) fBCM[i].get()->PrintValue();
  QwMessage << "HaloMonitor" << QwLog::endl;
  for (size_t i = 0; i < fHaloMonitor.size();       i++) fHaloMonitor[i].PrintValue();
  QwMessage << "BCM combo" << QwLog::endl;
  for (size_t i = 0; i < fBCMCombo.size();  i++) fBCMCombo[i].get()->PrintValue();
  QwMessage << "BPM combo" << QwLog::endl;
  for (size_t i = 0; i < fBPMCombo.size();  i++) fBPMCombo[i].get()->PrintValue();
  QwMessage << "Energy " << QwLog::endl;
  for (size_t i = 0; i < fECalculator.size();  i++) fECalculator[i].PrintValue();
 
}
 
//*****************************************************************//
void QwBeamLine::AccumulateRunningSum(VQwSubsystem* value1, Int_t count, Int_t ErrorMask)
{
  if (Compare(value1)) {
    QwBeamLine* value = dynamic_cast<QwBeamLine*>(value1);
 
    for (size_t i = 0; i < fClock.size();       i++)
      fClock[i].get()->AccumulateRunningSum(*(value->fClock[i].get()), count, ErrorMask);
    for (size_t i = 0; i < fStripline.size(); i++)
      fStripline[i].get()->AccumulateRunningSum(*(value->fStripline[i].get()), count, ErrorMask);
    for (size_t i = 0; i < fCavity.size(); i++)
      fCavity[i].AccumulateRunningSum(value->fCavity[i], count, ErrorMask);
    for (size_t i = 0; i < fBCM.size();       i++)
      fBCM[i].get()->AccumulateRunningSum(*(value->fBCM[i].get()), count, ErrorMask);
    for (size_t i = 0; i < fBCMCombo.size();  i++)
      fBCMCombo[i].get()->AccumulateRunningSum(*(value->fBCMCombo[i].get()), count, ErrorMask);
    for (size_t i = 0; i < fBPMCombo.size();  i++)
      fBPMCombo[i].get()->AccumulateRunningSum(*(value->fBPMCombo[i].get()), count, ErrorMask);
    for (size_t i = 0; i < fECalculator.size();  i++)
      fECalculator[i].AccumulateRunningSum(value->fECalculator[i], count, ErrorMask);
    for (size_t i = 0; i < fQPD.size();  i++)
      fQPD[i].AccumulateRunningSum(value->fQPD[i], count, ErrorMask);
    for (size_t i = 0; i < fLinearArray.size();  i++)
      fLinearArray[i].AccumulateRunningSum(value->fLinearArray[i], count, ErrorMask);
    for (size_t i = 0; i <fHaloMonitor.size();  i++)
      fHaloMonitor[i].AccumulateRunningSum(value->fHaloMonitor[i], count, ErrorMask);
 
  }
}
 
//*****************************************************************//
void QwBeamLine::DeaccumulateRunningSum(VQwSubsystem* value1, Int_t ErrorMask){
    if (Compare(value1)) {
    QwBeamLine* value = dynamic_cast<QwBeamLine*>(value1);
    for (size_t i = 0; i < fClock.size(); i++)
      fClock[i].get()->DeaccumulateRunningSum(*(value->fClock[i].get()), ErrorMask);
    for (size_t i = 0; i < fStripline.size(); i++)
      fStripline[i].get()->DeaccumulateRunningSum(*(value->fStripline[i].get()), ErrorMask);
    for (size_t i = 0; i < fCavity.size(); i++)
      fCavity[i].DeaccumulateRunningSum(value->fCavity[i], ErrorMask);
    for (size_t i = 0; i < fBCM.size();       i++)
      fBCM[i].get()->DeaccumulateRunningSum(*(value->fBCM[i].get()), ErrorMask);
    for (size_t i = 0; i < fBCMCombo.size();  i++)
      fBCMCombo[i].get()->DeaccumulateRunningSum(*(value->fBCMCombo[i].get()), ErrorMask);
    for (size_t i = 0; i < fBPMCombo.size();  i++)
      fBPMCombo[i].get()->DeaccumulateRunningSum(*(value->fBPMCombo[i].get()), ErrorMask);
    for (size_t i = 0; i < fQPD.size();  i++)
      fQPD[i].DeaccumulateRunningSum(value->fQPD[i], ErrorMask);
    for (size_t i = 0; i < fLinearArray.size();  i++)
      fLinearArray[i].DeaccumulateRunningSum(value->fLinearArray[i], ErrorMask);
    for (size_t i = 0; i < fECalculator.size();  i++)
      fECalculator[i].DeaccumulateRunningSum(value->fECalculator[i], ErrorMask);
    for (size_t i = 0; i <fHaloMonitor.size();  i++)
      fHaloMonitor[i].DeaccumulateRunningSum(value->fHaloMonitor[i], ErrorMask);
 
  }
};
 
//*****************************************************************//
Bool_t QwBeamLine::Compare(VQwSubsystem *value)
{
  //  std::cout<<" Here in QwBeamLine::Compare \n";
 
  Bool_t res=kTRUE;
  if(typeid(*value)!=typeid(*this))
    {
      res=kFALSE;
      //      std::cout<<" types are not ok \n";
      //      std::cout<<" this is bypassed just for now but should be fixed eventually \n";
    }
  else
    {
      QwBeamLine* input = dynamic_cast<QwBeamLine*>(value);
      if(input->fStripline.size()!=fStripline.size())
    {
      //      std::cout<<" not the same number of striplines \n";
      res=kFALSE;
    }
      else if(input->fBCM.size()!=fBCM.size())
      {
        res=kFALSE;
      //      std::cout<<" not the same number of bcms \n";
      }
      else if(input->fHaloMonitor.size()!=fHaloMonitor.size())
      {
        res=kFALSE;
      //      std::cout<<" not the same number of halomonitors \n";
      }
      else if(input->fClock.size()!=fClock.size()){
        res=kFALSE;
      //      std::cout<<" not the same number of halomonitors \n";
      }else if(input->fBCMCombo.size()!=fBCMCombo.size()){
    res=kFALSE;
      }else if(input->fBPMCombo.size()!=fBPMCombo.size()){
    res=kFALSE;
      }else if(input->fLinearArray.size()!=fLinearArray.size()){
    res=kFALSE;
      }else if(input->fECalculator.size()!=fECalculator.size()){
    res=kFALSE;
      }else if(input->fCavity.size()!=fCavity.size()){
    res=kFALSE;
      }else if(input->fQPD.size()!=fQPD.size()){
    res=kFALSE;
      }
 
    }
  return res;
}
 
 
//*****************************************************************//
void  QwBeamLine::ConstructHistograms(TDirectory *folder, TString &prefix)
{
 
  //  std::cout<<" here is QwBeamLine::ConstructHistogram with prefix ="<<prefix<<"\n";
  for(size_t i=0;i<fClock.size();i++)
      fClock[i].get()->ConstructHistograms(folder,prefix);
 
  for(size_t i=0;i<fStripline.size();i++)
      fStripline[i].get()->ConstructHistograms(folder,prefix);
 
  for(size_t i=0;i<fQPD.size();i++)
      fQPD[i].ConstructHistograms(folder,prefix);
 
  for(size_t i=0;i<fLinearArray.size();i++)
      fLinearArray[i].ConstructHistograms(folder,prefix);
 
  for(size_t i=0;i<fCavity.size();i++)
      fCavity[i].ConstructHistograms(folder,prefix);
 
  for(size_t i=0;i<fBCM.size();i++)
      fBCM[i].get()->ConstructHistograms(folder,prefix);
 
  for(size_t i=0;i<fHaloMonitor.size();i++)
      fHaloMonitor[i].ConstructHistograms(folder,prefix);
 
  for(size_t i=0;i<fBCMCombo.size();i++)
      fBCMCombo[i].get()->ConstructHistograms(folder,prefix);
 
  for(size_t i=0;i<fBPMCombo.size();i++)
      fBPMCombo[i].get()->ConstructHistograms(folder,prefix);
 
  for(size_t i=0;i<fECalculator.size();i++)
      fECalculator[i].ConstructHistograms(folder,prefix);
  return;
}
 
//*****************************************************************//
void  QwBeamLine::FillHistograms()
{
  for(size_t i=0;i<fClock.size();i++)
    fClock[i].get()->FillHistograms();
  for(size_t i=0;i<fStripline.size();i++)
    fStripline[i].get()->FillHistograms();
  for(size_t i=0;i<fQPD.size();i++)
    fQPD[i].FillHistograms();
  for(size_t i=0;i<fLinearArray.size();i++)
    fLinearArray[i].FillHistograms();
  for(size_t i=0;i<fCavity.size();i++)
    fCavity[i].FillHistograms();
  for(size_t i=0;i<fBCM.size();i++)
    fBCM[i].get()->FillHistograms();
  for(size_t i=0;i<fHaloMonitor.size();i++)
    fHaloMonitor[i].FillHistograms();
  for(size_t i=0;i<fBCMCombo.size();i++)
    fBCMCombo[i].get()->FillHistograms();
  for(size_t i=0;i<fBPMCombo.size();i++)
    fBPMCombo[i].get()->FillHistograms();
  for(size_t i=0;i<fECalculator.size();i++)
    fECalculator[i].FillHistograms();
 
  return;
}
 
 
//*****************************************************************//
void QwBeamLine::ConstructBranchAndVector(TTree *tree, TString & prefix, QwRootTreeBranchVector &values)
{
 
  for(size_t i = 0; i < fClock.size(); i++)
    fClock[i].get()->ConstructBranchAndVector(tree, prefix, values);
  for(size_t i = 0; i < fStripline.size(); i++)
    fStripline[i].get()->ConstructBranchAndVector(tree, prefix, values);
  for(size_t i = 0; i < fQPD.size(); i++)
    fQPD[i].ConstructBranchAndVector(tree, prefix, values);
  for(size_t i = 0; i < fLinearArray.size(); i++)
    fLinearArray[i].ConstructBranchAndVector(tree, prefix, values);
  for(size_t i = 0; i < fCavity.size(); i++)
    fCavity[i].ConstructBranchAndVector(tree, prefix, values);
  for(size_t i = 0; i < fBCM.size(); i++)
    fBCM[i].get()->ConstructBranchAndVector(tree, prefix, values);
  for(size_t i = 0; i < fHaloMonitor.size(); i++)
    fHaloMonitor[i].ConstructBranchAndVector(tree, prefix, values);
  for(size_t i = 0; i <fBCMCombo.size();i++)
    fBCMCombo[i].get()->ConstructBranchAndVector(tree, prefix, values);
  for(size_t i = 0; i <fBPMCombo.size();i++)
    fBPMCombo[i].get()->ConstructBranchAndVector(tree, prefix, values);
  for(size_t i = 0; i <fECalculator.size();i++)
    fECalculator[i].ConstructBranchAndVector(tree, prefix, values);
 
  return;
}
 
//*****************************************************************//
void QwBeamLine::ConstructBranch(TTree *tree, TString & prefix)
{
  for(size_t i = 0; i < fClock.size(); i++)
    fClock[i].get()->ConstructBranch(tree, prefix);
  for(size_t i = 0; i < fStripline.size(); i++)
    fStripline[i].get()->ConstructBranch(tree, prefix);
  for(size_t i = 0; i < fQPD.size(); i++)
    fQPD[i].ConstructBranch(tree, prefix);
  for(size_t i = 0; i < fLinearArray.size(); i++)
    fLinearArray[i].ConstructBranch(tree, prefix);
  for(size_t i = 0; i < fBCM.size(); i++)
    fBCM[i].get()->ConstructBranch(tree, prefix);
  for(size_t i = 0; i <fCavity.size(); i++)
    fStripline[i].get()->ConstructBranch(tree, prefix);
  for(size_t i = 0; i < fHaloMonitor.size(); i++)
    fHaloMonitor[i].ConstructBranch(tree, prefix);
  for(size_t i = 0; i <fBCMCombo.size();i++)
    fBCMCombo[i].get()->ConstructBranch(tree, prefix);
  for(size_t i = 0; i <fBPMCombo.size();i++)
    fBPMCombo[i].get()->ConstructBranch(tree, prefix);
  for(size_t i = 0; i <fECalculator.size();i++)
    fECalculator[i].ConstructBranch(tree, prefix);
 
 
  return;
}
 
//*****************************************************************//
void QwBeamLine::ConstructBranch(TTree *tree, TString & prefix, QwParameterFile& trim_file)
{
  TString tmp,varname,varvalue;
  tmp="QwBCM";
  std::unique_ptr<QwParameterFile> nextmodule;
  trim_file.RewindToFileStart();
 
 
  tmp="QwBPMStripline";
  trim_file.RewindToFileStart();
  if (trim_file.FileHasModuleHeader(tmp)){
    nextmodule=trim_file.ReadUntilNextModule();//This section contains sub modules and or channels to be included in the tree
    for(size_t i = 0; i < fStripline.size(); i++)
      fStripline[i].get()->ConstructBranch(tree, prefix,*nextmodule);
 
  }
 
  tmp="QwQPD";
  trim_file.RewindToFileStart();
  if (trim_file.FileHasModuleHeader(tmp)){
    nextmodule=trim_file.ReadUntilNextModule();//This section contains sub modules and or channels to be included in the tree
    for(size_t i = 0; i < fQPD.size(); i++)
      fQPD[i].ConstructBranch(tree, prefix,*nextmodule);
  }
 
  tmp="QwLinearDiodeArray";
  trim_file.RewindToFileStart();
  if (trim_file.FileHasModuleHeader(tmp)){
    nextmodule=trim_file.ReadUntilNextModule();//This section contains sub modules and or channels to be included in the tree
    for(size_t i = 0; i < fLinearArray.size(); i++)
      fLinearArray[i].ConstructBranch(tree, prefix,*nextmodule);
  }
 
  tmp="QwBPMCavity";
  trim_file.RewindToFileStart();
  if (trim_file.FileHasModuleHeader(tmp)){
    nextmodule=trim_file.ReadUntilNextModule();//This section contains sub modules and or channels to be included in the tree
    for(size_t i = 0; i < fCavity.size(); i++)
      fCavity[i].ConstructBranch(tree, prefix,*nextmodule);
 
  }
 
  tmp="QwBCM";
  trim_file.RewindToFileStart();
  if (trim_file.FileHasModuleHeader(tmp)){
    nextmodule=trim_file.ReadUntilNextModule();//This section contains sub modules and or channels to be included in the tree
    for(size_t i = 0; i < fBCM.size(); i++)
      fBCM[i].get()->ConstructBranch(tree, prefix,*nextmodule);
  }
 
  tmp="QwClock";
  trim_file.RewindToFileStart();
  if (trim_file.FileHasModuleHeader(tmp)){
    nextmodule=trim_file.ReadUntilNextModule();//This section contains sub modules and or channels to be included in the tree
    for(size_t i = 0; i < fClock.size(); i++)
      fClock[i].get()->ConstructBranch(tree, prefix,*nextmodule);
  }
 
  tmp="QwHaloMonitor";
  trim_file.RewindToFileStart();
  if (trim_file.FileHasModuleHeader(tmp)){
    nextmodule=trim_file.ReadUntilNextModule();//This section contains sub modules and or channels to be included in the tree
    for(size_t i = 0; i < fHaloMonitor.size(); i++)
      fHaloMonitor[i].ConstructBranch(tree, prefix,*nextmodule);
  }
 
 
  tmp="QwCombinedBCM";
  trim_file.RewindToFileStart();
  if (trim_file.FileHasModuleHeader(tmp)){
    nextmodule=trim_file.ReadUntilNextModule();//This section contains sub modules and or channels to be included in the tree
    for(size_t i = 0; i <fBCMCombo.size();i++)
      fBCMCombo[i].get()->ConstructBranch(tree, prefix,*nextmodule);
  }
 
 
  tmp="QwCombinedBPM";
  trim_file.RewindToFileStart();
  if (trim_file.FileHasModuleHeader(tmp)){
    nextmodule=trim_file.ReadUntilNextModule();//This section contains sub modules and or channels to be included in the tree
    for(size_t i = 0; i <fBPMCombo.size();i++)
      fBPMCombo[i].get()->ConstructBranch(tree, prefix,*nextmodule);
  }
 
  tmp="QwEnergyCalculator";
  trim_file.RewindToFileStart();
  if (trim_file.FileHasModuleHeader(tmp)){
    nextmodule=trim_file.ReadUntilNextModule();//This section contains sub modules and or channels to be included in the tree
    for(size_t i = 0; i <fECalculator.size();i++)
      fECalculator[i].ConstructBranch(tree, prefix,*nextmodule);
  }
 
  return;
}
 
//*****************************************************************//
void QwBeamLine::FillTreeVector(QwRootTreeBranchVector &values) const
{
  for(size_t i = 0; i < fClock.size(); i++)
    fClock[i].get()->FillTreeVector(values);
  for(size_t i = 0; i < fStripline.size(); i++)
    fStripline[i].get()->FillTreeVector(values);
  for(size_t i = 0; i < fQPD.size(); i++)
    fQPD[i].FillTreeVector(values);
  for(size_t i = 0; i < fLinearArray.size(); i++)
    fLinearArray[i].FillTreeVector(values);
  for(size_t i = 0; i < fCavity.size(); i++)
    fCavity[i].FillTreeVector(values);
  for(size_t i = 0; i < fBCM.size(); i++)
    fBCM[i].get()->FillTreeVector(values);
  for(size_t i = 0; i < fHaloMonitor.size(); i++)
    fHaloMonitor[i].FillTreeVector(values);
  for(size_t i = 0; i < fBCMCombo.size(); i++)
    fBCMCombo[i].get()->FillTreeVector(values);
  for(size_t i = 0; i < fBPMCombo.size(); i++)
    fBPMCombo[i].get()->FillTreeVector(values);
  for(size_t i = 0; i < fECalculator.size(); i++){
    fECalculator[i].FillTreeVector(values);}
  return;
}
 
#ifdef HAS_RNTUPLE_SUPPORT
//*****************************************************************//
void QwBeamLine::ConstructNTupleAndVector(std::unique_ptr<ROOT::RNTupleModel>& model, TString& prefix, std::vector<Double_t>& values, std::vector<std::shared_ptr<Double_t>>& fieldPtrs)
{
  for(size_t i = 0; i < fClock.size(); i++)
    fClock[i].get()->ConstructNTupleAndVector(model, prefix, values, fieldPtrs);
  for(size_t i = 0; i < fStripline.size(); i++)
    fStripline[i].get()->ConstructNTupleAndVector(model, prefix, values, fieldPtrs);
  for(size_t i = 0; i < fQPD.size(); i++)
    fQPD[i].ConstructNTupleAndVector(model, prefix, values, fieldPtrs);
  for(size_t i = 0; i < fLinearArray.size(); i++)
    fLinearArray[i].ConstructNTupleAndVector(model, prefix, values, fieldPtrs);
  for(size_t i = 0; i < fCavity.size(); i++)
    fCavity[i].ConstructNTupleAndVector(model, prefix, values, fieldPtrs);
  for(size_t i = 0; i < fBCM.size(); i++)
    fBCM[i].get()->ConstructNTupleAndVector(model, prefix, values, fieldPtrs);
  for(size_t i = 0; i < fHaloMonitor.size(); i++)
    fHaloMonitor[i].ConstructNTupleAndVector(model, prefix, values, fieldPtrs);
  for(size_t i = 0; i <fBCMCombo.size();i++)
    fBCMCombo[i].get()->ConstructNTupleAndVector(model, prefix, values, fieldPtrs);
  for(size_t i = 0; i <fBPMCombo.size();i++)
    fBPMCombo[i].get()->ConstructNTupleAndVector(model, prefix, values, fieldPtrs);
  for(size_t i = 0; i <fECalculator.size();i++)
    fECalculator[i].ConstructNTupleAndVector(model, prefix, values, fieldPtrs);
 
  return;
}
 
//*****************************************************************//
void QwBeamLine::FillNTupleVector(std::vector<Double_t>& values) const
{
  for(size_t i = 0; i < fClock.size(); i++)
    fClock[i].get()->FillNTupleVector(values);
  for(size_t i = 0; i < fStripline.size(); i++)
    fStripline[i].get()->FillNTupleVector(values);
  for(size_t i = 0; i < fQPD.size(); i++)
    fQPD[i].FillNTupleVector(values);
  for(size_t i = 0; i < fLinearArray.size(); i++)
    fLinearArray[i].FillNTupleVector(values);
  for(size_t i = 0; i < fCavity.size(); i++)
    fCavity[i].FillNTupleVector(values);
  for(size_t i = 0; i < fBCM.size(); i++)
    fBCM[i].get()->FillNTupleVector(values);
  for(size_t i = 0; i < fHaloMonitor.size(); i++)
    fHaloMonitor[i].FillNTupleVector(values);
  for(size_t i = 0; i < fBCMCombo.size(); i++)
    fBCMCombo[i].get()->FillNTupleVector(values);
  for(size_t i = 0; i < fBPMCombo.size(); i++)
    fBPMCombo[i].get()->FillNTupleVector(values);
  for(size_t i = 0; i < fECalculator.size(); i++){
    fECalculator[i].FillNTupleVector(values);}
  return;
}
#endif
 
 
//*****************************************************************//
void  QwBeamLine::PrintInfo() const
{
  std::cout<<"Name of the subsystem ="<<fSystemName<<"\n";
  std::cout<<"there are "<<fClock.size()<<" clock \n";
  std::cout<<"there are "<<fStripline.size()<<" striplines \n";
  std::cout<<"there are "<<fQPD.size()<<" QPDs \n";
  std::cout<<"there are "<<fLinearArray.size()<<" LinearArrays \n";
  std::cout<<"there are "<<fCavity.size()<<" cavities \n";
  std::cout<<"there are "<<fBCM.size()<<" bcm \n";
  std::cout<<"there are "<<fHaloMonitor.size()<<" halomonitors \n";
  std::cout<<"there are "<<fBCMCombo.size()<<" combined bcms \n";
  std::cout<<"there are "<<fBPMCombo.size()<<" combined bpms \n";
  std::cout<<"there are "<<fECalculator.size()<<" energy calculators \n";
  std::cout<<" Printing Running AVG and other channel info for BCMs"<<std::endl;
  for(size_t i=0;i<fBCM.size();i++)
    fBCM[i].get()->PrintInfo();
  for(size_t i=0;i<fHaloMonitor.size();i++)
    fHaloMonitor[i].PrintInfo();
  return;
}
 
 
//*****************************************************************//
void  QwBeamLine::PrintDetectorID() const
{
  for (size_t i=0;i<fBeamDetectorID.size();i++)
    {
      std::cout<<"============================="<<std::endl;
      std::cout<<" Detector ID="<<i<<std::endl;
      fBeamDetectorID[i].Print();
    }
  return;
}
 
 
 
//*****************************************************************//
void  QwBeamLine::CopyTemplatedDataElements(const VQwSubsystem *source)
{
  const QwBeamLine* input = dynamic_cast<const QwBeamLine*>(source);
 
  this->fClock.reserve(input->fClock.size());
  for(size_t i=0;i<input->fClock.size();i++) {
    this->fClock.push_back(VQwClock_ptr(VQwClock::Create(*(input->fClock[i].get()))));
  }
 
  this->fStripline.reserve(input->fStripline.size());
  for(size_t i=0;i<input->fStripline.size();i++) {
    this->fStripline.push_back(VQwBPM_ptr(
        VQwBPM::CreateStripline(*(input->fStripline[i].get()))));
  }
 
  this->fBCM.reserve(input->fBCM.size());
  for(size_t i=0;i<input->fBCM.size();i++) {
    this->fBCM.push_back(VQwBCM_ptr(
        VQwBCM::Create(*(input->fBCM[i].get()))));
  }
 
  this->fBCMCombo.reserve(input->fBCMCombo.size());
  for(size_t i=0;i<input->fBCMCombo.size();i++) {
    this->fBCMCombo.push_back(VQwBCM_ptr(
        VQwBCM::CreateCombo(*(
            input->fBCMCombo[i].get()))));
  }
 
  this->fBPMCombo.reserve(input->fBPMCombo.size());
  for(size_t i=0;i<input->fBPMCombo.size();i++){
    this->fBPMCombo.push_back(VQwBPM_ptr(
        VQwBPM::CreateCombo(*(input->fBPMCombo[i].get()))));
  }
}
 
//*****************************************************************//
#ifdef __USE_DATABASE__
void QwBeamLine::FillDB(QwParityDB *db, TString datatype)
{
 
  Bool_t local_print_flag = false;
 
  if(local_print_flag) {
    QwMessage << " --------------------------------------------------------------- " << QwLog::endl;
    QwMessage << "                         QwBeamLine::FillDB                      " << QwLog::endl;
    QwMessage << " --------------------------------------------------------------- " << QwLog::endl;
  }
 
  std::vector<QwDBInterface> interface;
  std::vector<QwParitySchema::beam_row> entrylist;
 
  UInt_t analysis_id = db->GetAnalysisID();
 
  TString measurement_type_bcm;
  TString measurement_type_bpm;
  TString measurement_type_halo;
 
  measurement_type_bcm =
    QwDBInterface::DetermineMeasurementTypeID(datatype,"q");
  measurement_type_bpm =
    QwDBInterface::DetermineMeasurementTypeID(datatype,"p",kTRUE);
  measurement_type_halo =
    QwDBInterface::DetermineMeasurementTypeID(datatype);
 
  UInt_t i,j;
  i = j = 0;
  // try to access BCM mean and its error
  // there are 2 different types BCM data we have at the moment
  // Yield and Asymmetry
  if(local_print_flag)  QwMessage <<  QwColor(Qw::kGreen) << "Beam Current Monitors" <<QwLog::endl;
 
  for(i=0; i< fBCM.size(); i++) {
    interface.clear();
    interface = fBCM[i].get()->GetDBEntry();
    for (j=0; j<interface.size(); j++){
      interface.at(j).SetAnalysisID( analysis_id );
      interface.at(j).SetMonitorID( db );
      interface.at(j).SetMeasurementTypeID( measurement_type_bcm );
      interface.at(j).PrintStatus( local_print_flag );
      interface.at(j).AddThisEntryToList( entrylist );
    }
  }
 
  ///   try to access BPM mean and its error
  if(local_print_flag) QwMessage <<  QwColor(Qw::kGreen) << "Beam Position Monitors" <<QwLog::endl;
  for(i=0; i< fStripline.size(); i++) {
    //    fStripline[i].MakeBPMList();
    interface.clear();
    interface = fStripline[i].get()->GetDBEntry();
    for (j=0; j<interface.size()-5; j++){
      interface.at(j).SetAnalysisID( analysis_id ) ;
      interface.at(j).SetMonitorID( db );
      interface.at(j).SetMeasurementTypeID( measurement_type_bpm );
      interface.at(j).PrintStatus( local_print_flag);
      interface.at(j).AddThisEntryToList( entrylist );
    }
    // effective charge (last 4 elements)  need to be saved as measurement_type_bcm
    for (j=interface.size()-5; j<interface.size(); j++){
      interface.at(j).SetAnalysisID( analysis_id ) ;
      interface.at(j).SetMonitorID( db );
      interface.at(j).SetMeasurementTypeID( measurement_type_bcm );
      interface.at(j).PrintStatus( local_print_flag);
      interface.at(j).AddThisEntryToList( entrylist );
    }
  }
 
 
  ///   try to access CombinedBPM means and errors
  if(local_print_flag) QwMessage <<  QwColor(Qw::kGreen) << "Combined Beam Position Monitors" <<QwLog::endl;
  for(i=0; i< fBPMCombo.size(); i++) {
    //    fBPMCombo[i].MakeBPMComboList();
    interface.clear();
    interface = fBPMCombo[i].get()->GetDBEntry();
    for (j=0; j<interface.size()-5; j++){
      interface.at(j).SetAnalysisID( analysis_id ) ;
      interface.at(j).SetMonitorID( db );
      interface.at(j).SetMeasurementTypeID( measurement_type_bpm );
      interface.at(j).PrintStatus( local_print_flag);
      interface.at(j).AddThisEntryToList( entrylist );
    }
    // effective charge (last element) need to be saved as measurement_type_bcm
    for (j=interface.size()-5; j<interface.size(); j++){
      interface.at(j).SetAnalysisID( analysis_id ) ;
      interface.at(j).SetMonitorID( db );
      interface.at(j).SetMeasurementTypeID( measurement_type_bcm );
      interface.at(j).PrintStatus( local_print_flag);
      interface.at(j).AddThisEntryToList( entrylist );
    }
  }
 
  ///   try to access CombinedBCM means and errors
  if(local_print_flag) QwMessage <<  QwColor(Qw::kGreen) << "Combined Beam Current Monitors" <<QwLog::endl;
 
  for(i=0; i< fBCMCombo.size(); i++) {
    interface.clear();
    interface = fBCMCombo[i].get()->GetDBEntry();
    for (j=0; j<interface.size(); j++){
      interface.at(j).SetAnalysisID( analysis_id );
      interface.at(j).SetMonitorID( db );
      interface.at(j).SetMeasurementTypeID( measurement_type_bcm );
      interface.at(j).PrintStatus( local_print_flag );
      interface.at(j).AddThisEntryToList( entrylist );
    }
  }
 
  ///   try to access Energy Calculator mean and its error
  if(local_print_flag)  QwMessage <<  QwColor(Qw::kGreen) << "Energy Calculator" <<QwLog::endl;
 
  for(i=0; i< fECalculator.size(); i++) {
    interface.clear();
    interface = fECalculator[i].GetDBEntry();
    for (j=0; j<interface.size(); j++){
      interface.at(j).SetAnalysisID( analysis_id );
      interface.at(j).SetMonitorID( db );
      interface.at(j).SetMeasurementTypeID( measurement_type_bpm );
      interface.at(j).PrintStatus( local_print_flag );
      interface.at(j).AddThisEntryToList( entrylist );
    }
  }
 
 
  ///   try to access QPD mean and its error
  if(local_print_flag) QwMessage <<  QwColor(Qw::kGreen) << "Quadrant PhotoDiodes" <<QwLog::endl;
  for(i=0; i< fQPD.size(); i++) {
    //    fQPD[i].MakeQPDList();
    interface.clear();
    interface = fQPD[i].GetDBEntry();
    for (j=0; j<interface.size()-5; j++){
      interface.at(j).SetAnalysisID( analysis_id ) ;
      interface.at(j).SetMonitorID( db );
      interface.at(j).SetMeasurementTypeID( measurement_type_bpm );
      interface.at(j).PrintStatus( local_print_flag);
      interface.at(j).AddThisEntryToList( entrylist );
    }
    // effective charge need (last element) to be saved as measurement_type_bcm
    for (j=interface.size()-5; j<interface.size(); j++){
      interface.at(j).SetAnalysisID( analysis_id ) ;
      interface.at(j).SetMonitorID( db );
      interface.at(j).SetMeasurementTypeID( measurement_type_bcm );
      interface.at(j).PrintStatus( local_print_flag);
      interface.at(j).AddThisEntryToList( entrylist );
    }
  }
 
  ///   try to access LinearArray mean and its error
  if(local_print_flag) QwMessage <<  QwColor(Qw::kGreen) << "Linear PhotoDiode Array" <<QwLog::endl;
  for(i=0; i< fLinearArray.size(); i++) {
    //    fLinearArray[i].MakeLinearArrayList();
    interface.clear();
    interface = fLinearArray[i].GetDBEntry();
    for (j=0; j<interface.size()-5; j++){
      interface.at(j).SetAnalysisID( analysis_id ) ;
      interface.at(j).SetMonitorID( db );
      interface.at(j).SetMeasurementTypeID( measurement_type_bpm );
      interface.at(j).PrintStatus( local_print_flag);
      interface.at(j).AddThisEntryToList( entrylist );
    }
    for (j=interface.size()-5; j<interface.size(); j++){
      interface.at(j).SetAnalysisID( analysis_id ) ;
      interface.at(j).SetMonitorID( db );
      interface.at(j).SetMeasurementTypeID( measurement_type_bcm );
      interface.at(j).PrintStatus( local_print_flag);
      interface.at(j).AddThisEntryToList( entrylist );
    }
  }
 
  ///   try to access cavity bpm mean and its error
  if(local_print_flag) QwMessage <<  QwColor(Qw::kGreen) << "Cavity Monitors" <<QwLog::endl;
  for(i=0; i< fCavity.size(); i++) {
    //    fCavity[i].MakeBPMCavityList();
    interface.clear();
    interface = fCavity[i].GetDBEntry();
    for (j=0; j<interface.size()-5; j++){
      interface.at(j).SetAnalysisID( analysis_id ) ;
      interface.at(j).SetMonitorID( db );
      interface.at(j).SetMeasurementTypeID( measurement_type_bpm );
      interface.at(j).PrintStatus( local_print_flag);
      interface.at(j).AddThisEntryToList( entrylist );
    }
    for (j=interface.size()-5; j<interface.size(); j++){
      interface.at(j).SetAnalysisID( analysis_id ) ;
      interface.at(j).SetMonitorID( db );
      interface.at(j).SetMeasurementTypeID( measurement_type_bcm );
      interface.at(j).PrintStatus( local_print_flag);
      interface.at(j).AddThisEntryToList( entrylist );
    }
  }
 
  // try to access halo mean and its error
  if(local_print_flag)  QwMessage <<  QwColor(Qw::kGreen) << "Halo Monitors" <<QwLog::endl;
 
  for(i=0; i< fHaloMonitor.size(); i++) {
    interface.clear();
    interface = fHaloMonitor[i].GetDBEntry();
    for (j=0; j<interface.size(); j++){
      interface.at(j).SetAnalysisID( analysis_id );
      interface.at(j).SetMonitorID( db );
      interface.at(j).SetMeasurementTypeID( measurement_type_halo );
      interface.at(j).PrintStatus( local_print_flag );
      interface.at(j).AddThisEntryToList( entrylist );
    }
  }
 
  if(local_print_flag){
    QwMessage << QwColor(Qw::kGreen)   << "Starting database insertion"
              << QwColor(Qw::kNormal)  << QwLog::endl;
  }
 
  // Check the entrylist size, if it isn't zero, start to query..
  if( entrylist.size() ) {
    auto c = db->GetScopedConnection();
    for (const auto& entry: entrylist) {
      c->QueryExecute(entry.insert_into());
    }
  } else {
    QwMessage << "QwBeamLine::FillDB :: This is the case when the entrlylist contains nothing in "<< datatype.Data() << QwLog::endl;
  }
}
 
 
void QwBeamLine::FillErrDB(QwParityDB *db, TString datatype)
{
 
  Bool_t local_print_flag = false;
 
  if(local_print_flag) {
    QwMessage << " --------------------------------------------------------------- " << QwLog::endl;
    QwMessage << "                      QwBeamLine::FillErrDB                      " << QwLog::endl;
    QwMessage << " --------------------------------------------------------------- " << QwLog::endl;
  }
 
  std::vector<QwErrDBInterface> interface;
  std::vector<QwParitySchema::beam_errors_row> entrylist;
 
  UInt_t analysis_id = db->GetAnalysisID();
 
  UInt_t i,j;
  i = j = 0;
 
  if(local_print_flag)  QwMessage <<  QwColor(Qw::kGreen) << "Beam Current Monitors" <<QwLog::endl;
  for(i=0; i< fBCM.size(); i++) {
    interface.clear();
    interface = fBCM[i].get()->GetErrDBEntry();
    for (j=0; j<interface.size(); j++){
      interface.at(j).SetAnalysisID( analysis_id );
      interface.at(j).SetMonitorID( db );
      interface.at(j).PrintStatus( local_print_flag );
      interface.at(j).AddThisEntryToList( entrylist );
    }
    if(local_print_flag) printf("\n");
  }
 
  ///   try to access BPM mean and its error
  if(local_print_flag) QwMessage <<  QwColor(Qw::kGreen) << "Beam Position Monitors" <<QwLog::endl;
  for(i=0; i< fStripline.size(); i++) {
    interface.clear();
    interface = fStripline[i].get()->GetErrDBEntry();
    for (j=0; j<interface.size(); j++){
      interface.at(j).SetAnalysisID( analysis_id ) ;
      interface.at(j).SetMonitorID( db );
      interface.at(j).PrintStatus( local_print_flag );
      interface.at(j).AddThisEntryToList( entrylist );
    }
    if(local_print_flag) printf("\n");
  }
 
  if(local_print_flag) QwMessage <<  QwColor(Qw::kGreen) << "Combined Beam Position Monitors" <<QwLog::endl;
  for(i=0; i< fBPMCombo.size(); i++) {
    interface.clear();
    interface = fBPMCombo[i].get()->GetErrDBEntry();
    for (j=0; j<interface.size(); j++){
      interface.at(j).SetAnalysisID( analysis_id ) ;
      interface.at(j).SetMonitorID( db );
      interface.at(j).PrintStatus( local_print_flag );
      interface.at(j).AddThisEntryToList( entrylist );
    }
    if(local_print_flag) printf("\n");
  }
 
  ///   try to access CombinedBCM means and errors
  if(local_print_flag) QwMessage <<  QwColor(Qw::kGreen) << "Combined Beam Current Monitors" <<QwLog::endl;
  for(i=0; i< fBCMCombo.size(); i++) {
    interface.clear();
    interface = fBCMCombo[i].get()->GetErrDBEntry();
    for (j=0; j<interface.size(); j++){
      interface.at(j).SetAnalysisID( analysis_id ) ;
      interface.at(j).SetMonitorID( db );
      interface.at(j).PrintStatus( local_print_flag );
      interface.at(j).AddThisEntryToList( entrylist );
    }
    if(local_print_flag) printf("\n");
  }
 
  ///   try to access Energy Calculator mean and its error
  if(local_print_flag)  QwMessage <<  QwColor(Qw::kGreen) << "Energy Calculator" <<QwLog::endl;
  for(i=0; i< fECalculator.size(); i++) {
    interface.clear();
    interface = fECalculator[i].GetErrDBEntry();
    for (j=0; j<interface.size(); j++){
      interface.at(j).SetAnalysisID( analysis_id ) ;
      interface.at(j).SetMonitorID( db );
      interface.at(j).PrintStatus( local_print_flag );
      interface.at(j).AddThisEntryToList( entrylist );
    }
    if(local_print_flag) printf("\n");
  }
 
  ///   try to access QPD mean and its error
  if(local_print_flag) QwMessage <<  QwColor(Qw::kGreen) << "Quadrant PhotoDiodes" <<QwLog::endl;
  for(i=0; i< fQPD.size(); i++) {
    interface.clear();
    interface = fQPD[i].GetErrDBEntry();
    for (j=0; j<interface.size(); j++){
      interface.at(j).SetAnalysisID( analysis_id ) ;
      interface.at(j).SetMonitorID( db );
      interface.at(j).PrintStatus( local_print_flag );
      interface.at(j).AddThisEntryToList( entrylist );
    }
    if(local_print_flag) printf("\n");
  }
 
  ///   try to access LinearArray mean and its error
  if(local_print_flag) QwMessage <<  QwColor(Qw::kGreen) << "Linear PhotoDiode Array" <<QwLog::endl;
  for(i=0; i< fLinearArray.size(); i++) {
    interface.clear();
    interface = fLinearArray[i].GetErrDBEntry();
    for (j=0; j<interface.size(); j++){
      interface.at(j).SetAnalysisID( analysis_id ) ;
      interface.at(j).SetMonitorID( db );
      interface.at(j).PrintStatus( local_print_flag );
      interface.at(j).AddThisEntryToList( entrylist );
    }
    if(local_print_flag) printf("\n");
  }
 
  ///   try to access cavity bpm mean and its error
  if(local_print_flag) QwMessage <<  QwColor(Qw::kGreen) << "Cavity Monitors" <<QwLog::endl;
  for(i=0; i< fCavity.size(); i++) {
    interface.clear();
    interface = fCavity[i].GetErrDBEntry();
    for (j=0; j<interface.size(); j++){
      interface.at(j).SetAnalysisID( analysis_id ) ;
      interface.at(j).SetMonitorID( db );
      interface.at(j).PrintStatus( local_print_flag );
      interface.at(j).AddThisEntryToList( entrylist );
    }
    if(local_print_flag) printf("\n");
  }
 
  // try to access halo mean and its error
  if(local_print_flag)  QwMessage <<  QwColor(Qw::kGreen) << "Halo Monitors" <<QwLog::endl;
  for(i=0; i< fHaloMonitor.size(); i++) {
    interface.clear();
    interface = fHaloMonitor[i].GetErrDBEntry();
    for (j=0; j<interface.size(); j++){
      interface.at(j).SetAnalysisID( analysis_id ) ;
      interface.at(j).SetMonitorID( db );
      interface.at(j).PrintStatus( local_print_flag );
      interface.at(j).AddThisEntryToList( entrylist );
    }
    if(local_print_flag) printf("\n");
  }
 
 
  if(local_print_flag){
    QwMessage << QwColor(Qw::kGreen)   << "Starting error database insertion"
              << QwColor(Qw::kNormal)  << QwLog::endl;
  }
 
  // Check the entrylist size, if it isn't zero, start to query..
  if (entrylist.size()) {
    auto c = db->GetScopedConnection();
    for (const auto& entry: entrylist) {
      c->QueryExecute(entry.insert_into());
    }
  } else {
    QwMessage << "QwBeamLine::FillErrDB :: This is the case when the entrlylist contains nothing in "<< datatype.Data() << QwLog::endl;
  }
}
#endif //__USE_DATABASE__
 
void QwBeamLine::WritePromptSummary(QwPromptSummary *ps, TString type)
{
  Bool_t local_print_flag = false;
  Bool_t local_add_element= type.Contains("yield");
 
 
 
  if(local_print_flag) {
    printf("---------------------------------------------------------------\n");
    printf("QwBeamLine::WritePromptSummary()  Type : %12s\n", type.Data());
    printf("---------------------------------------------------------------\n");
  }
 
  const VQwHardwareChannel* tmp_channel = 0;
  TString  element_name        = "";
  Double_t element_value       = 0.0;
  Double_t element_value_err   = 0.0;
  Double_t element_value_width = 0.0;
 
  PromptSummaryElement *local_ps_element = NULL;
  Bool_t local_add_these_elements= false;
 
  // Add BCMs
  for (size_t i = 0; i < fBCM.size();  i++)
    {
      tmp_channel=GetChannel(kQwBCM, i,"");
      element_name        = tmp_channel->GetElementName();
      element_value       = 0.0;
      element_value_err   = 0.0;
      element_value_width = 0.0;
 
      local_add_these_elements =  element_name.EqualTo("bcm_an_us")||element_name.EqualTo("bcm_an_ds")||element_name.EqualTo("bcm_an_ds3")|| (element_name.Contains("cav4") ); // Need to change this to add other BCMs in summary
 
      if(local_add_these_elements && local_add_element){
        ps->AddElement(new PromptSummaryElement(element_name));
      }
 
 
      local_ps_element=ps->GetElementByName(element_name);
 
 
      if(local_ps_element) {
        element_value       = tmp_channel->GetValue();
        element_value_err   = tmp_channel->GetValueError();
        element_value_width = tmp_channel->GetValueWidth();
 
        local_ps_element->Set(type, element_value, element_value_err, element_value_width);
      }
 
      if( local_print_flag && local_ps_element) {
        printf("Type %12s, Element %32s, value %12.4e error %8.4e  width %12.4e\n",
          type.Data(), element_name.Data(), element_value, element_value_err, element_value_width);
      }
    }
 
    // Add combined BCMs
    for (size_t i = 0; i < fBCMCombo.size();  i++) {
      tmp_channel = GetChannel(kQwCombinedBCM, i,"");
      // Need to change this to add other BCMs in summary
      static const TString kTargetBCMName = "bcm_target";
      local_add_these_elements =  element_name.EqualTo(kTargetBCMName);
      element_value_err   = 0.0;
      element_value_width = 0.0;
 
      // Need to change this to add other BCMs in summary
      local_add_these_elements =  element_name.EqualTo("bcm_target");
 
      if(local_add_these_elements && local_add_element){
        ps->AddElement(new PromptSummaryElement(element_name));
      }
 
 
      local_ps_element=ps->GetElementByName(element_name);
 
 
      if(local_ps_element) {
        element_value       = tmp_channel->GetValue();
        element_value_err   = tmp_channel->GetValueError();
        element_value_width = tmp_channel->GetValueWidth();
 
        local_ps_element->Set(type, element_value, element_value_err, element_value_width);
      }
 
      if( local_print_flag && local_ps_element) {
        printf("Type %12s, Element %32s, value %12.4e error %8.4e  width %12.4e\n",
          type.Data(), element_name.Data(), element_value, element_value_err, element_value_width);
      }
    }
 
  // Add BPM Cavity
  for (size_t i = 0; i < fCavity.size();  i++)
    {
      tmp_channel=GetChannel(kQwBPMCavity, i,"ef");
      element_name        = tmp_channel->GetElementName();
      element_value       = 0.0;
      element_value_err   = 0.0;
      element_value_width = 0.0;
 
      local_add_these_elements =  (element_name.EqualTo("cav4bQ")||element_name.EqualTo("cav4cQ")|| element_name.EqualTo("cav4dQ") ); // Need to change this to add other cavities in summary
 
      if(local_add_these_elements && local_add_element){
        ps->AddElement(new PromptSummaryElement(element_name));
      }
 
 
      local_ps_element=ps->GetElementByName(element_name);
 
 
      if(local_ps_element) {
        element_value       = tmp_channel->GetValue();
        element_value_err   = tmp_channel->GetValueError();
        element_value_width = tmp_channel->GetValueWidth();
 
        local_ps_element->Set(type, element_value, element_value_err, element_value_width);
      }
 
      if( local_print_flag && local_ps_element) {
        printf("Type %12s, Element %32s, value %12.4e error %8.4e  width %12.4e\n",
           type.Data(), element_name.Data(), element_value, element_value_err, element_value_width);
      }
    }
 
//////
 
 
 
  char property[2][6]={"x","y"};
 
  local_ps_element=NULL;
  local_add_these_elements=false;
 
  // Add BPM striplines
  for(size_t i=0; i< fStripline.size(); i++)
     {
     for (Int_t j=0;j<2;j++){
      tmp_channel= GetChannel(kQwBPMStripline, i,property[j]);
      element_name= tmp_channel->GetElementName();
      element_value       = 0.0;
      element_value_err   = 0.0;
      element_value_width = 0.0;
 
      // List of BPM names to include in the summary (move to config if needed)
      static const std::vector<TString> kSummaryBPMNames = {
        "bpm4", "bpm18", "bpm14", "bpm12"
      };
      local_add_these_elements = false;
      for (const auto& bpm_name : kSummaryBPMNames) {
        if (element_name.Contains(bpm_name)) {
          local_add_these_elements = true;
          break;
        }
      }
 
      if( local_add_these_elements && local_add_element){
        ps->AddElement(new PromptSummaryElement(element_name));
      }
 
      local_ps_element=ps->GetElementByName(element_name);
 
 
 
      if(local_ps_element) {
        element_value       = tmp_channel->GetValue();
        element_value_err   = tmp_channel->GetValueError();
        element_value_width = tmp_channel->GetValueWidth();
        local_ps_element->Set(type, element_value, element_value_err, element_value_width);
      }
 
      if( local_print_flag && local_ps_element) {
        printf("Type %12s, Element %32s, value %12.4e error %8.4e  width %12.4e\n",
          type.Data(), element_name.Data(), element_value, element_value_err, element_value_width);
      }
 
      }
    }
 
  // Add combined BPM striplines
  for(size_t i=0; i< fBPMCombo.size(); i++)
     {
     for (Int_t j=0;j<2;j++){
      tmp_channel= GetChannel(kQwCombinedBPM, i, property[j]);
      element_name= tmp_channel->GetElementName();
      element_value       = 0.0;
      element_value_err   = 0.0;
      element_value_width = 0.0;
 
      local_add_these_elements=element_name.Contains("bpm4")||element_name.Contains("bpm18")||element_name.Contains("bpm14")||element_name.Contains("bpm12"); //Need to change this to add other stripline BPMs in summary
 
      if( local_add_these_elements && local_add_element){
        ps->AddElement(new PromptSummaryElement(element_name));
      }
 
      local_ps_element=ps->GetElementByName(element_name);
 
 
 
      if(local_ps_element) {
        element_value       = tmp_channel->GetValue();
        element_value_err   = tmp_channel->GetValueError();
        element_value_width = tmp_channel->GetValueWidth();
        local_ps_element->Set(type, element_value, element_value_err, element_value_width);
      }
 
      if( local_print_flag && local_ps_element) {
        printf("Type %12s, Element %32s, value %12.4e error %8.4e  width %12.4e\n",
          type.Data(), element_name.Data(), element_value, element_value_err, element_value_width);
      }
 
      }
    }
 
};