<StimulationSet>

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

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

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

<!-- pulses 4-20 = half of pulse 3, faster injection to decrease stochasticity-->
<InjectionStim specieID="DaDbuf" injectionSite="sa1[0].pointA">
        <onset>               21450          </onset>
        <duration>            15              </duration>
        <rate>                3           </rate>
        <period>              50             </period>
        <end>                 22299          </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>               31300          </onset>
        <duration>            15              </duration>
        <rate>                12            </rate>
        <period>              50             </period>
        <end>                 31320          </end>
</InjectionStim>

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

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

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

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

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

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

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

 <!-- rebound burst: develops more slowly -->
<InjectionStim specieID="ACh" injectionSite="sa1[0].pointA">
        <onset>               22400          </onset>
        <duration>            800            </duration>
        <rate>                0            </rate>
        <period>              900           </period>
        <end>                 23299          </end>
        <numTrains>           20             </numTrains>
        <interTrainInterval>  9101          </interTrainInterval>  
</InjectionStim>
<!-- Train 1 -->
<InjectionStim specieID="GluGbuf" injectionSite="sa1[0].pointA">
        <onset>               20000          </onset>
        <duration>            5              </duration>
        <rate>                1.125           </rate>
        <period>              50             </period>
        <end>                 20999          </end>
        <numTrains>           20             </numTrains>
        <interTrainInterval>  9001          </interTrainInterval>   
</InjectionStim>

<InjectionStim specieID="Ca" injectionSite="sa1[0].pointA">
        <onset>               20000          </onset>
        <duration>            3              </duration>
        <rate>                156            </rate>
        <period>              50             </period>
        <end>                 20999          </end>
        <numTrains>           20             </numTrains>
        <interTrainInterval>  9001          </interTrainInterval>  
</InjectionStim>

<InjectionStim specieID="Ca" injectionSite="dend:submembrane">
        <onset>               20000          </onset>
        <duration>            3              </duration>
        <rate>                156            </rate>
        <period>              50             </period>
        <end>                 20999          </end>
        <numTrains>           20             </numTrains>
        <interTrainInterval>  9001          </interTrainInterval>  
</InjectionStim>


</StimulationSet>