CommonData/interface/ElectronTree.h

#ifndef ElectronTree_H
#define ElectronTree_H

#include "TFile.h"
#include "TTree.h"
#include "TError.h"
#include <cmath>
#include "assert.h"

namespace citana
{
  class ElectronTree {

    public:

      /// bit map
      /// DON'T CHANGE ORDER

      //*******************************************
      //=== ElectronTriggerBits  ====
      //*******************************************
      enum ElectronTriggerBits { kEleTrigger_Ele                                    = 0x000001, 
                                 kEleTrigger_Ele_CaloIdL_CaloIsoVL                  = 0x000002, 
                                 kEleTrigger_Ele_CaloIdT_TrkIdVL_CaloIsoVL_TrkIsoVL = 0x000004 
      };

      /// variables
      Float_t                 fWeight;
      UInt_t                  fRunNumber;
      UInt_t                  fLumiSectionNumber;
      UInt_t                  fEventNumber;
      Bool_t                  fEleEventNumberParity;
      Float_t                 fElePt; 
      Float_t                 fEleEta; 
      Float_t                 fElePhi; 
      Float_t                 fEleSCEt; 
      Float_t                 fEleSCEta; 
      Float_t                 fEleSCPhi; 
      Float_t                 fEleEcalEnergy; 
      Bool_t                  fEleIsEcalDriven;
      UInt_t                  fEleTriggerBit;
      Float_t                 fRho; 
      UInt_t                  fNVertices; 

      // Conversion and IP
      Float_t                 fEleD0; 
      Float_t                 fEleDZ; 
      Float_t                 fEleIP3d; 
      Float_t                 fEleIP3dSig; 
      Bool_t                  fEleMatchedConversion;
      Float_t                 fEleConvDCot;
      Float_t                 fEleConvDist;
      UInt_t                  fEleNMissHits;

      // E/P variables
      Float_t                 fEleNBrem; 
      Float_t                 fEleFBrem; 
      Float_t                 fEleEOverP; 
      Float_t                 fEleESeedClusterOverPIn; 
      Float_t                 fEleESeedClusterOverPout; 
      Float_t                 fEleEEleClusterOverPout; 
      Float_t                 fEleOneOverEMinusOneOverP; 

      // track cluster matching
      Float_t                 fEleDEtaIn; 
      Float_t                 fEleDPhiIn; 
      Float_t                 fEledEtaCalo;
      Float_t                 fEledPhiCalo;

      //shower shape
      Float_t                 fEleSigmaIEtaIEta; 
      Float_t                 fEleSigmaIPhiIPhi; 
      Float_t                 fEleSigmaIEtaIPhi;
      Float_t                 fEleSCEtaWidth;
      Float_t                 fEleSCPhiWidth;
      Float_t                 fEleR9;
      Float_t                 fElePreShowerOverRaw;
      Float_t                 fEleHoverE; 

      //track quality
      Float_t                 fEleGsfTrackChi2OverNdof;    
      Float_t                 fEleKFTrackChi2OverNDoF;
      Int_t                   fEleKFTrackNHits;
      Int_t                   fEleKFTrackNLayersWithMeasurement;
      Float_t                 fEleOneMinusSeedE1x5OverE5x5;

      //Isolation Variables
      Float_t                 fElePFMVA;
      Float_t                 fEleTrkIso03; 
      Float_t                 fEleEMIso03; 
      Float_t                 fEleHadIso03; 
      Float_t                 fEleTrkIso04; 
      Float_t                 fEleEMIso04; 
      Float_t                 fEleHadIso04;     
      Float_t                 fElePFIso04; 
      Float_t                 fChargedIso_DR0p0To0p1;
      Float_t                 fChargedIso_DR0p1To0p2;
      Float_t                 fChargedIso_DR0p2To0p3;
      Float_t                 fChargedIso_DR0p3To0p4;
      Float_t                 fChargedIso_DR0p4To0p5;
      Float_t                 fGammaIso_DR0p0To0p1;
      Float_t                 fGammaIso_DR0p1To0p2;
      Float_t                 fGammaIso_DR0p2To0p3;
      Float_t                 fGammaIso_DR0p3To0p4;
      Float_t                 fGammaIso_DR0p4To0p5;
      Float_t                 fNeutralHadronIso_DR0p0To0p1;
      Float_t                 fNeutralHadronIso_DR0p1To0p2;
      Float_t                 fNeutralHadronIso_DR0p2To0p3;
      Float_t                 fNeutralHadronIso_DR0p3To0p4;
      Float_t                 fNeutralHadronIso_DR0p4To0p5;

      Bool_t                  fElePassTriggerDenominator;

      //Regression Variables
      Bool_t                  fIsEB;           
      Bool_t                  fIsEE;           
      Float_t                 fSCRawEnergy;
      Float_t                 fE5x5;
      Float_t                 fNClusters;
      Float_t                 fEtaSeed;
      Float_t                 fPhiSeed;
      Float_t                 fESeed;
      Float_t                 fE3x3Seed;
      Float_t                 fE5x5Seed;
      Float_t                 fEMaxSeed;
      Float_t                 fE2ndSeed;
      Float_t                 fETopSeed;
      Float_t                 fEBottomSeed;
      Float_t                 fELeftSeed;
      Float_t                 fERightSeed;
      Float_t                 fE2x5MaxSeed;
      Float_t                 fE2x5TopSeed;
      Float_t                 fE2x5BottomSeed;
      Float_t                 fE2x5LeftSeed;
      Float_t                 fE2x5RightSeed;
      Float_t                 fIEtaSeed;
      Float_t                 fIPhiSeed;
      Float_t                 fEtaCrySeed;
      Float_t                 fPhiCrySeed;
      Float_t                 fGsfTrackPIn;
      Int_t                   fCharge;
      Float_t                 fGeneratedEnergy;
      Float_t                 fGeneratedEnergyStatus1;
      Float_t                 fGeneratedEnergyStatus3;

    public:
      /// this is the main element
      TTree *tree_;
      TFile *f_;
  
      /// hold the names of variables to facilitate things (filled during Init)
      std::vector<std::string> variables_;

      /// default constructor  
      ElectronTree()  {};
      /// default destructor
      ~ElectronTree(){ 
        if (f_) f_->Close();  
      };
    
