BmnTriggersCheck.cxx

Go to the documentation of this file.

#include "BmnTriggersCheck.h"
 
#include "BmnTrigDigit.h"
 
BmnTriggersCheck::BmnTriggersCheck(Bool_t isExp, Int_t runP, Int_t runN)
{
    fIsExp = isExp;
    fRunPeriod = runP;
    fRunNumber = runN;
 
    fT0Branch = "T0";
    fVetoBranch = "VC";
    fBC2Branch = "BC2";
    fBC1Branch = "BC1";
    fBC3Branch = "BC3";
    fBDBranch = "BD";
    fBmnEventHeaderBranchName = "EventHeader";
    fBmnEvQualityBranchName = "BmnEventQuality";
}
 
InitStatus BmnTriggersCheck::Init()
{
    FairRootManager* ioman = FairRootManager::Instance();
 
    fBmnEventHeader = (TClonesArray*)ioman->GetObject(fBmnEventHeaderBranchName);
    fT0Array = (TClonesArray*)ioman->GetObject(fT0Branch.Data());
    fVetoArray = (TClonesArray*)ioman->GetObject(fVetoBranch.Data());
    fBC2Array = (TClonesArray*)ioman->GetObject(fBC2Branch.Data());
    fBDArray = (TClonesArray*)ioman->GetObject(fBDBranch.Data());
    fBC1Array = (TClonesArray*)ioman->GetObject(fBC1Branch.Data());
    fBC3Array = (TClonesArray*)ioman->GetObject(fBC3Branch.Data());
 
    fBmnEvQuality = new TClonesArray(fBmnEvQualityBranchName);
    ioman->Register(fBmnEvQualityBranchName, "QUALITY", fBmnEvQuality, true);
 
    return kSUCCESS;
}
 
void BmnTriggersCheck::Exec(Option_t* opt)
{
    if (!fIsExp)
        return;
 
    fBmnEvQuality->Delete();
 
    BmnEventQuality* evQual = new ((*fBmnEvQuality)[fBmnEvQuality->GetEntriesFast()]) BmnEventQuality("GOOD");
 
    Int_t cntr2 = 0;
    Bool_t time2 = 0;
    for (Int_t iDig = 0; iDig < fBC2Array->GetEntriesFast(); ++iDig) {
        BmnTrigDigit* dig = (BmnTrigDigit*)fBC2Array->At(iDig);
        if (dig->GetMod() == 0) {
            cntr2++;
            //  if (dig->GetTime() < 420 || dig->GetTime() > 600) time2 =1;
            // cout<<" time BC2 "<<dig->GetTime()<<" time2 "<<time2<<endl;
        }
    }
    if (cntr2 > 2 || time2 == 1)
        evQual->SetIsGoodEvent("BAD");
 
    Int_t cntr1 = 0;
    Bool_t time1 = 0;
    for (Int_t iDig = 0; iDig < fBC1Array->GetEntriesFast(); ++iDig) {
        BmnTrigDigit* dig = (BmnTrigDigit*)fBC1Array->At(iDig);
        if (dig->GetMod() == 0) {
            cntr1++;
            // if (dig->GetTime() < 460 ||  dig->GetTime() > 600) time1 =1;
            // cout<<" time BC1 "<<dig->GetTime()<<" time1 "<<time1<<endl;
        }
    }
    if (cntr1 > 2 || time1 == 1)
        evQual->SetIsGoodEvent("BAD");
 
    Int_t cntr3 = 0;
    Bool_t time3 = 0;
    for (Int_t iDig = 0; iDig < fBC3Array->GetEntriesFast(); ++iDig) {
        BmnTrigDigit* dig = (BmnTrigDigit*)fBC3Array->At(iDig);
        if (dig->GetMod() == 0) {
            cntr1++;
            //  if (dig->GetTime() < 470 || dig->GetTime() > 620) time3 =1;
            // cout<<" time BC3 "<<dig->GetTime()<<" time3 "<<time3<<endl;
        }
    }
    if (cntr3 > 2 || time3 == 1)
        evQual->SetIsGoodEvent("BAD");
 
    // if (fVetoArray != NULL){
    Int_t cntrv = 0;
    for (Int_t iDig = 0; iDig < fVetoArray->GetEntriesFast(); ++iDig) {
        BmnTrigDigit* dig = (BmnTrigDigit*)fVetoArray->At(iDig);
        if (dig->GetMod() == 0)
            cntrv++;
    }
    if (cntrv >= 1)
        evQual->SetIsGoodEvent("BAD");
 
    // } else cout<<" fVetoArray not found "<<endl;
 
    //  BmnEventQuality* evQual = new ((*fBmnEvQuality)[fBmnEvQuality->GetEntriesFast()]) BmnEventQuality("GOOD");
 
    //     const Int_t kEnergies = 4;
    //     const Int_t kTargets = 5;
    //     Double_t energies[kEnergies] = {3.5, 4., 4.5, 5.14};
    //     TString targets[kTargets] = {"C", "Al", "Cu", "Pb", "C2H4"};
 
    //     BmnEventHeader* evHeader = (BmnEventHeader*) fBmnEventHeader->UncheckedAt(0);
    //     if (!evHeader)
    //         return;
 
    //     UniRun* runInfo = UniRun::GetRun(fRunPeriod, fRunNumber);
 
    //     BmnTriggerType trigType = evHeader->GetTrig();
    //     if (trigType == kBMNMINBIAS) {
    //         // Setup BD-threshold
    //         // To be read subsequently from the UniDB directly when processing RUN7, FIXME
    //         TString target = *runInfo->GetTargetParticle();
    //         Double_t energy = *runInfo->GetEnergy();
    //         Int_t bdThresh;
    //         if (Abs(energy - energies[0]) < FLT_EPSILON) // 3.5 GeV/n, all targets
    //             bdThresh = 2;
    //         else if (Abs(energy - energies[1]) < FLT_EPSILON && (target == targets[0] || target == targets[4])) // 4
    //         GeV/n, C, C2H4
    //             bdThresh = 2;
    //         else if (Abs(energy - energies[1]) < FLT_EPSILON && target != targets[0] && target != targets[4]) // 4
    //         GeV/n, Al, Cu, Pb
    //             bdThresh = 3;
    //         else if (Abs(energy - energies[2]) < FLT_EPSILON && (target == targets[0] || target == targets[4]))
    //         // 4.5 GeV/n, C, C2H4
    //             bdThresh = 2;
    //         else if (Abs(energy - energies[2]) < FLT_EPSILON && target != targets[0] && target != targets[4]) // 4
    //         GeV/n, Al, Cu, Pb
    //             bdThresh = 3;
    //         else if (Abs(energy - energies[3]) < FLT_EPSILON) // 5.14 GeV/n, all targets
    //             bdThresh = 2;
    //         else
    //             bdThresh = 0;
 
    //         if (fT0Array->GetEntriesFast() != 1 || fBC2Array->GetEntriesFast() != 1 || fVetoArray->GetEntriesFast()
    //         != 0 || fBDArray->GetEntriesFast() < bdThresh)
    //             evQual->SetIsGoodEvent("BAD");
    //     }
 
    //     else if (trigType == kBMNBEAM) {
    //         if (fVetoArray->GetEntriesFast() > 1 || fBDArray->GetEntriesFast() > 0)
    //             evQual->SetIsGoodEvent("BAD");
    //     }
 
    //     else
    //         evQual->SetIsGoodEvent("BAD");
}
 
void BmnTriggersCheck::Finish() {}