In This Package:

AdQuality.cc

Go to the documentation of this file.

#include "AdQuality.h"

#include "Context/Context.h"
#include "Conventions/Detectors.h"
#include "DataSvc/ICableSvc.h"
#include "Event/ReadoutHeader.h"
#include "Event/Readout.h"
#include "Event/ReadoutPmtCrate.h"
#include "Event/ReadoutPmtChannel.h"

#include "GaudiKernel/ITHistSvc.h"

#include "TH2F.h"
#include "TH1F.h"

using namespace std;
using namespace DayaBay;

class AdQualityPlots {
    DayaBay::Detector detector;
    MsgStream& log;
public:

    // direct access
    TH2F* hitMap2d;
    TH1F* hitMap1d;

    AdQualityPlots(DayaBay::Detector det, MsgStream& msg) : detector(det), log(msg) { }

    StatusCode book(ITHistSvc* hsvc, string filepath) {

        {
            string name = detector.detName() + "_HitMap2D";
            string title = "Hit map for " + detector.detName();
            hitMap2d = new TH2F(name.c_str(),title.c_str(),
                                24,0,24, 8,0,8);
            
            name = filepath+"/"+detector.detName()+"/hitMap2D";
            if (hsvc->regHist(name,hitMap2d).isFailure()) {
                log << MSG::ERROR << "Could not register " << name << endreq;
                delete hitMap2d; hitMap2d = 0;
                return StatusCode::FAILURE;
            }
        }

        {
            string name = detector.detName() + "_HitMap1D";
            string title = "Hit map for " + detector.detName();
            hitMap1d = new TH1F(name.c_str(), title.c_str(),
                                8*24,0,8*24);
            hitMap1d->Sumw2();
            name = filepath+"/"+detector.detName()+"/hitMap1D";
            if (hsvc->regHist(name,hitMap1d).isFailure()) {
                log << MSG::ERROR << "Could not register " << name << endreq;
                delete hitMap1d; hitMap1d = 0;
                return StatusCode::FAILURE;
            }
        }

        return StatusCode::SUCCESS;
    }

    StatusCode fill(ICableSvc* cable, const ServiceMode &sm, 
                    const ReadoutPmtCrate* ro) 
    {

        // std::map<DayaBay::FeeChannelId,DayaBay::ReadoutPmtChannel> 
        const ReadoutPmtCrate::PmtChannelReadouts& 
            channelMap = ro->channelReadout();
        ReadoutPmtCrate::PmtChannelReadouts::const_iterator
            it, done = channelMap.end();
        for (it=channelMap.begin(); it != done; ++it) {

            float charge = it->second.sumAdc();

            DayaBay::AdPmtSensor pmt = cable->adPmtSensor(it->first,sm);

            hitMap1d->Fill(pmt.column()-1.0 + 24*(pmt.ring()-1) + 0.5, charge);
            hitMap2d->Fill(pmt.column()-0.5, pmt.ring()-0.5,charge);
        }
        return StatusCode::SUCCESS;
    }

};


AdQuality::AdQuality(const std::string& name, ISvcLocator* pSvcLocator)
    : GaudiAlgorithm(name,pSvcLocator)
    , m_hsvc(0)
{
    declareProperty("Location",m_location=DayaBay::ReadoutHeaderLocation::Default,
                    "Location in the TES to ReadoutHeaders");
    declareProperty("FilePath",m_filepath="/file1/adqual/",
                    "File path with with to register histograms.");
}

AdQuality::~AdQuality()
{
}

StatusCode AdQuality::initialize()
{
    this->GaudiAlgorithm::initialize();

    // Get the histogram service
    if ( service("THistSvc", m_hsvc).isFailure()) {
        error() << " No THistSvc available." << endreq;
        return StatusCode::FAILURE;
    } 
    // Get the cable service
    if ( service("StaticCableSvc", m_cable).isFailure()) {
        error() << " No CableService available." << endreq;
        return StatusCode::FAILURE;
    } 

    return StatusCode::SUCCESS;
}

StatusCode AdQuality::execute()
{
    DayaBay::ReadoutHeader* roh = get<DayaBay::ReadoutHeader>(m_location);

    const DayaBay::ReadoutPmtCrate* ro = 
        dynamic_cast<const DayaBay::ReadoutPmtCrate*>(roh->readout());
    if (!ro) {
        debug() << "Got a non-PmtCrate readout, skipping..." << endreq;
        return StatusCode::SUCCESS;
    }

    const Context &ctx = roh->context();
    DayaBay::Detector det(ctx.GetSite(),ctx.GetDetId());
    AdQualityPlots* plots = m_plots[det];
    if (!plots) {
        plots = new AdQualityPlots(det,info());
        m_plots[det] = plots;
        plots->book(m_hsvc,m_filepath);
    }

    return plots->fill(m_cable,ServiceMode(ctx,0),ro);
}

StatusCode AdQuality::finalize()
{
    return this->GaudiAlgorithm::finalize();
}

---