Updates Leticia
authormorsch <morsch@f7af4fe6-9843-0410-8265-dc069ae4e863>
Mon, 18 Jun 2012 15:45:20 +0000 (15:45 +0000)
committermorsch <morsch@f7af4fe6-9843-0410-8265-dc069ae4e863>
Mon, 18 Jun 2012 15:45:20 +0000 (15:45 +0000)
PWGJE/AliAnalysisTaskJetCore.cxx
PWGJE/AliAnalysisTaskJetCore.h

index a163fe0..6f7b635 100644 (file)
@@ -80,6 +80,9 @@ fAngStructCloseTracks(0),
 fCheckMethods(0),
 fDoEventMixing(0), 
 fFlagPhiBkg(0),
+fFlagRandom(0),
+fRPAngle(0),
+fNRPBins(3),
 fJetEtaMin(-.5),
 fJetEtaMax(.5),
 fNevents(0),
@@ -107,8 +110,9 @@ fh2JetCoreMethod1C30(0x0),
 fh2JetCoreMethod2C30(0x0),
 fh2JetCoreMethod1C60(0x0),
 fh2JetCoreMethod2C60(0x0),
-fh3JetTrackC10(0x0),
+fh3JetTrackC3060(0x0),
 fh3JetTrackC20(0x0),
+fh3JetTrackC4080(0x0), 
 fh2AngStructpt1C10(0x0),
 fh2AngStructpt2C10(0x0),
 fh2AngStructpt3C10(0x0),
@@ -125,75 +129,14 @@ fh2AngStructpt1C60(0x0),
 fh2AngStructpt2C60(0x0),
 fh2AngStructpt3C60(0x0),
 fh2AngStructpt4C60(0x0),
-fh2JetsumHT3R2a(0x0),
-fh2JetsumHT3R2ap(0x0),
-fh2JetsumHT3R4a(0x0),
-fh2JetsumHT3R4ap(0x0),
-fh2JetsumHT3R6a(0x0),
-fh2JetsumHT3R6ap(0x0),
-fh2JetsumHT3R8a(0x0),
-fh2JetsumHT3R8ap(0x0),
-fh2JetsumHT3R10a(0x0),
-fh2JetsumHT3R10ap(0x0),
-fh2JetsumHT3R2aa(0x0),
-fh2JetsumHT3R2aap(0x0),
-fh2JetsumHT3R4aa(0x0),
-fh2JetsumHT3R4aap(0x0),
-fh2JetsumHT3R6aa(0x0),
-fh2JetsumHT3R6aap(0x0),
-fh2JetsumHT3R8aa(0x0),
-fh2JetsumHT3R8aap(0x0),
-fh2JetsumHT3R10aa(0x0),
-fh2JetsumHT3R10aap(0x0),
-fh2JetsumHT3R2aaa(0x0),
-fh2JetsumHT3R2aaap(0x0),
-fh2JetsumHT3R4aaa(0x0),
-fh2JetsumHT3R4aaap(0x0),
-fh2JetsumHT3R6aaa(0x0),
-fh2JetsumHT3R6aaap(0x0),
-fh2JetsumHT3R8aaa(0x0),
-fh2JetsumHT3R8aaap(0x0),
-fh2JetsumHT3R10aaa(0x0),
-fh2JetsumHT3R10aaap(0x0),
-fh2JetsumHT3R2b(0x0),
-fh2JetsumHT3R2bp(0x0),
-fh2JetsumHT3R4b(0x0),
-fh2JetsumHT3R4bp(0x0),
-fh2JetsumHT3R6b(0x0),
-fh2JetsumHT3R6bp(0x0),
-fh2JetsumHT3R8b(0x0),
-fh2JetsumHT3R8bp(0x0),
-fh2JetsumHT3R10b(0x0),
-fh2JetsumHT3R10bp(0x0),
-fh2JetsumHT3R2bb(0x0),
-fh2JetsumHT3R2bbp(0x0),
-fh2JetsumHT3R4bb(0x0),
-fh2JetsumHT3R4bbp(0x0),
-fh2JetsumHT3R6bb(0x0),
-fh2JetsumHT3R6bbp(0x0),
-fh2JetsumHT3R8bb(0x0),
-fh2JetsumHT3R8bbp(0x0),
-fh2JetsumHT3R10bb(0x0),
-fh2JetsumHT3R10bbp(0x0),
-fh2JetsumHT3R2bbb(0x0),
-fh2JetsumHT3R2bbbp(0x0),
-fh2JetsumHT3R4bbb(0x0),
-fh2JetsumHT3R4bbbp(0x0),
-fh2JetsumHT3R6bbb(0x0),
-fh2JetsumHT3R6bbbp(0x0),
-fh2JetsumHT3R8bbb(0x0),
-fh2JetsumHT3R8bbbp(0x0),
-fh2JetsumHT3R10bbb(0x0),
-fh2JetsumHT3R10bbbp(0x0),
 fh2Ntriggers(0x0),
 fh2JetDensity(0x0),
 fh2JetDensityA4(0x0),
-fh3spectriggeredC10(0x0),
+fh2RPJets(0x0),
+fh3spectriggeredC4080(0x0),
 fh3spectriggeredC20(0x0),
