[u/mrichter/AliRoot.git] / HLT / PHOS / AliHLTPHOSSanityInspector.cxx

// $Id$

/**************************************************************************
 * This file is property of and copyright by the ALICE HLT Project        * 
 * All rights reserved.                                                   *
 *                                                                        *
 * Primary Authors: Oystein Djuvsland                                     *
 *                                                                        *
 * Permission to use, copy, modify and distribute this software and its   *
 * documentation strictly for non-commercial purposes is hereby granted   *
 * without fee, provided that the above copyright notice appears in all   *
 * copies and that both the copyright notice and this permission notice   *
 * appear in the supporting documentation. The authors make no claims     *
 * about the suitability of this software for any purpose. It is          * 
 * provided "as is" without express or implied warranty.                  *
 **************************************************************************/

/** 
 * @file   AliHLTPHOSSanityInspector.cxx
 * @author Oystein Djuvsland
 * @date 
 * @brief  Sanity inspector for PHOS HLT 
 */

// see header file for class documentation
// or
// refer to README to build package
// or
// visit http://web.ift.uib.no/~kjeks/doc/alice-hlt

#include "AliHLTPHOSSanityInspector.h"
#include "Rtypes.h"

ClassImp(AliHLTPHOSSanityInspector);


AliHLTPHOSSanityInspector::AliHLTPHOSSanityInspector() : 
  //  AliHLTPHOSBase(),
  fMaxDifference(120)
{
  //See header file for documentation
}


AliHLTPHOSSanityInspector::~AliHLTPHOSSanityInspector()
{
  //See header file for documentation
}


Int_t  
AliHLTPHOSSanityInspector::CheckInsanity(const UInt_t* data, const Int_t N) const
{
   //See header file for documentation

  for(Int_t i = 1; i < N; i++)
    {
      if((((Int_t)data[i] - (Int_t)data[i-1]) > fMaxDifference) || (((Int_t)data[i-1] - (Int_t)data[i]) > fMaxDifference))
	return 1;
    }
  return 0;
}

/*
Int_t
AliHLTPHOSSanityInspector::CheckInsanity(Int_t* data, Int_t N)
{
   //See header file for documentation
  for(Int_t i = 1; i < N; i++)
  {
    if((((Int_t)data[i] - (Int_t)data[i-1]) > fMaxDifference) || (((Int_t)data[i-1] - (Int_t)data[i]) > fMaxDifference))
      return 1;
  }
  return 0;
}
*/


Int_t 
AliHLTPHOSSanityInspector::CheckAndHealInsanity(UInt_t* data, Int_t N)
  //
{
   //See header file for documentation

  Int_t crazyness = 0;

  if(N > 3)
    {
      //Require a stable start 
      if((((Int_t)data[0] - (Int_t)data[1]) > fMaxDifference) || (((Int_t)data[1] - (Int_t)data[0]) > fMaxDifference))
	return -1;
      if((((Int_t)data[1] - (Int_t)data[2]) > fMaxDifference) || (((Int_t)data[2] - (Int_t)data[1]) > fMaxDifference))
	return -1;


      for(Int_t i = 2; i < N - 3; i++)
	{
	  if((((Int_t)data[i] - (Int_t)data[i+1]) > fMaxDifference) || (((Int_t)data[i+1] - (Int_t)data[i]) > fMaxDifference))
	    {
	      i++;
	      if((((Int_t)data[i] -(Int_t)data[i+1]) > fMaxDifference) || (((Int_t)data[i+1] - (Int_t)data[i]) > fMaxDifference))
		{
		  i++;
		  if((((Int_t)data[i] - (Int_t)data[i+1]) > fMaxDifference) || (((Int_t)data[i+1] - (Int_t)data[i]) > fMaxDifference))
		    {
		      return -2;  //Too crazy
		    }
		  data[i-1] = ((Int_t)data[i] + (Int_t)data[i-2])/2;
		  crazyness++;
		}
	      else 
		return -3;    //Two spikes in a row? 
	    }
	}
      
      
      if((((Int_t)data[N - 3] -(Int_t) data[N - 2]) > fMaxDifference) || 
	 (((Int_t)data[N - 2] - (Int_t)data[N - 3]) > fMaxDifference))
	{
	  if((((Int_t)data[N - 2] - (Int_t)data[N - 1]) > fMaxDifference) || 
	     (((Int_t)data[N - 1] - (Int_t)data[N - 2]) > fMaxDifference))
	    {
	      data[N - 2] = ((Int_t)data[N - 3] +  (Int_t)data[N - 1])/2;
	      return crazyness++;
	    }
	  return -4;

	}
      
      if((((Int_t)data[N - 2] - (Int_t)data[N - 1]) > fMaxDifference) || 
	 (((Int_t)data[N - 1] - (Int_t)data[N - 2]) > fMaxDifference))
	{
	  //	  (Int_t)data[N - 3] = (Int_t)data[N - 4] -(Int_t) data[N - 5] + (Int_t)data[N-4];
	  data[N - 1] = data[N - 2];
	  return crazyness++;
	}
      
    }
  
  return crazyness;
  
}


