missing initialisation
authorkleinb <kleinb@f7af4fe6-9843-0410-8265-dc069ae4e863>
Wed, 16 May 2012 09:52:29 +0000 (09:52 +0000)
committerkleinb <kleinb@f7af4fe6-9843-0410-8265-dc069ae4e863>
Wed, 16 May 2012 09:52:29 +0000 (09:52 +0000)
PWGJE/AliAnalysisTaskJetSpectrum2.cxx
PWGJE/macros/AddTaskJetSpectrum2.C
PWGJE/macros/AnalysisTrainPWGJets.C
PWGJE/macros/ConfigLegoTrainPWGJE.C

index 7d67132..560441c 100644 (file)
@@ -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;
index a39c926..7cc4692 100644 (file)
@@ -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
index 6190c8a..935dbbe 100644 (file)
@@ -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
index 3cfc32a..ba3a9ae 100644 (file)
@@ -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);