coverity and rule checker fixes (B. Bathen)
authorkleinb <kleinb@f7af4fe6-9843-0410-8265-dc069ae4e863>
Wed, 25 Jan 2012 14:31:13 +0000 (14:31 +0000)
committerkleinb <kleinb@f7af4fe6-9843-0410-8265-dc069ae4e863>
Wed, 25 Jan 2012 14:31:13 +0000 (14:31 +0000)
JETAN/AliAnalysisTaskFastEmbedding.cxx
JETAN/AliAnalysisTaskFastEmbedding.h
PWGJE/AliAnalysisTaskJetResponse.cxx
PWGJE/AliAnalysisTaskJetResponse.h
PWGJE/AliAnalysisTaskJetResponseV2.cxx
PWGJE/AliAnalysisTaskJetResponseV2.h

index c092b06..60fa320 100644 (file)
@@ -548,7 +548,7 @@ void AliAnalysisTaskFastEmbedding::Init()
 //__________________________________________________________________________
 Bool_t AliAnalysisTaskFastEmbedding::UserNotify()
 {
-
+   // User defined Notify(), called once
    if(fDebug > 1) Printf("AliAnalysisTaskFastEmbedding::UserNotify()");
 
    // get total nb of events in tree (of this subjob)
@@ -564,6 +564,7 @@ Bool_t AliAnalysisTaskFastEmbedding::UserNotify()
 //__________________________________________________________________________
 void AliAnalysisTaskFastEmbedding::UserExec(Option_t *)
 {
+   // Analysis, called once per event
    if(fDebug > 1) Printf("AliAnalysisTaskFastEmbedding::UserExec()");
 
    if(!fAODout){
@@ -962,6 +963,7 @@ void AliAnalysisTaskFastEmbedding::UserExec(Option_t *)
 //__________________________________________________________________________
 void AliAnalysisTaskFastEmbedding::Terminate(Option_t *)
 {
+   // terminate
    if(fDebug > 1) Printf("AliAnalysisTaskFastEmbedding::Terminate()");
 
    if(fAODfile && fAODfile->IsOpen())
@@ -972,6 +974,8 @@ void AliAnalysisTaskFastEmbedding::Terminate(Option_t *)
 //__________________________________________________________________________
 Int_t AliAnalysisTaskFastEmbedding::GetJobID()
 {
+   // gives back the alien subjob id, if available, else -1
+
    Int_t id=-1;
 
    const char* env = gSystem->Getenv("ALIEN_PROC_ID"); // GRID
@@ -989,8 +993,9 @@ Int_t AliAnalysisTaskFastEmbedding::GetJobID()
 }
 
 //__________________________________________________________________________
-
-Int_t AliAnalysisTaskFastEmbedding::SelectAODfile(){
+Int_t AliAnalysisTaskFastEmbedding::SelectAODfile()
+{
+   // select an AOD file from fAODPathArray
 
    Int_t nFiles = fAODPathArray->GetEntries();
 
@@ -999,6 +1004,7 @@ Int_t AliAnalysisTaskFastEmbedding::SelectAODfile(){
    Float_t tmpProp = -1.;
    Int_t pendingEvents = fInputEntries-Entry();
    //Printf("input entries %d, entry %d, pending events %d", fInputEntries, (Int_t)Entry(), pendingEvents);
+   if(fAODEntriesArray){
    while(rndm->Rndm()>tmpProp){
       n = rndm->Integer(nFiles);
       fAODEntries = fAODEntriesArray->At(n);
@@ -1006,6 +1012,11 @@ Int_t AliAnalysisTaskFastEmbedding::SelectAODfile(){
       tmpProp = fAODEntriesMax ? (Float_t)tmpEntries/fAODEntriesMax : 1.;
    }
    fAODEntry = (Int_t)(rndm->Uniform(fAODEntries));
+   } else {
+      AliWarning("Number of entries in extra AODs not set!");
+      n = rndm->Integer(nFiles);
+      fAODEntry = 0;
+   }
 
    AliInfo(Form("Select AOD file %d", n));
    TObjString *objStr = (TObjString*) fAODPathArray->At(n);
@@ -1019,14 +1030,15 @@ Int_t AliAnalysisTaskFastEmbedding::SelectAODfile(){
 }
 
 //__________________________________________________________________________
-
-Int_t AliAnalysisTaskFastEmbedding::OpenAODfile(Int_t trial){
+Int_t AliAnalysisTaskFastEmbedding::OpenAODfile(Int_t trial)
+{
+   // select and open an AOD file from fAODPathArray
 
    if(fAODPathArray) fFileId = SelectAODfile();
    if(fFileId<0) return -1;
 
    TDirectory *owd = gDirectory;
-   if (fAODfile) fAODfile->Close();
+   if (fAODfile && fAODfile->IsOpen()) fAODfile->Close();
    fAODfile = TFile::Open(fAODPath.Data(),"TIMEOUT=180");
    owd->cd();
    if(!fAODfile){
@@ -1108,6 +1120,8 @@ Int_t AliAnalysisTaskFastEmbedding::OpenAODfile(Int_t trial){
 //____________________________________________________________________________
 Float_t AliAnalysisTaskFastEmbedding::GetPtHard(Bool_t bSet, Float_t newValue){
 
+   // static stored, available for other tasks in train
+
    static Float_t ptHard = -1.;
    if(bSet) ptHard = newValue;
 
@@ -1117,6 +1131,8 @@ Float_t AliAnalysisTaskFastEmbedding::GetPtHard(Bool_t bSet, Float_t newValue){
 //____________________________________________________________________________
 Int_t AliAnalysisTaskFastEmbedding::GetPtHardBin(Double_t ptHard){
 
+   // returns pt hard bin (for LHC10e14, LHC11a1x, LHC11a2x simulations)
+
    const Int_t nBins = 10;
    Double_t binLimits[nBins] = { 5., 11., 21., 36., 57., 84., 117., 156., 200., 249. }; // lower limits
 
index d17a43c..92f1611 100644 (file)
@@ -60,9 +60,9 @@ public:
    void SetJetBranch(TString name) {fJetBranch = name;}
 
    void SetEmbedMode(Int_t m) {fEmbedMode = m;}
-   Int_t GetEmbedMode() {return fEmbedMode;} 
+   Int_t GetEmbedMode() const {return fEmbedMode;} 
    void SetEvtSelecMode(Int_t s) {fEvtSelecMode = s;}
-   Int_t GetEvtSelecMode() {return fEvtSelecMode;}
+   Int_t GetEvtSelecMode() const {return fEvtSelecMode;}
 
    void SetEvtSelJetPtRange(Float_t minPt, Float_t maxPt) {fEvtSelMinJetPt = minPt; fEvtSelMaxJetPt = maxPt;}
    void SetEvtSelJetEtaRange(Float_t minEta, Float_t maxEta) {fEvtSelMinJetEta = minEta; fEvtSelMaxJetEta = maxEta;}
@@ -127,8 +127,8 @@ private:
    Int_t fAODEntry; // entry of extra AOD
    Int_t fCountEvents; // count processed events in this file
 
-   Int_t fEmbedMode;
-   Int_t fEvtSelecMode;
+   Int_t fEmbedMode;     // embedding mode: kAODFull=0, kAODJetTracks=1, kAODJet4Mom=2, kToyTracks=3
+   Int_t fEvtSelecMode;  // event selection criterion: kEventsAll=0, kEventsJetPt=1
 
    // event selection from AOD
    Float_t fEvtSelMinJetPt;       // minimum pt of the leading jet
@@ -190,7 +190,7 @@ private:
    Int_t OpenAODfile(Int_t trial = 0);
 
 
-   ClassDef(AliAnalysisTaskFastEmbedding, 5);
+   ClassDef(AliAnalysisTaskFastEmbedding, 6);
 };
 
 #endif
index bb05a55..bbf0816 100644 (file)
@@ -1,3 +1,25 @@
+/**************************************************************************
+ * Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved. *
+ *                                                                        *
+ * Author: The ALICE Off-line Project.                                    *
+ * Contributors are mentioned in the code where appropriate.              *
+ *                                                                        *
+ * Permission to use, copy, modify and distribute this software and its   *
+ * documentation strictly for non-commercial purposes is hereby granted   *
+ * without fee, provided that the above copyright notice appears in all   *
+ * copies and that both the copyright notice and this permission notice   *
+ * appear in the supporting documentation. The authors make no claims     *
+ * about the suitability of this software for any purpose. It is          *
+ * provided "as is" without express or implied warranty.                  *
+ **************************************************************************/
+
+//
+// task compares jets in two branches,
+// written for analysis of jet embedding in HI events
+//
+// newer class: AliAnalysisTaskJetResponseV2
+//
+
 #include "TChain.h"
 #include "TTree.h"
 #include "TMath.h"
index 8ff02ab..5d42b5b 100644 (file)
@@ -1,6 +1,16 @@
 #ifndef ALIANALYSISTASKJETRESPONSE_H
 #define ALIANALYSISTASKJETRESPONSE_H
 
+/* Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved. *
+ * See cxx source for full Copyright notice                               */
+
+//
+// task compares jets in two branches,
+// written for analysis of jet embedding in HI events
+//
+// newer class: AliAnalysisTaskJetResponseV2
+//
+
 class TH1F;
 class TH2F;
 class TH3F;
index d9f5dbf..9276a38 100644 (file)
@@ -1,3 +1,24 @@
+/**************************************************************************
+ * Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved. *
+ *                                                                        *
+ * Author: The ALICE Off-line Project.                                    *
+ * Contributors are mentioned in the code where appropriate.              *
+ *                                                                        *
+ * Permission to use, copy, modify and distribute this software and its   *
+ * documentation strictly for non-commercial purposes is hereby granted   *
+ * without fee, provided that the above copyright notice appears in all   *
+ * copies and that both the copyright notice and this permission notice   *
+ * appear in the supporting documentation. The authors make no claims     *
+ * about the suitability of this software for any purpose. It is          *
+ * provided "as is" without express or implied warranty.                  *
+ **************************************************************************/
+
+//
+// task compares jets in two branches,
+// written for analysis of jet embedding in HI events
+//
+
+
 #include "TChain.h"
 #include "TTree.h"
 #include "TMath.h"
@@ -716,6 +737,7 @@ void AliAnalysisTaskJetResponseV2::Terminate(const Option_t *)
 
 Int_t AliAnalysisTaskJetResponseV2::GetNInputTracks()
 {
+   // number of tracks in the event, obtained via jet finder
 
    Int_t nInputTracks = 0;
 
@@ -751,6 +773,8 @@ Int_t AliAnalysisTaskJetResponseV2::GetNInputTracks()
 
 THnSparse* AliAnalysisTaskJetResponseV2::NewTHnSparseF(const char* name, UInt_t entries, UInt_t opt)
 {
+   // generate new THnSparseF, axes are defined in GetDimParams()
+
    Int_t count = 0;
    UInt_t tmp = entries;
    while(tmp!=0){
@@ -784,6 +808,7 @@ THnSparse* AliAnalysisTaskJetResponseV2::NewTHnSparseF(const char* name, UInt_t
 
 void AliAnalysisTaskJetResponseV2::GetDimParams(Int_t iEntry, Bool_t hr, TString &label, Int_t &nbins, Double_t &xmin, Double_t &xmax)
 {
+   // stores label and binning of axis for THnSparse
 
    const Double_t pi = TMath::Pi();
    
@@ -1064,6 +1089,8 @@ void AliAnalysisTaskJetResponseV2::GetDimParams(Int_t iEntry, Bool_t hr, TString
 //____________________________________________________________________________
 Int_t AliAnalysisTaskJetResponseV2::GetPtHardBin(Double_t ptHard){
 
+   // returns pt hard bin (for LHC10e14, LHC11a1x, LHC11a2x simulations)
+
    const Int_t nBins = 10;
    Double_t binLimits[nBins] = { 5., 11., 21., 36., 57., 84., 117., 156., 200., 249. }; // lower limits
    
@@ -1078,6 +1105,8 @@ Int_t AliAnalysisTaskJetResponseV2::GetPtHardBin(Double_t ptHard){
 //____________________________________________________________________________
 Double_t AliAnalysisTaskJetResponseV2::GetPt(AliAODJet *j, Int_t mode=0){
 
+   // returns jet pt, also negative pt after background subtraction if available
+
    Double_t pt = 0.;
 
    if(fKeepJets && mode==1){  // background subtracted pt, can be negative
index b847223..74d2356 100644 (file)
@@ -1,6 +1,14 @@
 #ifndef ALIANALYSISTASKJETRESPONSEV2_H
 #define ALIANALYSISTASKJETRESPONSEV2_H
 
+/* Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved. *
+ * See cxx source for full Copyright notice                               */
+
+//
+// task compares jets in two branches,
+// written for analysis of jet embedding in HI events
+//
+
 class TH1F;
 class TH2F;
 class TH3F;
@@ -87,7 +95,7 @@ private:
    TString fJetBranchName[2]; //  name of jet branches to compare
    TList *fListJets[2];       //! jet lists
 
-   TString fBackgroundBranch;
+   TString fBackgroundBranch;  // branch of external background
 
    // event selection
    Bool_t fIsPbPb;         // is Pb-Pb (fast embedding) or p-p (detector response)