Int_t 
AliHLTPHOSSanityInspector::CheckAndHealInsanity(Int_t* data, Int_t N)
{
   //See header file for documentation

  Int_t crazyness = 0;

  if(N > 3)
    {
      if((((Int_t)data[0] - (Int_t)data[1]) > fMaxDifference) || (((Int_t)data[1] - (Int_t)data[0]) > fMaxDifference))
	return -1;
      if((((Int_t)data[1] - (Int_t)data[2]) > fMaxDifference) || (((Int_t)data[2] - (Int_t)data[1]) > fMaxDifference))
	return -1;


      for(Int_t i = 2; i < N - 3; i++)
	{
	  if((((Int_t)data[i] - (Int_t)data[i+1]) > fMaxDifference) || (((Int_t)data[i+1] - (Int_t)data[i]) > fMaxDifference))
	    {
	      i++;
	      if((((Int_t)data[i] -(Int_t)data[i+1]) > fMaxDifference) || (((Int_t)data[i+1] - (Int_t)data[i]) > fMaxDifference))
		{
		  i++;
		  if((((Int_t)data[i] - (Int_t)data[i+1]) > fMaxDifference) || (((Int_t)data[i+1] - (Int_t)data[i]) > fMaxDifference))
		    {
		      return -2;  //Too crazy
		    }
		  data[i-1] = ((Int_t)data[i] + (Int_t)data[i-2])/2;
		  crazyness++;
		}
	      else 
		return -3;    //Two spikes in a row? 
	    }
	}
      
      
      if((((Int_t)data[N - 3] -(Int_t) data[N - 2]) > fMaxDifference) || 
	 (((Int_t)data[N - 2] - (Int_t)data[N - 3]) > fMaxDifference))
	{
	  if((((Int_t)data[N - 2] - (Int_t)data[N - 1]) > fMaxDifference) || 
	     (((Int_t)data[N - 1] - (Int_t)data[N - 2]) > fMaxDifference))
	    {
	      data[N - 2] = ((Int_t)data[N - 3] +  (Int_t)data[N - 1])/2;
	      return crazyness++;
	    }
	  return -4;

	}
      
      if((((Int_t)data[N - 2] - (Int_t)data[N - 1]) > fMaxDifference) || 
	 (((Int_t)data[N - 1] - (Int_t)data[N - 2]) > fMaxDifference))
	{
	  //	  (Int_t)data[N - 3] = (Int_t)data[N - 4] -(Int_t) data[N - 5] + (Int_t)data[N-4];
	  data[N - 1] = data[N - 2];
	  return crazyness++;
	}
      
    }
  
  return crazyness;
  
}
Commit	Line	Data
1b41ab20	1	// $Id$
1b41ab20	2
3ed012c1	3	/**************************************************************************
	4	* This file is property of and copyright by the ALICE HLT Project *
	5	* All rights reserved. *
	6	* *
	7	* Primary Authors: Oystein Djuvsland *
	8	* *
	9	* Permission to use, copy, modify and distribute this software and its *
	10	* documentation strictly for non-commercial purposes is hereby granted *
	11	* without fee, provided that the above copyright notice appears in all *
	12	* copies and that both the copyright notice and this permission notice *
	13	* appear in the supporting documentation. The authors make no claims *
	14	* about the suitability of this software for any purpose. It is *
	15	* provided "as is" without express or implied warranty. *
	16	**************************************************************************/
	17
