Back to home page

Project CMSSW displayed by LXR

 
 

    


File indexing completed on 2024-04-06 11:59:22

0001 //////////////////////////////////////////////////////////
0002 // This class has been automatically generated on
0003 // Wed Aug 11 09:24:19 2010 by ROOT version 5.22/00d
0004 // from TTree tree/tree
0005 //////////////////////////////////////////////////////////
0006 
0007 #ifndef TreeAnalysisHcalNZS_h
0008 #define TreeAnalysisHcalNZS_h
0009 
0010 #include <cstdlib>
0011 #include <iostream>
0012 #include <fstream>
0013 #include <vector>
0014 #include <string>
0015 #include <cmath>
0016 #include <iomanip>
0017 
0018 #include <TROOT.h>
0019 #include <TChain.h>
0020 #include <TFile.h>
0021 #include "TDirectory.h"
0022 #include "TH1F.h"
0023 #include "TH2.h"
0024 #include "TProfile.h"
0025 #include "TChain.h"
0026 #include "TString.h"
0027 
0028 class TreeAnalysisHcalNZS {
0029 
0030 public :
0031 
0032   std::string ecalCharIso;
0033   std::string hcalCharIso;
0034   std::string dataType;
0035   std::string L1Seed;
0036   double      ebNeutIso, eeNeutIso, hhNeutIso;
0037   double      dRL1Jet;
0038   int         GoodPVCut;
0039 
0040   //declaration of histograms
0041   static const int NEtaBins = 12;
0042   static const int NPBins   = 15; 
0043 
0044   double genPartPBins[NPBins+1], genPartEtaBins[NEtaBins+1];
0045   TFile *fout;
0046 
0047   TProfile *h_HcalMeanEneVsEta;
0048 
0049   TProfile *h_trackPCaloE11x11H3x3[NEtaBins+1], *h_trackPCaloE9x9H3x3[NEtaBins+1], *h_trackPCaloE7x7H3x3[NEtaBins+1];
0050 
0051   TH1F *h_NPV[2], *h_PVx[2], *h_PVy[2], *h_PVr[2], *h_PVz[2], *h_PVNDOF[2];
0052   TH2F *h_PVNTracksSumPt[2];
0053   TH1F *h_PVNTracks[2], *h_PVTracksSumPt[2], *h_PVNTracksWt[2], *h_PVTracksSumPtWt[2];
0054 
0055   TH1F *h_trackP[6],  *h_trackEta[6], *h_trackPhi[6];
0056   TH1F *h_trackPt[6], *h_trackDxy[6], *h_trackDz[6], *h_trackChiSq[6];
0057 
0058   TH1F *h_trackPhi_2_2[NEtaBins];
0059   TH1F *h_trackPhi_3_3[NEtaBins];
0060 
0061   TH1F *h_trackPhi_3_Inner[NEtaBins];
0062   TH1F *h_trackPhi_3_Outer[NEtaBins];
0063 
0064   TH2F *h_eECAL11x11VsHCAL3x3[NEtaBins];
0065 
0066   TH1F *h_meanTrackP[NPBins][NEtaBins];
0067 
0068   TH1F *h_drTrackL1[NPBins][NEtaBins];
0069     
0070   TH1F *h_eECAL7x7_Frac[NPBins][NEtaBins],
0071        *h_eECAL9x9_Frac[NPBins][NEtaBins],
0072        *h_eECAL11x11_Frac[NPBins][NEtaBins],
0073        *h_eECAL15x15_Frac[NPBins][NEtaBins],
0074        *h_eECAL21x21_Frac[NPBins][NEtaBins],
0075        *h_eECAL25x25_Frac[NPBins][NEtaBins],
0076        *h_eECAL31x31_Frac[NPBins][NEtaBins];
0077 
0078   TH1F *hh_eECAL7x7_Frac[NEtaBins],
0079        *hh_eECAL9x9_Frac[NEtaBins],
0080        *hh_eECAL11x11_Frac[NEtaBins],
0081        *hh_eECAL15x15_Frac[NEtaBins],
0082        *hh_eECAL21x21_Frac[NEtaBins],
0083        *hh_eECAL25x25_Frac[NEtaBins],
0084        *hh_eECAL31x31_Frac[NEtaBins];
0085 
0086   TH1F *h_eHCAL3x3_Frac[NPBins][NEtaBins],
0087        *h_eHCAL5x5_Frac[NPBins][NEtaBins],
0088        *h_eHCAL7x7_Frac[NPBins][NEtaBins];
0089   TH1F *h_eHCAL3x3_Frac_20Sig[NPBins][NEtaBins],
0090        *h_eHCAL5x5_Frac_20Sig[NPBins][NEtaBins],
0091        *h_eHCAL7x7_Frac_20Sig[NPBins][NEtaBins];
0092   TH1F *h_eHCAL3x3MIP_Frac[NPBins][NEtaBins],
0093        *h_eHCAL5x5MIP_Frac[NPBins][NEtaBins],
0094        *h_eHCAL7x7MIP_Frac[NPBins][NEtaBins];
0095   TH1F *h_eHCAL3x3MIP_Frac_20Sig[NPBins][NEtaBins],
0096        *h_eHCAL5x5MIP_Frac_20Sig[NPBins][NEtaBins],
0097        *h_eHCAL7x7MIP_Frac_20Sig[NPBins][NEtaBins];
0098 
0099   TH1F *hh_eHCAL3x3_Frac[NEtaBins],
0100        *hh_eHCAL5x5_Frac[NEtaBins],
0101        *hh_eHCAL7x7_Frac[NEtaBins];
0102   TH1F *hh_eHCAL3x3_Frac_20Sig[NEtaBins],
0103        *hh_eHCAL5x5_Frac_20Sig[NEtaBins],
0104        *hh_eHCAL7x7_Frac_20Sig[NEtaBins];
0105 
0106   TH1F *h_eHCAL3x3_eECAL11x11_response[NPBins][NEtaBins],
0107        *h_eHCAL5x5_eECAL11x11_response[NPBins][NEtaBins],
0108        *h_eHCAL7x7_eECAL11x11_response[NPBins][NEtaBins];
0109   TH1F *h_eHCAL3x3_eECAL11x11_responseMIP[NPBins][NEtaBins],
0110        *h_eHCAL5x5_eECAL11x11_responseMIP[NPBins][NEtaBins],
0111        *h_eHCAL7x7_eECAL11x11_responseMIP[NPBins][NEtaBins];
0112   TH1F *h_eHCAL3x3_eECAL11x11_responseInteract[NPBins][NEtaBins],
0113        *h_eHCAL5x5_eECAL11x11_responseInteract[NPBins][NEtaBins],
0114        *h_eHCAL7x7_eECAL11x11_responseInteract[NPBins][NEtaBins];
0115 
0116   TH1F *h_eHCAL3x3_eECAL9x9_response[NPBins][NEtaBins],
0117        *h_eHCAL5x5_eECAL9x9_response[NPBins][NEtaBins],
0118        *h_eHCAL7x7_eECAL9x9_response[NPBins][NEtaBins];
0119   TH1F *h_eHCAL3x3_eECAL9x9_responseMIP[NPBins][NEtaBins],
0120        *h_eHCAL5x5_eECAL9x9_responseMIP[NPBins][NEtaBins],
0121        *h_eHCAL7x7_eECAL9x9_responseMIP[NPBins][NEtaBins];
0122   TH1F *h_eHCAL3x3_eECAL9x9_responseInteract[NPBins][NEtaBins],
0123        *h_eHCAL5x5_eECAL9x9_responseInteract[NPBins][NEtaBins],
0124        *h_eHCAL7x7_eECAL9x9_responseInteract[NPBins][NEtaBins];
0125 
0126   TH1F *h_eHCAL3x3_eECAL7x7_response[NPBins][NEtaBins],
0127        *h_eHCAL5x5_eECAL7x7_response[NPBins][NEtaBins],
0128        *h_eHCAL7x7_eECAL7x7_response[NPBins][NEtaBins];
0129   TH1F *h_eHCAL3x3_eECAL7x7_responseMIP[NPBins][NEtaBins],
0130        *h_eHCAL5x5_eECAL7x7_responseMIP[NPBins][NEtaBins],
0131        *h_eHCAL7x7_eECAL7x7_responseMIP[NPBins][NEtaBins];
0132   TH1F *h_eHCAL3x3_eECAL7x7_responseInteract[NPBins][NEtaBins],
0133        *h_eHCAL5x5_eECAL7x7_responseInteract[NPBins][NEtaBins],
0134        *h_eHCAL7x7_eECAL7x7_responseInteract[NPBins][NEtaBins];
0135 
0136 
0137   //===
0138   TH1F *h_diff_e15x15e11x11[NPBins][NEtaBins],
0139     *h_diff_e15x15e11x11_20Sig[NPBins][NEtaBins];
0140 
0141   TH1F *h_diff_h7x7h5x5[NPBins][NEtaBins];
0142   TH1F *h_maxNearP31x31[NPBins][NEtaBins];
0143         
0144   TH1F *h_eECAL7x7_Frac_20Sig[NPBins][NEtaBins],
0145        *h_eECAL9x9_Frac_20Sig[NPBins][NEtaBins],
0146        *h_eECAL11x11_Frac_20Sig[NPBins][NEtaBins];
0147   TH1F *h_eHCAL3x3_eECAL11x11_response_20Sig[NPBins][NEtaBins];
0148   TH1F *h_eHCAL3x3_eECAL11x11_responseMIP_20Sig[NPBins][NEtaBins];
0149   TH1F *h_eHCAL3x3_eECAL11x11_responseInteract_20Sig[NPBins][NEtaBins];
0150   TH1F *h_eHCAL3x3_eECAL9x9_response_20Sig[NPBins][NEtaBins];
0151   TH1F *h_eHCAL3x3_eECAL9x9_responseMIP_20Sig[NPBins][NEtaBins];
0152   TH1F *h_eHCAL3x3_eECAL9x9_responseInteract_20Sig[NPBins][NEtaBins];
0153   TH1F *h_eHCAL3x3_eECAL7x7_response_20Sig[NPBins][NEtaBins];
0154   TH1F *h_eHCAL3x3_eECAL7x7_responseMIP_20Sig[NPBins][NEtaBins];
0155   TH1F *h_eHCAL3x3_eECAL7x7_responseInteract_20Sig[NPBins][NEtaBins];
0156   TH1F *h_eHCAL5x5_eECAL11x11_response_20Sig[NPBins][NEtaBins];
0157   TH1F *h_eHCAL5x5_eECAL11x11_responseMIP_20Sig[NPBins][NEtaBins];
0158   TH1F *h_eHCAL5x5_eECAL11x11_responseInteract_20Sig[NPBins][NEtaBins];
0159 
0160   TH1F *hh_eECAL7x7_Frac_20Sig[NEtaBins],
0161        *hh_eECAL9x9_Frac_20Sig[NEtaBins],
0162        *hh_eECAL11x11_Frac_20Sig[NEtaBins];
0163 
0164   //================================
0165    TChain          *fChain;   //!pointer to the analyzed TTree or TChain
0166    Int_t           fCurrent; //!current Tree number in a TChain
0167 
0168    // Declaration of leaf types
0169    Int_t                    t_EvtNo;
0170    Int_t                    t_RunNo;
0171    Int_t                    t_Lumi;
0172    Int_t                    t_Bunch;
0173    std::vector<double>      *PVx;
0174    std::vector<double>      *PVy;
0175    std::vector<double>      *PVz;
0176    std::vector<int>         *PVisValid;
0177    std::vector<int>         *PVndof;
0178    std::vector<int>         *PVNTracks;
0179    std::vector<int>         *PVNTracksWt;
0180    std::vector<double>      *t_PVTracksSumPt;
0181    std::vector<double>      *t_PVTracksSumPtWt;
0182    std::vector<int>         *PVNTracksHP;
0183    std::vector<int>         *PVNTracksHPWt;
0184    std::vector<double>      *t_PVTracksSumPtHP;
0185    std::vector<double>      *t_PVTracksSumPtHPWt;
0186 
0187    Int_t                     t_L1Decision[128];
0188    std::vector<std::string> *t_L1AlgoNames;
0189    std::vector<int>         *t_L1PreScale;
0190    std::vector<double>      *t_L1CenJetPt;
0191    std::vector<double>      *t_L1CenJetEta;
0192    std::vector<double>      *t_L1CenJetPhi;
0193    std::vector<double>      *t_L1FwdJetPt;
0194    std::vector<double>      *t_L1FwdJetEta;
0195    std::vector<double>      *t_L1FwdJetPhi;
0196    std::vector<double>      *t_L1TauJetPt;
0197    std::vector<double>      *t_L1TauJetEta;
0198    std::vector<double>      *t_L1TauJetPhi;
0199    std::vector<double>      *t_L1MuonPt;
0200    std::vector<double>      *t_L1MuonEta;
0201    std::vector<double>      *t_L1MuonPhi;
0202    std::vector<double>      *t_L1IsoEMPt;
0203    std::vector<double>      *t_L1IsoEMEta;
0204    std::vector<double>      *t_L1IsoEMPhi;
0205    std::vector<double>      *t_L1NonIsoEMPt;
0206    std::vector<double>      *t_L1NonIsoEMEta;
0207    std::vector<double>      *t_L1NonIsoEMPhi;
0208    std::vector<double>      *t_L1METPt;
0209    std::vector<double>      *t_L1METEta;
0210    std::vector<double>      *t_L1METPhi;
0211 
0212    std::vector<double>      *t_jetPt;
0213    std::vector<double>      *t_jetEta;
0214    std::vector<double>      *t_jetPhi;
0215    std::vector<double>      *t_nTrksJetCalo;
0216    std::vector<double>      *t_nTrksJetVtx;
0217 
0218    std::vector<double>      *t_trackPAll;
0219    std::vector<double>      *t_trackEtaAll;
0220    std::vector<double>      *t_trackPhiAll;
0221    std::vector<double>      *t_trackPtAll;
0222    std::vector<double>      *t_trackDxyAll;
0223    std::vector<double>      *t_trackDzAll;
0224    std::vector<double>      *t_trackDxyPVAll;
0225    std::vector<double>      *t_trackDzPVAll;
0226    std::vector<double>      *t_trackChiSqAll;
0227 
0228    std::vector<double>      *t_trackP;
0229    std::vector<double>      *t_trackPt;
0230    std::vector<double>      *t_trackEta;
0231    std::vector<double>      *t_trackPhi;
0232    std::vector<double>      *t_trackEcalEta;
0233    std::vector<double>      *t_trackEcalPhi;
0234    std::vector<double>      *t_trackHcalEta;
0235    std::vector<double>      *t_trackHcalPhi;
0236    std::vector<int>         *t_trackNOuterHits;
0237    std::vector<int>         *t_NLayersCrossed;
0238    std::vector<int>         *t_trackHitsTOB;
0239    std::vector<int>         *t_trackHitsTEC;
0240    std::vector<int>         *t_trackHitInMissTOB;
0241    std::vector<int>         *t_trackHitInMissTEC;
0242    std::vector<int>         *t_trackHitInMissTIB;
0243    std::vector<int>         *t_trackHitInMissTID;
0244    std::vector<int>         *t_trackHitInMissTIBTID;
0245    std::vector<int>         *t_trackHitOutMissTOB;
0246    std::vector<int>         *t_trackHitOutMissTEC;
0247    std::vector<int>         *t_trackHitOutMissTIB;
0248    std::vector<int>         *t_trackHitOutMissTID;
0249    std::vector<int>         *t_trackHitOutMissTOBTEC;
0250    std::vector<int>         *t_trackHitInMeasTOB;
0251    std::vector<int>         *t_trackHitInMeasTEC;
0252    std::vector<int>         *t_trackHitInMeasTIB;
0253    std::vector<int>         *t_trackHitInMeasTID;
0254    std::vector<int>         *t_trackHitOutMeasTOB;
0255    std::vector<int>         *t_trackHitOutMeasTEC;
0256    std::vector<int>         *t_trackHitOutMeasTIB;
0257    std::vector<int>         *t_trackHitOutMeasTID;
0258    std::vector<double>      *t_trackOutPosOutHitDr;
0259    std::vector<double>      *t_trackL;
0260    std::vector<double>      *t_trackDxy;
0261    std::vector<double>      *t_trackDxyBS;
0262    std::vector<double>      *t_trackDz;
0263    std::vector<double>      *t_trackDzBS;
0264    std::vector<double>      *t_trackDxyPV;
0265    std::vector<double>      *t_trackDzPV;
0266    std::vector<double>      *t_trackChiSq;
0267    std::vector<int>         *t_trackPVIdx;
0268 
0269    std::vector<double>      *t_maxNearP31x31;
0270    std::vector<double>      *t_maxNearP21x21;
0271    std::vector<int>         *t_ecalSpike11x11;
0272    std::vector<double>      *t_e7x7;
0273    std::vector<double>      *t_e9x9;
0274    std::vector<double>      *t_e11x11;
0275    std::vector<double>      *t_e15x15;
0276    std::vector<double>      *t_e7x7_10Sig;
0277    std::vector<double>      *t_e9x9_10Sig;
0278    std::vector<double>      *t_e11x11_10Sig;
0279    std::vector<double>      *t_e15x15_10Sig;
0280    std::vector<double>      *t_e7x7_15Sig;
0281    std::vector<double>      *t_e9x9_15Sig;
0282    std::vector<double>      *t_e11x11_15Sig;
0283    std::vector<double>      *t_e15x15_15Sig;
0284    std::vector<double>      *t_e7x7_20Sig;
0285    std::vector<double>      *t_e9x9_20Sig;
0286    std::vector<double>      *t_e11x11_20Sig;
0287    std::vector<double>      *t_e15x15_20Sig;
0288    std::vector<double>      *t_e7x7_25Sig;
0289    std::vector<double>      *t_e9x9_25Sig;
0290    std::vector<double>      *t_e11x11_25Sig;
0291    std::vector<double>      *t_e15x15_25Sig;
0292    std::vector<double>      *t_e7x7_30Sig;
0293    std::vector<double>      *t_e9x9_30Sig;
0294    std::vector<double>      *t_e11x11_30Sig;
0295    std::vector<double>      *t_e15x15_30Sig;
0296 
0297    std::vector<double>      *t_maxNearHcalP3x3;
0298    std::vector<double>      *t_maxNearHcalP5x5;
0299    std::vector<double>      *t_maxNearHcalP7x7;
0300    std::vector<double>      *t_h3x3;
0301    std::vector<double>      *t_h5x5;
0302    std::vector<double>      *t_h7x7;
0303    std::vector<double>      *t_h3x3Sig;
0304    std::vector<double>      *t_h5x5Sig;
0305    std::vector<double>      *t_h7x7Sig;
0306    std::vector<int>         *t_infoHcal;
0307    Int_t                     t_nTracks;
0308 
0309    // List of branches
0310    TBranch                  *b_t_EvtNo;
0311    TBranch                  *b_t_RunNo;
0312    TBranch                  *b_t_Lumi; 
0313    TBranch                  *b_t_Bunch;
0314    TBranch                  *b_PVx;
0315    TBranch                  *b_PVy;
0316    TBranch                  *b_PVz;
0317    TBranch                  *b_PVisValid;
0318    TBranch                  *b_PVndof; 
0319    TBranch                  *b_PVNTracks;
0320    TBranch                  *b_PVNTracksWt;
0321    TBranch                  *b_t_PVTracksSumPt;
0322    TBranch                  *b_t_PVTracksSumPtWt;
0323    TBranch                  *b_PVNTracksHP;
0324    TBranch                  *b_PVNTracksHPWt;
0325    TBranch                  *b_t_PVTracksSumPtHP;
0326    TBranch                  *b_t_PVTracksSumPtHPWt;
0327 
0328    TBranch                  *b_t_L1Decision;
0329    TBranch                  *b_t_L1AlgoNames;
0330    TBranch                  *b_t_L1PreScale;
0331    TBranch                  *b_t_L1CenJetPt;
0332    TBranch                  *b_t_L1CenJetEta;
0333    TBranch                  *b_t_L1CenJetPhi;
0334    TBranch                  *b_t_L1FwdJetPt;
0335    TBranch                  *b_t_L1FwdJetEta;
0336    TBranch                  *b_t_L1FwdJetPhi;
0337    TBranch                  *b_t_L1TauJetPt;
0338    TBranch                  *b_t_L1TauJetEta;
0339    TBranch                  *b_t_L1TauJetPhi;
0340    TBranch                  *b_t_L1MuonPt;
0341    TBranch                  *b_t_L1MuonEta;
0342    TBranch                  *b_t_L1MuonPhi;
0343    TBranch                  *b_t_L1IsoEMPt;
0344    TBranch                  *b_t_L1IsoEMEta;
0345    TBranch                  *b_t_L1IsoEMPhi;
0346    TBranch                  *b_t_L1NonIsoEMPt;
0347    TBranch                  *b_t_L1NonIsoEMEta;
0348    TBranch                  *b_t_L1NonIsoEMPhi;
0349    TBranch                  *b_t_L1METPt;
0350    TBranch                  *b_t_L1METEta;
0351    TBranch                  *b_t_L1METPhi;
0352 
0353    TBranch                  *b_t_jetPt;
0354    TBranch                  *b_t_jetEta;
0355    TBranch                  *b_t_jetPhi;
0356    TBranch                  *b_t_nTrksJetCalo;
0357    TBranch                  *b_t_nTrksJetVtx;
0358 
0359    TBranch                  *b_t_trackPAll;
0360    TBranch                  *b_t_trackEtaAll;
0361    TBranch                  *b_t_trackPhiAll;
0362    TBranch                  *b_t_trackPtAll;
0363    TBranch                  *b_t_trackDxyAll;
0364    TBranch                  *b_t_trackDzAll;
0365    TBranch                  *b_t_trackDxyPVAll;
0366    TBranch                  *b_t_trackDzPVAll;
0367    TBranch                  *b_t_trackChiSqAll;
0368 
0369    TBranch                  *b_t_trackP;
0370    TBranch                  *b_t_trackPt;
0371    TBranch                  *b_t_trackEta;
0372    TBranch                  *b_t_trackPhi;
0373    TBranch                  *b_t_trackEcalEta;
0374    TBranch                  *b_t_trackEcalPhi;
0375    TBranch                  *b_t_trackHcalEta;
0376    TBranch                  *b_t_trackHcalPhi;
0377    TBranch                  *b_t_trackNOuterHits;
0378    TBranch                  *b_t_NLayersCrossed;
0379    TBranch                  *b_t_trackHitsTOB;
0380    TBranch                  *b_t_trackHitsTEC;
0381    TBranch                  *b_t_trackHitInMissTOB;
0382    TBranch                  *b_t_trackHitInMissTEC;
0383    TBranch                  *b_t_trackHitInMissTIB;
0384    TBranch                  *b_t_trackHitInMissTID;
0385    TBranch                  *b_t_trackHitInMissTIBTID;
0386    TBranch                  *b_t_trackHitOutMissTOB;
0387    TBranch                  *b_t_trackHitOutMissTEC;
0388    TBranch                  *b_t_trackHitOutMissTIB;
0389    TBranch                  *b_t_trackHitOutMissTID;
0390    TBranch                  *b_t_trackHitOutMissTOBTEC;
0391    TBranch                  *b_t_trackHitInMeasTOB;
0392    TBranch                  *b_t_trackHitInMeasTEC;
0393    TBranch                  *b_t_trackHitInMeasTIB;
0394    TBranch                  *b_t_trackHitInMeasTID;
0395    TBranch                  *b_t_trackHitOutMeasTOB;
0396    TBranch                  *b_t_trackHitOutMeasTEC;
0397    TBranch                  *b_t_trackHitOutMeasTIB;
0398    TBranch                  *b_t_trackHitOutMeasTID;
0399    TBranch                  *b_t_trackOutPosOutHitDr;
0400    TBranch                  *b_t_trackL;
0401    TBranch                  *b_t_trackDxy;
0402    TBranch                  *b_t_trackDxyBS;
0403    TBranch                  *b_t_trackDz;
0404    TBranch                  *b_t_trackDzBS;
0405    TBranch                  *b_t_trackDxyPV;
0406    TBranch                  *b_t_trackDzPV; 
0407    TBranch                  *b_t_trackChiSq;
0408    TBranch                  *b_t_trackPVIdx;
0409 
0410    TBranch                  *b_t_maxNearP31x31;
0411    TBranch                  *b_t_maxNearP21x21;
0412    TBranch                  *b_t_ecalSpike11x11;
0413    TBranch                  *b_t_e7x7;
0414    TBranch                  *b_t_e9x9;
0415    TBranch                  *b_t_e11x11;
0416    TBranch                  *b_t_e15x15;
0417    TBranch                  *b_t_e7x7_10Sig;
0418    TBranch                  *b_t_e9x9_10Sig;
0419    TBranch                  *b_t_e11x11_10Sig;
0420    TBranch                  *b_t_e15x15_10Sig;
0421    TBranch                  *b_t_e7x7_15Sig;
0422    TBranch                  *b_t_e9x9_15Sig;
0423    TBranch                  *b_t_e11x11_15Sig;
0424    TBranch                  *b_t_e15x15_15Sig;
0425    TBranch                  *b_t_e7x7_20Sig;
0426    TBranch                  *b_t_e9x9_20Sig;
0427    TBranch                  *b_t_e11x11_20Sig;
0428    TBranch                  *b_t_e15x15_20Sig;
0429    TBranch                  *b_t_e7x7_25Sig;
0430    TBranch                  *b_t_e9x9_25Sig;
0431    TBranch                  *b_t_e11x11_25Sig;
0432    TBranch                  *b_t_e15x15_25Sig;
0433    TBranch                  *b_t_e7x7_30Sig;
0434    TBranch                  *b_t_e9x9_30Sig;
0435    TBranch                  *b_t_e11x11_30Sig;
0436    TBranch                  *b_t_e15x15_30Sig;
0437 
0438    TBranch                  *b_t_maxNearHcalP3x3;
0439    TBranch                  *b_t_maxNearHcalP5x5;
0440    TBranch                  *b_t_maxNearHcalP7x7;
0441    TBranch                  *b_t_h3x3;
0442    TBranch                  *b_t_h5x5;
0443    TBranch                  *b_t_h7x7;
0444    TBranch                  *b_t_h3x3Sig;
0445    TBranch                  *b_t_h5x5Sig;
0446    TBranch                  *b_t_h7x7Sig;
0447    TBranch                  *b_t_infoHcal;
0448    TBranch                  *b_t_nTracks; 
0449 
0450    TreeAnalysisHcalNZS(TChain *tree, const char *outFileName);
0451    virtual ~TreeAnalysisHcalNZS();
0452 
0453    virtual Int_t    Cut(Long64_t entry);
0454    virtual Int_t    GetEntry(Long64_t entry);
0455    virtual Long64_t LoadTree(Long64_t entry);
0456    virtual void     Init(TChain *tree);
0457    virtual void     Loop(int cut=1);
0458    virtual Bool_t   Notify();
0459    virtual void     Show(Long64_t entry = -1);
0460 
0461    double DeltaPhi(double v1, double v2);
0462    double DeltaR(double eta1, double phi1, double eta2, double phi2);
0463    void   BookHistograms(const char *outFileName);
0464 };
0465 
0466 #endif