#include "Gate1D.hpp"
#include <iostream>
#include <cmath>

Gate1D::Gate1D(const std::string &config_dir,
	       const std::string &name, 
	       const bool no_files)
  : Processor(name, config_dir, no_files, std::vector<int>())
{
  // make all the in/outputs, parameters and coefficients 
  // known by name.
  NAME_PARAMETER(num_channels);
  init();
  read_parameters();
  
  NAME_CONDITION_ARRAY(reference_gate, Window1D, parameter(num_channels));

  NAME_INPUT_ARRAY(value, parameter(num_channels));
  NAME_INPUT_ARRAY(value2test, parameter(num_channels));
  
  NAME_OUTPUT_ARRAY(GatedValue, parameter(num_channels));
  

  
  init();
  read_conditions();
}

Gate1D::~Gate1D()
{
}


void Gate1D::process(prespec::viscon::Interface &viscon_interface, int trigger)
{

int chnr = parameter(num_channels);		  

  for (int i = 0; i<chnr; i++)
	{

double left_gate = condition_array(reference_gate, i, 0);
double right_gate = condition_array(reference_gate, i, 1);

  if(input_array_value(value, i) == 0 && input_array_value(value2test, i) == 0 )
   return;

	if (input_array_value(value2test, i) > left_gate && input_array_value(value2test, i) < right_gate)
	{
	double girdi = input_array_value(value, i);
	fill_output_array(GatedValue, i, girdi);
	}
    
  }
  
}

go4analysis: prespec/process/Processor.cpp:665: const double& prespec::process::Processor::input_array_value(int, int): Assertion `input_arrays_[channel].size() > i' failed.

cp Gate1D.hpp MyNewProcessor.hpp; cp Gate1D.cpp MyNewProcessor.cpp

#include "process/MyNewProcessor.hpp"

ELDER_MODULE_PROCESSOR_BEGIN;
...
ELDER_MODULE_PROCESSOR_ADD(MyNewProcessor);
...
ELDER_MODULE_PROCESSOR_END;

libprespecUTILS_la_SOURCES = ...    \
                             ...    \
                             plugins/UTILS/process/MyNewProcessor.cpp \
                             plugins/UTILS/process/MyNewProcessor.hpp \
                             ... 

Assertion `input_arrays_[channel].size() > i' failed

NAME_INPUT_ARRAY(my_array, 4) // 4 different indices allowed: 0,1,2,3

       |entry0 | entry1 | entry2 | entry3 |
-------------------------------------------
index  |  1    |    3   |   2    |   3    |
value  |  3.2  |  100.5 |   3.3  |  234.5 |
-------------------------------------------

       |entry0 | entry1 | entry2 | entry3 | entry4 | entry5 |
-------------------------------------------------------------
index  |  3    |    0   |   1    |   2    |   2    |   1    |
value  |  3.2  |  100.5 |   3.3  |  234.5 |  150.1 |   3.14 |
--------------------------------------------------------------

       |
--------
index  |
value  |
--------

int N = input_array_size(my_array); // N can be any positive integer or zero
for (int i = 0; i < N; ++i)
{
  int    index = input_array_index(my_array, i); // index refers to the index of the i-th array entry
  double value = input_array_value(my_array, i); // value refers to the value of the i-th array entry


  // do something with the value and the index
  // .. for example testing
}

input_array_index(my_array, 5); // this will cause an assertion failure if the number of entries is smaller than 6

