#include <TH2F.h>
#include <THnSparse.h>
#include <TSystem.h>
+#include <TObjArray.h>
#include <TInterpreter.h>
#include <TList.h>
#include <TLorentzVector.h>
ClassImp(AliAnalysisTaskJetBackgroundSubtract)
AliAnalysisTaskJetBackgroundSubtract::~AliAnalysisTaskJetBackgroundSubtract(){
+ //
+ // destructor
+ //
delete fJBArray;
delete fOutJetArrayList;
delete fInJetArrayList;
fAODOut(0x0),
fAODIn(0x0),
fAODExtension(0x0),
- fJBArray(new TObjArray()),
+ fJBArray(0),
fBackgroundBranch(""),
fNonStdFile(""),
fReplaceString1("B0"),
fAODOut(0x0),
fAODIn(0x0),
fAODExtension(0x0),
- fJBArray(new TObjArray()),
+ fJBArray(0),
fBackgroundBranch(""),
fNonStdFile(""),
fReplaceString1("B0"),
Bool_t AliAnalysisTaskJetBackgroundSubtract::Notify()
{
+ //
+ // exec with every new file
//
fAODIn = dynamic_cast<AliAODEvent*>(InputEvent());
fInJetArrayList->Clear();
fOutJetArrayList->Clear();
- for(int iJB = 0;iJB<fJBArray->GetEntries();iJB++){
+ for(int iJB = 0;iJB<(fJBArray?fJBArray->GetEntries():0);iJB++){
TObjString *ostr = (TObjString*)fJBArray->At(iJB);
if(!fInJetArrayList)fInJetArrayList =new TList();
if(!fOutJetArrayList)fOutJetArrayList =new TList();
- for(int iJB = 0;iJB<fJBArray->GetEntries();iJB++){
+ for(int iJB = 0;iJB<(fJBArray?fJBArray->GetEntries():0);iJB++){
TObjString *ostr = (TObjString*)fJBArray->At(iJB);
TString newName(ostr->GetString().Data());
if(!newName.Contains(fReplaceString1.Data())){
const Double_t xminDelta[nSparseBinsDelta] = {-120.5, 0, 0, 0};
const Double_t xmaxDelta[nSparseBinsDelta] = { 120.5, 1.0, 100,5000};
- for(int iJB = 0;iJB<fJBArray->GetEntries();iJB++){
+ for(int iJB = 0;iJB<(fJBArray?fJBArray->GetEntries():0);iJB++){
TObjString *ostr = (TObjString*)fJBArray->At(iJB);
TString oldName(ostr->GetString().Data());
TString newName(ostr->GetString().Data());
if(fBackgroundBranch.Length()==0)
AliError(Form("%s:%d No BackgroundBranch defined",(char*)__FILE__,__LINE__));
- if(fJBArray->GetEntries()==0)
+ if((fJBArray?fJBArray->GetEntries():0)==0)
AliError(Form("%s:%d No Jet Branches defined defined",(char*)__FILE__,__LINE__));
}
void AliAnalysisTaskJetBackgroundSubtract::UserExec(Option_t */*option*/)
{
+ //
+ // Execute for every selected event
+ //
if (fDebug > 1) printf("AnalysisTaskJetBackgroundSubtract::UserExec() \n");
ResetOutJets();
- if(fBackgroundBranch.Length()==0||fJBArray->GetEntries()==0){
+ if(fBackgroundBranch.Length()==0||(fJBArray?fJBArray->GetEntries():0)==0){
if(fDebug)Printf("%s:%d No background subtraction done",(char*)__FILE__,__LINE__);
PostData(1,fHistList);
}
Double_t AliAnalysisTaskJetBackgroundSubtract::RecalcRho(TClonesArray* bkgClusters,Double_t meanarea){
+ //
+ // recalc rhoo
+ //
+
Double_t ptarea=0.;
Int_t count=0;
Double_t rho=0.;
- const Double_t Rlimit2=0.8*0.8; //2*jet radius.
+ const Double_t rLimit2=0.8*0.8; //2*jet radius.
TClonesArray* jarray=0;
for(int iJB = 0;iJB<fInJetArrayList->GetEntries();iJB++){
(first->Phi()-clus->Phi())*(first->Phi()-clus->Phi());
Double_t distance2= (second->Eta()-clus->Eta())*(second->Eta()-clus->Eta())+
(second->Phi()-clus->Phi())*(second->Phi()-clus->Phi());
- if((distance1<Rlimit2)||(distance2<Rlimit2)) continue;
+ if((distance1<rLimit2)||(distance2<rLimit2)) continue;
ptarea=ptarea+clus->Pt()/clus->EffectiveAreaCharged();
count=count+1;}
if(count!=0) rho=ptarea/count;
return rho;
}
- Double_t AliAnalysisTaskJetBackgroundSubtract::RhoRC(TClonesArray* bkgClustersRC){
+Double_t AliAnalysisTaskJetBackgroundSubtract::RhoRC(TClonesArray* bkgClustersRC){
+
+ //
+ // calc rho from random cones
+ //
Double_t ptarea=0.;
Int_t count=0;
Double_t rho=0.;
- const Double_t Rlimit2=0.8*0.8; //2*jet radius.
+ const Double_t rLimit2=0.8*0.8; //2*jet radius.
TClonesArray* jarray=0;
for(int iJB = 0;iJB<fInJetArrayList->GetEntries();iJB++){
TObjString *ostr = (TObjString*)fInJetArrayList->At(iJB);
(first->Phi()-clus->Phi())*(first->Phi()-clus->Phi());
Double_t distance2= (second->Eta()-clus->Eta())*(second->Eta()-clus->Eta())+
(second->Phi()-clus->Phi())*(second->Phi()-clus->Phi());
- if((distance1<Rlimit2)||(distance2<Rlimit2)) continue;
+ if((distance1<rLimit2)||(distance2<rLimit2)) continue;
ptarea=ptarea+clus->Pt()/clus->EffectiveAreaCharged();
count=count+1;}
if(count!=0) rho=ptarea/count; }
void AliAnalysisTaskJetBackgroundSubtract::ResetOutJets(){
+ //
+ // Reset the output jets
+ //
+
if(!fOutJetArrayList)return;
for(int iJB = 0;iJB<fOutJetArrayList->GetEntries();iJB++){
TClonesArray* jarray = (TClonesArray*)fOutJetArrayList->At(iJB);
void AliAnalysisTaskJetBackgroundSubtract::PrintAODContents(){
+
+ //
+ // guess from the name what this function does
+ //
+
if(fAODIn){
Printf("%s:%d >>>>>> Input",(char*)__FILE__,__LINE__);
fAODIn->Print();
}
Int_t AliAnalysisTaskJetBackgroundSubtract::MultFromJetRefs(TClonesArray* jets){
+ //
+ // calculate multiplicty based on jet references
+ //
+
if(!jets)return 0;
Int_t refMult = 0;
return refMult;
}
+
+void AliAnalysisTaskJetBackgroundSubtract::AddJetBranch(const char* c){
+ if(!fJBArray)fJBArray = new TObjArray();
+ fJBArray->Add(new TObjString(c));
+}