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

 /**************************************************************************
 * 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.                  *
 **************************************************************************/

#include "AliHLTPHOSMIPCounter.h"
#include "AliHLTPHOSDigitContainerDataStruct.h"
#include "TH2I.h"

AliHLTPHOSMIPCounter::AliHLTPHOSMIPCounter()
    : AliHLTPHOSBase(),
    fMIPCountEvent ( 0 ),
    fMIPCountTotal ( 0 ),
    fMIPRate ( 0 ),
    fLowerBound ( 0 ),
    fUpperBound ( 0 ),
    fUpperStartTime ( 0 ),
    fLowerStartTime ( 0 ),
    fZeroThreshold ( 0 ),
    fChannelHistPtr ( 0 )
{ 

}


AliHLTPHOSMIPCounter::~AliHLTPHOSMIPCounter()
{
}

Int_t
AliHLTPHOSMIPCounter::CountMIPs(AliHLTPHOSDigitContainerDataStruct* digitContainerPtr)
{
  fMIPCountEvent = 0;
  Bool_t IsMIP = true;
  Int_t *dataPtr = 0;
  AliHLTPHOSDigitDataStruct *digitPtr;
  for(UInt_t i = 0; i < digitContainerPtr->fNDigits; i++)
  {
    digitPtr = &(digitContainerPtr->fDigitDataStruct[i]);
    dataPtr = digitPtr->fData;
    if(digitPtr->fCrazyness != 0)
    {
      continue;
    }
    if(digitPtr->fAmplitude < fLowerBound || digitPtr->fAmplitude > fUpperBound)
    {
      continue;
    }
    for(Int_t time = (Int_t)(digitPtr->fTime - 2); time < (digitPtr->fTime - 3); time++)
    {
      if((Float_t)dataPtr[time] < (digitPtr->fAmplitude - (digitPtr->fAmplitude)/10))
      {
	IsMIP = false;
	break;
      }
    }
    if(!IsMIP)
      continue;
    for(Int_t sample = 0; sample < fLowerStartTime; sample++)
      {
	if(dataPtr[sample] > fZeroThreshold || dataPtr[sample] < -fZeroThreshold)
	{
	  IsMIP = false;
	  break;
	}
      }
      if(dataPtr[(Int_t)fUpperStartTime + 3] < fZeroThreshold)
	IsMIP = false;
    if(IsMIP)
    {
      fMIPCountEvent++;
      fChannelHistPtr->Fill(digitPtr->fX, digitPtr->fZ);
    }
  }
  fMIPCountTotal += fMIPCountEvent;
  return fMIPCountEvent; 
}
Commit	Line	Data
78cc8fa3	1	/**************************************************************************
	2	* This file is property of and copyright by the ALICE HLT Project *
	3	* All rights reserved. *
	4	* *
	5	* Primary Authors: Oystein Djuvsland *
	6	* *
	7	* Permission to use, copy, modify and distribute this software and its *
	8	* documentation strictly for non-commercial purposes is hereby granted *
	9	* without fee, provided that the above copyright notice appears in all *
	10	* copies and that both the copyright notice and this permission notice *
	11	* appear in the supporting documentation. The authors make no claims *
	12	* about the suitability of this software for any purpose. It is *
	13	* provided "as is" without express or implied warranty. *
	14	**************************************************************************/
	15
	16	#include "AliHLTPHOSMIPCounter.h"
	17	#include "AliHLTPHOSDigitContainerDataStruct.h"
	18	#include "TH2I.h"
	19
	20	AliHLTPHOSMIPCounter::AliHLTPHOSMIPCounter()
	21	: AliHLTPHOSBase(),
	22	fMIPCountEvent ( 0 ),
	23	fMIPCountTotal ( 0 ),
	24	fMIPRate ( 0 ),
	25	fLowerBound ( 0 ),
	26	fUpperBound ( 0 ),
	27	fUpperStartTime ( 0 ),
	28	fLowerStartTime ( 0 ),
	29	fZeroThreshold ( 0 ),
	30	fChannelHistPtr ( 0 )
	31	{
	32
	33	}
	34
	35
	36	AliHLTPHOSMIPCounter::~AliHLTPHOSMIPCounter()
	37	{
	38	}
	39
	40	Int_t
	41	AliHLTPHOSMIPCounter::CountMIPs(AliHLTPHOSDigitContainerDataStruct* digitContainerPtr)
	42	{
	43	fMIPCountEvent = 0;
	44	Bool_t IsMIP = true;
	45	Int_t *dataPtr = 0;
	46	AliHLTPHOSDigitDataStruct *digitPtr;
1804b020	47	for(UInt_t i = 0; i < digitContainerPtr->fNDigits; i++)
78cc8fa3	48	{
	49	digitPtr = &(digitContainerPtr->fDigitDataStruct[i]);
	50	dataPtr = digitPtr->fData;
	51	if(digitPtr->fCrazyness != 0)
	52	{
	53	continue;
	54	}
	55	if(digitPtr->fAmplitude < fLowerBound \|\| digitPtr->fAmplitude > fUpperBound)
	56	{
	57	continue;
	58	}
	59	for(Int_t time = (Int_t)(digitPtr->fTime - 2); time < (digitPtr->fTime - 3); time++)
	60	{
	61	if((Float_t)dataPtr[time] < (digitPtr->fAmplitude - (digitPtr->fAmplitude)/10))
	62	{
	63	IsMIP = false;
	64	break;
	65	}
	66	}
	67	if(!IsMIP)
	68	continue;
	69	for(Int_t sample = 0; sample < fLowerStartTime; sample++)
	70	{
	71	if(dataPtr[sample] > fZeroThreshold \|\| dataPtr[sample] < -fZeroThreshold)
	72	{
	73	IsMIP = false;
	74	break;
	75	}
	76	}
	77	if(dataPtr[(Int_t)fUpperStartTime + 3] < fZeroThreshold)
	78	IsMIP = false;
	79	if(IsMIP)
	80	{
	81	fMIPCountEvent++;
	82	fChannelHistPtr->Fill(digitPtr->fX, digitPtr->fZ);
	83	}
	84	}
	85	fMIPCountTotal += fMIPCountEvent;
	86	return fMIPCountEvent;
	87	}