//___________________________________________________________________
Float_t AliJetBkg::BkgFastJet(){
-
+// Return background
AliFastJetHeaderV1 *header = (AliFastJetHeaderV1*)fHeader;
Bool_t debug = header->GetDebug(); // debug option
}
//___________________________________________________________________
Float_t AliJetBkg::BkgChargedFastJet(){
-
+// Background for charged jets
AliFastJetHeaderV1 *header = (AliFastJetHeaderV1*)fHeader;
Bool_t debug = header->GetDebug(); // debug option
Float_t rCone=0.4;
if(nJ==1) {
- AliAODJet *jettmp = dynamic_cast<AliAODJet*>(fAODJets->At(0));
+ AliAODJet *jettmp = (AliAODJet*)(fAODJets->At(0));
jeteta=jettmp->Eta();
jetphi=jettmp->Phi();
acc=EmcalAcceptance(jeteta,jetphi,rCone);
if(nJ>=2) {
- AliAODJet *jettmp = dynamic_cast<AliAODJet*>(fAODJets->At(0));
- AliAODJet *jettmp1 = dynamic_cast<AliAODJet*>(fAODJets->At(1));
+ AliAODJet *jettmp = (AliAODJet*)(fAODJets->At(0));
+ AliAODJet *jettmp1 = (AliAODJet*)(fAODJets->At(1));
jeteta=jettmp->Eta();
jetphi=jettmp->Phi();
jeteta1=jettmp1->Eta();
}
-Bool_t AliJetBkg::EmcalAcceptance(const Float_t eta, const Float_t phi, const Float_t radius){
-
+Bool_t AliJetBkg::EmcalAcceptance(const Float_t eta, const Float_t phi, const Float_t radius) const{
+// Apply emcal acceptance cuts
Float_t meanPhi=190./180.*TMath::Pi()-110./180.*TMath::Pi()/2;
Float_t deltaphi=110./180.*TMath::Pi();
Float_t phicut=deltaphi/2.-radius;