Updates for k_T algo wiht different distance paramater
authorkleinb <kleinb@f7af4fe6-9843-0410-8265-dc069ae4e863>
Sat, 30 Jan 2010 11:13:22 +0000 (11:13 +0000)
committerkleinb <kleinb@f7af4fe6-9843-0410-8265-dc069ae4e863>
Sat, 30 Jan 2010 11:13:22 +0000 (11:13 +0000)
PWG4/JetTasks/AliAnalysisTaskJetSpectrum2.cxx
PWG4/JetTasks/AliAnalysisTaskJetSpectrum2.h
PWG4/macros/AddTaskJetSpectrum2.C
PWG4/macros/AnalysisTrainPWG4Jets.C

index fc3732b..9529e09 100644 (file)
@@ -69,7 +69,9 @@ AliAnalysisTaskJetSpectrum2::AliAnalysisTaskJetSpectrum2(): AliAnalysisTaskSE(),
   fhnCorrelation(0x0),
   fBranchRec("jets"),
   fBranchGen(""),
-  fUseAODInput(kFALSE),
+  fUseAODJetInput(kFALSE),
+  fUseAODTrackInput(kFALSE),
+  fUseAODMCInput(kFALSE),
   fUseGlobalSelection(kFALSE),
   fUseExternalWeightOnly(kFALSE),
   fLimitGenJetEta(kFALSE),
@@ -91,12 +93,16 @@ AliAnalysisTaskJetSpectrum2::AliAnalysisTaskJetSpectrum2(): AliAnalysisTaskSE(),
   fh1SumPtTrackRec(0x0),  
   fh1SumPtTrackAreaRec(0x0),  
   fh1PtJetsRecIn(0x0),
+  fh1PtJetsLeadingRecIn(0x0),
   fh1PtTracksRecIn(0x0),
+  fh1PtTracksLeadingRecIn(0x0),
   fh1PtTracksGenIn(0x0),
   fh2NRecJetsPt(0x0),
   fh2NRecTracksPt(0x0),
   fh2JetsLeadingPhiEta(0x0),
+  fh2JetsLeadingPhiPt(0x0),
   fh2TracksLeadingPhiEta(0x0),
