dd/d77/L1PromptAnalysis_8h_source.html

//-------------------------------------------------

//

//   \class L1PromptAnalysis

//

//   I. Mikulec            HEPHY Vienna

//

//--------------------------------------------------

#ifndef L1_PROMPT_ANALYSIS_H

#define L1_PROMPT_ANALYSIS_H


//---------------

// C++ Headers --

//---------------


#include <memory>

#include <string>


//----------------------

// Base Class Headers --

//----------------------

#include "FWCore/Framework/interface/EDAnalyzer.h"


#include "FWCore/Framework/interface/Frameworkfwd.h"

#include "FWCore/Framework/interface/MakerMacros.h"

#include "FWCore/Framework/interface/Event.h"

#include "DataFormats/Common/interface/Handle.h"

#include "FWCore/ParameterSet/interface/ParameterSet.h"


//------------------------------------

// Collaborating Class Declarations --

//------------------------------------


#include "DataFormats/L1GlobalMuonTrigger/interface/L1MuRegionalCand.h"

#include "DataFormats/L1GlobalMuonTrigger/interface/L1MuGMTReadoutCollection.h"

#include "FWCore/Utilities/interface/InputTag.h"


class TFile;

class TTree;


//              ---------------------

//              -- Class Interface --

//              ---------------------


const int MAXGEN = 20;

const int MAXRPC = 12;

const int MAXDTBX = 12;

const int MAXCSC = 12;

const int MAXGMT = 12;

const int MAXGT = 12;

const int MAXRCTREG = 400;

const int MAXDTPH = 50;

const int MAXDTTH = 50;

const int MAXDTTR = 50;


class L1PromptAnalysis : public edm::EDAnalyzer {

public:

  // constructor

  explicit L1PromptAnalysis(const edm::ParameterSet&);

  virtual ~L1PromptAnalysis();


  // fill tree

  virtual void analyze(const edm::Event&, const edm::EventSetup&);

  void book();


  virtual void beginJob();

  virtual void endJob();


private:

  //GENERAL block

  int runn;

  int eventn;

  int lumi;

  int bx;

  uint64_t orbitn;

  uint64_t timest;


  // Generator info

  float weight;

  float pthat;


  // simulation block

  int ngen;

  float pxgen[MAXGEN];

  float pygen[MAXGEN];

  float pzgen[MAXGEN];

  float ptgen[MAXGEN];

  float etagen[MAXGEN];

  float phigen[MAXGEN];

  int chagen[MAXGEN];

  float vxgen[MAXGEN];

  float vygen[MAXGEN];

  float vzgen[MAXGEN];

  int pargen[MAXGEN];


  // GMT data

  int bxgmt;


  //DTBX Trigger block

  int ndt;

  int bxd[MAXDTBX];

  float ptd[MAXDTBX];

  int chad[MAXDTBX];

  float etad[MAXDTBX];

  int etafined[MAXDTBX];

  float phid[MAXDTBX];

  int quald[MAXDTBX];

  int dwd[MAXDTBX];

  int chd[MAXDTBX];


  //CSC Trigger block

  int ncsc;

  int bxc[MAXCSC];

  float ptc[MAXCSC];

  int chac[MAXCSC];

  float etac[MAXCSC];

  float phic[MAXCSC];

  int qualc[MAXCSC];

  int dwc[MAXCSC];


  //RPCb Trigger

  int nrpcb;

  int bxrb[MAXRPC];

  float ptrb[MAXRPC];

  int charb[MAXRPC];

  float etarb[MAXRPC];

  float phirb[MAXRPC];

  int qualrb[MAXRPC];

  int dwrb[MAXRPC];


  //RPCf Trigger

  int nrpcf;

  int bxrf[MAXRPC];

  float ptrf[MAXRPC];

  int charf[MAXRPC];

  float etarf[MAXRPC];

  float phirf[MAXRPC];

  int qualrf[MAXRPC];

  int dwrf[MAXRPC];


  //Global Muon Trigger

  int ngmt;

  int bxg[MAXGMT];

  float ptg[MAXGMT];

  int chag[MAXGMT];

  float etag[MAXGMT];

  float phig[MAXGMT];

  int qualg[MAXGMT];

  int detg[MAXGMT];

  int rankg[MAXGMT];

  int isolg[MAXGMT];

  int mipg[MAXGMT];

  int dwg[MAXGMT];

  int idxRPCb[MAXGMT];

  int idxRPCf[MAXGMT];

  int idxDTBX[MAXGMT];

