ATLAS Offline Software
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 
5 import os, sys, ROOT
6 
7 class MyHistoFiller(object):
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())
DCubeHistograms.MyHistoFiller.__sector_sel
__sector_sel
Definition: DCubeHistograms.py:265
DCubeHistograms.MyHistoFiller.__init__
def __init__(self, chamber_name="", eta_sel=None, sector_sel=None)
Definition: DCubeHistograms.py:262
DCubeHistograms.MyHistoFiller
Definition: DCubeHistograms.py:7
DCubeHistograms.MyHistoFiller.fill
def fill(self, TTree, n)
Definition: DCubeHistograms.py:274
DCubeHistograms.MyHistoFiller.__chamber_name
__chamber_name
Definition: DCubeHistograms.py:263
pickleTool.object
object
Definition: pickleTool.py:29
DCubeHistograms.MyHistoFiller.write
def write(self, outdir)
Definition: DCubeHistograms.py:593
DCubeHistograms.MyHistoFiller.Eta
def Eta(eta)
Definition: DCubeHistograms.py:268
DCubeHistograms.MyHistoFiller.__eta_sel
__eta_sel
Definition: DCubeHistograms.py:264
Generated on Thu Oct 30 2025 21:08:52 for ATLAS Offline Software by   doxygen 1.8.18