2374af72	18	/**
	19	* @file AliHLTPHOSSanityInspector.cxx
	20	* @author Oystein Djuvsland
	21	* @date
	22	* @brief Sanity inspector for PHOS HLT
	23	*/
	24
	25	// see header file for class documentation
	26	// or
	27	// refer to README to build package
	28	// or
	29	// visit http://web.ift.uib.no/~kjeks/doc/alice-hlt
	30
3ed012c1	31	#include "AliHLTPHOSSanityInspector.h"
3ed012c1	32	#include "Rtypes.h"
	33
	34	ClassImp(AliHLTPHOSSanityInspector);
	35
	36
	37	AliHLTPHOSSanityInspector::AliHLTPHOSSanityInspector() :
9f050726	38	// AliHLTPHOSBase(),
3ed012c1	39	fMaxDifference(120)
3ed012c1	40	{
2374af72	41	//See header file for documentation
3ed012c1	42	}
	43
	44
	45	AliHLTPHOSSanityInspector::~AliHLTPHOSSanityInspector()
	46	{
2374af72	47	//See header file for documentation
3ed012c1	48	}
	49
	50
939c67e7	51
e304ea31	52	Int_t
939c67e7	53	AliHLTPHOSSanityInspector::CheckInsanity(const UInt_t* data, const Int_t N) const
3ed012c1	54	{
2374af72	55	//See header file for documentation
2374af72	56
3ed012c1	57	for(Int_t i = 1; i < N; i++)
	58	{
	59	if((((Int_t)data[i] - (Int_t)data[i-1]) > fMaxDifference) \|\| (((Int_t)data[i-1] - (Int_t)data[i]) > fMaxDifference))
	60	return 1;
	61	}
	62	return 0;
	63	}
	64
2374af72	65	/*
	66	Int_t
	67	AliHLTPHOSSanityInspector::CheckInsanity(Int_t* data, Int_t N)
	68	{
	69	//See header file for documentation
	70	for(Int_t i = 1; i < N; i++)
	71	{
	72	if((((Int_t)data[i] - (Int_t)data[i-1]) > fMaxDifference) \|\| (((Int_t)data[i-1] - (Int_t)data[i]) > fMaxDifference))
	73	return 1;
	74	}
	75	return 0;
	76	}
	77	*/
3ed012c1	78
	79
	80	Int_t
	81	AliHLTPHOSSanityInspector::CheckAndHealInsanity(UInt_t* data, Int_t N)
939c67e7	82	//
3ed012c1	83	{
2374af72	84	//See header file for documentation
2374af72	85
3ed012c1	86	Int_t crazyness = 0;
939c67e7	87
3ed012c1	88	if(N > 3)
	89	{
	90	//Require a stable start
	91	if((((Int_t)data[0] - (Int_t)data[1]) > fMaxDifference) \|\| (((Int_t)data[1] - (Int_t)data[0]) > fMaxDifference))
	92	return -1;
	93	if((((Int_t)data[1] - (Int_t)data[2]) > fMaxDifference) \|\| (((Int_t)data[2] - (Int_t)data[1]) > fMaxDifference))
	94	return -1;
939c67e7	95
939c67e7	96
3ed012c1	97	for(Int_t i = 2; i < N - 3; i++)
	98	{
	99	if((((Int_t)data[i] - (Int_t)data[i+1]) > fMaxDifference) \|\| (((Int_t)data[i+1] - (Int_t)data[i]) > fMaxDifference))
	100	{
	101	i++;
	102	if((((Int_t)data[i] -(Int_t)data[i+1]) > fMaxDifference) \|\| (((Int_t)data[i+1] - (Int_t)data[i]) > fMaxDifference))
	103	{
	104	i++;
	105	if((((Int_t)data[i] - (Int_t)data[i+1]) > fMaxDifference) \|\| (((Int_t)data[i+1] - (Int_t)data[i]) > fMaxDifference))
	106	{
	107	return -2; //Too crazy
	108	}
	109	data[i-1] = ((Int_t)data[i] + (Int_t)data[i-2])/2;
	110	crazyness++;
	111	}
	112	else
	113	return -3; //Two spikes in a row?
	114	}
	115	}
	116
	117
	118
	119	if((((Int_t)data[N - 3] -(Int_t) data[N - 2]) > fMaxDifference) \|\|
	120	(((Int_t)data[N - 2] - (Int_t)data[N - 3]) > fMaxDifference))
	121	{
	122	if((((Int_t)data[N - 2] - (Int_t)data[N - 1]) > fMaxDifference) \|\|
	123	(((Int_t)data[N - 1] - (Int_t)data[N - 2]) > fMaxDifference))
	124	{
	125	data[N - 2] = ((Int_t)data[N - 3] + (Int_t)data[N - 1])/2;
	126	return crazyness++;
	127	}
	128	return -4;
	129
	130	}
	131
	132	if((((Int_t)data[N - 2] - (Int_t)data[N - 1]) > fMaxDifference) \|\|
	133	(((Int_t)data[N - 1] - (Int_t)data[N - 2]) > fMaxDifference))
	134	{
	135	// (Int_t)data[N - 3] = (Int_t)data[N - 4] -(Int_t) data[N - 5] + (Int_t)data[N-4];
	136	data[N - 1] = data[N - 2];
	137	return crazyness++;
	138	}
	139
	140	}
	141
	142	return crazyness;
	143
	144	}