-fh3spectriggeredC3060(0x0),
-fh3specbiased(0x0),
-fh3spectot(0x0),
-fh3spectotb(0x0)
+fh3spectriggeredC3060(0x0)
+
  
 {
    // default Constructor
@@ -243,6 +186,9 @@ fAngStructCloseTracks(0),
 fCheckMethods(0),
 fDoEventMixing(0),
 fFlagPhiBkg(0),
+fFlagRandom(0),
+fRPAngle(0),
+fNRPBins(3),
 fJetEtaMin(-.5),
 fJetEtaMax(.5),
 fNevents(0),
@@ -270,8 +216,9 @@ fh2JetCoreMethod1C30(0x0),
 fh2JetCoreMethod2C30(0x0),
 fh2JetCoreMethod1C60(0x0),
 fh2JetCoreMethod2C60(0x0),
-fh3JetTrackC10(0x0),
+fh3JetTrackC3060(0x0),
 fh3JetTrackC20(0x0),
+fh3JetTrackC4080(0x0), 
 fh2AngStructpt1C10(0x0),
 fh2AngStructpt2C10(0x0),
 fh2AngStructpt3C10(0x0),
@@ -288,75 +235,14 @@ fh2AngStructpt1C60(0x0),
 fh2AngStructpt2C60(0x0),
 fh2AngStructpt3C60(0x0),
 fh2AngStructpt4C60(0x0),    
-fh2JetsumHT3R2a(0x0),
-fh2JetsumHT3R2ap(0x0),
-fh2JetsumHT3R4a(0x0),
-fh2JetsumHT3R4ap(0x0),
-fh2JetsumHT3R6a(0x0),
-fh2JetsumHT3R6ap(0x0),
-fh2JetsumHT3R8a(0x0),
-fh2JetsumHT3R8ap(0x0),
-fh2JetsumHT3R10a(0x0),
-fh2JetsumHT3R10ap(0x0),
-fh2JetsumHT3R2aa(0x0),
-fh2JetsumHT3R2aap(0x0),
-fh2JetsumHT3R4aa(0x0),
-fh2JetsumHT3R4aap(0x0),
-fh2JetsumHT3R6aa(0x0),
-fh2JetsumHT3R6aap(0x0),
-fh2JetsumHT3R8aa(0x0),
-fh2JetsumHT3R8aap(0x0),
-fh2JetsumHT3R10aa(0x0),
-fh2JetsumHT3R10aap(0x0),
-fh2JetsumHT3R2aaa(0x0),
-fh2JetsumHT3R2aaap(0x0),
-fh2JetsumHT3R4aaa(0x0),
-fh2JetsumHT3R4aaap(0x0),
-fh2JetsumHT3R6aaa(0x0),
-fh2JetsumHT3R6aaap(0x0),
-fh2JetsumHT3R8aaa(0x0),
-fh2JetsumHT3R8aaap(0x0),
-fh2JetsumHT3R10aaa(0x0),
-fh2JetsumHT3R10aaap(0x0),
-fh2JetsumHT3R2b(0x0),
-fh2JetsumHT3R2bp(0x0),
-fh2JetsumHT3R4b(0x0),
-fh2JetsumHT3R4bp(0x0),
-fh2JetsumHT3R6b(0x0),
-fh2JetsumHT3R6bp(0x0),
-fh2JetsumHT3R8b(0x0),
-fh2JetsumHT3R8bp(0x0),
-fh2JetsumHT3R10b(0x0),
-fh2JetsumHT3R10bp(0x0),
-fh2JetsumHT3R2bb(0x0),
-fh2JetsumHT3R2bbp(0x0),
-fh2JetsumHT3R4bb(0x0),
-fh2JetsumHT3R4bbp(0x0),
-fh2JetsumHT3R6bb(0x0),
-fh2JetsumHT3R6bbp(0x0),
-fh2JetsumHT3R8bb(0x0),
-fh2JetsumHT3R8bbp(0x0),
-fh2JetsumHT3R10bb(0x0),
-fh2JetsumHT3R10bbp(0x0),
-fh2JetsumHT3R2bbb(0x0),
-fh2JetsumHT3R2bbbp(0x0),
-fh2JetsumHT3R4bbb(0x0),
-fh2JetsumHT3R4bbbp(0x0),
-fh2JetsumHT3R6bbb(0x0),
-fh2JetsumHT3R6bbbp(0x0),
-fh2JetsumHT3R8bbb(0x0),
-fh2JetsumHT3R8bbbp(0x0),
-fh2JetsumHT3R10bbb(0x0),
-fh2JetsumHT3R10bbbp(0x0),
 fh2Ntriggers(0x0),
 fh2JetDensity(0x0),
 fh2JetDensityA4(0x0),
-fh3spectriggeredC10(0x0),
+fh2RPJets(0x0),
+fh3spectriggeredC4080(0x0),
 fh3spectriggeredC20(0x0),
-fh3spectriggeredC3060(0x0),
-fh3specbiased(0x0),
-fh3spectot(0x0),
-fh3spectotb(0x0)
+fh3spectriggeredC3060(0x0)
+
  {
    // Constructor
 
@@ -468,9 +354,9 @@ void AliAnalysisTaskJetCore::UserCreateOutputObjects()
     fh2JetCoreMethod1C60 = new TH2F("JetCoreMethod1C60","",150, 0., 150.,100, 0., 1.5);
     fh2JetCoreMethod2C60 = new TH2F("JetCoreMethod2C60","",150, 0., 150.,100, 0., 1.5);}
 
-    fh3JetTrackC10=new TH3F("JetTrackC10","",50,0,50,150,0.,150.,32,0,3.2);
-    fh3JetTrackC20=new TH3F("JetTrackC20","",50,0,50,150,0.,150.,32,0,3.2);
-   
+    fh3JetTrackC3060=new TH3F("JetTrackC3060","",50,0,50,150,0.,150.,34,0,3.4);
+    fh3JetTrackC20=new TH3F("JetTrackC20","",50,0,50,150,0.,150.,34,0,3.4);
+    fh3JetTrackC4080=new TH3F("JetTrackC4080","",50,0,50,150,0.,150.,34,0,3.4);
     if(fAngStructCloseTracks>0){
     fh2AngStructpt1C10 = new TH2F("Ang struct pt1 C10","",15,0.,1.5,150,0.,10.);
     fh2AngStructpt2C10 = new TH2F("Ang struct pt2 C10","",15,0.,1.5,150,0.,10.);
@@ -489,84 +375,18 @@ void AliAnalysisTaskJetCore::UserCreateOutputObjects()
     fh2AngStructpt3C60 = new TH2F("Ang struct pt3 C60","",15,0.,1.5,150,0.,10.);
     fh2AngStructpt4C60 = new TH2F("Ang struct pt4 C60","",15,0.,1.5,150,0.,10.); }
 
-    fh2JetsumHT3R2a = new TH2F("Pt sum R02 HT0 TT10","",20,0.,200.,100,0.,10.);
-    fh2JetsumHT3R2ap = new TH2F("Pt sum R02 HT0 TT10 p","",20,0.,200.,100,0.,10.);
-    fh2JetsumHT3R4a = new TH2F("Pt sum R04 HT0 TT10","",20,0.,200.,100,0.,10.);
-    fh2JetsumHT3R4ap = new TH2F("Pt sum R04 HT0 TT10 p","",20,0.,200.,100,0.,10.);
-    fh2JetsumHT3R6a = new TH2F("Pt sum R06 HT0 TT10","",20,0.,200.,100,0.,10.);
-    fh2JetsumHT3R6ap = new TH2F("Pt sum R06 HT0 TT10 p","",20,0.,200.,100,0.,10.);
-    fh2JetsumHT3R8a = new TH2F("Pt sum R08 HT0 TT10","",20,0.,200.,100,0.,10.);
-    fh2JetsumHT3R8ap = new TH2F("Pt sum R08 HT0 TT10 p","",20,0.,200.,100,0.,10.);
-    fh2JetsumHT3R10a = new TH2F("Pt sum R10 HT0 TT10","",20,0.,200.,100,0.,10.);
-    fh2JetsumHT3R10ap = new TH2F("Pt sum R10 HT0 TT10 p","",20,0.,200.,100,0.,10.);  
-
-    fh2JetsumHT3R2aa = new TH2F("Pt sum R02 HT0 TT20","",20,0.,200.,100,0.,10.);
-    fh2JetsumHT3R2aap = new TH2F("Pt sum R02 HT0 TT20 p","",20,0.,200.,100,0.,10.);
-    fh2JetsumHT3R4aa = new TH2F("Pt sum R04 HT0 TT20","",20,0.,200.,100,0.,10.);
-    fh2JetsumHT3R4aap = new TH2F("Pt sum R04 HT0 TT20 p","",20,0.,200.,100,0.,10.);
-    fh2JetsumHT3R6aa = new TH2F("Pt sum R06 HT0 TT20","",20,0.,200.,100,0.,10.);
-    fh2JetsumHT3R6aap = new TH2F("Pt sum R06 HT0 TT20 p","",20,0.,200.,100,0.,10.);
-    fh2JetsumHT3R8aa = new TH2F("Pt sum R08 HT0 TT20","",20,0.,200.,100,0.,10.);
-    fh2JetsumHT3R8aap = new TH2F("Pt sum R08 HT0 TT20 p","",20,0.,200.,100,0.,10.);
-    fh2JetsumHT3R10aa = new TH2F("Pt sum R10 HT0 TT20","",20,0.,200.,100,0.,10.);
-    fh2JetsumHT3R10aap = new TH2F("Pt sum R10 HT0 TT20 p","",20,0.,200.,100,0.,10.);  
-
-    fh2JetsumHT3R2aaa = new TH2F("Pt sum R02 HT0 TT0","",20,0.,200.,100,0.,10.);
-    fh2JetsumHT3R2aaap = new TH2F("Pt sum R02 HT0 TT0 p","",20,0.,200.,100,0.,10.);
-    fh2JetsumHT3R4aaa = new TH2F("Pt sum R04 HT0 TT0","",20,0.,200.,100,0.,10.);
-    fh2JetsumHT3R4aaap = new TH2F("Pt sum R04 HT0 TT0 p","",20,0.,200.,100,0.,10.);
-    fh2JetsumHT3R6aaa = new TH2F("Pt sum R06 HT0 TT0","",20,0.,200.,100,0.,10.);
-    fh2JetsumHT3R6aaap = new TH2F("Pt sum R06 HT0 TT0 p","",20,0.,200.,100,0.,10.);
-    fh2JetsumHT3R8aaa = new TH2F("Pt sum R08 HT0 TT0","",20,0.,200.,100,0.,10.);
-    fh2JetsumHT3R8aaap = new TH2F("Pt sum R08 HT0 TT0 p","",20,0.,200.,100,0.,10.);
-    fh2JetsumHT3R10aaa = new TH2F("Pt sum R10 HT0 TT0","",20,0.,200.,100,0.,10.);
-    fh2JetsumHT3R10aaap = new TH2F("Pt sum R10 HT0 TT0 p","",20,0.,200.,100,0.,10.);  
-
-    fh2JetsumHT3R2b = new TH2F("Pt sum R02 HT6 TT10","",20,0.,200.,100,0.,10.);
-    fh2JetsumHT3R2bp = new TH2F("Pt sum R02 HT6 TT10 p","",20,0.,200.,100,0.,10.);
-    fh2JetsumHT3R4b = new TH2F("Pt sum R04 HT6 TT10","",20,0.,200.,100,0.,10.);
-    fh2JetsumHT3R4bp = new TH2F("Pt sum R04 HT6 TT10 p","",20,0.,200.,100,0.,10.);
-    fh2JetsumHT3R6b = new TH2F("Pt sum R06 HT6 TT10","",20,0.,200.,100,0.,10.);
-    fh2JetsumHT3R6bp = new TH2F("Pt sum R06 HT6 TT10 p","",20,0.,200.,100,0.,10.);
-    fh2JetsumHT3R8b = new TH2F("Pt sum R08 HT6 TT10","",20,0.,200.,100,0.,10.);
-    fh2JetsumHT3R8bp = new TH2F("Pt sum R08 HT6 TT10 p","",20,0.,200.,100,0.,10.);
-    fh2JetsumHT3R10b = new TH2F("Pt sum R10 HT6 TT10","",20,0.,200.,100,0.,10.);
-    fh2JetsumHT3R10bp = new TH2F("Pt sum R10 HT6 TT10 p","",20,0.,200.,100,0.,10.);  
-
-    fh2JetsumHT3R2bb = new TH2F("Pt sum R02 HT6 TT20","",20,0.,200.,100,0.,10.);
-    fh2JetsumHT3R2bbp = new TH2F("Pt sum R02 HT6 TT20 p","",20,0.,200.,100,0.,10.);
-    fh2JetsumHT3R4bb = new TH2F("Pt sum R04 HT6 TT20","",20,0.,200.,100,0.,10.);
-    fh2JetsumHT3R4bbp = new TH2F("Pt sum R04 HT6 TT20 p","",20,0.,200.,100,0.,10.);
-    fh2JetsumHT3R6bb = new TH2F("Pt sum R06 HT6 TT20","",20,0.,200.,100,0.,10.);
-    fh2JetsumHT3R6bbp = new TH2F("Pt sum R06 HT6 TT20 p","",20,0.,200.,100,0.,10.);
-    fh2JetsumHT3R8bb = new TH2F("Pt sum R08 HT6 TT20","",20,0.,200.,100,0.,10.);
-    fh2JetsumHT3R8bbp = new TH2F("Pt sum R08 HT6 TT20 p","",20,0.,200.,100,0.,10.);
-    fh2JetsumHT3R10bb = new TH2F("Pt sum R10 HT6 TT20","",20,0.,200.,100,0.,10.);
-    fh2JetsumHT3R10bbp = new TH2F("Pt sum R10 HT6 TT20 p","",20,0.,200.,100,0.,10.);  
-
-    fh2JetsumHT3R2bbb = new TH2F("Pt sum R02 HT6 TT0","",20,0.,200.,100,0.,10.);
-    fh2JetsumHT3R2bbbp = new TH2F("Pt sum R02 HT6 TT0 p","",20,0.,200.,100,0.,10.);
-    fh2JetsumHT3R4bbb = new TH2F("Pt sum R04 HT6 TT0","",20,0.,200.,100,0.,10.);
-    fh2JetsumHT3R4bbbp = new TH2F("Pt sum R04 HT6 TT0 p","",20,0.,200.,100,0.,10.);
-    fh2JetsumHT3R6bbb = new TH2F("Pt sum R06 HT6 TT0","",20,0.,200.,100,0.,10.);
-    fh2JetsumHT3R6bbbp = new TH2F("Pt sum R06 HT6 TT0 p","",20,0.,200.,100,0.,10.);
-    fh2JetsumHT3R8bbb = new TH2F("Pt sum R08 HT6 TT0","",20,0.,200.,100,0.,10.);
-    fh2JetsumHT3R8bbbp = new TH2F("Pt sum R08 HT6 TT0 p","",20,0.,200.,100,0.,10.);
-    fh2JetsumHT3R10bbb = new TH2F("Pt sum R10 HT6 TT0","",20,0.,200.,100,0.,10.);
-    fh2JetsumHT3R10bbbp = new TH2F("Pt sum R10 HT6 TT0 p","",20,0.,200.,100,0.,10.);  
-
    
 
     fh2Ntriggers=new TH2F("# of triggers","",10,0.,100.,50,0.,50.);
     fh2JetDensity=new TH2F("Jet density vs centrality A>0.4","",100,0.,4000.,100,0.,5.);
     fh2JetDensityA4=new TH2F("Jet density vs multiplicity A>0.4","",100,0.,4000.,100,0.,5.);
-    fh3spectriggeredC10 = new TH3F("Triggered spectrumC10","",5,0.,1.,140,-80.,200.,50,0.,50.);
+    fh2RPJets=new TH2F("RPJet","",3,0.,3.,150,0.,150.); 
+    fh3spectriggeredC4080 = new TH3F("Triggered spectrumC4080","",5,0.,1.,140,-80.,200.,50,0.,50.);
     fh3spectriggeredC20 = new TH3F("Triggered spectrumC20","",5,0.,1.,140,-80.,200.,50,0.,50.);
     fh3spectriggeredC3060 = new TH3F("Triggered spectrumC3060","",5,0.,1.,140,-80.,200.,50,0.,50.);
 
-    fh3specbiased = new TH3F("Biased spectrum","",10,0,100,50,0.,200.,50,0.,50.);
-    fh3spectot = new TH3F("Total spectrum 0-10","",50,0.,200.,50,0.,50.,50,0.,50.);
-    fh3spectotb = new TH3F("Total spectrum 30-60","",50,0.,200.,50,0.,50.,50,0.,50.);    
+    
+    
    fOutputList->Add(fHistEvtSelection);
 
    fOutputList->Add(fhnDeltaR);
@@ -586,8 +406,9 @@ void AliAnalysisTaskJetCore::UserCreateOutputObjects()
       fOutputList->Add(fh2JetCoreMethod1C60);
       fOutputList->Add(fh2JetCoreMethod2C60);}
       
