coltrims_analysis.cpp

Go to the documentation of this file.

//Copyright (C) 2011 Lutz Foucar

/**
 * @file coltrims_analysis.cpp file contains the processor specific for
 *                             coltrims analysis
 *
 * @author Lutz Foucar
 */

#include <stdexcept>
#include <cmath>
#include <algorithm>

#include <QtCore/QString>

#include "coltrims_analysis.h"

#include "acqiris_detectors_helper.h"
#include "cass_event.h"
#include "acqiris_device.hpp"
#include "cass.h"
#include "convenience_functions.h"
#include "cass_settings.h"
#include "log.h"

using namespace cass;
using namespace ACQIRIS;
using namespace std;


//----------------Electron Energy----------------------------------------------
pp5000::pp5000(const name_t &name)
  : Processor(name)

{
  loadSettings(0);
}

void pp5000::loadSettings(size_t)
{
  CASSSettings settings;
  settings.beginGroup("Processor");
  settings.beginGroup(QString::fromStdString(name()));
  _detector = loadDelayDet(settings,5000,name());
  _particle = loadParticle(settings,_detector,5000,name());
  setupGeneral();
  if (!setupCondition())
    return;
  createHistList(set1DHist(name()));
  Log::add(Log::INFO,"Processor '" + name() +
           "' extracts the total momentum of the particle '" + _particle +
           "' of detector '" + _detector +
           "' and then converts it to Energy assuming that it is an electron. " +
           "Condition is '" + _condition->name() + "'");
}

void pp5000::process(const CASSEvent& evt, result_t &result)
{
  DetectorBackend &rawdet(
        HelperAcqirisDetectors::instance(_detector)->detector(evt));
  DelaylineDetector &det (dynamic_cast<DelaylineDetector&>(rawdet));
  Particle &particle(det.particles()[_particle]);
  particleHits_t::iterator it(particle.hits().begin());
  for (; it != particle.hits().end(); ++it)
  {
    double p = (*it)[roh];
    double e_energy = p*p*13.6;
    result.histogram(e_energy);
  }
}



//----------------PIPIPICO-------------------------------------------------------
pp5001::pp5001(const name_t &name)
  :cass:: Processor(name)
{
  loadSettings(0);
}

void pp5001::loadSettings(size_t)
{
  CASSSettings settings;
  settings.beginGroup("Processor");
  settings.beginGroup(QString::fromStdString(name()));
  _detector = settings.value("Detector","blubb").toString().toStdString();
  setupGeneral();
  if (!setupCondition())
    return;
  createHistList(set2DHist(name()));
  HelperAcqirisDetectors::instance(_detector)->loadSettings();
  Log::add(Log::INFO,"Processor '" + name() +
           "' create a tripple coincidence Histogram of particles of '" + _detector +
           "'. Condition is '" + _condition->name() + "'");
}

void pp5001::process(const CASSEvent& evt, result_t &result)
{
  DetectorBackend &rawdet(
        HelperAcqirisDetectors::instance(_detector)->detector(evt));
  DelaylineDetector &det01 (dynamic_cast<DelaylineDetector&>(rawdet));
  SignalProducer::signals_t::const_iterator it01(det01.mcp().output().begin());
  SignalProducer::signals_t::const_iterator end(det01.mcp().output().end());
  for (; it01 != end;++it01)
  {
    SignalProducer::signals_t::const_iterator it02(it01+1);
    for (; it02 != end; ++it02)
    {

      SignalProducer::signals_t::const_iterator it03(it02+1);
      for (; it03 != end; ++it03)
      {
        result.histogram(make_pair((*it01)[ACQIRIS::time]+(*it02)[ACQIRIS::time],
                                   (*it03)[ACQIRIS::time]));
      }
    }
  }
}