ab38011b	145
939c67e7	146
939c67e7	147
3ed012c1	148	Int_t
ab38011b	149	AliHLTPHOSSanityInspector::CheckAndHealInsanity(Int_t* data, Int_t N)
3ed012c1	150	{
2374af72	151	//See header file for documentation
2374af72	152
3ed012c1	153	Int_t crazyness = 0;
939c67e7	154
3ed012c1	155	if(N > 3)
3ed012c1	156	{
939c67e7	157	if((((Int_t)data[0] - (Int_t)data[1]) > fMaxDifference) \|\| (((Int_t)data[1] - (Int_t)data[0]) > fMaxDifference))
	158	return -1;
	159	if((((Int_t)data[1] - (Int_t)data[2]) > fMaxDifference) \|\| (((Int_t)data[2] - (Int_t)data[1]) > fMaxDifference))
	160	return -1;
	161
	162
	163	for(Int_t i = 2; i < N - 3; i++)
3ed012c1	164	{
3ed012c1	165	if((((Int_t)data[i] - (Int_t)data[i+1]) > fMaxDifference) \|\| (((Int_t)data[i+1] - (Int_t)data[i]) > fMaxDifference))
939c67e7	166	{
	167	i++;
	168	if((((Int_t)data[i] -(Int_t)data[i+1]) > fMaxDifference) \|\| (((Int_t)data[i+1] - (Int_t)data[i]) > fMaxDifference))
	169	{
	170	i++;
	171	if((((Int_t)data[i] - (Int_t)data[i+1]) > fMaxDifference) \|\| (((Int_t)data[i+1] - (Int_t)data[i]) > fMaxDifference))
	172	{
	173	return -2; //Too crazy
	174	}
	175	data[i-1] = ((Int_t)data[i] + (Int_t)data[i-2])/2;
	176	crazyness++;
	177	}
	178	else
	179	return -3; //Two spikes in a row?
	180	}
3ed012c1	181	}
3ed012c1	182
	183
	184
939c67e7	185	if((((Int_t)data[N - 3] -(Int_t) data[N - 2]) > fMaxDifference) \|\|
	186	(((Int_t)data[N - 2] - (Int_t)data[N - 3]) > fMaxDifference))
	187	{
	188	if((((Int_t)data[N - 2] - (Int_t)data[N - 1]) > fMaxDifference) \|\|
	189	(((Int_t)data[N - 1] - (Int_t)data[N - 2]) > fMaxDifference))
	190	{
	191	data[N - 2] = ((Int_t)data[N - 3] + (Int_t)data[N - 1])/2;
	192	return crazyness++;
	193	}
	194	return -4;
3ed012c1	195
939c67e7	196	}
3ed012c1	197
939c67e7	198	if((((Int_t)data[N - 2] - (Int_t)data[N - 1]) > fMaxDifference) \|\|
	199	(((Int_t)data[N - 1] - (Int_t)data[N - 2]) > fMaxDifference))
	200	{
3ed012c1	201	// (Int_t)data[N - 3] = (Int_t)data[N - 4] -(Int_t) data[N - 5] + (Int_t)data[N-4];
939c67e7	202	data[N - 1] = data[N - 2];
	203	return crazyness++;
	204	}
3ed012c1	205
939c67e7	206	}
3ed012c1	207
	208	return crazyness;
	209
	210	}