<StimulationSet>

    <!-- the injectionSite must match a regionLabel defined in the morphology -->
   <!-- 20 Glu and Ca pulses / sec for 1 sec * 20 trains = 400 pulses -->
<!--first burst of Da - full strength -->
<!-- pulse 1 to produce 0.5 uM alone -->
<InjectionStim specieID="DaDbuf" injectionSite="sa1[0].pointA">
        <onset>               30000          </onset>
        <duration>            15              </duration>
        <rate>                18            </rate>
        <period>              50             </period>
        <end>                 30020          </end>
</InjectionStim>

<!-- pulse 2 = half of pulse 1 -->
<InjectionStim specieID="DaDbuf" injectionSite="sa1[0].pointA">
        <onset>               30050          </onset>
        <duration>            15              </duration>
        <rate>                9.0           </rate>
        <period>              50             </period>
        <end>                 30070          </end>
</InjectionStim>

<!-- pulse 3 = half of pulse 2  -->
<InjectionStim specieID="DaDbuf" injectionSite="sa1[0].pointA">
        <onset>               30100          </onset>
        <duration>            15              </duration>
        <rate>                6.3           </rate>
        <period>              50             </period>
        <end>                 30120          </end>
</InjectionStim>

<!-- pulses 4-20 = 2/3 of ctrl Da, to produce smaller peak than theta, because 40 pulses of 20hz are needed to produce same amount as 40 theta pulses-->
<InjectionStim specieID="DaDbuf" injectionSite="sa1[0].pointA">
        <onset>               30150          </onset>
        <duration>            15              </duration>
        <rate>                3           </rate>
        <period>              50             </period>
        <end>                 30999          </end>
</InjectionStim>

<!--later bursts of Da - partial strength (2/3d) -->
<!-- pulse 1 to produce 0.5 uM alone -->
<InjectionStim specieID="DaDbuf" injectionSite="sa1[0].pointA">
        <onset>               40000          </onset>
        <duration>            15              </duration>
        <rate>                12            </rate>
        <period>              50             </period>
        <end>                 40020          </end>
</InjectionStim>

<!-- pulse 2 = half of pulse 1 -->
<InjectionStim specieID="DaDbuf" injectionSite="sa1[0].pointA">
        <onset>               40050          </onset>
        <duration>            15              </duration>
        <rate>                6.0           </rate>
        <period>              50             </period>
        <end>                 40070          </end>
</InjectionStim>

<!-- pulse 3 = half of pulse 2  -->
<InjectionStim specieID="DaDbuf" injectionSite="sa1[0].pointA">
        <onset>               40100          </onset>
        <duration>            15              </duration>
        <rate>                4.2           </rate>
        <period>              50             </period>
        <end>                 40120          </end>
</InjectionStim>

<!-- pulses 4-20 = 2/3 of ctrl Da, faster injection to decrease stochasticity-->
<InjectionStim specieID="DaDbuf" injectionSite="sa1[0].pointA">
        <onset>               40150          </onset>
        <duration>            15              </duration>
        <rate>                2           </rate>
        <period>              50             </period>
        <end>                 40999          </end>
</InjectionStim>

<!-- third and subsequent trains even smaller -->
<InjectionStim specieID="DaDbuf" injectionSite="sa1[0].pointA">
        <onset>               50000          </onset>
        <duration>            15              </duration>
        <rate>                2.0            </rate>
        <period>              50             </period>
        <end>                 50999          </end>
</InjectionStim>

<!-- fourth and subsequent trains smaller still -->
<InjectionStim specieID="DaDbuf" injectionSite="sa1[0].pointA">
        <onset>               60000          </onset>
        <duration>            15              </duration>
        <rate>                1.0            </rate>
        <period>              50             </period>
        <end>                 60999          </end>
        <numTrains>           5             </numTrains>
        <interTrainInterval>  9001          </interTrainInterval>   
</InjectionStim>

<!-- basal ACh value prior to and after stim  -->
<InjectionStim specieID="ACh" injectionSite="dend:submembrane">
        <onset>               0          </onset>
        <duration>            30000      </duration>
        <rate>                0.2         </rate>
</InjectionStim>
<InjectionStim specieID="ACh" injectionSite="dend:submembrane">
        <onset>               110000          </onset>
        <duration>            200000      </duration>
        <rate>                0.2         </rate>
</InjectionStim>
<!-- basal ACh stim paused during Da stim-->
<InjectionStim specieID="ACh" injectionSite="dend:submembrane">
        <onset>               33200          </onset>
        <duration>             6800      </duration>
        <rate>                0.2         </rate>
        <period>              10000             </period>
        <end>                 110000          </end>
</InjectionStim>
<!-- burst during stim : Assume that ACh neurons can follow 20 Hz -->
<InjectionStim specieID="ACh" injectionSite="sa1[0].pointA">
        <onset>               30000          </onset>
        <duration>            15              </duration>
        <rate>                3.0           </rate>
        <period>              50             </period>
        <end>                 30999          </end>
        <numTrains>           8             </numTrains>
        <interTrainInterval>  9001          </interTrainInterval>  
</InjectionStim>

<!-- Should this stimulation be eliminated???-->
<InjectionStim specieID="ACh" injectionSite="dend:submembrane">
        <onset>               30000          </onset>
        <duration>            15              </duration>
        <rate>                3.0           </rate>
        <period>              50             </period>
        <end>                 30999          </end>
        <numTrains>           8             </numTrains>
        <interTrainInterval>  9001          </interTrainInterval>  
</InjectionStim>

<InjectionStim specieID="GluGbuf" injectionSite="sa1[0].pointA">
        <onset>               30000          </onset>
        <duration>            5              </duration>
        <rate>                1.125           </rate>
        <period>              50             </period>
        <end>                 30999          </end>
        <numTrains>           8             </numTrains>
        <interTrainInterval>  9001          </interTrainInterval>  
</InjectionStim>

<InjectionStim specieID="Ca" injectionSite="sa1[0].pointA">
        <onset>               30000          </onset>
        <duration>            3              </duration>
        <rate>                312            </rate>
        <period>              50             </period>
        <end>                 30999          </end>
        <numTrains>           8             </numTrains>
        <interTrainInterval>  9001          </interTrainInterval>  
</InjectionStim>

<InjectionStim specieID="Ca" injectionSite="dend:submembrane">
        <onset>               30000          </onset>
        <duration>            3              </duration>
        <rate>                156            </rate>
        <period>              50             </period>
        <end>                 30999          </end>
        <numTrains>           8             </numTrains>
        <interTrainInterval>  9001          </interTrainInterval>  
</InjectionStim>


</StimulationSet>