<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>               20000          </onset>
        <duration>            15              </duration>
        <rate>                14.4            </rate>
        <period>              50             </period>
        <end>                 20020          </end>
</InjectionStim>

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

<!-- pulse 3 = half of pulse 2  -->
<InjectionStim specieID="DaDbuf" injectionSite="sa1[0].pointA">
        <onset>               20100          </onset>
        <duration>            15              </duration>
        <rate>                5.04           </rate>
        <period>              50             </period>
        <end>                 20120          </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>               20150          </onset>
        <duration>            15              </duration>
        <rate>                2.4           </rate>
        <period>              50             </period>
        <end>                 20999          </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>               30000          </onset>
        <duration>            15              </duration>
        <rate>                9.6            </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>                4.8           </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>                3.36           </rate>
        <period>              50             </period>
        <end>                 30120          </end>
</InjectionStim>

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

<!-- third and subsequent trains even smaller -->
<InjectionStim specieID="DaDbuf" injectionSite="sa1[0].pointA">
        <onset>               40000          </onset>
        <duration>            15              </duration>
        <rate>                1.6            </rate>
        <period>              50             </period>
        <end>                 40999          </end>
        <numTrains>           18             </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>