histo switch
authorloizides <loizides@f7af4fe6-9843-0410-8265-dc069ae4e863>
Mon, 21 May 2012 19:46:51 +0000 (19:46 +0000)
committerloizides <loizides@f7af4fe6-9843-0410-8265-dc069ae4e863>
Mon, 21 May 2012 19:46:51 +0000 (19:46 +0000)
PWGGA/EMCALJetTasks/AliAnalysisTaskRho.cxx
PWGGA/EMCALJetTasks/AliAnalysisTaskRho.h
PWGGA/EMCALJetTasks/AliHadCorrTask.cxx
PWGGA/EMCALJetTasks/AliHadCorrTask.h
PWGGA/EMCALJetTasks/macros/AddTaskHadCorr.C
PWGGA/EMCALJetTasks/macros/AddTaskRho.C

index 54e564e..b4dff30 100644 (file)
@@ -57,6 +57,7 @@ AliAnalysisTaskRho::AliAnalysisTaskRho() :
 {
   // Constructor
 }
+
 //________________________________________________________________________
 AliAnalysisTaskRho::AliAnalysisTaskRho(const char *name) :
   AliAnalysisTaskRhoBase(name),
@@ -92,10 +93,49 @@ AliAnalysisTaskRho::AliAnalysisTaskRho(const char *name) :
 {
   // Constructor
 
-  DefineOutput(1, TList::Class());
 }
 
 //________________________________________________________________________
+AliAnalysisTaskRho::AliAnalysisTaskRho(const char *name, Bool_t histo) :
+  AliAnalysisTaskRhoBase(name),
+  fTracksName("tracks"),
+  fJetsName("KtJets"),
+  fClustersName("caloClusters"),
+  fRhoScaledName(""),
+  fPhiMin(0),
+  fPhiMax(0),
+  fEtaMin(0),
+  fEtaMax(0),
+  fAreaCut(0),
+  fNExclLeadJets(0),
+  fScaleFunction(0),
+  fCreateHisto(histo),
+  fOutputList(0),
+  fHistCentrality(0),
+  fHistJetPt(0),
+  fHistJetArea(0),
+  fHistRhovsCent(0),
+  fHistDeltaRhovsCent(0),
+  fHistDeltaRhoScalevsCent(0),
+  fHistJetPtvsCent(0),
+  fHistJetAreavsCent(0),
+  fHistNjetvsCent(0),
+  fHistRhovsNtrack(0),
+  fHistDeltaRhovsNtrack(0),
+  fHistDeltaRhoScalevsNtrack(0),
+  fHistJetPtvsNtrack(0),
+  fHistJetAreavsNtrack(0),
+  fHistNjetvsNtrack(0),
+  fRhoScaled(0)
+{
+  // Constructor
+
+  if (fCreateHisto)
+    DefineOutput(1, TList::Class());
+}
+
+
+//________________________________________________________________________
 void AliAnalysisTaskRho::UserCreateOutputObjects()
 {
   AliAnalysisTaskRhoBase::UserCreateOutputObjects();
@@ -104,6 +144,9 @@ void AliAnalysisTaskRho::UserCreateOutputObjects()
   fRhoScaledName += "_Scaled";
   fRhoScaled = new TParameter<Double_t>(fRhoScaledName, 0);  
 
+  if (!fCreateHisto)
+    return;
+
   OpenFile(1);
   fOutputList = new TList();
   fOutputList->SetOwner();
index db14abe..bf9e87e 100644 (file)
@@ -19,6 +19,7 @@ class AliAnalysisTaskRho : public AliAnalysisTaskRhoBase {
  public:
   AliAnalysisTaskRho();
   AliAnalysisTaskRho(const char *name);
+  AliAnalysisTaskRho(const char *name, Bool_t histo);
   virtual ~AliAnalysisTaskRho() {}
   
   virtual void           UserCreateOutputObjects();
@@ -32,7 +33,6 @@ class AliAnalysisTaskRho : public AliAnalysisTaskRhoBase {
   void                   SetScaleFunction(TF1* sf)                             { fScaleFunction = sf   ; }
   void                   SetTracksName(const char *n)                          { fTracksName    = n    ; }
   void                   SetExcludeLeadJets(UInt_t n)                          { fNExclLeadJets = n    ; }
-  void                   SetCreateHistograms(Bool_t h = kTRUE)                 { fCreateHisto   = h    ; }
   
  protected:
   virtual Double_t       GetScaleFactor(Double_t cent);
index 3ab26df..dbf561c 100644 (file)
@@ -106,8 +106,52 @@ AliHadCorrTask::AliHadCorrTask(const char *name) :
   } 
 
   fBranchNames="ESD:AliESDRun.,AliESDHeader.,PrimaryVertex.";
+}
+
+//________________________________________________________________________
+AliHadCorrTask::AliHadCorrTask(const char *name, Bool_t histo) : 
+  AliAnalysisTaskSE(name),
+  fTracksName("Tracks"),
+  fCaloName("CaloClusters"),
+  fOutCaloName("CaloClustersCorr"),
+  fPhiMatch(0.05),
+  fEtaMatch(0.025),
+  fDoTrackClus(0),
+  fHadCorr(0),
+  fMinPt(0.15),
+  fCreateHisto(histo),
+  fOutClusters(0),
+  fOutputList(0),
+  fHistNclusvsCent(0),
+  fHistNclusMatchvsCent(0),
+  fHistEbefore(0),
+  fHistEafter(0),
+  fHistEoPCent(0),
+  fHistNMatchCent(0),
+  fHistNMatchCent_trk(0),
+  fHistCentrality(0)
 
-  DefineOutput(1, TList::Class());
+{
+  // Standard constructor.
+
+   for(Int_t i=0; i<8; i++) {
+    if (i<4) {
+      fHistMatchEvsP[i]   = 0;
+      fHistMatchdRvsEP[i] = 0;
+      for(Int_t j=0; j<3; j++) {
+       fHistEsubPch[i][j] = 0;
+       fHistEsubPchRat[i][j] = 0;
+      }
+    }
+    for(Int_t j=0; j<9; j++) {
+      fHistMatchEtaPhi[i][j] = 0;
+    }
+  } 
+
+  fBranchNames="ESD:AliESDRun.,AliESDHeader.,PrimaryVertex.";
+
+  if (fCreateHisto)
+    DefineOutput(1, TList::Class());
 }
 
 //________________________________________________________________________
@@ -385,6 +429,9 @@ void AliHadCorrTask::UserCreateOutputObjects()
   }
   fOutClusters->SetName(fOutCaloName);
 