  int idxCSC[MAXGMT];


  // GT info

  uint64_t gttw1[3];

  uint64_t gttw2[3];

  uint64_t gttt[3];


  //PSB info

  int nele;

  int bxel[MAXGT];

  float rankel[MAXGT];

  float phiel[MAXGT];

  float etael[MAXGT];


  int njet;

  int bxjet[MAXGT];

  float rankjet[MAXGT];

  float phijet[MAXGT];

  float etajet[MAXGT];


  //GCT

  edm::InputTag gctCenJetsSource_;

  edm::InputTag gctForJetsSource_;

  edm::InputTag gctTauJetsSource_;

  edm::InputTag gctEnergySumsSource_;

  edm::InputTag gctIsoEmSource_;

  edm::InputTag gctNonIsoEmSource_;


  int gctIsoEmSize;

  float gctIsoEmEta[4];

  float gctIsoEmPhi[4];

  float gctIsoEmRnk[4];

  int gctNonIsoEmSize;

  float gctNonIsoEmEta[4];

  float gctNonIsoEmPhi[4];

  float gctNonIsoEmRnk[4];

  int gctCJetSize;

  float gctCJetEta[4];

  float gctCJetPhi[4];

  float gctCJetRnk[4];

  int gctFJetSize;

  float gctFJetEta[4];

  float gctFJetPhi[4];

  float gctFJetRnk[4];

  int gctTJetSize;

  float gctTJetEta[4];

  float gctTJetPhi[4];

  float gctTJetRnk[4];

  float gctEtMiss;

  float gctEtMissPhi;

  float gctEtHad;

  float gctEtTot;

  int gctHFRingEtSumSize;

  float gctHFRingEtSumEta[4];

  float gctHFBitCountsSize;

  float gctHFBitCountsEta[4];

  //  RCT


  edm::InputTag rctSource_;

  int rctRegSize;

  float rctRegEta[MAXRCTREG];

  float rctRegPhi[MAXRCTREG];

  float rctRegRnk[MAXRCTREG];

  int rctRegVeto[MAXRCTREG];

  int rctRegBx[MAXRCTREG];

  int rctRegOverFlow[MAXRCTREG];

  int rctRegMip[MAXRCTREG];

  int rctRegFGrain[MAXRCTREG];

  int rctEmSize;

  int rctIsIsoEm[MAXRCTREG];

  float rctEmEta[MAXRCTREG];

  float rctEmPhi[MAXRCTREG];

  float rctEmRnk[MAXRCTREG];

  int rctEmBx[MAXRCTREG];


  // DTTF

  edm::InputTag dttfSource_;


  int dttf_phSize;

  int dttf_phBx[MAXDTPH];

  int dttf_phWh[MAXDTPH];

  int dttf_phSe[MAXDTPH];

  int dttf_phSt[MAXDTPH];

  float dttf_phAng[MAXDTPH];

  float dttf_phBandAng[MAXDTPH];

  int dttf_phCode[MAXDTPH];

  float dttf_phX[MAXDTPH];

  float dttf_phY[MAXDTPH];


  int dttf_thSize;

  int dttf_thBx[MAXDTTH];

  int dttf_thWh[MAXDTTH];

  int dttf_thSe[MAXDTTH];

  int dttf_thSt[MAXDTTH];

  float dttf_thX[MAXDTTH];

  float dttf_thY[MAXDTTH];

  float dttf_thTheta[MAXDTTH][7];

  int dttf_thCode[MAXDTTH][7];


  int dttf_trSize;

  int dttf_trBx[MAXDTTR];

  int dttf_trTag[MAXDTTR];

  int dttf_trQual[MAXDTTR];

  int dttf_trPtPck[MAXDTTR];

  float dttf_trPtVal[MAXDTTR];

  int dttf_trPhiPck[MAXDTTR];

  float dttf_trPhiVal[MAXDTTR];

  int dttf_trPhiGlob[MAXDTTR];

  int dttf_trChPck[MAXDTTR];

  int dttf_trWh[MAXDTTR];

  int dttf_trSc[MAXDTTR];


  TFile* m_file;

  TTree* m_tree;


  edm::InputTag m_GMTInputTag;

  edm::InputTag m_GTEvmInputTag;

  edm::InputTag m_GTInputTag;

  edm::InputTag m_GeneratorInputTag;

  edm::InputTag m_SimulationInputTag;


  bool m_PhysVal;

  bool verbose_;

  std::string m_outfilename;

};


#endif