<p>
<b> Single pi (e=50GeV) scan of ECAL+HCAL </b>
<p>

<font color="#aa0000"> >>> click on thumbnail to get a full-size plot </font>

<p>
<hr>
<br>
Simulation details: <br>
 - 50k events
 - no field <br>
 - no vertex smearing <br>
 - no material in front of Calorimetry <br>
 - no Triggers <br>
<br>
<hr>


Plotted in this row is the HCAL subdetectors RecHits energy (mostly noise with single-pion contributin on the right side <br> 
Plots ====>>> HB >>> HE >>> HO >>> HF

<p>

<A HREF= "RecHits_energy_HB.gif">
<img src="RecHits_energy_HB.gif" width=24%> </A>
<A HREF= "RecHits_energy_HE.gif">
<img src="RecHits_energy_HE.gif" width=24%> </A>
<A HREF= "RecHits_energy_HO.gif">
<img src="RecHits_energy_HO.gif" width=24%> </A>
<A HREF= "RecHits_energy_HF.gif">
<img src="RecHits_energy_HF.gif" width=24%> </A>


<hr>

Plotted in this row is the RecHits timing vs energy <br>
Plots ====>>> HB >>> HE >>> HF 

<p>
<A HREF= "RecHitsTask_timing_vs_energy_profile_HB.gif">
<img src="RecHitsTask_timing_vs_energy_profile_HB.gif" width=24%> </A>
<A HREF= "RecHitsTask_timing_vs_energy_profile_HE.gif">
<img src="RecHitsTask_timing_vs_energy_profile_HE.gif" width=24%> </A>
<A HREF= "RecHitsTask_timing_vs_energy_profile_HF.gif">
<img src="RecHitsTask_timing_vs_energy_profile_HF.gif" width=24%> </A>


<hr>

Plotted in this row is the RecHits energy collected in R=0.3 : EM, HAD and EM+HAD without HO <br> 
Plots ====>>> Em component >>> Had component (no HO) >>> Total energy (Em+Had) >>>

<p>

<A HREF= "RecHitsTask_emean_vs_ieta_E.gif">
<img src="RecHitsTask_emean_vs_ieta_E.gif" width=24%> </A>
<A HREF= "RecHitsTask_emean_vs_ieta_H.gif">
<img src="RecHitsTask_emean_vs_ieta_H.gif" width=24%> </A>
<A HREF= "RecHitsTask_emean_vs_ieta_EH.gif">
<img src="RecHitsTask_emean_vs_ieta_EH.gif" width=24%> </A>


<hr>

Plotted in the row below is the same quantities as above, <br>
but for ONE single CaloTower in the event, which is closest to the MC direction of the pion
<br>
Also added is a plot of the single CaloTower scale ratio (new / previous).
The idea is to avoid any clustering (which masks transition) and to pinpoint each ieta ring in the most simplified and clear way.  
<p>
<pre>
Plots ====>>> Em component >>> Had component (no HO) >>> Total energy (Em+Had) >>> Ratio of total energy
</pre>


<p>
<A HREF= "CaloTowersTask_emean_vs_ieta_E1.gif">
<img src="CaloTowersTask_emean_vs_ieta_E1.gif" width=24%> </A>
<A HREF= "CaloTowersTask_emean_vs_ieta_H1.gif">
<img src="CaloTowersTask_emean_vs_ieta_H1.gif" width=24%> </A>
<A HREF= "CaloTowersTask_emean_vs_ieta_EH1.gif">
<img src="CaloTowersTask_emean_vs_ieta_EH1.gif" width=24%> </A>
<A HREF= "Ratio.gif">
<img src="Ratio.gif" width=24%> </A>
<hr>

Finally: number of CaloTowers in HB, HE and HF (separately) <br>

<p>
<A HREF= "N_calotowers_HB.gif">
<img src="N_calotowers_HB.gif" width=32%> </A>
<A HREF= "N_calotowers_HE.gif">
<img src="N_calotowers_HE.gif" width=32%> </A>
<A HREF= "N_calotowers_HF.gif">
<img src="N_calotowers_HF.gif" width=32%> </A>
<hr>