+  if (!fCreateHisto)
+    return;
+
   OpenFile(1);
   fOutputList = new TList();
   fOutputList->SetOwner();
index a6670cb..86168b9 100644 (file)
@@ -14,6 +14,7 @@ class AliHadCorrTask : public AliAnalysisTaskSE {
  public:
   AliHadCorrTask();
   AliHadCorrTask(const char *name); 
+  AliHadCorrTask(const char *name, Bool_t histo); 
   virtual ~AliHadCorrTask();
 
   void         UserCreateOutputObjects();
@@ -28,7 +29,6 @@ class AliHadCorrTask : public AliAnalysisTaskSE {
   void         SetPhiMatch(Double_t phi)               { fPhiMatch       = phi; }
   void         SetTracksName(const char *n)            { fTracksName     = n;   }
   void         SetTrackClus(Int_t c)                   { fDoTrackClus    = c;   }
-  void         SetCreateHistograms(Bool_t h = kTRUE)   { fCreateHisto    = h;   }
 
  protected:
   Int_t        GetCentBin(Double_t cent) const;
index 6d775b5..23e9bc9 100644 (file)
@@ -8,7 +8,8 @@ AliHadCorrTask* AddTaskHadCorr(
   const Double_t minPt       = 0.15,
   const Double_t phiMatch    = 0.050,
   const Double_t etaMatch    = 0.025,
-  const char *outputname     = "AnalysisResults.root"
+  const char *outputname     = "AnalysisResults.root",
+  const Bool_t   histo       = kFALSE
 )
 {  
   // Get the pointer to the existing analysis manager via the static access method.
@@ -33,7 +34,7 @@ AliHadCorrTask* AddTaskHadCorr(
   //-------------------------------------------------------
 
   TString name(Form("HadCorr_%s", outClusName));
-  AliHadCorrTask *hcor = new AliHadCorrTask(name);
+  AliHadCorrTask *hcor = new AliHadCorrTask(name, histo);
   hcor->SetTracksName(nTracks);
   hcor->SetClusName(nClusters);
   hcor->SetOutClusName(outClusName);
index 0eb4349..365da93 100644 (file)
@@ -10,7 +10,8 @@ AliAnalysisTaskRho* AddTaskRho(
    const Double_t minEta      = -0.3,
    const Double_t maxEta      = 0.3,
    const Double_t minArea     = 0.0,
-   const UInt_t   exclJets    = 0
+   const UInt_t   exclJets    = 0,
+   const Bool_t   histo       = kFALSE
 )
 {  
   // Get the pointer to the existing analysis manager via the static access method.
@@ -35,7 +36,7 @@ AliAnalysisTaskRho* AddTaskRho(
   //-------------------------------------------------------
 
   TString name(Form("Rho_%s", nJets));
-  AliAnalysisTaskRho *rhotask = new AliAnalysisTaskRho(name);
+  AliAnalysisTaskRho *rhotask = new AliAnalysisTaskRho(name, histo);
   rhotask->SetJetsName(nJets);
   rhotask->SetTracksName(nTracks);
   rhotask->SetRhoName(nRho);