Skip to content

Modifiy JetAnalyzer to include more processes at once #144

New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

Sign up for GitHub

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

Draft
davzuliani wants to merge 5 commits into
base: master
Choose a base branch
from
Draft

Modifiy JetAnalyzer to include more processes at once #144

Show file tree
Hide file tree
Changes from all commits
Commits
Show all changes
5 commits
Select commit Hold shift + click to select a range
686014b
Add 4quarks data
gianelle Jan 14, 2025
b0408c1
initial commit
Apr 10, 2025
fb13d4a
update JetAnalyzer tool
Apr 14, 2025
c6971b6
add cc and qq
Apr 14, 2025
b0a8130
Âfinished jetanalyzer, adding also dijet
Apr 17, 2025
File filter

Filter by extension

Filter by extension
Conversations
Failed to load comments.
Loading
Jump to
Jump to file
Failed to load files.
Loading
Diff view
Diff view
79 changes: 72 additions & 7 deletions Calorimetry/include/JetAnalyzer.h
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -54,12 +54,13 @@ class JetAnalyzer : public Processor {
protected:

std::string m_inputMCParticleCollection="";
std::string m_inputRECOParticleCollection="";
std::string m_genjetColName="";
std::string m_recojetColName="";
// std::string m_inputRECOParticleCollection="";
// std::string m_genjetColName="";
// std::string m_recojetColName="";


std::string m_rootFileName="";
std::string m_processName="";

// Run and event counters
int m_eventNumber=0;
Expand All @@ -78,6 +79,70 @@ class JetAnalyzer : public Processor {
float m_d2_mcPz=0;
//particle level information sums --> visible and invisible energy

//#######NEW#####################
int m_H1_mcPDGID = 0;
float m_H1_mcE = 0;
float m_H1_mcPx = 0;
float m_H1_mcPy = 0;
float m_H1_mcPz = 0;

int m_H2_mcPDGID = 0;
float m_H2_mcE = 0;
float m_H2_mcPx = 0;
float m_H2_mcPy = 0;
float m_H2_mcPz = 0;

int m_b1_mcPDGID = 0;
float m_b1_mcE = 0;
float m_b1_mcPx = 0;
float m_b1_mcPy = 0;
float m_b1_mcPz = 0;

int m_b2_mcPDGID = 0;
float m_b2_mcE = 0;
float m_b2_mcPx = 0;
float m_b2_mcPy = 0;
float m_b2_mcPz = 0;

int m_b3_mcPDGID = 0;
float m_b3_mcE = 0;
float m_b3_mcPx = 0;
float m_b3_mcPy = 0;
float m_b3_mcPz = 0;

int m_b4_mcPDGID = 0;
float m_b4_mcE = 0;
float m_b4_mcPx = 0;
float m_b4_mcPy = 0;
float m_b4_mcPz = 0;

int m_b5_mcPDGID = 0;
float m_b5_mcE = 0;
float m_b5_mcPx = 0;
float m_b5_mcPy = 0;
float m_b5_mcPz = 0;

int m_b6_mcPDGID = 0;
float m_b6_mcE = 0;
float m_b6_mcPx = 0;
float m_b6_mcPy = 0;
float m_b6_mcPz = 0;

int m_b7_mcPDGID = 0;
float m_b7_mcE = 0;
float m_b7_mcPx = 0;
float m_b7_mcPy = 0;
float m_b7_mcPz = 0;

int m_b8_mcPDGID = 0;
float m_b8_mcE = 0;
float m_b8_mcPx = 0;
float m_b8_mcPy = 0;
float m_b8_mcPz = 0;


//#######ENDNEW#####################

//only neutrinos
float m_E_trueInv=0;
float m_px_trueInv=0;
Expand All @@ -90,20 +155,20 @@ class JetAnalyzer : public Processor {
float m_py_trueVis=0;
float m_pz_trueVis=0;

//all reconstructed particles
/* //all reconstructed particles
float m_E_totPFO=0;
float m_px_totPFO=0;
float m_py_totPFO=0;
float m_pz_totPFO=0;

*/
//if extended parton quantities are saved: dibosons and their decay products
std::vector<float>* m_trueME_E=NULL;
std::vector<float>* m_trueME_Px=NULL;
std::vector<float>* m_trueME_Py=NULL;
std::vector<float>* m_trueME_Pz=NULL;
std::vector<int>* m_trueME_PDGID=NULL;

std::vector<float>* m_genJet_E=NULL;
/* std::vector<float>* m_genJet_E=NULL;
std::vector<float>* m_genJet_Px=NULL;
std::vector<float>* m_genJet_Py=NULL;
std::vector<float>* m_genJet_Pz=NULL;
Expand All @@ -112,7 +177,7 @@ class JetAnalyzer : public Processor {
std::vector<float>* m_recoJet_Px=NULL;
std::vector<float>* m_recoJet_Py=NULL;
std::vector<float>* m_recoJet_Pz=NULL;

*/

// Call to get collections
void getCollection(LCCollection*&, std::string, LCEvent*);
Expand Down
Loading