// $Id$
+//
+// Scale task.
+//
+//
-#include "TChain.h"
-#include "TTree.h"
-#include "TList.h"
-#include "TH1.h"
-#include "TH2.h"
+#include <TChain.h>
#include <TClonesArray.h>
-#include <TParticle.h>
+#include <TH1.h>
+#include <TH2.h>
+#include <TList.h>
#include <TLorentzVector.h>
-
-#include "AliAnalysisTask.h"
+#include <TParticle.h>
+#include <TTree.h>
#include "AliAnalysisManager.h"
-#include "AliVCluster.h"
-
+#include "AliAnalysisTask.h"
+#include "AliCentrality.h"
#include "AliESDEvent.h"
#include "AliESDInputHandler.h"
-#include "AliCentrality.h"
-
+#include "AliVCluster.h"
#include "AliAnalysisTaskScale.h"
ClassImp(AliAnalysisTaskScale)
//________________________________________________________________________
AliAnalysisTaskScale::AliAnalysisTaskScale(const char *name)
- : AliAnalysisTaskSE(name), fESD(0), fOutputList(0), fHistCentrality(0), fHistPtTPCvsCent(0), fHistPtEMCALvsCent(0), fHistEtvsCent(0), fHistScalevsCent(0), fHistDeltaScalevsCent(0), fHistPtTPCvsNtrack(0), fHistPtEMCALvsNtrack(0), fHistEtvsNtrack(0), fHistScalevsNtrack(0), fHistDeltaScalevsNtrack(0),
+ : AliAnalysisTaskSE(name),
fTracksName("tracks"),
- fClustersName("clusters")
-
+ fClustersName("clusters"),
+ fESD(0),
+ fOutputList(0),
+ fHistCentrality(0),
+ fHistPtTPCvsCent(0),
+ fHistPtEMCALvsCent(0),
+ fHistEtvsCent(0),
+ fHistScalevsCent(0),
+ fHistDeltaScalevsCent(0),
+ fHistPtTPCvsNtrack(0),
+ fHistPtEMCALvsNtrack(0),
+ fHistEtvsNtrack(0),
+ fHistScalevsNtrack(0),
+ fHistDeltaScalevsNtrack(0)
{
- // Constructor
+ // Constructor.
DefineInput(0, TChain::Class());
DefineOutput(1, TList::Class());
//________________________________________________________________________
void AliAnalysisTaskScale::UserCreateOutputObjects()
{
+ // Create my user objects.
+
OpenFile(1);
fOutputList = new TList();
fOutputList->SetOwner();
fHistScalevsCent = new TH2F("fHistScalevsCent","rho vs cent",101,-1,100,400,0,4.0);
fHistDeltaScalevsCent = new TH2F("fHistDeltaScalevsCent","rho vs cent",101,-1,100,400,-2.0,2.0);
-
fHistPtTPCvsNtrack = new TH2F("fHistPtTPCvsNtrack","rho vs cent",500,0,2500,500,0,1000);
fHistPtEMCALvsNtrack = new TH2F("fHistPtEMCALvsNtrack","rho vs cent",500,0,2500,500,0,1000);
fHistEtvsNtrack = new TH2F("fHistEtvsNtrack","rho vs cent",500,0,2500,500,0,1000);
fOutputList->Add(fHistEtvsNtrack);
fOutputList->Add(fHistScalevsNtrack);
fOutputList->Add(fHistDeltaScalevsNtrack);
-
PostData(1, fOutputList);
}
//________________________________________________________________________
void AliAnalysisTaskScale::UserExec(Option_t *)
{
+ // Execute on each event.
+
fESD = dynamic_cast<AliESDEvent*>(InputEvent());
if (!fESD) {
printf("ERROR: fESD not available\n");
return;
}
- AliAnalysisManager *am = AliAnalysisManager::GetAnalysisManager();
TClonesArray *tracks = 0;
TClonesArray *clusters = 0;
AliCentrality *centrality = 0;
AliError(Form("Pointer to clusters %s == 0", fClustersName.Data() ));
return;
}
-
centrality = dynamic_cast<AliCentrality*>(l->FindObject("Centrality"));
float fCent = centrality->GetCentralityPercentile("V0M");
float Et = 0;
InputEvent()->GetPrimaryVertex()->GetXYZ(vertex);
const Int_t Nclus = clusters->GetEntries();
- for (Int_t iClus = 0, iN = 0; iClus < Nclus; ++iClus) {
+ for (Int_t iClus = 0; iClus < Nclus; ++iClus) {
AliVCluster *c = dynamic_cast<AliVCluster*>(clusters->At(iClus));
if (!c->IsEMCAL())
//________________________________________________________________________
void AliAnalysisTaskScale::Terminate(Option_t *)
-{}
-
-
-
-
+{
+ // Nothing to be done for the moment.
+}
class TH1F;
class TH1I;
class TH2F;
-//class TH3F;
class AliESDEvent;
#include "AliAnalysisTaskSE.h"
class AliAnalysisTaskScale : public AliAnalysisTaskSE {
public:
- AliAnalysisTaskScale() : AliAnalysisTaskSE(), fESD(0), fOutputList(0), fHistCentrality(0), fHistPtTPCvsCent(0), fHistPtEMCALvsCent(0), fHistEtvsCent(0), fHistScalevsCent(0), fHistDeltaScalevsCent(0), fHistPtTPCvsNtrack(0), fHistPtEMCALvsNtrack(0), fHistEtvsNtrack(0), fHistScalevsNtrack(0), fHistDeltaScalevsNtrack(0){}
+ AliAnalysisTaskScale() : AliAnalysisTaskSE(), fTracksName("tracks"), fClustersName("clusters"), fESD(0),
+ fOutputList(0), fHistCentrality(0), fHistPtTPCvsCent(0), fHistPtEMCALvsCent(0), fHistEtvsCent(0),
+ fHistScalevsCent(0), fHistDeltaScalevsCent(0), fHistPtTPCvsNtrack(0), fHistPtEMCALvsNtrack(0),
+ fHistEtvsNtrack(0), fHistScalevsNtrack(0), fHistDeltaScalevsNtrack(0) {}
AliAnalysisTaskScale(const char *name);
virtual ~AliAnalysisTaskScale() {}
virtual void UserCreateOutputObjects();
virtual void UserExec(Option_t *option);
virtual void Terminate(Option_t *);
- virtual void SetTracksName(const char *n) { fTracksName = n; }
- virtual void SetClustersName(const char *n) { fClustersName = n; }
-
+ virtual void SetTracksName(const char *n) { fTracksName = n; }
+ virtual void SetClustersName(const char *n) { fClustersName = n; }
private:
TString fTracksName; // name of track collection
- TString fJetsName; // name of jet collection
- TString fClustersName; // name of clusters collection
- AliESDEvent *fESD; //! ESD object
- TList *fOutputList; //! Output list
- TH1F *fHistCentrality;
-
- TH2F *fHistPtTPCvsCent;
- TH2F *fHistPtEMCALvsCent;
- TH2F *fHistEtvsCent;
- TH2F *fHistScalevsCent;
- TH2F *fHistDeltaScalevsCent;
-
- TH2F *fHistPtTPCvsNtrack;
- TH2F *fHistPtEMCALvsNtrack;
- TH2F *fHistEtvsNtrack;
- TH2F *fHistScalevsNtrack;
- TH2F *fHistDeltaScalevsNtrack;
-
-
+ TString fClustersName; // name of clusters collection
+ AliESDEvent *fESD; //!ESD object
+ TList *fOutputList; //!Output list
+ TH1F *fHistCentrality; //!
+ TH2F *fHistPtTPCvsCent; //!
+ TH2F *fHistPtEMCALvsCent; //!
+ TH2F *fHistEtvsCent; //!
+ TH2F *fHistScalevsCent; //!
+ TH2F *fHistDeltaScalevsCent; //!
+ TH2F *fHistPtTPCvsNtrack; //!
+ TH2F *fHistPtEMCALvsNtrack; //!
+ TH2F *fHistEtvsNtrack; //!
+ TH2F *fHistScalevsNtrack; //!
+ TH2F *fHistDeltaScalevsNtrack; //!
AliAnalysisTaskScale(const AliAnalysisTaskScale&); // not implemented
AliAnalysisTaskScale& operator=(const AliAnalysisTaskScale&); // not implemented
ClassDef(AliAnalysisTaskScale, 1); // example of analysis
};
-
#endif