File indexing completed on 2024-09-27 07:03:27
0001
0002
0003 #ifndef EXCLUDE_DELPHES
0004
0005
0006
0007
0008
0009 #pragma once
0010
0011
0012 #include "Kinematics.h"
0013
0014
0015 #include <classes/DelphesClasses.h>
0016 #include <fastjet/ClusterSequence.hh>
0017 #ifdef INCLUDE_CENTAURO
0018 #include <fastjet/plugins/Centauro/Centauro.hh>
0019 #endif
0020
0021 using std::map;
0022 using std::cout;
0023 using std::cerr;
0024 using std::endl;
0025
0026 class KinematicsJets : public Kinematics
0027 {
0028 public:
0029 KinematicsJets(Double_t enEleBeam, Double_t enIonBeam, Double_t crossAng);
0030 ~KinematicsJets();
0031
0032
0033
0034
0035 void GetJets(
0036 TObjArrayIter itEFlowTrack, TObjArrayIter itEFlowPhoton,
0037 TObjArrayIter itEFlowNeutralHadron, TObjArrayIter itParticle,
0038 int jetAlgo, double jetRadius, double jetMinPt
0039 );
0040 void CalculateJetKinematics(fastjet::PseudoJet jet);
0041 void CalculateJetResolution(double deltaRCut);
0042 #ifdef INCLUDE_CENTAURO
0043 void GetBreitFrameJets(
0044 TObjArrayIter itEFlowTrack, TObjArrayIter itEFlowPhoton,
0045 TObjArrayIter itEFlowNeutralHadron, TObjArrayIter itParticle
0046 );
0047 void CalculateBreitJetKinematics(fastjet::PseudoJet jet);
0048 #endif
0049
0050
0051
0052
0053 int jetAlgo;
0054 double jetRad, jetMinPt;
0055
0056 std::vector<fastjet::PseudoJet> jetsRec, jetsTrue;
0057 std::vector<fastjet::PseudoJet> breitJetsRec, breitJetsTrue;
0058 std::map<double, int> jetConstituents;
0059 fastjet::ClusterSequence csRec;
0060 fastjet::ClusterSequence csTrue;
0061
0062 Double_t zjet, pTjet, qTjet, mTjet, etajet, phijet, mjet, ejet;
0063 Double_t deltaRjet;
0064 int matchStatusjet;
0065 Double_t pTmtjet, mTmtjet, etamtjet, phimtjet, mmtjet, emtjet;
0066 std::vector<double> jperp;
0067 std::vector<double> zhad_jet;
0068
0069 Double_t quarkpT;
0070
0071 ClassDef(KinematicsJets,1);
0072 };
0073
0074 #endif