+  fh2TracksLeadingPhiPt(0x0),
   fHistList(0x0)  
 {
   for(int i = 0;i < kMaxStep*2;++i){
@@ -117,7 +123,9 @@ AliAnalysisTaskJetSpectrum2::AliAnalysisTaskJetSpectrum2(const char* name):
   fhnCorrelation(0x0),
   fBranchRec("jets"),
   fBranchGen(""),
-  fUseAODInput(kFALSE),
+  fUseAODJetInput(kFALSE),
+  fUseAODTrackInput(kFALSE),
+  fUseAODMCInput(kFALSE),
   fUseGlobalSelection(kFALSE),
   fUseExternalWeightOnly(kFALSE),
   fLimitGenJetEta(kFALSE),
@@ -139,12 +147,16 @@ AliAnalysisTaskJetSpectrum2::AliAnalysisTaskJetSpectrum2(const char* name):
   fh1SumPtTrackRec(0x0),  
   fh1SumPtTrackAreaRec(0x0),  
   fh1PtJetsRecIn(0x0),
+  fh1PtJetsLeadingRecIn(0x0),
   fh1PtTracksRecIn(0x0),
+  fh1PtTracksLeadingRecIn(0x0),
   fh1PtTracksGenIn(0x0),
   fh2NRecJetsPt(0x0),
   fh2NRecTracksPt(0x0),
   fh2JetsLeadingPhiEta(0x0),
+  fh2JetsLeadingPhiPt(0x0),
   fh2TracksLeadingPhiEta(0x0),
+  fh2TracksLeadingPhiPt(0x0),
   fHistList(0x0)
 {
 
@@ -272,17 +284,23 @@ void AliAnalysisTaskJetSpectrum2::UserCreateOutputObjects()
   fh1SumPtTrackAreaRec = new TH1F("fh1SumPtTrackAreaRec","Sum Rec track P_T #eta <0.9 / 1.8 * 2 * 0.4*0.4;p_{T,sum} (GeV/c)",nBinPt,binLimitsPt);
   
   fh1PtJetsRecIn  = new TH1F("fh1PtJetsRecIn","Rec jets P_T;p_{T} (GeV/c)",nBinPt,binLimitsPt);
+  fh1PtJetsLeadingRecIn = new TH1F("fh1PtJetsLeadingRecIn","Rec jets P_T;p_{T} (GeV/c)",nBinPt,binLimitsPt);
   fh1PtTracksRecIn  = new TH1F("fh1PtTracksRecIn","Rec tracks P_T #eta < 0.9;p_{T} (GeV/c)",nBinPt,binLimitsPt);
+  fh1PtTracksLeadingRecIn  = new TH1F("fh1PtTracksLeadingRecIn","Rec tracks P_T #eta < 0.9;p_{T} (GeV/c)",nBinPt,binLimitsPt);
   fh1PtTracksGenIn  = new TH1F("fh1PtTracksGenIn","gen tracks P_T #eta < 0.9;p_{T} (GeV/c)",nBinPt,binLimitsPt);
 
   fh2NRecJetsPt = new TH2F("fh2NRecJetsPt","Number of jets above threshhold;p_{T,cut} (GeV/c);N_{jets}",nBinPt,binLimitsPt,50,-0.5,49.5);
   fh2NRecTracksPt = new TH2F("fh2NRecTracksPt","Number of tracks above threshhold;p_{T,cut} (GeV/c);N_{tracks}",nBinPt,binLimitsPt,50,-0.5,49.5);
   // 
 
-  fh2JetsLeadingPhiEta = new TH2F("fh2JetsLeadingPhiEta","delta eta vs delta phi to leading;#Delta#phi;#Delta#eta",
-                               nBinPhi,binLimitsPhi,nBinEta,binLimitsEta);
-  fh2TracksLeadingPhiEta = new TH2F("fh2TracksLeadingPhiEta","delta eta vs delta phi to leading;#Delta#phi;#Delta#eta",
+  fh2JetsLeadingPhiEta = new TH2F("fh2JetsLeadingPhiEta","delta eta vs delta phi to leading jet;#Delta#phi;#Delta#eta",
                                nBinPhi,binLimitsPhi,nBinEta,binLimitsEta);
+  fh2JetsLeadingPhiPt = new TH2F("fh2JetsLeadingPhiPt","leading p_T vs delta phi to leading jet;#Delta#phi;p_{T} (GeV/c)",
+                               nBinPhi,binLimitsPhi,nBinPt,binLimitsPt);
+  fh2TracksLeadingPhiEta = new TH2F("fh2TracksLeadingPhiEta","delta eta vs delta phi to leading track;#Delta#phi;#Delta#eta",
+                                   nBinPhi,binLimitsPhi,nBinEta,binLimitsEta);
+  fh2TracksLeadingPhiPt = new TH2F("fh2TracksLeadingPhiPt","leading p_T vs delta phi to leading jet;#Delta#phi;p_{T} (GeV/c)",
+                                nBinPhi,binLimitsPhi,nBinPt,binLimitsPt);
 
   for(int ij = 0;ij < kMaxJets;++ij){
     fh1PtRecIn[ij] = new TH1F(Form("fh1PtRecIn_j%d",ij),"rec p_T input ;p_{T,rec}",nBinPt,binLimitsPt);
@@ -319,7 +337,9 @@ void AliAnalysisTaskJetSpectrum2::UserCreateOutputObjects()
       fHistList->Add(fh1PtTracksGenIn);
     }
     fHistList->Add(fh1PtJetsRecIn);
+    fHistList->Add(fh1PtJetsLeadingRecIn);
     fHistList->Add(fh1PtTracksRecIn);
+    fHistList->Add(fh1PtTracksLeadingRecIn);
     fHistList->Add(fh1NRecJets);
     fHistList->Add(fh1PtTrackRec);
     fHistList->Add(fh1SumPtTrackRec);
@@ -327,7 +347,9 @@ void AliAnalysisTaskJetSpectrum2::UserCreateOutputObjects()
     fHistList->Add(fh2NRecJetsPt);
     fHistList->Add(fh2NRecTracksPt);
     fHistList->Add(fh2JetsLeadingPhiEta );
+    fHistList->Add(fh2JetsLeadingPhiPt );
     fHistList->Add(fh2TracksLeadingPhiEta);
+    fHistList->Add(fh2TracksLeadingPhiPt);
     for(int i = 0;i<kMaxStep*2;++i)fHistList->Add(fhnJetContainer[i]);
     for(int ij = 0;ij<kMaxJets;++ij){
       fHistList->Add( fh1PtRecIn[ij]);
@@ -382,10 +404,10 @@ void AliAnalysisTaskJetSpectrum2::UserExec(Option_t */*option*/)
   // Execute analysis for current event
   //
   AliESDEvent *fESD = 0;
-  if(fUseAODInput){    
+  if(fUseAODJetInput){    
     fAOD = dynamic_cast<AliAODEvent*>(InputEvent());
     if(!fAOD){
-      Printf("%s:%d AODEvent not found in Input Manager %d",(char*)__FILE__,__LINE__,fUseAODInput);
+      Printf("%s:%d AODEvent not found in Input Manager %d",(char*)__FILE__,__LINE__,fUseAODJetInput);
       return;
     }
     // fethc the header
@@ -499,8 +521,13 @@ void AliAnalysisTaskJetSpectrum2::UserExec(Option_t */*option*/)
   TList recParticles;
   TList genParticles;
 
-  GetListOfTracks(&recParticles,fTrackTypeRec);
-  GetListOfTracks(&genParticles,fTrackTypeGen);
+
+
+
+  Int_t nT = GetListOfTracks(&recParticles,fTrackTypeRec);
+  if(fDebug>2)Printf("%s:%d Selected Rec tracks: %d %d",(char*)__FILE__,__LINE__,nT,recParticles.GetEntries());
+  nT = GetListOfTracks(&genParticles,fTrackTypeGen);
+  if(fDebug>2)Printf("%s:%d Selected Gen tracks: %d %d",(char*)__FILE__,__LINE__,nT,genParticles.GetEntries());
 
 
   if (fDebug > 10)Printf("%s:%d",(char*)__FILE__,__LINE__);
@@ -581,10 +608,13 @@ void AliAnalysisTaskJetSpectrum2::UserExec(Option_t */*option*/)
     Float_t phi = leading->Phi();
     if(phi<0)phi+=TMath::Pi()*2.;    
     Float_t eta = leading->Eta();
-    while((tmpRec = (AliAODJet*)(recIter->Next()))){
+     while((tmpRec = (AliAODJet*)(recIter->Next()))){
       Float_t tmpPt = tmpRec->Pt();
       fh1PtJetsRecIn->Fill(tmpPt);
-      if(tmpRec==leading)continue;
+      if(tmpRec==leading){
+       fh1PtJetsLeadingRecIn->Fill(tmpPt);
+       continue;
+      }
       // correlation
       Float_t tmpPhi =  tmpRec->Phi();
 
@@ -594,6 +624,7 @@ void AliAnalysisTaskJetSpectrum2::UserExec(Option_t */*option*/)
       if(dPhi<(-1.*TMath::Pi()))dPhi = dPhi + 2.*TMath::Pi();      
       Float_t dEta = eta - tmpRec->Eta();
       fh2JetsLeadingPhiEta->Fill(dPhi,dEta);
+      fh2JetsLeadingPhiPt->Fill(dPhi,pt);
     }  
     delete recIter;
   }
@@ -626,7 +657,10 @@ void AliAnalysisTaskJetSpectrum2::UserExec(Option_t */*option*/)
     while((tmpRec = (AliVParticle*)(recIter->Next()))){
       Float_t tmpPt = tmpRec->Pt();
       fh1PtTracksRecIn->Fill(tmpPt);
-      if(tmpRec==leading)continue;
+      if(tmpRec==leading){
+       fh1PtTracksLeadingRecIn->Fill(tmpPt);
+       continue;
+      }
       // correlation
       Float_t tmpPhi =  tmpRec->Phi();
 
@@ -636,6 +670,7 @@ void AliAnalysisTaskJetSpectrum2::UserExec(Option_t */*option*/)
       if(dPhi<(-1.*TMath::Pi()))dPhi = dPhi + 2.*TMath::Pi();      
       Float_t dEta = eta - tmpRec->Eta();
       fh2TracksLeadingPhiEta->Fill(dPhi,dEta);
+      fh2TracksLeadingPhiPt->Fill(dPhi,pt);
     }  
     delete recIter;
   }
@@ -916,12 +951,13 @@ void AliAnalysisTaskJetSpectrum2::Terminate(Option_t */*option*/)
 
 Int_t  AliAnalysisTaskJetSpectrum2::GetListOfTracks(TList *list,Int_t type){
 
+  if(fDebug>2)Printf("%s:%d Selecting tracks with %d",(char*)__FILE__,__LINE__,type);
 
   Int_t iCount = 0;
-  if(type==kTrackAODIn||type==kTrackAODOut){
+  if(type==kTrackAOD){
     AliAODEvent *aod = 0;
-    if(type==kTrackAODIn)aod = dynamic_cast<AliAODEvent*>(InputEvent());
-    else if (type==kTrackAODOut) aod = AODEvent();
+    if(fUseAODTrackInput)aod = dynamic_cast<AliAODEvent*>(InputEvent());
+    else aod = AODEvent();
     if(!aod){
       return iCount;
     }
@@ -951,9 +987,10 @@ Int_t  AliAnalysisTaskJetSpectrum2::GetListOfTracks(TList *list,Int_t type){
       }
     }
   }
-  else if (type == kTrackAODMCCharged || type == kTrackAODMCAll ) {
-    AliAODEvent *aod = dynamic_cast<AliAODEvent*>(InputEvent());
-    if(!aod) aod = AODEvent();
+  else if (type == kTrackAODMCCharged || type == kTrackAODMCAll || type == kTrackAODMCChargedAcceptance) {
+    AliAODEvent *aod = 0;
+    if(fUseAODMCInput)dynamic_cast<AliAODEvent*>(InputEvent());
+    else aod = AODEvent();
     if(!aod)return iCount;
     TClonesArray *tca = dynamic_cast<TClonesArray*>(aod->FindListObject(AliAODMCParticle::StdBranchName()));
     if(!tca)return iCount;
@@ -964,9 +1001,15 @@ Int_t  AliAnalysisTaskJetSpectrum2::GetListOfTracks(TList *list,Int_t type){
        list->Add(part);
        iCount++;
       }
-      else if (type == kTrackAODMCCharged){
+      else if (type == kTrackAODMCCharged || type == kTrackAODMCChargedAcceptance ){
        if(part->Charge()==0)continue;
-       list->Add(part);
+       if(kTrackAODMCCharged){
+         list->Add(part);
+       }
+       else {
+         if(TMath::Abs(part->Eta())>0.9)continue;
+         list->Add(part);
+       }
        iCount++;
       }
     }
index 1501fa5..ace7cfa 100644 (file)
@@ -45,7 +45,9 @@ class AliAnalysisTaskJetSpectrum2 : public AliAnalysisTaskSE
     virtual void SetUseGlobalSelection(Bool_t b){fUseGlobalSelection = b;}
     virtual void SetExternalWeight(Float_t f){fExternalWeight = f;}
     virtual void SetUseExternalWeightOnly(Bool_t b){fUseExternalWeightOnly = b;}
-    virtual void SetAODInput(Bool_t b){fUseAODInput = b;}
+    virtual void SetAODJetInput(Bool_t b){fUseAODJetInput = b;}
+    virtual void SetAODTrackInput(Bool_t b){fUseAODTrackInput = b;}
+    virtual void SetAODMCInput(Bool_t b){fUseAODMCInput = b;}
     virtual void SetLimitGenJetEta(Bool_t b){fLimitGenJetEta = b;}
     virtual void SetRecEtaWindow(Float_t f){fRecEtaWindow = f;}
     virtual void SetAnalysisType(Int_t i){fAnalysisType = i;}
@@ -64,7 +66,7 @@ class AliAnalysisTaskJetSpectrum2 : public AliAnalysisTaskSE
     // AOD reading -> MC from AOD
     // ESD reading -> MC from Kinematics
     // this has to match with our selection of input events
-    enum {kTrackUndef = 0, kTrackAODIn, kTrackAODOut, kTrackKineAll,kTrackKineCharged, kTrackAODMCAll, kTrackAODMCCharged };
+    enum {kTrackUndef = 0, kTrackAOD, kTrackKineAll,kTrackKineCharged, kTrackAODMCAll, kTrackAODMCCharged, kTrackAODMCChargedAcceptance};
     enum {kAnaMC =  0x1, kAnaMCESD = 0x2};
     enum {kMaxJets = 4};
     enum {kMaxCorrelation =  3};
@@ -96,7 +98,9 @@ class AliAnalysisTaskJetSpectrum2 : public AliAnalysisTaskSE
     TString       fBranchRec;  // AOD branch name for reconstructed
     TString       fBranchGen;  // AOD brnach for genereated
 
-    Bool_t        fUseAODInput;           // take jet from input AOD not from ouptu AOD
+    Bool_t        fUseAODJetInput;        // take jet from input AOD not from ouptu AOD
+    Bool_t        fUseAODTrackInput;      // take track from input AOD not from ouptu AOD
+    Bool_t        fUseAODMCInput;         // take MC from input AOD not from ouptu AOD
     Bool_t        fUseGlobalSelection;    // Limit the eta of the generated jets
     Bool_t        fUseExternalWeightOnly; // use only external weight
     Bool_t        fLimitGenJetEta;        // Limit the eta of the generated jets
@@ -124,14 +128,18 @@ class AliAnalysisTaskJetSpectrum2 : public AliAnalysisTaskSE
     TH1F*         fh1PtGenIn[kMaxJets];  // Jet pt with corellated generated jet    
 
     TH1F*         fh1PtJetsRecIn;  // Jet pt for all jets
+    TH1F*         fh1PtJetsLeadingRecIn;  // Jet pt for all jets
     TH1F*         fh1PtTracksRecIn;  // track pt for all tracks
+    TH1F*         fh1PtTracksLeadingRecIn;  // track pt for all tracks
     TH1F*         fh1PtTracksGenIn;  // track pt for all tracks
 
 
     TH2F*         fh2NRecJetsPt;            // Number of found jets above threshold
     TH2F*         fh2NRecTracksPt;          // Number of found tracks above threshold
-    TH2F*         fh2JetsLeadingPhiEta;     // jet corrlation with leading jet
-    TH2F*         fh2TracksLeadingPhiEta;   // jet corrlation with leading jet
+    TH2F*         fh2JetsLeadingPhiEta;     // jet correlation with leading jet
+    TH2F*         fh2JetsLeadingPhiPt;      // jet correlation with leading jet
+    TH2F*         fh2TracksLeadingPhiEta;   // track correlation with leading track
+    TH2F*         fh2TracksLeadingPhiPt;    // track correlation with leading track
 
     TH2F*         fh2PhiPt[kMaxJets];    // delta phi correlation of tracks with the jet      
     TH2F*         fh2PhiEta[kMaxJets];   // eta   phi correlation of tracks with the jet      
index 55c11fa..0c44110 100644 (file)
@@ -10,14 +10,18 @@ AliAnalysisTaskJetSpectrum2 *AddTaskJetSpectrum2Delta(UInt_t filterMask = 16,Boo
     if(iFlag&(1<<2))js = AddTaskJetSpectrum2("jetsAOD_FASTJET04","jetsAODMC_FASTJET04",filterMask,iPhysicsSelection);\r
     if(iFlag&(1<<3))js = AddTaskJetSpectrum2("jetsAOD_FASTJET04","jetsAODMC2_FASTJET04",filterMask,iPhysicsSelection);\r
 \r
-    if(iFlag&(1<<4))js = AddTaskJetSpectrum2("jetsAOD_FASTKT04","jetsAODMC_FASTKT04",filterMask,iPhysicsSelection);\r
+    if(iFlag&(1<<4)){\r
+      js = AddTaskJetSpectrum2("jetsAOD_FASTKT04","jetsAODMC_FASTKT04",filterMask,iPhysicsSelection);\r
+    }\r
     if(iFlag&(1<<5))js = AddTaskJetSpectrum2("jetsAOD_FASTKT04","jetsAODMC2_FASTKT04",filterMask,iPhysicsSelection);\r
-\r
     if(iFlag&(1<<6))js = AddTaskJetSpectrum2("jetsAOD_UA107","jetsAODMC_UA107",filterMask,iPhysicsSelection);\r
   }\r
+\r
   if(iFlag&(1<<7))js = AddTaskJetSpectrum2("jets","jetsAOD_FASTJET04",filterMask,iPhysicsSelection);\r
   if(iFlag&(1<<8))js = AddTaskJetSpectrum2("jetsAOD_FASTJET04","",filterMask,iPhysicsSelection);\r
   if(iFlag&(1<<9))js = AddTaskJetSpectrum2("jetsAOD_FASTKT04","",filterMask,iPhysicsSelection);\r
+  \r
+\r
   if(iFlag&(1<<10)){\r
     js = AddTaskJetSpectrum2("jetsAOD_UA107","",filterMask,iPhysicsSelection);\r
     js->SetRecEtaWindow(0.2);\r
@@ -46,8 +50,10 @@ AliAnalysisTaskJetSpectrum2 *AddTaskJetSpectrum2(char* bRec,char* bGen ,UInt_t f
    }\r
 \r
    TString type = mgr->GetInputEventHandler()->GetDataType();\r
-   TString typeMC(bGen);\r
-   typeMC.ToUpper();\r
+   TString typeRec(bRec);\r
+   TString typeGen(bGen);\r
+   typeGen.ToUpper();\r
+   typeRec.ToUpper();\r
    // Create the task and configure it.\r
    //===========================================================================\r
    \r
@@ -63,25 +69,42 @@ AliAnalysisTaskJetSpectrum2 *AddTaskJetSpectrum2(char* bRec,char* bGen ,UInt_t f
    //   pwg4spec->SetUseGlobalSelection(kTRUE); \r
 \r
    if(type == "AOD"){\r
-     pwg4spec->SetTrackTypeRec(AliAnalysisTaskJetSpectrum2::kTrackAODIn);\r
-     pwg4spec->SetAODInput(kTRUE);\r
+     // Assume all jet are produced already\r
+     pwg4spec->SetAODJetInput(kTRUE);\r
+     pwg4spec->SetAODTrackInput(kTRUE);\r
+     pwg4spec->SetAODMCInput(kTRUE);\r
+   }\r
+\r
+   if(typeRec.Contains("AODMC2b")){// work down from the top AODMC2b -> AODMC2 -> AODMC -> AOD\r
+     pwg4spec->SetTrackTypeRec(AliAnalysisTaskJetSpectrum2::kTrackAODMCChargedAcceptance);\r
+   }\r
+   else if (typeRec.Contains("AODMC2")){\r
+     pwg4spec->SetTrackTypeRec(AliAnalysisTaskJetSpectrum2::kTrackAODMCCharged);\r
+   }\r
+   else if (typeRec.Contains("AODMC")){\r
+     pwg4spec->SetTrackTypeRec(AliAnalysisTaskJetSpectrum2::kTrackAODMCall);\r
+   }\r
+   else if (typeRec.Contains("AOD")) {\r
+     pwg4spec->SetTrackTypeRec(AliAnalysisTaskJetSpectrum2::kTrackAOD);\r
    }\r
-   else pwg4spec->SetTrackTypeRec(AliAnalysisTaskJetSpectrum2::kTrackAODOut);\r
 \r
-   if(typeMC.Contains("AODMC2")){\r
+   if(typeGen.Contains("AODMC2b")){// work down from the top AODMC2b -> AODMC2 -> AODMC -> AOD\r
+     pwg4spec->SetTrackTypeGen(AliAnalysisTaskJetSpectrum2::kTrackAODMCChargedAcceptance);\r
+   }\r
+   else if (typeGen.Contains("AODMC2")){\r
      pwg4spec->SetTrackTypeGen(AliAnalysisTaskJetSpectrum2::kTrackAODMCCharged);\r
    }\r
-   else if (typeMC.Contains("AODMC2")){\r
-     pwg4spec->SetTrackTypeGen(AliAnalysisTaskJetSpectrum2::kTrackAODMCAll);\r
+   else if (typeGen.Contains("AODMC")){\r
+     pwg4spec->SetTrackTypeGen(AliAnalysisTaskJetSpectrum2::kTrackAODMCall);\r
    }\r
-   else if (typeMC.Contains("AOD")) {\r
-     if(type == "AOD"){\r
-       pwg4spec->SetTrackTypeGen(AliAnalysisTaskJetSpectrum2::kTrackAODIn);\r
-     }\r
-     else{\r
-       pwg4spec->SetTrackTypeGen(AliAnalysisTaskJetSpectrum2::kTrackAODOut);\r
-     }\r
+   else if (typeGen.Contains("AOD")) {\r
+     pwg4spec->SetTrackTypeGen(AliAnalysisTaskJetSpectrum2::kTrackAOD);\r
    }\r
+\r
+\r
+\r
+\r
+\r
    if(iPhysicsSelection)pwg4spec->SelectCollisionCandidates();\r
 \r
    mgr->AddTask(pwg4spec);\r
index 82c0f32..b1625a3 100644 (file)
@@ -366,10 +366,19 @@ void AnalysisTrainPWG4Jets(const char *analysis_mode="local",
        AddTaskJetSpectrum2Delta(kHighPtFilterMask,kUseAODMC,iPhysicsSelection,selection);
      }
      if(iJETAN&3&&!kFillAOD&!iAODanalysis){
+       AddTaskJetSpectrum2("jetsAOD_FASTKT00","",kHighPtFilterMask,iPhysicsSelection);
        AddTaskJetSpectrum2("jetsAOD_FASTKT01","",kHighPtFilterMask,iPhysicsSelection);
        AddTaskJetSpectrum2("jetsAOD_FASTKT02","",kHighPtFilterMask,iPhysicsSelection);
        AddTaskJetSpectrum2("jetsAOD_FASTKT06","",kHighPtFilterMask,iPhysicsSelection);
        AddTaskJetSpectrum2("jetsAOD_FASTKT08","",kHighPtFilterMask,iPhysicsSelection);
+       if(kIsMC){
+        AddTaskJetSpectrum2("jetsAODMC2b_FASTKT00","",kHighPtFilterMask,iPhysicsSelection);
+        AddTaskJetSpectrum2("jetsAODMC2b_FASTKT00","",kHighPtFilterMask,iPhysicsSelection);
+        AddTaskJetSpectrum2("jetsAODMC2b_FASTKT01","",kHighPtFilterMask,iPhysicsSelection);
+        AddTaskJetSpectrum2("jetsAODMC2b_FASTKT02","",kHighPtFilterMask,iPhysicsSelection);
+        AddTaskJetSpectrum2("jetsAODMC2b_FASTKT06","",kHighPtFilterMask,iPhysicsSelection);
+        AddTaskJetSpectrum2("jetsAODMC2b_FASTKT08","",kHighPtFilterMask,iPhysicsSelection);
+       }
      }
      if(kUseMC)taskjetSpectrum->SetAnalysisType( AliAnalysisTaskJetSpectrum2::kAnaMCESD);
      taskjetSpectrum->SetDebugLevel(0);