ATLAS Offline Software
Loading...
Searching...
No Matches
DCubeHistograms.py
Go to the documentation of this file.
1# Copyright (C) 2002-2025 CERN for the benefit of the ATLAS collaboration
2
3# this script can be used to create DCube histograms from the output ntuples of NSWPRDValAlg
4
5import os, sys, ROOT
6
8
9
11 truthPtHist = ROOT.TH1F("truthPt","truthPt;p_{T} [GeV]",100,0,100)
12 truthEtaHist = ROOT.TH1F("truthEta","truthEta;#eta",100,-2.8,2.8)
13 truthPhiHist = ROOT.TH1F("truthPhi","truthPhi;#phi",100,-3.15,3.15)
14
16 rpcLocalX = ROOT.TH1F("rpcLocalX","rpcLocalX;RPC_SIM_hitLocalPositionX",100,-1.05,1.05)
17 rpcLocalY = ROOT.TH1F("rpcLocalY","rpcLocalY;RPC_SIM_hitLocalPositionY",100,-1500,1500)
18 rpcStationEta = ROOT.TH1F("rpcStationEta","rpcStationEta;RPC_SIM_stationEta",16,-8,8)
19 rpcStationPhi = ROOT.TH1F("rpcStationPhi","rpcStationPhi;RPC_SIM_stationPhi",9,0,9)
20 rpcGlobalX = ROOT.TH1F("rpcGlobalX","rpcGlobalX;RPC_SIM_hitGlobalPositionX",100,-15000,15000)
21 rpcGlobalY = ROOT.TH1F("rpcGlobalY","rpcGlobalY;RPC_SIM_hitGlobalPositionY",100,-15000,15000)
22 rpcGlobalZ = ROOT.TH1F("rpcGlobalZ","rpcGlobalZ;RPC_SIM_hitGlobalPositionZ",100,-15000,15000)
23 rpcGasGap = ROOT.TH1F("rpcGasGap","rpcGasGap;RPC_SIM_GasGap",4,0,4)
24
25 # RPC Digit Histograms (all stations)
26 rpcDigitlocalX = ROOT.TH1F("rpcDigitlocalX","rpcDigitlocalX;Digits_RPC_localPosX [mm]",100,-1200,1200)
27 rpcDigitlocalY = ROOT.TH1F("rpcDigitlocalY","rpcDigitlocalY;Digits_RPC_localPosY [mm]",100,-1200,1200)
28 rpcDigitStationEta = ROOT.TH1F("rpcDigitStationEta","rpcDigitStationEta;Digits_RPC_stationEta",20,-10,10)
29 rpcDigitStationPhi = ROOT.TH1F("rpcDigitStationPhi","rpcDigitStationPhi;Digits_RPC_stationPhi",10,0,10)
30 rpcDigitglobalX = ROOT.TH1F("rpcDigitglobalX","rpcDigitglobalX;Digits_RPC_globalPosX [mm]",100,-13000,12000)
31 rpcDigitglobalY = ROOT.TH1F("rpcDigitglobalY","rpcDigitglobalY;Digits_RPC_globalPosY [mm]",100,-15000,15000)
32 rpcDigitglobalZ = ROOT.TH1F("rpcDigitglobalZ","rpcDigitglobalZ;Digits_RPC_globalPosZ [mm]",100,-15000,15000)
33 rpcDigitGasGap = ROOT.TH1F("rpcDigitGasGap","rpcDigitGasGap;Digits_RPC_GasGap",5,0,5)
34 rpcDigitdoubletR = ROOT.TH1F("rpcDigitdoubletR","rpcDigitdoubletR;Digits_RPC_doubletR",4,0,4)
35 rpcDigitdoubletZ = ROOT.TH1F("rpcDigitdoubletZ","rpcDigitdoubletZ;Digits_RPC_doubletZ",4,0,4)
36 rpcDigitdoubletPhi = ROOT.TH1F("rpcDigitdoubletPhi","rpcDigitdoubletPhi;Digits_RPC_doubletPhi",4,0,4)
37 rpcDigitTime = ROOT.TH1F("rpcDigitTime","rpcDigitTime;Digits_RPC_time",100,0,150)
38 rpcDigitStrip = ROOT.TH1F("rpcDigitStrip","rpcDigitStrip;Digits_RPC_strip",10,0,80)
39 rpcDigitStripNumber = ROOT.TH1F("rpcDigitStripNumber","rpcDigitStripNumber;Digits_RPC_stripNumber",10,0,80)
40 rpcDigitMeasuresPhi = ROOT.TH1F("rpcDigitMeasuresPhi","rpcDigitMeasuresPhi;Digits_RPC_measuresPhi",2,0,2)
41 rpcDigitglobalXY = ROOT.TH2F("rpcDigitglobalXY","rpcDigitglobalXY;Digits_RPC_globalPosX [mm];Digits_RPC_globalPosY [mm]",100,-15000,15000,100,-15000,15000)
42 rpcDigitglobalZY = ROOT.TH2F("rpcDigitglobalZY","rpcDigitglobalZY;Digits_RPC_globalPosZ [mm];Digits_RPC_globalPosY [mm]",100,-15000,15000,100,-15000,15000)
43 rpcDigitglobalZX = ROOT.TH2F("rpcDigitglobalZX","rpcDigitglobalZX;Digits_RPC_globalPosZ [mm];Digits_RPC_globalPosX [mm]",100,-15000,15000,100,-15000,15000)
44
45 # RPC Digit Histograms (BIS RPCs only)
46 bis_rpcDigitlocalX = ROOT.TH1F("bis_rpcDigitlocalX","bis_rpcDigitlocalX;Digits_bis_RPC_localPosX [mm]",100,-1000,1000)
47 bis_rpcDigitlocalY = ROOT.TH1F("bis_rpcDigitlocalY","bis_rpcDigitlocalY;Digits_bis_RPC_localPosY [mm]",100,-1000,1000)
48 bis7_rpcDigitlocalX = ROOT.TH1F("bis7_rpcDigitlocalX","bis7_rpcDigitlocalX;Digits_bis7_RPC_localPosX [mm]",100,-1000,1000)
49 bis7_rpcDigitlocalY = ROOT.TH1F("bis7_rpcDigitlocalY","bis7_rpcDigitlocalY;Digits_bis7_RPC_localPosY [mm]",100,-1000,1000)
50 bis8_rpcDigitlocalX = ROOT.TH1F("bis8_rpcDigitlocalX","bis8_rpcDigitlocalX;Digits_bis8_RPC_localPosX [mm]",100,-1000,1000)
51 bis8_rpcDigitlocalY = ROOT.TH1F("bis8_rpcDigitlocalY","bis8_rpcDigitlocalY;Digits_bis8_RPC_localPosY [mm]",100,-1000,1000)
52 bis_rpcDigitStationEta = ROOT.TH1F("bis_rpcDigitStationEta","bis_rpcDigitStationEta;Digits_bis_RPC_stationEta",20,-10,10)
53 bis_rpcDigitStationPhi = ROOT.TH1F("bis_rpcDigitStationPhi","bis_rpcDigitStationPhi;Digits_bis_RPC_stationPhi",10,0,10)
54 bis_rpcDigitglobalX = ROOT.TH1F("bis_rpcDigitglobalX","bis_rpcDigitglobalX;Digits_bis_RPC_globalPosX [mm]",100,-5000,5000)
55 bis_rpcDigitglobalY = ROOT.TH1F("bis_rpcDigitglobalY","bis_rpcDigitglobalY;Digits_bis_RPC_globalPosY [mm]",100,-5000,5000)
56 bis_rpcDigitglobalZ = ROOT.TH1F("bis_rpcDigitglobalZ","bis_rpcDigitglobalZ;Digits_bis_RPC_globalPosZ [mm]",100,-10000,10000)
57 bis_rpcDigitGasGap = ROOT.TH1F("bis_rpcDigitGasGap","bis_rpcDigitGasGap;Digits_bis_RPC_GasGap",5,0,5)
58 bis_rpcDigitdoubletR = ROOT.TH1F("bis_rpcDigitdoubletR","bis_rpcDigitdoubletR;Digits_bis_RPC_doubletR",4,0,4)
59 bis_rpcDigitdoubletZ = ROOT.TH1F("bis_rpcDigitdoubletZ","bis_rpcDigitdoubletZ;Digits_bis_RPC_doubletZ",4,0,4)
60 bis_rpcDigitdoubletPhi = ROOT.TH1F("bis_rpcDigitdoubletPhi","bis_rpcDigitdoubletPhi;Digits_bis_RPC_doubletPhi",4,0,4)
61 bis_rpcDigitTime = ROOT.TH1F("bis_rpcDigitTime","bis_rpcDigitTime;Digits_bis_RPC_time [ns]",80,0,40)
62 bis7_rpcDigitTime = ROOT.TH1F("bis7_rpcDigitTime","bis7_rpcDigitTime;Digits_bis7_RPC_time [ns]",80,0,40)
63 bis8_rpcDigitTime = ROOT.TH1F("bis8_rpcDigitTime","bis8_rpcDigitTime;Digits_bis8_RPC_time [ns]",80,0,40)
64 bis_rpcDigitStrip = ROOT.TH1F("bis_rpcDigitStrip","bis_rpcDigitStrip;Digits_bis_RPC_strip",10,0,10)
65 bis_rpcDigitStripNumber = ROOT.TH1F("bis_rpcDigitStripNumber","bis_rpcDigitStripNumber;Digits_bis_RPC_stripNumber",10,0,80)
66 bis_rpcDigitMeasuresPhi = ROOT.TH1F("bis_rpcDigitMeasuresPhi","bis_rpcDigitMeasuresPhi;Digits_bis_RPC_measuresPhi",2,0,2)
67 bis_rpcDigitglobalXY = ROOT.TH2F("bis_rpcDigitglobalXY","bis_rpcDigitglobalXY;Digits_bis_RPC_globalPosX [mm];Digits_bis_RPC_globalPosY [mm]",100,-5000,5000,100,-5000,5000)
68 bis_rpcDigitglobalZY = ROOT.TH2F("bis_rpcDigitglobalZY","bis_rpcDigitglobalZY;Digits_bis_RPC_globalPosZ [mm];Digits_bis_RPC_globalPosY [mm]",100,-15000,15000,100,-5000,5000)
69 bis_rpcDigitglobalZX = ROOT.TH2F("bis_rpcDigitglobalZX","bis_rpcDigitglobalZX;Digits_bis_RPC_globalPosZ [mm];Digits_bis_RPC_globalPosX [mm]",100,-15000,15000,100,-5000,5000)
70
71
72
74 mdtLocalX = ROOT.TH1F("mdtLocalX","mdtLocalX;MDT_Sim_hitLocalPositionX",100,-15,15)
75 mdtLocalY = ROOT.TH1F("mdtLocalY","mdtLocalY;MDT_Sim_hitLocalPositionY",100,-16,16)
76 mdtLocalZ = ROOT.TH1F("mdtLocalZ","mdtLocalZ;MDT_Sim_hitLocalPositionZ",100,-2500,2500)
77 mdtStationEta = ROOT.TH1F("mdtStationEta","mdtStationEta;MDT_Sim_stationEta",16,-8,8)
78 mdtStationPhi = ROOT.TH1F("mdtStationPhi","mdtStationPhi;MDT_Sim_stationPhi",9,0,9)
79 mdtGlobalX = ROOT.TH1F("mdtGlobalX","mdtGlobalX;MDT_Sim_hitGlobalPositionX",100,-15000,15000)
80 mdtGlobalY = ROOT.TH1F("mdtGlobalY","mdtGlobalY;MDT_Sim_hitGlobalPositionY",100,-15000,15000)
81 mdtGlobalZ = ROOT.TH1F("mdtGlobalZ","mdtGlobalZ;MDT_Sim_hitGlobalPositionZ",100,-26000,26000)
82 mdtTube = ROOT.TH1F("mdtTube","mdtTube;MDT_Sim_tube",100,0,100)
83
84 # MDT Digit Histograms (all stations)
85 mdtDigitStationEta = ROOT.TH1F("mdtDigitStationEta","mdtDigitStationEta;Digits_MDT_stationEta",20,-10,10)
86 mdtDigitStationPhi = ROOT.TH1F("mdtDigitStationPhi","mdtDigitStationPhi;Digits_MDT_stationPhi",10,0,10)
87 mdtDigitglobalX = ROOT.TH1F("mdtDigitglobalX","mdtDigitglobalX;Digits_MDT_globalPosX [mm]",100,-15000,15000)
88 mdtDigitglobalY = ROOT.TH1F("mdtDigitglobalY","mdtDigitglobalY;Digits_MDT_globalPosY [mm]",100,-15000,15000)
89 mdtDigitglobalZ = ROOT.TH1F("mdtDigitglobalZ","mdtDigitglobalZ;Digits_MDT_globalPosZ [mm]",100,-25000,25000)
90 mdtDigitNumberOfMultilayers = ROOT.TH1F("mdtDigitNumberOfMultilayers","mdtDigitNumberOfMultilayers;Digits_MDT_numberOfMultilayers",4,0,4)
91 mdtDigitTube = ROOT.TH1F("mdtDigitTube","mdtDigitTube;Digits_MDT_tube",100,0,100)
92 mdtDigitTubeLayer = ROOT.TH1F("mdtDigitTubeLayer","mdtDigitTubeLayer;Digits_MDT_tubeLayer",6,0,6)
93 mdtDigitMultilayer = ROOT.TH1F("mdtDigitMultilayer","mdtDigitMultilayer;Digits_MDT_Multilayer",4,0,4)
94 mdtDigitTime = ROOT.TH1F("mdtDigitTime","mdtDigitTime;Digits_MDT_time [ns]",100,0,2000)
95 mdtDigitCharge = ROOT.TH1F("mdtDigitCharge","mdtDigitCharge;Digits_MDT_charge [e]",100,0,600)
96 mdtDigitglobalXY = ROOT.TH2F("mdtDigitglobalXY","mdtDigitglobalXY;Digits_MDT_globalPosX [mm];Digits_MDT_globalPosY [mm]",100,-15000,15000,100,-15000,15000)
97 mdtDigitglobalZY = ROOT.TH2F("mdtDigitglobalZY","mdtDigitglobalZY;Digits_MDT_globalPosZ [mm];Digits_MDT_globalPosY [mm]",100,-15000,15000,100,-15000,15000)
98 mdtDigitglobalZX = ROOT.TH2F("mdtDigitglobalZX","mdtDigitglobalZX;Digits_MDT_globalPosZ [mm];Digits_MDT_globalPosX [mm]",100,-15000,15000,100,-15000,15000)
99 mdtDigitlocalTubePosX = ROOT.TH1F("mdtDigitlocalTubePosX","mdtDigitlocalTubePosX;Digits_MDT_localTubePosX [mm]",100,-100,500)
100 mdtDigitlocalTubePosY = ROOT.TH1F("mdtDigitlocalTubePosY","mdtDigitlocalTubePosY;Digits_MDT_localTubePosY [mm]",100,-100,500)
101 mdtDigitlocalTubePosZ = ROOT.TH1F("mdtDigitlocalTubePosZ","mdtDigitlocalTubePosZ;Digits_MDT_localTubePosZ [mm]",100,-1200,1200)
102
103 # MDT Digit Histograms (BIS mdtDigits only)
104 bis_mdtDigitStationEta = ROOT.TH1F("bis_mdtDigitStationEta","bis_mdtDigitStationEta;Digits_bis_MDT_stationEta",20,-10,10)
105 bis_mdtDigitStationPhi = ROOT.TH1F("bis_mdtDigitStationPhi","bis_mdtDigitStationPhi;Digits_bis_MDT_stationPhi",10,0,10)
106 bis_mdtDigitglobalX = ROOT.TH1F("bis_mdtDigitglobalX","bis_mdtDigitglobalX;Digits_bis_MDT_globalPosX [mm]",100,-5000,5000)
107 bis_mdtDigitglobalY = ROOT.TH1F("bis_mdtDigitglobalY","bis_mdtDigitglobalY;Digits_bis_MDT_globalPosY [mm]",100,-5000,5000)
108 bis_mdtDigitglobalZ = ROOT.TH1F("bis_mdtDigitglobalZ","bis_mdtDigitglobalZ;Digits_bis_MDT_globalPosZ [mm]",100,-10000,10000)
109 bis_mdtDigitNumberOfMultilayers = ROOT.TH1F("bis_mdtDigitNumberOfMultilayers","bis_mdtDigitNumberOfMultilayers;Digits_bis_MDT_numberOfMultilayers",4,0,4)
110 bis_mdtDigitTube = ROOT.TH1F("bis_mdtDigitTube","bis_mdtDigitTube;Digits_bis_MDT_tube",100,0,100)
111 bis_mdtDigitTubeLayer = ROOT.TH1F("bis_mdtDigitTubeLayer","bis_mdtDigitTubeLayer;Digits_bis_MDT_tubeLayer",6,0,6)
112 bis_mdtDigitMultilayer = ROOT.TH1F("bis_mdtDigitMultilayer","bis_mdtDigitMultilayer;Digits_bis_MDT_Multilayer",4,0,4)
113 bis_mdtDigitTime = ROOT.TH1F("bis_mdtDigitTime","bis_mdtDigitTime;Digits_bis_MDT_time [ns]",100,0,2000)
114 bis_mdtDigitCharge = ROOT.TH1F("bis_mdtDigitCharge","bis_mdtDigitCharge;Digits_bis_MDT_charge [e]",100,0,600)
115 bis_mdtDigitChannel = ROOT.TH1F("bis_mdtDigitChannel","bis_mdtDigitChannel;Digits_bis_MDT_channel",10,0,10)
116 bis_mdtDigitglobalXY = ROOT.TH2F("bis_mdtDigitglobalXY","bis_mdtDigitglobalXY;Digits_bis_MDT_globalPosX [mm];Digits_bis_MDT_globalPosY [mm]",100,-5000,5000,100,-5000,5000)
117 bis_mdtDigitglobalZY = ROOT.TH2F("bis_mdtDigitglobalZY","bis_mdtDigitglobalZY;Digits_bis_MDT_globalPosZ [mm];Digits_bis_MDT_globalPosY [mm]",100,-15000,15000,100,-5000,5000)
118 bis_mdtDigitglobalZX = ROOT.TH2F("bis_mdtDigitglobalZX","bis_mdtDigitglobalZX;Digits_bis_MDT_globalPosZ [mm];Digits_bis_MDT_globalPosX [mm]",100,-15000,15000,100,-5000,5000)
119 bis_mdtDigitlocalTubePosX = ROOT.TH1F("bis_mdtDigitlocalTubePosX","bis_mdtDigitlocalTubePosX;Digits_MDT_localTubePosX [mm]",120,-60,60)
120 bis_mdtDigitlocalTubePosY = ROOT.TH1F("bis_mdtDigitlocalTubePosY","bis_mdtDigitlocalTubePosY;Digits_MDT_localTubePosY [mm]",960,-900,900)
121 bis_mdtDigitlocalTubePosZ = ROOT.TH1F("bis_mdtDigitlocalTubePosZ","bis_mdtDigitlocalTubePosZ;Digits_MDT_localTubePosZ [mm]",960,-900,900)
122
123 bis_mdtDigitlocalTubePosXY = ROOT.TH2F("bis_mdtDigitlocalTubePosXY","bis_mdtDigitlocalTubePosXY;Digits_MDT_localTubePosX [mm];Digits_MDT_localTubePosY [mm]",60,-30,30,100,-100,100)
124 bis_mdtDigitlocalTubePosZY = ROOT.TH2F("bis_mdtDigitlocalTubePosZY","bis_mdtDigitlocalTubePosZY;Digits_MDT_localTubePosZ [mm];Digits_MDT_localTubePosY [mm]",100,-1200,1200,100,-100,100)
125 bis_mdtDigitlocalTubePosZX = ROOT.TH2F("bis_mdtDigitlocalTubePosZX","bis_mdtDigitlocalTubePosZX;Digits_MDT_localTubePosZ [mm];Digits_MDT_localTubePosX [mm]",48,1,180,60,-30,30)
126
127
128
130 cscGlobalX = ROOT.TH1F("cscGlobalX","cscGlobalX;CSC_Sim_hitGlobalPositionX",100,-1500,2200)
131 cscGlobalY = ROOT.TH1F("cscGlobalY","cscGlobalY;CSC_Sim_hitGlobalPositionY",100,-2000,2000)
132 cscGlobalZ = ROOT.TH1F("cscGlobalZ","cscGlobalZ;CSC_Sim_hitGlobalPositionZ",100,-8000,-8000)
133 cscWireLayer = ROOT.TH1F("cscWireLayer","cscWireLayer;CSC_Sim_wireLayer",5,0,5)
134 cscStrip = ROOT.TH1F("cscStrip","cscStrip;CSC_Sim_strip",3,0,3)
135 cscglobalTime = ROOT.TH1F("cscglobalTime","cscglobalTime;CSC_Sim_globalTime",100,0,60)
136 cscKineticEnergy = ROOT.TH1F("cscKineticEnergy","cscKineticEnergy;CSC_Sim_kineticEnergy",100,0,400000)
137 cscDepositEnergy = ROOT.TH1F("cscDepositEnergy","cscDeposityEnergy;CSC_Sim_depositEnergy",100,0,0.008)
138 cscSimStationEta = ROOT.TH1F("cscSimStationEta","cscSimStationEta;CSC_Sim_stationEta",4,-3,1)
139 cscSimStationPhi = ROOT.TH1F("cscSimStationPhi","cscSimStationPhi;CSC_Sim_stationPhi",9,0,9)
140 cscSimChamberLayer = ROOT.TH1F("cscSimChamberLayer","cscSimChamberLayer;CSC_Sim_chamberLayer",4,0,4)
141 # CSC Digit Histograms
142 CSCDigitStationEta = ROOT.TH1F("CSCDigitStationEta","CSCDigitStationEta;Digits_CSC_stationEta",4,-3,1)
143 CSCDigitStationPhi = ROOT.TH1F("CSCDigitStationPhi","CSCDigitStationPhi;Digits_CSC_stationPhi",9,0,9)
144 CSCDigitWireLayer = ROOT.TH1F("CSCDigitWireLayer","CSCDigitWireLayer;Digits_CSC_wireLayer",5,0,5)
145 CSCDigitlocalPosX = ROOT.TH1F("CSCDigitlocalPosX","CSCDigitlocalPosX;Digits_CSC_localPosX",100,0,100)
146 CSCDigitglobalPosX = ROOT.TH1F("CSCDigitglobalPosX","CSCDigitglobalPosX;Digits_CSC_globalPosX",100,-1600,1600)
147 CSCDigitglobalPosY = ROOT.TH1F("CSCDigitglobalPosY","CSCDigitglobalPosY;Digits_CSC_globalPosY",100,-1600,1600)
148 CSCDigitglobalPosZ = ROOT.TH1F("CSCDigitglobalPosZ","CSCDigitglobalPosZ;Digits_CSC_globalPosZ",100,-8000,-8000)
149 # CSC SDO Histograms
150 CSCSDOStationEta = ROOT.TH1F("CSCSDOStationEta","CSCSDOStationEta;SDO_CSC_stationEta",4,-3,1)
151 CSCSDOStationPhi = ROOT.TH1F("CSCSDOStationPhi","CSCSDOStationPhi;SDO_CSC_stationPhi",9,0,9)
152 CSCSDOChamberLayer = ROOT.TH1F("CSCSDOChamberLayer","CSCSDOChamberLayer;SDO_CSC_chamberLayer",4,0,4)
153 CSCSDOWireLayer = ROOT.TH1F("CSCSDOWireLayer","CSCSDOWireLayer;SDO_CSC_wireLayer",5,0,5)
154 CSCSDOStrip = ROOT.TH1F("CSCSDOStrip","CSCSDOStrip;SDO_CSC_strip",100,0,200)
155 CSCSDOBarcode = ROOT.TH1F("CSCSDOBarcode","CSCSDOBarcode;SDO_CSC_barcode",100,0,250000)
156 CSCSDOGlobalX = ROOT.TH1F("CSCSDOGlobalX","CSCSDOGlobalX;SDO_CSC_globalPosX",100,-2000,2000)
157 CSCSDOGlobalY = ROOT.TH1F("CSCSDOGlobalY","CSCSDOGlobalY;SDO_CSC_globalPosY",100,-2000,2000)
158 CSCSDOGlobalZ = ROOT.TH1F("CSCSDOGlobalZ","CSCSDOGlobalZ;SDO_CSC_globalPosZ",100,-8000,8000)
159 CSCSDOLocalX = ROOT.TH1F("CSCSDOLocalX","CSCSDOLocalX;SDO_CSC_localPosX",100,-1000,1000)
160 CSCSDOLocalY = ROOT.TH1F("CSCSDOLocalY","CSCSDOLocalY;SDO_CSC_localPosY",100,-1000,1000)
161 # CSC RDO Histograms
162 CSCRDOStationEta = ROOT.TH1F("CSCRDOStationEta","CSCRDOStationEta;RDO_CSC_stationEta",4,-3,1)
163 CSCRDOStationPhi = ROOT.TH1F("CSCRDOStationPhi","CSCRDOStationPhi;RDO_CSC_stationPhi",9,0,9)
164 CSCRDOChamberLayer = ROOT.TH1F("CSCRDOChamberLayer","CSCRDOChamberLayer;RDO_CSC_chamberLayer",4,0,4)
165 CSCRDOWireLayer = ROOT.TH1F("CSCRDOWireLayer","CSCRDOWireLayer;RDO_CSC_wireLayer",5,0,5)
166 CSCRDOStrip = ROOT.TH1F("CSCRDOStrip","CSCRDOStrip;RDO_CSC_strip",100,0,200)
167 CSCRDOmeasurePhi = ROOT.TH1F("CSCRDOmeasurePhi","CSCRDOmeasurePhi;RDO_CSC_measurePhi",2,0,2)
168 CSCRDOTime = ROOT.TH1F("CSCRDOTime","CSCRDOTime;RDO_CSC_time",4,-2,2)
169 CSCRDOGlobalX = ROOT.TH1F("CSCRDOGlobalX","CSCRDOGlobalX;RDO_CSC_globalPosX",100,-2500,2500)
170 CSCRDOGlobalY = ROOT.TH1F("CSCRDOGlobalY","CSCRDOGlobalY;RDO_CSC_globalPosY",100,-2500,2500)
171 CSCRDOGlobalZ = ROOT.TH1F("CSCRDOGlobalZ","CSCRDOGlobalZ;RDO_CSC_globalPosZ",100,-8000,8000)
172 CSCRDOLocalX = ROOT.TH1F("CSCRDOLocalX","CSCRDOLocalX;RDO_CSC_localPosX",100,0,100)
173 CSCRDOLocalY = ROOT.TH1F("CSCRDOLocalY","CSCRDOLocalY;RDO_CSC_localPosY",100,0,100)
174 # CSC PRD Histograms
175 CSCPRDStationEta = ROOT.TH1F("CSCPRDStationEta","CSCPRDStationEta;PRD_CSC_stationEta",4,-3,1)
176 CSCPRDStationPhi = ROOT.TH1F("CSCPRDStationPhi","CSCPRDStationPhi;PRD_CSC_stationPhi",9,0,9)
177 CSCPRDChamberLayer = ROOT.TH1F("CSCPRDChamberLayer","CSCPRDChamberLayer;PRD_CSC_chamberLayer",4,0,4)
178 CSCPRDWireLayer = ROOT.TH1F("CSCPRDWireLayer","CSCPRDWireLayer;PRD_CSC_wireLayer",5,0,5)
179 CSCPRDStrip = ROOT.TH1F("CSCPRDStrip","CSCPRDStrip;PRD_CSC_strip",100,0,200)
180 CSCPRDmeasurePhi = ROOT.TH1F("CSCPRDmeasurePhi","CSCPRDmeasurePhi;PRD_CSC_measurePhi",2,0,2)
181 CSCPRDTime = ROOT.TH1F("CSCPRDTime","CSCPRDTime;PRD_CSC_time",100,-100,1000)
182 CSCPRDCharge = ROOT.TH1F("CSCPRDCharge","CSCPRDCharge;PRD_CSC_charge",100,0,10000000)
183 CSCPRDGlobalX = ROOT.TH1F("CSCPRDGlobalX","CSCPRDGlobalX;PRD_CSC_globalPosX",100,-2500,2500)
184 CSCPRDGlobalY = ROOT.TH1F("CSCPRDGlobalY","CSCPRDGlobalY;PRD_CSC_globalPosY",100,-2500,2500)
185 CSCPRDGlobalZ = ROOT.TH1F("CSCPRDGlobalZ","CSCPRDGlobalZ;PRD_CSC_globalPosZ",100,-8000,8000)
186 CSCPRDLocalX = ROOT.TH1F("CSCPRDLocalX","CSCPRDLocalX;PRD_CSC_localPosX",100,-1000,1000)
187
189 tgcLocalX = ROOT.TH1F("tgcLocalX","tgcLocalX;TGC_Sim_hitLocalPositionX",100,-1.5,1.5)
190 tgcLocalY = ROOT.TH1F("tgcLocalY","tgcLocalY;TGC_Sim_hitLocalPositionY",100,-800,800)
191 tgcLocalZ = ROOT.TH1F("tgcLocalZ","tgcLocalZ;TGC_Sim_hitLocalPositionZ",100,-1200,1200)
192 tgcStationEta = ROOT.TH1F("tgcStationEta","tgcStationEta;TGC_Sim_stationEta",12,-6,6)
193 tgcStationPhi = ROOT.TH1F("tgcStationPhi","tgcStationPhi;TGC_Sim_stationPhi",50,0,50)
194 tgcGlobalX = ROOT.TH1F("tgcGlobalX","tgcGlobalX;TGC_Sim_hitGlobalPositionX",100,-12000,12000)
195 tgcGlobalY = ROOT.TH1F("tgcGlobalY","tgcGlobalY;TGC_Sim_hitGlobalPositionY",100,-12000,12000)
196 tgcGlobalZ = ROOT.TH1F("tgcGlobalZ","tgcGlobalZ;TGC_Sim_hitGlobalPositionZ",100,-20000,20000)
197 tgcGasGap = ROOT.TH1F("tgcGasGap","tgcGasGap;TGC_Sim_GasGap",4,0,4)
198 tgcChannel = ROOT.TH1F("tgcChannel","tgcChannel;TGC_Sim_channel",3,0,3)
199 tgcGlobalTime = ROOT.TH1F("tgcGlobalTime","tgcGlobalTime;TGC_sim_globalTime",100,0,500)
200 tgcKineticEnergy = ROOT.TH1F("tgcKineticEnergy","tgcKineticEnergy;TGC_Sim_kineticEnergy",100,0,400000)
201 tgcDepositEnergy = ROOT.TH1F("tgcDepositEnergy","tgcDepositEnergy;TGC_Sim_depositEnergy",100,0,0.002)
202 # TGC Digit Histograms
203 TGCDigitStationEta = ROOT.TH1F("TGCDigitStationEta","TGCDigitStationEta;Digits_TGC_stationEta",12,-6,6)
204 TGCDigitStationPhi = ROOT.TH1F("TGCDigitStationPhi","TGCDigitStationPhi;Digits_TGC_stationPhi",50,0,50)
205 TGCDigitGasGap = ROOT.TH1F("TGCDigitGasGap","TGCDigitGasGap;Digits_TGC_GasGap",4,0,4)
206 TGCDigitChannel = ROOT.TH1F("TGCDigitChannel","TGCDigitChannel;Digits_TGC_channel",120,0,120)
207 TGCDigitlocalX = ROOT.TH1F("TGCDigitlocalX","TGCDigitlocalX;Digits_TGC_localPosX",100,-1500,1500)
208 TGCDigitlocalY = ROOT.TH1F("TGCDigitlocalY","TGCDigitlocalY;Digits_TGC_localPosY",100,-1500,1500)
209 TGCDigitglobalX = ROOT.TH1F("TGCDigitglobalX","TGCDigitglobalX;Digits_TGC_globalPosX",100,-12000,12000)
210 TGCDigitglobalY = ROOT.TH1F("TGCDigitglobalY","TGCDigitglobalY;Digits_TGC_globalPosY",100,-12000,12000)
211 TGCDigitglobalZ = ROOT.TH1F("TGCDigitglobalZ","TGCDigitglobalZ;Digits_TGC_globalPosZ",100,-18000,18000)
212 # TGC SDO Histograms
213 TGCSDOStationEta = ROOT.TH1F("TGCSDOStationEta","TGCSDOStationEta;SDO_TGC_stationEta",12,-6,6)
214 TGCSDOStationPhi = ROOT.TH1F("TGCSDOStationPhi","TGCSDOStationPhi;SDO_TGC_stationPhi",50,0,50)
215 TGCSDOGasGap = ROOT.TH1F("TGCSDOGasGap","TGCSDOGasGap;SDO_TGC_GasGap",4,0,4)
216 TGCSDOChannel = ROOT.TH1F("TGCSDOChannel","TGCSDOChannel;SDO_TGC_channel",100,0,150)
217 TGCSDOWord = ROOT.TH1F("TGCSDOWord","TGCSDOWord;SDO_TGC_word",4,-2,2)
218 TGCSDOBarcode = ROOT.TH1F("TGCSDOBarcode","TGCSDOBarcode;SDO_TGC_barcode",100,0,200000)
219 TGCSDOGlobalX = ROOT.TH1F("TGCSDOGlobalX","TGCSDOGlobalX;SDO_TGC_globalPosX",100,-15000,15000)
220 TGCSDOGlobalY = ROOT.TH1F("TGCSDOGlobalY","TGCSDOGlobalY;SDO_TGC_globalPosY",100,-15000,15000)
221 TGCSDOGlobalZ = ROOT.TH1F("TGCSDOGlobalZ","TGCSDOGlobalZ;SDO_TGC_globalPosZ",100,-20000,20000)
222 TGCSDOLocalX = ROOT.TH1F("TGCSDOLocalX","TGCSDOLocalX;SDO_TGC_localPosX",100,0,100)
223 TGCSDOLocalY = ROOT.TH1F("TGCSDOLocalY","TGCSDOLocalY;SDO_TGC_localPosY",100,0,100)
224 TGCSDOGlotime = ROOT.TH1F("TGCSDOGlotime","TGCSDOGlotime;SDO_TGC_global_time",100,0,100)
225 # TGC RDO Histograms
226 TGCRDOStationEta = ROOT.TH1F("TGCRDOStationEta","TGCRDOStationEta;RDO_TGC_stationEta",12,-6,6)
227 TGCRDOStationPhi = ROOT.TH1F("TGCRDOStationPhi","TGCRDOStationPhi;RDO_TGC_stationPhi",50,0,50)
228 TGCRDOGasGap = ROOT.TH1F("TGCRDOGasGap","TGCRDOGasGap;RDO_TGC_GasGap",3,0,3)
229 TGCRDOChannel = ROOT.TH1F("TGCRDOChannel","TGCRDOChannel;RDO_TGC_channel",3,0,3)
230 TGCRDOGlobalX = ROOT.TH1F("TGCRDOGlobalX","TGCRDOGlobalX;RDO_TGC_globalPosX",100,-10000,10000)
231 TGCRDOGlobalY = ROOT.TH1F("TGCRDOGlobalY","TGCRDOGlobalY;RDO_TGC_globalPosY",100,-10000,10000)
232 TGCRDOGlobalZ = ROOT.TH1F("TGCRDOGlobalZ","TGCRDOGlobalZ;RDO_TGC_globalPosZ",100,-16000,16000)
233 TGCRDOLocalX = ROOT.TH1F("TGCRDOLocalX","TGCRDOLocalX;RDO_TGC_localPosX",100,-1200,-400)
234 # TGC PRD Histograms
235 TGCPRDStationEta = ROOT.TH1F("TGCPRDStationEta","TGCPRDStationEta;PRD_TGC_stationEta",12,-6,6)
236 TGCPRDStationPhi = ROOT.TH1F("TGCPRDStationPhi","TGCPRDStationPhi;PRD_TGC_stationPhi",50,0,50)
237 TGCPRDGasGap = ROOT.TH1F("TGCPRDGasGap","TGCPRDGasGap;PRD_TGC_GasGap",4,0,4)
238 TGCPRDChannel = ROOT.TH1F("TGCPRDChannel","TGCPRDChannel;PRD_TGC_channel",100,0,150)
239 TGCPRDisStrip = ROOT.TH1F("TGCPRDisStrip","TGCPRDisStrip;PRD_TGC_isStrip",2,0,2)
240 TGCPRDGlobalX = ROOT.TH1F("TGCPRDGlobalX","TGCPRDGlobalX;PRD_TGC_globalPosX",100,-15000,15000)
241 TGCPRDGlobalY = ROOT.TH1F("TGCPRDGlobalY","TGCPRDGlobalY;PRD_TGC_globalPosY",100,-15000,15000)
242 TGCPRDGlobalZ = ROOT.TH1F("TGCPRDGlobalZ","TGCPRDGlobalZ;PRD_TGC_globalPosZ",100,-18000,18000)
243 TGCPRDLocalX = ROOT.TH1F("TGCPRDLocalX","TGCPRDLocalX;PRD_TGC_localPosX",100,-1500,1500)
244
246 mmGlobalX = ROOT.TH1F("mmGlobalX","mmGlobalX;MM_Sim_hitGlobalPositionX",100,-5000,5000)
247 mmGlobalY = ROOT.TH1F("mmGlobalY","mmGlobalY;MM_Sim_hitGlobalPositionY",100,-5000,5000)
248 mmGlobalZ = ROOT.TH1F("mmGlobalZ","mmGlobalZ;MM_Sim_hitGlobalPositionZ",100,7000,8000)
249 mmStationEta = ROOT.TH1F("mmStationEta","mmStationEta;Hits_MM_off_stationEta",3,0,3)
250 mmStationPhi = ROOT.TH1F("mmStationPhi","mmStationPhi;Hits_MM_off_stationPhi",9,0,9)
251
253 stgcGlobalX = ROOT.TH1F("stgcGlobalX","stgcGlobalX;sTGC_Sim_hitGlobalPositionX",100,-5000,5000)
254 stgcGlobalY = ROOT.TH1F("stgcGlobalY","stgcGlobalY;sTGC_Sim_hitGlobalPositionY",100,-5000,5000)
255 stgcGlobalZ = ROOT.TH1F("stgcGlobalZ","stgcGlobalZ;sTGC_Sim_hitGlobalPositionZ",100,6000,8000)
256 stgcStationEta = ROOT.TH1F("stgcStationEta","stgcStationEta;Hits_sTGC_off_stationEta",4,0,4)
257 stgcStationPhi = ROOT.TH1F("stgcStationPhi","stgcStationPhi;Hits_sTGC_off_stationPhi",9,0,9)
258
259
260
261
262 def __init__(self, chamber_name="", eta_sel = None, sector_sel = None):
263 self.__chamber_name = chamber_name
264 self.__eta_sel = eta_sel
265 self.__sector_sel = sector_sel
266
267 # For negative eta x, ord() have returned values of (256 + x).
268 def Eta(eta):
269 if eta > 9:
270 return (eta - 256)
271 else:
272 return eta
273
274 def fill(self, TTree, n):
275
276 if self.__chamber_name == "TruthInfo":
277 MyHistoFiller.truthPtHist.Fill(TTree.MuEntry_ParticlePt[n]*0.001)
278 MyHistoFiller.truthEtaHist.Fill(TTree.MuEntry_ParticleEta[n])
279 MyHistoFiller.truthPhiHist.Fill(TTree.MuEntry_ParticlePhi[n])
280
281 if self.__chamber_name == "RPC_Sim":
282 if not (self.__eta_sel(TTree) and self.__sector_sel(TTree)):
283 return
284 else:
285 MyHistoFiller.rpcLocalX.Fill(TTree.RPC_SIM_localPostionX[n])
286 MyHistoFiller.rpcLocalY.Fill(TTree.RPC_SIM_localPostionY[n])
287 MyHistoFiller.rpcStationEta.Fill(MyHistoFiller.Eta(ord(TTree.RPC_SIM_stationEta[n])))
288 MyHistoFiller.rpcStationPhi.Fill(ord(TTree.RPC_SIM_stationPhi[n]))
289 MyHistoFiller.rpcGlobalX.Fill(TTree.RPC_SIM_GlobalPositionX[n])
290 MyHistoFiller.rpcGlobalY.Fill(TTree.RPC_SIM_GlobalPositionY[n])
291 MyHistoFiller.rpcGlobalZ.Fill(TTree.RPC_SIM_GlobalPositionZ[n])
292 MyHistoFiller.rpcGasGap.Fill(ord(TTree.RPC_SIM_gasGap[n]))
293
294 if self.__chamber_name == "RPC_Digit":
295 if not (self.__eta_sel(TTree) and self.__sector_sel(TTree)):
296 return
297 else:
298 MyHistoFiller.rpcDigitlocalX.Fill(TTree.Digits_RPC_localPosX[n])
299 MyHistoFiller.rpcDigitlocalY.Fill(TTree.Digits_RPC_localPosY[n])
300 MyHistoFiller.rpcDigitStationEta.Fill(MyHistoFiller.Eta(ord(TTree.Digits_RPC_stationEta[n])))
301 MyHistoFiller.rpcDigitStationPhi.Fill(ord(TTree.Digits_RPC_stationPhi[n]))
302 MyHistoFiller.rpcDigitglobalX.Fill(TTree.Digits_RPC_globalPosX[n])
303 MyHistoFiller.rpcDigitglobalY.Fill(TTree.Digits_RPC_globalPosY[n])
304 MyHistoFiller.rpcDigitglobalZ.Fill(TTree.Digits_RPC_globalPosZ[n])
305 MyHistoFiller.rpcDigitGasGap.Fill(ord(TTree.Digits_RPC_gasGap[n]))
306 MyHistoFiller.rpcDigitdoubletR.Fill(ord(TTree.Digits_RPC_doubletR[n]))
307 MyHistoFiller.rpcDigitdoubletZ.Fill(ord(TTree.Digits_RPC_doubletZ[n]))
308 MyHistoFiller.rpcDigitdoubletPhi.Fill(ord(TTree.Digits_RPC_doubletPhi[n]))
309 MyHistoFiller.rpcDigitMeasuresPhi.Fill(TTree.Digits_RPC_measuresPhi[n])
310 MyHistoFiller.rpcDigitTime.Fill(TTree.Digits_RPC_time[n])
311 MyHistoFiller.rpcDigitStrip.Fill(ord(TTree.Digits_RPC_strip[n]))
312 MyHistoFiller.rpcDigitStripNumber.Fill(TTree.Digits_RPC_stripNumber[n])
313 MyHistoFiller.rpcDigitglobalXY.Fill(TTree.Digits_RPC_globalPosX[n],TTree.Digits_RPC_globalPosY[n])
314 MyHistoFiller.rpcDigitglobalZY.Fill(TTree.Digits_RPC_globalPosZ[n],TTree.Digits_RPC_globalPosY[n])
315 MyHistoFiller.rpcDigitglobalZX.Fill(TTree.Digits_RPC_globalPosZ[n],TTree.Digits_RPC_globalPosX[n])
316
317 # Cut on BIS78 modules only -- remove stationEta condition to also plot BIS16 in the future
318 if TTree.Digits_RPC_stationIndex[n]=="BIS" and MyHistoFiller.Eta(ord(TTree.Digits_RPC_stationEta[n]))>6:
319 #use measuresPhi=0 to extract local eta/Y coordinate from localX
320 if TTree.Digits_RPC_measuresPhi[n]==0:
321 MyHistoFiller.bis_rpcDigitlocalY.Fill(TTree.Digits_RPC_localPosX[n])
322 #use measuresPhi=1 to extract local phi/X coordinate from localX
323 elif TTree.Digits_RPC_measuresPhi[n]==1:
324 MyHistoFiller.bis_rpcDigitlocalX.Fill(TTree.Digits_RPC_localPosX[n])
325
326 #use globalZ<6700 to get BIS7
327 if TTree.Digits_RPC_globalPosZ[n]<6700:
328 MyHistoFiller.bis7_rpcDigitTime.Fill(TTree.Digits_RPC_time[n])
329 if TTree.Digits_RPC_measuresPhi[n]==0:
330 MyHistoFiller.bis7_rpcDigitlocalY.Fill(TTree.Digits_RPC_localPosX[n])
331 elif TTree.Digits_RPC_measuresPhi[n]==1:
332 MyHistoFiller.bis7_rpcDigitlocalX.Fill(TTree.Digits_RPC_localPosX[n])
333 else:
334 #use globalZ>6700 to get BIS8
335 MyHistoFiller.bis8_rpcDigitTime.Fill(TTree.Digits_RPC_time[n])
336 if TTree.Digits_RPC_measuresPhi[n]==0:
337 MyHistoFiller.bis8_rpcDigitlocalY.Fill(TTree.Digits_RPC_localPosX[n])
338 elif TTree.Digits_RPC_measuresPhi[n]==1:
339 MyHistoFiller.bis8_rpcDigitlocalX.Fill(TTree.Digits_RPC_localPosX[n])
340
341 MyHistoFiller.bis_rpcDigitStationEta.Fill(MyHistoFiller.Eta(ord(TTree.Digits_RPC_stationEta[n])))
342 MyHistoFiller.bis_rpcDigitStationPhi.Fill(ord(TTree.Digits_RPC_stationPhi[n]))
343 MyHistoFiller.bis_rpcDigitglobalX.Fill(TTree.Digits_RPC_globalPosX[n])
344 MyHistoFiller.bis_rpcDigitglobalY.Fill(TTree.Digits_RPC_globalPosY[n])
345 MyHistoFiller.bis_rpcDigitglobalZ.Fill(TTree.Digits_RPC_globalPosZ[n])
346 MyHistoFiller.bis_rpcDigitGasGap.Fill(ord(TTree.Digits_RPC_gasGap[n]))
347 MyHistoFiller.bis_rpcDigitdoubletR.Fill(ord(TTree.Digits_RPC_doubletR[n]))
348 MyHistoFiller.bis_rpcDigitdoubletZ.Fill(ord(TTree.Digits_RPC_doubletZ[n]))
349 MyHistoFiller.bis_rpcDigitdoubletPhi.Fill(ord(TTree.Digits_RPC_doubletPhi[n]))
350 MyHistoFiller.bis_rpcDigitMeasuresPhi.Fill(TTree.Digits_RPC_measuresPhi[n])
351 MyHistoFiller.bis_rpcDigitTime.Fill(TTree.Digits_RPC_time[n])
352 MyHistoFiller.bis_rpcDigitStrip.Fill(ord(TTree.Digits_RPC_strip[n]))
353 MyHistoFiller.bis_rpcDigitStripNumber.Fill(TTree.Digits_RPC_stripNumber[n])
354 MyHistoFiller.bis_rpcDigitglobalXY.Fill(TTree.Digits_RPC_globalPosX[n],TTree.Digits_RPC_globalPosY[n])
355 MyHistoFiller.bis_rpcDigitglobalZY.Fill(TTree.Digits_RPC_globalPosZ[n],TTree.Digits_RPC_globalPosY[n])
356 MyHistoFiller.bis_rpcDigitglobalZX.Fill(TTree.Digits_RPC_globalPosZ[n],TTree.Digits_RPC_globalPosX[n])
357
358
359 if self.__chamber_name == "MDT_Sim":
360 if not (self.__eta_sel(TTree) and self.__sector_sel(TTree)):
361 return
362 else:
363 MyHistoFiller.mdtLocalX.Fill(TTree.MDT_Sim_hitLocalPositionX[n])
364 MyHistoFiller.mdtLocalY.Fill(TTree.MDT_Sim_hitLocalPositionY[n])
365 MyHistoFiller.mdtLocalZ.Fill(TTree.MDT_Sim_hitLocalPositionZ[n])
366 MyHistoFiller.mdtStationEta.Fill(MyHistoFiller.Eta(ord(TTree.MDT_Sim_stationEta[n])))
367 MyHistoFiller.mdtStationPhi.Fill(ord(TTree.MDT_Sim_stationPhi[n]))
368 MyHistoFiller.mdtGlobalX.Fill(TTree.MDT_Sim_hitGlobalPositionX[n])
369 MyHistoFiller.mdtGlobalY.Fill(TTree.MDT_Sim_hitGlobalPositionY[n])
370 MyHistoFiller.mdtGlobalZ.Fill(TTree.MDT_Sim_hitGlobalPositionZ[n])
371 MyHistoFiller.mdtTube.Fill(ord(TTree.MDT_Sim_tube[n]))
372
373 if self.__chamber_name == "MDT_Digit":
374 if not (self.__eta_sel(TTree) and self.__sector_sel(TTree)):
375 return
376 else:
377 MyHistoFiller.mdtDigitlocalTubePosX.Fill(TTree.Digits_MDT_localTubePosX[n])
378 MyHistoFiller.mdtDigitlocalTubePosY.Fill(TTree.Digits_MDT_localTubePosY[n])
379 MyHistoFiller.mdtDigitlocalTubePosZ.Fill(TTree.Digits_MDT_localTubePosZ[n])
380 MyHistoFiller.mdtDigitStationEta.Fill(MyHistoFiller.Eta(ord(TTree.Digits_MDT_stationEta[n])))
381 MyHistoFiller.mdtDigitStationPhi.Fill(ord(TTree.Digits_MDT_stationPhi[n]))
382 MyHistoFiller.mdtDigitglobalX.Fill(TTree.Digits_MDT_globalPosX[n])
383 MyHistoFiller.mdtDigitglobalY.Fill(TTree.Digits_MDT_globalPosY[n])
384 MyHistoFiller.mdtDigitglobalZ.Fill(TTree.Digits_MDT_globalPosZ[n])
385 MyHistoFiller.mdtDigitNumberOfMultilayers.Fill(TTree.Digits_MDT_numberOfMultilayers[n])
386 MyHistoFiller.mdtDigitTube.Fill(ord(TTree.Digits_MDT_tube[n]))
387 MyHistoFiller.mdtDigitTubeLayer.Fill(ord(TTree.Digits_MDT_tubeLayer[n]))
388 MyHistoFiller.mdtDigitMultilayer.Fill(ord(TTree.Digits_MDT_multiLayer[n]))
389 MyHistoFiller.mdtDigitTime.Fill(TTree.Digits_MDT_time[n])
390 MyHistoFiller.mdtDigitCharge.Fill(TTree.Digits_MDT_charge[n])
391 MyHistoFiller.mdtDigitglobalXY.Fill(TTree.Digits_MDT_globalPosX[n],TTree.Digits_MDT_globalPosY[n])
392 MyHistoFiller.mdtDigitglobalZY.Fill(TTree.Digits_MDT_globalPosZ[n],TTree.Digits_MDT_globalPosY[n])
393 MyHistoFiller.mdtDigitglobalZX.Fill(TTree.Digits_MDT_globalPosZ[n],TTree.Digits_MDT_globalPosX[n])
394 # Cut on BIS78 modules only
395 if TTree.Digits_MDT_stationIndex[n]=="BIS" and MyHistoFiller.Eta(ord(TTree.Digits_MDT_stationEta[n]))>6:
396 MyHistoFiller.bis_mdtDigitlocalTubePosX.Fill(TTree.Digits_MDT_localTubePosX[n])
397 MyHistoFiller.bis_mdtDigitlocalTubePosY.Fill(TTree.Digits_MDT_localTubePosY[n])
398 MyHistoFiller.bis_mdtDigitlocalTubePosZ.Fill(TTree.Digits_MDT_localTubePosZ[n])
399 MyHistoFiller.bis_mdtDigitStationEta.Fill(MyHistoFiller.Eta(ord(TTree.Digits_MDT_stationEta[n])))
400 MyHistoFiller.bis_mdtDigitStationPhi.Fill(ord(TTree.Digits_MDT_stationPhi[n]))
401 MyHistoFiller.bis_mdtDigitglobalX.Fill(TTree.Digits_MDT_globalPosX[n])
402 MyHistoFiller.bis_mdtDigitglobalY.Fill(TTree.Digits_MDT_globalPosY[n])
403 MyHistoFiller.bis_mdtDigitglobalZ.Fill(TTree.Digits_MDT_globalPosZ[n])
404 MyHistoFiller.bis_mdtDigitNumberOfMultilayers.Fill(TTree.Digits_MDT_numberOfMultilayers[n])
405 MyHistoFiller.bis_mdtDigitTube.Fill(ord(TTree.Digits_MDT_tube[n]))
406 MyHistoFiller.bis_mdtDigitTubeLayer.Fill(ord(TTree.Digits_MDT_tubeLayer[n]))
407 MyHistoFiller.bis_mdtDigitMultilayer.Fill(ord(TTree.Digits_MDT_multiLayer[n]))
408 MyHistoFiller.bis_mdtDigitTime.Fill(TTree.Digits_MDT_time[n])
409 MyHistoFiller.bis_mdtDigitCharge.Fill(TTree.Digits_MDT_charge[n])
410 MyHistoFiller.bis_mdtDigitglobalXY.Fill(TTree.Digits_MDT_globalPosX[n],TTree.Digits_MDT_globalPosY[n])
411 MyHistoFiller.bis_mdtDigitglobalZY.Fill(TTree.Digits_MDT_globalPosZ[n],TTree.Digits_MDT_globalPosY[n])
412 MyHistoFiller.bis_mdtDigitglobalZX.Fill(TTree.Digits_MDT_globalPosZ[n],TTree.Digits_MDT_globalPosX[n])
413 MyHistoFiller.bis_mdtDigitlocalTubePosXY.Fill(TTree.Digits_MDT_localTubePosX[n],TTree.Digits_MDT_localTubePosY[n])
414 MyHistoFiller.bis_mdtDigitlocalTubePosZY.Fill(TTree.Digits_MDT_localTubePosZ[n],TTree.Digits_MDT_localTubePosY[n])
415 MyHistoFiller.bis_mdtDigitlocalTubePosZX.Fill(TTree.Digits_MDT_localTubePosZ[n],TTree.Digits_MDT_localTubePosX[n])
416
417
418 if self.__chamber_name == "CSC_Sim":
419 if not (self.__eta_sel(TTree) and self.__sector_sel(TTree)):
420 return
421 else:
422 MyHistoFiller.cscGlobalX.Fill(TTree.CSC_Sim_hitGlobalPositionX[n])
423 MyHistoFiller.cscGlobalY.Fill(TTree.CSC_Sim_hitGlobalPositionY[n])
424 MyHistoFiller.cscGlobalZ.Fill(TTree.CSC_Sim_hitGlobalPositionZ[n])
425 MyHistoFiller.cscWireLayer.Fill(ord(TTree.CSC_Sim_wireLayer[n]))
426 MyHistoFiller.cscStrip.Fill(ord(TTree.CSC_Sim_strip[n]))
427 MyHistoFiller.cscglobalTime.Fill(TTree.CSC_Sim_globalTime[n])
428 MyHistoFiller.cscKineticEnergy.Fill(TTree.CSC_Sim_kineticEnergy[n])
429 MyHistoFiller.cscDepositEnergy.Fill(TTree.CSC_Sim_depositEnergy[n])
430 MyHistoFiller.cscSimStationEta.Fill(MyHistoFiller.Eta(ord(TTree.CSC_Sim_stationEta[n])))
431 MyHistoFiller.cscSimStationPhi.Fill(ord(TTree.CSC_Sim_stationPhi[n]))
432 MyHistoFiller.cscSimChamberLayer.Fill(ord(TTree.CSC_Sim_chamberLayer[n]))
433
434 if self.__chamber_name == "CSC_Digit":
435 if not (self.__eta_sel(TTree) and self.__sector_sel(TTree)):
436 return
437 else:
438 MyHistoFiller.CSCDigitStationEta.Fill(MyHistoFiller.Eta(ord(TTree.Digits_CSC_stationEta[n])))
439 MyHistoFiller.CSCDigitStationPhi.Fill(ord(TTree.Digits_CSC_stationPhi[n]))
440 MyHistoFiller.CSCDigitWireLayer.Fill(ord(TTree.Digits_CSC_wireLayer[n]))
441 MyHistoFiller.CSCDigitlocalPosX.Fill(TTree.Digits_CSC_localPosX[n])
442 MyHistoFiller.CSCDigitglobalPosX.Fill(TTree.Digits_CSC_globalPosX[n])
443 MyHistoFiller.CSCDigitglobalPosY.Fill(TTree.Digits_CSC_globalPosY[n])
444 MyHistoFiller.CSCDigitglobalPosZ.Fill(TTree.Digits_CSC_globalPosZ[n])
445
446 if self.__chamber_name == "CSC_SDO":
447 if not (self.__eta_sel(TTree) and self.__sector_sel(TTree)):
448 return
449 else:
450 MyHistoFiller.CSCSDOStationEta.Fill(MyHistoFiller.Eta(ord(TTree.SDO_CSC_stationEta[n])))
451 MyHistoFiller.CSCSDOStationPhi.Fill(ord(TTree.SDO_CSC_stationPhi[n]))
452 MyHistoFiller.CSCSDOChamberLayer.Fill(ord(TTree.SDO_CSC_chamberLayer[n]))
453 MyHistoFiller.CSCSDOWireLayer.Fill(ord(TTree.SDO_CSC_wireLayer[n]))
454 MyHistoFiller.CSCSDOStrip.Fill(ord(TTree.SDO_CSC_strip[n]))
455 MyHistoFiller.CSCSDOBarcode.Fill(TTree.SDO_CSC_barcode[n])
456 MyHistoFiller.CSCSDOGlobalX.Fill(TTree.SDO_CSC_globalPosX[n])
457 MyHistoFiller.CSCSDOGlobalY.Fill(TTree.SDO_CSC_globalPosY[n])
458 MyHistoFiller.CSCSDOGlobalZ.Fill(TTree.SDO_CSC_globalPosZ[n])
459 MyHistoFiller.CSCSDOLocalX.Fill(TTree.SDO_CSC_localPosX[n])
460 MyHistoFiller.CSCSDOLocalY.Fill(TTree.SDO_CSC_localPosY[n])
461
462 if self.__chamber_name == "CSC_RDO":
463 if not (self.__eta_sel(TTree) and self.__sector_sel(TTree)):
464 return
465 else:
466 MyHistoFiller.CSCRDOStationEta.Fill(MyHistoFiller.Eta(ord(TTree.RDO_CSC_stationEta[n])))
467 MyHistoFiller.CSCRDOStationPhi.Fill(ord(TTree.RDO_CSC_stationPhi[n]))
468 MyHistoFiller.CSCRDOChamberLayer.Fill(ord(TTree.RDO_CSC_chamberLayer[n]))
469 MyHistoFiller.CSCRDOWireLayer.Fill(ord(TTree.RDO_CSC_wireLayer[n]))
470 MyHistoFiller.CSCRDOStrip.Fill(ord(TTree.RDO_CSC_strip[n]))
471 MyHistoFiller.CSCRDOmeasurePhi.Fill(TTree.RDO_CSC_measuresPhi[n])
472 MyHistoFiller.CSCRDOTime.Fill(TTree.RDO_CSC_time[n])
473 MyHistoFiller.CSCRDOGlobalX.Fill(TTree.RDO_CSC_globalPosX[n])
474 MyHistoFiller.CSCRDOGlobalY.Fill(TTree.RDO_CSC_globalPosY[n])
475 MyHistoFiller.CSCRDOGlobalZ.Fill(TTree.RDO_CSC_globalPosZ[n])
476 MyHistoFiller.CSCRDOLocalX.Fill(TTree.RDO_CSC_localPosX[n])
477 MyHistoFiller.CSCRDOLocalY.Fill(TTree.RDO_CSC_localPosY[n])
478
479 if self.__chamber_name == "CSC_PRD":
480 if not (self.__eta_sel(TTree) and self.__sector_sel(TTree)):
481 return
482 else:
483 MyHistoFiller.CSCPRDStationEta.Fill(MyHistoFiller.Eta(ord(TTree.PRD_CSC_stationEta[n])))
484 MyHistoFiller.CSCPRDStationPhi.Fill(ord(TTree.PRD_CSC_stationPhi[n]))
485 MyHistoFiller.CSCPRDChamberLayer.Fill(ord(TTree.PRD_CSC_chamberLayer[n]))
486 MyHistoFiller.CSCPRDWireLayer.Fill(ord(TTree.PRD_CSC_wireLayer[n]))
487 MyHistoFiller.CSCPRDStrip.Fill(ord(TTree.PRD_CSC_strip[n]))
488 MyHistoFiller.CSCPRDmeasurePhi.Fill(TTree.PRD_CSC_measuresPhi[n])
489 MyHistoFiller.CSCPRDTime.Fill(TTree.PRD_CSC_time[n])
490 MyHistoFiller.CSCPRDCharge.Fill(TTree.PRD_CSC_charge[n])
491 MyHistoFiller.CSCPRDGlobalX.Fill(TTree.PRD_CSC_globalPosX[n])
492 MyHistoFiller.CSCPRDGlobalY.Fill(TTree.PRD_CSC_globalPosY[n])
493 MyHistoFiller.CSCPRDGlobalZ.Fill(TTree.PRD_CSC_globalPosZ[n])
494 MyHistoFiller.CSCPRDLocalX.Fill(TTree.PRD_CSC_localPosX[n])
495
496 if self.__chamber_name == "TGC_Sim":
497 if not (self.__eta_sel(TTree) and self.__sector_sel(TTree)):
498 return
499 else:
500 MyHistoFiller.tgcLocalX.Fill(TTree.TGC_Sim_hitLocalPositionX[n])
501 MyHistoFiller.tgcLocalY.Fill(TTree.TGC_Sim_hitLocalPositionY[n])
502 MyHistoFiller.tgcLocalZ.Fill(TTree.TGC_Sim_hitLocalPositionZ[n])
503 MyHistoFiller.tgcStationEta.Fill(MyHistoFiller.Eta(ord(TTree.TGC_Sim_stationEta[n])))
504 MyHistoFiller.tgcStationPhi.Fill(ord(TTree.TGC_Sim_stationPhi[n]))
505 MyHistoFiller.tgcGlobalX.Fill(TTree.TGC_Sim_hitGlobalPositionX[n])
506 MyHistoFiller.tgcGlobalY.Fill(TTree.TGC_Sim_hitGlobalPositionY[n])
507 MyHistoFiller.tgcGlobalZ.Fill(TTree.TGC_Sim_hitGlobalPositionZ[n])
508 MyHistoFiller.tgcGasGap.Fill(ord(TTree.TGC_Sim_gasGap[n]))
509 MyHistoFiller.tgcChannel.Fill(ord(TTree.TGC_Sim_channel[n]))
510 MyHistoFiller.tgcGlobalTime.Fill(TTree.TGC_Sim_globalTime[n])
511 MyHistoFiller.tgcKineticEnergy.Fill(TTree.TGC_Sim_kineticEnergy[n])
512 MyHistoFiller.tgcDepositEnergy.Fill(TTree.TGC_Sim_depositEnergy[n])
513
514 if self.__chamber_name == "TGC_Digit":
515 if not (self.__eta_sel(TTree) and self.__sector_sel(TTree)):
516 return
517 else:
518 MyHistoFiller.TGCDigitStationEta.Fill(MyHistoFiller.Eta(ord(TTree.Digits_TGC_stationEta[n])))
519 MyHistoFiller.TGCDigitStationPhi.Fill(ord(TTree.Digits_TGC_stationPhi[n]))
520 MyHistoFiller.TGCDigitGasGap.Fill(ord(TTree.Digits_TGC_gasGap[n]))
521 MyHistoFiller.TGCDigitChannel.Fill(ord(TTree.Digits_TGC_channel[n]))
522 MyHistoFiller.TGCDigitlocalX.Fill(TTree.Digits_TGC_localPosX[n])
523 MyHistoFiller.TGCDigitlocalY.Fill(TTree.Digits_TGC_localPosY[n])
524 MyHistoFiller.TGCDigitglobalX.Fill(TTree.Digits_TGC_globalPosX[n])
525 MyHistoFiller.TGCDigitglobalY.Fill(TTree.Digits_TGC_globalPosY[n])
526 MyHistoFiller.TGCDigitglobalZ.Fill(TTree.Digits_TGC_globalPosZ[n])
527
528 if self.__chamber_name == "TGC_SDO":
529 if not (self.__eta_sel(TTree) and self.__sector_sel(TTree)):
530 return
531 else:
532 MyHistoFiller.TGCSDOStationEta.Fill(MyHistoFiller.Eta(ord(TTree.SDO_TGC_stationEta[n])))
533 MyHistoFiller.TGCSDOStationPhi.Fill(ord(TTree.SDO_TGC_stationPhi[n]))
534 MyHistoFiller.TGCSDOGasGap.Fill(ord(TTree.SDO_TGC_gasGap[n]))
535 MyHistoFiller.TGCSDOChannel.Fill(ord(TTree.SDO_TGC_channel[n]))
536 MyHistoFiller.TGCSDOWord.Fill(TTree.SDO_TGC_word[n])
537 MyHistoFiller.TGCSDOBarcode.Fill(TTree.SDO_TGC_barcode[n])
538 MyHistoFiller.TGCSDOGlobalX.Fill(TTree.SDO_TGC_globalPosX[n])
539 MyHistoFiller.TGCSDOGlobalY.Fill(TTree.SDO_TGC_globalPosY[n])
540 MyHistoFiller.TGCSDOGlobalZ.Fill(TTree.SDO_TGC_globalPosZ[n])
541 MyHistoFiller.TGCSDOLocalX.Fill(TTree.SDO_TGC_localPosX[n])
542 MyHistoFiller.TGCSDOLocalY.Fill(TTree.SDO_TGC_localPosY[n])
543 MyHistoFiller.TGCSDOGlotime.Fill(TTree.SDO_TGC_global_time[n])
544
545 if self.__chamber_name == "TGC_RDO":
546 if not (self.__eta_sel(TTree) and self.__sector_sel(TTree)):
547 return
548 else:
549 MyHistoFiller.TGCRDOStationEta.Fill(MyHistoFiller.Eta(ord(TTree.RDO_TGC_stationEta[n])))
550 MyHistoFiller.TGCRDOStationPhi.Fill(ord(TTree.RDO_TGC_stationPhi[n]))
551 MyHistoFiller.TGCRDOGasGap.Fill(ord(TTree.RDO_TGC_gasGap[n]))
552 MyHistoFiller.TGCRDOChannel.Fill(ord(TTree.RDO_TGC_channel[n]))
553 MyHistoFiller.TGCRDOGlobalX.Fill(TTree.RDO_TGC_globalPosX[n])
554 MyHistoFiller.TGCRDOGlobalY.Fill(TTree.RDO_TGC_globalPosY[n])
555 MyHistoFiller.TGCRDOGlobalZ.Fill(TTree.RDO_TGC_globalPosZ[n])
556 MyHistoFiller.TGCRDOLocalX.Fill(TTree.RDO_TGC_localPosX[n])
557
558 if self.__chamber_name == "TGC_PRD":
559 if not (self.__eta_sel(TTree) and self.__sector_sel(TTree)):
560 return
561 else:
562 MyHistoFiller.TGCPRDStationEta.Fill(MyHistoFiller.Eta(ord(TTree.PRD_TGC_stationEta[n])))
563 MyHistoFiller.TGCPRDStationPhi.Fill(ord(TTree.PRD_TGC_stationPhi[n]))
564 MyHistoFiller.TGCPRDGasGap.Fill(ord(TTree.PRD_TGC_gasGap[n]))
565 MyHistoFiller.TGCPRDChannel.Fill(ord(TTree.PRD_TGC_channel[n]))
566 MyHistoFiller.TGCPRDisStrip.Fill(TTree.PRD_TGC_isStrip[n])
567 MyHistoFiller.TGCPRDGlobalX.Fill(TTree.PRD_TGC_globalPosX[n])
568 MyHistoFiller.TGCPRDGlobalY.Fill(TTree.PRD_TGC_globalPosY[n])
569 MyHistoFiller.TGCPRDGlobalZ.Fill(TTree.PRD_TGC_globalPosZ[n])
570 MyHistoFiller.TGCPRDLocalX.Fill(TTree.PRD_TGC_localPosX[n])
571
572 if self.__chamber_name == "MM_Sim":
573 if not (self.__eta_sel(TTree) and self.__sector_sel(TTree)):
574 return
575 else:
576 MyHistoFiller.mmGlobalX.Fill(TTree.Hits_MM_GlobalPositionX[n])
577 MyHistoFiller.mmGlobalY.Fill(TTree.Hits_MM_GlobalPositionY[n])
578 MyHistoFiller.mmGlobalZ.Fill(TTree.Hits_MM_GlobalPositionZ[n])
579 MyHistoFiller.mmStationEta.Fill(MyHistoFiller.Eta(ord(TTree.Hits_MM_off_stationEta[n])))
580 MyHistoFiller.mmStationPhi.Fill(ord(TTree.Hits_MM_off_stationPhi[n]))
581
582 if self.__chamber_name == "sTGC_Sim":
583 if not (self.__eta_sel(TTree) and self.__sector_sel(TTree)):
584 return
585 else:
586 MyHistoFiller.stgcGlobalX.Fill(TTree.Hits_sTGC_GlobalPositionX[n])
587 MyHistoFiller.stgcGlobalY.Fill(TTree.Hits_sTGC_GlobalPositionY[n])
588 MyHistoFiller.stgcGlobalZ.Fill(TTree.Hits_sTGC_GlobalPositionZ[n])
589 MyHistoFiller.stgcStationEta.Fill(MyHistoFiller.Eta(ord(TTree.Hits_sTGC_off_stationEta[n])))
590 MyHistoFiller.stgcStationPhi.Fill(ord(TTree.Hits_sTGC_off_stationPhi[n]))
591
592
593 def write(self, outdir):
594
595 if self.__chamber_name == "TruthInfo":
596 outdir.WriteTObject(MyHistoFiller.truthPtHist, MyHistoFiller.truthPtHist.GetName())
597 outdir.WriteTObject(MyHistoFiller.truthEtaHist, MyHistoFiller.truthEtaHist.GetName())
598 outdir.WriteTObject(MyHistoFiller.truthPhiHist, MyHistoFiller.truthPhiHist.GetName())
599
600 if self.__chamber_name == "RPC_Sim":
601 outdir.WriteTObject(MyHistoFiller.rpcLocalX, MyHistoFiller.rpcLocalX.GetName())
602 outdir.WriteTObject(MyHistoFiller.rpcLocalY, MyHistoFiller.rpcLocalY.GetName())
603 outdir.WriteTObject(MyHistoFiller.rpcStationEta, MyHistoFiller.rpcStationEta.GetName())
604 outdir.WriteTObject(MyHistoFiller.rpcStationPhi, MyHistoFiller.rpcStationPhi.GetName())
605 outdir.WriteTObject(MyHistoFiller.rpcGlobalX, MyHistoFiller.rpcGlobalX.GetName())
606 outdir.WriteTObject(MyHistoFiller.rpcGlobalY, MyHistoFiller.rpcGlobalY.GetName())
607 outdir.WriteTObject(MyHistoFiller.rpcGlobalZ, MyHistoFiller.rpcGlobalZ.GetName())
608 outdir.WriteTObject(MyHistoFiller.rpcGasGap, MyHistoFiller.rpcGasGap.GetName())
609
610 if self.__chamber_name == "RPC_Digit":
611 outdir.WriteTObject(MyHistoFiller.rpcDigitlocalX, MyHistoFiller.rpcDigitlocalX.GetName())
612 outdir.WriteTObject(MyHistoFiller.rpcDigitlocalY, MyHistoFiller.rpcDigitlocalY.GetName())
613 outdir.WriteTObject(MyHistoFiller.rpcDigitStationEta, MyHistoFiller.rpcDigitStationEta.GetName())
614 outdir.WriteTObject(MyHistoFiller.rpcDigitStationPhi, MyHistoFiller.rpcDigitStationPhi.GetName())
615 outdir.WriteTObject(MyHistoFiller.rpcDigitglobalX, MyHistoFiller.rpcDigitglobalX.GetName())
616 outdir.WriteTObject(MyHistoFiller.rpcDigitglobalY, MyHistoFiller.rpcDigitglobalY.GetName())
617 outdir.WriteTObject(MyHistoFiller.rpcDigitglobalZ, MyHistoFiller.rpcDigitglobalZ.GetName())
618 outdir.WriteTObject(MyHistoFiller.rpcDigitGasGap, MyHistoFiller.rpcDigitGasGap.GetName())
619 outdir.WriteTObject(MyHistoFiller.rpcDigitStrip, MyHistoFiller.rpcDigitStrip.GetName())
620 outdir.WriteTObject(MyHistoFiller.rpcDigitStripNumber, MyHistoFiller.rpcDigitStripNumber.GetName())
621 outdir.WriteTObject(MyHistoFiller.rpcDigitMeasuresPhi, MyHistoFiller.rpcDigitMeasuresPhi.GetName())
622 outdir.WriteTObject(MyHistoFiller.rpcDigitTime, MyHistoFiller.rpcDigitTime.GetName())
623 outdir.WriteTObject(MyHistoFiller.rpcDigitdoubletR, MyHistoFiller.rpcDigitdoubletR.GetName())
624 outdir.WriteTObject(MyHistoFiller.rpcDigitdoubletZ, MyHistoFiller.rpcDigitdoubletZ.GetName())
625 outdir.WriteTObject(MyHistoFiller.rpcDigitdoubletPhi, MyHistoFiller.rpcDigitdoubletPhi.GetName())
626 outdir.WriteTObject(MyHistoFiller.rpcDigitglobalXY, MyHistoFiller.rpcDigitglobalXY.GetName())
627 outdir.WriteTObject(MyHistoFiller.rpcDigitglobalZY, MyHistoFiller.rpcDigitglobalZY.GetName())
628 outdir.WriteTObject(MyHistoFiller.rpcDigitglobalZX, MyHistoFiller.rpcDigitglobalZX.GetName())
629 # BIS rpcDigit
630 outdir.WriteTObject(MyHistoFiller.bis_rpcDigitlocalX, MyHistoFiller.bis_rpcDigitlocalX.GetName())
631 outdir.WriteTObject(MyHistoFiller.bis_rpcDigitlocalY, MyHistoFiller.bis_rpcDigitlocalY.GetName())
632 outdir.WriteTObject(MyHistoFiller.bis_rpcDigitStationEta, MyHistoFiller.bis_rpcDigitStationEta.GetName())
633 outdir.WriteTObject(MyHistoFiller.bis_rpcDigitStationPhi, MyHistoFiller.bis_rpcDigitStationPhi.GetName())
634 outdir.WriteTObject(MyHistoFiller.bis_rpcDigitglobalX, MyHistoFiller.bis_rpcDigitglobalX.GetName())
635 outdir.WriteTObject(MyHistoFiller.bis_rpcDigitglobalY, MyHistoFiller.bis_rpcDigitglobalY.GetName())
636 outdir.WriteTObject(MyHistoFiller.bis_rpcDigitglobalZ, MyHistoFiller.bis_rpcDigitglobalZ.GetName())
637 outdir.WriteTObject(MyHistoFiller.bis_rpcDigitGasGap, MyHistoFiller.bis_rpcDigitGasGap.GetName())
638 outdir.WriteTObject(MyHistoFiller.bis_rpcDigitStrip, MyHistoFiller.bis_rpcDigitStrip.GetName())
639 outdir.WriteTObject(MyHistoFiller.bis_rpcDigitStripNumber, MyHistoFiller.bis_rpcDigitStripNumber.GetName())
640 outdir.WriteTObject(MyHistoFiller.bis_rpcDigitMeasuresPhi, MyHistoFiller.bis_rpcDigitMeasuresPhi.GetName())
641 outdir.WriteTObject(MyHistoFiller.bis_rpcDigitTime, MyHistoFiller.bis_rpcDigitTime.GetName())
642 outdir.WriteTObject(MyHistoFiller.bis7_rpcDigitTime, MyHistoFiller.bis7_rpcDigitTime.GetName())
643 outdir.WriteTObject(MyHistoFiller.bis8_rpcDigitTime, MyHistoFiller.bis8_rpcDigitTime.GetName())
644 outdir.WriteTObject(MyHistoFiller.bis_rpcDigitdoubletR, MyHistoFiller.bis_rpcDigitdoubletR.GetName())
645 outdir.WriteTObject(MyHistoFiller.bis_rpcDigitdoubletZ, MyHistoFiller.bis_rpcDigitdoubletZ.GetName())
646 outdir.WriteTObject(MyHistoFiller.bis_rpcDigitdoubletPhi, MyHistoFiller.bis_rpcDigitdoubletPhi.GetName())
647 outdir.WriteTObject(MyHistoFiller.bis_rpcDigitglobalXY, MyHistoFiller.bis_rpcDigitglobalXY.GetName())
648 outdir.WriteTObject(MyHistoFiller.bis_rpcDigitglobalZY, MyHistoFiller.bis_rpcDigitglobalZY.GetName())
649 outdir.WriteTObject(MyHistoFiller.bis_rpcDigitglobalZX, MyHistoFiller.bis_rpcDigitglobalZX.GetName())
650
651 if self.__chamber_name == "MDT_Sim":
652 outdir.WriteTObject(MyHistoFiller.mdtLocalX, MyHistoFiller.mdtLocalX.GetName())
653 outdir.WriteTObject(MyHistoFiller.mdtLocalY, MyHistoFiller.mdtLocalY.GetName())
654 outdir.WriteTObject(MyHistoFiller.mdtLocalZ, MyHistoFiller.mdtLocalZ.GetName())
655 outdir.WriteTObject(MyHistoFiller.mdtStationEta, MyHistoFiller.mdtStationEta.GetName())
656 outdir.WriteTObject(MyHistoFiller.mdtStationPhi, MyHistoFiller.mdtStationPhi.GetName())
657 outdir.WriteTObject(MyHistoFiller.mdtGlobalX, MyHistoFiller.mdtGlobalX.GetName())
658 outdir.WriteTObject(MyHistoFiller.mdtGlobalY, MyHistoFiller.mdtGlobalY.GetName())
659 outdir.WriteTObject(MyHistoFiller.mdtGlobalZ, MyHistoFiller.mdtGlobalZ.GetName())
660 outdir.WriteTObject(MyHistoFiller.mdtTube, MyHistoFiller.mdtTube.GetName())
661
662 if self.__chamber_name == "MDT_Digit":
663 outdir.WriteTObject(MyHistoFiller.mdtDigitlocalTubePosX, MyHistoFiller.mdtDigitlocalTubePosX.GetName())
664 outdir.WriteTObject(MyHistoFiller.mdtDigitlocalTubePosY, MyHistoFiller.mdtDigitlocalTubePosY.GetName())
665 outdir.WriteTObject(MyHistoFiller.mdtDigitlocalTubePosZ, MyHistoFiller.mdtDigitlocalTubePosZ.GetName())
666 outdir.WriteTObject(MyHistoFiller.mdtDigitStationEta, MyHistoFiller.mdtDigitStationEta.GetName())
667 outdir.WriteTObject(MyHistoFiller.mdtDigitStationPhi, MyHistoFiller.mdtDigitStationPhi.GetName())
668 outdir.WriteTObject(MyHistoFiller.mdtDigitglobalX, MyHistoFiller.mdtDigitglobalX.GetName())
669 outdir.WriteTObject(MyHistoFiller.mdtDigitglobalY, MyHistoFiller.mdtDigitglobalY.GetName())
670 outdir.WriteTObject(MyHistoFiller.mdtDigitglobalZ, MyHistoFiller.mdtDigitglobalZ.GetName())
671 outdir.WriteTObject(MyHistoFiller.mdtDigitNumberOfMultilayers, MyHistoFiller.mdtDigitNumberOfMultilayers.GetName())
672 outdir.WriteTObject(MyHistoFiller.mdtDigitTube, MyHistoFiller.mdtDigitTube.GetName())
673 outdir.WriteTObject(MyHistoFiller.mdtDigitTubeLayer, MyHistoFiller.mdtDigitTubeLayer.GetName())
674 outdir.WriteTObject(MyHistoFiller.mdtDigitMultilayer, MyHistoFiller.mdtDigitMultilayer.GetName())
675 outdir.WriteTObject(MyHistoFiller.mdtDigitTime, MyHistoFiller.mdtDigitTime.GetName())
676 outdir.WriteTObject(MyHistoFiller.mdtDigitCharge, MyHistoFiller.mdtDigitCharge.GetName())
677 outdir.WriteTObject(MyHistoFiller.mdtDigitglobalXY, MyHistoFiller.mdtDigitglobalXY.GetName())
678 outdir.WriteTObject(MyHistoFiller.mdtDigitglobalZY, MyHistoFiller.mdtDigitglobalZY.GetName())
679 outdir.WriteTObject(MyHistoFiller.mdtDigitglobalZX, MyHistoFiller.mdtDigitglobalZX.GetName())
680 # BIS mdtDigit
681 outdir.WriteTObject(MyHistoFiller.bis_mdtDigitlocalTubePosX, MyHistoFiller.bis_mdtDigitlocalTubePosX.GetName())
682 outdir.WriteTObject(MyHistoFiller.bis_mdtDigitlocalTubePosY, MyHistoFiller.bis_mdtDigitlocalTubePosY.GetName())
683 outdir.WriteTObject(MyHistoFiller.bis_mdtDigitlocalTubePosZ, MyHistoFiller.bis_mdtDigitlocalTubePosZ.GetName())
684 outdir.WriteTObject(MyHistoFiller.bis_mdtDigitStationEta, MyHistoFiller.bis_mdtDigitStationEta.GetName())
685 outdir.WriteTObject(MyHistoFiller.bis_mdtDigitStationPhi, MyHistoFiller.bis_mdtDigitStationPhi.GetName())
686 outdir.WriteTObject(MyHistoFiller.bis_mdtDigitglobalX, MyHistoFiller.bis_mdtDigitglobalX.GetName())
687 outdir.WriteTObject(MyHistoFiller.bis_mdtDigitglobalY, MyHistoFiller.bis_mdtDigitglobalY.GetName())
688 outdir.WriteTObject(MyHistoFiller.bis_mdtDigitglobalZ, MyHistoFiller.bis_mdtDigitglobalZ.GetName())
689 outdir.WriteTObject(MyHistoFiller.bis_mdtDigitNumberOfMultilayers, MyHistoFiller.bis_mdtDigitNumberOfMultilayers.GetName())
690 outdir.WriteTObject(MyHistoFiller.bis_mdtDigitTube, MyHistoFiller.bis_mdtDigitTube.GetName())
691 outdir.WriteTObject(MyHistoFiller.bis_mdtDigitTubeLayer, MyHistoFiller.bis_mdtDigitTubeLayer.GetName())
692 outdir.WriteTObject(MyHistoFiller.bis_mdtDigitMultilayer, MyHistoFiller.bis_mdtDigitMultilayer.GetName())
693 outdir.WriteTObject(MyHistoFiller.bis_mdtDigitTime, MyHistoFiller.bis_mdtDigitTime.GetName())
694 outdir.WriteTObject(MyHistoFiller.bis_mdtDigitCharge, MyHistoFiller.bis_mdtDigitCharge.GetName())
695 outdir.WriteTObject(MyHistoFiller.bis_mdtDigitChannel, MyHistoFiller.bis_mdtDigitChannel.GetName())
696 outdir.WriteTObject(MyHistoFiller.bis_mdtDigitglobalXY, MyHistoFiller.bis_mdtDigitglobalXY.GetName())
697 outdir.WriteTObject(MyHistoFiller.bis_mdtDigitglobalZY, MyHistoFiller.bis_mdtDigitglobalZY.GetName())
698 outdir.WriteTObject(MyHistoFiller.bis_mdtDigitglobalZX, MyHistoFiller.bis_mdtDigitglobalZX.GetName())
699 outdir.WriteTObject(MyHistoFiller.bis_mdtDigitlocalTubePosXY, MyHistoFiller.bis_mdtDigitlocalTubePosXY.GetName())
700 outdir.WriteTObject(MyHistoFiller.bis_mdtDigitlocalTubePosZY, MyHistoFiller.bis_mdtDigitlocalTubePosZY.GetName())
701 outdir.WriteTObject(MyHistoFiller.bis_mdtDigitlocalTubePosZX, MyHistoFiller.bis_mdtDigitlocalTubePosZX.GetName())
702
703
704 if self.__chamber_name == "CSC_Sim":
705 outdir.WriteTObject(MyHistoFiller.cscGlobalX, MyHistoFiller.cscGlobalX.GetName())
706 outdir.WriteTObject(MyHistoFiller.cscGlobalY, MyHistoFiller.cscGlobalY.GetName())
707 outdir.WriteTObject(MyHistoFiller.cscGlobalZ, MyHistoFiller.cscGlobalZ.GetName())
708 outdir.WriteTObject(MyHistoFiller.cscWireLayer, MyHistoFiller.cscWireLayer.GetName())
709 outdir.WriteTObject(MyHistoFiller.cscStrip, MyHistoFiller.cscStrip.GetName())
710 outdir.WriteTObject(MyHistoFiller.cscglobalTime, MyHistoFiller.cscglobalTime.GetName())
711 outdir.WriteTObject(MyHistoFiller.cscKineticEnergy, MyHistoFiller.cscKineticEnergy.GetName())
712 outdir.WriteTObject(MyHistoFiller.cscDepositEnergy, MyHistoFiller.cscDepositEnergy.GetName())
713 outdir.WriteTObject(MyHistoFiller.cscSimStationEta, MyHistoFiller.cscSimStationEta.GetName())
714 outdir.WriteTObject(MyHistoFiller.cscSimStationPhi, MyHistoFiller.cscSimStationPhi.GetName())
715 outdir.WriteTObject(MyHistoFiller.cscSimChamberLayer, MyHistoFiller.cscSimChamberLayer.GetName())
716
717
718 if self.__chamber_name == "CSC_Digit":
719 outdir.WriteTObject(MyHistoFiller.CSCDigitStationEta, MyHistoFiller.CSCDigitStationEta.GetName())
720 outdir.WriteTObject(MyHistoFiller.CSCDigitStationPhi, MyHistoFiller.CSCDigitStationPhi.GetName())
721 outdir.WriteTObject(MyHistoFiller.CSCDigitWireLayer, MyHistoFiller.CSCDigitWireLayer.GetName())
722 outdir.WriteTObject(MyHistoFiller.CSCDigitlocalPosX, MyHistoFiller.CSCDigitlocalPosX.GetName())
723 outdir.WriteTObject(MyHistoFiller.CSCDigitglobalPosX, MyHistoFiller.CSCDigitglobalPosX.GetName())
724 outdir.WriteTObject(MyHistoFiller.CSCDigitglobalPosY, MyHistoFiller.CSCDigitglobalPosY.GetName())
725 outdir.WriteTObject(MyHistoFiller.CSCDigitglobalPosZ, MyHistoFiller.CSCDigitglobalPosZ.GetName())
726
727 if self.__chamber_name == "CSC_SDO":
728 outdir.WriteTObject(MyHistoFiller.CSCSDOStationEta, MyHistoFiller.CSCSDOStationEta.GetName())
729 outdir.WriteTObject(MyHistoFiller.CSCSDOStationPhi, MyHistoFiller.CSCSDOStationPhi.GetName())
730 outdir.WriteTObject(MyHistoFiller.CSCSDOChamberLayer, MyHistoFiller.CSCSDOChamberLayer.GetName())
731 outdir.WriteTObject(MyHistoFiller.CSCSDOWireLayer, MyHistoFiller.CSCSDOWireLayer.GetName())
732 outdir.WriteTObject(MyHistoFiller.CSCSDOStrip, MyHistoFiller.CSCSDOStrip.GetName())
733 outdir.WriteTObject(MyHistoFiller.CSCSDOBarcode, MyHistoFiller.CSCSDOBarcode.GetName())
734 outdir.WriteTObject(MyHistoFiller.CSCSDOGlobalX, MyHistoFiller.CSCSDOGlobalX.GetName())
735 outdir.WriteTObject(MyHistoFiller.CSCSDOGlobalY, MyHistoFiller.CSCSDOGlobalY.GetName())
736 outdir.WriteTObject(MyHistoFiller.CSCSDOGlobalZ, MyHistoFiller.CSCSDOGlobalZ.GetName())
737 outdir.WriteTObject(MyHistoFiller.CSCSDOLocalX, MyHistoFiller.CSCSDOLocalX.GetName())
738 outdir.WriteTObject(MyHistoFiller.CSCSDOLocalY, MyHistoFiller.CSCSDOLocalY.GetName())
739
740 if self.__chamber_name == "CSC_RDO":
741 outdir.WriteTObject(MyHistoFiller.CSCRDOStationEta, MyHistoFiller.CSCRDOStationEta.GetName())
742 outdir.WriteTObject(MyHistoFiller.CSCRDOStationPhi, MyHistoFiller.CSCRDOStationPhi.GetName())
743 outdir.WriteTObject(MyHistoFiller.CSCRDOChamberLayer, MyHistoFiller.CSCRDOChamberLayer.GetName())
744 outdir.WriteTObject(MyHistoFiller.CSCRDOWireLayer, MyHistoFiller.CSCRDOWireLayer.GetName())
745 outdir.WriteTObject(MyHistoFiller.CSCRDOStrip, MyHistoFiller.CSCRDOStrip.GetName())
746 outdir.WriteTObject(MyHistoFiller.CSCRDOmeasurePhi, MyHistoFiller.CSCRDOmeasurePhi.GetName())
747 outdir.WriteTObject(MyHistoFiller.CSCRDOTime, MyHistoFiller.CSCRDOTime.GetName())
748 outdir.WriteTObject(MyHistoFiller.CSCRDOGlobalX, MyHistoFiller.CSCRDOGlobalX.GetName())
749 outdir.WriteTObject(MyHistoFiller.CSCRDOGlobalY, MyHistoFiller.CSCRDOGlobalY.GetName())
750 outdir.WriteTObject(MyHistoFiller.CSCRDOGlobalZ, MyHistoFiller.CSCRDOGlobalZ.GetName())
751 outdir.WriteTObject(MyHistoFiller.CSCRDOLocalX, MyHistoFiller.CSCRDOLocalX.GetName())
752 outdir.WriteTObject(MyHistoFiller.CSCRDOLocalY, MyHistoFiller.CSCRDOLocalY.GetName())
753
754 if self.__chamber_name == "CSC_PRD":
755 outdir.WriteTObject(MyHistoFiller.CSCPRDStationEta, MyHistoFiller.CSCPRDStationEta.GetName())
756 outdir.WriteTObject(MyHistoFiller.CSCPRDStationPhi, MyHistoFiller.CSCPRDStationPhi.GetName())
757 outdir.WriteTObject(MyHistoFiller.CSCPRDChamberLayer, MyHistoFiller.CSCPRDChamberLayer.GetName())
758 outdir.WriteTObject(MyHistoFiller.CSCPRDWireLayer, MyHistoFiller.CSCPRDWireLayer.GetName())
759 outdir.WriteTObject(MyHistoFiller.CSCPRDStrip, MyHistoFiller.CSCPRDStrip.GetName())
760 outdir.WriteTObject(MyHistoFiller.CSCPRDmeasurePhi, MyHistoFiller.CSCPRDmeasurePhi.GetName())
761 outdir.WriteTObject(MyHistoFiller.CSCPRDTime, MyHistoFiller.CSCPRDTime.GetName())
762 outdir.WriteTObject(MyHistoFiller.CSCPRDCharge, MyHistoFiller.CSCPRDCharge.GetName())
763 outdir.WriteTObject(MyHistoFiller.CSCPRDGlobalX, MyHistoFiller.CSCPRDGlobalX.GetName())
764 outdir.WriteTObject(MyHistoFiller.CSCPRDGlobalY, MyHistoFiller.CSCPRDGlobalY.GetName())
765 outdir.WriteTObject(MyHistoFiller.CSCPRDGlobalZ, MyHistoFiller.CSCPRDGlobalZ.GetName())
766 outdir.WriteTObject(MyHistoFiller.CSCPRDLocalX, MyHistoFiller.CSCPRDLocalX.GetName())
767
768
769 if self.__chamber_name == "TGC_Sim":
770 outdir.WriteTObject(MyHistoFiller.tgcLocalX, MyHistoFiller.tgcLocalX.GetName())
771 outdir.WriteTObject(MyHistoFiller.tgcLocalY, MyHistoFiller.tgcLocalY.GetName())
772 outdir.WriteTObject(MyHistoFiller.tgcLocalZ, MyHistoFiller.tgcLocalZ.GetName())
773 outdir.WriteTObject(MyHistoFiller.tgcStationEta, MyHistoFiller.tgcStationEta.GetName())
774 outdir.WriteTObject(MyHistoFiller.tgcStationPhi, MyHistoFiller.tgcStationPhi.GetName())
775 outdir.WriteTObject(MyHistoFiller.tgcGlobalX, MyHistoFiller.tgcGlobalX.GetName())
776 outdir.WriteTObject(MyHistoFiller.tgcGlobalY, MyHistoFiller.tgcGlobalY.GetName())
777 outdir.WriteTObject(MyHistoFiller.tgcGlobalZ, MyHistoFiller.tgcGlobalZ.GetName())
778 outdir.WriteTObject(MyHistoFiller.tgcGasGap, MyHistoFiller.tgcGasGap.GetName())
779 outdir.WriteTObject(MyHistoFiller.tgcChannel, MyHistoFiller.tgcChannel.GetName())
780 outdir.WriteTObject(MyHistoFiller.tgcGlobalTime, MyHistoFiller.tgcGlobalTime.GetName())
781 outdir.WriteTObject(MyHistoFiller.tgcKineticEnergy, MyHistoFiller.tgcKineticEnergy.GetName())
782 outdir.WriteTObject(MyHistoFiller.tgcDepositEnergy, MyHistoFiller.tgcDepositEnergy.GetName())
783
784 if self.__chamber_name == "TGC_Digit":
785 outdir.WriteTObject(MyHistoFiller.TGCDigitStationEta, MyHistoFiller.TGCDigitStationEta.GetName())
786 outdir.WriteTObject(MyHistoFiller.TGCDigitStationPhi, MyHistoFiller.TGCDigitStationPhi.GetName())
787 outdir.WriteTObject(MyHistoFiller.TGCDigitGasGap, MyHistoFiller.TGCDigitGasGap.GetName())
788 outdir.WriteTObject(MyHistoFiller.TGCDigitChannel, MyHistoFiller.TGCDigitChannel.GetName())
789 outdir.WriteTObject(MyHistoFiller.TGCDigitlocalX, MyHistoFiller.TGCDigitlocalX.GetName())
790 outdir.WriteTObject(MyHistoFiller.TGCDigitlocalY, MyHistoFiller.TGCDigitlocalY.GetName())
791 outdir.WriteTObject(MyHistoFiller.TGCDigitglobalX, MyHistoFiller.TGCDigitglobalX.GetName())
792 outdir.WriteTObject(MyHistoFiller.TGCDigitglobalY, MyHistoFiller.TGCDigitglobalY.GetName())
793 outdir.WriteTObject(MyHistoFiller.TGCDigitglobalZ, MyHistoFiller.TGCDigitglobalZ.GetName())
794
795 if self.__chamber_name == "TGC_SDO":
796 outdir.WriteTObject(MyHistoFiller.TGCSDOStationEta, MyHistoFiller.TGCSDOStationEta.GetName())
797 outdir.WriteTObject(MyHistoFiller.TGCSDOStationPhi, MyHistoFiller.TGCSDOStationPhi.GetName())
798 outdir.WriteTObject(MyHistoFiller.TGCSDOGasGap, MyHistoFiller.TGCSDOGasGap.GetName())
799 outdir.WriteTObject(MyHistoFiller.TGCSDOChannel, MyHistoFiller.TGCSDOChannel.GetName())
800 outdir.WriteTObject(MyHistoFiller.TGCSDOWord, MyHistoFiller.TGCSDOWord.GetName())
801 outdir.WriteTObject(MyHistoFiller.TGCSDOBarcode, MyHistoFiller.TGCSDOBarcode.GetName())
802 outdir.WriteTObject(MyHistoFiller.TGCSDOGlobalX, MyHistoFiller.TGCSDOGlobalX.GetName())
803 outdir.WriteTObject(MyHistoFiller.TGCSDOGlobalY, MyHistoFiller.TGCSDOGlobalY.GetName())
804 outdir.WriteTObject(MyHistoFiller.TGCSDOGlobalZ, MyHistoFiller.TGCSDOGlobalZ.GetName())
805 outdir.WriteTObject(MyHistoFiller.TGCSDOLocalX, MyHistoFiller.TGCSDOLocalX.GetName())
806 outdir.WriteTObject(MyHistoFiller.TGCSDOLocalY, MyHistoFiller.TGCSDOLocalY.GetName())
807 outdir.WriteTObject(MyHistoFiller.TGCSDOGlotime, MyHistoFiller.TGCSDOGlotime.GetName())
808
809 if self.__chamber_name == "TGC_RDO":
810 outdir.WriteTObject(MyHistoFiller.TGCRDOStationEta, MyHistoFiller.TGCRDOStationEta.GetName())
811 outdir.WriteTObject(MyHistoFiller.TGCRDOStationPhi, MyHistoFiller.TGCRDOStationPhi.GetName())
812 outdir.WriteTObject(MyHistoFiller.TGCRDOGasGap, MyHistoFiller.TGCRDOGasGap.GetName())
813 outdir.WriteTObject(MyHistoFiller.TGCRDOChannel, MyHistoFiller.TGCRDOChannel.GetName())
814 outdir.WriteTObject(MyHistoFiller.TGCRDOGlobalX, MyHistoFiller.TGCRDOGlobalX.GetName())
815 outdir.WriteTObject(MyHistoFiller.TGCRDOGlobalY, MyHistoFiller.TGCRDOGlobalY.GetName())
816 outdir.WriteTObject(MyHistoFiller.TGCRDOGlobalZ, MyHistoFiller.TGCRDOGlobalZ.GetName())
817 outdir.WriteTObject(MyHistoFiller.TGCRDOLocalX, MyHistoFiller.TGCRDOLocalX.GetName())
818
819 if self.__chamber_name == "TGC_PRD":
820 outdir.WriteTObject(MyHistoFiller.TGCPRDStationEta, MyHistoFiller.TGCPRDStationEta.GetName())
821 outdir.WriteTObject(MyHistoFiller.TGCPRDStationPhi, MyHistoFiller.TGCPRDStationPhi.GetName())
822 outdir.WriteTObject(MyHistoFiller.TGCPRDGasGap, MyHistoFiller.TGCPRDGasGap.GetName())
823 outdir.WriteTObject(MyHistoFiller.TGCPRDChannel, MyHistoFiller.TGCPRDChannel.GetName())
824 outdir.WriteTObject(MyHistoFiller.TGCPRDisStrip, MyHistoFiller.TGCPRDisStrip.GetName())
825 outdir.WriteTObject(MyHistoFiller.TGCPRDGlobalX, MyHistoFiller.TGCPRDGlobalX.GetName())
826 outdir.WriteTObject(MyHistoFiller.TGCPRDGlobalY, MyHistoFiller.TGCPRDGlobalY.GetName())
827 outdir.WriteTObject(MyHistoFiller.TGCPRDGlobalZ, MyHistoFiller.TGCPRDGlobalZ.GetName())
828 outdir.WriteTObject(MyHistoFiller.TGCPRDLocalX, MyHistoFiller.TGCPRDLocalX.GetName())
829
830
831 if self.__chamber_name == "MM_Sim":
832 outdir.WriteTObject(MyHistoFiller.mmGlobalX, MyHistoFiller.mmGlobalX.GetName())
833 outdir.WriteTObject(MyHistoFiller.mmGlobalY, MyHistoFiller.mmGlobalY.GetName())
834 outdir.WriteTObject(MyHistoFiller.mmGlobalZ, MyHistoFiller.mmGlobalZ.GetName())
835 outdir.WriteTObject(MyHistoFiller.mmStationEta, MyHistoFiller.mmStationEta.GetName())
836 outdir.WriteTObject(MyHistoFiller.mmStationPhi, MyHistoFiller.mmStationPhi.GetName())
837
838 if self.__chamber_name == "sTGC_Sim":
839 outdir.WriteTObject(MyHistoFiller.stgcGlobalX, MyHistoFiller.stgcGlobalX.GetName())
840 outdir.WriteTObject(MyHistoFiller.stgcGlobalY, MyHistoFiller.stgcGlobalY.GetName())
841 outdir.WriteTObject(MyHistoFiller.stgcGlobalZ, MyHistoFiller.stgcGlobalZ.GetName())
842 outdir.WriteTObject(MyHistoFiller.stgcStationEta, MyHistoFiller.stgcStationEta.GetName())
843 outdir.WriteTObject(MyHistoFiller.stgcStationPhi, MyHistoFiller.stgcStationPhi.GetName())
@ Eta
Definition RPCdef.h:8
__init__(self, chamber_name="", eta_sel=None, sector_sel=None)