def generate(morphologies=[16,17,18,19]): for m in morphologies: # try: db = Database( db_name='postgresql://hal08.g-node.pri/morphjongleur', exec_role='morphjokey_admin', exec_path='mitsubachi') # must be mapped before Object is created #mapping = mrj.orm.experiment.Mapper( db.engine ) mapping = mrj.orm.morphology.Mapper( db.engine ) mapping.orm_map() morphology = db.load_morphology( m ) print morphology iclamp = mrj.model.clamp.IClamp(morphology.biggest) r = mrj.model.experiment.RecordingPoint(morphology.compartments[len(morphology.compartments)-1]) experiment = mrj.model.experiment.Experiment( morphology, clamps=[iclamp], recording_points=r, neuron_passive_parameter=Neuron_passive_parameter(Ra=35.4,g=0.001) ) db.store( experiment ) print experiment experiment.run_simulation() voltage_trace = experiment.get_voltage_trace() db.store( voltage_trace ) print voltage_trace
def experiment(morphology, compartment, frequency, amplitude): vclamp = mrj.model.clamp.VClamp(compartment=morphology.biggest) iclamp = mrj.util.pattern_generator.SinusClamp(compartment=compartment, amplitude=amplitude, frequency=frequency, duration=3./frequency ) neuron_passive_parameter = mrj.model.neuron_passive.Neuron_passive_parameter(Ra=35.4,g=0.001) experiment = mrj.model.experiment.Experiment( morphology, vclamps=[vclamp], iclamps=[iclamp], neuron_passive_parameter=neuron_passive_parameter, duration=3./frequency, dt=1e-4, description = "SinusExperiment %s %i %f" % (morphology.name, compartment.compartment_id, frequency), ) experiment.run_simulation() voltage_trace = experiment.get_voltage_trace(delay=2./frequency)# nach Einschwinphase return mrj.model.sinusresult.SinusResult(voltage_trace, compartment, frequency, amplitude)