for( //...
{
  int    index = input_array_index(my_array, i);
  double value = input_array_value(my_array, i);
  //...

  // the following works if 'reference_gate' is defined 
  //   with 'NAME_CONDITION_ARRAY(reference_gate, Window1D, parameter(num_channels))'
  if (condition_valid(value, reference_gate, index)) 
  {
    //...
  }
}

#include "ArraysWithCondition.hpp"
#include <iostream>
#include <cmath>

ArraysWithCondition::ArraysWithCondition(const std::string &config_dir,
	       const std::string &name, 
	       const bool no_files)
  : Processor(name, config_dir, no_files, std::vector<int>())
{
  // make all the in/outputs, parameters and coefficients 
  // known by name.


  NAME_INPUT_ARRAY(value1, 7);
  NAME_INPUT_ARRAY(value2, 7);
  
  NAME_OUTPUT_ARRAY(GatedValue, 7);
  
  NAME_CONDITION_ARRAY(gate, Window1D, 7);
  
  init();
  read_conditions();
}

ArraysWithCondition::~ArraysWithCondition()
{
}

void ArraysWithCondition::process(prespec::viscon::Interface &viscon_interface, int trigger)
{
  

	for(int i = 0; i < input_array_size(value1); i++)
	{
		int index = input_array_index(value1, i);
		double val2 = input_array_value(value2, index);
		double val1 = input_array_value(value1, index);
		double left_gate = condition_array(gate, index, 0);
		double right_gate = condition_array(gate, index, 1);
		if(val2 < right_gate && val2 > left_gate)
		{
			fill_output_array(GatedValue, index, val1);

		}

	}

  
  
}

#ifndef UTILS_ArraysWithCondition_HPP_
#define UTILS_ArraysWithCondition_HPP_

#include <prespec/process/Processor.hpp>


//! @brief A processor that evaluates the information coming from a
//!        thin circular plastic membrane detectors, surrounded by 
//!        many photomultipliers. These detectors are used in the
//!        Lycca time of flight subsystem.  WTF!!!!! wrong doc!
//! 
//! It evaluates the time of the particle impact on the membrane, given
//! the particle position (deduced by other tracking detectors).
//! If this information is missing, the processor tries to deduce the 
//! impact position from the time information. That is possible with an
//! accuracy of between (6-10 milimeters).
class ArraysWithCondition : public prespec::process::Processor
{
public:
  ArraysWithCondition(const std::string &config_dir, 
	 const std::string &name, 
	 const bool no_files = false);
  ~ArraysWithCondition();
  
  virtual void process(prespec::viscon::Interface &viscon_interface, int trigger);
  
  double analyze(double x_particle, double y_particle, bool find_position = false);
  
  // the inputs to this processor
  enum InputArray
    {
      value1,
      value2,
    };
  
  // the outputs of this processor
  enum OutputArray
    {
      GatedValue,
    };
  
 
   enum Conditions
    {
      gate,
    };
   
};


#endif

processor HighLevel/particle/incoming/MusicPileUp	UTILS.ArraysWithCondition
	value1[0:6] <- MuiscPileupCor.music[0:6]
	value2[0:6] <- MuiscPileupCor.music[8:14]

	display GatedValue
end

gate[0]		25	160
gate[1]		25	160
gate[2]		25	160
gate[3]		25	160
gate[4]		25	160
gate[5]		25	160
gate[6]		25	160

prespec/process/Processor.cpp:665: const double& prespec::process::Processor::input_array_value(int, int): Assertion `input_arrays_[channel].size() > i' failed.

// check first if both arrays have the same number of entries
if (input_array_size(value1) == input_array_size(value2)) 
{
  for (int i = 0; i < input_array_size(value1); ++i)
  {
    // At this point, your code has a problem, because you use "index" instead of "i".
    // The second parameter in the functions input_array_index/input_array_value
    // must not exceed the size of the array (i.e. the number of entries of the array) 
    // Otherwise you'll get an assertion failure.                                                             
    int    index_of_ith_entry_in_value1 = input_array_index(value1, i); 
    double value_of_ith_entry_in_value1 = input_array_value(value1, i); 
    int    index_of_ith_entry_in_value2 = input_array_index(value2, i); 
    double value_of_ith_entry_in_value2 = input_array_value(value2, i); 

    // eventually, you want to check if the indices are the same
    if (index_of_ith_entry_in_value1 == index_of_ith_entry_in_value2)
    {
      if (condition_valid(value_of_ith_entry_in_value2, gate, index_of_ith_entry_in_value1))
      {
        fill_output_array(GatedValue, index_of_ith_entry_in_value2, value_of_ith_entry_in_value1)
      }
    }
  }
}