Warning removal
[u/mrichter/AliRoot.git] / TPC / AliTPCcalibBase.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
17 ///////////////////////////////////////////////////////////////////////////////
18 //                                                                           //
19 //  Base class for the calibration components using 
20 //  as input TPCseeds and ESDs
21 //  Event loop outside of the component
22 //
23 //
24 // Base functionality to be implemeneted by component 
25 /* 
26    //In some cases only one of this function to be implemented
27    virtual void     Process(AliESDEvent *event)
28    virtual void     Process(AliTPCseed *track)
29    //
30    virtual Long64_t Merge(TCollection *li);
31    virtual void     Analyze()
32    void             Terminate();
33 */
34 // Functionality provided by base class for Algorith debuging:
35 //  TTreeSRedirector * cstream =  GetDebugStreamer() - get debug streamer which can be use for numerical debugging
36 //                      
37
38
39
40 //  marian.ivanov@cern.ch
41 // 
42 #include "AliTPCcalibBase.h"
43 #include "TSystem.h"
44 #include "TFile.h"
45 #include "TTreeStream.h"
46 #include "AliLog.h"
47 #include "TTimeStamp.h"
48 #include "AliESDEvent.h"
49
50
51 ClassImp(AliTPCcalibBase)
52
53 AliTPCcalibBase::AliTPCcalibBase():
54     TNamed(),
55     fDebugStreamer(0),
56     fStreamLevel(0),   
57     fRun(0),                  //!  current Run number
58     fEvent(0),                //!  current Event number
59     fTime(0),                 //!  current Time
60     fTrigger(0),              //! current trigger type
61     fMagF(0),                 //! current magnetic field
62     fTriggerMaskReject(-1),   //trigger mask - reject trigger
63     fTriggerMaskAccept(-1),   //trigger mask - accept trigger
64     fHasLaser(kFALSE),                    //flag the laser is overlayed with given event 
65     fRejectLaser(kTRUE),                 //flag- reject laser
66     fTriggerClass(),
67     fDebugLevel(0)
68 {
69   //
70   // Constructor
71   //
72 }
73
74 AliTPCcalibBase::AliTPCcalibBase(const char * name, const char * title):
75   TNamed(name,title),
76   fDebugStreamer(0),
77   fStreamLevel(0),   
78   fRun(0),                  //!  current Run number
79   fEvent(0),                //!  current Event number
80   fTime(0),                 //!  current Time
81   fTrigger(0),              //! current trigger type
82   fMagF(0),                 //! current magnetic field
83   fTriggerMaskReject(-1),   //trigger mask - reject trigger
84   fTriggerMaskAccept(-1),   //trigger mask - accept trigger
85   fHasLaser(kFALSE),                    //flag the laser is overlayed with given event 
86   fRejectLaser(kTRUE),                 //flag- reject laser
87   fTriggerClass(),
88   fDebugLevel(0)
89 {
90   //
91   // Constructor
92   //
93 }
94
95 AliTPCcalibBase::AliTPCcalibBase(const AliTPCcalibBase&calib):
96   TNamed(calib),
97   fDebugStreamer(0),
98   fStreamLevel(calib.fStreamLevel),
99   fRun(0),                  //!  current Run number
100   fEvent(0),                //!  current Event number
101   fTime(0),                 //!  current Time
102   fTrigger(0),              //! current trigger type
103   fMagF(0),                 //! current magnetic field
104   fTriggerMaskReject(calib.fTriggerMaskReject),   //trigger mask - reject trigger
105   fTriggerMaskAccept(calib.fTriggerMaskAccept),   //trigger mask - accept trigger
106   fHasLaser(calib.fHasLaser),                    //flag the laser is overlayed with given event
107   fRejectLaser(calib.fRejectLaser),                 //flag- reject laser
108   fTriggerClass(calib.fTriggerClass),
109   fDebugLevel(calib.fDebugLevel)
110 {
111   //
112   // copy constructor
113   //
114 }
115
116 AliTPCcalibBase &AliTPCcalibBase::operator=(const AliTPCcalibBase&calib){
117   //
118   //
119   //
120   ((TNamed *)this)->operator=(calib);
121   fDebugStreamer=0;
122   fStreamLevel=calib.fStreamLevel;
123   fDebugLevel=calib.fDebugLevel;
124   return *this;
125 }
126
127
128 AliTPCcalibBase::~AliTPCcalibBase() {
129   //
130   // destructor
131   //
132   if (fDebugLevel>0) printf("AliTPCcalibBase::~AliTPCcalibBase\n");
133   if (fDebugStreamer) delete fDebugStreamer;
134   fDebugStreamer=0;
135 }
136
137 void  AliTPCcalibBase::Terminate(){
138   //
139   //
140   //
141   if (fDebugLevel>0) printf("AliTPCcalibBase::Terminate\n");
142   if (fDebugStreamer) delete fDebugStreamer;
143   fDebugStreamer = 0;
144   return;
145 }
146
147 TTreeSRedirector *AliTPCcalibBase::GetDebugStreamer(){
148   //
149   // Get Debug streamer
150   // In case debug streamer not yet initialized and StreamLevel>0 create new one
151   //
152   if (fStreamLevel==0) return 0;
153   if (fDebugStreamer) return fDebugStreamer;
154   TString dsName;
155   dsName=GetName();
156   dsName+="Debug.root";
157   dsName.ReplaceAll(" ",""); 
158   fDebugStreamer = new TTreeSRedirector(dsName.Data());
159   return fDebugStreamer;
160 }
161
162
163 void    AliTPCcalibBase::UpdateEventInfo(AliESDEvent * event){
164   //
165   //
166   //
167   fRun     = event->GetRunNumber();
168   fEvent   = event->GetEventNumberInFile();
169   fTime    = event->GetTimeStamp();
170   fTrigger = event->GetTriggerMask();
171   fMagF    = event->GetMagneticField();
172   fTriggerClass = event->GetFiredTriggerClasses().Data();
173   fHasLaser = HasLaser(event); 
174   
175 }
176
177
178 Bool_t AliTPCcalibBase::HasLaser(AliESDEvent *event){
179   //
180   //
181   //
182   // Thresholds more than 8 tracks with small dip angle
183   
184   const Int_t kMinLaserTracks = 8;
185   const Float_t kThrLaser       = 0.3;
186   const Float_t kLaserTgl       = 0.01;
187
188   Int_t ntracks = event->GetNumberOfTracks();
189   if (ntracks<kMinLaserTracks) return kFALSE;
190   Float_t nlaser=0;
191   Float_t nall=0;
192   for (Int_t i=0;i<ntracks;++i) {
193     AliESDtrack *track=event->GetTrack(i);
194     if (!track) continue;
195     if (track->GetTPCNcls()<=0) continue; 
196     nall++;
197     if (TMath::Abs(track->GetTgl())<kLaserTgl) nlaser++;
198   }
199   if (nlaser>kMinLaserTracks) return kTRUE;
200   if (nall>0 && nlaser/nall>kThrLaser) return kTRUE;
201   return kFALSE;
202 }
203
204
205
206 Bool_t AliTPCcalibBase::AcceptTrigger(){
207   //
208   // Apply trigger mask - Don't do calibration for non proper triggers
209   // 
210   if (fTriggerMaskReject==(Int_t)fTrigger) return kFALSE;
211   if (fTriggerMaskAccept>0 && fTriggerMaskAccept!=(Int_t)fTrigger) return kFALSE;
212   if (fHasLaser && fRejectLaser) return kFALSE;
213   return kTRUE;
214 }
215
216
217 void AliTPCcalibBase::RegisterDebugOutput(const char *path){
218   //
219   // store  - copy debug output to the destination position
220   // currently ONLY for local copy
221   if (fDebugLevel>0) printf("AliTPCcalibBase::RegisterDebugOutput(%s)\n",path);
222   if (fStreamLevel==0) return;
223   TString dsName;
224   dsName=GetName();
225   dsName+="Debug.root";
226   dsName.ReplaceAll(" ",""); 
227   TString dsName2=path;
228   gSystem->MakeDirectory(dsName2.Data());
229   dsName2+=gSystem->HostName();
230   gSystem->MakeDirectory(dsName2.Data());
231   dsName2+="/";
232   TTimeStamp s;
233   dsName2+=Int_t(s.GetNanoSec());
234   dsName2+="/";
235   gSystem->MakeDirectory(dsName2.Data());
236   dsName2+=dsName;
237   AliInfo(Form("copy %s\t%s\n",dsName.Data(),dsName2.Data()));
238   printf("copy %s\t%s\n",dsName.Data(),dsName2.Data());
239   TFile::Cp(dsName.Data(),dsName2.Data());
240 }