HarpoSimChamber.cxx

Go to the documentation of this file.

#include "HarpoSimChamber.h"
#include "HarpoDebug.h"
#include "TpcSimIonisationTrack.h"

#include "TRandom3.h"
#include "TMath.h"

ClassImp(HarpoSimChamber)


HarpoSimChamber::HarpoSimChamber() : TObject()
{
  fReadout = (HarpoSimReadout*)NULL;
  fDriftChamber = (TpcSimDriftChamber*)NULL;
  fDigi = (HarpoSimElectronics*)NULL;
}

HarpoSimChamber::~HarpoSimChamber()
{
}

HarpoSimEvent* HarpoSimChamber::ProcessEvent(TpcSimMCEvent* /* mcev */)
{
  
  HarpoSimEvent *event = new HarpoSimEvent();
  
  // StartEvent();
  // Int_t nTr = mcev->GetNtracks();
  // //cout << nTr << endl;
   
  // event->SetMCEvent(mcev);
  
  // //event->SetMCEvent(ev);
  // for(Int_t iTr = 0; iTr<nTr; iTr++){
  
  //   //cout << "  Generate Track " << iTr << endl;
  //   TpcSimIonisationTrack* tr = fIonisationChamber->GenerateTrack(mcev->GetTrack(iTr));
  //   //cout << "  Add Track in Event" << endl;
  //   event->AddIonisationTrack(iTr,tr);
    
  //   fHarpoReadout->Reset();
  //   fHarpoChamber->ProcessTrack(tr);
    
  //   HarpoSimStrips* strips = fHarpoReadout->GetStrips();
  //   strips->SetTrackId(iTr);
  //   event->AddStrips(iTr,strips);
    
  //   fHarpoElectronics->ProcessStrips(strips);
  // }
  // //cout << "Finish Event" << endl;
  // FinishEvent();
  // //cout << "GetDccMapX" << endl;
  // HarpoDccMap* mapX = fHarpoElectronics->GetDccMapX();
  // event->SetDccMapX(mapX);
  // //  cout << "GetDccMapY" << endl;
  // HarpoDccMap* mapY = fHarpoElectronics->GetDccMapX();
  // event->SetDccMapY(mapY);
  
  return event;
  
}
  
  

Int_t HarpoSimChamber::ProcessTrack(TpcSimIonisationTrack* tr)
{
  
  if(!fDriftChamber){
    return -1;
  }
  if(!fReadout){
    return -1;
  }

  const Int_t npoints = tr->GetNpoints();
  Double_t xinit= 0, yinit =0 , zinit = 0;
  Double_t xfin = 0, yfin = 0, zfin = 0;
  Int_t nTot = 0;
  Double_t mean = 0;
  for(Int_t i = 0; i<npoints; i++){
    if(gHarpoDebug>2) Info("ProcessTrack","Point %d",i);
    TpcSimIonisationPoint* point = tr->GetPoint(i);
    
    Double_t x0 = point->GetX();
    Double_t y0 = point->GetY();
    Double_t z0 = point->GetZ();
    Double_t t0 = point->GetT();
    if(i==0){
      xinit = x0;
      yinit = y0;
      zinit = z0;
    }
    Double_t length = TMath::Sqrt((x0-xfin)*(x0-xfin) + (y0-yfin)*(y0-yfin) + (z0-zfin)*(z0-zfin));
    xfin = x0;
    yfin = y0;
    zfin = z0;
    Int_t nElectrons = point->GetNe();
    nTot += nElectrons; 

    mean += nElectrons*1./length;

    for(Int_t j = 0; j<nElectrons; j++){
      Double_t x, y, z, t;
      Int_t test = fDriftChamber->DriftElectron(x0,y0,z0,t0,x,y,z,t);
      if(test<1) continue;
     
      hSimReadoutRegion region = hSimReadoutRegionDrift; 
      if(z0-z<4) region = hSimReadoutRegionTransfer2;
      if(z0-z<2) region = hSimReadoutRegionTransfer1;
      fReadout->ProcessElectron(x,y,t,region);
    }
  }
  
  Double_t length2 = TMath::Sqrt((xinit-xfin)*(xinit-xfin) + (yinit-yfin)*(yinit-yfin) + (zinit-zfin)*(zinit-zfin));
  // if(length)
  //   Info("ProcessTrack","%d electrons / %.3g mm = %.3g e/cm",nTot,length,10*nTot/length);
  if(npoints && length2){
    mean /= npoints;
    if(gHarpoDebug>1) 
      Info("ProcessTrack","%d electrons / %.3g mm = %.3g e/cm           %.3g e/cm",nTot,length2,10*nTot/length2, 10*mean);
  }

  return 1;
  
}


Int_t HarpoSimChamber::StartEvent()
{

  if(gHarpoDebug>1) Info("StartEvent","Reset readout");

  if(fReadout)
    fReadout->Reset();
  if(gHarpoDebug>1) Info("StartEvent","Reset electronics");
  if(fDigi)
    fDigi->Reset();

  return 1;
}


Int_t HarpoSimChamber::FinishEvent()
{
  fDigi->Noise();
  fDigi->ZeroSuppression();
  // fDigi->AGETmultiplicity();
  return 1;
}