BmnGemDigitQa.h

Go to the documentation of this file.

// BmnGemDigitQa: example task for DIGI files - strip signal spectra of the GEM stations
class BmnGemDigitQa : public FairTask
{
  public:
    BmnGemDigitQa(Bool_t isExp = kFALSE)
        : FairTask("BmnGemDigitQa")
        , fIsExp(isExp)
    {}
 
    InitStatus Init() override
    {
        // the GEM digits have different branch names in experimental and simulated files
        FairRootManager* ioman = FairRootManager::Instance();
        fDigits = (TClonesArray*)ioman->GetObject(fIsExp ? "GEM" : "BmnGemStripDigit");
        if (!fDigits) {
            LOG(error) << "BmnGemDigitQa::Init(): GEM digit branch not found, the task is deactivated";
            return kERROR;
        }
        fHistSignal = new TH2D("hGemSignal", "GEM strip signal;station;signal, ADC counts", 8, 0, 8, 100, 0, 1000);
        return kSUCCESS;
    }
 
    void Exec(Option_t* /*option*/) override
    {
        for (Int_t i = 0; i < fDigits->GetEntriesFast(); i++) {
            BmnStripDigit* digit = (BmnStripDigit*)fDigits->UncheckedAt(i);
            fHistSignal->Fill(digit->GetStation(), digit->GetStripSignal());
        }
    }
 
    void Finish() override
    {
        // write the results to the output file of the macro
        FairRootManager::Instance()->GetSink()->WriteObject(fHistSignal, "hGemSignal");
    }
 
  private:
    Bool_t fIsExp;
    TClonesArray* fDigits = nullptr;
    TH2D* fHistSignal = nullptr;
 
    ClassDefOverride(BmnGemDigitQa, 1);
};