missing initialisation

diff --git a/PWGJE/AliAnalysisTaskJetSpectrum2.cxx b/PWGJE/AliAnalysisTaskJetSpectrum2.cxx
index 7d671324d206c7d4a6b7211c1f521ff6014c6b65..560441cacc73b6f91a6ac4f930b7f3edf4f37164 100644 (file)
--- a/PWGJE/AliAnalysisTaskJetSpectrum2.cxx
+++ b/PWGJE/AliAnalysisTaskJetSpectrum2.cxx
@@ -112,6 +112,10 @@ AliAnalysisTaskJetSpectrum2::AliAnalysisTaskJetSpectrum2():
   fMinJetPt(0),
   fMinTrackPt(0.15),
   fDeltaPhiWindow(90./180.*TMath::Pi()),
+  fAcceptancePhiMin(0x0),
+  fAcceptancePhiMax(0x0),
+  fAcceptanceEtaMin(0x0),
+  fAcceptanceEtaMax(0x0),
   fCentrality(100),
   fRPAngle(0),
   fMultRec(0),
@@ -212,6 +216,10 @@ AliAnalysisTaskJetSpectrum2::AliAnalysisTaskJetSpectrum2(const char* name):
   fMinJetPt(0),
   fMinTrackPt(0.15),
   fDeltaPhiWindow(90./180.*TMath::Pi()),
+  fAcceptancePhiMin(0x0),
+  fAcceptancePhiMax(0x0),
+  fAcceptanceEtaMin(0x0),
+  fAcceptanceEtaMax(0x0),
   fCentrality(100),
   fRPAngle(0),
   fMultRec(0),
@@ -1017,7 +1025,7 @@ void AliAnalysisTaskJetSpectrum2::FillJetHistos(TList &jetsList,TList &particles
       var2[3] = etaJet;
       var2[4] = phiJet;
       var2[5] = jet->EffectiveAreaCharged();
-      var2[7] = var2[8];
+      var2[7] = var1[8];
       if(ij<kMaxJets){
        if(leadTrack)fh2LTrackPtJetPt[iType][ij]->Fill(leadTrack->Pt(),ptJet);
        var1[0] = ij;

diff --git a/PWGJE/macros/AddTaskJetSpectrum2.C b/PWGJE/macros/AddTaskJetSpectrum2.C
index a39c9264888a3e173a6f7c4e2d040d0f5505f0a1..7cc469235a99003f4e33e4ab40c991960b9b2c81 100644 (file)
--- a/PWGJE/macros/AddTaskJetSpectrum2.C
+++ b/PWGJE/macros/AddTaskJetSpectrum2.C
@@ -181,6 +181,23 @@ Bool_t ConfigWithFlagsJetSpectrum2(){
      Printf("%s%d: kNTrigger not defined",(char*)__FILE__,__LINE__); return kFALSE; \r
    }\r
 \r
+\r
+  Int_t nAcceptance = AliAnalysisManager::GetGlobalInt("kNAcceptanceSpec",bGood1);\r
+  \r
+  if(bGood1){\r
+    jetspec->SetNAcceptance(nAcceptance);\r
+    for(int ia = 0;ia < nAcceptance;ia++){\r
+      jetspec->SetAcceptance(ia,\r
+                            AliAnalysisManager::GetGlobalDbl(Form("kAcceptancePhiMinSpec%d",ia),bGood1),\r
+                            AliAnalysisManager::GetGlobalDbl(Form("kAcceptancePhiMaxSpec0%d",ia),bGood2),\r
+                            AliAnalysisManager::GetGlobalDbl(Form("kAcceptanceEtaMinSpec%d",ia),bGood1),\r
+                            AliAnalysisManager::GetGlobalDbl(Form("kAcceptanceEtaMaxSpec0%d",ia),bGood2));\r
+    }\r
+   }\r
+   else {\r
+     Printf("%s%d: kNAcceptance not defined",(char*)__FILE__,__LINE__); \r
+   }\r
+\r
      \r
    AliAnalysisManager::GetGlobalInt("kPhysicsSelectionFlag",bGood1);if(bGood1){\r
      jetspec->SelectCollisionCandidates(AliAnalysisManager::GetGlobalInt("kPhysicsSelectionFlag",bGood1));\r

diff --git a/PWGJE/macros/AnalysisTrainPWGJets.C b/PWGJE/macros/AnalysisTrainPWGJets.C
index 6190c8ae9b296c70a586143dcbeb257a6ecc3220..935dbbec8a73f69e6256a910e4ce5cd4e766e1fd 100644 (file)
--- a/PWGJE/macros/AnalysisTrainPWGJets.C
+++ b/PWGJE/macros/AnalysisTrainPWGJets.C
@@ -339,7 +339,7 @@ void AnalysisTrainPWGJets(const char *analysis_mode="local",
    }   
 
    gROOT->LoadMacro("$ALICE_ROOT/PWGJE/macros/ConfigLegoTrainPWGJE.C");
-   ConfigLegoTrainPWGJE(1013);
+   ConfigLegoTrainPWGJE(1108);
    // Create input handler (input container created automatically)
    if (iAODanalysis) {
    // AOD input handler

diff --git a/PWGJE/macros/ConfigLegoTrainPWGJE.C b/PWGJE/macros/ConfigLegoTrainPWGJE.C
index 3cfc32a9fc86bd089ec310cb616ee87d6a15c2e4..ba3a9aec066834ffd4d494616ec998dd701a0aba 100644 (file)
--- a/PWGJE/macros/ConfigLegoTrainPWGJE.C
+++ b/PWGJE/macros/ConfigLegoTrainPWGJE.C
@@ -50,6 +50,27 @@ void ConfigLegoTrainPWGJE(int iFlag = 0){
     AliAnalysisManager::SetGlobalStr("kTriggerName3","kSemiCentral");    
     AliAnalysisManager::SetGlobalStr("kTriggerName4","kEMCEJE");    
 
+
+    // aceptance windows for Spectrum task
+
+    Float_t fid = 0.2;
+    Float_t fEMCPhiLo = 2.17 - 1.87/2. + fid; // some fiducial region
+    Float_t fEMCPhiUp = 2.17 + 1.87/2. - fid; // some fiducial region
+    Float_t fEMCEtaLo = - 0.7 + fid;
+    Float_t fEMCEtaUp =   0.7 - fid;
+    AliAnalysisManager::SetGlobalInt("kNAcceptanceSpec",2);    
+    AliAnalysisManager::SetGlobalDbl("kAcceptancePhiMinSpec0",fEMCPhiLo);
+    AliAnalysisManager::SetGlobalDbl("kAcceptancePhiMaxSpec0",fEMCPhiUp);
+    AliAnalysisManager::SetGlobalDbl("kAcceptanceEtaMinSpec0",fEMCEtaLo);
+    AliAnalysisManager::SetGlobalDbl("kAcceptanceEtaMaxSpec0",fEMCEtaUp);
+
+    // iroc 13 
+    AliAnalysisManager::SetGlobalDbl("kAcceptancePhiMinSpec1",4.7-0.4);
+    AliAnalysisManager::SetGlobalDbl("kAcceptancePhiMaxSpec1",4.7+0.4);
+    AliAnalysisManager::SetGlobalDbl("kAcceptanceEtaMinSpec1",-0.9);
+    AliAnalysisManager::SetGlobalDbl("kAcceptanceEtaMaxSpec1",0);
+
+
     // jet selection
     AliAnalysisManager::SetGlobalDbl("kJetEtaWindow",0.5);