Back to home page

Project CMSSW displayed by LXR

 
 

    


File indexing completed on 2024-04-06 12:26:03

0001 // This function makes plots from root file output of CSCValidation, input_file
0002 // Pass it to root e.g. see makePlots.sh script
0003 // It requires local file "myFunctions.C"
0004 // Original author: Andy Kubik (NWU)
0005 // - Please contact CSC DPG for current status - 12.07.2022
0006 
0007 // - includes rechit and segment position plots
0008 // - does not include segment time plots - MUST BE ADDED
0009 
0010 #include "myFunctions.C"
0011 
0012 void makePlots(std::string input_file) {
0013   extern TFile* OpenFiles(std::string path);
0014   extern void makeCSCOccupancy(std::string histoname, TFile * f1, std::string histotitle, std::string savename);
0015   extern void Draw2DTempPlot(
0016       std::string histo, TFile * f1, bool includeME11, std::string savename, bool hasLabels = false);
0017   extern void make1DPlot(
0018       std::string histoname, TFile * f1, std::string histotitle, int statoption, std::string savename);
0019   extern void make1DPlot(std::string histoname,
0020                          TFile * f1,
0021                          std::string histotitle,
0022                          std::string xtitle,
0023                          std::string ytitle,
0024                          int statoption,
0025                          std::string savename);
0026   extern void printEmptyChambers(std::string histoname, std::string oname, TFile * f);
0027   extern void GlobalPosfromTree(
0028       std::string graphname, TFile * f1, int endcap, int station, std::string type, std::string savename);
0029   extern void NikolaiPlots(TFile * f_in, int flag);
0030   extern void makeEffGif(std::string histoname, TFile * f1, std::string histotitle, std::string savename);
0031   extern void Draw2DEfficiency(std::string histo, TFile * f1, std::string title, std::string savename);
0032   extern void make2DPlot(
0033       std::string histoname, TFile * f1, std::string histotitle, int statoption, std::string savename);
0034   extern void make2DPlot(std::string histoname,
0035                          TFile * f1,
0036                          std::string histotitle,
0037                          std::string xtitle,
0038                          std::string ytitle,
0039                          int statoption,
0040                          std::string savename);
0041   extern void makeProfile(
0042       std::string histoname, TFile * f1, std::string histotitle, int statoption, std::string savename);
0043   extern void makeProfile(std::string histoname,
0044                           TFile * f1,
0045                           std::string histotitle,
0046                           std::string xtitle,
0047                           std::string ytitle,
0048                           int statoption,
0049                           std::string savename);
0050 
0051   //  gROOT->Reset();
0052   //  gROOT->ProcessLine(".L myFunctions.C");
0053 
0054   TFile* f1;
0055   f1 = OpenFiles(input_file);
0056 
0057   //produce occupancy plots
0058   makeCSCOccupancy("GeneralHists/hCSCOccupancy", f1, "CSC Occupancy", "CSCOccupancy.png");
0059   Draw2DTempPlot("Digis/hOWires", f1, true, "hOWires.png");
0060   Draw2DTempPlot("Digis/hOStrips", f1, true, "hOStrips.png");
0061   Draw2DTempPlot("recHits/hORecHits", f1, true, "hORecHits.png");
0062   Draw2DTempPlot("Segments/hOSegments", f1, true, "hOSegments.png");
0063   make1DPlot("Digis/hOStripSerial", f1, "Strip Occupancy by Chamber Serial", 10, "hOStripSerial.png");
0064   make1DPlot("Digis/hOWireSerial", f1, "Wire Occupancy by Chamber Serial", 10, "hOWireSerial.png");
0065   make1DPlot("Segments/hOSegmentsSerial", f1, "Segment Occupancy by Chamber Serial", 10, "hOSegmentsSerial.png");
0066   make1DPlot("recHits/hORecHitsSerial", f1, "recHit Occupancy by Chamber Serial", 10, "hORecHitsSerial.png");
0067 
0068   //Print dead chamber lists
0069   printEmptyChambers("Digis/hOWires", "wire digis", f1);
0070   printEmptyChambers("Digis/hOStrips", "strip digis", f1);
0071   printEmptyChambers("recHits/hORecHits", "rechits", f1);
0072 
0073   //Make global position graphs from trees
0074   GlobalPosfromTree("Global recHit positions ME+1", f1, 1, 1, "rechit", "rHglobal_station_+1.png");
0075   GlobalPosfromTree("Global recHit positions ME+2", f1, 1, 2, "rechit", "rHglobal_station_+2.png");
0076   GlobalPosfromTree("Global recHit positions ME+3", f1, 1, 3, "rechit", "rHglobal_station_+3.png");
0077   GlobalPosfromTree("Global recHit positions ME+4", f1, 1, 4, "rechit", "rHglobal_station_+4.png");
0078   GlobalPosfromTree("Global recHit positions ME-1", f1, 2, 1, "rechit", "rHglobal_station_-1.png");
0079   GlobalPosfromTree("Global recHit positions ME-2", f1, 2, 2, "rechit", "rHglobal_station_-2.png");
0080   GlobalPosfromTree("Global recHit positions ME-3", f1, 2, 3, "rechit", "rHglobal_station_-3.png");
0081   GlobalPosfromTree("Global recHit positions ME-4", f1, 2, 4, "rechit", "rHglobal_station_-4.png");
0082   GlobalPosfromTree("Global Segment positions ME+1", f1, 1, 1, "segment", "Sglobal_station_+1.png");
0083   GlobalPosfromTree("Global Segment positions ME+2", f1, 1, 2, "segment", "Sglobal_station_+2.png");
0084   GlobalPosfromTree("Global Segment positions ME+3", f1, 1, 3, "segment", "Sglobal_station_+3.png");
0085   GlobalPosfromTree("Global Segment positions ME+4", f1, 1, 4, "segment", "Sglobal_station_+4.png");
0086   GlobalPosfromTree("Global Segment positions ME-1", f1, 2, 1, "segment", "Sglobal_station_-1.png");
0087   GlobalPosfromTree("Global Segment positions ME-2", f1, 2, 2, "segment", "Sglobal_station_-2.png");
0088   GlobalPosfromTree("Global Segment positions ME-3", f1, 2, 3, "segment", "Sglobal_station_-3.png");
0089   GlobalPosfromTree("Global Segment positions ME-4", f1, 2, 4, "segment", "Sglobal_station_-4.png");
0090 
0091   //Nikolai's plots
0092   NikolaiPlots(f1, 1);
0093   NikolaiPlots(f1, 2);
0094   NikolaiPlots(f1, 3);
0095   NikolaiPlots(f1, 4);
0096 
0097   //produce number of X per event plots
0098   make1DPlot("Digis/hStripNFired", f1, "Fired Strips per Event", 1110, "Digis_hStripNFired.png");
0099   make1DPlot("Digis/hWirenGroupsTotal", f1, "Fired Wires per Event", 1110, "Digis_hWirenGroupsTotal.png");
0100   make1DPlot("recHits/hRHnrechits", f1, "RecHits per Event", 1110, "recHits_hRHnrechits.png");
0101   make1DPlot("Segments/hSnSegments", f1, "Segments per Event", 1110, "Segments_hSnSegments.png");
0102 
0103   //efficiency plots
0104   makeEffGif("Efficiency/hRHSTE", f1, "RecHit Efficiency", "Efficiency_hRHEff.png");
0105   makeEffGif("Efficiency/hSSTE", f1, "Segment Efficiency", "Efficiency_hSEff.png");
0106   Draw2DEfficiency("Efficiency/hRHSTE2", f1, "RecHit Efficiency 2D", "Efficiency_hRHEff2.png");
0107   Draw2DEfficiency("Efficiency/hSSTE2", f1, "Segment Efficiency 2D", "Efficiency_hSEff2.png");
0108   Draw2DEfficiency("Efficiency/hWireSTE2", f1, "Wire Efficiency 2D", "Efficiency_hWireEff2.png");
0109   Draw2DEfficiency("Efficiency/hStripSTE2", f1, "Strip Efficiency 2D", "Efficiency_hStripEff2.png");
0110   Draw2DTempPlot("Efficiency/hSensitiveAreaEvt", f1, false, "Efficiency_hEvts2.png");
0111 
0112   //produce wire timing plots
0113   make1DPlot("Digis/hWireTBin+11", f1, "Wire TimeBin Fired ME+1/1", 1110, "Digis_hWireTBin+11.png");
0114   make1DPlot("Digis/hWireTBin+12", f1, "Wire TimeBin Fired ME+1/2", 1110, "Digis_hWireTBin+12.png");
0115   make1DPlot("Digis/hWireTBin+13", f1, "Wire TimeBin Fired ME+1/3", 1110, "Digis_hWireTBin+13.png");
0116   make1DPlot("Digis/hWireTBin+21", f1, "Wire TimeBin Fired ME+2/1", 1110, "Digis_hWireTBin+21.png");
0117   make1DPlot("Digis/hWireTBin+22", f1, "Wire TimeBin Fired ME+2/2", 1110, "Digis_hWireTBin+22.png");
0118   make1DPlot("Digis/hWireTBin+31", f1, "Wire TimeBin Fired ME+3/1", 1110, "Digis_hWireTBin+31.png");
0119   make1DPlot("Digis/hWireTBin+32", f1, "Wire TimeBin Fired ME+3/2", 1110, "Digis_hWireTBin+32.png");
0120   make1DPlot("Digis/hWireTBin+41", f1, "Wire TimeBin Fired ME+4/1", 1110, "Digis_hWireTBin+41.png");
0121   make1DPlot("Digis/hWireTBin-11", f1, "Wire TimeBin Fired ME-1/1", 1110, "Digis_hWireTBin-11.png");
0122   make1DPlot("Digis/hWireTBin-12", f1, "Wire TimeBin Fired ME-1/2", 1110, "Digis_hWireTBin-12.png");
0123   make1DPlot("Digis/hWireTBin-13", f1, "Wire TimeBin Fired ME-1/3", 1110, "Digis_hWireTBin-13.png");
0124   make1DPlot("Digis/hWireTBin-21", f1, "Wire TimeBin Fired ME-2/1", 1110, "Digis_hWireTBin-21.png");
0125   make1DPlot("Digis/hWireTBin-22", f1, "Wire TimeBin Fired ME-2/2", 1110, "Digis_hWireTBin-22.png");
0126   make1DPlot("Digis/hWireTBin-31", f1, "Wire TimeBin Fired ME-3/1", 1110, "Digis_hWireTBin-31.png");
0127   make1DPlot("Digis/hWireTBin-32", f1, "Wire TimeBin Fired ME-3/2", 1110, "Digis_hWireTBin-32.png");
0128   make1DPlot("Digis/hWireTBin-41", f1, "Wire TimeBin Fired ME-4/1", 1110, "Digis_hWireTBin-41.png");
0129 
0130   //produce pedestal noise plots
0131   make1DPlot("PedestalNoise/hStripPedME+11",
0132              f1,
0133              "Pedestal Noise Distribution ME+1/1b",
0134              1110,
0135              "PedestalNoise_hStripPedME+11.png");
0136   make1DPlot("PedestalNoise/hStripPedME+14",
0137              f1,
0138              "Pedestal Noise Distribution ME+1/1a",
0139              1110,
0140              "PedestalNoise_hStripPedME+11a.png");
0141   make1DPlot("PedestalNoise/hStripPedME+12",
0142              f1,
0143              "Pedestal Noise Distribution ME+1/2",
0144              1110,
0145              "PedestalNoise_hStripPedME+12.png");
0146   make1DPlot("PedestalNoise/hStripPedME+13",
0147              f1,
0148              "Pedestal Noise Distribution ME+1/3",
0149              1110,
0150              "PedestalNoise_hStripPedME+13.png");
0151   make1DPlot("PedestalNoise/hStripPedME+21",
0152              f1,
0153              "Pedestal Noise Distribution ME+2/1",
0154              1110,
0155              "PedestalNoise_hStripPedME+21.png");
0156   make1DPlot("PedestalNoise/hStripPedME+22",
0157              f1,
0158              "Pedestal Noise Distribution ME+2/2",
0159              1110,
0160              "PedestalNoise_hStripPedME+22.png");
0161   make1DPlot("PedestalNoise/hStripPedME+31",
0162              f1,
0163              "Pedestal Noise Distribution ME+3/1",
0164              1110,
0165              "PedestalNoise_hStripPedME+31.png");
0166   make1DPlot("PedestalNoise/hStripPedME+32",
0167              f1,
0168              "Pedestal Noise Distribution ME+3/2",
0169              1110,
0170              "PedestalNoise_hStripPedME+32.png");
0171   make1DPlot("PedestalNoise/hStripPedME+41",
0172              f1,
0173              "Pedestal Noise Distribution ME+4/1",
0174              1110,
0175              "PedestalNoise_hStripPedME+41.png");
0176   make1DPlot("PedestalNoise/hStripPedME-11",
0177              f1,
0178              "Pedestal Noise Distribution ME-1/1b",
0179              1110,
0180              "PedestalNoise_hStripPedME-11.png");
0181   make1DPlot("PedestalNoise/hStripPedME-14",
0182              f1,
0183              "Pedestal Noise Distribution ME-1/1a",
0184              1110,
0185              "PedestalNoise_hStripPedME-11a.png");
0186   make1DPlot("PedestalNoise/hStripPedME-12",
0187              f1,
0188              "Pedestal Noise Distribution ME-1/2",
0189              1110,
0190              "PedestalNoise_hStripPedME-12.png");
0191   make1DPlot("PedestalNoise/hStripPedME-13",
0192              f1,
0193              "Pedestal Noise Distribution ME-1/3",
0194              1110,
0195              "PedestalNoise_hStripPedME-13.png");
0196   make1DPlot("PedestalNoise/hStripPedME-21",
0197              f1,
0198              "Pedestal Noise Distribution ME-2/1",
0199              1110,
0200              "PedestalNoise_hStripPedME-21.png");
0201   make1DPlot("PedestalNoise/hStripPedME-22",
0202              f1,
0203              "Pedestal Noise Distribution ME-2/2",
0204              1110,
0205              "PedestalNoise_hStripPedME-22.png");
0206   make1DPlot("PedestalNoise/hStripPedME-31",
0207              f1,
0208              "Pedestal Noise Distribution ME-3/1",
0209              1110,
0210              "PedestalNoise_hStripPedME-31.png");
0211   make1DPlot("PedestalNoise/hStripPedME-32",
0212              f1,
0213              "Pedestal Noise Distribution ME-3/2",
0214              1110,
0215              "PedestalNoise_hStripPedME-32.png");
0216   make1DPlot("PedestalNoise/hStripPedME-41",
0217              f1,
0218              "Pedestal Noise Distribution ME-4/1",
0219              1110,
0220              "PedestalNoise_hStripPedME-41.png");
0221 
0222   // resolution
0223   make1DPlot("Resolution/hSResid+11",
0224              f1,
0225              "Expected Position from Fit - Reconstructed, ME+1/1b",
0226              1110,
0227              "Resolution_hSResid+11.png");
0228   make1DPlot("Resolution/hSResid+12",
0229              f1,
0230              "Expected Position from Fit - Reconstructed, ME+1/2",
0231              1110,
0232              "Resolution_hSResid+12.png");
0233   make1DPlot("Resolution/hSResid+13",
0234              f1,
0235              "Expected Position from Fit - Reconstructed, ME+1/3",
0236              1110,
0237              "Resolution_hSResid+13.png");
0238   make1DPlot("Resolution/hSResid+14",
0239              f1,
0240              "Expected Position from Fit - Reconstructed, ME+1/1a",
0241              1110,
0242              "Resolution_hSResid+11a.png");
0243   make1DPlot("Resolution/hSResid+21",
0244              f1,
0245              "Expected Position from Fit - Reconstructed, ME+2/1",
0246              1110,
0247              "Resolution_hSResid+21.png");
0248   make1DPlot("Resolution/hSResid+22",
0249              f1,
0250              "Expected Position from Fit - Reconstructed, ME+2/2",
0251              1110,
0252              "Resolution_hSResid+22.png");
0253   make1DPlot("Resolution/hSResid+31",
0254              f1,
0255              "Expected Position from Fit - Reconstructed, ME+3/1",
0256              1110,
0257              "Resolution_hSResid+31.png");
0258   make1DPlot("Resolution/hSResid+32",
0259              f1,
0260              "Expected Position from Fit - Reconstructed, ME+3/2",
0261              1110,
0262              "Resolution_hSResid+32.png");
0263   make1DPlot("Resolution/hSResid+41",
0264              f1,
0265              "Expected Position from Fit - Reconstructed, ME+4/1",
0266              1110,
0267              "Resolution_hSResid+41.png");
0268   make1DPlot("Resolution/hSResid-11",
0269              f1,
0270              "Expected Position from Fit - Reconstructed, ME-1/1b",
0271              1110,
0272              "Resolution_hSResid-11.png");
0273   make1DPlot("Resolution/hSResid-12",
0274              f1,
0275              "Expected Position from Fit - Reconstructed, ME-1/2",
0276              1110,
0277              "Resolution_hSResid-12.png");
0278   make1DPlot("Resolution/hSResid-13",
0279              f1,
0280              "Expected Position from Fit - Reconstructed, ME-1/3",
0281              1110,
0282              "Resolution_hSResid-13.png");
0283   make1DPlot("Resolution/hSResid-14",
0284              f1,
0285              "Expected Position from Fit - Reconstructed, ME-1/1a",
0286              1110,
0287              "Resolution_hSResid-11a.png");
0288   make1DPlot("Resolution/hSResid-21",
0289              f1,
0290              "Expected Position from Fit - Reconstructed, ME-2/1",
0291              1110,
0292              "Resolution_hSResid-21.png");
0293   make1DPlot("Resolution/hSResid-22",
0294              f1,
0295              "Expected Position from Fit - Reconstructed, ME-2/2",
0296              1110,
0297              "Resolution_hSResid-22.png");
0298   make1DPlot("Resolution/hSResid-31",
0299              f1,
0300              "Expected Position from Fit - Reconstructed, ME-3/1",
0301              1110,
0302              "Resolution_hSResid-31.png");
0303   make1DPlot("Resolution/hSResid-32",
0304              f1,
0305              "Expected Position from Fit - Reconstructed, ME-3/2",
0306              1110,
0307              "Resolution_hSResid-32.png");
0308   make1DPlot("Resolution/hSResid-41",
0309              f1,
0310              "Expected Position from Fit - Reconstructed, ME-4/1",
0311              1110,
0312              "Resolution_hSResid-41.png");
0313 
0314   // rechit strip position
0315   make1DPlot("recHits/hRHstpos+11", f1, "Strip Position (ME+1/1b)", 1110, "recHits_hRHstpos+11.png");
0316   make1DPlot("recHits/hRHstpos+14", f1, "Strip Position (ME+1/1a)", 1110, "recHits_hRHstpos+11a.png");
0317   make1DPlot("recHits/hRHstpos+12", f1, "Strip Position (ME+1/2)", 1110, "recHits_hRHstpos+12.png");
0318   make1DPlot("recHits/hRHstpos+13", f1, "Strip Position (ME+1/3)", 1110, "recHits_hRHstpos+13.png");
0319   make1DPlot("recHits/hRHstpos+21", f1, "Strip Position (ME+2/1)", 1110, "recHits_hRHstpos+21.png");
0320   make1DPlot("recHits/hRHstpos+22", f1, "Strip Position (ME+2/2)", 1110, "recHits_hRHstpos+22.png");
0321   make1DPlot("recHits/hRHstpos+31", f1, "Strip Position (ME+3/1)", 1110, "recHits_hRHstpos+31.png");
0322   make1DPlot("recHits/hRHstpos+32", f1, "Strip Position (ME+3/2)", 1110, "recHits_hRHstpos+32.png");
0323   make1DPlot("recHits/hRHstpos+41", f1, "Strip Position (ME+4/1)", 1110, "recHits_hRHstpos+41.png");
0324   make1DPlot("recHits/hRHstpos-11", f1, "Strip Position (ME-1/1b)", 1110, "recHits_hRHstpos-11.png");
0325   make1DPlot("recHits/hRHstpos-14", f1, "Strip Position (ME-1/1a)", 1110, "recHits_hRHstpos-11a.png");
0326   make1DPlot("recHits/hRHstpos-12", f1, "Strip Position (ME-1/2)", 1110, "recHits_hRHstpos-12.png");
0327   make1DPlot("recHits/hRHstpos-13", f1, "Strip Position (ME-1/3)", 1110, "recHits_hRHstpos-13.png");
0328   make1DPlot("recHits/hRHstpos-21", f1, "Strip Position (ME-2/1)", 1110, "recHits_hRHstpos-21.png");
0329   make1DPlot("recHits/hRHstpos-22", f1, "Strip Position (ME-2/2)", 1110, "recHits_hRHstpos-22.png");
0330   make1DPlot("recHits/hRHstpos-31", f1, "Strip Position (ME-3/1)", 1110, "recHits_hRHstpos-31.png");
0331   make1DPlot("recHits/hRHstpos-32", f1, "Strip Position (ME-3/2)", 1110, "recHits_hRHstpos-32.png");
0332   make1DPlot("recHits/hRHstpos-41", f1, "Strip Position (ME-4/1)", 1110, "recHits_hRHstpos-41.png");
0333 
0334   // rechit timing
0335   make1DPlot("recHits/hRHTiming+11", f1, "RecHit Timing ME+1/1b", 1110, "recHits_hRHTiming+11.png");
0336   make1DPlot("recHits/hRHTiming+14", f1, "RecHit Timing ME+1/1a", 1110, "recHits_hRHTiming+11a.png");
0337   make1DPlot("recHits/hRHTiming+12", f1, "RecHit Timing ME+1/2", 1110, "recHits_hRHTiming+12.png");
0338   make1DPlot("recHits/hRHTiming+13", f1, "RecHit Timing ME+1/3", 1110, "recHits_hRHTiming+13.png");
0339   make1DPlot("recHits/hRHTiming+21", f1, "RecHit Timing ME+2/1", 1110, "recHits_hRHTiming+21.png");
0340   make1DPlot("recHits/hRHTiming+22", f1, "RecHit Timing ME+2/2", 1110, "recHits_hRHTiming+22.png");
0341   make1DPlot("recHits/hRHTiming+31", f1, "RecHit Timing ME+3/1", 1110, "recHits_hRHTiming+31.png");
0342   make1DPlot("recHits/hRHTiming+32", f1, "RecHit Timing ME+3/2", 1110, "recHits_hRHTiming+32.png");
0343   make1DPlot("recHits/hRHTiming+41", f1, "RecHit Timing ME+4/1", 1110, "recHits_hRHTiming+41.png");
0344   make1DPlot("recHits/hRHTiming-11", f1, "RecHit Timing ME-1/1b", 1110, "recHits_hRHTiming-11.png");
0345   make1DPlot("recHits/hRHTiming-14", f1, "RecHit Timing ME-1/1a", 1110, "recHits_hRHTiming-11a.png");
0346   make1DPlot("recHits/hRHTiming-12", f1, "RecHit Timing ME-1/2", 1110, "recHits_hRHTiming-12.png");
0347   make1DPlot("recHits/hRHTiming-13", f1, "RecHit Timing ME-1/3", 1110, "recHits_hRHTiming-13.png");
0348   make1DPlot("recHits/hRHTiming-21", f1, "RecHit Timing ME-2/1", 1110, "recHits_hRHTiming-21.png");
0349   make1DPlot("recHits/hRHTiming-22", f1, "RecHit Timing ME-2/2", 1110, "recHits_hRHTiming-22.png");
0350   make1DPlot("recHits/hRHTiming-31", f1, "RecHit Timing ME-3/1", 1110, "recHits_hRHTiming-31.png");
0351   make1DPlot("recHits/hRHTiming-32", f1, "RecHit Timing ME-3/2", 1110, "recHits_hRHTiming-32.png");
0352   make1DPlot("recHits/hRHTiming-41", f1, "RecHit Timing ME-4/1", 1110, "recHits_hRHTiming-41.png");
0353 
0354   // rechit charge
0355   make1DPlot("recHits/hRHSumQ+11", f1, "Sum 3x3 RecHit Charge ME+1/1b", 1110, "recHits_hRHSumQ+11.png");
0356   make1DPlot("recHits/hRHSumQ+14", f1, "Sum 3x3 RecHit Charge ME+1/1a", 1110, "recHits_hRHSumQ+11a.png");
0357   make1DPlot("recHits/hRHSumQ+12", f1, "Sum 3x3 RecHit Charge ME+1/2", 1110, "recHits_hRHSumQ+12.png");
0358   make1DPlot("recHits/hRHSumQ+13", f1, "Sum 3x3 RecHit Charge ME+1/3", 1110, "recHits_hRHSumQ+13.png");
0359   make1DPlot("recHits/hRHSumQ+21", f1, "Sum 3x3 RecHit Charge ME+2/1", 1110, "recHits_hRHSumQ+21.png");
0360   make1DPlot("recHits/hRHSumQ+22", f1, "Sum 3x3 RecHit Charge ME+2/2", 1110, "recHits_hRHSumQ+22.png");
0361   make1DPlot("recHits/hRHSumQ+31", f1, "Sum 3x3 RecHit Charge ME+3/1", 1110, "recHits_hRHSumQ+31.png");
0362   make1DPlot("recHits/hRHSumQ+32", f1, "Sum 3x3 RecHit Charge ME+3/2", 1110, "recHits_hRHSumQ+32.png");
0363   make1DPlot("recHits/hRHSumQ+41", f1, "Sum 3x3 RecHit Charge ME+4/1", 1110, "recHits_hRHSumQ+41.png");
0364   make1DPlot("recHits/hRHSumQ-11", f1, "Sum 3x3 RecHit Charge ME-1/1b", 1110, "recHits_hRHSumQ-11.png");
0365   make1DPlot("recHits/hRHSumQ-14", f1, "Sum 3x3 RecHit Charge ME-1/1a", 1110, "recHits_hRHSumQ-11a.png");
0366   make1DPlot("recHits/hRHSumQ-12", f1, "Sum 3x3 RecHit Charge ME-1/2", 1110, "recHits_hRHSumQ-12.png");
0367   make1DPlot("recHits/hRHSumQ-13", f1, "Sum 3x3 RecHit Charge ME-1/3", 1110, "recHits_hRHSumQ-13.png");
0368   make1DPlot("recHits/hRHSumQ-21", f1, "Sum 3x3 RecHit Charge ME-2/1", 1110, "recHits_hRHSumQ-21.png");
0369   make1DPlot("recHits/hRHSumQ-22", f1, "Sum 3x3 RecHit Charge ME-2/2", 1110, "recHits_hRHSumQ-22.png");
0370   make1DPlot("recHits/hRHSumQ-31", f1, "Sum 3x3 RecHit Charge ME-3/1", 1110, "recHits_hRHSumQ-31.png");
0371   make1DPlot("recHits/hRHSumQ-32", f1, "Sum 3x3 RecHit Charge ME-3/2", 1110, "recHits_hRHSumQ-32.png");
0372   make1DPlot("recHits/hRHSumQ-41", f1, "Sum 3x3 RecHit Charge ME-4/1", 1110, "recHits_hRHSumQ-41.png");
0373 
0374   make1DPlot("recHits/hRHRatioQ+11", f1, "Charge Ratio (Ql_Qr)/Qt ME+1/1b", 1110, "recHits_hRHRatioQ+11.png");
0375   make1DPlot("recHits/hRHRatioQ+14", f1, "Charge Ratio (Ql_Qr)/Qt ME+1/1a", 1110, "recHits_hRHRatioQ+11a.png");
0376   make1DPlot("recHits/hRHRatioQ+12", f1, "Charge Ratio (Ql_Qr)/Qt ME+1/2", 1110, "recHits_hRHRatioQ+12.png");
0377   make1DPlot("recHits/hRHRatioQ+13", f1, "Charge Ratio (Ql_Qr)/Qt ME+1/3", 1110, "recHits_hRHRatioQ+13.png");
0378   make1DPlot("recHits/hRHRatioQ+21", f1, "Charge Ratio (Ql_Qr)/Qt ME+2/1", 1110, "recHits_hRHRatioQ+21.png");
0379   make1DPlot("recHits/hRHRatioQ+22", f1, "Charge Ratio (Ql_Qr)/Qt ME+2/2", 1110, "recHits_hRHRatioQ+22.png");
0380   make1DPlot("recHits/hRHRatioQ+31", f1, "Charge Ratio (Ql_Qr)/Qt ME+3/1", 1110, "recHits_hRHRatioQ+31.png");
0381   make1DPlot("recHits/hRHRatioQ+32", f1, "Charge Ratio (Ql_Qr)/Qt ME+3/2", 1110, "recHits_hRHRatioQ+32.png");
0382   make1DPlot("recHits/hRHRatioQ+41", f1, "Charge Ratio (Ql_Qr)/Qt ME+4/1", 1110, "recHits_hRHRatioQ+41.png");
0383   make1DPlot("recHits/hRHRatioQ-11", f1, "Charge Ratio (Ql_Qr)/Qt ME-1/1b", 1110, "recHits_hRHRatioQ-11.png");
0384   make1DPlot("recHits/hRHRatioQ-14", f1, "Charge Ratio (Ql_Qr)/Qt ME-1/1a", 1110, "recHits_hRHRatioQ-11a.png");
0385   make1DPlot("recHits/hRHRatioQ-12", f1, "Charge Ratio (Ql_Qr)/Qt ME-1/2", 1110, "recHits_hRHRatioQ-12.png");
0386   make1DPlot("recHits/hRHRatioQ-13", f1, "Charge Ratio (Ql_Qr)/Qt ME-1/3", 1110, "recHits_hRHRatioQ-13.png");
0387   make1DPlot("recHits/hRHRatioQ-21", f1, "Charge Ratio (Ql_Qr)/Qt ME-2/1", 1110, "recHits_hRHRatioQ-21.png");
0388   make1DPlot("recHits/hRHRatioQ-22", f1, "Charge Ratio (Ql_Qr)/Qt ME-2/2", 1110, "recHits_hRHRatioQ-22.png");
0389   make1DPlot("recHits/hRHRatioQ-31", f1, "Charge Ratio (Ql_Qr)/Qt ME-3/1", 1110, "recHits_hRHRatioQ-31.png");
0390   make1DPlot("recHits/hRHRatioQ-32", f1, "Charge Ratio (Ql_Qr)/Qt ME-3/2", 1110, "recHits_hRHRatioQ-32.png");
0391   make1DPlot("recHits/hRHRatioQ-41", f1, "Charge Ratio (Ql_Qr)/Qt ME-4/1", 1110, "recHits_hRHRatioQ-41.png");
0392 
0393   //hits on a segment
0394   make1DPlot("Segments/hSnHits+11", f1, "N Hits on Segments ME+1/1b", 1110, "Segments_hSnHits+11.png");
0395   make1DPlot("Segments/hSnHits+14", f1, "N Hits on Segments ME+1/1a", 1110, "Segments_hSnHits+11a.png");
0396   make1DPlot("Segments/hSnHits+12", f1, "N Hits on Segments ME+1/2", 1110, "Segments_hSnHits+12.png");
0397   make1DPlot("Segments/hSnHits+13", f1, "N Hits on Segments ME+1/3", 1110, "Segments_hSnHits+13.png");
0398   make1DPlot("Segments/hSnHits+21", f1, "N Hits on Segments ME+2/1", 1110, "Segments_hSnHits+21.png");
0399   make1DPlot("Segments/hSnHits+22", f1, "N Hits on Segments ME+2/2", 1110, "Segments_hSnHits+22.png");
0400   make1DPlot("Segments/hSnHits+31", f1, "N Hits on Segments ME+3/1", 1110, "Segments_hSnHits+31.png");
0401   make1DPlot("Segments/hSnHits+32", f1, "N Hits on Segments ME+3/2", 1110, "Segments_hSnHits+32.png");
0402   make1DPlot("Segments/hSnHits+41", f1, "N Hits on Segments ME+4/1", 1110, "Segments_hSnHits+41.png");
0403   make1DPlot("Segments/hSnHits-11", f1, "N Hits on Segments ME-1/1b", 1110, "Segments_hSnHits-11.png");
0404   make1DPlot("Segments/hSnHits-14", f1, "N Hits on Segments ME-1/1a", 1110, "Segments_hSnHits-11a.png");
0405   make1DPlot("Segments/hSnHits-12", f1, "N Hits on Segments ME-1/2", 1110, "Segments_hSnHits-12.png");
0406   make1DPlot("Segments/hSnHits-13", f1, "N Hits on Segments ME-1/3", 1110, "Segments_hSnHits-13.png");
0407   make1DPlot("Segments/hSnHits-21", f1, "N Hits on Segments ME-2/1", 1110, "Segments_hSnHits-21.png");
0408   make1DPlot("Segments/hSnHits-22", f1, "N Hits on Segments ME-2/2", 1110, "Segments_hSnHits-22.png");
0409   make1DPlot("Segments/hSnHits-31", f1, "N Hits on Segments ME-3/1", 1110, "Segments_hSnHits-31.png");
0410   make1DPlot("Segments/hSnHits-32", f1, "N Hits on Segments ME-3/2", 1110, "Segments_hSnHits-32.png");
0411   make1DPlot("Segments/hSnHits-41", f1, "N Hits on Segments ME-4/1", 1110, "Segments_hSnHits-41.png");
0412 
0413   // segment chi2
0414   make1DPlot("Segments/hSChiSq+11", f1, "Segment Chi2/ndof ME+1/1b", 1110, "Segments_hSChiSq+11.png");
0415   make1DPlot("Segments/hSChiSq+14", f1, "Segment Chi2/ndof ME+1/1a", 1110, "Segments_hSChiSq+11a.png");
0416   make1DPlot("Segments/hSChiSq+12", f1, "Segment Chi2/ndof ME+1/2", 1110, "Segments_hSChiSq+12.png");
0417   make1DPlot("Segments/hSChiSq+13", f1, "Segment Chi2/ndof ME+1/3", 1110, "Segments_hSChiSq+13.png");
0418   make1DPlot("Segments/hSChiSq+21", f1, "Segment Chi2/ndof ME+2/1", 1110, "Segments_hSChiSq+21.png");
0419   make1DPlot("Segments/hSChiSq+22", f1, "Segment Chi2/ndof ME+2/2", 1110, "Segments_hSChiSq+22.png");
0420   make1DPlot("Segments/hSChiSq+31", f1, "Segment Chi2/ndof ME+3/1", 1110, "Segments_hSChiSq+31.png");
0421   make1DPlot("Segments/hSChiSq+32", f1, "Segment Chi2/ndof ME+3/2", 1110, "Segments_hSChiSq+32.png");
0422   make1DPlot("Segments/hSChiSq+41", f1, "Segment Chi2/ndof ME+4/1", 1110, "Segments_hSChiSq+41.png");
0423   make1DPlot("Segments/hSChiSq-11", f1, "Segment Chi2/ndof ME-1/1b", 1110, "Segments_hSChiSq-11.png");
0424   make1DPlot("Segments/hSChiSq-14", f1, "Segment Chi2/ndof ME-1/1a", 1110, "Segments_hSChiSq-11a.png");
0425   make1DPlot("Segments/hSChiSq-12", f1, "Segment Chi2/ndof ME-1/2", 1110, "Segments_hSChiSq-12.png");
0426   make1DPlot("Segments/hSChiSq-13", f1, "Segment Chi2/ndof ME-1/3", 1110, "Segments_hSChiSq-13.png");
0427   make1DPlot("Segments/hSChiSq-21", f1, "Segment Chi2/ndof ME-2/1", 1110, "Segments_hSChiSq-21.png");
0428   make1DPlot("Segments/hSChiSq-22", f1, "Segment Chi2/ndof ME-2/2", 1110, "Segments_hSChiSq-22.png");
0429   make1DPlot("Segments/hSChiSq-31", f1, "Segment Chi2/ndof ME-3/1", 1110, "Segments_hSChiSq-31.png");
0430   make1DPlot("Segments/hSChiSq-32", f1, "Segment Chi2/ndof ME-3/2", 1110, "Segments_hSChiSq-32.png");
0431   make1DPlot("Segments/hSChiSq-41", f1, "Segment Chi2/ndof ME-4/1", 1110, "Segments_hSChiSq-41.png");
0432 
0433   //miscellaneous
0434   make1DPlot("Segments/hSGlobalPhi", f1, "Segment Global Phi", 1110, "Segments_hSGlobalPhi.png");
0435   make1DPlot("Segments/hSGlobalTheta", f1, "Segment Global Theta", 1110, "Segments_hSGlobalTheta.png");
0436 }