-      fOutputList->Add(fh3JetTrackC10);
+      fOutputList->Add(fh3JetTrackC3060);
       fOutputList->Add(fh3JetTrackC20);
+      fOutputList->Add(fh3JetTrackC4080); 
             
      
 
@@ -611,80 +432,18 @@ void AliAnalysisTaskJetCore::UserCreateOutputObjects()
        fOutputList->Add(fh2AngStructpt4C60);}  
 
 
-       fOutputList->Add(fh2JetsumHT3R2a);
-        fOutputList->Add(fh2JetsumHT3R2ap);
-               fOutputList->Add(fh2JetsumHT3R4a);
-        fOutputList->Add(fh2JetsumHT3R4ap);
-        fOutputList->Add(fh2JetsumHT3R6a);
-        fOutputList->Add(fh2JetsumHT3R6ap);
-               fOutputList->Add(fh2JetsumHT3R8a);
-        fOutputList->Add(fh2JetsumHT3R8ap);
-        fOutputList->Add(fh2JetsumHT3R10a);
-        fOutputList->Add(fh2JetsumHT3R10ap);
-       fOutputList->Add(fh2JetsumHT3R2aa);
-        fOutputList->Add(fh2JetsumHT3R2aap);
-               fOutputList->Add(fh2JetsumHT3R4aa);
-        fOutputList->Add(fh2JetsumHT3R4aap);
-        fOutputList->Add(fh2JetsumHT3R6aa);
-        fOutputList->Add(fh2JetsumHT3R6aap);
-               fOutputList->Add(fh2JetsumHT3R8aa);
-        fOutputList->Add(fh2JetsumHT3R8aap);
-        fOutputList->Add(fh2JetsumHT3R10aa);
-        fOutputList->Add(fh2JetsumHT3R10aap);
-        fOutputList->Add(fh2JetsumHT3R2aaa);
-        fOutputList->Add(fh2JetsumHT3R2aaap);
-               fOutputList->Add(fh2JetsumHT3R4aaa);
-        fOutputList->Add(fh2JetsumHT3R4aaap);
-        fOutputList->Add(fh2JetsumHT3R6aaa);
-        fOutputList->Add(fh2JetsumHT3R6aaap);
-               fOutputList->Add(fh2JetsumHT3R8aaa);
-        fOutputList->Add(fh2JetsumHT3R8aaap);
-        fOutputList->Add(fh2JetsumHT3R10aaa);
-        fOutputList->Add(fh2JetsumHT3R10aaap);
-
-        fOutputList->Add(fh2JetsumHT3R2b);
-        fOutputList->Add(fh2JetsumHT3R2bp);
-               fOutputList->Add(fh2JetsumHT3R4b);
-        fOutputList->Add(fh2JetsumHT3R4bp);
-        fOutputList->Add(fh2JetsumHT3R6b);
-        fOutputList->Add(fh2JetsumHT3R6bp);
-               fOutputList->Add(fh2JetsumHT3R8b);
-        fOutputList->Add(fh2JetsumHT3R8bp);
-        fOutputList->Add(fh2JetsumHT3R10b);
-        fOutputList->Add(fh2JetsumHT3R10bp);
-       fOutputList->Add(fh2JetsumHT3R2bb);
-        fOutputList->Add(fh2JetsumHT3R2bbp);
-               fOutputList->Add(fh2JetsumHT3R4bb);
-        fOutputList->Add(fh2JetsumHT3R4bbp);
-        fOutputList->Add(fh2JetsumHT3R6bb);
-        fOutputList->Add(fh2JetsumHT3R6bbp);
-               fOutputList->Add(fh2JetsumHT3R8bb);
-        fOutputList->Add(fh2JetsumHT3R8bbp);
-        fOutputList->Add(fh2JetsumHT3R10bb);
-        fOutputList->Add(fh2JetsumHT3R10bbp);
-        fOutputList->Add(fh2JetsumHT3R2bbb);
-        fOutputList->Add(fh2JetsumHT3R2bbbp);
-               fOutputList->Add(fh2JetsumHT3R4bbb);
-        fOutputList->Add(fh2JetsumHT3R4bbbp);
-        fOutputList->Add(fh2JetsumHT3R6bbb);
-        fOutputList->Add(fh2JetsumHT3R6bbbp);
-               fOutputList->Add(fh2JetsumHT3R8bbb);
-        fOutputList->Add(fh2JetsumHT3R8bbbp);
-        fOutputList->Add(fh2JetsumHT3R10bbb);
-        fOutputList->Add(fh2JetsumHT3R10bbbp);
+
 
  
        fOutputList->Add(fh2Ntriggers);
         fOutputList->Add(fh2JetDensity);
         fOutputList->Add(fh2JetDensityA4);
