fgMCEvGen changed to fMCEvGen.

[u/mrichter/AliRoot.git] / EMCAL / AliEMCALJetFinderPlots.h

diff --git a/EMCAL/AliEMCALJetFinderPlots.h b/EMCAL/AliEMCALJetFinderPlots.h
index 343d608b23ab94ccf8651893550fd57d61565a0c..1dc54705af5a941b4be6d9428f1e1f8ae8788c96 100755 (executable)
--- a/EMCAL/AliEMCALJetFinderPlots.h
+++ b/EMCAL/AliEMCALJetFinderPlots.h
@@ -4,6 +4,7 @@
 /* Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved. *
  *  *  * See cxx source for full Copyright notice     */
 
+
 /* $Id$ */
 
 //_________________________________________________________________________
@@ -30,6 +31,8 @@ class AliEMCALJetFinderPlots : public TObject
        void SetNominalEnergy(Float_t energy){fNominalEnergy = energy;}
        void SetDebug(Int_t debug){fDebug = debug;}
        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;}        
@@ -45,13 +48,35 @@ class AliEMCALJetFinderPlots : public TObject
        TH1F* GetPhiDiff(){return hPhiDiff;}
        TH2F* GetEtaPhiSpread(){return hEtaPhiSpread;}
        TH1F* GetNJets(){return hNJets;}
-       
+
+       //========================== 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;}
+
        private:
        void InitPlots();
        Int_t                           fDebug;         // Debug value
        Float_t                         fConeRadius;    // Cone radius to be used in filling
        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.);
@@ -68,10 +93,36 @@ class AliEMCALJetFinderPlots : public TObject
        TH2F                            *hEtaPhiSpread; // ("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);
+  
+       //============================== 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 
+                                                       //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);
+
+
+
+       
        Bool_t                          fInitialised;
                
        
-       ClassDef(AliEMCALJetFinderPlots,2)
+       ClassDef(AliEMCALJetFinderPlots,3)
        
 };
 #endif