      /// initialize varibles and fill list of available variables
      void InitVariables() {
        fWeight                        = 0.0;
        fRunNumber                     = 0.0;
        fLumiSectionNumber             = 0.0;
        fEventNumber                   = 0.0;
        fEleEventNumberParity          = 0.0;
        fElePt                         = 0.0;
        fEleEta                        = 0.0;
        fElePhi                        = 0.0;
        fEleSCEt                       = 0.0;
        fEleSCEta                      = 0.0;
        fEleSCPhi                      = 0.0;
        fEleEcalEnergy                 = 0.0;
        fEleIsEcalDriven               = 0.0;
        fEleTriggerBit                 = 0.0;
        fRho                           = 0.0;
        fNVertices                     = 0.0;
        fEleD0                         = 0.0;
        fEleDZ                         = 0.0;
        fEleIP3d                       = 0.0;
        fEleIP3dSig                    = 0.0;
        fEleMatchedConversion          = 0.0;
        fEleConvDCot                   = 0.0;
        fEleConvDist                   = 0.0;
        fEleNMissHits                  = 0.0;
        fEleNBrem                      = 0.0;
        fEleFBrem                      = 0.0;
        fEleEOverP                     = 0.0;
        fEleESeedClusterOverPIn        = 0.0;
        fEleESeedClusterOverPout       = 0.0;
        fEleEEleClusterOverPout        = 0.0;
        fEleOneOverEMinusOneOverP      = 0.0;
        fEleDEtaIn                     = 0.0;
        fEleDPhiIn                     = 0.0;
        fEledEtaCalo                   = 0.0;
        fEledPhiCalo                   = 0.0;
        fEleSigmaIEtaIEta              = 0.0;
        fEleSigmaIPhiIPhi              = 0.0;
        fEleSigmaIEtaIPhi              = 0.0;
        fEleSCEtaWidth                 = 0.0;
        fEleSCPhiWidth                 = 0.0;
        fEleR9                         = 0.0;
        fElePreShowerOverRaw           = 0.0;
        fEleHoverE                     = 0.0;
        fEleGsfTrackChi2OverNdof       = 0.0;
        fEleKFTrackChi2OverNDoF        = 0.0;
        fEleKFTrackNHits               = 0.0;
        fEleKFTrackNLayersWithMeasurement = 0.0;
        fEleOneMinusSeedE1x5OverE5x5   = 0.0;
        fElePFMVA                      = 0.0;
        fEleTrkIso03                   = 0.0;
        fEleEMIso03                    = 0.0;
        fEleHadIso03                   = 0.0;
        fEleTrkIso04                   = 0.0;
        fEleEMIso04                    = 0.0;
        fEleHadIso04                   = 0.0;
        fElePFIso04                    = 0.0;
        fChargedIso_DR0p0To0p1         = 0.0;
        fChargedIso_DR0p1To0p2         = 0.0;
        fChargedIso_DR0p2To0p3         = 0.0;
        fChargedIso_DR0p3To0p4         = 0.0;
        fChargedIso_DR0p4To0p5         = 0.0;
        fGammaIso_DR0p0To0p1           = 0.0;
        fGammaIso_DR0p1To0p2           = 0.0;
        fGammaIso_DR0p2To0p3           = 0.0;
        fGammaIso_DR0p3To0p4           = 0.0;
        fGammaIso_DR0p4To0p5           = 0.0;
        fNeutralHadronIso_DR0p0To0p1   = 0.0;
        fNeutralHadronIso_DR0p1To0p2   = 0.0;
        fNeutralHadronIso_DR0p2To0p3   = 0.0;
        fNeutralHadronIso_DR0p3To0p4   = 0.0;
        fNeutralHadronIso_DR0p4To0p5   = 0.0;
        fElePassTriggerDenominator     = 0.0;
        fIsEB                          = 0.0;           
        fIsEE                          = 0.0;           
        fSCRawEnergy                   = 0.0;
        fE5x5                          = 0.0;
        fNClusters                     = 0.0;
        fEtaSeed                       = 0.0;
        fPhiSeed                       = 0.0;
        fESeed                         = 0.0;
        fE3x3Seed                      = 0.0;
        fE5x5Seed                      = 0.0;
        fEMaxSeed                      = 0.0;
        fE2ndSeed                      = 0.0;
        fETopSeed                      = 0.0;
        fEBottomSeed                   = 0.0;
        fELeftSeed                     = 0.0;
        fERightSeed                    = 0.0;
        fE2x5MaxSeed                   = 0.0;
        fE2x5TopSeed                   = 0.0;
        fE2x5BottomSeed                = 0.0;
        fE2x5LeftSeed                  = 0.0;
        fE2x5RightSeed                 = 0.0;
        fIEtaSeed                      = 0.0;
        fIPhiSeed                      = 0.0;
        fEtaCrySeed                    = 0.0;
        fPhiCrySeed                    = 0.0;
        fGsfTrackPIn                   = 0.0;
        fCharge                        = 0.0;
        fGeneratedEnergy               = 0.0;
        fGeneratedEnergyStatus1        = 0.0;
        fGeneratedEnergyStatus3        = 0.0;


      }
    
      /// load a ElectronTree
      void LoadTree(const char* file){
        f_ = TFile::Open(file);
        assert(f_);
        tree_ = dynamic_cast<TTree*>(f_->Get("Electrons"));
        assert(tree_);
      }
    
