open2c / polychrom

Chromosome polymer simulations and in silico capture experiments
MIT License
44 stars 16 forks source link

question #72

Open Jincz-lang opened 1 month ago

Jincz-lang commented 1 month ago

radius = sim_object.conlen radiusMult nbCutOffDist = radius repul_energy = ( "rsc12 (rsc2 - 1.0) REPe / emin12 + REPe;" "rsc12 = rsc4 rsc4 rsc4;" "rsc4 = rsc2 rsc2;" "rsc2 = rsc rsc;" "rsc = r / REPsigma rmin12; What is REPsigma, rmin12,rsc?

henrik-dahl-pinholt commented 1 month ago

I don't think I completely understand your question. Do you mean "what do the parameters represent mathematically"? Or do you mean "where are they defined"?

If you mean "what do the parameters represent mathematically" you can see the potential (assuming you are referring to the polynomial_repulsive potential) plotted here with the mathematical formula written out. I believe the parameters you ask about are represented in that equation as:

If you mean "where are they defined", then you can find their definition a little further down in the code using the function addGlobalParameter

I hope that clarifies things a bit. If not, don't hesitate to ask clarifying questions.