fConeRadius = 0.3;
fDebug = 0;
fOutput=0;
- hFragmFcn=0;// = new TH1F("hFragmFcn","Fragmentation Function",100,0,1);
- hPartonFragmFcn=0;// = new TH1F("hPartonFragmFcn","Fragmentation Function",100,0,1);
- hPartonJT=0;// = new TH1F("hPartonJT","Track Momentum Perpendicular to Parton Axis",100,0.,10.);
- hPartonPL=0;// = new TH1F("hPartonPL","Track Momentum Parallel to Parton Axis ",100,0.,100.);
- hJetJT=0;// = new TH1F("hJetJT","Track Momentum Perpendicular to Jet Axis",100,0.,10.);
- hJetPL=0;// = new TH1F("hJetPL","Track Momentum Parallel to Jet Axis ",100,0.,100.);
- hJetEt=0;// = new TH1F("hJetEt","E_{T}^{reco}",250,0.,250.);
- hJetEta=0;// = new TH1F("hJetEta","#eta_{jet}^{reco}",180,-0.9,0.9);
- hJetPhi=0;// = new TH1F("hJetPhi","#phi_{jet}^{reco}",62,0.,3.1);
- hPartonEta=0;// = new TH1F("hPartonEta","#eta_{Parton}",180,-0.9,0.9);
- hPartonPhi=0;// = new TH1F("hPartonPhi","#phi_{Parton}",62,0.,3.1);
- hEtaDiff=0;// = new TH1F("hEtaDiff","#eta_{jet}^{reco}-#eta_{jet}^{input}",100,-0.5,0.5);
- hPhiDiff=0;// = new TH1F("hPhiDiff","#phi_{jet}^{reco}-#phi_{jet}^{input}",100,-0.5,0.5);
- hNJets=0;// = new TH1F("hNJets","N Reconstructed jets",11,-0.5,10.5);
- hEtaPhiSpread=0;
-
-hFragmFcn2=0; // ("hFragmFcn2","Fragmentation Function",100,0,1);
-hPartonFragmFcn2=0;// ("hFragmFcn2","Parton Fragmentation Function",100,0,1);
-hPartonJT2=0; // ("hPartonJT2","Track Momentum Perpendicular to Parton Axis",100,0.,10.);
-hPartonPL2=0; // ("hPartonPL2","Track Momentum Parallel to Parton Axis ",100,0.,100.);
-hJetJT2=0; // ("hJetJT2","Track Momentum Perpendicular to Jet Axis",100,0.,10.);
-hJetPL2=0; // ("hJetPL2","Track Momentum Parallel to Jet Axis ",100,0.,100.);
-hJetEt2=0; // ("hJetEt2","E_{T}^{reco}",250,0.,250.);
-hJetEta2=0; // ("hJetEta2","#eta_{jet}^{reco}",180,-0.9,0.9);
-hJetPhi2=0; // ("hJetPhi2","#phi_{jet}^{reco}",62,0.,3.1);
-hPartonEta2=0; // ("hPartonEta2","#eta_{Parton}",180,-0.9,0.9);
-hPartonPhi2=0; // ("hPartonPhi2","#phi_{Parton}",62,0.,3.1);
-hEtaDiff2=0; // ("hEtaDiff2","#eta_{jet}^{reco}-#eta_{jet}^{input}",100,-0.5,0.5);
-hPhiDiff2=0; // ("hPhiDiff2","#phi_{jet}^{reco}-#phi_{jet}^{input}",100,-0.5,0.5);
-hEtaPhiSpread2=0; // ("hEtaPhiSpread2","#eta - #phi Distribution
+ fhFragmFcn=0;// = new TH1F("hFragmFcn","Fragmentation Function",100,0,1);
+ fhPartonFragmFcn=0;// = new TH1F("hPartonFragmFcn","Fragmentation Function",100,0,1);
+ fhPartonJT=0;// = new TH1F("hPartonJT","Track Momentum Perpendicular to Parton Axis",100,0.,10.);
+ fhPartonPL=0;// = new TH1F("hPartonPL","Track Momentum Parallel to Parton Axis ",100,0.,100.);
+ fhJetJT=0;// = new TH1F("hJetJT","Track Momentum Perpendicular to Jet Axis",100,0.,10.);
+ fhJetPL=0;// = new TH1F("hJetPL","Track Momentum Parallel to Jet Axis ",100,0.,100.);
+ fhJetEt=0;// = new TH1F("hJetEt","E_{T}^{reco}",250,0.,250.);
+ fhJetEta=0;// = new TH1F("hJetEta","#eta_{jet}^{reco}",180,-0.9,0.9);
+ fhJetPhi=0;// = new TH1F("hJetPhi","#phi_{jet}^{reco}",62,0.,3.1);
+ fhPartonEta=0;// = new TH1F("hPartonEta","#eta_{Parton}",180,-0.9,0.9);
+ fhPartonPhi=0;// = new TH1F("hPartonPhi","#phi_{Parton}",62,0.,3.1);
+ fhEtaDiff=0;// = new TH1F("hEtaDiff","#eta_{jet}^{reco}-#eta_{jet}^{input}",100,-0.5,0.5);
+ fhPhiDiff=0;// = new TH1F("hPhiDiff","#phi_{jet}^{reco}-#phi_{jet}^{input}",100,-0.5,0.5);
+ fhNJets=0;// = new TH1F("hNJets","N Reconstructed jets",11,-0.5,10.5);
+ fhEtaPhiSpread=0;
+
+fhFragmFcn2=0; // ("hFragmFcn2","Fragmentation Function",100,0,1);
+fhPartonFragmFcn2=0;// ("hFragmFcn2","Parton Fragmentation Function",100,0,1);
+fhPartonJT2=0; // ("hPartonJT2","Track Momentum Perpendicular to Parton Axis",100,0.,10.);
+fhPartonPL2=0; // ("hPartonPL2","Track Momentum Parallel to Parton Axis ",100,0.,100.);
+fhJetJT2=0; // ("hJetJT2","Track Momentum Perpendicular to Jet Axis",100,0.,10.);
+fhJetPL2=0; // ("hJetPL2","Track Momentum Parallel to Jet Axis ",100,0.,100.);
+fhJetEt2=0; // ("hJetEt2","E_{T}^{reco}",250,0.,250.);
+fhJetEta2=0; // ("hJetEta2","#eta_{jet}^{reco}",180,-0.9,0.9);
+fhJetPhi2=0; // ("hJetPhi2","#phi_{jet}^{reco}",62,0.,3.1);
+fhPartonEta2=0; // ("hPartonEta2","#eta_{Parton}",180,-0.9,0.9);
+fhPartonPhi2=0; // ("hPartonPhi2","#phi_{Parton}",62,0.,3.1);
+fhEtaDiff2=0; // ("hEtaDiff2","#eta_{jet}^{reco}-#eta_{jet}^{input}",100,-0.5,0.5);
+fhPhiDiff2=0; // ("hPhiDiff2","#phi_{jet}^{reco}-#phi_{jet}^{input}",100,-0.5,0.5);
+fhEtaPhiSpread2=0; // ("hEtaPhiSpread2","#eta - #phi Distribution
//of Reconstructed Jets",192,-0.7,0.7,288,pi/3,pi);
-hNJets2=0; // ("hNJets2","N Reconstructed jets",11,-0.5,10.5);
-hJetEtSecond2=0; //("hJetEtSecond2","E_{T}^{reco}",250,0.,250.);
-hJetEtRatio2=0; //("hJetEtRatio2","Ratio of Second Highest to Highest",100,0,1);
-hEtaPhiDist2=0; //("hEtaPhiDist2","Angular Distance Between First and Second",100,0,3);
+fhNJets2=0; // ("hNJets2","N Reconstructed jets",11,-0.5,10.5);
+fhJetEtSecond2=0; //("hJetEtSecond2","E_{T}^{reco}",250,0.,250.);
+fhJetEtRatio2=0; //("hJetEtRatio2","Ratio of Second Highest to Highest",100,0,1);
+fhEtaPhiDist2=0; //("hEtaPhiDist2","Angular Distance Between First and Second",100,0,3);
}
void AliEMCALJetFinderPlots::InitPlots()
{
//========================= CASE 1 =======================================
- hFragmFcn = new TH1F("hFragmFcn","Fragmentation Function",100,0,1);
- hFragmFcn->Sumw2();
- hPartonFragmFcn = new TH1F("hPartonFragmFcn","Parton Fragmentation Function",100,0,1);
- hPartonFragmFcn->Sumw2();
- hPartonJT = new TH1F("hPartonJT","Track Momentum Perpendicular to Parton Axis",100,0.,10.);
- hPartonJT->Sumw2();
- hPartonPL = new TH1F("hPartonPL","Track Momentum Parallel to Parton Axis ",100,0.,100.);
- hPartonPL->Sumw2();
- hJetJT = new TH1F("hJetJT","Track Momentum Perpendicular to Jet Axis",100,0.,10.);
- hJetJT->Sumw2();
- hJetPL = new TH1F("hJetPL","Track Momentum Parallel to Jet Axis ",100,0.,100.);
- hJetPL->Sumw2();
- hJetEt = new TH1F("hJetEt","E_{T}^{reco}",250,0.,250.);
- hJetEt->Sumw2();
- hJetEta = new TH1F("hJetEta","#eta_{jet}^{reco}",180,-0.9,0.9);
- hJetEta->Sumw2();
- hJetPhi = new TH1F("hJetPhi","#phi_{jet}^{reco}",62,0.,3.1);
- hJetPhi->Sumw2();
- hPartonEta = new TH1F("hPartonEta","#eta_{Parton}",180,-0.9,0.9);
- hPartonEta->Sumw2();
- hPartonPhi = new TH1F("hPartonPhi","#phi_{Parton}",62,0.,3.1);
- hPartonPhi->Sumw2();
- hEtaDiff = new TH1F("hEtaDiff","#eta_{jet}^{reco}-#eta_{jet}^{input}",100,-0.5,0.5);
- hEtaDiff->Sumw2();
- hPhiDiff = new TH1F("hPhiDiff","#phi_{jet}^{reco}-#phi_{jet}^{input}",100,-0.5,0.5);
- hPhiDiff->Sumw2();
- hNJets = new TH1F("hNJets","N Reconstructed jets",11,-0.5,10.5);
- hNJets->Sumw2();
- hEtaPhiSpread = new TH2F("hEtaPhiSpread","#eta - #phi Distribution of Reconstructed Jets",100,-0.5,0.5,100,-0.5,0.5);
- hEtaPhiSpread->Sumw2();
- hNJets->SetXTitle("N_{jets}^{reco}/event");
- hNJets->SetYTitle("N_{events}");
+ fhFragmFcn = new TH1F("hFragmFcn","Fragmentation Function",100,0,1);
+ fhFragmFcn->Sumw2();
+ fhPartonFragmFcn = new TH1F("hPartonFragmFcn","Parton Fragmentation Function",100,0,1);
+ fhPartonFragmFcn->Sumw2();
+ fhPartonJT = new TH1F("hPartonJT","Track Momentum Perpendicular to Parton Axis",100,0.,10.);
+ fhPartonJT->Sumw2();
+ fhPartonPL = new TH1F("hPartonPL","Track Momentum Parallel to Parton Axis ",100,0.,100.);
+ fhPartonPL->Sumw2();
+ fhJetJT = new TH1F("hJetJT","Track Momentum Perpendicular to Jet Axis",100,0.,10.);
+ fhJetJT->Sumw2();
+ fhJetPL = new TH1F("hJetPL","Track Momentum Parallel to Jet Axis ",100,0.,100.);
+ fhJetPL->Sumw2();
+ fhJetEt = new TH1F("hJetEt","E_{T}^{reco}",250,0.,250.);
+ fhJetEt->Sumw2();
+ fhJetEta = new TH1F("hJetEta","#eta_{jet}^{reco}",180,-0.9,0.9);
+ fhJetEta->Sumw2();
+ fhJetPhi = new TH1F("hJetPhi","#phi_{jet}^{reco}",62,0.,3.1);
+ fhJetPhi->Sumw2();
+ fhPartonEta = new TH1F("hPartonEta","#eta_{Parton}",180,-0.9,0.9);
+ fhPartonEta->Sumw2();
+ fhPartonPhi = new TH1F("hPartonPhi","#phi_{Parton}",62,0.,3.1);
+ fhPartonPhi->Sumw2();
+ fhEtaDiff = new TH1F("hEtaDiff","#eta_{jet}^{reco}-#eta_{jet}^{input}",100,-0.5,0.5);
+ fhEtaDiff->Sumw2();
+ fhPhiDiff = new TH1F("hPhiDiff","#phi_{jet}^{reco}-#phi_{jet}^{input}",100,-0.5,0.5);
+ fhPhiDiff->Sumw2();
+ fhNJets = new TH1F("hNJets","N Reconstructed jets",11,-0.5,10.5);
+ fhNJets->Sumw2();
+ fhEtaPhiSpread = new TH2F("hEtaPhiSpread","#eta - #phi Distribution of Reconstructed Jets",100,-0.5,0.5,100,-0.5,0.5);
+ fhEtaPhiSpread->Sumw2();
+ fhNJets->SetXTitle("N_{jets}^{reco}/event");
+ fhNJets->SetYTitle("N_{events}");
//Jet properties
- hJetEt->SetFillColor(16);
- hJetEt->SetXTitle("E_{T}^{reco}");
+ fhJetEt->SetFillColor(16);
+ fhJetEt->SetXTitle("E_{T}^{reco}");
- hJetEta->SetFillColor(16);
- hJetEta->SetXTitle("#eta_{jet}^{reco}");
+ fhJetEta->SetFillColor(16);
+ fhJetEta->SetXTitle("#eta_{jet}^{reco}");
- hJetPhi->SetFillColor(16);
- hJetPhi->SetXTitle("#phi_{jet}^{reco}");
+ fhJetPhi->SetFillColor(16);
+ fhJetPhi->SetXTitle("#phi_{jet}^{reco}");
- hPartonEta->SetFillColor(16);
- hPartonEta->SetXTitle("#eta_{parton}");
+ fhPartonEta->SetFillColor(16);
+ fhPartonEta->SetXTitle("#eta_{parton}");
- hPartonPhi->SetFillColor(16);
- hPartonPhi->SetXTitle("#phi_{parton}");
+ fhPartonPhi->SetFillColor(16);
+ fhPartonPhi->SetXTitle("#phi_{parton}");
- hPartonPL->SetXTitle("p (GeV/c)");
- hPartonJT->SetXTitle("p (GeV/c)");
+ fhPartonPL->SetXTitle("p (GeV/c)");
+ fhPartonJT->SetXTitle("p (GeV/c)");
- hPartonFragmFcn->SetXTitle("Z = p_{T}^{Chg}/E_{T}^{parton}");
+ fhPartonFragmFcn->SetXTitle("Z = p_{T}^{Chg}/E_{T}^{parton}");
//Jet component properties
- hJetPL->SetXTitle("p (GeV/c)");
- hJetJT->SetXTitle("p (GeV/c)");
- hFragmFcn->SetXTitle("Z = p_{T}^{Chg}/E_{T}^{reco}");
- hPartonFragmFcn->SetXTitle("Z = p_{T}^{Chg}/E_{T}^{reco}");
+ fhJetPL->SetXTitle("p (GeV/c)");
+ fhJetJT->SetXTitle("p (GeV/c)");
+ fhFragmFcn->SetXTitle("Z = p_{T}^{Chg}/E_{T}^{reco}");
+ fhPartonFragmFcn->SetXTitle("Z = p_{T}^{Chg}/E_{T}^{reco}");
- hEtaDiff->SetXTitle("#eta_{jet}^{reco}-#eta_{jet}^{input}");
- hPhiDiff->SetXTitle("#phi_{jet}^{reco}-#phi_{jet}^{input}");
- hEtaPhiSpread->SetXTitle("#eta");
- hEtaPhiSpread->SetYTitle("#phi");
+ fhEtaDiff->SetXTitle("#eta_{jet}^{reco}-#eta_{jet}^{input}");
+ fhPhiDiff->SetXTitle("#phi_{jet}^{reco}-#phi_{jet}^{input}");
+ fhEtaPhiSpread->SetXTitle("#eta");
+ fhEtaPhiSpread->SetYTitle("#phi");
//======================= CASE 2 ======================================
-hFragmFcn2 = new TH1F("hFragmFcn2","Fragmentation Function",100,0,1);
-hFragmFcn2->Sumw2();
-hPartonFragmFcn2 = new TH1F("hPartonFragmFcn2","Parton Fragmentation Function",100,0,1);
-hPartonFragmFcn2->Sumw2();
-hPartonJT2 = new TH1F("hPartonJT2","Track Momentum Perpendicular to Parton Axis",100,0.,10.);
-hPartonJT2->Sumw2();
-hPartonPL2 = new TH1F("hPartonPL2","Track Momentum Parallel to Parton Axis ",100,0.,100.);
-hPartonPL2->Sumw2();
-hJetJT2 = new TH1F("hJetJT2","Track Momentum Perpendicular to Jet Axis",100,0.,10.);
-hJetJT2->Sumw2();
-hJetPL2 = new TH1F("hJetPL2","Track Momentum Parallel to Jet Axis ",100,0.,100.);
-hJetPL2->Sumw2();
-hJetEt2 = new TH1F("hJetEt2","E_{T}^{reco}",250,0.,250.);
-hJetEt2->Sumw2();
-hJetEta2 = new TH1F("hJetEta2","#eta_{jet}^{reco}",180,-0.9,0.9);
-hJetEta2->Sumw2();
-hJetPhi2 = new TH1F("hJetPhi2","#phi_{jet}^{reco}",62,0.,3.1);
-hJetPhi2->Sumw2();
-hPartonEta2 = new TH1F("hPartonEta2","#eta_{Parton}",180,-0.9,0.9);
-hPartonEta2->Sumw2();
-hPartonPhi2 = new TH1F("hPartonPhi2","#phi_{Parton}",62,0.,3.1);
-hPartonPhi2->Sumw2();
-hEtaDiff2 = new TH1F("hEtaDiff2","#eta_{jet}^{reco}-#eta_{jet}^{input}",100,-0.5,0.5);
-hEtaDiff2->Sumw2();
-hPhiDiff2 = new TH1F("hPhiDiff2","#phi_{jet}^{reco}-#phi_{jet}^{input}",100,-0.5,0.5);
-hPhiDiff2->Sumw2();
-hEtaPhiSpread2 = new TH2F("hEtaPhiSpread2","#eta - #phi Distribution of Reconstructed Jets",100,-0.5,0.5,100,-0.5,0.5);
-hEtaPhiSpread2->Sumw2();
-hNJets2 = new TH1F("hNJets2","N Reconstructed jets",11,-0.5,10.5);
-hNJets2->Sumw2();
-hJetEtSecond2 = new TH1F("hJetEtSecond2","E_{T}^{reco}",250,0.,250.);
-hJetEtSecond2->Sumw2();
-hJetEtRatio2 = new TH1F("hJetEtRatio2","Ratio of Second Highest to Highest",100,0,1);
-hJetEtRatio2->Sumw2();
-hEtaPhiDist2 = new TH1F("hEtaPhiDist2","Angular Distance Between First and Second",100,0,3);
-hEtaPhiDist2->Sumw2();
+fhFragmFcn2 = new TH1F("hFragmFcn2","Fragmentation Function",100,0,1);
+fhFragmFcn2->Sumw2();
+fhPartonFragmFcn2 = new TH1F("hPartonFragmFcn2","Parton Fragmentation Function",100,0,1);
+fhPartonFragmFcn2->Sumw2();
+fhPartonJT2 = new TH1F("hPartonJT2","Track Momentum Perpendicular to Parton Axis",100,0.,10.);
+fhPartonJT2->Sumw2();
+fhPartonPL2 = new TH1F("hPartonPL2","Track Momentum Parallel to Parton Axis ",100,0.,100.);
+fhPartonPL2->Sumw2();
+fhJetJT2 = new TH1F("hJetJT2","Track Momentum Perpendicular to Jet Axis",100,0.,10.);
+fhJetJT2->Sumw2();
+fhJetPL2 = new TH1F("hJetPL2","Track Momentum Parallel to Jet Axis ",100,0.,100.);
+fhJetPL2->Sumw2();
+fhJetEt2 = new TH1F("hJetEt2","E_{T}^{reco}",250,0.,250.);
+fhJetEt2->Sumw2();
+fhJetEta2 = new TH1F("hJetEta2","#eta_{jet}^{reco}",180,-0.9,0.9);
+fhJetEta2->Sumw2();
+fhJetPhi2 = new TH1F("hJetPhi2","#phi_{jet}^{reco}",62,0.,3.1);
+fhJetPhi2->Sumw2();
+fhPartonEta2 = new TH1F("hPartonEta2","#eta_{Parton}",180,-0.9,0.9);
+fhPartonEta2->Sumw2();
+fhPartonPhi2 = new TH1F("hPartonPhi2","#phi_{Parton}",62,0.,3.1);
+fhPartonPhi2->Sumw2();
+fhEtaDiff2 = new TH1F("hEtaDiff2","#eta_{jet}^{reco}-#eta_{jet}^{input}",100,-0.5,0.5);
+fhEtaDiff2->Sumw2();
+fhPhiDiff2 = new TH1F("hPhiDiff2","#phi_{jet}^{reco}-#phi_{jet}^{input}",100,-0.5,0.5);
+fhPhiDiff2->Sumw2();
+fhEtaPhiSpread2 = new TH2F("hEtaPhiSpread2","#eta - #phi Distribution of Reconstructed Jets",100,-0.5,0.5,100,-0.5,0.5);
+fhEtaPhiSpread2->Sumw2();
+fhNJets2 = new TH1F("hNJets2","N Reconstructed jets",11,-0.5,10.5);
+fhNJets2->Sumw2();
+fhJetEtSecond2 = new TH1F("hJetEtSecond2","E_{T}^{reco}",250,0.,250.);
+fhJetEtSecond2->Sumw2();
+fhJetEtRatio2 = new TH1F("hJetEtRatio2","Ratio of Second Highest to Highest",100,0,1);
+fhJetEtRatio2->Sumw2();
+fhEtaPhiDist2 = new TH1F("hEtaPhiDist2","Angular Distance Between First and Second",100,0,3);
+fhEtaPhiDist2->Sumw2();
fInitialised = kTRUE;
AliEMCALJetFinderPlots::~AliEMCALJetFinderPlots(){
-delete hFragmFcn;// = new TH1F("hFragmFcn","Fragmentation Function",100,0,1);
-delete hPartonFragmFcn;// = new TH1F("hFragmFcn","Fragmentation Function",100,0,1);
-delete hPartonJT;// = new TH1F("hPartonJT","Track Momentum Perpendicular to Parton Axis",100,0.,10.);
-delete hPartonPL;// = new TH1F("hPartonPL","Track Momentum Parallel to Parton Axis ",100,0.,100.);
-delete hJetJT;// = new TH1F("hJetJT","Track Momentum Perpendicular to Jet Axis",100,0.,10.);
-delete hJetPL;// = new TH1F("hJetPL","Track Momentum Parallel to Jet Axis ",100,0.,100.);
-delete hJetEt;// = new TH1F("hJetEt","E_{T}^{reco}",250,0.,250.);
-delete hJetEta;// = new TH1F("hJetEta","#eta_{jet}^{reco}",180,-0.9,0.9);
-delete hJetPhi;// = new TH1F("hJetPhi","#phi_{jet}^{reco}",62,0.,3.1);
-delete hPartonEta;// = new TH1F("hPartonEta","#eta_{Parton}",180,-0.9,0.9);
-delete hPartonPhi;// = new TH1F("hPartonPhi","#phi_{Parton}",62,0.,3.1);
-delete hEtaDiff;// = new TH1F("hEtaDiff","#eta_{jet}^{reco}-#eta_{jet}^{input}",100,-0.5,0.5);
-delete hPhiDiff;// = new TH1F("hPhiDiff","#phi_{jet}^{reco}-#phi_{jet}^{input}",100,-0.5,0.5);
-delete hNJets;// = new TH1F("hNJets","N Reconstructed jets",11,-0.5,10.5);
-delete hEtaPhiSpread;
-
- delete hFragmFcn2; // ("hFragmFcn2","Fragmentation Function",100,0,1);
- delete hPartonFragmFcn2;// ("hFragmFcn2","Parton Fragmentation Function",100,0,1);
- delete hPartonJT2; // ("hPartonJT2","Track Momentum Perpendicular to Parton Axis",100,0.,10.);
- delete hPartonPL2; // ("hPartonPL2","Track Momentum Parallel to Parton Axis ",100,0.,100.);
- delete hJetJT2; // ("hJetJT2","Track Momentum Perpendicular to Jet Axis",100,0.,10.);
- delete hJetPL2; // ("hJetPL2","Track Momentum Parallel to Jet Axis ",100,0.,100.);
- delete hJetEt2; // ("hJetEt2","E_{T}^{reco}",250,0.,250.);
- delete hJetEta2; // ("hJetEta2","#eta_{jet}^{reco}",180,-0.9,0.9);
- delete hJetPhi2; // ("hJetPhi2","#phi_{jet}^{reco}",62,0.,3.1);
- delete hPartonEta2; // ("hPartonEta2","#eta_{Parton}",180,-0.9,0.9);
- delete hPartonPhi2; // ("hPartonPhi2","#phi_{Parton}",62,0.,3.1);
- delete hEtaDiff2; // ("hEtaDiff2","#eta_{jet}^{reco}-#eta_{jet}^{input}",100,-0.5,0.5);
- delete hPhiDiff2; // ("hPhiDiff2","#phi_{jet}^{reco}-#phi_{jet}^{input}",100,-0.5,0.5);
- delete hEtaPhiSpread2; // ("hEtaPhiSpread2","#eta - #phi Distribution
+delete fhFragmFcn;// = new TH1F("hFragmFcn","Fragmentation Function",100,0,1);
+delete fhPartonFragmFcn;// = new TH1F("hFragmFcn","Fragmentation Function",100,0,1);
+delete fhPartonJT;// = new TH1F("hPartonJT","Track Momentum Perpendicular to Parton Axis",100,0.,10.);
+delete fhPartonPL;// = new TH1F("hPartonPL","Track Momentum Parallel to Parton Axis ",100,0.,100.);
+delete fhJetJT;// = new TH1F("hJetJT","Track Momentum Perpendicular to Jet Axis",100,0.,10.);
+delete fhJetPL;// = new TH1F("hJetPL","Track Momentum Parallel to Jet Axis ",100,0.,100.);
+delete fhJetEt;// = new TH1F("hJetEt","E_{T}^{reco}",250,0.,250.);
+delete fhJetEta;// = new TH1F("hJetEta","#eta_{jet}^{reco}",180,-0.9,0.9);
+delete fhJetPhi;// = new TH1F("hJetPhi","#phi_{jet}^{reco}",62,0.,3.1);
+delete fhPartonEta;// = new TH1F("hPartonEta","#eta_{Parton}",180,-0.9,0.9);
+delete fhPartonPhi;// = new TH1F("hPartonPhi","#phi_{Parton}",62,0.,3.1);
+delete fhEtaDiff;// = new TH1F("hEtaDiff","#eta_{jet}^{reco}-#eta_{jet}^{input}",100,-0.5,0.5);
+delete fhPhiDiff;// = new TH1F("hPhiDiff","#phi_{jet}^{reco}-#phi_{jet}^{input}",100,-0.5,0.5);
+delete fhNJets;// = new TH1F("hNJets","N Reconstructed jets",11,-0.5,10.5);
+delete fhEtaPhiSpread;
+
+ delete fhFragmFcn2; // ("hFragmFcn2","Fragmentation Function",100,0,1);
+ delete fhPartonFragmFcn2;// ("hFragmFcn2","Parton Fragmentation Function",100,0,1);
+ delete fhPartonJT2; // ("hPartonJT2","Track Momentum Perpendicular to Parton Axis",100,0.,10.);
+ delete fhPartonPL2; // ("hPartonPL2","Track Momentum Parallel to Parton Axis ",100,0.,100.);
+ delete fhJetJT2; // ("hJetJT2","Track Momentum Perpendicular to Jet Axis",100,0.,10.);
+ delete fhJetPL2; // ("hJetPL2","Track Momentum Parallel to Jet Axis ",100,0.,100.);
+ delete fhJetEt2; // ("hJetEt2","E_{T}^{reco}",250,0.,250.);
+ delete fhJetEta2; // ("hJetEta2","#eta_{jet}^{reco}",180,-0.9,0.9);
+ delete fhJetPhi2; // ("hJetPhi2","#phi_{jet}^{reco}",62,0.,3.1);
+ delete fhPartonEta2; // ("hPartonEta2","#eta_{Parton}",180,-0.9,0.9);
+ delete fhPartonPhi2; // ("hPartonPhi2","#phi_{Parton}",62,0.,3.1);
+ delete fhEtaDiff2; // ("hEtaDiff2","#eta_{jet}^{reco}-#eta_{jet}^{input}",100,-0.5,0.5);
+ delete fhPhiDiff2; // ("hPhiDiff2","#phi_{jet}^{reco}-#phi_{jet}^{input}",100,-0.5,0.5);
+ delete fhEtaPhiSpread2; // ("hEtaPhiSpread2","#eta - #phi Distribution
//of Reconstructed Jets",192,-0.7,0.7,288,pi/3,pi);
- delete hNJets2; // ("hNJets2","N Reconstructed jets",11,-0.5,10.5);
- delete hJetEtSecond2; //("hJetEtSecond2","E_{T}^{reco}",250,0.,250.);
- delete hJetEtRatio2; //("hJetEtRatio2","Ratio of Second Highest to Highest",100,0,1);
- delete hEtaPhiDist2; //("hEtaPhiDist2","Angular Distance Between First and Second",100,0,3);
+ delete fhNJets2; // ("hNJets2","N Reconstructed jets",11,-0.5,10.5);
+ delete fhJetEtSecond2; //("hJetEtSecond2","E_{T}^{reco}",250,0.,250.);
+ delete fhJetEtRatio2; //("hJetEtRatio2","Ratio of Second Highest to Highest",100,0,1);
+ delete fhEtaPhiDist2; //("hEtaPhiDist2","Angular Distance Between First and Second",100,0,3);
if (!fInitialised) InitPlots();
fOutput = output;
if (!fOutput) return;
- hNJets->Fill(fOutput->GetNJets());
+ fhNJets->Fill(fOutput->GetNJets());
if (fOutput->GetNJets()>1)
{
//========================= CASE 2 ===========================
Int_t nPartons = fOutput->GetNPartons();
- hNJets2->Fill(fOutput->GetNJets());
+ fhNJets2->Fill(fOutput->GetNJets());
AliEMCALParton* parton;
AliEMCALJet* jethighest=0;
AliEMCALJet* jetsecond=0;
}
// End finding highest and second highest and continue
- hJetEt2->Fill(jethighest->Energy());
- hJetEta2->Fill(jethighest->Eta() );
- hJetPhi2->Fill(jethighest->Phi() );
+ fhJetEt2->Fill(jethighest->Energy());
+ fhJetEta2->Fill(jethighest->Eta() );
+ fhJetPhi2->Fill(jethighest->Phi() );
if (nPartons ==0) return;
parton = fOutput->GetParton(0);
- hPartonEta2->Fill( parton->Eta() );
- hPartonPhi2->Fill( parton->Phi() );
+ fhPartonEta2->Fill( parton->Eta() );
+ fhPartonPhi2->Fill( parton->Phi() );
//hJetEtDiff->Fill( jet->Energy() - parton->Energy() );
- hEtaDiff2->Fill( jethighest->Eta() - parton->Eta() );
- hPhiDiff2->Fill( jethighest->Phi() - parton->Phi() );
- hEtaPhiSpread2->Fill(jethighest->Eta()-parton->Eta(),jethighest->Phi() - parton->Phi());
- hJetEtSecond2->Fill(jetsecond->Energy());
- hJetEtRatio2->Fill(jetsecond->Energy()/jethighest->Energy());
- hEtaPhiDist2->Fill( TMath::Sqrt((jethighest->Eta() - jetsecond->Eta())*(jethighest->Eta() - jetsecond->Eta())
+ fhEtaDiff2->Fill( jethighest->Eta() - parton->Eta() );
+ fhPhiDiff2->Fill( jethighest->Phi() - parton->Phi() );
+ fhEtaPhiSpread2->Fill(jethighest->Eta()-parton->Eta(),jethighest->Phi() - parton->Phi());
+ fhJetEtSecond2->Fill(jetsecond->Energy());
+ fhJetEtRatio2->Fill(jetsecond->Energy()/jethighest->Energy());
+ fhEtaPhiDist2->Fill( TMath::Sqrt((jethighest->Eta() - jetsecond->Eta())*(jethighest->Eta() - jetsecond->Eta())
+ (jethighest->Phi() - jetsecond->Phi())*(jethighest->Phi() - jetsecond->Phi()) ));
/*
Float_t *pt,*phi,*eta;
Double_t srp = sin(parton->Phi());
Double_t alpha = acos(crp*ctp*srt*stt+srp*stp*srt*stt+crt*ctt);
Double_t correctp = pt[iT]/stt;
- hPartonPL2->Fill( correctp*cos(alpha));
+ fhPartonPL2->Fill( correctp*cos(alpha));
if ( (parton->Eta()-eta[iT])*(parton->Eta()-eta[iT]) +
(parton->Phi()-phi[iT])*(parton->Phi()-phi[iT]) < 0.2*0.2 )
- hPartonJT2->Fill( correctp*sin(alpha));
+ fhPartonJT2->Fill( correctp*sin(alpha));
if (fNominalEnergy == 0.0) {
- hPartonFragmFcn2->Fill( correctp*sin(tt)/parton->Energy() );
+ fhPartonFragmFcn2->Fill( correctp*sin(tt)/parton->Energy() );
}else
{
- hPartonFragmFcn2->Fill(correctp*sin(tt)/fNominalEnergy);
+ fhPartonFragmFcn2->Fill(correctp*sin(tt)/fNominalEnergy);
}
}// loop over tracks
Double_t srp = sin(jethighest->Phi());
Double_t alpha = acos(crp*ctp*srt*stt+srp*stp*srt*stt+crt*ctt);
Double_t correctp = pt[iT]/stt;
- hJetPL2->Fill( correctp*cos(alpha));
+ fhJetPL2->Fill( correctp*cos(alpha));
if ( (jethighest->Eta()-eta[iT])*(jethighest->Eta()-eta[iT]) +
(jethighest->Phi()-phi[iT])*(jethighest->Phi()-phi[iT]) < 0.2*0.2 )
- hJetJT2->Fill( correctp*sin(alpha));
+ fhJetJT2->Fill( correctp*sin(alpha));
if (fNominalEnergy==0.0){
- hFragmFcn2->Fill( correctp*sin(tt)/parton->Energy() );
+ fhFragmFcn2->Fill( correctp*sin(tt)/parton->Energy() );
} else
{
- hFragmFcn2->Fill( correctp*sin(tt)/fNominalEnergy );
+ fhFragmFcn2->Fill( correctp*sin(tt)/fNominalEnergy );
}
}// loop over tracks
}
AliEMCALParton* parton;
AliEMCALJet* jet;
jet = fOutput->GetJet(0);
- hJetEt->Fill(jet->Energy());
- hJetEta->Fill(jet->Eta() );
- hJetPhi->Fill(jet->Phi() );
+ fhJetEt->Fill(jet->Energy());
+ fhJetEta->Fill(jet->Eta() );
+ fhJetPhi->Fill(jet->Phi() );
if (nPartons ==0) return;
parton = fOutput->GetParton(0);
- hPartonEta->Fill( parton->Eta() );
- hPartonPhi->Fill( parton->Phi() );
+ fhPartonEta->Fill( parton->Eta() );
+ fhPartonPhi->Fill( parton->Phi() );
//hJetEtDiff->Fill( jet->Energy() - parton->Energy() );
- hEtaDiff->Fill( jet->Eta() - parton->Eta() );
- hPhiDiff->Fill( jet->Phi() - parton->Phi() );
- hEtaPhiSpread->Fill(jet->Eta()-parton->Eta(),jet->Phi() - parton->Phi());
+ fhEtaDiff->Fill( jet->Eta() - parton->Eta() );
+ fhPhiDiff->Fill( jet->Phi() - parton->Phi() );
+ fhEtaPhiSpread->Fill(jet->Eta()-parton->Eta(),jet->Phi() - parton->Phi());
/*
Float_t *pt,*phi,*eta;
Int_t *pdg;
Double_t srp = sin(parton->Phi());
Double_t alpha = acos(crp*ctp*srt*stt+srp*stp*srt*stt+crt*ctt);
Double_t correctp = pt[iT]/stt;
- hPartonPL->Fill( correctp*cos(alpha));
+ fhPartonPL->Fill( correctp*cos(alpha));
if ( (parton->Eta()-eta[iT])*(parton->Eta()-eta[iT]) +
(parton->Phi()-phi[iT])*(parton->Phi()-phi[iT]) < 0.2*0.2 )
- hPartonJT->Fill( correctp*sin(alpha));
+ fhPartonJT->Fill( correctp*sin(alpha));
if (fNominalEnergy == 0.0) {
- hPartonFragmFcn->Fill( correctp*sin(tt)/parton->Energy() );
+ fhPartonFragmFcn->Fill( correctp*sin(tt)/parton->Energy() );
}else
{
- hPartonFragmFcn->Fill(correctp*sin(tt)/fNominalEnergy);
+ fhPartonFragmFcn->Fill(correctp*sin(tt)/fNominalEnergy);
}
}// loop over tracks
Double_t srp = sin(jet->Phi());
Double_t alpha = acos(crp*ctp*srt*stt+srp*stp*srt*stt+crt*ctt);
Double_t correctp = pt[iT]/stt;
- hJetPL->Fill( correctp*cos(alpha));
+ fhJetPL->Fill( correctp*cos(alpha));
if ( (jet->Eta()-eta[iT])*(jet->Eta()-eta[iT]) +
(jet->Phi()-phi[iT])*(jet->Phi()-phi[iT]) < 0.2*0.2 )
- hJetJT->Fill( correctp*sin(alpha));
+ fhJetJT->Fill( correctp*sin(alpha));
if (fNominalEnergy==0.0){
- hFragmFcn->Fill( correctp*sin(tt)/parton->Energy() );
+ fhFragmFcn->Fill( correctp*sin(tt)/parton->Energy() );
} else
{
- hFragmFcn->Fill( correctp*sin(tt)/fNominalEnergy );
+ fhFragmFcn->Fill( correctp*sin(tt)/fNominalEnergy );
}
}// loop over tracks
}
void FillFromOutput(AliEMCALJetFinderOutput* output);
//========================== CASE 1 ========================
// Only consider events with only 1 jet
- TH1F* GetFragmFcn(){return hFragmFcn;}
- TH1F* GetPartonFragmFcn(){return hPartonFragmFcn;}
- TH1F* GetJetJT(){return hJetJT;}
- TH1F* GetPartonJT(){return hPartonJT;}
- TH1F* GetJetPL(){return hJetPL;}
- TH1F* GetPartonPL(){return hPartonPL;}
- TH1F* GetJetEt(){return hJetEt;}
- TH1F* GetJetEta(){return hJetEta;}
- TH1F* GetPartonEta(){return hPartonEta;}
- TH1F* GetPartonPhi(){return hPartonPhi;}
- TH1F* GetJetPhi(){return hJetPhi;}
- TH1F* GetEtaDiff(){return hEtaDiff;}
- TH1F* GetPhiDiff(){return hPhiDiff;}
- TH2F* GetEtaPhiSpread(){return hEtaPhiSpread;}
- TH1F* GetNJets(){return hNJets;}
+ TH1F* GetFragmFcn(){return fhFragmFcn;}
+ TH1F* GetPartonFragmFcn(){return fhPartonFragmFcn;}
+ TH1F* GetJetJT(){return fhJetJT;}
+ TH1F* GetPartonJT(){return fhPartonJT;}
+ TH1F* GetJetPL(){return fhJetPL;}
+ TH1F* GetPartonPL(){return fhPartonPL;}
+ TH1F* GetJetEt(){return fhJetEt;}
+ TH1F* GetJetEta(){return fhJetEta;}
+ TH1F* GetPartonEta(){return fhPartonEta;}
+ TH1F* GetPartonPhi(){return fhPartonPhi;}
+ TH1F* GetJetPhi(){return fhJetPhi;}
+ TH1F* GetEtaDiff(){return fhEtaDiff;}
+ TH1F* GetPhiDiff(){return fhPhiDiff;}
+ TH2F* GetEtaPhiSpread(){return fhEtaPhiSpread;}
+ TH1F* GetNJets(){return fhNJets;}
//========================== CASE 2 ========================
// Only consider events with at least 2 jets
- TH1F* GetFragmFcn2(){return hFragmFcn2;}
- TH1F* GetPartonFragmFcn2(){return hPartonFragmFcn2;}
- TH1F* GetJetJT2(){return hJetJT2;}
- TH1F* GetPartonJT2(){return hPartonJT2;}
- TH1F* GetJetPL2(){return hJetPL2;}
- TH1F* GetPartonPL2(){return hPartonPL2;}
- TH1F* GetJetEt2(){return hJetEt2;}
- TH1F* GetJetEta2(){return hJetEta2;}
- TH1F* GetPartonEta2(){return hPartonEta2;}
- TH1F* GetPartonPhi2(){return hPartonPhi2;}
- TH1F* GetJetPhi2(){return hJetPhi2;}
- TH1F* GetEtaDiff2(){return hEtaDiff2;}
- TH1F* GetPhiDiff2(){return hPhiDiff2;}
- TH2F* GetEtaPhiSpread2(){return hEtaPhiSpread2;}
- TH1F* GetNJets2(){return hNJets2;}
- TH1F* GetJetEtSecond2(){return hJetEtSecond2;}
- TH1F* GetJetEtRatio2(){return hJetEtRatio2;}
- TH1F* GetEtaPhiDist2(){return hEtaPhiDist2;}
+ TH1F* GetFragmFcn2(){return fhFragmFcn2;}
+ TH1F* GetPartonFragmFcn2(){return fhPartonFragmFcn2;}
+ TH1F* GetJetJT2(){return fhJetJT2;}
+ TH1F* GetPartonJT2(){return fhPartonJT2;}
+ TH1F* GetJetPL2(){return fhJetPL2;}
+ TH1F* GetPartonPL2(){return fhPartonPL2;}
+ TH1F* GetJetEt2(){return fhJetEt2;}
+ TH1F* GetJetEta2(){return fhJetEta2;}
+ TH1F* GetPartonEta2(){return fhPartonEta2;}
+ TH1F* GetPartonPhi2(){return fhPartonPhi2;}
+ TH1F* GetJetPhi2(){return fhJetPhi2;}
+ TH1F* GetEtaDiff2(){return fhEtaDiff2;}
+ TH1F* GetPhiDiff2(){return fhPhiDiff2;}
+ TH2F* GetEtaPhiSpread2(){return fhEtaPhiSpread2;}
+ TH1F* GetNJets2(){return fhNJets2;}
+ TH1F* GetJetEtSecond2(){return fhJetEtSecond2;}
+ TH1F* GetJetEtRatio2(){return fhJetEtRatio2;}
+ TH1F* GetEtaPhiDist2(){return fhEtaPhiDist2;}
private:
void InitPlots();
Float_t fNominalEnergy; // Force a nominal energy - specifically for 80+20 jets
AliEMCALJetFinderOutput* fOutput; // Output object to be analysed
//===================== CASE 1 ===========================================
- TH1F *hFragmFcn; // ("hFragmFcn","Fragmentation Function",100,0,1);
- TH1F *hPartonFragmFcn;// ("hFragmFcn","Parton Fragmentation Function",100,0,1);
- TH1F *hPartonJT; // ("hPartonJT","Track Momentum Perpendicular to Parton Axis",100,0.,10.);
- TH1F *hPartonPL; // ("hPartonPL","Track Momentum Parallel to Parton Axis ",100,0.,100.);
- TH1F *hJetJT; // ("hJetJT","Track Momentum Perpendicular to Jet Axis",100,0.,10.);
- TH1F *hJetPL; // ("hJetPL","Track Momentum Parallel to Jet Axis ",100,0.,100.);
- TH1F *hJetEt; // ("hJetEt","E_{T}^{reco}",250,0.,250.);
- TH1F *hJetEta; // ("hJetEta","#eta_{jet}^{reco}",180,-0.9,0.9);
- TH1F *hJetPhi; // ("hJetPhi","#phi_{jet}^{reco}",62,0.,3.1);
- TH1F *hPartonEta; // ("hPartonEta","#eta_{Parton}",180,-0.9,0.9);
- TH1F *hPartonPhi; // ("hPartonPhi","#phi_{Parton}",62,0.,3.1);
- TH1F *hEtaDiff; // ("hEtaDiff","#eta_{jet}^{reco}-#eta_{jet}^{input}",100,-0.5,0.5);
- TH1F *hPhiDiff; // ("hPhiDiff","#phi_{jet}^{reco}-#phi_{jet}^{input}",100,-0.5,0.5);
- TH2F *hEtaPhiSpread; // ("hEtaPhiSpread","#eta - #phi Distribution
+ TH1F *fhFragmFcn; // ("hFragmFcn","Fragmentation Function",100,0,1);
+ TH1F *fhPartonFragmFcn;// ("hFragmFcn","Parton Fragmentation Function",100,0,1);
+ TH1F *fhPartonJT; // ("hPartonJT","Track Momentum Perpendicular to Parton Axis",100,0.,10.);
+ TH1F *fhPartonPL; // ("hPartonPL","Track Momentum Parallel to Parton Axis ",100,0.,100.);
+ TH1F *fhJetJT; // ("hJetJT","Track Momentum Perpendicular to Jet Axis",100,0.,10.);
+ TH1F *fhJetPL; // ("hJetPL","Track Momentum Parallel to Jet Axis ",100,0.,100.);
+ TH1F *fhJetEt; // ("hJetEt","E_{T}^{reco}",250,0.,250.);
+ TH1F *fhJetEta; // ("hJetEta","#eta_{jet}^{reco}",180,-0.9,0.9);
+ TH1F *fhJetPhi; // ("hJetPhi","#phi_{jet}^{reco}",62,0.,3.1);
+ TH1F *fhPartonEta; // ("hPartonEta","#eta_{Parton}",180,-0.9,0.9);
+ TH1F *fhPartonPhi; // ("hPartonPhi","#phi_{Parton}",62,0.,3.1);
+ TH1F *fhEtaDiff; // ("hEtaDiff","#eta_{jet}^{reco}-#eta_{jet}^{input}",100,-0.5,0.5);
+ TH1F *fhPhiDiff; // ("hPhiDiff","#phi_{jet}^{reco}-#phi_{jet}^{input}",100,-0.5,0.5);
+ TH2F *fhEtaPhiSpread; // ("hEtaPhiSpread","#eta - #phi Distribution
//of Reconstructed Jets",192,-0.7,0.7,288,pi/3,pi);
- TH1F *hNJets; // ("hNJets","N Reconstructed jets",11,-0.5,10.5);
+ TH1F *fhNJets; // ("hNJets","N Reconstructed jets",11,-0.5,10.5);
//============================== CASE 2 ============================================
- TH1F *hFragmFcn2; // ("hFragmFcn2","Fragmentation Function",100,0,1);
- TH1F *hPartonFragmFcn2;// ("hFragmFcn2","Parton Fragmentation Function",100,0,1);
- TH1F *hPartonJT2; // ("hPartonJT2","Track Momentum Perpendicular to Parton Axis",100,0.,10.);
- TH1F *hPartonPL2; // ("hPartonPL2","Track Momentum Parallel to Parton Axis ",100,0.,100.);
- TH1F *hJetJT2; // ("hJetJT2","Track Momentum Perpendicular to Jet Axis",100,0.,10.);
- TH1F *hJetPL2; // ("hJetPL2","Track Momentum Parallel to Jet Axis ",100,0.,100.);
- TH1F *hJetEt2; // ("hJetEt2","E_{T}^{reco}",250,0.,250.);
- TH1F *hJetEta2; // ("hJetEta2","#eta_{jet}^{reco}",180,-0.9,0.9);
- TH1F *hJetPhi2; // ("hJetPhi2","#phi_{jet}^{reco}",62,0.,3.1);
- TH1F *hPartonEta2; // ("hPartonEta2","#eta_{Parton}",180,-0.9,0.9);
- TH1F *hPartonPhi2; // ("hPartonPhi2","#phi_{Parton}",62,0.,3.1);
- TH1F *hEtaDiff2; // ("hEtaDiff2","#eta_{jet}^{reco}-#eta_{jet}^{input}",100,-0.5,0.5);
- TH1F *hPhiDiff2; // ("hPhiDiff2","#phi_{jet}^{reco}-#phi_{jet}^{input}",100,-0.5,0.5);
- TH2F *hEtaPhiSpread2; // ("hEtaPhiSpread2","#eta - #phi Distribution
+ TH1F *fhFragmFcn2; // ("hFragmFcn2","Fragmentation Function",100,0,1);
+ TH1F *fhPartonFragmFcn2;// ("hFragmFcn2","Parton Fragmentation Function",100,0,1);
+ TH1F *fhPartonJT2; // ("hPartonJT2","Track Momentum Perpendicular to Parton Axis",100,0.,10.);
+ TH1F *fhPartonPL2; // ("hPartonPL2","Track Momentum Parallel to Parton Axis ",100,0.,100.);
+ TH1F *fhJetJT2; // ("hJetJT2","Track Momentum Perpendicular to Jet Axis",100,0.,10.);
+ TH1F *fhJetPL2; // ("hJetPL2","Track Momentum Parallel to Jet Axis ",100,0.,100.);
+ TH1F *fhJetEt2; // ("hJetEt2","E_{T}^{reco}",250,0.,250.);
+ TH1F *fhJetEta2; // ("hJetEta2","#eta_{jet}^{reco}",180,-0.9,0.9);
+ TH1F *fhJetPhi2; // ("hJetPhi2","#phi_{jet}^{reco}",62,0.,3.1);
+ TH1F *fhPartonEta2; // ("hPartonEta2","#eta_{Parton}",180,-0.9,0.9);
+ TH1F *fhPartonPhi2; // ("hPartonPhi2","#phi_{Parton}",62,0.,3.1);
+ TH1F *fhEtaDiff2; // ("hEtaDiff2","#eta_{jet}^{reco}-#eta_{jet}^{input}",100,-0.5,0.5);
+ TH1F *fhPhiDiff2; // ("hPhiDiff2","#phi_{jet}^{reco}-#phi_{jet}^{input}",100,-0.5,0.5);
+ TH2F *fhEtaPhiSpread2; // ("hEtaPhiSpread2","#eta - #phi Distribution
//of Reconstructed Jets",192,-0.7,0.7,288,pi/3,pi);
- TH1F *hNJets2; // ("hNJets2","N Reconstructed jets",11,-0.5,10.5);
- TH1F *hJetEtSecond2; //("hJetEtSecond2","E_{T}^{reco}",250,0.,250.);
- TH1F *hJetEtRatio2; //("hJetEtRatio2","Ratio of Second Highest to Highest",100,0,1);
- TH1F *hEtaPhiDist2; //("hEtaPhiDist2","Angular Distance Between First and Second",100,0,3);
+ TH1F *fhNJets2; // ("hNJets2","N Reconstructed jets",11,-0.5,10.5);
+ TH1F *fhJetEtSecond2; //("hJetEtSecond2","E_{T}^{reco}",250,0.,250.);
+ TH1F *fhJetEtRatio2; //("hJetEtRatio2","Ratio of Second Highest to Highest",100,0,1);
+ TH1F *fhEtaPhiDist2; //("hEtaPhiDist2","Angular Distance Between First and Second",100,0,3);