changed binning, vertex range, leading, subleading
[u/mrichter/AliRoot.git] / PWGJE / AliAnalysisTaskJetCore.h
index 98674f8fdd4ae210d9ef64bea6ac45f15a63595f..efcdcf351f497989cb194947b6aef80ce856660a 100644 (file)
@@ -1,7 +1,19 @@
 #ifndef ALIANALYSISTASKJETCORE_H
 #define ALIANALYSISTASKJETCORE_H
 
+/* Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved. *
+ * See cxx source for full Copyright notice                               */
+
+// **************************************
+// This task computes several jet observables like 
+// the fraction of energy in inner and outer coronnas,
+// the distance from track to jet axis and a 
+// correlation strength distribution of particles inside jets.    
+// Author: lcunquei@cern.ch
+// *******************************************
+
 class TH1F;
+class TH1I;
 class TH2F;
 class TH3F;
 class THnSparse;
@@ -120,102 +132,244 @@ private:
    TH2F  *fh2JetSelection;                    //! jet selection statistic, with 
      
 
-   TH2F      *fh2JetCoreMethod1C10;            //energy fraction in the jet core for
-   TH2F      *fh2JetCoreMethod2C10;            //different centralities and 3 methods 
-   TH2F      *fh2JetCoreMethod3C10;
-   TH2F      *fh2JetCoreMethod1C20;
-   TH2F      *fh2JetCoreMethod2C20;
-   TH2F      *fh2JetCoreMethod3C20;
-   TH2F      *fh2JetCoreMethod1C30;
-   TH2F      *fh2JetCoreMethod2C30;
-   TH2F      *fh2JetCoreMethod3C30;
-   TH2F      *fh2JetCoreMethod1C60;
-   TH2F      *fh2JetCoreMethod2C60;
-   TH2F      *fh2JetCoreMethod3C60;
-   TH2F      *fh2SumPtInC10;                 //energy fraction in external coronas
-   TH2F      *fh2SumPtInC20;
-   TH2F      *fh2SumPtInC30;
-   TH2F      *fh2SumPtInC60;
-   TH2F      *fh2SumPtOutC10;
-   TH2F      *fh2SumPtOutC10b;   
-   TH2F      *fh2SumPtOutC20;
-   TH2F      *fh2SumPtOutC30;
-   TH2F      *fh2SumPtOutC60;
-   TH2F      *fh2SumPtInC10bkg;        //energy fraction in external coronas expected
-   TH2F      *fh2SumPtInC20bkg;        //only from background
-   TH2F      *fh2SumPtInC30bkg;
-   TH2F      *fh2SumPtInC60bkg;
-   TH2F      *fh2SumPtOutC10bkg;
-   TH2F      *fh2SumPtOutC20bkg;
-   TH2F      *fh2SumPtOutC30bkg;
-   TH2F      *fh2SumPtOutC60bkg;
+   TH2F      *fh2JetCoreMethod1C10;          //Energy fraction in the core C10 method 1
+   TH2F      *fh2JetCoreMethod2C10;          //Energy fraction in the core C10 method 2
+   TH2F      *fh2JetCoreMethod3C10;          //Energy fraction in the core C10 method 3
+   TH2F      *fh2JetCoreMethod1C20;          //Energy fraction in the core C20 method 1 
+   TH2F      *fh2JetCoreMethod2C20;          //Energy fraction in the core C20 method 2
+   TH2F      *fh2JetCoreMethod3C20;          //Energy fraction in the core C20 method 3
+   TH2F      *fh2JetCoreMethod1C30;          //Energy fraction in the core C30 method 1
+   TH2F      *fh2JetCoreMethod2C30;          //Energy fraction in the core C30 method 2
+   TH2F      *fh2JetCoreMethod3C30;          //Energy fraction in the core C30 method 3
+   TH2F      *fh2JetCoreMethod1C60;          //Energy fraction in the core C60 method 1
+   TH2F      *fh2JetCoreMethod2C60;          //Energy fraction in the core C60 method 2
+   TH2F      *fh2JetCoreMethod3C60;          //Energy fraction in the core C60 method 3
+   TH2F      *fh2JetCoreMethod3C10lead;          //Energy fraction in the core C30 method 3
+   TH2F      *fh2JetCoreMethod3C20lead;          //Energy fraction in the core C60 method 1
+   TH2F      *fh2JetCoreMethod3C30lead;          //Energy fraction in the core C60 method 2
+   TH2F      *fh2JetCoreMethod3C60lead;          //Energy fraction in the core C60 method 3
+   TH2F      *fh2JetCoreMethod3C10sublead;          //Energy fraction in the core C30 method 3
+   TH2F      *fh2JetCoreMethod3C20sublead;          //Energy fraction in the core C60 method 1
+   TH2F      *fh2JetCoreMethod3C30sublead;          //Energy fraction in the core C60 method 2
+   TH2F      *fh2JetCoreMethod3C60sublead;          //Energy fraction in the core C60 method 3
 
