Trigger class for ZDC

diff --git a/ZDC/AliZDCTrigger.cxx b/ZDC/AliZDCTrigger.cxx
new file mode 100644 (file)
index 0000000..00b96fc
--- /dev/null
+++ b/ZDC/AliZDCTrigger.cxx
@@ -0,0 +1,200 @@
+/**************************************************************************
+ * Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved. *
+ *                                                                         *
+ * Author: The ALICE Off-line Project.                                    *
+ * Contributors are mentioned in the code where appropriate.              *
+ *                                                                        *
+ * Permission to use, copy, modify and distribute this software and its   *
+ * documentation strictly for non-commercial purposes is hereby granted   *
+ * without fee, provided that the above copyright notice appears in all   *
+ * copies and that both the copyright notice and this permission notice   *
+ * appear in the supporting documentation. The authors make no claims     *
+ * about the suitability of this software for any purpose. It is          *
+ * provided "as is" without express or implied warranty.                  *
+ **************************************************************************/
+
+#include "AliLog.h"
+#include "AliRun.h"
+#include "AliLoader.h"
+#include "AliRunLoader.h"
+#include "AliTriggerInput.h"
+
+#include "AliZDC.h"
+#include "AliZDCDigit.h"
+#include "AliZDCTrigger.h"
+
+//________________________________________________________________
+ClassImp(AliZDCTrigger)
+
+//________________________________________________________________
+AliZDCTrigger::AliZDCTrigger() : AliTriggerDetector() 
+{
+   SetName("ZDC");
+   CreateInputs();
+   //
+   SetZNMinCut(0);
+   SetZDCMinCut(0);
+   SetZEMMinCut(0);
+   SetZDCLeftEMDCuts(0,0);
+   SetZDCRightEMDCuts(0,0);
+   SetZDCMBCut(0);
+   SetZDCCentrCut(0);
+   SetZDCSemiCentrCut(0);
+   SetZEMCentrCut(0);
+
+}
+
+//________________________________________________________________
+void AliZDCTrigger::CreateInputs()
+{
+   // inputs 
+   
+   // Do not create inputs again!!
+   if( fInputs.GetEntriesFast() > 0 ) return;
+   
+   fInputs.AddLast(new AliTriggerInput("ZDC_1_L1",   "ZDC Minimum Bias", 0x01));
+   fInputs.AddLast(new AliTriggerInput("ZDC_2_L1",   "ZDC Central", 0x02));
+   fInputs.AddLast(new AliTriggerInput("ZDC_3_L1",   "ZDC Semi-central", 0x04));
+   fInputs.AddLast(new AliTriggerInput("ZDC_EMD_L1", "ZDC EMD events", 0x08));
+}
+
+//________________________________________________________________
+void AliZDCTrigger::Trigger()
+{
+
+
+   AliRunLoader *runLoader = gAlice->GetRunLoader();
+
+   AliLoader *ZDCLoader = runLoader->GetLoader("ZDCLoader");
+   ZDCLoader->LoadDigits("READ");
+   AliZDCDigit digit;
+   AliZDCDigit* pdigit = &digit;
+   TTree* TD = ZDCLoader->TreeD();
+   if (!TD) cerr<<"AliZDCTrigger: digits tree not found\n";
+   TD->SetBranchAddress("ZDC", &pdigit);
+   //
+   Float_t ZNSignal[2], ZPSignal[2], ZDCSumSignal[2], ZEMSumSignal[2];
+   for(Int_t iDigit=0; iDigit<TD->GetEntries(); iDigit++){
+      TD->GetEntry(iDigit);
+      //
+      if(digit.GetSector(0)==1)
+         for(Int_t i=0; i<2; i++){
+           ZNSignal[i] += digit.GetADCValue(i);
+           ZDCSumSignal[i] += digit.GetADCValue(i);
+         }
+      else if(digit.GetSector(0)==2)
+         for(Int_t i=0; i<2; i++){
+           ZPSignal[i] += digit.GetADCValue(i);
+           ZDCSumSignal[i] += digit.GetADCValue(i);
+         }
+      else if(digit.GetSector(0)==3)
+         for(Int_t i=0; i<2; i++) ZEMSumSignal[i] += digit.GetADCValue(i);
+   }
+   // *******************************************************************
+   if(ZNSignal[1]>fZDCLeftEMDCuts[0] && ZNSignal[1]<fZDCLeftEMDCuts[1] && 
+       ZEMSumSignal[1]<fZEMMinCut){ // *** n EMD event
+       SetInput("ZDC_EMD_L1");
+   }
+   // *******************************************************************
+   if(ZDCSumSignal[0]>fZDCMBCut) // *** ZDC minimum bias
+     SetInput("ZDC_1_L1");
+   // *******************************************************************
+   if(ZDCSumSignal[0]>fZDCMinCut && ZDCSumSignal[0]<fZDCMBCut) 
+     // *** ZDC central (0-10%)
+     SetInput("ZDC_2_L1");
+   // *******************************************************************
+   if(ZDCSumSignal[0]>fZDCCentrCut && ZDCSumSignal[0]<fZDCSemiCentrCut
+          && ZEMSumSignal[0]>fZEMCentrCut) 
+     // *** ZDC semi-central (10-40%)
+     SetInput("ZDC_3_L1");
+   
+}
+
+//________________________________________________________________
+void AliZDCTrigger::SetZNMinCut(Float_t ZNMinCut) 
+{
+  if(ZNMinCut)  fZNMinCut = ZNMinCut;
+  else  fZNMinCut = 400.;
+}
+
+//________________________________________________________________
+void AliZDCTrigger::SetZDCMinCut(Float_t ZDCMinCut) 
+{
+  if(ZDCMinCut)  fZDCMinCut = ZDCMinCut;
+  else  fZDCMinCut = 800.;
+}
+
+//________________________________________________________________
+void AliZDCTrigger::SetZEMMinCut(Float_t ZEMMinCut) 
+{
+  if(ZEMMinCut)  fZEMMinCut = ZEMMinCut;
+  else  fZEMMinCut = 80.;
+}
+//________________________________________________________________
+void AliZDCTrigger::SetZDCLeftEMDCuts(Float_t* ZDCLeftEMDCuts) 
+{
+  if(ZDCLeftEMDCuts) for(int j=0; j<2; j++) fZDCLeftEMDCuts[j] = ZDCLeftEMDCuts[j];
+  else{
+    fZDCLeftEMDCuts[0] = 600.;
+    fZDCLeftEMDCuts[1] = 1000.;
+  }
+}
+//________________________________________________________________
+void AliZDCTrigger::SetZDCLeftEMDCuts(Float_t ZDCLeftEMDCutInf, 
+       Float_t ZDCLeftEMDCutSup) 
+{
+  if(ZDCLeftEMDCutInf && ZDCLeftEMDCutSup){
+    fZDCLeftEMDCuts[0]=ZDCLeftEMDCutInf; 
+    fZDCLeftEMDCuts[1]=ZDCLeftEMDCutSup;
+  }    
+  else{
+    fZDCLeftEMDCuts[0] = 600.;
+    fZDCLeftEMDCuts[1] = 1000.;
+  }
+}
+//________________________________________________________________
+void AliZDCTrigger::SetZDCRightEMDCuts(Float_t* ZDCRightEMDCuts) 
+{
+  if(ZDCRightEMDCuts) for(int j=0; j<2; j++) fZDCRightEMDCuts[j] = ZDCRightEMDCuts[j];
+  else{
+    fZDCRightEMDCuts[0] = 600.;
+    fZDCRightEMDCuts[1] = 1000.;
+  }
+}
+//________________________________________________________________
+void AliZDCTrigger::SetZDCRightEMDCuts(Float_t ZDCRightEMDCutInf, 
+       Float_t ZDCRightEMDCutSup) 
+{
+  if(ZDCRightEMDCutInf && ZDCRightEMDCutSup){
+    fZDCRightEMDCuts[0]=ZDCRightEMDCutInf; 
+    fZDCRightEMDCuts[1]=ZDCRightEMDCutSup;
+  }    
+  else{
+    fZDCRightEMDCuts[0] = 600.;
+    fZDCRightEMDCuts[1] = 1000.;
+  }
+}
+//________________________________________________________________
+void AliZDCTrigger::SetZDCMBCut(Float_t ZDCMBCut) 
+{
+  if(ZDCMBCut) fZDCMBCut = ZDCMBCut;
+  else fZDCMBCut = 800.;
+}
+//________________________________________________________________
+void AliZDCTrigger::SetZDCCentrCut(Float_t ZDCCentrCut) 
+{
+  if(ZDCCentrCut) fZDCCentrCut = ZDCCentrCut;
+  else fZDCCentrCut = 10000.;
+}
+//________________________________________________________________
+void AliZDCTrigger::SetZDCSemiCentrCut(Float_t ZDCSemiCentrCut) 
+{
+  if(ZDCSemiCentrCut) fZDCSemiCentrCut = ZDCSemiCentrCut;
+  else fZDCSemiCentrCut = 18500.;
+}
+//________________________________________________________________
+void AliZDCTrigger::SetZEMCentrCut(Float_t ZEMCentrCut) 
+{
+  if(ZEMCentrCut) fZEMCentrCut = ZEMCentrCut;
+  else fZEMCentrCut = 210.;
+}

