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     fDebugLevel(0)
  67 {
  68   //
  69   // Constructor
  70   //
  71 }
  72
  73 AliTPCcalibBase::AliTPCcalibBase(const AliTPCcalibBase&calib):
  74   TNamed(calib),
  75   fDebugStreamer(0),
  76   fStreamLevel(calib.fStreamLevel),
  77   fRun(0),                  //!  current Run number
  78   fEvent(0),                //!  current Event number
  79   fTime(0),                 //!  current Time
  80   fTrigger(0),              //! current trigger type
  81   fMagF(0),                 //! current magnetic field
  82   fTriggerMaskReject(calib.fTriggerMaskReject),   //trigger mask - reject trigger
  83   fTriggerMaskAccept(calib.fTriggerMaskAccept),   //trigger mask - accept trigger
  84   fHasLaser(calib.fHasLaser),                    //flag the laser is overlayed with given event
  85   fRejectLaser(calib.fRejectLaser),                 //flag- reject laser
  86   fDebugLevel(calib.fDebugLevel)
  87 {
  88   //
  89   // copy constructor
  90   //
  91 }
  92
  93 AliTPCcalibBase &AliTPCcalibBase::operator=(const AliTPCcalibBase&calib){
  94   //
  95   //
  96   //
  97   ((TNamed *)this)->operator=(calib);
  98   fDebugStreamer=0;
  99   fStreamLevel=calib.fStreamLevel;
 100   fDebugLevel=calib.fDebugLevel;
 101   return *this;
 102 }
 103
 104
 105 AliTPCcalibBase::~AliTPCcalibBase() {
 106   //
 107   // destructor
 108   //
 109   if (fDebugLevel>0) printf("AliTPCcalibBase::~AliTPCcalibBase\n");
 110   if (fDebugStreamer) delete fDebugStreamer;
 111   fDebugStreamer=0;
 112 }
 113
 114 void  AliTPCcalibBase::Terminate(){
 115   //
 116   //
 117   //
 118   if (fDebugLevel>0) printf("AliTPCcalibBase::Terminate\n");
 119   if (fDebugStreamer) delete fDebugStreamer;
 120   fDebugStreamer = 0;
 121   return;
 122 }
 123
 124 TTreeSRedirector *AliTPCcalibBase::GetDebugStreamer(){
 125   //
 126   // Get Debug streamer
 127   // In case debug streamer not yet initialized and StreamLevel>0 create new one
 128   //
 129   if (fStreamLevel==0) return 0;
 130   if (fDebugStreamer) return fDebugStreamer;
 131   TString dsName;
 132   dsName=GetName();
 133   dsName+="Debug.root";
 134   dsName.ReplaceAll(" ","");
 135   fDebugStreamer = new TTreeSRedirector(dsName.Data());
 136   return fDebugStreamer;
 137 }
 138
 139
 140 void    AliTPCcalibBase::UpdateEventInfo(AliESDEvent * event){
 141   //
 142   //
 143   //
 144   fRun     = event->GetRunNumber();
 145   fEvent   = event->GetEventNumberInFile();
 146   fTime    = event->GetTimeStamp();
 147   fTrigger = event->GetTriggerMask();
 148   fMagF    = event->GetMagneticField();
 149   fHasLaser = HasLaser(event);
 150 }
 151
 152
 153 Bool_t AliTPCcalibBase::HasLaser(AliESDEvent *event){
 154   //
 155   //
 156   //
 157   // Thresholds more than 8 tracks with small dip angle
 158
 159   const Int_t kMinLaserTracks = 8;
 160   const Float_t kThrLaser       = 0.3;
 161   const Float_t kLaserTgl       = 0.01;
 162
 163   Int_t ntracks = event->GetNumberOfTracks();
 164   if (ntracks<kMinLaserTracks) return kFALSE;
 165   Float_t nlaser=0;
 166   Float_t nall=0;
 167   for (Int_t i=0;i<ntracks;++i) {
 168     AliESDtrack *track=event->GetTrack(i);
 169     if (!track) continue;
 170     if (track->GetTPCNcls()<=0) continue;
 171     nall++;
 172     if (TMath::Abs(track->GetTgl())<kLaserTgl) nlaser++;
 173   }
 174   if (nlaser>kMinLaserTracks) return kTRUE;
 175   if (nall>0 && nlaser/nall>kThrLaser) return kTRUE;
 176   return kFALSE;
 177 }
 178
 179
 180
 181 Bool_t AliTPCcalibBase::AcceptTrigger(){
 182   //
 183   // Apply trigger mask - Don't do calibration for non proper triggers
 184   //
 185   if (fTriggerMaskReject==(Int_t)fTrigger) return kFALSE;
 186   if (fTriggerMaskAccept>0 && fTriggerMaskAccept!=(Int_t)fTrigger) return kFALSE;
 187   if (fHasLaser && fRejectLaser) return kFALSE;
 188   return kTRUE;
 189 }
 190
 191
 192 void AliTPCcalibBase::RegisterDebugOutput(const char *path){
 193   //
 194   // store  - copy debug output to the destination position
 195   // currently ONLY for local copy
 196   if (fDebugLevel>0) printf("AliTPCcalibBase::RegisterDebugOutput(%s)\n",path);
 197   if (fStreamLevel==0) return;
 198   TString dsName;
 199   dsName=GetName();
 200   dsName+="Debug.root";
 201   dsName.ReplaceAll(" ","");
 202   TString dsName2=path;
 203   gSystem->MakeDirectory(dsName2.Data());
 204   dsName2+=gSystem->HostName();
 205   gSystem->MakeDirectory(dsName2.Data());
 206   dsName2+="/";
 207   TTimeStamp s;
 208   dsName2+=Int_t(s.GetNanoSec());
 209   dsName2+="/";
 210   gSystem->MakeDirectory(dsName2.Data());
 211   dsName2+=dsName;
 212   AliInfo(Form("copy %s\t%s\n",dsName.Data(),dsName2.Data()));
 213   printf("copy %s\t%s\n",dsName.Data(),dsName2.Data());
 214   TFile::Cp(dsName.Data(),dsName2.Data());
 215 }