VZERO/AliVZEROTriggerMask.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 #include <Riostream.h>
  18 #include <TGeoManager.h>
  19 #include <TGeoMatrix.h>
  20 #include <TGeoPhysicalNode.h>
  21 #include <AliGeomManager.h>
  22 #include "AliCDBManager.h"
  23 #include "AliCDBStorage.h"
  24 #include "AliCDBEntry.h"
  25 #include "AliLog.h"
  26 #include "AliVZEROTriggerMask.h"
  27 #include "AliVZEROdigit.h"
  28
  29 //______________________________________________________________________
  30 ClassImp(AliVZEROTriggerMask)
  31
  32 //______________________________________________________________________
  33
  34 AliVZEROTriggerMask::AliVZEROTriggerMask()
  35   :TObject(),
  36    fAdcThresHold(0.0),
  37    fTimeWindowWidthBBA(50.0),
  38    fTimeWindowWidthBGA(20.0),
  39    fTimeWindowWidthBBC(50.0),
  40    fTimeWindowWidthBGC(20.0),
  41    fBBtriggerV0A(0),
  42    fBGtriggerV0A(0),
  43    fBBtriggerV0C(0),
  44    fBGtriggerV0C(0)
  45
  46 {
  47    SetAdcThreshold();
  48 }
  49
  50 //________________________________________________________________________________
  51 Double_t AliVZEROTriggerMask::GetZPosition(const char* symname){
  52 // Get the global z coordinate of the given V0 alignable volume
  53 //
  54   Double_t *tr;
  55   TGeoPNEntry *pne = gGeoManager->GetAlignableEntry(symname);
  56   if (!pne) return 0;
  57
  58
  59   TGeoPhysicalNode *pnode = pne->GetPhysicalNode();
  60   if(pnode){
  61           TGeoHMatrix* hm = pnode->GetMatrix();
  62            tr = hm->GetTranslation();
  63   }else{
  64           const char* path = pne->GetTitle();
  65           if(!gGeoManager->cd(path)){
  66                   AliErrorClass(Form("Volume path %s not valid!",path));
  67                   return 0;
  68           }
  69          tr = gGeoManager->GetCurrentMatrix()->GetTranslation();
  70   }
  71   return tr[2];
  72
  73 }
  74
  75 //________________________________________________________________________________
  76
  77
  78 void AliVZEROTriggerMask::FillMasks(TTree* vzeroDigitsTree,
  79                                 TClonesArray* vzeroDigits)
  80 {
  81   const Double_t LightSpeed = 2.9979245800; // cm/100 ps
  82   Float_t Z_V0A = TMath::Abs(GetZPosition("VZERO/V0A"));
  83   Float_t Z_V0C = TMath::Abs(GetZPosition("VZERO/V0C"));
  84
  85   // distance in time units from nominal vertex to V0
  86   Float_t v0a_dist = Z_V0A/LightSpeed; // 100 of picoseconds
  87   Float_t v0c_dist = Z_V0C/LightSpeed; // 100 of picoseconds
  88   Float_t bunch_separation = 1000.0; // 100 of picoseconds
  89
  90   // mask
  91   UInt_t one=1;
  92
  93   // loop over vzero entries
  94   Int_t nEntries = (Int_t)vzeroDigitsTree->GetEntries();
  95   for (Int_t e=0; e<nEntries; e++) {
  96     vzeroDigitsTree->GetEvent(e);
  97
  98     Int_t nDigits = vzeroDigits->GetEntriesFast();
  99
 100     for (Int_t d=0; d<nDigits; d++) {
 101       //      vzeroDigitsTree->GetEvent(d);
 102       AliVZEROdigit* digit = (AliVZEROdigit*)vzeroDigits->At(d);
 103
 104       Int_t   PMNumber   = digit->PMNumber();
 105       Float_t adc        = digit->ADC();
 106       Float_t tdc        = digit->Time(); // in 100 of picoseconds
 107
 108       if (adc>fAdcThresHold) {
 109         if (PMNumber<32) { // in V0C
 110           if (tdc>(v0c_dist-fTimeWindowWidthBBC/2.0) &&
 111               tdc<(v0c_dist+fTimeWindowWidthBBC/2.0))
 112             fBBtriggerV0C+=(one<<PMNumber);
 113           if (tdc>(bunch_separation-v0c_dist-fTimeWindowWidthBGC/2.0) &&
 114               tdc<(bunch_separation-v0c_dist+fTimeWindowWidthBGC/2.0))
 115            fBGtriggerV0C+=(one<<PMNumber);
 116         }
 117         if (PMNumber>31) { // in V0A
 118           Int_t shift = PMNumber-32;
 119           if (tdc>(v0a_dist-fTimeWindowWidthBBA/2.0) &&
 120               tdc<(v0a_dist+fTimeWindowWidthBBA/2.0))
 121             fBBtriggerV0A+=(one<<shift);
 122           if (tdc>(bunch_separation-v0a_dist-fTimeWindowWidthBGA/2.0) &&
 123               tdc<(bunch_separation-v0a_dist+fTimeWindowWidthBGA/2.0))
 124             fBGtriggerV0A+=(one<<shift);
 125         }
 126       }
 127     } // end of loop over digits
 128   } // end of loop over events in digits tree
 129 }