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   fDebugLevel(calib.fDebugLevel)
 109 {
 110   //
 111   // copy constructor
 112   //
 113 }
 114
 115 AliTPCcalibBase &AliTPCcalibBase::operator=(const AliTPCcalibBase&calib){
 116   //
 117   //
 118   //
 119   ((TNamed *)this)->operator=(calib);
 120   fDebugStreamer=0;
 121   fStreamLevel=calib.fStreamLevel;
 122   fDebugLevel=calib.fDebugLevel;
 123   return *this;
 124 }
 125
 126
 127 AliTPCcalibBase::~AliTPCcalibBase() {
 128   //
 129   // destructor
 130   //
 131   if (fDebugLevel>0) printf("AliTPCcalibBase::~AliTPCcalibBase\n");
 132   if (fDebugStreamer) delete fDebugStreamer;
 133   fDebugStreamer=0;
 134 }
 135
 136 void  AliTPCcalibBase::Terminate(){
 137   //
 138   //
 139   //
 140   if (fDebugLevel>0) printf("AliTPCcalibBase::Terminate\n");
 141   if (fDebugStreamer) delete fDebugStreamer;
 142   fDebugStreamer = 0;
 143   return;
 144 }
 145
 146 TTreeSRedirector *AliTPCcalibBase::GetDebugStreamer(){
 147   //
 148   // Get Debug streamer
 149   // In case debug streamer not yet initialized and StreamLevel>0 create new one
 150   //
 151   if (fStreamLevel==0) return 0;
 152   if (fDebugStreamer) return fDebugStreamer;
 153   TString dsName;
 154   dsName=GetName();
 155   dsName+="Debug.root";
 156   dsName.ReplaceAll(" ","");
 157   fDebugStreamer = new TTreeSRedirector(dsName.Data());
 158   return fDebugStreamer;
 159 }
 160
 161
 162 void    AliTPCcalibBase::UpdateEventInfo(AliESDEvent * event){
 163   //
 164   //
 165   //
 166   fRun     = event->GetRunNumber();
 167   fEvent   = event->GetEventNumberInFile();
 168   fTime    = event->GetTimeStamp();
 169   fTrigger = event->GetTriggerMask();
 170   fMagF    = event->GetMagneticField();
 171   fTriggerClass = event->GetFiredTriggerClasses().Data();
 172   fHasLaser = HasLaser(event);
 173
 174 }
 175
 176
 177 Bool_t AliTPCcalibBase::HasLaser(AliESDEvent *event){
 178   //
 179   //
 180   //
 181   // Thresholds more than 8 tracks with small dip angle
 182
 183   const Int_t kMinLaserTracks = 8;
 184   const Float_t kThrLaser       = 0.3;
 185   const Float_t kLaserTgl       = 0.01;
 186
 187   Int_t ntracks = event->GetNumberOfTracks();
 188   if (ntracks<kMinLaserTracks) return kFALSE;
 189   Float_t nlaser=0;
 190   Float_t nall=0;
 191   for (Int_t i=0;i<ntracks;++i) {
 192     AliESDtrack *track=event->GetTrack(i);
 193     if (!track) continue;
 194     if (track->GetTPCNcls()<=0) continue;
 195     nall++;
 196     if (TMath::Abs(track->GetTgl())<kLaserTgl) nlaser++;
 197   }
 198   if (nlaser>kMinLaserTracks) return kTRUE;
 199   if (nall>0 && nlaser/nall>kThrLaser) return kTRUE;
 200   return kFALSE;
 201 }
 202
 203
 204
 205 Bool_t AliTPCcalibBase::AcceptTrigger(){
 206   //
 207   // Apply trigger mask - Don't do calibration for non proper triggers
 208   //
 209   if (fTriggerMaskReject==(Int_t)fTrigger) return kFALSE;
 210   if (fTriggerMaskAccept>0 && fTriggerMaskAccept!=(Int_t)fTrigger) return kFALSE;
 211   if (fHasLaser && fRejectLaser) return kFALSE;
 212   return kTRUE;
 213 }
 214
 215
 216 void AliTPCcalibBase::RegisterDebugOutput(const char *path){
 217   //
 218   // store  - copy debug output to the destination position
 219   // currently ONLY for local copy
 220   if (fDebugLevel>0) printf("AliTPCcalibBase::RegisterDebugOutput(%s)\n",path);
 221   if (fStreamLevel==0) return;
 222   TString dsName;
 223   dsName=GetName();
 224   dsName+="Debug.root";
 225   dsName.ReplaceAll(" ","");
 226   TString dsName2=path;
 227   gSystem->MakeDirectory(dsName2.Data());
 228   dsName2+=gSystem->HostName();
 229   gSystem->MakeDirectory(dsName2.Data());
 230   dsName2+="/";
 231   TTimeStamp s;
 232   dsName2+=Int_t(s.GetNanoSec());
 233   dsName2+="/";
 234   gSystem->MakeDirectory(dsName2.Data());
 235   dsName2+=dsName;
 236   AliInfo(Form("copy %s\t%s\n",dsName.Data(),dsName2.Data()));
 237   printf("copy %s\t%s\n",dsName.Data(),dsName2.Data());
 238   TFile::Cp(dsName.Data(),dsName2.Data());
 239 }