Switched on vertex constrained again for Jet SPectrum, limit some printout in three...
authorkleinb <kleinb@f7af4fe6-9843-0410-8265-dc069ae4e863>
Mon, 11 Jan 2010 17:25:24 +0000 (17:25 +0000)
committerkleinb <kleinb@f7af4fe6-9843-0410-8265-dc069ae4e863>
Mon, 11 Jan 2010 17:25:24 +0000 (17:25 +0000)
PWG4/JetTasks/AliAnalysisTaskJetServices.cxx
PWG4/JetTasks/AliAnalysisTaskJetSpectrum2.cxx
PWG4/JetTasks/AliAnalysisTaskThreeJets.cxx

index 83d0bac..83dc75c 100644 (file)
@@ -351,10 +351,15 @@ void AliAnalysisTaskJetServices::UserExec(Option_t */*option*/)
        Float_t zvtx = vtxESD->GetZ();
        Float_t yvtx = vtxESD->GetY();
        Float_t xvtx = vtxESD->GetX();
+       Float_t r2   = yvtx*yvtx+xvtx*xvtx;  
        fh2ESDTriggerVtx->Fill(it,zvtx);
-       if(TMath::Abs(zvtx)<fZVtxCut&&esdTrig&&TMath::Abs(xvtx)<0.5&&TMath::Abs(yvtx)<0.5){
+       if(TMath::Abs(zvtx)<fZVtxCut&&esdTrig&&r2<1){
          fh2ESDTriggerCount->Fill(it,kTriggeredVertexIn);
-         if(cand)fh2ESDTriggerCount->Fill(it,kSelectedALICEVertexIn);
+         if(cand){
+           fh2ESDTriggerCount->Fill(it,kSelectedALICEVertexIn);
+           fh2ESDTriggerCount->Fill(it,kSelected);
+           AliAnalysisHelperJetTasks::Selected(kTRUE,kTRUE);// select this event
+         }
          // here we select based on ESD info...
        }
       }
index ff52940..a5f53f4 100644 (file)
@@ -613,6 +613,7 @@ void AliAnalysisTaskJetSpectrum2::UserExec(Option_t */*option*/)
     if(ptRec>10.&&fDebug>0){
       // need to cast to int, otherwise the printf overwrites
       Printf("Jet found in Event %d with p_T, %E",(int)Entry(),ptRec);
+      Printf("%s read event, %d",fInputHandler->GetTree()->GetCurrentFile()->GetName(),fInputHandler->GetTree()->GetReadEntry());
       fAOD->GetHeader()->Print();
       Printf("TriggerClasses: %s",fAOD->GetFiredTriggerClasses().Data());
       for(int it = 0;it < fAOD->GetNumberOfTracks();++it){
@@ -701,7 +702,7 @@ void AliAnalysisTaskJetSpectrum2::MakeJetContainer(){
   // link it
   //
   const Int_t kNvar   = 3 ; //number of variables on the grid:pt,eta, phi
-  const Double_t kPtmin = 10.0, kPtmax = 260.; // we do not want to have empty bins at the beginning...
+  const Double_t kPtmin = 5.0, kPtmax = 105.; // we do not want to have empty bins at the beginning...
   const Double_t kEtamin = -3.0, kEtamax = 3.0;
   const Double_t kPhimin = 0., kPhimax = 2. * TMath::Pi();
 
@@ -724,7 +725,7 @@ void AliAnalysisTaskJetSpectrum2::MakeJetContainer(){
 
   //values for bin lower bounds
   //  for(Int_t i=0; i<=iBin[0]; i++) binEdges[0][i]=(Double_t)TMath::Power(10,TMath::Log10(kPtmin) + (TMath::Log10(kPtmax)-TMath::Log10(kPtmin))/iBin[0]*(Double_t)i);  
-  for(Int_t i=0; i<=iBin[0]; i++) binEdges[0][i]=(Double_t)kPtmin  + (kPtmax-kPtmin)/iBin[0]*(Double_t)i;
+  for(Int_t i=0; i<=iBin[0]; i++) binEdges[0][i]=(Double_t)kPtmin  + (kPtmax-kPtmin)/(Double_t)iBin[0]*(Double_t)i;
   for(Int_t i=0; i<=iBin[1]; i++) binEdges[1][i]=(Double_t)kEtamin  + (kEtamax-kEtamin)/iBin[1]*(Double_t)i;
   for(Int_t i=0; i<=iBin[2]; i++) binEdges[2][i]=(Double_t)kPhimin  + (kPhimax-kPhimin)/iBin[2]*(Double_t)i;
 
index 46391d5..ca6ec55 100644 (file)
@@ -510,7 +510,7 @@ void AliAnalysisTaskThreeJets::UserExec(Option_t * )
   
   // reconstructed jets
   nRecJets = aodRecJets->GetEntries(); 
-  Printf("--- Jets found in bRec: %d", nRecJets);
+  if(fDebug)Printf("--- Jets found in bRec: %d", nRecJets);
   nRecJets = TMath::Min(nRecJets, kMaxJets);
   
   for(int ir = 0;ir < nRecJets;++ir)
@@ -915,7 +915,7 @@ void AliAnalysisTaskThreeJets::UserExec(Option_t * )
   Int_t tag1 = 0;
 
   AliAODJet recSelJets[kMaxJets];
-  Printf("---- Number of reco jets: %d\n",nRecJets);
+  if(fDebug)Printf("---- Number of reco jets: %d\n",nRecJets);
   for(Int_t i = 0; i < nRecJets; i++)
     {
       if(nRecJets == 1)