#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));
+}
#include "AliAnalysisTaskSE.h"
#include "TObjString.h"
#include "TString.h"
-#include "TObjArray.h"
////////////////
class TProfile;
class TRandom3;
class TRefArray;
-
+class TObjArray;
class AliAnalysisTaskJetBackgroundSubtract : public AliAnalysisTaskSE
{
virtual Bool_t Notify();
// Task specific methods...
- virtual void AddJetBranch(const char* c){fJBArray->Add(new TObjString(c));}
+ virtual void AddJetBranch(const char* c);
virtual void SetSubtractionMethod(Int_t i){fSubtraction = i;}
virtual Int_t GetSubtractionMethod(){return fSubtraction;}
virtual void SetKeepJets(Bool_t b = kTRUE){fKeepJets = b;}
ClassImp(AliAnalysisTaskJetCluster)
AliAnalysisTaskJetCluster::~AliAnalysisTaskJetCluster(){
+ //
+ // Destructor
+ //
delete fRef;
delete fRandom;
fh2TracksLeadingJetPhiPtWRan(0x0),
fHistList(0x0)
{
+ //
+ // Constructor
+ //
+
for(int i = 0;i<3;i++){
fh1BiARandomCones[i] = 0;
fh1BiARandomConesRan[i] = 0;
fh2TracksLeadingJetPhiPtWRan(0x0),
fHistList(0x0)
{
+ //
+ // named ctor
+ //
for(int i = 0;i<3;i++){
fh1BiARandomCones[i] = 0;
fh1BiARandomConesRan[i] = 0;
void AliAnalysisTaskJetCluster::Terminate(Option_t */*option*/)
{
-// Terminate analysis
-//
- if (fDebug > 1) printf("AnalysisJetCluster: Terminate() \n");
+ //
+ // Terminate analysis
+ //
+ if (fDebug > 1) printf("AnalysisJetCluster: Terminate() \n");
}
Int_t AliAnalysisTaskJetCluster::GetListOfTracks(TList *list,Int_t type){
+ //
+ // get list of tracks/particles for different types
+ //
+
if(fDebug>2)Printf("%s:%d Selecting tracks with %d",(char*)__FILE__,__LINE__,type);
Int_t iCount = 0;
// Fast jet
- Double_t fRparam;
+ Double_t fRparam; // fastjet distance parameter
fastjet::JetAlgorithm fAlgorithm; //fastjet::kt_algorithm
fastjet::Strategy fStrategy; //= fastjet::Best;
fastjet::RecombinationScheme fRecombScheme; // = fastjet::BIpt_scheme;
- fastjet::AreaType fAreaType;
- Double_t fGhostArea;
- Int_t fActiveAreaRepeats;
- Double_t fGhostEtamax;
+ fastjet::AreaType fAreaType; // fastjet area type
+ Double_t fGhostArea; // fasjet ghost area
+ Int_t fActiveAreaRepeats; // fast jet active area repeats
+ Double_t fGhostEtamax; // fast jet ghost area
TClonesArray *fTCAJetsOut; //! TCA of output jets
TClonesArray *fTCAJetsOutRan; //! TCA of output jets in randomized event
+//
+// Taks for jet finding in ALICE
+//
#ifndef ALIANALYSISTASKJETS_H
#define ALIANALYSISTASKJETS_H
+//
+// CDF jet finding
+//
#ifndef ALIANALYSISTASKJETSCDF_H
#define ALIANALYSISTASKJETSCDF_H
#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));
+}
#include "AliAnalysisTaskSE.h"
#include "TObjString.h"
#include "TString.h"
-#include "TObjArray.h"
////////////////
class TProfile;
class TRandom3;
class TRefArray;
-
+class TObjArray;
class AliAnalysisTaskJetBackgroundSubtract : public AliAnalysisTaskSE
{
virtual Bool_t Notify();
// Task specific methods...
- virtual void AddJetBranch(const char* c){fJBArray->Add(new TObjString(c));}
+ virtual void AddJetBranch(const char* c);
virtual void SetSubtractionMethod(Int_t i){fSubtraction = i;}
virtual Int_t GetSubtractionMethod(){return fSubtraction;}
virtual void SetKeepJets(Bool_t b = kTRUE){fKeepJets = b;}
ClassImp(AliAnalysisTaskJetCluster)
AliAnalysisTaskJetCluster::~AliAnalysisTaskJetCluster(){
+ //
+ // Destructor
+ //
delete fRef;
delete fRandom;
fh2TracksLeadingJetPhiPtWRan(0x0),
fHistList(0x0)
{
+ //
+ // Constructor
+ //
+
for(int i = 0;i<3;i++){
fh1BiARandomCones[i] = 0;
fh1BiARandomConesRan[i] = 0;
fh2TracksLeadingJetPhiPtWRan(0x0),
fHistList(0x0)
{
+ //
+ // named ctor
+ //
for(int i = 0;i<3;i++){
fh1BiARandomCones[i] = 0;
fh1BiARandomConesRan[i] = 0;
void AliAnalysisTaskJetCluster::Terminate(Option_t */*option*/)
{
-// Terminate analysis
-//
- if (fDebug > 1) printf("AnalysisJetCluster: Terminate() \n");
+ //
+ // Terminate analysis
+ //
+ if (fDebug > 1) printf("AnalysisJetCluster: Terminate() \n");
}
Int_t AliAnalysisTaskJetCluster::GetListOfTracks(TList *list,Int_t type){
+ //
+ // get list of tracks/particles for different types
+ //
+
if(fDebug>2)Printf("%s:%d Selecting tracks with %d",(char*)__FILE__,__LINE__,type);
Int_t iCount = 0;
// Fast jet
- Double_t fRparam;
+ Double_t fRparam; // fastjet distance parameter
fastjet::JetAlgorithm fAlgorithm; //fastjet::kt_algorithm
fastjet::Strategy fStrategy; //= fastjet::Best;
fastjet::RecombinationScheme fRecombScheme; // = fastjet::BIpt_scheme;
- fastjet::AreaType fAreaType;
- Double_t fGhostArea;
- Int_t fActiveAreaRepeats;
- Double_t fGhostEtamax;
+ fastjet::AreaType fAreaType; // fastjet area type
+ Double_t fGhostArea; // fasjet ghost area
+ Int_t fActiveAreaRepeats; // fast jet active area repeats
+ Double_t fGhostEtamax; // fast jet ghost area
TClonesArray *fTCAJetsOut; //! TCA of output jets
TClonesArray *fTCAJetsOutRan; //! TCA of output jets in randomized event