33 addUsageLine(
"Analyze a PDB file and produce different statistics");
37 addParamsLine(
" distance_histogram <fileOut> <Nnearest=3> <MaxDistance=-1>: Compute the distance histogram between");
39 addParamsLine(
" : The maximum distance of the neighbours may be limited (in Angstroms)");
40 addExampleLine(
"Compute the histogram of interatomic distances",
false);
41 addExampleLine(
"xmipp_pdb_analysis -i mypdb.pdb --operation distance_histogram distance.hist");
48 if (
op==
"distance_histogram")
61 <<
"PDB: " <<
fn_pdb << std::endl
62 <<
"Operation: " <<
op << std::endl;
63 if (
op==
"distance_histogram")
64 std::cout <<
"Output histogram: " <<
fn_hist << std::endl
65 <<
"Nnearest: " <<
Nnearest << std::endl
71 if (
op==
"distance_histogram")
void read(const FileName &fnPDB)
Read phantom from either a PDB of CIF file.
double getDoubleParam(const char *param, int arg=0)
const char * getParam(const char *param, int arg=0)
void distanceHistogramPDB(const PDBPhantom &phantomPDB, size_t Nnearest, double maxDistance, int Nbins, Histogram1D &hist)
void addExampleLine(const char *example, bool verbatim=true)
int verbose
Verbosity level.
void write(const FileName &fn, MDLabel=MDL_X, MDLabel=MDL_COUNT)
void addUsageLine(const char *line, bool verbatim=false)
int getIntParam(const char *param, int arg=0)
void addParamsLine(const String &line)