-
-       fOutputList->Add(fh3spectriggeredC10);
+        fOutputList->Add(fh2RPJets);
+       fOutputList->Add(fh3spectriggeredC4080);
        fOutputList->Add(fh3spectriggeredC20); 
        fOutputList->Add(fh3spectriggeredC3060);   
 
-       fOutputList->Add(fh3specbiased);
-       fOutputList->Add(fh3spectot);
-       fOutputList->Add(fh3spectotb); 
+     
    // =========== Switch on Sumw2 for all histos ===========
    for (Int_t i=0; i<fOutputList->GetEntries(); ++i) {
       TH1 *h1 = dynamic_cast<TH1*>(fOutputList->At(i));
@@ -823,8 +582,11 @@ void AliAnalysisTaskJetCore::UserExec(Option_t *)
     }
    
    Float_t rho = 0;
-   if(externalBackground)rho = externalBackground->GetBackground(0);
 
+   if(fFlagRandom==0){
+     if(externalBackground)rho = externalBackground->GetBackground(0);}
+   if(fFlagRandom==1){
+      if(externalBackground)rho = externalBackground->GetBackground(2);}
 
    // fetch jets
    TClonesArray *aodJets[2];
@@ -872,19 +634,12 @@ void AliAnalysisTaskJetCore::UserExec(Option_t *)
    Double_t phismall=0.;
          
   