diff --git a/ZDC/AliZDCTrigger.h b/ZDC/AliZDCTrigger.h
new file mode 100644 (file)
index 0000000..f533e0f
--- /dev/null
+++ b/ZDC/AliZDCTrigger.h
@@ -0,0 +1,71 @@
+/* Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved. *
+ * See cxx source for full Copyright notice                               */
+
+#ifndef ALIZDCTRIGGER_H
+#define ALIZDCTRIGGER_H
+
+/// \ingroup sim
+/// \class AliZDCTrigger
+/// \brief ZDC trigger class
+///
+/////////////////////////////////////////////////
+///  ZDC Trigger Detector Class               //
+/////////////////////////////////////////////////
+
+#include "AliTriggerDetector.h"
+
+class AliZDCTrigger : public AliTriggerDetector
+{
+ public:
+   AliZDCTrigger();  // constructor
+   virtual ~AliZDCTrigger(){}  // destructor
+   virtual void    CreateInputs();
+   virtual void    Trigger();
+
+   // Print method
+   virtual void Print(Option_t *) const {
+     printf("\t AliZDCTrigger: fZNMinCut = %1.0f, fZDCMinCut = %1.0f, fZEMMinCut= %1.0f \n"
+     "fZDCLeftEMDCuts = [%1.0f, %1.0f], fZDCRightEMDCuts = [%1.0f, %1.0f], fZDCMBCut = %1.0f\n"
+     "fZDCCentrCut = %1.0f, fZDCSemiCentrCut = %1.0f, fZEMCentrCut = %1.0f\n\n",
+     fZNMinCut,fZDCMinCut,fZEMMinCut,fZDCLeftEMDCuts[0],fZDCLeftEMDCuts[1],
+     fZDCRightEMDCuts[0],fZDCRightEMDCuts[1],fZDCMBCut,fZDCCentrCut,fZDCSemiCentrCut,
+     fZEMCentrCut);
+   }
+ 
+ protected:
+   
+   // Setters   
+   void SetZNMinCut(Float_t ZNMinCut);
+   void SetZDCMinCut(Float_t ZDCMinCut);
+   void SetZEMMinCut(Float_t ZEMMinCut);
+   void SetZDCLeftEMDCuts(Float_t *ZDCLeftEMDCuts);
+   void SetZDCLeftEMDCuts(Float_t ZDCLeftEMDCutInf, Float_t ZDCLeftEMDCutSup);
+   void SetZDCRightEMDCuts(Float_t *ZDCRightEMDCuts);
+   void SetZDCRightEMDCuts(Float_t ZDCRightEMDCutInf, Float_t  ZDCRightEMDCutSup);
+   void SetZDCMBCut(Float_t ZDCMBCut);
+   void SetZDCCentrCut(Float_t ZDCCentrCuts);
+   void SetZDCSemiCentrCut(Float_t ZDCSemiCentrCut);
+   void SetZEMCentrCut(Float_t ZEMCentrCut);
+
+   // Data member
+   Float_t fZNMinCut;   
+   Float_t fZDCMinCut;   
+   Float_t fZEMMinCut;   
+   Float_t fZDCLeftEMDCuts[2];
+   Float_t fZDCRightEMDCuts[2];
+   Float_t fZDCMBCut;
+   Float_t fZDCCentrCut;
+   Float_t fZDCSemiCentrCut;
+   Float_t fZEMCentrCut;   
+    
+  ClassDef(AliZDCTrigger,1)  // ZDC Trigger Detector class
+};
+#endif
+
+
+
+
+
+
+
+