/*******************************************************************************
* Gen.cpp Signal Generator classes
* T.Barnaby, BEAM Ltd, 2006-09-12
*******************************************************************************
*/
#include <Gen.h>
#include <BArray.h>
#include <TmsLib.h>
#include <math.h>
#define AVG_PRINT 0
const double pllFrequency = 125e6; // The PLL frequency
const int defaultFrefGain = 4096; // The default PLL FREF gain
const int defaultPllGain = 9; // The default PLL feedback gain (In shifts to the right)
using namespace Tms;
BString PNameValueList::getValue(BString name){
BString* s;
if(s = find(name)){
return *s;
}
return "";
}
void PNameValueList::setValue(BString name, BString value){
BString* s;
if(s = find(name)){
*s = value;
return;
}
else {
append(PNameValue(name, value));
// BList<GenParam>::append(GenParam(name, value));
}
}
Gen::Gen(){
osampleRate = 125000000.0;
ofref = 437000.0;
ophaseTableFile = 0;
ouseMsbFref = 0;
opllInitialFrequencyDelay = 0;
ophaseDelay = 0;
ophaseDelayAll = -1;
}
Gen::~Gen(){
}
BError Gen::init(PNameValueList& params){
BError err;
BString s;
oparams = params;
if((s = params.getValue("fileName")) != "")
ofileName = s;
if((s = params.getValue("sampleRate")) != "")
osampleRate = atof(s);
if((s = params.getValue("fref")) != "")
ofref = atof(s);
if((s = params.getValue("phaseTableFile")) != "")
ophaseTableFile = 1;
if((s = params.getValue("useMsbFref")) != "")
ouseMsbFref = atoi(s);
if((s = params.getValue("pllInitialFrequencyDelay")) != "")
opllInitialFrequencyDelay = atoi(s);
if((s = params.getValue("phaseDelay")) != "")
ophaseDelay = atoi(s);
if((s = params.getValue("phaseDelayAll")) != "")
ophaseDelayAll = atoi(s);
if((s = params.getValue("name")) != "")
oname = s;
if((s = params.getValue("info")) != "")
ocycleParam.info = s + " ";
return ofile.open(ofileName, "w");
}
BList<BString> Gen::getTypes(){
BList<BString> sl;
sl.append("Beam1: Beam with 4 particle bunches at harmonic 8.");
sl.append("Beam2: Beam with 4 particle bunches at harmonic 8. Beam then moves to 4 particle bunches at harmonic 16.");
sl.append("Beam3: Beam with 4 particle bunches at harmonic 8. Beam then moves to 8 particle bunches at harmonic 16.");
sl.append("Beam3Delat: Beam with 4 particle bunches at harmonic 8. Beam then moves to 8 particle bunches at harmonic 16. Has Delay");
sl.append("Test1: Beam with 4 particle bunches at harmonic 8. Beam then moves to 8 particle bunches at harmonic 16.");
return sl;
}
BError Gen::generate(BString cycleType){
BError err;
if(cycleType == "Beam1")
err = genBeam1();
else if(cycleType == "Beam2")
err = genBeam2();
else if(cycleType == "Beam3")
err = genBeam3();
else if(cycleType == "Beam3Delay")
err = genBeam3Delay();
else if(cycleType == "Test1")
err = genTest1();
else
err.set(1, "Unknown cycle type");
if(!err && ophaseTableFile){
BFile f;
unsigned int s;
unsigned int n;
TmsPhase p;
for(s = 0; s < ocycleParam.stateTable.size(); s++){
f.open(BString("phaseTable") + s + ".txt", "w");
for(n = 0; n < ocycleParam.stateTable[s].phaseTable.size(); n++){
p.value = ocycleParam.stateTable[s].phaseTable[n];
f.printf("%d ", p.lo1);
f.printf("%d ", p.lo2);
f.printf("%d ", p.gate);
f.printf("%d ", p.blr);
f.printf("%d ", p.meanFilter1);
f.printf("%d\n", p.meanFilter2);
}
f.close();
}
}
return err;
}
void Gen::setPhaseDelay(){
BUInt32 n;
ocycleParam.getdefaultPickupPositions(ocycleParam.frefPhaseDelay);
if(ophaseDelayAll >= 0){
for(n = 0; n < ocycleParam.frefPhaseDelay.size(); n++){
ocycleParam.frefPhaseDelay[n] = ophaseDelayAll;
}
}
else {
for(n = 0; n < ocycleParam.frefPhaseDelay.size(); n++){
ocycleParam.frefPhaseDelay[n] += ophaseDelay;
}
}
}
BError Gen::genBeam1(){
BError err;
CycleParamState st;
BList<CycleParamState> cpl;
BUInt num = 0;
// Set up base parameters
ocycleParam.cycleType = oname;
ocycleParam.info += "Beam with 4 particle bunches at harmonic 8.";
ocycleParam.channel = 0;
ocycleParam.pllCycleStartFrequency = int(pow(2.0, 32) / (pllFrequency / ofref));
ocycleParam.pllInitialFrequency = int(pow(2.0, 32) / (pllFrequency / ofref));
ocycleParam.pllInitialFrequencyDelay = opllInitialFrequencyDelay;
ocycleParam.pllFrefGain = defaultFrefGain;
ocycleParam.pllGain = defaultPllGain;
ocycleParam.pllDdsMinimum = 0;
ocycleParam.pllDdsMaximum = 0;
// Set up Positions
setPhaseDelay();
// Set up States
ocycleParam.getDefaultState(st);
st.setNext(num++, CyclePeriodAll, 0, 1, 1, 0, 0, 1);
st.state.calStart = 1;
st.state.injection = 3;
st.lo1Harmonic = 8;
st.lo1Phase = 0.0;
st.lo2Harmonic = 8;
st.lo2Phase = 0.0;
st.bunchMask = 0x3C;
st.mean1Mask = 0x3C;
st.mean2Mask = 0x04;
cpl.append(st);
st.setNext(num++, CyclePeriodCalibration, 1, 0, 1, 0, 0, 1);
st.state.calStop = 2;
cpl.append(st);
st.setNext(num++, CyclePeriodAll, 0, 1, 1, 0, 0, 1);
st.state.injection = 3;
cpl.append(st);
st.setNext(num++, CyclePeriodEvent0, 1, 0, 1, 0, 0, 1);
cpl.append(st);
if(err = ocycleParam.setStates(cpl))
return err;
ofile.writeString(ocycleParam.getString());
ofile.close();
return err;
}
BError Gen::genBeam2(){
BError err;
CycleParamState st;
BList<CycleParamState> cpl;
BUInt num = 0;
// Set up base parameters
ocycleParam.cycleType = oname;
ocycleParam.info += "Beam with 4 particle bunches at harmonic 8. Beam then moves to 4 particle bunches at harmonic 16.";
ocycleParam.channel = 0;
ocycleParam.pllCycleStartFrequency = int(pow(2.0, 32) / (pllFrequency / ofref));
ocycleParam.pllInitialFrequency = int(pow(2.0, 32) / (pllFrequency / ofref));
ocycleParam.pllInitialFrequencyDelay = opllInitialFrequencyDelay;
ocycleParam.pllFrefGain = defaultFrefGain;
ocycleParam.pllGain = defaultPllGain;
ocycleParam.pllDdsMinimum = 0;
ocycleParam.pllDdsMaximum = 0;
// Set up Positions
setPhaseDelay();
// Set up States
ocycleParam.getDefaultState(st);
st.setNext(num++, CyclePeriodAll, 0, 1, 1, 0, 0, 1);
st.state.calStart = 1;
st.state.injection = 3;
st.lo1Harmonic = 8;
st.lo1Phase = 0.0;
st.lo2Harmonic = 8;
st.lo2Phase = 0.0;
st.bunchMask = 0x3C;
st.mean1Mask = 0x3C;
st.mean2Mask = 0x04;
cpl.append(st);
st.setNext(num++, CyclePeriodCalibration, 1, 0, 1, 0, 0, 1);
st.state.calStop = 2;
cpl.append(st);
st.setNext(num++, CyclePeriodAll, 0, 1, 1, 0, 0, 1);
st.state.injection = 3;
cpl.append(st);
st.setNext(num++, CyclePeriodEvent0, 1, 0, 1, 0, 0, 1);
st.state.hchange = 4;
cpl.append(st);
st.setNext(num++, CyclePeriodEvent1, 1, 0, 1, 0, 0, 1);
st.lo1Harmonic = 16;
st.lo1Phase = 0.0;
st.bunchMask = 0x03C;
st.mean1Mask = 0x03C;
st.mean2Mask = 0x004;
cpl.append(st);
if(err = ocycleParam.setStates(cpl))
return err;
ofile.writeString(ocycleParam.getString());
ofile.close();
return err;
}
BError Gen::genBeam3(){
BError err;
CycleParamState st;
BList<CycleParamState> cpl;
BUInt num = 0;
// Set up base parameters
ocycleParam.cycleType = oname;
ocycleParam.info += "Beam with 4 particle bunches at harmonic 8. Beam then moves to 4 particle bunches at harmonic 16.";
ocycleParam.channel = 0;
ocycleParam.pllCycleStartFrequency = int(pow(2.0, 32) / (pllFrequency / ofref));
ocycleParam.pllInitialFrequency = int(pow(2.0, 32) / (pllFrequency / ofref));
ocycleParam.pllInitialFrequencyDelay = opllInitialFrequencyDelay;
ocycleParam.pllFrefGain = defaultFrefGain;
ocycleParam.pllGain = defaultPllGain;
ocycleParam.pllDdsMinimum = 0;
ocycleParam.pllDdsMaximum = 0;
// Set up Positions
setPhaseDelay();
// Set up States
ocycleParam.getDefaultState(st);
st.setNext(num++, CyclePeriodAll, 0, 1, 1, 0, 0, 1);
st.state.calStart = 1;
st.state.injection = 3;
st.lo1Harmonic = 8;
st.lo1Phase = 0.0;
st.lo2Harmonic = 8;
st.lo2Phase = 0.0;
st.bunchMask = 0x3C;
st.mean1Mask = 0x3C;
st.mean2Mask = 0x04;
cpl.append(st);
st.setNext(num++, CyclePeriodCalibration, 1, 0, 1, 0, 0, 1);
st.state.calStop = 2;
cpl.append(st);
st.setNext(num++, CyclePeriodAll, 0, 1, 1, 0, 0, 1);
st.state.injection = 3;
cpl.append(st);
st.setNext(num++, CyclePeriodEvent0, 1, 0, 1, 0, 0, 1);
st.state.hchange = 4;
cpl.append(st);
st.setNext(num++, CyclePeriodEvent1, 1, 0, 1, 0, 0, 1);
st.lo1Harmonic = 16;
st.lo1Phase = 0.0;
st.bunchMask = 0x03C0;
st.mean1Mask = 0x03C0;
st.mean2Mask = 0x0040;
cpl.append(st);
if(err = ocycleParam.setStates(cpl))
return err;
ofile.writeString(ocycleParam.getString());
ofile.close();
return err;
}
BError Gen::genBeam3Delay(){
BError err;
CycleParamState st;
BList<CycleParamState> cpl;
BUInt num = 0;
// Set up base parameters
ocycleParam.cycleType = oname;
ocycleParam.info += "Beam with 4 particle bunches at harmonic 8. Beam then moves to 4 particle bunches at harmonic 16. With delay state.";
ocycleParam.channel = 0;
ocycleParam.pllCycleStartFrequency = int(pow(2.0, 32) / (pllFrequency / ofref));
ocycleParam.pllInitialFrequency = int(pow(2.0, 32) / (pllFrequency / ofref));
ocycleParam.pllInitialFrequencyDelay = opllInitialFrequencyDelay;
ocycleParam.pllFrefGain = defaultFrefGain;
ocycleParam.pllGain = defaultPllGain;
ocycleParam.pllDdsMinimum = 0;
ocycleParam.pllDdsMaximum = 0;
ocycleParam.stateDelay = 5;
// Set up Positions
setPhaseDelay();
// Set up States
ocycleParam.getDefaultState(st);
st.setNext(num++, CyclePeriodAll, 0, 1, 1, 0, 0, 1);
st.state.calStart = 1;
st.state.injection = 3;
st.lo1Harmonic = 8;
st.lo1Phase = 0.0;
st.lo2Harmonic = 8;
st.lo2Phase = 0.0;
st.bunchMask = 0x3C;
st.mean1Mask = 0x3C;
st.mean2Mask = 0x04;
cpl.append(st);
st.setNext(num++, CyclePeriodCalibration, 1, 0, 1, 0, 0, 1);
st.state.calStop = 2;
cpl.append(st);
st.setNext(num++, CyclePeriodAll, 0, 1, 1, 0, 0, 1);
st.state.injection = 3;
cpl.append(st);
st.setNext(num++, CyclePeriodEvent0, 1, 0, 1, 0, 0, 1);
st.state.delay = 4;
st.lo1Harmonic = 8;
st.lo1Phase = 0.0;
cpl.append(st);
st.setNext(num++, CyclePeriodEvent0, 1, 0, 1, 0, 0, 1);
st.state.hchange = 5;
st.lo1Harmonic = 8;
st.lo1Phase = 0.0;
cpl.append(st);
st.setNext(num++, CyclePeriodEvent1, 1, 0, 1, 0, 0, 1);
st.lo1Harmonic = 16;
st.lo1Phase = 0.0;
st.bunchMask = 0x03C0;
st.mean1Mask = 0x03C0;
st.mean2Mask = 0x0040;
cpl.append(st);
if(err = ocycleParam.setStates(cpl))
return err;
ofile.writeString(ocycleParam.getString());
ofile.close();
return err;
}
BError Gen::genTest1(){
BError err;
CycleParamState st;
BList<CycleParamState> cpl;
BUInt num = 0;
// Set up base parameters
ocycleParam.cycleType = oname;
ocycleParam.info += "Beam with 4 particle bunches at harmonic 8. Beam then moves to 4 particle bunches at harmonic 16.";
ocycleParam.channel = 0;
ocycleParam.pllCycleStartFrequency = int(pow(2.0, 32) / (pllFrequency / ofref));
ocycleParam.pllInitialFrequency = int(pow(2.0, 32) / (pllFrequency / ofref));
ocycleParam.pllInitialFrequencyDelay = opllInitialFrequencyDelay;
ocycleParam.pllFrefGain = defaultFrefGain;
ocycleParam.pllGain = defaultPllGain;
ocycleParam.pllDdsMinimum = 0;
ocycleParam.pllDdsMaximum = 0;
// Set up Positions
setPhaseDelay();
// Set up States
ocycleParam.getDefaultState(st);
st.setNext(num++, CyclePeriodAll, 0, 1, 1, 0, 0, 1);
st.state.calStart = 1;
st.state.injection = 3;
st.lo1Harmonic = 8;
st.lo1Phase = 0.0;
st.lo2Harmonic = 8;
st.lo2Phase = 0.0;
st.bunchMask = 0x3C;
st.mean1Mask = 0x3C;
st.mean2Mask = 0x04;
cpl.append(st);
st.setNext(num++, CyclePeriodCalibration, 1, 0, 1, 0, 0, 1);
st.state.calStop = 2;
cpl.append(st);
st.setNext(num++, CyclePeriodAll, 0, 1, 1, 0, 0, 1);
st.state.injection = 3;
cpl.append(st);
st.setNext(num++, CyclePeriodEvent0, 1, 0, 1, 0, 0, 1);
st.state.hchange = 4;
cpl.append(st);
st.setNext(num++, CyclePeriodEvent1, 0, 1, 1, 0, 0, 1);
st.lo1Harmonic = 16;
st.lo1Phase = 0.0;
st.bunchMask = 0x03C0;
st.mean1Mask = 0x03C0;
st.mean2Mask = 0x0040;
cpl.append(st);
if(err = ocycleParam.setStates(cpl))
return err;
ofile.writeString(ocycleParam.getString());
ofile.close();
return err;
}
