void RecalibrationModel::InitializeFromJSON(Json::Value &recal_params, string &my_block_key, bool spam_enabled) { // this needs to signal when it fails in some way int flowStart, flowEnd, flowSpan, xMin, xMax, xSpan, yMin, yMax, ySpan, max_hp_calibrated; flowStart = rGetParamsInt(recal_params[my_block_key],"flowStart",0); //cout << flowStart << endl; flowEnd = rGetParamsInt(recal_params[my_block_key],"flowEnd",0); //cout << flowEnd << endl; flowSpan = rGetParamsInt(recal_params[my_block_key],"flowSpan",0); //cout << flowSpan << endl; xMin = rGetParamsInt(recal_params[my_block_key],"xMin",0); //cout << xMin << endl; xMax =rGetParamsInt(recal_params[my_block_key],"xMax",0); //cout << xMax << endl; xSpan =rGetParamsInt(recal_params[my_block_key],"xSpan",0); //cout << xSpan << endl; yMin = rGetParamsInt(recal_params[my_block_key],"yMin",0); //cout << yMin << endl; yMax =rGetParamsInt(recal_params[my_block_key],"yMax",0); //cout << yMax << endl; ySpan = rGetParamsInt(recal_params[my_block_key],"ySpan",0); //cout << ySpan << endl; max_hp_calibrated = rGetParamsInt(recal_params[my_block_key],"max_hp_calibrated",0); stratification.SetupRegion(xMin, xMax, xSpan, yMin, yMax, ySpan); //calculate number of partitions and initialize the stratifiedAs and stratifiedBs SetupStratification(flowStart,flowEnd, flowSpan,xMin,xMax,xSpan,yMin,yMax,ySpan,max_hp_calibrated); // stratification setup done // now iterate and obtain each line from the JSON int iter_size = recal_params[my_block_key]["modelParameters"].size(); for (int i_item=0; i_item<iter_size; i_item++) { // extract my single line //model_file >> flowBase >> flowStart >> flowEnd >> xMin >> xMax >> yMin >> yMax >> refHP >> paramA >> paramB; //flowBase is a special extraction flowStart = rGetParamsInt(recal_params[my_block_key]["modelParameters"][i_item],"flowStart",0); flowEnd = rGetParamsInt(recal_params[my_block_key]["modelParameters"][i_item],"flowEnd",0); xMin = rGetParamsInt(recal_params[my_block_key]["modelParameters"][i_item],"xMin",0); xMax =rGetParamsInt(recal_params[my_block_key]["modelParameters"][i_item],"xMax",0); yMin = rGetParamsInt(recal_params[my_block_key]["modelParameters"][i_item],"yMin",0); yMax =rGetParamsInt(recal_params[my_block_key]["modelParameters"][i_item],"yMax",0); int refHP; refHP = rGetParamsInt(recal_params[my_block_key]["modelParameters"][i_item],"refHP",0); float paramA, paramB; paramA = rGetParamsDbl(recal_params[my_block_key]["modelParameters"][i_item],"paramA",1.0); paramB = rGetParamsDbl(recal_params[my_block_key]["modelParameters"][i_item],"paramB",0.0); //string flowBase = recal_params[my_block_key]["modelParameters"][i_item]["flowBase"].asCString(); char flowBase = (char) rGetParamsInt(recal_params[my_block_key]["modelParameters"][i_item],"flowBase",0); int nucInd = NuctoInt(flowBase); // all set with the values int offsetRegion = stratification.OffsetRegion(xMin,yMin); FillIndexes(offsetRegion,nucInd, refHP, flowStart, flowEnd, paramA, paramB); } // now we're done! if (spam_enabled) printf("Recalibration: enabled (using recalibration comment %s)\n\n", my_block_key.c_str()); // if something bad happened above, how do we find out? is_enabled_ = true; }
void RecalibrationModel::InitializeFromJSON(Json::Value &recal_params, string &my_block_key, bool spam_enabled, int over_flow_protect) { // this needs to signal when it fails in some way verbose_ = spam_enabled; int flowStart, flowEnd, flowSpan, xMin, xMax, xSpan, yMin, yMax, ySpan, max_hp_calibrated; flowStart = rGetParamsInt(recal_params[my_block_key],"flowStart",0); //cout << flowStart << endl; flowEnd = rGetParamsInt(recal_params[my_block_key],"flowEnd",0); if (over_flow_protect>flowEnd) flowEnd = over_flow_protect; // allocate for combined bam files that need to access flows "past the end" without crashing //cout << flowEnd << endl; flowSpan = rGetParamsInt(recal_params[my_block_key],"flowSpan",0); //cout << flowSpan << endl; xMin = rGetParamsInt(recal_params[my_block_key],"xMin",0); //cout << xMin << endl; xMax =rGetParamsInt(recal_params[my_block_key],"xMax",0); //cout << xMax << endl; xSpan =rGetParamsInt(recal_params[my_block_key],"xSpan",0); //cout << xSpan << endl; yMin = rGetParamsInt(recal_params[my_block_key],"yMin",0); //cout << yMin << endl; yMax =rGetParamsInt(recal_params[my_block_key],"yMax",0); //cout << yMax << endl; ySpan = rGetParamsInt(recal_params[my_block_key],"ySpan",0); //cout << ySpan << endl; max_hp_calibrated = rGetParamsInt(recal_params[my_block_key],"max_hp_calibrated",0); stratification.SetupChipRegions(xMin, xMax, xSpan, yMin, yMax, ySpan); stratification.SetupFlowRegions(flowStart,flowEnd, flowSpan); //calculate number of partitions and initialize the stratifiedAs and stratifiedBs SetupStratification(flowStart,flowEnd, flowSpan,xMin,xMax,xSpan,yMin,yMax,ySpan,max_hp_calibrated); // stratification setup done // now iterate and obtain each line from the JSON int iter_size = recal_params[my_block_key]["modelParameters"].size(); for (int i_item=0; i_item<iter_size; i_item++) { // extract my single line //model_file >> flowBase >> flowStart >> flowEnd >> xMin >> xMax >> yMin >> yMax >> refHP >> paramA >> paramB; //flowBase is a special extraction flowStart = rGetParamsInt(recal_params[my_block_key]["modelParameters"][i_item],"flowStart",0); flowEnd = rGetParamsInt(recal_params[my_block_key]["modelParameters"][i_item],"flowEnd",0); xMin = rGetParamsInt(recal_params[my_block_key]["modelParameters"][i_item],"xMin",0); xMax =rGetParamsInt(recal_params[my_block_key]["modelParameters"][i_item],"xMax",0); yMin = rGetParamsInt(recal_params[my_block_key]["modelParameters"][i_item],"yMin",0); yMax =rGetParamsInt(recal_params[my_block_key]["modelParameters"][i_item],"yMax",0); int refHP; refHP = rGetParamsInt(recal_params[my_block_key]["modelParameters"][i_item],"refHP",0); float paramA, paramB; paramA = rGetParamsDbl(recal_params[my_block_key]["modelParameters"][i_item],"paramA",1.0); paramB = rGetParamsDbl(recal_params[my_block_key]["modelParameters"][i_item],"paramB",0.0); //string flowBase = recal_params[my_block_key]["modelParameters"][i_item]["flowBase"].asCString(); char flowBase = (char) rGetParamsInt(recal_params[my_block_key]["modelParameters"][i_item],"flowBase",0); int nucInd = NuctoInt(flowBase); // all set with the values int offsetRegion = stratification.OffsetRegion(xMin,yMin); // note we only fill >in< flows fit by the recalibration model FillIndexes(offsetRegion,nucInd, refHP, flowStart, flowEnd, paramA, paramB); } // if something bad happened above, how do we find out? - We do some basic sanity checks CheckArraySize(); // now we're done! if (verbose_) cout << "Recalibration: enabled (using recalibration comment " << my_block_key << ") in a " << stratification.xCuts << 'x' << stratification.yCuts << 'x' << stratification.flowCuts << " grid." <<endl << endl; is_enabled_ = true; }