ZDC/AliZDCTrigger.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 //      Trigger class for ZDC
  19 //
  20 // ****************************************************************
  21
  22 #include "AliLog.h"
  23 #include "AliRun.h"
  24 #include "AliLoader.h"
  25 #include "AliRunLoader.h"
  26 #include "AliTriggerInput.h"
  27
  28 #include "AliZDC.h"
  29 #include "AliZDCDigit.h"
  30 #include "AliZDCTrigger.h"
  31
  32 //________________________________________________________________
  33 ClassImp(AliZDCTrigger)
  34
  35 //________________________________________________________________
  36 AliZDCTrigger::AliZDCTrigger() : AliTriggerDetector()
  37 {
  38    // Constructor
  39    SetName("ZDC");
  40    CreateInputs();
  41    //
  42    SetZDCLeftMinCut(0);
  43    SetZDCRightMinCut(0);
  44    SetZEMMinCut(0);
  45    SetZDCLeftEMDCuts(0,0);
  46    SetZDCRightEMDCuts(0,0);
  47    SetZDCLeftMBCut(0);
  48    SetZDCRightMBCut(0);
  49    SetZDCLeftCentrCut(0);
  50    SetZDCRightCentrCut(0);
  51    SetZDCLeftSemiCentrCut(0);
  52    SetZDCRightSemiCentrCut(0);
  53    SetZEMCentrCut(0);
  54
  55 }
  56
  57 //________________________________________________________________
  58 void AliZDCTrigger::CreateInputs()
  59 {
  60    // Trigger inputs
  61
  62    // Do not create inputs again!!
  63    if( fInputs.GetEntriesFast() > 0 ) return;
  64
  65    fInputs.AddLast(new AliTriggerInput("ZDC_1_L1",   "ZDC Minimum Bias", 0x01));
  66    fInputs.AddLast(new AliTriggerInput("ZDC_2_L1",   "ZDC Semi-central", 0x02));
  67    fInputs.AddLast(new AliTriggerInput("ZDC_3_L1",   "ZDC Central", 0x04));
  68    fInputs.AddLast(new AliTriggerInput("ZDC_EMD_L1", "ZDC EMD events", 0x08));
  69 }
  70
  71 //________________________________________________________________
  72 void AliZDCTrigger::Trigger()
  73 {
  74
  75    // Trigger selection
  76    //
  77    AliRunLoader *runLoader = gAlice->GetRunLoader();
  78
  79    AliLoader *aZDCLoader = runLoader->GetLoader("ZDCLoader");
  80    aZDCLoader->LoadDigits("READ");
  81    AliZDCDigit digit;
  82    AliZDCDigit* pdigit = &digit;
  83    TTree* tD = aZDCLoader->TreeD();
  84    if (!tD) cerr<<"AliZDCTrigger: digits tree not found\n";
  85    tD->SetBranchAddress("ZDC", &pdigit);
  86    //
  87    Float_t signalZNLeft[2], signalZPLeft[2], signalZDCLeftSum[2];
  88    Float_t signalZNRight[2], signalZPRight[2], signalZDCRightSum[2];
  89    Float_t signalZEMSum[2];
  90    for(Int_t iDigit=0; iDigit<tD->GetEntries(); iDigit++){
  91       tD->GetEntry(iDigit);
  92       //
  93       // *** ZDC LEFT
  94       if(digit.GetSector(0)==1)
  95          for(Int_t i=0; i<2; i++){ //0=high range; 1=low range
  96             signalZNLeft[i] += digit.GetADCValue(i);
  97             signalZDCLeftSum[i] += digit.GetADCValue(i);
  98           }
  99       else if(digit.GetSector(0)==2)
 100          for(Int_t i=0; i<2; i++){
 101             signalZPLeft[i] += digit.GetADCValue(i);
 102             signalZDCLeftSum[i] += digit.GetADCValue(i);
 103           }
 104       else if(digit.GetSector(0)==3)
 105          for(Int_t i=0; i<2; i++) signalZEMSum[i] += digit.GetADCValue(i);
 106       // *** ZDC RIGHT
 107       else if(digit.GetSector(0)==4)
 108          for(Int_t i=0; i<2; i++){ //0=high range; 1=low range
 109             signalZNRight[i] += digit.GetADCValue(i);
 110             signalZDCRightSum[i] += digit.GetADCValue(i);
 111           }
 112       else if(digit.GetSector(0)==5)
 113          for(Int_t i=0; i<2; i++){
 114             signalZPRight[i] += digit.GetADCValue(i);
 115             signalZDCRightSum[i] += digit.GetADCValue(i);
 116           }
 117    }
 118    // *******************************************************************
 119    if(signalZDCLeftSum[1]>fZDCLeftMBCut && signalZDCRightSum[1]>fZDCRightMBCut)
 120        // *** ZDC minimum bias trigger
 121        SetInput("ZDC_1_L1");
 122    // *******************************************************************
 123    if(signalZDCLeftSum[1]>fZDCLeftCentrCut && signalZDCLeftSum[1]<fZDCLeftSemiCentrCut &&
 124       signalZDCRightSum[1]>fZDCRightCentrCut && signalZDCRightSum[1]<fZDCRightSemiCentrCut
 125       && signalZEMSum[1]>fZEMCentrCut)
 126        // *** ZDC semi-central (10-40%)
 127        SetInput("ZDC_2_L1");
 128    // *******************************************************************
 129    if(signalZDCLeftSum[1]>fZDCLeftMinCut && signalZDCLeftSum[1]<fZDCLeftCentrCut &&
 130       signalZDCRightSum[1]>fZDCRightMinCut && signalZDCRightSum[1]<fZDCRightCentrCut &&
 131       signalZEMSum[1]>fZEMCentrCut)
 132        // *** ZDC central (0-10%)
 133        SetInput("ZDC_3_L1");
 134    // *******************************************************************
 135    if(signalZNLeft[0]>fZDCLeftEMDCuts[0] && signalZNLeft[0]<fZDCLeftEMDCuts[1] &&
 136       signalZNRight[0]>fZDCRightEMDCuts[0] && signalZNRight[0]<fZDCRightEMDCuts[1] &&
 137       signalZEMSum[1]<fZEMMinCut){ // *** 1n EMD trigger
 138         SetInput("ZDC_EMD_L1");
 139    }
 140
 141 }
 142
 143 //________________________________________________________________
 144 void AliZDCTrigger::SetZDCLeftMinCut(Float_t ZDCLeftMinCut)
 145 {
 146   // Set default cut values for ZDC trigger
 147   //
 148   if(ZDCLeftMinCut)  fZDCLeftMinCut = ZDCLeftMinCut;
 149   else  fZDCLeftMinCut = 800.;
 150 }
 151 //________________________________________________________________
 152 void AliZDCTrigger::SetZDCRightMinCut(Float_t ZDCRightMinCut)
 153 {
 154   // Set default cut values for ZDC trigger
 155   //
 156   if(ZDCRightMinCut)  fZDCRightMinCut = ZDCRightMinCut;
 157   else  fZDCRightMinCut = 800.;
 158 }
 159
 160 //________________________________________________________________
 161 void AliZDCTrigger::SetZEMMinCut(Float_t ZEMMinCut)
 162 {
 163   // Set default cut values for ZDC trigger
 164   //
 165   if(ZEMMinCut)  fZEMMinCut = ZEMMinCut;
 166   else  fZEMMinCut = 80.;
 167 }
 168 //________________________________________________________________
 169 void AliZDCTrigger::SetZDCLeftEMDCuts(Float_t* ZDCLeftEMDCuts)
 170 {
 171   // Set default cut values for ZDC trigger
 172   //
 173   if(ZDCLeftEMDCuts) for(int j=0; j<2; j++) fZDCLeftEMDCuts[j] = ZDCLeftEMDCuts[j];
 174   else{
 175     fZDCLeftEMDCuts[0] = 600.;
 176     fZDCLeftEMDCuts[1] = 1000.;
 177   }
 178 }
 179 //________________________________________________________________
 180 void AliZDCTrigger::SetZDCLeftEMDCuts(Float_t ZDCLeftEMDCutInf,
 181         Float_t ZDCLeftEMDCutSup)
 182 {
 183   // Set default cut values for ZDC trigger
 184   //
 185   if(ZDCLeftEMDCutInf && ZDCLeftEMDCutSup){
 186     fZDCLeftEMDCuts[0]=ZDCLeftEMDCutInf;
 187     fZDCLeftEMDCuts[1]=ZDCLeftEMDCutSup;
 188   }
 189   else{
 190     fZDCLeftEMDCuts[0] = 600.;
 191     fZDCLeftEMDCuts[1] = 1000.;
 192   }
 193 }
 194 //________________________________________________________________
 195 void AliZDCTrigger::SetZDCRightEMDCuts(Float_t* ZDCRightEMDCuts)
 196 {
 197   // Set default cut values for ZDC trigger
 198   //
 199   if(ZDCRightEMDCuts) for(int j=0; j<2; j++) fZDCRightEMDCuts[j] = ZDCRightEMDCuts[j];
 200   else{
 201     fZDCRightEMDCuts[0] = 600.;
 202     fZDCRightEMDCuts[1] = 1000.;
 203   }
 204 }
 205 //________________________________________________________________
 206 void AliZDCTrigger::SetZDCRightEMDCuts(Float_t ZDCRightEMDCutInf,
 207         Float_t ZDCRightEMDCutSup)
 208 {
 209   // Set default cut values for ZDC trigger
 210   //
 211   if(ZDCRightEMDCutInf && ZDCRightEMDCutSup){
 212     fZDCRightEMDCuts[0]=ZDCRightEMDCutInf;
 213     fZDCRightEMDCuts[1]=ZDCRightEMDCutSup;
 214   }
 215   else{
 216     fZDCRightEMDCuts[0] = 600.;
 217     fZDCRightEMDCuts[1] = 1000.;
 218   }
 219 }
 220 //________________________________________________________________
 221 void AliZDCTrigger::SetZDCLeftMBCut(Float_t ZDCLeftMBCut)
 222 {
 223   // Set default cut values for ZDC trigger
 224   //
 225   if(ZDCLeftMBCut) fZDCLeftMBCut = ZDCLeftMBCut;
 226   else fZDCLeftMBCut = 800.;
 227 }
 228 //________________________________________________________________
 229 void AliZDCTrigger::SetZDCRightMBCut(Float_t ZDCRightMBCut)
 230 {
 231   // Set default cut values for ZDC trigger
 232   //
 233   if(ZDCRightMBCut) fZDCRightMBCut = ZDCRightMBCut;
 234   else fZDCRightMBCut = 800.;
 235 }
 236 //________________________________________________________________
 237 void AliZDCTrigger::SetZDCLeftCentrCut(Float_t ZDCLeftCentrCut)
 238 {
 239   // Set default cut values for ZDC trigger
 240   //
 241   if(ZDCLeftCentrCut) fZDCLeftCentrCut = ZDCLeftCentrCut;
 242   else fZDCLeftCentrCut = 10000.;
 243 }
 244 //________________________________________________________________
 245 void AliZDCTrigger::SetZDCRightCentrCut(Float_t ZDCRightCentrCut)
 246 {
 247   // Set default cut values for ZDC trigger
 248   //
 249   if(ZDCRightCentrCut) fZDCRightCentrCut = ZDCRightCentrCut;
 250   else fZDCRightCentrCut = 10000.;
 251 }
 252 //________________________________________________________________
 253 void AliZDCTrigger::SetZDCLeftSemiCentrCut(Float_t ZDCLeftSemiCentrCut)
 254 {
 255   // Set default cut values for ZDC trigger
 256   //
 257   if(ZDCLeftSemiCentrCut) fZDCLeftSemiCentrCut = ZDCLeftSemiCentrCut;
 258   else fZDCLeftSemiCentrCut = 18500.;
 259 }
 260 //________________________________________________________________
 261 void AliZDCTrigger::SetZDCRightSemiCentrCut(Float_t ZDCRightSemiCentrCut)
 262 {
 263   // Set default cut values for ZDC trigger
 264   //
 265   if(ZDCRightSemiCentrCut) fZDCRightSemiCentrCut = ZDCRightSemiCentrCut;
 266   else fZDCRightSemiCentrCut = 18500.;
 267 }
 268 //________________________________________________________________
 269 void AliZDCTrigger::SetZEMCentrCut(Float_t ZEMCentrCut)
 270 {
 271   // Set default cut values for ZDC trigger
 272   //
 273   if(ZEMCentrCut) fZEMCentrCut = ZEMCentrCut;
 274   else fZEMCentrCut = 210.;
 275 }