reduced number of track cuts for QA task, changed binning for centrality
authorkleinb <kleinb@f7af4fe6-9843-0410-8265-dc069ae4e863>
Thu, 7 Jul 2011 12:34:52 +0000 (12:34 +0000)
committerkleinb <kleinb@f7af4fe6-9843-0410-8265-dc069ae4e863>
Thu, 7 Jul 2011 12:34:52 +0000 (12:34 +0000)
PWG4/JetTasks/AliAnalysisTaskJetServices.cxx
PWG4/JetTasks/AliAnalysisTaskJetSpectrum2.cxx
PWG4/macros/AddTaskPWG4HighPtTrackQA.C
PWG4/macros/AnalysisTrainPWG4Jets.C
PWG4/macros/ConfigTrainGrid.C

index cd4f2ff..f68feb4 100644 (file)
@@ -262,10 +262,10 @@ void AliAnalysisTaskJetServices::UserCreateOutputObjects()
   }
   
 
-  fh1CentralityESD = new TH1F("fh1CentralityESD","cent",102,-0.5,101.5);
+  fh1CentralityESD = new TH1F("fh1CentralityESD","cent",103,-1,102);
   fHistList->Add(fh1CentralityESD);
   
-  fh1Centrality = new TH1F("fh1Centrality","cent",102,-0.5,101.5);
+  fh1Centrality = new TH1F("fh1Centrality","cent",103,-1,102);
   fHistList->Add(fh1Centrality);
 
   fh1RP = new TH1F("fh1RP","RP;#Psi",440, -1.*TMath::Pi(), 2.*TMath::Pi());
index 85cfb11..04675f0 100644 (file)
@@ -414,7 +414,7 @@ void AliAnalysisTaskJetSpectrum2::UserCreateOutputObjects()
   fh1RP = new TH1F("fh1RP","RP;#Psi",440, -1.*TMath::Pi(), 2.*TMath::Pi());
   fHistList->Add(fh1RP);
 
-  fh1Centrality = new TH1F("fh1Centrality","cent;cent (%)",101,-0.5,100.5);
+  fh1Centrality = new TH1F("fh1Centrality","cent;cent (%)",103,-1,102);
   fHistList->Add(fh1Centrality);
 
   fh2MultRec = new TH2F("fh2MultRec","multiplicity rec;# tracks;# jetrefs",500,0,5000,500,0.,5000);
index 6c8533b..5ea70d6 100644 (file)
@@ -35,6 +35,21 @@ void AddTaskPWG4HighPtTrackQAAll(char *prodType = "LHC10h",Bool_t isPbPb=kTRUE,
 
 }
 
+void AddTaskPWG4HighPtTrackQAAllReduced(char *prodType = "LHC10h",Bool_t isPbPb=kTRUE, Int_t iAODanalysis = 0) 
+{    
+
+  int cent = 10;
+  
+  if(isPbPb) {
+    for(cent=0; cent<4; cent++) {
+      AliPWG4HighPtTrackQA *taskTrackQA00 = AddTaskPWG4HighPtTrackQA(prodType,isPbPb,iAODanalysis,cent,0,0);
+      AliPWG4HighPtTrackQA *taskTrackQA01 = AddTaskPWG4HighPtTrackQA(prodType,isPbPb,iAODanalysis,cent,0,1);
+      AliPWG4HighPtTrackQA *taskTrackQA02 = AddTaskPWG4HighPtTrackQA(prodType,isPbPb,iAODanalysis,cent,0,2);
+    }
+  }
+
+}
+
 void AddTaskPWG4HighPtTrackQAAOD(char *prodType = "LHC10h",Bool_t isPbPb=kTRUE, Int_t iAODanalysis = 1, Int_t filterBit) 
 {   
   AliPWG4HighPtTrackQA *taskTrackQA = AddTaskPWG4HighPtTrackQA(prodType,isPbPb,iAODanalysis,0,0,0);
index 61c1032..a81476f 100644 (file)
@@ -1203,7 +1203,8 @@ void AnalysisTrainPWG4Jets(const char *analysis_mode="local",
 
    if(iPWG4PtTrackQA){
      gROOT->LoadMacro("$ALICE_ROOT/PWG4/macros/AddTaskPWG4HighPtTrackQA.C");
-     AddTaskPWG4HighPtTrackQAAll(kGridDataSet.Data(),kIsPbPb,iAODanalysis);
+     if(iPWG4PtTrackQA&2)AddTaskPWG4HighPtTrackQAAll(kGridDataSet.Data(),kIsPbPb,iAODanalysis);
+     else AddTaskPWG4HighPtTrackQAAllReduced(kGridDataSet.Data(),kIsPbPb,iAODanalysis);
    }
 
    if(iPWG4PtQATPC){
index 8eff021..918f552 100644 (file)
@@ -39,7 +39,7 @@
   // bextra == 0 4 plus
   // bextra == 1 large pass1 split..
   // bextra == 2 3 plus
-  Int_t bRun = 802; Int_t bExtra = 0;  char* cDate = "110704a";
+  Int_t bRun = 802; Int_t bExtra = 0;  char* cDate = "110706a";
   iAODanalysis = 0; 
   // 1 == Read Jets and tracks form the input AOD
   // needs the jet branchnames set explicitly
       
       iPWG4PtQAMC     = 0;  // 
       iPWG4PtQATPC   = 0;  // 
-      iPWG4PtSpectra   = 1;  //  
+      iPWG4PtSpectra   = 0;  //  
       iPWG4PtTrackQA    = 1;
       iPWG4JetSpectrum = 1; 
       iPWG4JetServices  = 1; // !!!!!!!!!!!