6dd1ab7a5a7aec3fef37fa04bd0cf68943d36bb5
[u/mrichter/AliRoot.git] / PWG / CaloTrackCorrBase / 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 //Root system
25 #include <TH1.h>
26 #include <TH1F.h>
27 //#include "Riostream.h"
28
29 //Analysis system
30 #include "AliAnaScale.h" 
31
32 //______________________________________________________________________________
33 AliAnaScale::AliAnaScale() : 
34   fDebug(0),
35   fScale(1.0),
36   fInputList(0x0), 
37   fOutputList(0x0),
38   fSumw2(0),
39   fhCount() 
40 {
41   //Default constructor
42 }
43 //______________________________________________________________________________
44 AliAnaScale::AliAnaScale(const char *name) : 
45   AliAnalysisTask(name,""),  
46   fDebug(0),
47   fScale(1.0), 
48   fInputList(0x0), 
49   fOutputList(0x0), 
50   fSumw2(0),
51   fhCount(0) 
52 {
53   // Constructor.
54   // Called only after the event loop
55   SetPostEventLoop(kTRUE);
56   // Input slot #0 
57   DefineInput(0,  TList::Class()) ; 
58   // Output slot 
59   DefineOutput(0,  TList::Class()) ; 
60 }
61
62 //______________________________________________________________________________
63 AliAnaScale::~AliAnaScale()
64 {
65   // dtor
66   
67 }
68
69
70 //______________________________________________________________________________
71 void AliAnaScale::ConnectInputData(const Option_t*)
72 {
73   // Initialisation of branch container and histograms 
74     
75   if(fDebug > 1) printf("*** Initialization of %s \n", GetName()) ; 
76   fInputList     = dynamic_cast<TList*>(GetInputData(0)) ;  
77 }
78 //________________________________________________________________________
79 void AliAnaScale::CreateOutputObjects()
80 {  
81   // Create the outputs containers
82
83   fOutputList = new TList() ; 
84   fOutputList->SetName(GetName()) ; 
85
86   fhCount =new TH1F("hCount","count files",1,0,1);  
87   fOutputList->Add(fhCount);
88
89 }
90
91 //______________________________________________________________________________
92 void AliAnaScale::Exec(Option_t *) 
93 {
94   // Do the Scaling
95
96   if(fDebug > 0 ) printf(">>>>> Scaling factor %e, do Sumw2 %d <<<<< \n",fScale,fSumw2) ;
97   
98   const Int_t buffersize = 255;
99   char name[buffersize] ; 
100
101   TIter next(fInputList) ;      
102   TObject * h ; 
103   while ( (h = next()) ) { 
104     if(h){
105       if ( !strncmp(h->ClassName(),"TH",2) ) {
106       snprintf(name, buffersize, "%sScaled", h->GetName()) ; 
107       TH1 * hout = dynamic_cast<TH1*> (h->Clone(name)) ; 
108       if(hout){
109         if(fSumw2) hout->Sumw2();
110         hout->Scale(fScale) ;  
111         fOutputList->Add(hout) ;
112         }// casting not null
113       } 
114       else  fOutputList->Add(h) ; 
115     }
116   }
117   // number of files
118
119   //File scaled, needed for file merging on grid
120   fhCount->Fill(0);
121  
122   PostData(0, fOutputList);
123 }
124
125
126 //______________________________________________________________________________
127 void AliAnaScale::Init()
128 {
129   // Intialisation of parameters
130   if(fDebug > 0 )printf("No initialization in scale class \n") ;
131
132 }
133
134 //______________________________________________________________________________
135 //void AliAnaScale::Terminate(Option_t *)
136 //{
137 //  // Processing when the event loop is ended
138 //  
139 //
140 //}