#include "AliLog.h"
#include "AliRun.h"
#include "AliLoader.h"
#include "AliRunLoader.h"
#include "AliTriggerInput.h"
#include "AliT0Parameters.h"
#include "AliT0TriggerParameters.h"
#include <AliCDBManager.h>
#include <AliCDBEntry.h>
#include <AliCDBStorage.h>
#include "AliT0.h"
#include "AliT0digit.h"
#include "AliT0Trigger.h"
#include "Riostream.h"
ClassImp(AliT0Trigger)
AliT0Trigger::AliT0Trigger()
: AliTriggerDetector(),
fT0(0x0),
fDigits(0x0),
fTrigPar(0x0)
{
SetName("T0");
CreateInputs();
AliCDBManager *stor =AliCDBManager::Instance();
AliCDBEntry* fCalibentry = stor->Get("T0/Calib/TriggerParam");
if (fCalibentry)
fTrigPar = (AliT0TriggerParameters*)fCalibentry->GetObject();
else {
AliWarning(" No trigger parameters in CDB , use default");
}
}
void AliT0Trigger::CreateInputs()
{
if( fInputs.GetEntriesFast() > 0 ) return;
fInputs.AddLast( new AliTriggerInput( "0T0A", "T0", 0 ) );
fInputs.AddLast( new AliTriggerInput( "0T0C", "T0", 0 ) );
fInputs.AddLast( new AliTriggerInput( "0TVX", "T0", 0 ) );
fInputs.AddLast( new AliTriggerInput( "T0_Centr_L0", "T0", 0 ) );
fInputs.AddLast( new AliTriggerInput( "T0_SemiCentral_L0", "T0", 0 ) );
}
void AliT0Trigger::Trigger()
{
AliRunLoader* runLoader = AliRunLoader::Instance();
AliLoader * fT0Loader = runLoader->GetLoader("T0Loader");
fT0Loader->LoadDigits("READ");
TTree* treeD = fT0Loader->TreeD();
if (!treeD) {
AliError("no digits tree");
return;
}
fDigits = new AliT0digit();
TBranch *brDigits = treeD->GetBranch("T0");
if (brDigits) {
brDigits->SetAddress(&fDigits);
}else{
AliError("Branch T0 DIGIT not found");
return;
}
brDigits->GetEntry(0);
Int_t besttimeA = fDigits->BestTimeA();
Int_t besttimeC = fDigits->BestTimeC();
Int_t timeDiff = fDigits->TimeDiff();
Int_t sumMult= fDigits->SumMult();
Float_t timeWindowLow = fTrigPar->GetTimeWindowLow();
Float_t timeWindowHigh = fTrigPar->GetTimeWindowHigh();
Int_t ampCentr = fTrigPar->GetAmpCentr();
Int_t ampSemiCentr = fTrigPar->GetAmpSemiCentr();
if (besttimeA > 0 && besttimeA <99999) SetInput("0T0A");
if (besttimeC > 0 && besttimeC<99999) SetInput("0T0C");
if (timeDiff >timeWindowLow && timeDiff < timeWindowHigh) SetInput("0TVX");
if (sumMult > ampCentr) SetInput("T0_Centr_L0");
if (sumMult> ampSemiCentr && sumMult <= ampCentr) SetInput("T0_SemiCentral_L0");;
}