-  
+   TH2F      *fh2SumPtInC10;                  //energy fraction in inner corona C10
+   TH2F      *fh2SumPtInC20;                  //energy fraction in inner corona C20 
+   TH2F      *fh2SumPtInC30;                  //energy fraction in inner corona C30
+   TH2F      *fh2SumPtInC60;                  //energy fraction in inner corona C60
+   TH2F      *fh2SumPtInC10lead;              //energy fraction in inner corona C10 leading
+   TH2F      *fh2SumPtInC20lead;              //energy fraction in inner corona C20 leading
+   TH2F      *fh2SumPtInC30lead;              //energy fraction in inner corona C30 leading
+   TH2F      *fh2SumPtInC60lead;              //energy fraction in inner corona C60 leading
+   TH2F      *fh2SumPtInC10sublead;           //energy fraction in inner corona C10 subleading
+   TH2F      *fh2SumPtInC20sublead;           //energy fraction in inner corona C20 subleading
+   TH2F      *fh2SumPtInC30sublead;           //energy fraction in inner corona C30 subleading
+   TH2F      *fh2SumPtInC60sublead;           //energy fraction in inner corona C60 subleading
+   TH2F      *fh2SumPtOutC10;                 //energy fraction in outer corona C10
+   TH2F      *fh2SumPtOutC20;                 //energy fraction in outer corona C20 
+   TH2F      *fh2SumPtOutC30;                 //energy fraction in outer corona C30
+   TH2F      *fh2SumPtOutC60;                 //energy fraction in outer corona C60
+   TH2F      *fh2SumPtOutC10lead;              //energy fraction in outer corona C10 leading
+   TH2F      *fh2SumPtOutC20lead;              //energy fraction in outer corona C20 leading
+   TH2F      *fh2SumPtOutC30lead;              //energy fraction in outer corona C30 leading
+   TH2F      *fh2SumPtOutC60lead;              //energy fraction in outer corona C60 leading
+   TH2F      *fh2SumPtOutC10sublead;           //energy fraction in outer corona C10 subleading
+   TH2F      *fh2SumPtOutC20sublead;           //energy fraction in outer corona C20 subleading
+   TH2F      *fh2SumPtOutC30sublead;           //energy fraction in outer corona C30 subleading
+   TH2F      *fh2SumPtOutC60sublead;           //energy fraction in outer corona C60 subleading
+   TH2F      *fh2SumPtInC10bkg;        //expected from background inner C10
+   TH2F      *fh2SumPtInC20bkg;        //expected from background inner C20
+   TH2F      *fh2SumPtInC30bkg;        //expected from background inner C30
+   TH2F      *fh2SumPtInC60bkg;        //expected from background inner C60
+   TH2F      *fh2SumPtInC10bkglead;        //expected from background inner C10 lead
+   TH2F      *fh2SumPtInC20bkglead;        //expected from background inner C20 lead
+   TH2F      *fh2SumPtInC30bkglead;        //expected from background inner C30 lead
+   TH2F      *fh2SumPtInC60bkglead;        //expected from background inner C60 lead
+   TH2F      *fh2SumPtInC10bkgsublead;        //expected from background inner C10 sublead
+   TH2F      *fh2SumPtInC20bkgsublead;        //expected from background inner C20 sublead  
+   TH2F      *fh2SumPtInC30bkgsublead;        //expected from background inner C30  sublead
+   TH2F      *fh2SumPtInC60bkgsublead;        //expected from background inner C60 sublead
+
+   TH2F      *fh2SumPtOutC10bkg;       //expected from background outer C10
+   TH2F      *fh2SumPtOutC20bkg;       //expected from background outer C10
+   TH2F      *fh2SumPtOutC30bkg;       //expected from background outer C10
+   TH2F      *fh2SumPtOutC60bkg;       //expected from background outer C10
+   TH2F      *fh2SumPtOutC10bkglead;       //expected from background outer C10 lead
+   TH2F      *fh2SumPtOutC20bkglead;       //expected from background outer C10 lead
+   TH2F      *fh2SumPtOutC30bkglead;       //expected from background outer C10 lead
+   TH2F      *fh2SumPtOutC60bkglead;       //expected from background outer C10 lead
+   TH2F      *fh2SumPtOutC10bkgsublead;       //expected from background outer C10 sublead
+   TH2F      *fh2SumPtOutC20bkgsublead;       //expected from background outer C10 sublead
+   TH2F      *fh2SumPtOutC30bkgsublead;       //expected from background outer C10 sublead
+   TH2F      *fh2SumPtOutC60bkgsublead;       //expected from background outer C10 sublead  
    