-   // Double_t up1[15]={0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.};
-   // Double_t up2[15]={0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.};
-   // Double_t up3[15]={0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.};
-   // Double_t up4[15]={0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.};
-   // Double_t down1[15]={0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.};
-   // Double_t down2[15]={0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.};
-   // Double_t down3[15]={0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.};
-   // Double_t down4[15]={0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.};
+
    Int_t iCount=0; 
    Int_t trigJet=-1;
    Int_t trigBBTrack=-1;
    Int_t trigInTrack=-1;
-     
+   fRPAngle = aod->GetHeader()->GetEventplane();     
 
    AliVParticle *partback = (AliVParticle*)ParticleList.At(nT);     
    if(!partback){  
@@ -900,21 +655,26 @@ void AliAnalysisTaskJetCore::UserExec(Option_t *)
            phibig  = jetbig->Phi();
            ptbig   = jetbig->Pt();
            if(ptbig==0) continue; 
+           Int_t phiBin = GetPhiBin(phibig-fRPAngle);       
            areabig = jetbig->EffectiveAreaCharged();
            Double_t ptcorr=ptbig-rho*areabig;
           if((etabig<fJetEtaMin)||(etabig>fJetEtaMax)) continue;
+       
            if(areabig>=0.2) injet=injet+1;
            if(areabig>=0.4) injet4=injet4+1;   
            Double_t dphi=RelativePhi(partback->Phi(),phibig);  
-           if(centValue<10.) fh3JetTrackC10->Fill(partback->Pt(),ptcorr,TMath::Abs(dphi));
+           if(centValue>40. && centValue<80.) fh3JetTrackC4080->Fill(partback->Pt(),ptcorr,TMath::Abs(dphi));
            if(centValue<20.) fh3JetTrackC20->Fill(partback->Pt(),ptcorr,TMath::Abs(dphi));
+           if(centValue>30. && centValue<60.) fh3JetTrackC3060->Fill(partback->Pt(),ptcorr,TMath::Abs(dphi));
+
            if(fFlagPhiBkg!=0) if((TMath::Abs(dphi)<TMath::Pi()/2.-0.1)||(TMath::Abs(dphi)>TMath::Pi()/2.+0.1)) continue;
            if(fFlagPhiBkg==0) if(TMath::Abs(dphi)<TMath::Pi()-0.6) continue;
+                   if(centValue<20.) fh2RPJets->Fill(phiBin, ptcorr);
                    Double_t dismin=100.;
                    Double_t ptmax=-10.; 
                    Int_t index1=-1;
                    Int_t index2=-1;
-           if(centValue<10.)  fh3spectriggeredC10->Fill(jetbig->EffectiveAreaCharged(),ptcorr,partback->Pt());
+           if(centValue>40. && centValue<80.)  fh3spectriggeredC4080->Fill(jetbig->EffectiveAreaCharged(),ptcorr,partback->Pt());
            if(centValue<20.)  fh3spectriggeredC20->Fill(jetbig->EffectiveAreaCharged(),ptcorr,partback->Pt());
            if(centValue>30. && centValue<60.)  fh3spectriggeredC3060->Fill(jetbig->EffectiveAreaCharged(),ptcorr,partback->Pt());
 
@@ -931,9 +691,7 @@ void AliAnalysisTaskJetCore::UserExec(Option_t *)
                         ippt=ir;}}
                         leadtrack=(AliAODTrack*)(genTrackList->At(ippt));
                         if(!leadtrack) continue;
-                        fh3specbiased->Fill(centValue,ptcorr,leadtrack->Pt());
-                        if(centValue<10)fh3spectot->Fill(ptcorr,leadtrack->Pt(),partback->Pt());
-                        if(centValue>30. && centValue<60.)fh3spectotb->Fill(ptcorr,leadtrack->Pt(),partback->Pt()); 
+                       
                        //store one trigger info                   
                         if(iCount==0){                        
                          trigJet=i;
@@ -941,7 +699,7 @@ void AliAnalysisTaskJetCore::UserExec(Option_t *)
                           trigInTrack=ippt;
                           iCount=iCount+1;} 
 
-
+   
                  if(fCheckMethods){
                   for(Int_t j=0; j<fListJets[1]->GetEntries(); ++j){
                   AliAODJet* jetsmall = (AliAODJet*)(fListJets[1]->At(j));
@@ -969,200 +727,23 @@ void AliAnalysisTaskJetCore::UserExec(Option_t *)
                   if((centValue>20)&&(centValue<40)) fh2JetCoreMethod2C20->Fill(ptcorr,jetmethod2->Pt()/ptbig); 
                  if((centValue>30)&&(centValue<60)) fh2JetCoreMethod2C30->Fill(ptcorr,jetmethod2->Pt()/ptbig);
                  if(centValue>60) fh2JetCoreMethod2C60->Fill(ptcorr,jetmethod2->Pt()/ptbig); }}  
-                 Double_t sumpt2a=0.;
-                  Double_t sumpt2b=0.;
-                  Double_t sumpt4a=0.;
-                  Double_t sumpt4b=0.;
-                  Double_t sumpt6a=0.;
-                  Double_t sumpt6b=0.;
-                  Double_t sumpt8a=0.;
-                  Double_t sumpt8b=0.;
-                  Double_t sumpt10a=0.;
-                  Double_t sumpt10b=0.;
-                  Double_t sumpt2aaa=0.;
-                  Double_t sumpt2bbb=0.;
-                  Double_t sumpt4aaa=0.;
-                  Double_t sumpt4bbb=0.;
-                  Double_t sumpt6aaa=0.;
-                  Double_t sumpt6bbb=0.;
-                  Double_t sumpt8aaa=0.;
-                  Double_t sumpt8bbb=0.;
-                  Double_t sumpt10aaa=0.;
-                  Double_t sumpt10bbb=0.; 
-                  Double_t sumpt2ap=0.;
-                  Double_t sumpt2bp=0.;
-                  Double_t sumpt4ap=0.;
-                  Double_t sumpt4bp=0.;
-                  Double_t sumpt6ap=0.;
-                  Double_t sumpt6bp=0.;
-                  Double_t sumpt8ap=0.;
-                  Double_t sumpt8bp=0.;
-                  Double_t sumpt10ap=0.;
-                  Double_t sumpt10bp=0.;
-                  Double_t sumpt2aaap=0.;
-                  Double_t sumpt2bbbp=0.;
-                  Double_t sumpt4aaap=0.;
-                  Double_t sumpt4bbbp=0.;
-                  Double_t sumpt6aaap=0.;
-                  Double_t sumpt6bbbp=0.;
-                  Double_t sumpt8aaap=0.;
-                  Double_t sumpt8bbbp=0.;
-                  Double_t sumpt10aaap=0.;
-                  Double_t sumpt10bbbp=0.; 
-
-
-
+       
                   
