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...
}
}
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){
// 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();
//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;
// 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)
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)