310 likes | 414 Views
Electron selection parameterization PS, 1-10GeV Hlushchenko Olena 28.03.2014. Starting point. The analysis used PS runs in the momentum range 1-10GeV Both C herenkov A and B used CO2 gas One should take into account that SPS runs at 10-300GeV runs used He instead of CO2
E N D
Electron selection parameterizationPS, 1-10GeVHlushchenkoOlena28.03.2014
Starting point • The analysis used PS runs in the momentum range 1-10GeV • Both Cherenkov A and B used CO2 gas • One should take into account that SPS runs at 10-300GeV runs used He instead of CO2 • Pressure settings of the two Cherenkov detectors for all mixed hadron runs:
Electron selection • For the next selection only runs with Cherenkov pressure BELOW the pion threshold were taken. Some 9 GeV runs where excluded from the analysis as the pressure in the CherenkovB detector was set up above the pion threshold . • Runs with issues have been removed: 660067,660068,660052. • For all runs efficiency to have an electron signal with Cherenkov is close too 100% • For runs >=4 GeV the pressure in Cherenkov isabove the muon threshold pressure. • For those runs, when both Cherenkov detectors have a signal we can expect both: electrons and muons • In 1-3 GeV runs the pressure in Cherenkov is lower that is needed to have a muon signal-> pure electron selection [1] /// • In 3-10GeV runs we can have a rather pure electron selection by asking for two Cherenkov signals combined with the density cuts in the calorimeter
For the 4GeV runs we can also have some separation for muons, as the CherenkovA pressure is set up below the threshold pressure for 4 GeVmuonsand CherenkovB above.
DHCAL identification criteria, used to reject muons • There were no interaction layer • Number of center of all hits in the z-axis is higher then 20 • Total number of hits in all layers is bigger than 30 • There were no electronics errors (box events etc.) • The number of hits in last 4 layers is higher than 2 • There were less than 2,25 hits per layer in average
DHCAL identification criteria, used to select electrons • First interaction was in second or third layer • Number of center of all hits in the z-axis is less than 10 • There were no electronics errors (box events etc.) • The number of hits in last 12 layers is less than 2 • There were more than 5 hits per layer on average
DHCAL identification criteria, used to reject pions • First interaction was between 1st and 11th layer • Number of hits is greater than 20 • There were no electronics errors (box events etc.) • There were more than 2 hits per layer on average
Efficiency • This all gives us better approach for calculation of the efficiency to identify an electron with DHCAL on tag-and-probe method • Tag—pure electron events obtained with Cherenkov+muon cuts • Probe—reconstructed with DHCAL electrons • Efficiency: • - Number of the DHCAL and CHERENKOV signals for electrons • - Number for CHERENKOV signals for electron • PURITYCHEKFORCHERENKOVWITHDHCAL.PY
Electrons amount visualisation • For the next slides one can see, that starting from the 4 GeV beam momentum runs the muons start to apper in the Cherencov selection, where both Cherenkov worked. • Code ELECTRONIAYERNUMdISTRIB.PY
1 GeV run • No separation in 1st row • Pure electrons on 2nd row • Muons and pions in 3rd row
2GeV run • No separation in 1st row • Pure electrons on 2nd row • Muons and pions in 3rd row
3GeV run • No separation in 1st row • Pure electrons on 2nd row • Muons and pions in 3rd row
4 GeV run • No separation in 1st row • Electrons with muons on 2nd row • Muons and pions in 3rd row
5 GeV run • No separation in 1st row • Electrons with muons on 2nd row • Muons and pions in 3rd row
6 GeV run • No separation in 1st row • Electrons with muons on 2nd row • Muons and pions in 3rd row
7 GeV run • No separation in 1st row • Electrons with muons on 2nd row • Muons and pions in 3rd row
8GeV run • No separation in 1st row • Electrons with muons on 2nd row • Muons and pions in 3rd row
9GeV run • No separation in 1st row • Electrons with muons on 2nd row • Muons and pions in 3rd row
10 GeV run • No separation in 1st row • Electrons with muons on 2nd row • Muons and pions in 3rd row
Electron calibrated number of hits • The center points of energy branch for electrons peak. • Error bars is RMS • 2 fits: exponential, line • Electrons were selected with both Cherenkov gigving signal, layerEnergy>0 and hitDensity>2,25 for energies above 4 GeV (included)
Electron calibrated fitted number of hits • The center points of energy branch for electrons peak fitted with gaus. • Sigma error bars • 2 fits: exponential, line • Electrons were selected with both Cherenkov gigving signal, layerEnergy>0 and hitDensity>2,25 for energies above 4 GeV (included)
Electron hitdensity • The center points of hitDensity for electrons peak. • Error bars RMS • 3 fits: polynomial, exponential, line, fitted with Chi-squared • Electrons were selected with both Cherenkov gigving signal, layerEnergy>0 and hitDensity>2,25 for energies above 4 GeV (included)
What should be next • Repeating the analysis in higher energies for other particles • Parameterize the cuts for DHCAL selections for other types of particles • Check how did the efficiency change with applying new cuts
Appendix 1. Runs with pressure above the thredhold for Pions with Cherenkov detectors energy = 9.0 number=660054 A Presure(real)-Presure(thed)=1.71 B Presure(real)-Presure(thed)=1.71 number=660055 A Presure(real)-Presure(thed)=1.71 B Presure(real)-Presure(thed)=1.71 number=660056 A Presure(real)-Presure(thed)=1.71 B Presure(real)-Presure(thed)=1.71 number=660057 A Presure(real)-Presure(thed)=1.71 B Presure(real)-Presure(thed)=1.71 number=660058 A Presure(real)-Presure(thed)=1.71 B Presure(real)-Presure(thed)=1.71 number=660059 A Presure(real)-Presure(thed)=1.71 B Presure(real)-Presure(thed)=1.71 number=660060 A Presure(real)-Presure(thed)=1.71 B Presure(real)-Presure(thed)=1.71 number=660061 A Presure(real)-Presure(thed)=1.71 B Presure(real)-Presure(thed)=1.71 number=660066 A Presure(real)-Presure(thed)=1.96 B Presure(real)-Presure(thed)=1.96 number=660070 A Presure(real)-Presure(thed)=0.71 B Presure(real)-Presure(thed)=0.71 Runnumbers = [660054 , 660055, 660056, 660057 660058, 660059, 660060, 660061, 660066, 660070]
Appendix 2. Efficiency list • [0.004552106380884963, 0.07819047516454633, 0.18142964887287355, 0.2849973634925541, 0.20510182731203794, 0.25447461698450363, 0.29051572327044023, 0.2813719692489651, 0.25024173274028233, 0.2650893288266538]