-
-                 for(int it = 0;it<ParticleList.GetEntries();++it){
+         if(fDoEventMixing==0){
+        for(int it = 0;it<ParticleList.GetEntries();++it){
          AliVParticle *part = (AliVParticle*)ParticleList.At(it);
                  Double_t deltaR = jetbig->DeltaR(part);
           Double_t deltaEta = etabig-part->Eta();
-          if(partback->Pt()>14.){
-
-          if(centValue<10){
-           //for one more centrality
-               if(deltaR<0.2){if(leadtrack->Pt()>0.)sumpt2a=sumpt2a+part->Pt();
-                             if(leadtrack->Pt()>6) sumpt2b=sumpt2b+part->Pt();}
-              if(deltaR>=0.2 && deltaR<0.4){if(leadtrack->Pt()>0.)sumpt4a=sumpt4a+part->Pt();
-                                              if(leadtrack->Pt()>6.)sumpt4b=sumpt4b+part->Pt();}
-
-               if(deltaR>=0.4 && deltaR<0.6){if(leadtrack->Pt()>0.)sumpt6a=sumpt6a+part->Pt();
-                                            if(leadtrack->Pt()>6.) sumpt6b=sumpt6b+part->Pt();} 
-               if(deltaR>=0.6 && deltaR<0.8){if(leadtrack->Pt()>0.)sumpt8a=sumpt8a+part->Pt();
-                                            if(leadtrack->Pt()>6.) sumpt8b=sumpt8b+part->Pt();}                               
-               if(deltaR>=0.8 && deltaR<1.2){if(leadtrack->Pt()>0.)sumpt10a=sumpt10a+part->Pt();
-                if(leadtrack->Pt()>6.) sumpt10b=sumpt10b+part->Pt();}}
-
-               if(centValue>30. && centValue<60.){
-           //for one more centrality
-               if(deltaR<0.2){if(leadtrack->Pt()>0.)sumpt2ap=sumpt2ap+part->Pt();
-                             if(leadtrack->Pt()>6) sumpt2bp=sumpt2bp+part->Pt();}
-              if(deltaR>=0.2 && deltaR<0.4){if(leadtrack->Pt()>0.)sumpt4ap=sumpt4ap+part->Pt();
-                                              if(leadtrack->Pt()>6.)sumpt4bp=sumpt4bp+part->Pt();}
-
-               if(deltaR>=0.4 && deltaR<0.6){if(leadtrack->Pt()>0.)sumpt6ap=sumpt6ap+part->Pt();
-                                            if(leadtrack->Pt()>6.) sumpt6bp=sumpt6bp+part->Pt();} 
-               if(deltaR>=0.6 && deltaR<0.8){if(leadtrack->Pt()>0.)sumpt8ap=sumpt8ap+part->Pt();
-                                            if(leadtrack->Pt()>6.) sumpt8bp=sumpt8bp+part->Pt();}                               
-               if(deltaR>=0.8 && deltaR<1.2){if(leadtrack->Pt()>0.)sumpt10ap=sumpt10ap+part->Pt();
-                if(leadtrack->Pt()>6.) sumpt10bp=sumpt10bp+part->Pt();}}}
-
-              
-
-             if(partback->Pt()<1.){
-
-              if(centValue<10){
-           //for one more centrality
-               if(deltaR<0.2){if(leadtrack->Pt()>0.)sumpt2aaa=sumpt2aaa+part->Pt();
-                             if(leadtrack->Pt()>6.) sumpt2bbb=sumpt2bbb+part->Pt();}
-              if(deltaR>=0.2 && deltaR<0.4){if(leadtrack->Pt()>0.)sumpt4aaa=sumpt4aaa+part->Pt();
-                                              if(leadtrack->Pt()>6.)sumpt4bbb=sumpt4bbb+part->Pt();}
-
-               if(deltaR>=0.4 && deltaR<0.6){if(leadtrack->Pt()>0.)sumpt6aaa=sumpt6aaa+part->Pt();
-                                            if(leadtrack->Pt()>6.) sumpt6bbb=sumpt6bbb+part->Pt();} 
-               if(deltaR>=0.6 && deltaR<0.8){if(leadtrack->Pt()>0.)sumpt8aaa=sumpt8aaa+part->Pt();
-                                            if(leadtrack->Pt()>6.) sumpt8bbb=sumpt8bbb+part->Pt();}                               
-               if(deltaR>=0.8 && deltaR<1.2){if(leadtrack->Pt()>0.)sumpt10aaa=sumpt10aaa+part->Pt();
-                if(leadtrack->Pt()>6.) sumpt10bbb=sumpt10bbb+part->Pt();}}
-
-               if(centValue>30. && centValue<60.){
-           //for one more centrality
-               if(deltaR<0.2){if(leadtrack->Pt()>0.)sumpt2aaap=sumpt2aaap+part->Pt();
-                             if(leadtrack->Pt()>6) sumpt2bbbp=sumpt2bbbp+part->Pt();}
-              if(deltaR>=0.2 && deltaR<0.4){if(leadtrack->Pt()>0.)sumpt4aaap=sumpt4aaap+part->Pt();
-                                              if(leadtrack->Pt()>6.)sumpt4bbbp=sumpt4bbbp+part->Pt();}
-
-               if(deltaR>=0.4 && deltaR<0.6){if(leadtrack->Pt()>0.)sumpt6aaap=sumpt6aaap+part->Pt();
-                                            if(leadtrack->Pt()>6.) sumpt6bbbp=sumpt6bbbp+part->Pt();} 
-               if(deltaR>=0.6 && deltaR<0.8){if(leadtrack->Pt()>0.)sumpt8aaap=sumpt8aaap+part->Pt();
-                                            if(leadtrack->Pt()>6.) sumpt8bbbp=sumpt8bbbp+part->Pt();}                               
-               if(deltaR>=0.8 && deltaR<1.2){if(leadtrack->Pt()>0.)sumpt10aaap=sumpt10aaap+part->Pt();
-                if(leadtrack->Pt()>6.) sumpt10bbbp=sumpt10bbbp+part->Pt();}}}
-
  
-
-
-         if(fDoEventMixing==0){
+         
           Double_t deltaPhi=phibig-part->Phi();
           if(deltaPhi<-0.5*TMath::Pi()) deltaPhi+=2.*TMath::Pi();
           if(deltaPhi>3./2.*TMath::Pi()) deltaPhi-=2.*TMath::Pi();
                  Double_t jetEntries[8] = {centValue,ptcorr,part->Pt(),deltaR,deltaEta,deltaPhi,leadtrack->Pt(),partback->Pt()};                     fhnDeltaR->Fill(jetEntries);}
-
-                 }
-          //end of track loop
-            Double_t rhoin2=rho*TMath::Pi()*0.2*0.2; 
-            Double_t rhoin4=rho*TMath::Pi()*(0.4*0.4-0.2*0.2);
-            Double_t rhoin6=rho*TMath::Pi()*(0.6*0.6-0.4*0.4);
-            Double_t rhoin8=rho*TMath::Pi()*(0.8*0.8-0.6*0.6);
-            Double_t rhoin10=rho*TMath::Pi()*(1.2*1.2-0.8*0.8); 
-           
-              
-          if(rho!=0){
-
-          if(partback->Pt()>14.){
-           if(centValue<10.){
-           if(leadtrack->Pt()>0.){
-          fh2JetsumHT3R2a->Fill(ptcorr,sumpt2a/rhoin2);
-          fh2JetsumHT3R4a->Fill(ptcorr,sumpt4a/rhoin4);
-          fh2JetsumHT3R6a->Fill(ptcorr,sumpt6a/rhoin6);
-          fh2JetsumHT3R8a->Fill(ptcorr,sumpt8a/rhoin8);
-          fh2JetsumHT3R10a->Fill(ptcorr,sumpt10a/rhoin10);}
-           if(leadtrack->Pt()>6.){
-          fh2JetsumHT3R2b->Fill(ptcorr,sumpt2b/rhoin2);
-          fh2JetsumHT3R4b->Fill(ptcorr,sumpt4b/rhoin4);
-          fh2JetsumHT3R6b->Fill(ptcorr,sumpt6b/rhoin6);
-          fh2JetsumHT3R8b->Fill(ptcorr,sumpt8b/rhoin8);
-          fh2JetsumHT3R10b->Fill(ptcorr,sumpt10b/rhoin10);}}
-
-                   if(centValue>30 && centValue<60.){
-            if(leadtrack->Pt()>0.){
-          fh2JetsumHT3R2ap->Fill(ptcorr,sumpt2ap/rhoin2);
-          fh2JetsumHT3R4ap->Fill(ptcorr,sumpt4ap/rhoin4);
-          fh2JetsumHT3R6ap->Fill(ptcorr,sumpt6ap/rhoin6);
-          fh2JetsumHT3R8ap->Fill(ptcorr,sumpt8ap/rhoin8);
-          fh2JetsumHT3R10ap->Fill(ptcorr,sumpt10ap/rhoin10);}
-           if(leadtrack->Pt()>6.){
-          fh2JetsumHT3R2bp->Fill(ptcorr,sumpt2bp/rhoin2);
-          fh2JetsumHT3R4bp->Fill(ptcorr,sumpt4bp/rhoin4);
-          fh2JetsumHT3R6bp->Fill(ptcorr,sumpt6bp/rhoin6);
-          fh2JetsumHT3R8bp->Fill(ptcorr,sumpt8bp/rhoin8);
-          fh2JetsumHT3R10bp->Fill(ptcorr,sumpt10bp/rhoin10);}}}
-
-
-
-                        if(partback->Pt()<1.){
-           if(centValue<10.){
-             if(leadtrack->Pt()>0.){
-          fh2JetsumHT3R2aaa->Fill(ptcorr,sumpt2aaa/rhoin2);
-          fh2JetsumHT3R4aaa->Fill(ptcorr,sumpt4aaa/rhoin4);
-          fh2JetsumHT3R6aaa->Fill(ptcorr,sumpt6aaa/rhoin6);
-          fh2JetsumHT3R8aaa->Fill(ptcorr,sumpt8aaa/rhoin8);
-          fh2JetsumHT3R10aaa->Fill(ptcorr,sumpt10aaa/rhoin10);}
-           if(leadtrack->Pt()>6.){
-          fh2JetsumHT3R2bbb->Fill(ptcorr,sumpt2bbb/rhoin2);
-          fh2JetsumHT3R4bbb->Fill(ptcorr,sumpt4bbb/rhoin4);
-          fh2JetsumHT3R6bbb->Fill(ptcorr,sumpt6bbb/rhoin6);
-          fh2JetsumHT3R8bbb->Fill(ptcorr,sumpt8bbb/rhoin8);
-          fh2JetsumHT3R10bbb->Fill(ptcorr,sumpt10bbb/rhoin10);}}
-
-                   if(centValue>30 && centValue<60.){
-            if(leadtrack->Pt()>0.){
-          fh2JetsumHT3R2aaap->Fill(ptcorr,sumpt2aaap/rhoin2);
-          fh2JetsumHT3R4aaap->Fill(ptcorr,sumpt4aaap/rhoin4);
-          fh2JetsumHT3R6aaap->Fill(ptcorr,sumpt6aaap/rhoin6);
-          fh2JetsumHT3R8aaap->Fill(ptcorr,sumpt8aaap/rhoin8);
-          fh2JetsumHT3R10aaap->Fill(ptcorr,sumpt10aaap/rhoin10);}
-           if(leadtrack->Pt()>6.){
-          fh2JetsumHT3R2bbbp->Fill(ptcorr,sumpt2bbbp/rhoin2);
-          fh2JetsumHT3R4bbbp->Fill(ptcorr,sumpt4bbbp/rhoin4);
-          fh2JetsumHT3R6bbbp->Fill(ptcorr,sumpt6bbbp/rhoin6);
-          fh2JetsumHT3R8bbbp->Fill(ptcorr,sumpt8bbbp/rhoin8);
-          fh2JetsumHT3R10bbbp->Fill(ptcorr,sumpt10bbbp/rhoin10);}}}
-
-         }
+        
+        }
+          //end of track loop, we only do it if EM is switched off
+         
 
 
 
@@ -1484,6 +1065,17 @@ Double_t AliAnalysisTaskJetCore::RelativePhi(Double_t mphi,Double_t vphi){
   return dphi;//dphi in [-Pi, Pi]
 }
 
+Int_t AliAnalysisTaskJetCore::GetPhiBin(Double_t phi)
+{
+    Int_t phibin=-1;
+    if(!(TMath::Abs(phi)<=2*TMath::Pi())){AliError("phi w.r.t. RP out of defined range");return -1;}
+    Double_t phiwrtrp=TMath::ACos(TMath::Abs(TMath::Cos(phi)));
+    phibin=Int_t(fNRPBins*phiwrtrp/(0.5*TMath::Pi()));
+    if(phibin<0||phibin>=fNRPBins){AliError("Phi Bin not defined");}
+    return phibin;
+}
+
+
 
 
 THnSparse* AliAnalysisTaskJetCore::NewTHnSparseF(const char* name, UInt_t entries)
index 9288f21..56661f6 100644 (file)
@@ -38,7 +38,8 @@ public:
    virtual Int_t      GetNInputTracks();
      
    Double_t RelativePhi(Double_t angle1,Double_t angle2);     
-    virtual THnSparse* NewTHnSparseF(const char* name, UInt_t entries);
+   Int_t   GetPhiBin(Double_t phi);
+   virtual THnSparse* NewTHnSparseF(const char* name, UInt_t entries);
    virtual void       GetDimParams(Int_t iEntry,TString &label, Int_t &nbins, Double_t &xmin, Double_t &xmax);
    virtual AliVEvent::EOfflineTriggerTypes GetOfflineTrgMask() const { return fOfflineTrgMask; }
    virtual void     GetBranchNames(TString &branch1, TString &branch2) const { branch1 = fJetBranchName[0]; branch2 = fJetBranchName[1]; }
@@ -77,6 +78,7 @@ public:
    virtual void     SetCheckMethods(Int_t yesno){fCheckMethods=yesno;}
    virtual void     SetEventMixing(Int_t yesno){fDoEventMixing=yesno;}
    virtual void     SetFlagPhiBkg(Int_t yesno){fFlagPhiBkg=yesno;}
+   virtual void     SetFlagRandom(Int_t yesno){fFlagRandom=yesno;}
 
    virtual void     SetJetEtaMin(Float_t eta) { fJetEtaMin = eta; }
    virtual void     SetJetEtaMax(Float_t eta) { fJetEtaMax = eta; }
@@ -123,6 +125,9 @@ private:
    Int_t   fCheckMethods;     //to look into more detail into the core
    Int_t   fDoEventMixing;
    Int_t   fFlagPhiBkg;
+   Int_t   fFlagRandom;
+   Int_t   fRPAngle;
+   Int_t   fNRPBins;
    Float_t fJetEtaMin;        // lower bound on eta for found jets
    Float_t fJetEtaMax;        // upper bound on eta for found jets
    Int_t   fNevents;          // number of events
@@ -156,8 +161,9 @@ private:
    TH2F      *fh2JetCoreMethod2C30;          //Energy fraction in the core C30 method 2
    TH2F      *fh2JetCoreMethod1C60;          //Energy fraction in the core C60 method 1
    TH2F      *fh2JetCoreMethod2C60;          //Energy fraction in the core C60 method 2
-     TH3F*      fh3JetTrackC10;         //C10 pt2
-     TH3F*      fh3JetTrackC20;         //C10 pt2
+     TH3F*      fh3JetTrackC3060;           //C3060 pt2
+     TH3F*      fh3JetTrackC20;             //C10 pt2
+     TH3F*      fh3JetTrackC4080;           //C4080
      TH2F*      fh2AngStructpt1C10;         //Average 
      TH2F*      fh2AngStructpt2C10;         //C10 pt2
      TH2F*      fh2AngStructpt3C10;         //C10 pt3
@@ -175,77 +181,14 @@ private:
      TH2F*      fh2AngStructpt3C60;         //C60 pt3
      TH2F*      fh2AngStructpt4C60;         //C60 pt4
   
-     TH2F*      fh2JetsumHT3R2a;             //jet shape 02
-     TH2F*      fh2JetsumHT3R2ap;             //jet shape 02
-     TH2F*      fh2JetsumHT3R4a;             //jet shape 02
-     TH2F*      fh2JetsumHT3R4ap;             //jet shape 02
-     TH2F*      fh2JetsumHT3R6a;             //jet shape 02
-     TH2F*      fh2JetsumHT3R6ap;             //jet shape 02
-     TH2F*      fh2JetsumHT3R8a;             //jet shape 02
-     TH2F*      fh2JetsumHT3R8ap;             //jet shape 02
-     TH2F*      fh2JetsumHT3R10a;             //jet shape 12
-     TH2F*      fh2JetsumHT3R10ap;             //jet shape 12
-     TH2F*      fh2JetsumHT3R2aa;             //jet shape 02
-     TH2F*      fh2JetsumHT3R2aap;             //jet shape 02
-     TH2F*      fh2JetsumHT3R4aa;             //jet shape 02
-     TH2F*      fh2JetsumHT3R4aap;             //jet shape 02
-     TH2F*      fh2JetsumHT3R6aa;             //jet shape 02
-     TH2F*      fh2JetsumHT3R6aap;             //jet shape 02
-     TH2F*      fh2JetsumHT3R8aa;             //jet shape 02
-     TH2F*      fh2JetsumHT3R8aap;             //jet shape 02
-     TH2F*      fh2JetsumHT3R10aa;             //jet shape 12
-     TH2F*      fh2JetsumHT3R10aap;             //jet shape 12
-     TH2F*      fh2JetsumHT3R2aaa;             //jet shape 02
-     TH2F*      fh2JetsumHT3R2aaap;             //jet shape 02
-     TH2F*      fh2JetsumHT3R4aaa;             //jet shape 02
-     TH2F*      fh2JetsumHT3R4aaap;             //jet shape 02
-     TH2F*      fh2JetsumHT3R6aaa;             //jet shape 02
-     TH2F*      fh2JetsumHT3R6aaap;             //jet shape 02
-     TH2F*      fh2JetsumHT3R8aaa;             //jet shape 02
-     TH2F*      fh2JetsumHT3R8aaap;             //jet shape 02
-     TH2F*      fh2JetsumHT3R10aaa;             //jet shape 12
-     TH2F*      fh2JetsumHT3R10aaap;             //jet shape 12
-     TH2F*      fh2JetsumHT3R2b;             //jet shape 02
-     TH2F*      fh2JetsumHT3R2bp;             //jet shape 02
-     TH2F*      fh2JetsumHT3R4b;             //jet shape 02
-     TH2F*      fh2JetsumHT3R4bp;             //jet shape 02
-     TH2F*      fh2JetsumHT3R6b;             //jet shape 02
-     TH2F*      fh2JetsumHT3R6bp;             //jet shape 02
-     TH2F*      fh2JetsumHT3R8b;             //jet shape 02
-     TH2F*      fh2JetsumHT3R8bp;             //jet shape 02
-     TH2F*      fh2JetsumHT3R10b;             //jet shape 12
-     TH2F*      fh2JetsumHT3R10bp;             //jet shape 12
-     TH2F*      fh2JetsumHT3R2bb;             //jet shape 02
-     TH2F*      fh2JetsumHT3R2bbp;             //jet shape 02
-     TH2F*      fh2JetsumHT3R4bb;             //jet shape 02
-     TH2F*      fh2JetsumHT3R4bbp;             //jet shape 02
-     TH2F*      fh2JetsumHT3R6bb;             //jet shape 02
-     TH2F*      fh2JetsumHT3R6bbp;             //jet shape 02
-     TH2F*      fh2JetsumHT3R8bb;             //jet shape 02
-     TH2F*      fh2JetsumHT3R8bbp;             //jet shape 02
-     TH2F*      fh2JetsumHT3R10bb;             //jet shape 12
-     TH2F*      fh2JetsumHT3R10bbp;             //jet shape 12
-     TH2F*      fh2JetsumHT3R2bbb;             //jet shape 02
-     TH2F*      fh2JetsumHT3R2bbbp;             //jet shape 02
-     TH2F*      fh2JetsumHT3R4bbb;             //jet shape 02
-     TH2F*      fh2JetsumHT3R4bbbp;             //jet shape 02
-     TH2F*      fh2JetsumHT3R6bbb;             //jet shape 02
-     TH2F*      fh2JetsumHT3R6bbbp;             //jet shape 02
-     TH2F*      fh2JetsumHT3R8bbb;             //jet shape 02
-     TH2F*      fh2JetsumHT3R8bbbp;             //jet shape 02
-     TH2F*      fh2JetsumHT3R10bbb;             //jet shape 12
-     TH2F*      fh2JetsumHT3R10bbbp;             //jet shape 12
-    
-
      TH2F*      fh2Ntriggers;
      TH2F*      fh2JetDensity;
      TH2F*      fh2JetDensityA4;
-     TH3F*      fh3spectriggeredC10;           //triggered spectra
+     TH2F*      fh2RPJets;
+     TH3F*      fh3spectriggeredC4080;         //triggered spectra
      TH3F*      fh3spectriggeredC20;           //triggered spectra
-     TH3F*      fh3spectriggeredC3060;           //triggered spectra
-     TH3F*      fh3specbiased;              //biased spectra
-     TH3F*      fh3spectot;                 //the two combined  
-     TH3F*      fh3spectotb;                 //the two combined    
+     TH3F*      fh3spectriggeredC3060;         //triggered spectra
+
      
      Double_t            fTrigBuffer[10][7];      //!buffer for triggers