KVElasticScatter.h

/*
$Id: KVElasticScatter.h,v 1.6 2007/04/04 10:39:17 ebonnet Exp $
$Revision: 1.6 $
$Date: 2007/04/04 10:39:17 $
*/
 
 
#ifndef __KVELASTICSCATTER_H
#define __KVELASTICSCATTER_H
 
#include "TVector3.h"
#include "KVNameValueList.h"
#include "KVNucleus.h"
 
class TH1F;
class TList;
class KVGroup;
class KVDetector;
class KVTarget;
class KV2Body;
class TObjArray;
 
class KVElasticScatter {
    TH1F* fDepth;
    TH1F* fTheta;
    Int_t fBinE;
 
    Double_t fEnergy;
    KVNucleus fProj;
   KVNucleus fOutgoing;            //scattered nucleus
   Bool_t has_outgoing;
 
    KV2Body* fKinematics;
 
    KVDetector* fDetector;
    TList fAlignedDetectors;
    Int_t fNDets;
    KVTarget* fTarget;
    Bool_t fMultiLayer;
 
    TVector3 fBeamDirection;
    Int_t fIntLayer;
 
    Double_t fExx;
 
    TObjArray* fHistos;
    KVNameValueList fDetInd;
 
public:
 
   KVElasticScatter();
   virtual ~ KVElasticScatter();
 
   void SetRun(Int_t run);
   void SetProjectile(Int_t Z, Int_t A);
   void SetEnergy(Double_t E);
   void SetTargetScatteringLayer(const Char_t* name);
   void SetEbinning(Int_t nbins = 500);
   void SetDetector(const Char_t* det);
   void CalculateScattering(Int_t N);
   void SetTargetExcitedState(Double_t ex)
   {
      fExx = ex;
   };
   void SetOutGoing(Int_t Z, Int_t A)
   {
      fOutgoing.SetZandA(Z, A);
      has_outgoing = true;
   }
 
   TH1F* GetDepth()
   {
      return fDepth;
   };
   TH1F* GetEnergy()
   {
      return GetEnergy(fNDets - 1);
   };
   TH1F* GetTheta()
   {
      return fTheta;
   };
   TH1F* GetEnergy(const Char_t* type);
   TH1F* GetEnergy(Int_t index);
   Int_t GetNDets() const
   {
      return fNDets;
   };
   Int_t GetEbinning(void)
   {
      return fBinE;
   };
 
   ClassDef(KVElasticScatter, 1)        //Calculate elastic scattering spectra in multidetector arrays
};
 
#endif