fh2ESDTriggerRun->Fill(run,iTrig+2);
fh2VtxXY->Fill(xvtx,yvtx);
}
- if(TMath::Abs(zvtx)<fZVtxCut&&TMath::Abs(xvtx)<0.5&&TMath::Abs(yvtx)<0.5)fh2ESDTriggerRun->Fill(run,iTrig+3);
+ Float_t r2 = xvtx *xvtx + yvtx *yvtx;
+ if(TMath::Abs(zvtx)<fZVtxCut&&r2<1)fh2ESDTriggerRun->Fill(run,iTrig+3);
}
else{
fh2ESDTriggerRun->Fill(run,0);
Float_t yvtx = vtxAOD->GetY();
Float_t xvtx = vtxAOD->GetX();
fh2TriggerVtx->Fill(it,zvtx);
- if(TMath::Abs(zvtx)<fZVtxCut&&aodTrig&&TMath::Abs(xvtx)<0.5&&TMath::Abs(yvtx)<0.5){
+ Float_t r2 = yvtx*yvtx+xvtx*xvtx;
+ if(TMath::Abs(zvtx)<fZVtxCut&&aodTrig&&r2<1){
fh2TriggerCount->Fill(it,kTriggeredVertexIn);
}
}
fh2JetsLeadingPhiPt(0x0),
fh2TracksLeadingPhiEta(0x0),
fh2TracksLeadingPhiPt(0x0),
+ fh2TracksLeadingJetPhiPt(0x0),
fHistList(0x0)
{
for(int i = 0;i < kMaxStep*2;++i){
fh2JetsLeadingPhiPt(0x0),
fh2TracksLeadingPhiEta(0x0),
fh2TracksLeadingPhiPt(0x0),
+ fh2TracksLeadingJetPhiPt(0x0),
fHistList(0x0)
{
binLimitsEta[iEta] = -2.0;
}
else{
- binLimitsEta[iEta] = binLimitsEta[iEta-1] + 1/(Float_t)nBinEta + 0.1;
+ binLimitsEta[iEta] = binLimitsEta[iEta-1] + 0.1;
}
}
}
fh1Trials->Fill("#sum{ntrials}",fAvgTrials);
-
+ // if(fDebug>0)aodH->SetFillAOD(kFALSE);
if (fDebug > 10)Printf("%s:%d",(char*)__FILE__,__LINE__);
if((fAnalysisType&kAnaMCESD)==kAnaMCESD){
// this is the part we only use when we have MC information
Float_t phi = leading->Phi();
if(phi<0)phi+=TMath::Pi()*2.;
Float_t eta = leading->Eta();
+ pt = leading->Pt();
while((tmpRec = (AliAODJet*)(recIter->Next()))){
Float_t tmpPt = tmpRec->Pt();
fh1PtJetsRecIn->Fill(tmpPt);
Float_t phi = leading->Phi();
if(phi<0)phi+=TMath::Pi()*2.;
Float_t eta = leading->Eta();
+ pt = leading->Pt();
while((tmpRec = (AliVParticle*)(recIter->Next()))){
Float_t tmpPt = tmpRec->Pt();
fh1PtTracksRecIn->Fill(tmpPt);
fh1SumPtTrackAreaRec->Fill(sumPt*0.4*0.4/(2.*1.8),eventW);
fh1SumPtTrackRec->Fill(sumPt,eventW);
-
+
// loop over reconstructed jets
for(int ir = 0;ir < nRecJets;++ir){
Double_t etaRec = recJets[ir].Eta();
// 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());
+ // aodH->SetFillAOD(kTRUE);
fAOD->GetHeader()->Print();
Printf("TriggerClasses: %s",fAOD->GetFiredTriggerClasses().Data());
for(int it = 0;it < fAOD->GetNumberOfTracks();++it){