1 - Rename AliAnaMaker and AliAnaBaseClass to AliAnaPartCorrMaker and AliAnaPartCorrB...
[u/mrichter/AliRoot.git] / PWG4 / AliAnaScale.cxx
1 /**************************************************************************
2  * Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved. *
3  *                                                                        *
4  * Author: The ALICE Off-line Project.                                    *
5  * Contributors are mentioned in the code where appropriate.              *
6  *                                                                        *
7  * Permission to use, copy, modify and distribute this software and its   *
8  * documentation strictly for non-commercial purposes is hereby granted   *
9  * without fee, provided that the above copyright notice appears in all   *
10  * copies and that both the copyright notice and this permission notice   *
11  * appear in the supporting documentation. The authors make no claims     *
12  * about the suitability of this software for any purpose. It is          *
13  * provided "as is" without express or implied warranty.                  *
14  **************************************************************************/
15
16 /* $Id: */
17
18 //_________________________________________________________________________
19 // A basic analysis task to analyse photon detected by PHOS
20 //
21 //*-- Yves Schutz 
22 //////////////////////////////////////////////////////////////////////////////
23
24 #include <TFile.h> 
25 #include <TH1.h>
26 #include <TH1F.h>
27 #include <TH1I.h>
28 #include <TCanvas.h> 
29
30 #include "AliAnaScale.h" 
31 #include "AliAnalysisManager.h"
32
33 //______________________________________________________________________________
34 AliAnaScale::AliAnaScale() : 
35   fDebug(0),
36   fScale(1.0),
37   fInputList(0x0), 
38   fOutputList(0x0),
39   fSumw2(0),
40   fhCount() 
41 {
42   //Default constructor
43 }
44 //______________________________________________________________________________
45 AliAnaScale::AliAnaScale(const char *name) : 
46   AliAnalysisTask(name,""),  
47   fDebug(0),
48   fScale(1.0), 
49   fInputList(0x0), 
50   fOutputList(0x0), 
51   fSumw2(0),
52   fhCount(0) 
53 {
54   // Constructor.
55   // Called only after the event loop
56   SetPostEventLoop(kTRUE);
57   // Input slot #0 
58   DefineInput(0,  TList::Class()) ; 
59   // Output slot 
60   DefineOutput(0,  TList::Class()) ; 
61 }
62
63 //______________________________________________________________________________
64 AliAnaScale::~AliAnaScale()
65 {
66   // dtor
67   
68 }
69
70
71 //______________________________________________________________________________
72 void AliAnaScale::ConnectInputData(const Option_t*)
73 {
74   // Initialisation of branch container and histograms 
75     
76   if(fDebug > 1) printf("*** Initialization of %s \n", GetName()) ; 
77   fInputList     = dynamic_cast<TList*>(GetInputData(0)) ;  
78 }
79 //________________________________________________________________________
80 void AliAnaScale::CreateOutputObjects()
81 {  
82   // Create the outputs containers
83
84   fOutputList = new TList() ; 
85   fOutputList->SetName(GetName()) ; 
86
87   fhCount =new TH1F("hCount","count files",1,0,1);  
88   fOutputList->Add(fhCount);
89
90 }
91
92 //______________________________________________________________________________
93 void AliAnaScale::Exec(Option_t *) 
94 {
95   // Do the Scaling
96
97   if(fDebug > 0 ) printf(">>>>> Scaling factor %e, do Sumw2 %d <<<<< \n",fScale,fSumw2) ;
98
99   TIter next(fInputList) ;      
100   TObject * h ; 
101   while ( (h = next()) ) { 
102     if(h){
103       if ( strcmp(h->ClassName(),"TNtuple") ) {
104       char name[128] ; 
105       sprintf(name, "%sScaled", h->GetName()) ; 
106       TH1 * hout = dynamic_cast<TH1*> (h->Clone(name)) ; 
107      
108       if(fSumw2) hout->Sumw2();
109       hout->Scale(fScale) ;  
110       fOutputList->Add(hout) ; 
111       } 
112       else  fOutputList->Add(h) ; 
113     }
114   }
115   // number of files
116
117   //File scaled, needed for file merging on grid
118   fhCount->Fill(0);
119  
120   PostData(0, fOutputList);
121 }
122
123
124 //______________________________________________________________________________
125 void AliAnaScale::Init()
126 {
127   // Intialisation of parameters
128   if(fDebug > 0 )printf("No initialization in scale class \n") ;
129
130 }
131
132 //______________________________________________________________________________
133 void AliAnaScale::Terminate(Option_t *)
134 {
135   // Processing when the event loop is ended
136   
137
138 }