本文整理汇总了C++中Peak1D::setPosition方法的典型用法代码示例。如果您正苦于以下问题：C++ Peak1D::setPosition方法的具体用法？C++ Peak1D::setPosition怎么用？C++ Peak1D::setPosition使用的例子？那么, 这里精选的方法代码示例或许可以为您提供帮助。您也可以进一步了解该方法所在类Peak1D的用法示例。
在下文中一共展示了Peak1D::setPosition方法的4个代码示例，这些例子默认根据受欢迎程度排序。您可以为喜欢或者感觉有用的代码点赞，您的评价将有助于系统推荐出更棒的C++代码示例。
示例1: getCIDSpectrumLight_
void CompNovoIdentificationBase::getCIDSpectrumLight_(PeakSpectrum & spec, const String & sequence, DoubleReal prefix, DoubleReal suffix)
{
static DoubleReal h2o_mass = EmpiricalFormula("H2O").getMonoWeight();
Peak1D p;
DoubleReal b_pos(0.0 + prefix);
DoubleReal y_pos(h2o_mass + suffix);
for (Size i = 0; i != sequence.size() - 1; ++i)
{
char aa(sequence[i]);
b_pos += aa_to_weight_[aa];
char aa2(sequence[sequence.size() - i - 1]);
y_pos += aa_to_weight_[aa2];
if (b_pos > min_mz_ && b_pos < max_mz_)
{
p.setPosition(b_pos + Constants::PROTON_MASS_U);
p.setIntensity(1.0f);
spec.push_back(p);
}
if (y_pos > min_mz_ && y_pos < max_mz_)
{
p.setPosition(y_pos + Constants::PROTON_MASS_U);
p.setIntensity(1.0f);
spec.push_back(p);
}
}
spec.sortByPosition();
return;
}示例2: main_
//.........这里部分代码省略.........
for (Size s = 0; s < variable_modifications.size(); ++s)
{
if (mod.getOneLetterCode() == mdb->getModification(variable_modifications[s]).getOrigin() && variable_modifications[s] != mod.getModification())
{
variable_modifications_ok = false;
break;
}
}
}
}
}
if (variable_modifications_ok && fixed_modifications_ok)
{
PeptideIdentification& translocate_pid = *i;
librar.getPeptideIdentifications().push_back(translocate_pid);
librar.setPrecursors(s->getPrecursors());
//library entry transformation
for (UInt l = 0; l < s->size(); ++l)
{
Peak1D peak;
if ((*s)[l].getIntensity() > remove_peaks_below_threshold)
{
const String& info = (*s)[l].getMetaValue("MSPPeakInfo");
if (info[0] == '?')
{
peak.setIntensity(sqrt(0.2 * (*s)[l].getIntensity()));
}
else
{
peak.setIntensity(sqrt((*s)[l].getIntensity()));
}
peak.setMZ((*s)[l].getMZ());
peak.setPosition((*s)[l].getPosition());
librar.push_back(peak);
}
}
if (found != MSLibrary.end())
{
found->second.push_back(librar);
}
else
{
vector<PeakSpectrum> tmp;
tmp.push_back(librar);
MSLibrary.insert(make_pair(MZ_multi, tmp));
}
}
}
}
time_t end_build_time = time(NULL);
cout << "Time needed for preprocessing data: " << (end_build_time - start_build_time) << "\n";
//compare function
PeakSpectrumCompareFunctor* comparor = Factory<PeakSpectrumCompareFunctor>::create(compare_function);
//-------------------------------------------------------------
// calculations
//-------------------------------------------------------------
double score;
StringList::iterator in, out_file;
for (in = in_spec.begin(), out_file = out.begin(); in < in_spec.end(); ++in, ++out_file)
{
time_t start_time = time(NULL);
spectra.load(*in, query);
//Will hold valuable hits
vector<PeptideIdentification> peptide_ids;
vector<ProteinIdentification> protein_ids;示例3: sm
bool ShiftModel1D::finish()
{
if (!isValid())
{
return false;
}
if (!system_)
{
Log.info() << "No valid system found!" << std::endl;
return false;
}
// compute the shift model if necessary
if (compute_shifts_)
{
BALL::ShiftModel sm(parameters_.getFilename());
system_->apply(sm);
}
String element = "";
// Peter Bayer proposed as peak width
// for H 15Hz
// for N 10hz
// for C 5Hz
// peakwidth is meassured in ppm, since
// experiments were done in Hz, we convert the values
// according to the formular
//
// offset [Hz] = offset[ppm] * basic frequency
//
// for our prediction we assume a basic frequency of 700 MHz
float peakwidth = 0.0;
switch(type_)
{
case H:
case H_ON_BACKBONE:
element = "H";
//peakwidth = 0.02142; // Peter Bayers estimation
peakwidth = 0.0032; // this is the former BALL estimation
break;
case N:
case N_BACKBONE:
element = "N";
peakwidth = 0.01428;
break;
case C:
case C_BACKBONE:
element = "C";
peakwidth = 0.00714;
break;
}
int counter = 0;
if (element == "" )
return true;
for (BALL::ResidueIterator r_it = system_->beginResidue(); +r_it; ++r_it)
{
Atom* atom = NULL;
for (BALL::AtomIterator at_it = r_it->beginAtom(); +at_it; ++at_it)
{
if (hasType_(&(*at_it), type_))
{
counter++;
atom = &(*at_it);
// we have, get the shift
float shift = atom->getProperty(BALL::ShiftModule::PROPERTY__SHIFT).getFloat();
Peak1D peak;
float pos = shift;
peak.setPosition(pos);
peak.setWidth(peakwidth);
peak.setIntensity(peak.getIntensity()+1);
//setAtom();
peaks_.push_back(peak);
}
}
}
std::cout << "Number of considered atoms: "<< counter << std::endl;
return true;
}示例4: getCIDSpectrum_
void CompNovoIdentificationBase::getCIDSpectrum_(PeakSpectrum & spec, const String & sequence, Size charge, DoubleReal prefix, DoubleReal suffix)
{
static DoubleReal h2o_mass = EmpiricalFormula("H2O").getMonoWeight();
static DoubleReal nh3_mass = EmpiricalFormula("NH3").getMonoWeight();
static DoubleReal co_mass = EmpiricalFormula("CO").getMonoWeight();
Peak1D p;
DoubleReal b_pos(0 + prefix);
DoubleReal y_pos(h2o_mass + suffix);
bool b_H2O_loss(false), b_NH3_loss(false), y_NH3_loss(false);
for (Size i = 0; i != sequence.size() - 1; ++i)
{
char aa(sequence[i]);
b_pos += aa_to_weight_[aa];
char aa2(sequence[sequence.size() - i - 1]);
y_pos += aa_to_weight_[aa2];
for (Size z = 1; z <= charge && z < 3; ++z)
{
// b-ions
if (b_pos >= min_mz_ && b_pos <= max_mz_)
{
for (Size j = 0; j != max_isotope_; ++j)
{
if (z == 1 /*|| b_pos > MIN_DOUBLE_MZ*/)
{
p.setPosition((b_pos + (DoubleReal)z * Constants::PROTON_MASS_U + (DoubleReal)j + Constants::NEUTRON_MASS_U) / (DoubleReal)z);
p.setIntensity(isotope_distributions_[(Size)b_pos][j] * 0.8 / (z * z));
spec.push_back(p);
}
}
}
// b-ion losses
if (b_pos - h2o_mass > min_mz_ && b_pos - h2o_mass < max_mz_)
{
if (b_H2O_loss || aa == 'S' || aa == 'T' || aa == 'E' || aa == 'D')
{
b_H2O_loss = true;
p.setPosition((b_pos + z * Constants::PROTON_MASS_U - h2o_mass) / z);
p.setIntensity(0.02 / (DoubleReal)(z * z));
if (z == 1 /* || b_pos > MIN_DOUBLE_MZ*/)
{
spec.push_back(p);
}
}
if (b_NH3_loss || aa == 'Q' || aa == 'N' || aa == 'R' || aa == 'K')
{
b_NH3_loss = true;
p.setPosition((b_pos + z * Constants::PROTON_MASS_U - nh3_mass) / z);
p.setIntensity(0.02 / (DoubleReal)(z * z));
if (z == 1 /* || b_pos > MIN_DOUBLE_MZ*/)
{
spec.push_back(p);
}
}
}
// a-ions only for charge 1
if (z == 1)
{
if (b_pos - co_mass > min_mz_ && b_pos - co_mass < max_mz_)
{
// a-ions
p.setPosition((b_pos + z * Constants::PROTON_MASS_U - co_mass) / (DoubleReal)z);
p.setIntensity(0.1f);
spec.push_back(p);
}
}
if (y_pos > min_mz_ && y_pos < max_mz_)
{
// y-ions
for (Size j = 0; j != max_isotope_; ++j)
{
if (z == 1 /* || y_pos > MIN_DOUBLE_MZ*/)
{
p.setPosition((y_pos + (DoubleReal)z * Constants::PROTON_MASS_U + (DoubleReal)j * Constants::NEUTRON_MASS_U) / (DoubleReal)z);
p.setIntensity(isotope_distributions_[(Size)y_pos][j] / (DoubleReal) (z * z));
spec.push_back(p);
}
}
// H2O loss
p.setPosition((y_pos + z * Constants::PROTON_MASS_U - h2o_mass) / (DoubleReal)z);
p.setIntensity(0.1 / (DoubleReal)(z * z));
if (aa2 == 'Q') // pyroglutamic acid formation
{
p.setIntensity(0.5f);
}
if (z == 1 /* || y_pos > MIN_DOUBLE_MZ*/)
{
spec.push_back(p);
}
// NH3 loss
if (y_NH3_loss || aa2 == 'Q' || aa2 == 'N' || aa2 == 'R' || aa2 == 'K')
//.........这里部分代码省略.........