      /// create a ElectronTree
      void CreateTree(){
        tree_ = new TTree("Electrons","Electrons");
        f_ = 0;

        //book the branches
        tree_->Branch("weight",&fWeight,"weight/F");
        tree_->Branch("run",&fRunNumber,"run/i");
        tree_->Branch("lumi",&fLumiSectionNumber,"lumi/i");
        tree_->Branch("event",&fEventNumber,"event/i");
        tree_->Branch("EventNumberParity",&fEleEventNumberParity,"EventNumberParity/O"); 
        tree_->Branch("pt",&fElePt,"pt/F"); 
        tree_->Branch("eta",&fEleEta,"eta/F"); 
        tree_->Branch("phi",&fElePhi,"phi/F"); 
        tree_->Branch("scEt",&fEleSCEt,"scEt/F"); 
        tree_->Branch("scEta",&fEleSCEta,"scEta/F"); 
        tree_->Branch("scPhi",&fEleSCPhi,"scPhi/F"); 
        tree_->Branch("ecalenergy",&fEleEcalEnergy,"ecalenergy/F"); 
        tree_->Branch("ecaldriven",&fEleIsEcalDriven,"ecaldriven/O"); 
        tree_->Branch("triggerBit",&fEleTriggerBit,"triggerBit/i"); 
        tree_->Branch("rho",&fRho,"rho/F"); 
        tree_->Branch("vertices",&fNVertices,"vertices/i"); 
        tree_->Branch("d0",&fEleD0,"d0/F"); 
        tree_->Branch("dz",&fEleDZ,"dz/F"); 
        tree_->Branch("ip3d",&fEleIP3d,"ip3d/F"); 
        tree_->Branch("ip3ds",&fEleIP3dSig,"ip3ds/F"); 
        tree_->Branch("matchConv",&fEleMatchedConversion,"matchConv/O"); 
        tree_->Branch("dcot",&fEleConvDCot,"dcot/F"); 
        tree_->Branch("dist",&fEleConvDist,"dist/F"); 
        tree_->Branch("missHits",&fEleNMissHits,"NMissHits/i"); 
        tree_->Branch("nbrems",&fEleNBrem,"nbrems/F"); 
        tree_->Branch("fbrem",&fEleFBrem,"fbrem/F"); 
        tree_->Branch("EoP",&fEleEOverP,"EoP/F"); 
        tree_->Branch("EoPin",&fEleESeedClusterOverPIn,"EoPin/F"); 
        tree_->Branch("ESeedoPout",&fEleESeedClusterOverPout,"ESeedoPout/F"); 
        tree_->Branch("EEleoPout",&fEleEEleClusterOverPout,"EEleoPout/F"); 
//       tree_->Branch("eleEoPout",&fEleEEleClusterOverPout,"eleEoPout/F"); //old name of the branch
        tree_->Branch("IoEmIoP",&fEleOneOverEMinusOneOverP,"IoEmIoP/F");  
        tree_->Branch("deta",&fEleDEtaIn,"deta/F"); 
        tree_->Branch("dphi",&fEleDPhiIn,"dphi/F"); 
        tree_->Branch("detacalo",&fEledEtaCalo,"detacalo/F"); 
        tree_->Branch("dphicalo",&fEledPhiCalo,"dphicalo/F"); 
        tree_->Branch("see",&fEleSigmaIEtaIEta,"see/F"); 
        tree_->Branch("spp",&fEleSigmaIPhiIPhi,"spp/F"); 
        tree_->Branch("sep",&fEleSigmaIEtaIPhi,"sep/F"); 
        tree_->Branch("etawidth",&fEleSCEtaWidth,"etawidth/F"); 
        tree_->Branch("phiwidth",&fEleSCPhiWidth,"phiwidth/F"); 
        tree_->Branch("R9",&fEleR9,"R9/F"); 
        tree_->Branch("PreShowerOverRaw",&fElePreShowerOverRaw,"PreShowerOverRaw/F"); 
        tree_->Branch("HoE",&fEleHoverE,"HoE/F"); 
        tree_->Branch("gsfchi2",&fEleGsfTrackChi2OverNdof,"gsfchi2/F"); 
        tree_->Branch("kfchi2",&fEleKFTrackChi2OverNDoF,"kfchi2/F"); 
        tree_->Branch("kfhits",&fEleKFTrackNHits,"kfhits/I"); 
        tree_->Branch("kflayers",&fEleKFTrackNLayersWithMeasurement,"kflayers/I"); 
        tree_->Branch("OneMinusSeedE1x5OverE5x5",&fEleOneMinusSeedE1x5OverE5x5,"OneMinusSeedE1x5OverE5x5/F"); 
        tree_->Branch("PFMVA",&fElePFMVA,"PFMVA/F"); 
        tree_->Branch("trkIso03",&fEleTrkIso03,"trkIso03/F"); 
        tree_->Branch("ecalIso03",&fEleEMIso03,"ecalIso03/F"); 
        tree_->Branch("hcalIso03",&fEleHadIso03,"hcalIso03/F"); 
        tree_->Branch("trkIso04",&fEleTrkIso04,"trkIso04/F"); 
        tree_->Branch("ecalIso04",&fEleEMIso04,"ecalIso04/F"); 
        tree_->Branch("hcalIso04",&fEleHadIso04,"hcalIso04/F"); 
        tree_->Branch("pfIso04",&fElePFIso04,"pfIso04/F"); 
        tree_->Branch("ChargedIso_DR0p0To0p1",&fChargedIso_DR0p0To0p1,"ChargedIso_DR0p0To0p1/F");
        tree_->Branch("ChargedIso_DR0p1To0p2",&fChargedIso_DR0p1To0p2,"ChargedIso_DR0p1To0p2/F");
        tree_->Branch("ChargedIso_DR0p2To0p3",&fChargedIso_DR0p2To0p3,"ChargedIso_DR0p2To0p3/F");
        tree_->Branch("ChargedIso_DR0p3To0p4",&fChargedIso_DR0p3To0p4,"ChargedIso_DR0p3To0p4/F");
        tree_->Branch("ChargedIso_DR0p4To0p5",&fChargedIso_DR0p4To0p5,"ChargedIso_DR0p4To0p5/F");
        tree_->Branch("GammaIso_DR0p0To0p1",&fGammaIso_DR0p0To0p1,"GammaIso_DR0p0To0p1/F");
        tree_->Branch("GammaIso_DR0p1To0p2",&fGammaIso_DR0p1To0p2,"GammaIso_DR0p1To0p2/F");
        tree_->Branch("GammaIso_DR0p2To0p3",&fGammaIso_DR0p2To0p3,"GammaIso_DR0p2To0p3/F");
        tree_->Branch("GammaIso_DR0p3To0p4",&fGammaIso_DR0p3To0p4,"GammaIso_DR0p3To0p4/F");
        tree_->Branch("GammaIso_DR0p4To0p5",&fGammaIso_DR0p4To0p5,"GammaIso_DR0p4To0p5/F");
        tree_->Branch("NeutralHadronIso_DR0p0To0p1",&fNeutralHadronIso_DR0p0To0p1,"NeutralHadronIso_DR0p0To0p1/F");
        tree_->Branch("NeutralHadronIso_DR0p1To0p2",&fNeutralHadronIso_DR0p1To0p2,"NeutralHadronIso_DR0p1To0p2/F");
        tree_->Branch("NeutralHadronIso_DR0p2To0p3",&fNeutralHadronIso_DR0p2To0p3,"NeutralHadronIso_DR0p2To0p3/F");
        tree_->Branch("NeutralHadronIso_DR0p3To0p4",&fNeutralHadronIso_DR0p3To0p4,"NeutralHadronIso_DR0p3To0p4/F");
        tree_->Branch("NeutralHadronIso_DR0p4To0p5",&fNeutralHadronIso_DR0p4To0p5,"NeutralHadronIso_DR0p4To0p5/F");
        tree_->Branch("PassTriggerDenominator",&fElePassTriggerDenominator,"PassTriggerDenominator/O"); 
        tree_->Branch("IsEB",&fIsEB,"IsEB/O");
        tree_->Branch("IsEE",&fIsEE,"IsEE/O");
        tree_->Branch("SCRawEnergy",&fSCRawEnergy,"SCRawEnergy/F");
        tree_->Branch("E5x5",&fE5x5,"E5x5/F");
        tree_->Branch("NClusters",&fNClusters,"NClusters/F");
        tree_->Branch("EtaSeed",&fEtaSeed,"EtaSeed/F");
        tree_->Branch("PhiSeed",&fPhiSeed,"PhiSeed/F");
        tree_->Branch("ESeed",&fESeed,"ESeed/F");
        tree_->Branch("E3x3Seed",&fE3x3Seed,"E3x3Seed/F");
        tree_->Branch("E5x5Seed",&fE5x5Seed,"E5x5Seed/F");
        tree_->Branch("EMaxSeed",&fEMaxSeed,"EMaxSeed/F");
        tree_->Branch("E2ndSeed",&fE2ndSeed,"E2ndSeed/F");
        tree_->Branch("ETopSeed",&fETopSeed,"ETopSeed/F");
        tree_->Branch("EBottomSeed",&fEBottomSeed,"EBottomSeed/F");
        tree_->Branch("ELeftSeed",&fELeftSeed,"ELeftSeed/F");
        tree_->Branch("ERightSeed",&fERightSeed,"ERightSeed/F");
        tree_->Branch("E2x5MaxSeed",&fE2x5MaxSeed,"E2x5MaxSeed/F");
        tree_->Branch("E2x5TopSeed",&fE2x5TopSeed,"E2x5TopSeed/F");
        tree_->Branch("E2x5BottomSeed",&fE2x5BottomSeed,"E2x5BottomSeed/F");
        tree_->Branch("E2x5LeftSeed",&fE2x5LeftSeed,"E2x5LeftSeed/F");
        tree_->Branch("E2x5RightSeed",&fE2x5RightSeed,"E2x5RightSeed/F");
        tree_->Branch("IEtaSeed",&fIEtaSeed,"IEtaSeed/F");
        tree_->Branch("IPhiSeed",&fIPhiSeed,"IPhiSeed/F");
        tree_->Branch("EtaCrySeed",&fEtaCrySeed,"EtaCrySeed/F");
        tree_->Branch("PhiCrySeed",&fPhiCrySeed,"PhiCrySeed/F");
        tree_->Branch("GsfTrackPIn",&fGsfTrackPIn,"GsfTrackPIn/F");
        tree_->Branch("Charge",&fCharge,"Charge/I");
        tree_->Branch("GeneratedEnergy",&fGeneratedEnergy,"GeneratedEnergy/F");
        tree_->Branch("GeneratedEnergyStatus1",&fGeneratedEnergyStatus1,"GeneratedEnergyStatus1/F");
        tree_->Branch("GeneratedEnergyStatus3",&fGeneratedEnergyStatus3,"GeneratedEnergyStatus3/F");


      } 