-     TH2F*      fh2DeltaRC10pt1;            //Jet track R distance in 3 high jet pT bins
-     TH2F*      fh2DeltaRC20pt1;            //for different centralities. 
-     TH2F*      fh2DeltaRC30pt1;
-     TH2F*      fh2DeltaRC60pt1;
-     TH2F*      fh2DeltaRC10pt2;
-     TH2F*      fh2DeltaRC20pt2;
-     TH2F*      fh2DeltaRC30pt2;
-     TH2F*      fh2DeltaRC60pt2;
-     TH2F*      fh2DeltaRC10pt3;
-     TH2F*      fh2DeltaRC20pt3;
-     TH2F*      fh2DeltaRC30pt3;
-     TH2F*      fh2DeltaRC60pt3;
-     TH2F*      fh2DeltaRC10pt4;
-     TH2F*      fh2DeltaRC20pt4;
-     TH2F*      fh2DeltaRC30pt4;
-     TH2F*      fh2DeltaRC60pt4;
-     TH2F*      fh2DeltaEtaC10pt1;             //The same but eta distance
-     TH2F*      fh2DeltaEtaC20pt1;
-     TH2F*      fh2DeltaEtaC30pt1;
-     TH2F*      fh2DeltaEtaC60pt1;
-     TH2F*      fh2DeltaEtaC10pt2;
-     TH2F*      fh2DeltaEtaC20pt2;
-     TH2F*      fh2DeltaEtaC30pt2;
-     TH2F*      fh2DeltaEtaC60pt2;
-     TH2F*      fh2DeltaEtaC10pt3;
-     TH2F*      fh2DeltaEtaC20pt3;
-     TH2F*      fh2DeltaEtaC30pt3;
-     TH2F*      fh2DeltaEtaC60pt3;
-     TH2F*      fh2DeltaEtaC10pt4;
-     TH2F*      fh2DeltaEtaC20pt4;
-     TH2F*      fh2DeltaEtaC30pt4;
-     TH2F*      fh2DeltaEtaC60pt4;
-     TH2F*      fh2DeltaPhiC10pt1;             //The same but phi distance
-     TH2F*      fh2DeltaPhiC20pt1;
-     TH2F*      fh2DeltaPhiC30pt1;
-     TH2F*      fh2DeltaPhiC60pt1;
-     TH2F*      fh2DeltaPhiC10pt2;
-     TH2F*      fh2DeltaPhiC20pt2;
-     TH2F*      fh2DeltaPhiC30pt2;
-     TH2F*      fh2DeltaPhiC60pt2;
-     TH2F*      fh2DeltaPhiC10pt3;
-     TH2F*      fh2DeltaPhiC20pt3;
-     TH2F*      fh2DeltaPhiC30pt3;
-     TH2F*      fh2DeltaPhiC60pt3;
-     TH2F*      fh2DeltaPhiC10pt4;
-     TH2F*      fh2DeltaPhiC20pt4;
-     TH2F*      fh2DeltaPhiC30pt4;
-     TH2F*      fh2DeltaPhiC60pt4;
-     TH2F*      fh2AngStructpt1C10;           //Average two particle correlation function
-     TH2F*      fh2AngStructpt2C10; 
-     TH2F*      fh2AngStructpt3C10;
-     TH2F*      fh2AngStructpt4C10;     
-     TH2F*      fh2AngStructpt1C20;
-     TH2F*      fh2AngStructpt2C20;
-     TH2F*      fh2AngStructpt3C20;
-     TH2F*      fh2AngStructpt4C20;     
-     TH2F*      fh2AngStructpt1C30;
-     TH2F*      fh2AngStructpt2C30;
-     TH2F*      fh2AngStructpt3C30;
-     TH2F*      fh2AngStructpt4C30;     
-     TH2F*      fh2AngStructpt1C60;
-     TH2F*      fh2AngStructpt2C60;
-     TH2F*      fh2AngStructpt3C60;
-     TH2F*      fh2AngStructpt4C60;     
+     TH2F*      fh2DeltaRC10pt1;            //Jet track R distance:C10 pt1
+     TH2F*      fh2DeltaRC20pt1;            //C20 pt1 
+     TH2F*      fh2DeltaRC30pt1;            //C30 pt1
+     TH2F*      fh2DeltaRC60pt1;            //C60 pt1
+     TH2F*      fh2DeltaRC10pt2;            //C10 pt2   
+     TH2F*      fh2DeltaRC20pt2;            //C20 pt2
+     TH2F*      fh2DeltaRC30pt2;            //C30 pt2
+     TH2F*      fh2DeltaRC60pt2;            //C60 pt2 
+     TH2F*      fh2DeltaRC10pt3;            //C10 pt3
+     TH2F*      fh2DeltaRC20pt3;            //C20 pt3
+     TH2F*      fh2DeltaRC30pt3;            //C30 pt3
+     TH2F*      fh2DeltaRC60pt3;            //C60 pt3
+     TH2F*      fh2DeltaRC10pt4;            //C10 pt4
+     TH2F*      fh2DeltaRC20pt4;            //C20 pt4
+     TH2F*      fh2DeltaRC30pt4;            //C30 pt4
+     TH2F*      fh2DeltaRC60pt4;            //C60 pt4 
+     TH2F*      fh2DeltaEtaC10pt1;          //The same but eta distance:C10 pt1
+     TH2F*      fh2DeltaEtaC20pt1;          //C20 pt1
+     TH2F*      fh2DeltaEtaC30pt1;          //C30 pt1
+     TH2F*      fh2DeltaEtaC60pt1;          //C60 pt1
+     TH2F*      fh2DeltaEtaC10pt2;          //C10 pt2  
+     TH2F*      fh2DeltaEtaC20pt2;          //C20 pt2
+     TH2F*      fh2DeltaEtaC30pt2;          //C30 pt2
+     TH2F*      fh2DeltaEtaC60pt2;          //C60 pt2
+     TH2F*      fh2DeltaEtaC10pt3;          //C10 pt3
+     TH2F*      fh2DeltaEtaC20pt3;          //C20 pt3
+     TH2F*      fh2DeltaEtaC30pt3;          //C30 pt3
+     TH2F*      fh2DeltaEtaC60pt3;          //C60 pt3
+     TH2F*      fh2DeltaEtaC10pt4;          //C10 pt4
+     TH2F*      fh2DeltaEtaC20pt4;          //C20 pt4
+     TH2F*      fh2DeltaEtaC30pt4;          //C30 pt4
+     TH2F*      fh2DeltaEtaC60pt4;          //C60 pt4
+     TH2F*      fh2DeltaPhiC10pt1;          //The same but phi distance:C10 pt1
+     TH2F*      fh2DeltaPhiC20pt1;          //C20 pt1
+     TH2F*      fh2DeltaPhiC30pt1;          //C30 pt1
+     TH2F*      fh2DeltaPhiC60pt1;          //C60 pt1
+     TH2F*      fh2DeltaPhiC10pt2;          //C10 pt2
+     TH2F*      fh2DeltaPhiC20pt2;          //C20 pt2 
+     TH2F*      fh2DeltaPhiC30pt2;          //C30 pt2
+     TH2F*      fh2DeltaPhiC60pt2;          //C60 pt2 
+     TH2F*      fh2DeltaPhiC10pt3;          //C10 pt3
+     TH2F*      fh2DeltaPhiC20pt3;          //C20 pt3
+     TH2F*      fh2DeltaPhiC30pt3;          //C30 pt3
+     TH2F*      fh2DeltaPhiC60pt3;          //C60 pt3
+     TH2F*      fh2DeltaPhiC10pt4;          //C10 pt4
+     TH2F*      fh2DeltaPhiC20pt4;          //C20 pt4
+     TH2F*      fh2DeltaPhiC30pt4;          //C30 pt4
+     TH2F*      fh2DeltaPhiC60pt4;          //C60 pt4
+
+     TH2F*      fh2DeltaRC10pt1lead;            //Jet track R distance:C10 pt1
+     TH2F*      fh2DeltaRC20pt1lead;            //C20 pt1 
+     TH2F*      fh2DeltaRC30pt1lead;            //C30 pt1
+     TH2F*      fh2DeltaRC60pt1lead;            //C60 pt1
+     TH2F*      fh2DeltaRC10pt2lead;            //C10 pt2   
+     TH2F*      fh2DeltaRC20pt2lead;            //C20 pt2
+     TH2F*      fh2DeltaRC30pt2lead;            //C30 pt2
+     TH2F*      fh2DeltaRC60pt2lead;            //C60 pt2 
+     TH2F*      fh2DeltaRC10pt3lead;            //C10 pt3
+     TH2F*      fh2DeltaRC20pt3lead;            //C20 pt3
+     TH2F*      fh2DeltaRC30pt3lead;            //C30 pt3
+     TH2F*      fh2DeltaRC60pt3lead;            //C60 pt3
+     TH2F*      fh2DeltaRC10pt4lead;            //C10 pt4
+     TH2F*      fh2DeltaRC20pt4lead;            //C20 pt4
+     TH2F*      fh2DeltaRC30pt4lead;            //C30 pt4
+     TH2F*      fh2DeltaRC60pt4lead;            //C60 pt4 
+     TH2F*      fh2DeltaEtaC10pt1lead;          //The same but eta distance:C10 pt1
+     TH2F*      fh2DeltaEtaC20pt1lead;          //C20 pt1
+     TH2F*      fh2DeltaEtaC30pt1lead;          //C30 pt1
+     TH2F*      fh2DeltaEtaC60pt1lead;          //C60 pt1
+     TH2F*      fh2DeltaEtaC10pt2lead;          //C10 pt2  
+     TH2F*      fh2DeltaEtaC20pt2lead;          //C20 pt2
+     TH2F*      fh2DeltaEtaC30pt2lead;          //C30 pt2
+     TH2F*      fh2DeltaEtaC60pt2lead;          //C60 pt2
+     TH2F*      fh2DeltaEtaC10pt3lead;          //C10 pt3
+     TH2F*      fh2DeltaEtaC20pt3lead;          //C20 pt3
+     TH2F*      fh2DeltaEtaC30pt3lead;          //C30 pt3
+     TH2F*      fh2DeltaEtaC60pt3lead;          //C60 pt3
+     TH2F*      fh2DeltaEtaC10pt4lead;          //C10 pt4
+     TH2F*      fh2DeltaEtaC20pt4lead;          //C20 pt4
+     TH2F*      fh2DeltaEtaC30pt4lead;          //C30 pt4
+     TH2F*      fh2DeltaEtaC60pt4lead;          //C60 pt4
+     TH2F*      fh2DeltaPhiC10pt1lead;          //The same but phi distance:C10 pt1
+     TH2F*      fh2DeltaPhiC20pt1lead;          //C20 pt1
+     TH2F*      fh2DeltaPhiC30pt1lead;          //C30 pt1
+     TH2F*      fh2DeltaPhiC60pt1lead;          //C60 pt1
+     TH2F*      fh2DeltaPhiC10pt2lead;          //C10 pt2
+     TH2F*      fh2DeltaPhiC20pt2lead;          //C20 pt2 
+     TH2F*      fh2DeltaPhiC30pt2lead;          //C30 pt2
+     TH2F*      fh2DeltaPhiC60pt2lead;          //C60 pt2 
+     TH2F*      fh2DeltaPhiC10pt3lead;          //C10 pt3
+     TH2F*      fh2DeltaPhiC20pt3lead;          //C20 pt3
+     TH2F*      fh2DeltaPhiC30pt3lead;          //C30 pt3
+     TH2F*      fh2DeltaPhiC60pt3lead;          //C60 pt3
+     TH2F*      fh2DeltaPhiC10pt4lead;          //C10 pt4
+     TH2F*      fh2DeltaPhiC20pt4lead;          //C20 pt4
+     TH2F*      fh2DeltaPhiC30pt4lead;          //C30 pt4
+     TH2F*      fh2DeltaPhiC60pt4lead;          //C60 pt4
+
+     TH2F*      fh2DeltaRC10pt1sublead;            //Jet track R distance:C10 pt1
+     TH2F*      fh2DeltaRC20pt1sublead;            //C20 pt1 
+     TH2F*      fh2DeltaRC30pt1sublead;            //C30 pt1
+     TH2F*      fh2DeltaRC60pt1sublead;            //C60 pt1
+     TH2F*      fh2DeltaRC10pt2sublead;            //C10 pt2   
+     TH2F*      fh2DeltaRC20pt2sublead;            //C20 pt2
+     TH2F*      fh2DeltaRC30pt2sublead;            //C30 pt2
+     TH2F*      fh2DeltaRC60pt2sublead;            //C60 pt2 
+     TH2F*      fh2DeltaRC10pt3sublead;            //C10 pt3
+     TH2F*      fh2DeltaRC20pt3sublead;            //C20 pt3
+     TH2F*      fh2DeltaRC30pt3sublead;            //C30 pt3
+     TH2F*      fh2DeltaRC60pt3sublead;            //C60 pt3
+     TH2F*      fh2DeltaRC10pt4sublead;            //C10 pt4
+     TH2F*      fh2DeltaRC20pt4sublead;            //C20 pt4
+     TH2F*      fh2DeltaRC30pt4sublead;            //C30 pt4
+     TH2F*      fh2DeltaRC60pt4sublead;            //C60 pt4 
+     TH2F*      fh2DeltaEtaC10pt1sublead;          //The same but eta distance:C10 pt1
+     TH2F*      fh2DeltaEtaC20pt1sublead;          //C20 pt1
+     TH2F*      fh2DeltaEtaC30pt1sublead;          //C30 pt1
+     TH2F*      fh2DeltaEtaC60pt1sublead;          //C60 pt1
+     TH2F*      fh2DeltaEtaC10pt2sublead;          //C10 pt2  
+     TH2F*      fh2DeltaEtaC20pt2sublead;          //C20 pt2
+     TH2F*      fh2DeltaEtaC30pt2sublead;          //C30 pt2
+     TH2F*      fh2DeltaEtaC60pt2sublead;          //C60 pt2
+     TH2F*      fh2DeltaEtaC10pt3sublead;          //C10 pt3
+     TH2F*      fh2DeltaEtaC20pt3sublead;          //C20 pt3
+     TH2F*      fh2DeltaEtaC30pt3sublead;          //C30 pt3
+     TH2F*      fh2DeltaEtaC60pt3sublead;          //C60 pt3
+     TH2F*      fh2DeltaEtaC10pt4sublead;          //C10 pt4
+     TH2F*      fh2DeltaEtaC20pt4sublead;          //C20 pt4
+     TH2F*      fh2DeltaEtaC30pt4sublead;          //C30 pt4
+     TH2F*      fh2DeltaEtaC60pt4sublead;          //C60 pt4
+     TH2F*      fh2DeltaPhiC10pt1sublead;          //The same but phi distance:C10 pt1
+     TH2F*      fh2DeltaPhiC20pt1sublead;          //C20 pt1
+     TH2F*      fh2DeltaPhiC30pt1sublead;          //C30 pt1
+     TH2F*      fh2DeltaPhiC60pt1sublead;          //C60 pt1
+     TH2F*      fh2DeltaPhiC10pt2sublead;          //C10 pt2
+     TH2F*      fh2DeltaPhiC20pt2sublead;          //C20 pt2 
+     TH2F*      fh2DeltaPhiC30pt2sublead;          //C30 pt2
+     TH2F*      fh2DeltaPhiC60pt2sublead;          //C60 pt2 
+     TH2F*      fh2DeltaPhiC10pt3sublead;          //C10 pt3
+     TH2F*      fh2DeltaPhiC20pt3sublead;          //C20 pt3
+     TH2F*      fh2DeltaPhiC30pt3sublead;          //C30 pt3
+     TH2F*      fh2DeltaPhiC60pt3sublead;          //C60 pt3
+     TH2F*      fh2DeltaPhiC10pt4sublead;          //C10 pt4
+     TH2F*      fh2DeltaPhiC20pt4sublead;          //C20 pt4
+     TH2F*      fh2DeltaPhiC30pt4sublead;          //C30 pt4
+     TH2F*      fh2DeltaPhiC60pt4sublead;          //C60 pt4
+
+
+
+
+
+     TH2F*      fh2AngStructpt1C10;         //Average two particle correlation function:C10 pt1
+     TH2F*      fh2AngStructpt2C10;         //C10 pt2
+     TH2F*      fh2AngStructpt3C10;         //C10 pt3
+     TH2F*      fh2AngStructpt4C10;         //C10 pt4
+     TH2F*      fh2AngStructpt1C20;         //C20 pt1
+     TH2F*      fh2AngStructpt2C20;         //C20 pt2
+     TH2F*      fh2AngStructpt3C20;         //C20 pt3 
+     TH2F*      fh2AngStructpt4C20;         //C20 pt4
+     TH2F*      fh2AngStructpt1C30;         //C30 pt1
+     TH2F*      fh2AngStructpt2C30;         //C30 pt2
+     TH2F*      fh2AngStructpt3C30;         //C30 pt3
+     TH2F*      fh2AngStructpt4C30;         //C30 pt4
+     TH2F*      fh2AngStructpt1C60;         //C60 pt1
+     TH2F*      fh2AngStructpt2C60;         //C60 pt2
+     TH2F*      fh2AngStructpt3C60;         //C60 pt3
+     TH2F*      fh2AngStructpt4C60;         //C60 pt4