Fixed names of histograms
authormhorner <mhorner@f7af4fe6-9843-0410-8265-dc069ae4e863>
Wed, 21 Jan 2004 16:43:22 +0000 (16:43 +0000)
committermhorner <mhorner@f7af4fe6-9843-0410-8265-dc069ae4e863>
Wed, 21 Jan 2004 16:43:22 +0000 (16:43 +0000)
EMCAL/AliEMCALJetFinderPlots.cxx
EMCAL/AliEMCALJetFinderPlots.h

index bf180ee..072ae4f 100755 (executable)
@@ -38,152 +38,152 @@ AliEMCALJetFinderPlots::AliEMCALJetFinderPlots()
   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;        
   
@@ -191,41 +191,41 @@ hEtaPhiDist2->Sumw2();
 
 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);
 
 
 
@@ -236,12 +236,12 @@ void AliEMCALJetFinderPlots::FillFromOutput(AliEMCALJetFinderOutput* output)
 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; 
@@ -269,22 +269,22 @@ if (fOutput->GetNJets()>1)
   }
 
   // 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;
@@ -318,15 +318,15 @@ if (fOutput->GetNJets()>1)
       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
 
@@ -350,15 +350,15 @@ if (fOutput->GetNJets()>1)
          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
   }
@@ -372,19 +372,19 @@ if (fOutput->GetNJets()>1)
   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;
@@ -417,15 +417,15 @@ if (fOutput->GetNJets()>1)
       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
 
@@ -449,15 +449,15 @@ if (fOutput->GetNJets()>1)
          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
   }
index 1dc5470..5f0d126 100755 (executable)
@@ -33,42 +33,42 @@ class AliEMCALJetFinderPlots : public TObject
        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();
@@ -77,44 +77,44 @@ class AliEMCALJetFinderPlots : public TObject
        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);