      // initialze a ElectronTree
      void InitTree(){
        assert(tree_);
        // don't forget to set pointers to zero before you set address
        // or you will fully appreciate that "ROOT sucks" :)
        InitVariables();
        //Set branch address
        Int_t currentState = gErrorIgnoreLevel;

        tree_->SetBranchAddress("weight",&fWeight);
        tree_->SetBranchAddress("run",&fRunNumber);
        tree_->SetBranchAddress("lumi",&fLumiSectionNumber);
        tree_->SetBranchAddress("event",&fEventNumber);
        tree_->SetBranchAddress("EventNumberParity",&fEleEventNumberParity);
        tree_->SetBranchAddress("pt",&fElePt);
        tree_->SetBranchAddress("eta",&fEleEta);
        tree_->SetBranchAddress("phi",&fElePhi);
        tree_->SetBranchAddress("scEt",&fEleSCEt);
        tree_->SetBranchAddress("scEta",&fEleSCEta);
        tree_->SetBranchAddress("scPhi",&fEleSCPhi);
        tree_->SetBranchAddress("ecalenergy",&fEleEcalEnergy);
        tree_->SetBranchAddress("ecaldriven",&fEleIsEcalDriven);
        tree_->SetBranchAddress("triggerBit",&fEleTriggerBit);
        tree_->SetBranchAddress("rho",&fRho);
        tree_->SetBranchAddress("vertices",&fNVertices);
        tree_->SetBranchAddress("d0",&fEleD0);
        tree_->SetBranchAddress("dz",&fEleDZ);
        tree_->SetBranchAddress("ip3d",&fEleIP3d);
        tree_->SetBranchAddress("ip3ds",&fEleIP3dSig);
        tree_->SetBranchAddress("matchConv",&fEleMatchedConversion);
        tree_->SetBranchAddress("dcot",&fEleConvDCot);
        tree_->SetBranchAddress("dist",&fEleConvDist);
        tree_->SetBranchAddress("missHits",&fEleNMissHits);
        tree_->SetBranchAddress("nbrems",&fEleNBrem);
        tree_->SetBranchAddress("fbrem",&fEleFBrem);
        tree_->SetBranchAddress("EoP",&fEleEOverP);
        tree_->SetBranchAddress("EoPin",&fEleESeedClusterOverPIn);
        tree_->SetBranchAddress("ESeedoPout",&fEleESeedClusterOverPout);
        tree_->SetBranchAddress("EEleoPout",&fEleEEleClusterOverPout);
        tree_->SetBranchAddress("IoEmIoP",&fEleOneOverEMinusOneOverP);
        tree_->SetBranchAddress("deta",&fEleDEtaIn);
        tree_->SetBranchAddress("dphi",&fEleDPhiIn);
        tree_->SetBranchAddress("detacalo",&fEledEtaCalo);
        tree_->SetBranchAddress("dphicalo",&fEledPhiCalo);
        tree_->SetBranchAddress("see",&fEleSigmaIEtaIEta);
        tree_->SetBranchAddress("spp",&fEleSigmaIPhiIPhi);
        tree_->SetBranchAddress("sep",&fEleSigmaIEtaIPhi);
        tree_->SetBranchAddress("etawidth",&fEleSCEtaWidth);
        tree_->SetBranchAddress("phiwidth",&fEleSCPhiWidth);
        tree_->SetBranchAddress("R9",&fEleR9);
        tree_->SetBranchAddress("PreShowerOverRaw",&fElePreShowerOverRaw);
        tree_->SetBranchAddress("HoE",&fEleHoverE);
        tree_->SetBranchAddress("gsfchi2",&fEleGsfTrackChi2OverNdof);
        tree_->SetBranchAddress("kfchi2",&fEleKFTrackChi2OverNDoF);
        tree_->SetBranchAddress("kfhits",&fEleKFTrackNHits);
        tree_->SetBranchAddress("kflayers",&fEleKFTrackNLayersWithMeasurement);
        tree_->SetBranchAddress("OneMinusSeedE1x5OverE5x5",&fEleOneMinusSeedE1x5OverE5x5);
        tree_->SetBranchAddress("PFMVA",&fElePFMVA);
        tree_->SetBranchAddress("trkIso03",&fEleTrkIso03);
        tree_->SetBranchAddress("ecalIso03",&fEleEMIso03);
        tree_->SetBranchAddress("hcalIso03",&fEleHadIso03);
        tree_->SetBranchAddress("trkIso04",&fEleTrkIso04);
        tree_->SetBranchAddress("ecalIso04",&fEleEMIso04);
        tree_->SetBranchAddress("hcalIso04",&fEleHadIso04);
        tree_->SetBranchAddress("pfIso04",&fElePFIso04);
        tree_->SetBranchAddress("ChargedIso_DR0p0To0p1",&fChargedIso_DR0p0To0p1);
        tree_->SetBranchAddress("ChargedIso_DR0p1To0p2",&fChargedIso_DR0p1To0p2);
        tree_->SetBranchAddress("ChargedIso_DR0p2To0p3",&fChargedIso_DR0p2To0p3);
        tree_->SetBranchAddress("ChargedIso_DR0p3To0p4",&fChargedIso_DR0p3To0p4);
        tree_->SetBranchAddress("ChargedIso_DR0p4To0p5",&fChargedIso_DR0p4To0p5);
        tree_->SetBranchAddress("GammaIso_DR0p0To0p1",&fGammaIso_DR0p0To0p1);
        tree_->SetBranchAddress("GammaIso_DR0p1To0p2",&fGammaIso_DR0p1To0p2);
        tree_->SetBranchAddress("GammaIso_DR0p2To0p3",&fGammaIso_DR0p2To0p3);
        tree_->SetBranchAddress("GammaIso_DR0p3To0p4",&fGammaIso_DR0p3To0p4);
        tree_->SetBranchAddress("GammaIso_DR0p4To0p5",&fGammaIso_DR0p4To0p5);
        tree_->SetBranchAddress("NeutralHadronIso_DR0p0To0p1",&fNeutralHadronIso_DR0p0To0p1);
        tree_->SetBranchAddress("NeutralHadronIso_DR0p1To0p2",&fNeutralHadronIso_DR0p1To0p2);
        tree_->SetBranchAddress("NeutralHadronIso_DR0p2To0p3",&fNeutralHadronIso_DR0p2To0p3);
        tree_->SetBranchAddress("NeutralHadronIso_DR0p3To0p4",&fNeutralHadronIso_DR0p3To0p4);
        tree_->SetBranchAddress("NeutralHadronIso_DR0p4To0p5",&fNeutralHadronIso_DR0p4To0p5);
        tree_->SetBranchAddress("PassTriggerDenominator",&fElePassTriggerDenominator);
        tree_->SetBranchAddress("IsEB",&fIsEB);
        tree_->SetBranchAddress("IsEE",&fIsEE);
        tree_->SetBranchAddress("SCRawEnergy",&fSCRawEnergy);
        tree_->SetBranchAddress("E5x5",&fE5x5);
        tree_->SetBranchAddress("NClusters",&fNClusters);
        tree_->SetBranchAddress("EtaSeed",&fEtaSeed);
        tree_->SetBranchAddress("PhiSeed",&fPhiSeed);
        tree_->SetBranchAddress("ESeed",&fESeed);
        tree_->SetBranchAddress("E3x3Seed",&fE3x3Seed);
        tree_->SetBranchAddress("E5x5Seed",&fE5x5Seed);
        tree_->SetBranchAddress("EMaxSeed",&fEMaxSeed);
        tree_->SetBranchAddress("E2ndSeed",&fE2ndSeed);
        tree_->SetBranchAddress("ETopSeed",&fETopSeed);
        tree_->SetBranchAddress("EBottomSeed",&fEBottomSeed);
        tree_->SetBranchAddress("ELeftSeed",&fELeftSeed);
        tree_->SetBranchAddress("ERightSeed",&fERightSeed);
        tree_->SetBranchAddress("E2x5MaxSeed",&fE2x5MaxSeed);
        tree_->SetBranchAddress("E2x5TopSeed",&fE2x5TopSeed);
        tree_->SetBranchAddress("E2x5BottomSeed",&fE2x5BottomSeed);
        tree_->SetBranchAddress("E2x5LeftSeed",&fE2x5LeftSeed);
        tree_->SetBranchAddress("E2x5RightSeed",&fE2x5RightSeed);
        tree_->SetBranchAddress("IEtaSeed",&fIEtaSeed);
        tree_->SetBranchAddress("IPhiSeed",&fIPhiSeed);
        tree_->SetBranchAddress("EtaCrySeed",&fEtaCrySeed);
        tree_->SetBranchAddress("PhiCrySeed",&fPhiCrySeed);
        tree_->SetBranchAddress("GsfTrackPIn",&fGsfTrackPIn);
        tree_->SetBranchAddress("Charge",&fCharge);
        tree_->SetBranchAddress("GeneratedEnergy",&fGeneratedEnergy);
        tree_->SetBranchAddress("GeneratedEnergyStatus1",&fGeneratedEnergyStatus1);
        tree_->SetBranchAddress("GeneratedEnergyStatus3",&fGeneratedEnergyStatus3);

        gErrorIgnoreLevel = currentState;
      }

  }; 

}

