A model of cholinergic modulation in olfactory bulb (OB) and piriform
cortex (PC).

Licurgo de Almeida, Marco Idiart and Christiane Linster, 2013.

So far the model was used in the following publications:

* de Almeida, L., Idiart, M., & Linster, C. (2013). A model of
  cholinergic modulation in olfactory bulb and piriform
  cortex. Journal of neurophysiology, 109(5), 1360-1377.
* Devore, S., de Almeida, L., & Linster, C. (2014). Distinct roles of
  bulbar muscarinic and nicotinic receptors in olfactory
  discrimination learning. Journal of Neuroscience, 34(34):
  11244-11260.

All codes presented here were implemented using Matlab 2011a and might
not work on other similar programs (i.e. Octave or very old versions
of Matlab).

To create the Bulb network using the integrate and fire model use the
following command:

[OSN,Pg,M1,M2,Gr] = CreateBulb(tsim,ncell,mod);

where:

tsim: simulation time
ncells: number of cells in each cell group
mod is either true (for mod ON) or false (for mod OFF)
The output are a group of networks:
OSN: Olfactory Sensory Neurons;
Pg: PG cells;
M1: Apical compartment from Mitral cells;
M2: Soma (spiking) compartment from Mitral cells;
Gr: Granule cells;

Parameters can be changed inside SetBulbParam.m and information about
these parameters are found in the publications cited above.

To create Mitral cell output using the simplified version of the model
(see de Almeida et al., 2013 for reasoning about the differences) use
the following command:

M2 = CreateSimpleBulb(tsim,ncell);
Where:
tsim: simulation time
ncells: number of cells in each cell group
The output is a Mitral network (M2) only.
Parameters can be changed inside SetBulbSimpleParam.m and information
about these parameters are found in the publications cited above.

You can plot the activation of a group of cells with the command:
Raster(M2);

After creating a pool of Mitral cells, one can use these neuronal
activities to feed the PC network with the following command:
[Ff,Py,Fb] = CreateCortex(M2,mod);
where:
M2 is the Mitral cell soma compartment generated by CreateBulb.m or
CreateSimpleBulb.m and mod is either true (for modulation ON) or false
(for modulation OFF)

The output are a group of networks:
OSN: Olfactory Sensory Neurons;
Ff: Feedforward neurons;
Py: Pyramidal cells;
Fb: Feedback neurons;
Parameters can be changed inside SetCortexParam.m and information
about these parameters are found in the publications cited above.

Other functions for the objects are:

Frequency(M2);  Calculates the average frequency for each cell in the
                network
Sparseness(M2); Calculates the level of sparseness
Coherence(M2);  Calculates the level of coherence