#endif
Revision:	1.3
Committed:	Sun Jul 29 10:34:09 2012 UTC (12 years, 9 months ago) by sixie
Content type:	text/plain
Branch:	MAIN
Changes since 1.2:	+147 -23 lines
Log Message:	add variables for regression
#	User	Rev	Content
1	sixie	1.1	#ifndef ElectronTree_H
2			#define ElectronTree_H
3
4			#include "TFile.h"
5			#include "TTree.h"
6			#include "TError.h"
7			#include <cmath>
8			#include "assert.h"
9
10	sixie	1.2	namespace citana
11			{
12			class ElectronTree {
13
14			public:
15
16			/// bit map
17			/// DON'T CHANGE ORDER
18
19			//*******************************************
20			//=== ElectronTriggerBits ====
21			//*******************************************
22			enum ElectronTriggerBits { kEleTrigger_Ele = 0x000001,
23			kEleTrigger_Ele_CaloIdL_CaloIsoVL = 0x000002,
24			kEleTrigger_Ele_CaloIdT_TrkIdVL_CaloIsoVL_TrkIsoVL = 0x000004
25			};
26
27			/// variables
28			Float_t fWeight;
29			UInt_t fRunNumber;
30			UInt_t fLumiSectionNumber;
31			UInt_t fEventNumber;
32			Bool_t fEleEventNumberParity;
33			Float_t fElePt;
34			Float_t fEleEta;
35			Float_t fElePhi;
36			Float_t fEleSCEt;
37			Float_t fEleSCEta;
38			Float_t fEleSCPhi;
39			Float_t fEleEcalEnergy;
40			Bool_t fEleIsEcalDriven;
41			UInt_t fEleTriggerBit;
42			Float_t fRho;
43			UInt_t fNVertices;
44
45			// Conversion and IP
46			Float_t fEleD0;
47			Float_t fEleDZ;
48			Float_t fEleIP3d;
49			Float_t fEleIP3dSig;
50			Bool_t fEleMatchedConversion;
51			Float_t fEleConvDCot;
52			Float_t fEleConvDist;
53			UInt_t fEleNMissHits;
54
55			// E/P variables
56			Float_t fEleNBrem;
57			Float_t fEleFBrem;
58			Float_t fEleEOverP;
59			Float_t fEleESeedClusterOverPIn;
60			Float_t fEleESeedClusterOverPout;
61			Float_t fEleEEleClusterOverPout;
62			Float_t fEleOneOverEMinusOneOverP;
63
64			// track cluster matching
65			Float_t fEleDEtaIn;
66			Float_t fEleDPhiIn;
67			Float_t fEledEtaCalo;
68			Float_t fEledPhiCalo;
69
70			//shower shape
71			Float_t fEleSigmaIEtaIEta;
72			Float_t fEleSigmaIPhiIPhi;
73			Float_t fEleSigmaIEtaIPhi;
74			Float_t fEleSCEtaWidth;
75			Float_t fEleSCPhiWidth;
76			Float_t fEleR9;
77			Float_t fElePreShowerOverRaw;
78			Float_t fEleHoverE;
79
80			//track quality
81			Float_t fEleGsfTrackChi2OverNdof;
82			Float_t fEleKFTrackChi2OverNDoF;
83			Int_t fEleKFTrackNHits;
84			Int_t fEleKFTrackNLayersWithMeasurement;
85			Float_t fEleOneMinusSeedE1x5OverE5x5;
86
87			//Isolation Variables
88			Float_t fElePFMVA;
89			Float_t fEleTrkIso03;
90			Float_t fEleEMIso03;
91			Float_t fEleHadIso03;
92			Float_t fEleTrkIso04;
93			Float_t fEleEMIso04;
94			Float_t fEleHadIso04;
95			Float_t fElePFIso04;
96			Float_t fChargedIso_DR0p0To0p1;
97			Float_t fChargedIso_DR0p1To0p2;
98			Float_t fChargedIso_DR0p2To0p3;
99			Float_t fChargedIso_DR0p3To0p4;
100			Float_t fChargedIso_DR0p4To0p5;
101			Float_t fGammaIso_DR0p0To0p1;
102			Float_t fGammaIso_DR0p1To0p2;
103			Float_t fGammaIso_DR0p2To0p3;
104			Float_t fGammaIso_DR0p3To0p4;
105			Float_t fGammaIso_DR0p4To0p5;
106			Float_t fNeutralHadronIso_DR0p0To0p1;
107			Float_t fNeutralHadronIso_DR0p1To0p2;
108			Float_t fNeutralHadronIso_DR0p2To0p3;
109			Float_t fNeutralHadronIso_DR0p3To0p4;
110			Float_t fNeutralHadronIso_DR0p4To0p5;
111
112			Bool_t fElePassTriggerDenominator;
113
114	sixie	1.3	//Regression Variables
115			Bool_t fIsEB;
116			Bool_t fIsEE;
117			Float_t fSCRawEnergy;
118			Float_t fE5x5;
119			Float_t fNClusters;
120			Float_t fEtaSeed;
121			Float_t fPhiSeed;
122			Float_t fESeed;
123			Float_t fE3x3Seed;
124			Float_t fE5x5Seed;
125			Float_t fEMaxSeed;
126			Float_t fE2ndSeed;
127			Float_t fETopSeed;
128			Float_t fEBottomSeed;
129			Float_t fELeftSeed;
130			Float_t fERightSeed;
131			Float_t fE2x5MaxSeed;
132			Float_t fE2x5TopSeed;
133			Float_t fE2x5BottomSeed;
134			Float_t fE2x5LeftSeed;
135			Float_t fE2x5RightSeed;
136			Float_t fIEtaSeed;
137			Float_t fIPhiSeed;
138			Float_t fEtaCrySeed;
139			Float_t fPhiCrySeed;
140			Float_t fGsfTrackPIn;
141			Int_t fCharge;
142			Float_t fGeneratedEnergy;
143			Float_t fGeneratedEnergyStatus1;
144			Float_t fGeneratedEnergyStatus3;
145
146	sixie	1.2	public:
147			/// this is the main element
148			TTree *tree_;
149			TFile *f_;
150	sixie	1.1
151	sixie	1.2	/// hold the names of variables to facilitate things (filled during Init)
152			std::vector<std::string> variables_;
153	sixie	1.1
154	sixie	1.2	/// default constructor
155			ElectronTree() {};
156			/// default destructor
157			~ElectronTree(){
158			if (f_) f_->Close();
159			};
160	sixie	1.1
161	sixie	1.2	/// initialize varibles and fill list of available variables
162			void InitVariables() {
163			fWeight = 0.0;
164			fRunNumber = 0.0;
165			fLumiSectionNumber = 0.0;
166			fEventNumber = 0.0;
167			fEleEventNumberParity = 0.0;
168			fElePt = 0.0;
169	sixie	1.3	fEleEta = 0.0;
170			fElePhi = 0.0;
171			fEleSCEt = 0.0;
172	sixie	1.2	fEleSCEta = 0.0;
173			fEleSCPhi = 0.0;
174			fEleEcalEnergy = 0.0;
175	sixie	1.3	fEleIsEcalDriven = 0.0;
176	sixie	1.2	fEleTriggerBit = 0.0;
177			fRho = 0.0;
178			fNVertices = 0.0;
179			fEleD0 = 0.0;
180			fEleDZ = 0.0;
181	sixie	1.3	fEleIP3d = 0.0;
182	sixie	1.2	fEleIP3dSig = 0.0;
183			fEleMatchedConversion = 0.0;
184			fEleConvDCot = 0.0;
185			fEleConvDist = 0.0;
186			fEleNMissHits = 0.0;
187			fEleNBrem = 0.0;
188			fEleFBrem = 0.0;
189			fEleEOverP = 0.0;
190	sixie	1.3	fEleESeedClusterOverPIn = 0.0;
191			fEleESeedClusterOverPout = 0.0;
192			fEleEEleClusterOverPout = 0.0;
193			fEleOneOverEMinusOneOverP = 0.0;
194	sixie	1.2	fEleDEtaIn = 0.0;
195			fEleDPhiIn = 0.0;
196			fEledEtaCalo = 0.0;
197			fEledPhiCalo = 0.0;
198			fEleSigmaIEtaIEta = 0.0;
199			fEleSigmaIPhiIPhi = 0.0;
200	sixie	1.3	fEleSigmaIEtaIPhi = 0.0;
201	sixie	1.2	fEleSCEtaWidth = 0.0;
202			fEleSCPhiWidth = 0.0;
203			fEleR9 = 0.0;
204			fElePreShowerOverRaw = 0.0;
205			fEleHoverE = 0.0;
206	sixie	1.3	fEleGsfTrackChi2OverNdof = 0.0;
207	sixie	1.2	fEleKFTrackChi2OverNDoF = 0.0;
208	sixie	1.3	fEleKFTrackNHits = 0.0;
209	sixie	1.2	fEleKFTrackNLayersWithMeasurement = 0.0;
210	sixie	1.3	fEleOneMinusSeedE1x5OverE5x5 = 0.0;
211			fElePFMVA = 0.0;
212	sixie	1.2	fEleTrkIso03 = 0.0;
213			fEleEMIso03 = 0.0;
214			fEleHadIso03 = 0.0;
215			fEleTrkIso04 = 0.0;
216			fEleEMIso04 = 0.0;
217	sixie	1.3	fEleHadIso04 = 0.0;
218	sixie	1.2	fElePFIso04 = 0.0;
219			fChargedIso_DR0p0To0p1 = 0.0;
220			fChargedIso_DR0p1To0p2 = 0.0;
221			fChargedIso_DR0p2To0p3 = 0.0;
222			fChargedIso_DR0p3To0p4 = 0.0;
223			fChargedIso_DR0p4To0p5 = 0.0;
224			fGammaIso_DR0p0To0p1 = 0.0;
225			fGammaIso_DR0p1To0p2 = 0.0;
226			fGammaIso_DR0p2To0p3 = 0.0;
227			fGammaIso_DR0p3To0p4 = 0.0;
228			fGammaIso_DR0p4To0p5 = 0.0;
229	sixie	1.3	fNeutralHadronIso_DR0p0To0p1 = 0.0;
230			fNeutralHadronIso_DR0p1To0p2 = 0.0;
231			fNeutralHadronIso_DR0p2To0p3 = 0.0;
232			fNeutralHadronIso_DR0p3To0p4 = 0.0;
233			fNeutralHadronIso_DR0p4To0p5 = 0.0;
234			fElePassTriggerDenominator = 0.0;
235			fIsEB = 0.0;
236			fIsEE = 0.0;
237			fSCRawEnergy = 0.0;
238			fE5x5 = 0.0;
239			fNClusters = 0.0;
240			fEtaSeed = 0.0;
241			fPhiSeed = 0.0;
242			fESeed = 0.0;
243			fE3x3Seed = 0.0;
244			fE5x5Seed = 0.0;
245			fEMaxSeed = 0.0;
246			fE2ndSeed = 0.0;
247			fETopSeed = 0.0;
248			fEBottomSeed = 0.0;
249			fELeftSeed = 0.0;
250			fERightSeed = 0.0;
251			fE2x5MaxSeed = 0.0;
252			fE2x5TopSeed = 0.0;
253			fE2x5BottomSeed = 0.0;
254			fE2x5LeftSeed = 0.0;
255			fE2x5RightSeed = 0.0;
256			fIEtaSeed = 0.0;
257			fIPhiSeed = 0.0;
258			fEtaCrySeed = 0.0;
259			fPhiCrySeed = 0.0;
260			fGsfTrackPIn = 0.0;
261			fCharge = 0.0;
262			fGeneratedEnergy = 0.0;
263			fGeneratedEnergyStatus1 = 0.0;
264			fGeneratedEnergyStatus3 = 0.0;
265
266
267	sixie	1.2	}
268	sixie	1.1
269	sixie	1.2	/// load a ElectronTree
270			void LoadTree(const char* file){
271			f_ = TFile::Open(file);
272			assert(f_);
273			tree_ = dynamic_cast<TTree*>(f_->Get("Electrons"));
274			assert(tree_);
275			}
276	sixie	1.1
277	sixie	1.2	/// create a ElectronTree
278			void CreateTree(){
279			tree_ = new TTree("Electrons","Electrons");
280			f_ = 0;
281
282			//book the branches
283			tree_->Branch("weight",&fWeight,"weight/F");
284			tree_->Branch("run",&fRunNumber,"run/i");
285			tree_->Branch("lumi",&fLumiSectionNumber,"lumi/i");
286			tree_->Branch("event",&fEventNumber,"event/i");
287			tree_->Branch("EventNumberParity",&fEleEventNumberParity,"EventNumberParity/O");
288			tree_->Branch("pt",&fElePt,"pt/F");
289			tree_->Branch("eta",&fEleEta,"eta/F");
290			tree_->Branch("phi",&fElePhi,"phi/F");
291			tree_->Branch("scEt",&fEleSCEt,"scEt/F");
292			tree_->Branch("scEta",&fEleSCEta,"scEta/F");
293			tree_->Branch("scPhi",&fEleSCPhi,"scPhi/F");
294			tree_->Branch("ecalenergy",&fEleEcalEnergy,"ecalenergy/F");
295			tree_->Branch("ecaldriven",&fEleIsEcalDriven,"ecaldriven/O");
296			tree_->Branch("triggerBit",&fEleTriggerBit,"triggerBit/i");
297			tree_->Branch("rho",&fRho,"rho/F");
298			tree_->Branch("vertices",&fNVertices,"vertices/i");
299			tree_->Branch("d0",&fEleD0,"d0/F");
300			tree_->Branch("dz",&fEleDZ,"dz/F");
301			tree_->Branch("ip3d",&fEleIP3d,"ip3d/F");
302			tree_->Branch("ip3ds",&fEleIP3dSig,"ip3ds/F");
303			tree_->Branch("matchConv",&fEleMatchedConversion,"matchConv/O");
304			tree_->Branch("dcot",&fEleConvDCot,"dcot/F");
305			tree_->Branch("dist",&fEleConvDist,"dist/F");
306			tree_->Branch("missHits",&fEleNMissHits,"NMissHits/i");
307			tree_->Branch("nbrems",&fEleNBrem,"nbrems/F");
308			tree_->Branch("fbrem",&fEleFBrem,"fbrem/F");
309			tree_->Branch("EoP",&fEleEOverP,"EoP/F");
310			tree_->Branch("EoPin",&fEleESeedClusterOverPIn,"EoPin/F");
311			tree_->Branch("ESeedoPout",&fEleESeedClusterOverPout,"ESeedoPout/F");
312			tree_->Branch("EEleoPout",&fEleEEleClusterOverPout,"EEleoPout/F");
313	sixie	1.3	// tree_->Branch("eleEoPout",&fEleEEleClusterOverPout,"eleEoPout/F"); //old name of the branch
314			tree_->Branch("IoEmIoP",&fEleOneOverEMinusOneOverP,"IoEmIoP/F");
315	sixie	1.2	tree_->Branch("deta",&fEleDEtaIn,"deta/F");
316			tree_->Branch("dphi",&fEleDPhiIn,"dphi/F");
317			tree_->Branch("detacalo",&fEledEtaCalo,"detacalo/F");
318			tree_->Branch("dphicalo",&fEledPhiCalo,"dphicalo/F");
319			tree_->Branch("see",&fEleSigmaIEtaIEta,"see/F");
320			tree_->Branch("spp",&fEleSigmaIPhiIPhi,"spp/F");
321			tree_->Branch("sep",&fEleSigmaIEtaIPhi,"sep/F");
322			tree_->Branch("etawidth",&fEleSCEtaWidth,"etawidth/F");
323			tree_->Branch("phiwidth",&fEleSCPhiWidth,"phiwidth/F");
324			tree_->Branch("R9",&fEleR9,"R9/F");
325			tree_->Branch("PreShowerOverRaw",&fElePreShowerOverRaw,"PreShowerOverRaw/F");
326			tree_->Branch("HoE",&fEleHoverE,"HoE/F");
327			tree_->Branch("gsfchi2",&fEleGsfTrackChi2OverNdof,"gsfchi2/F");
328			tree_->Branch("kfchi2",&fEleKFTrackChi2OverNDoF,"kfchi2/F");
329			tree_->Branch("kfhits",&fEleKFTrackNHits,"kfhits/I");
330			tree_->Branch("kflayers",&fEleKFTrackNLayersWithMeasurement,"kflayers/I");
331			tree_->Branch("OneMinusSeedE1x5OverE5x5",&fEleOneMinusSeedE1x5OverE5x5,"OneMinusSeedE1x5OverE5x5/F");
332			tree_->Branch("PFMVA",&fElePFMVA,"PFMVA/F");
333			tree_->Branch("trkIso03",&fEleTrkIso03,"trkIso03/F");
334			tree_->Branch("ecalIso03",&fEleEMIso03,"ecalIso03/F");
335			tree_->Branch("hcalIso03",&fEleHadIso03,"hcalIso03/F");
336			tree_->Branch("trkIso04",&fEleTrkIso04,"trkIso04/F");
337			tree_->Branch("ecalIso04",&fEleEMIso04,"ecalIso04/F");
338			tree_->Branch("hcalIso04",&fEleHadIso04,"hcalIso04/F");
339			tree_->Branch("pfIso04",&fElePFIso04,"pfIso04/F");
340			tree_->Branch("ChargedIso_DR0p0To0p1",&fChargedIso_DR0p0To0p1,"ChargedIso_DR0p0To0p1/F");
341			tree_->Branch("ChargedIso_DR0p1To0p2",&fChargedIso_DR0p1To0p2,"ChargedIso_DR0p1To0p2/F");
342			tree_->Branch("ChargedIso_DR0p2To0p3",&fChargedIso_DR0p2To0p3,"ChargedIso_DR0p2To0p3/F");
343			tree_->Branch("ChargedIso_DR0p3To0p4",&fChargedIso_DR0p3To0p4,"ChargedIso_DR0p3To0p4/F");
344			tree_->Branch("ChargedIso_DR0p4To0p5",&fChargedIso_DR0p4To0p5,"ChargedIso_DR0p4To0p5/F");
345			tree_->Branch("GammaIso_DR0p0To0p1",&fGammaIso_DR0p0To0p1,"GammaIso_DR0p0To0p1/F");
346			tree_->Branch("GammaIso_DR0p1To0p2",&fGammaIso_DR0p1To0p2,"GammaIso_DR0p1To0p2/F");
347			tree_->Branch("GammaIso_DR0p2To0p3",&fGammaIso_DR0p2To0p3,"GammaIso_DR0p2To0p3/F");
348			tree_->Branch("GammaIso_DR0p3To0p4",&fGammaIso_DR0p3To0p4,"GammaIso_DR0p3To0p4/F");
349			tree_->Branch("GammaIso_DR0p4To0p5",&fGammaIso_DR0p4To0p5,"GammaIso_DR0p4To0p5/F");
350			tree_->Branch("NeutralHadronIso_DR0p0To0p1",&fNeutralHadronIso_DR0p0To0p1,"NeutralHadronIso_DR0p0To0p1/F");
351			tree_->Branch("NeutralHadronIso_DR0p1To0p2",&fNeutralHadronIso_DR0p1To0p2,"NeutralHadronIso_DR0p1To0p2/F");
352			tree_->Branch("NeutralHadronIso_DR0p2To0p3",&fNeutralHadronIso_DR0p2To0p3,"NeutralHadronIso_DR0p2To0p3/F");
353			tree_->Branch("NeutralHadronIso_DR0p3To0p4",&fNeutralHadronIso_DR0p3To0p4,"NeutralHadronIso_DR0p3To0p4/F");
354			tree_->Branch("NeutralHadronIso_DR0p4To0p5",&fNeutralHadronIso_DR0p4To0p5,"NeutralHadronIso_DR0p4To0p5/F");
355			tree_->Branch("PassTriggerDenominator",&fElePassTriggerDenominator,"PassTriggerDenominator/O");
356	sixie	1.3	tree_->Branch("IsEB",&fIsEB,"IsEB/O");
357			tree_->Branch("IsEE",&fIsEE,"IsEE/O");
358			tree_->Branch("SCRawEnergy",&fSCRawEnergy,"SCRawEnergy/F");
359			tree_->Branch("E5x5",&fE5x5,"E5x5/F");
360			tree_->Branch("NClusters",&fNClusters,"NClusters/F");
361			tree_->Branch("EtaSeed",&fEtaSeed,"EtaSeed/F");
362			tree_->Branch("PhiSeed",&fPhiSeed,"PhiSeed/F");
363			tree_->Branch("ESeed",&fESeed,"ESeed/F");
364			tree_->Branch("E3x3Seed",&fE3x3Seed,"E3x3Seed/F");
365			tree_->Branch("E5x5Seed",&fE5x5Seed,"E5x5Seed/F");
366			tree_->Branch("EMaxSeed",&fEMaxSeed,"EMaxSeed/F");
367			tree_->Branch("E2ndSeed",&fE2ndSeed,"E2ndSeed/F");
368			tree_->Branch("ETopSeed",&fETopSeed,"ETopSeed/F");
369			tree_->Branch("EBottomSeed",&fEBottomSeed,"EBottomSeed/F");
370			tree_->Branch("ELeftSeed",&fELeftSeed,"ELeftSeed/F");
371			tree_->Branch("ERightSeed",&fERightSeed,"ERightSeed/F");
372			tree_->Branch("E2x5MaxSeed",&fE2x5MaxSeed,"E2x5MaxSeed/F");
373			tree_->Branch("E2x5TopSeed",&fE2x5TopSeed,"E2x5TopSeed/F");
374			tree_->Branch("E2x5BottomSeed",&fE2x5BottomSeed,"E2x5BottomSeed/F");
375			tree_->Branch("E2x5LeftSeed",&fE2x5LeftSeed,"E2x5LeftSeed/F");
376			tree_->Branch("E2x5RightSeed",&fE2x5RightSeed,"E2x5RightSeed/F");
377			tree_->Branch("IEtaSeed",&fIEtaSeed,"IEtaSeed/F");
378			tree_->Branch("IPhiSeed",&fIPhiSeed,"IPhiSeed/F");
379			tree_->Branch("EtaCrySeed",&fEtaCrySeed,"EtaCrySeed/F");
380			tree_->Branch("PhiCrySeed",&fPhiCrySeed,"PhiCrySeed/F");
381			tree_->Branch("GsfTrackPIn",&fGsfTrackPIn,"GsfTrackPIn/F");
382			tree_->Branch("Charge",&fCharge,"Charge/I");
383			tree_->Branch("GeneratedEnergy",&fGeneratedEnergy,"GeneratedEnergy/F");
384			tree_->Branch("GeneratedEnergyStatus1",&fGeneratedEnergyStatus1,"GeneratedEnergyStatus1/F");
385			tree_->Branch("GeneratedEnergyStatus3",&fGeneratedEnergyStatus3,"GeneratedEnergyStatus3/F");
386	sixie	1.2
387
388			}
389
390			// initialze a ElectronTree
391			void InitTree(){
392			assert(tree_);
393			// don't forget to set pointers to zero before you set address
394			// or you will fully appreciate that "ROOT sucks" :)
395			InitVariables();
396			//Set branch address
397			Int_t currentState = gErrorIgnoreLevel;
398
399			tree_->SetBranchAddress("weight",&fWeight);
400			tree_->SetBranchAddress("run",&fRunNumber);
401			tree_->SetBranchAddress("lumi",&fLumiSectionNumber);
402			tree_->SetBranchAddress("event",&fEventNumber);
403			tree_->SetBranchAddress("EventNumberParity",&fEleEventNumberParity);
404			tree_->SetBranchAddress("pt",&fElePt);
405			tree_->SetBranchAddress("eta",&fEleEta);
406			tree_->SetBranchAddress("phi",&fElePhi);
407			tree_->SetBranchAddress("scEt",&fEleSCEt);
408			tree_->SetBranchAddress("scEta",&fEleSCEta);
409			tree_->SetBranchAddress("scPhi",&fEleSCPhi);
410			tree_->SetBranchAddress("ecalenergy",&fEleEcalEnergy);
411			tree_->SetBranchAddress("ecaldriven",&fEleIsEcalDriven);
412			tree_->SetBranchAddress("triggerBit",&fEleTriggerBit);
413			tree_->SetBranchAddress("rho",&fRho);
414			tree_->SetBranchAddress("vertices",&fNVertices);
415			tree_->SetBranchAddress("d0",&fEleD0);
416			tree_->SetBranchAddress("dz",&fEleDZ);
417			tree_->SetBranchAddress("ip3d",&fEleIP3d);
418			tree_->SetBranchAddress("ip3ds",&fEleIP3dSig);
419			tree_->SetBranchAddress("matchConv",&fEleMatchedConversion);
420			tree_->SetBranchAddress("dcot",&fEleConvDCot);
421			tree_->SetBranchAddress("dist",&fEleConvDist);
422			tree_->SetBranchAddress("missHits",&fEleNMissHits);
423			tree_->SetBranchAddress("nbrems",&fEleNBrem);
424			tree_->SetBranchAddress("fbrem",&fEleFBrem);
425			tree_->SetBranchAddress("EoP",&fEleEOverP);
426			tree_->SetBranchAddress("EoPin",&fEleESeedClusterOverPIn);
427			tree_->SetBranchAddress("ESeedoPout",&fEleESeedClusterOverPout);
428			tree_->SetBranchAddress("EEleoPout",&fEleEEleClusterOverPout);
429			tree_->SetBranchAddress("IoEmIoP",&fEleOneOverEMinusOneOverP);
430			tree_->SetBranchAddress("deta",&fEleDEtaIn);
431			tree_->SetBranchAddress("dphi",&fEleDPhiIn);
432			tree_->SetBranchAddress("detacalo",&fEledEtaCalo);
433			tree_->SetBranchAddress("dphicalo",&fEledPhiCalo);
434			tree_->SetBranchAddress("see",&fEleSigmaIEtaIEta);
435			tree_->SetBranchAddress("spp",&fEleSigmaIPhiIPhi);
436			tree_->SetBranchAddress("sep",&fEleSigmaIEtaIPhi);
437			tree_->SetBranchAddress("etawidth",&fEleSCEtaWidth);
438			tree_->SetBranchAddress("phiwidth",&fEleSCPhiWidth);
439			tree_->SetBranchAddress("R9",&fEleR9);
440			tree_->SetBranchAddress("PreShowerOverRaw",&fElePreShowerOverRaw);
441			tree_->SetBranchAddress("HoE",&fEleHoverE);
442			tree_->SetBranchAddress("gsfchi2",&fEleGsfTrackChi2OverNdof);
443			tree_->SetBranchAddress("kfchi2",&fEleKFTrackChi2OverNDoF);
444			tree_->SetBranchAddress("kfhits",&fEleKFTrackNHits);
445			tree_->SetBranchAddress("kflayers",&fEleKFTrackNLayersWithMeasurement);
446			tree_->SetBranchAddress("OneMinusSeedE1x5OverE5x5",&fEleOneMinusSeedE1x5OverE5x5);
447			tree_->SetBranchAddress("PFMVA",&fElePFMVA);
448			tree_->SetBranchAddress("trkIso03",&fEleTrkIso03);
449			tree_->SetBranchAddress("ecalIso03",&fEleEMIso03);
450			tree_->SetBranchAddress("hcalIso03",&fEleHadIso03);
451			tree_->SetBranchAddress("trkIso04",&fEleTrkIso04);
452			tree_->SetBranchAddress("ecalIso04",&fEleEMIso04);
453			tree_->SetBranchAddress("hcalIso04",&fEleHadIso04);
454			tree_->SetBranchAddress("pfIso04",&fElePFIso04);
455			tree_->SetBranchAddress("ChargedIso_DR0p0To0p1",&fChargedIso_DR0p0To0p1);
456			tree_->SetBranchAddress("ChargedIso_DR0p1To0p2",&fChargedIso_DR0p1To0p2);
457			tree_->SetBranchAddress("ChargedIso_DR0p2To0p3",&fChargedIso_DR0p2To0p3);
458			tree_->SetBranchAddress("ChargedIso_DR0p3To0p4",&fChargedIso_DR0p3To0p4);
459			tree_->SetBranchAddress("ChargedIso_DR0p4To0p5",&fChargedIso_DR0p4To0p5);
460			tree_->SetBranchAddress("GammaIso_DR0p0To0p1",&fGammaIso_DR0p0To0p1);
461			tree_->SetBranchAddress("GammaIso_DR0p1To0p2",&fGammaIso_DR0p1To0p2);
462			tree_->SetBranchAddress("GammaIso_DR0p2To0p3",&fGammaIso_DR0p2To0p3);
463			tree_->SetBranchAddress("GammaIso_DR0p3To0p4",&fGammaIso_DR0p3To0p4);
464			tree_->SetBranchAddress("GammaIso_DR0p4To0p5",&fGammaIso_DR0p4To0p5);
465			tree_->SetBranchAddress("NeutralHadronIso_DR0p0To0p1",&fNeutralHadronIso_DR0p0To0p1);
466			tree_->SetBranchAddress("NeutralHadronIso_DR0p1To0p2",&fNeutralHadronIso_DR0p1To0p2);
467			tree_->SetBranchAddress("NeutralHadronIso_DR0p2To0p3",&fNeutralHadronIso_DR0p2To0p3);
468			tree_->SetBranchAddress("NeutralHadronIso_DR0p3To0p4",&fNeutralHadronIso_DR0p3To0p4);
469			tree_->SetBranchAddress("NeutralHadronIso_DR0p4To0p5",&fNeutralHadronIso_DR0p4To0p5);
470			tree_->SetBranchAddress("PassTriggerDenominator",&fElePassTriggerDenominator);
471	sixie	1.3	tree_->SetBranchAddress("IsEB",&fIsEB);
472			tree_->SetBranchAddress("IsEE",&fIsEE);
473			tree_->SetBranchAddress("SCRawEnergy",&fSCRawEnergy);
474			tree_->SetBranchAddress("E5x5",&fE5x5);
475			tree_->SetBranchAddress("NClusters",&fNClusters);
476			tree_->SetBranchAddress("EtaSeed",&fEtaSeed);
477			tree_->SetBranchAddress("PhiSeed",&fPhiSeed);
478			tree_->SetBranchAddress("ESeed",&fESeed);
479			tree_->SetBranchAddress("E3x3Seed",&fE3x3Seed);
480			tree_->SetBranchAddress("E5x5Seed",&fE5x5Seed);
481			tree_->SetBranchAddress("EMaxSeed",&fEMaxSeed);
482			tree_->SetBranchAddress("E2ndSeed",&fE2ndSeed);
483			tree_->SetBranchAddress("ETopSeed",&fETopSeed);
484			tree_->SetBranchAddress("EBottomSeed",&fEBottomSeed);
485			tree_->SetBranchAddress("ELeftSeed",&fELeftSeed);
486			tree_->SetBranchAddress("ERightSeed",&fERightSeed);
487			tree_->SetBranchAddress("E2x5MaxSeed",&fE2x5MaxSeed);
488			tree_->SetBranchAddress("E2x5TopSeed",&fE2x5TopSeed);
489			tree_->SetBranchAddress("E2x5BottomSeed",&fE2x5BottomSeed);
490			tree_->SetBranchAddress("E2x5LeftSeed",&fE2x5LeftSeed);
491			tree_->SetBranchAddress("E2x5RightSeed",&fE2x5RightSeed);
492			tree_->SetBranchAddress("IEtaSeed",&fIEtaSeed);
493			tree_->SetBranchAddress("IPhiSeed",&fIPhiSeed);
494			tree_->SetBranchAddress("EtaCrySeed",&fEtaCrySeed);
495			tree_->SetBranchAddress("PhiCrySeed",&fPhiCrySeed);
496			tree_->SetBranchAddress("GsfTrackPIn",&fGsfTrackPIn);
497			tree_->SetBranchAddress("Charge",&fCharge);
498			tree_->SetBranchAddress("GeneratedEnergy",&fGeneratedEnergy);
499			tree_->SetBranchAddress("GeneratedEnergyStatus1",&fGeneratedEnergyStatus1);
500			tree_->SetBranchAddress("GeneratedEnergyStatus3",&fGeneratedEnergyStatus3);
501	sixie	1.1
502	sixie	1.2	gErrorIgnoreLevel = currentState;
503			}
504	sixie	1.1
505	sixie	1.2	};
506	sixie	1.1
507	sixie	1.2	}
508	sixie	1.1
509			#endif