MDT_ResponseTest.cxx File Reference

#include "MDT_Response/MDT_Response.h"
#include "CLHEP/Random/RandomEngine.h"
#include "CLHEP/Random/RanluxEngine.h"
#include "CLHEP/Random/RandPoisson.h"
#include "CLHEP/Random/RandFlat.h"
#include "TH1.h"
#include "TFile.h"
#include <cerrno>
#include <climits>
#include <cstdlib>
#include <iostream>
#include "CxxUtils/checker_macros.h"

Go to the source code of this file.

Functions
CLHEP::RandFlat	ranFlat (ranEngine)
double	generateFlat ()
double	generateRadius ()
double	generatePositionAlongTube ()
void	generateEvents (int eventMax)
int	main (int argc, char *argv[])

Variables
	ATLAS_NO_CHECK_FILE_THREAD_SAFETY
CLHEP::RanluxEngine	ranEngine
MDT_Response	response
double	tubeMax = 14.7
double	tubeLength = 5000.

Function Documentation

generateEvents()

void generateEvents

(

int

eventMax

)

Definition at line 46 of file MDT_ResponseTest.cxx.

                                    {
  TH1* adcHist = new TH1F("adc","adc",400,0,200);
  TH1* driftHist = new TH1F("driftTime","driftTime",1000,0,1000);
  TH1* radiusHit = new TH1F("radiusHit","radiusHit",100,0,tubeMax);
  TH1* radiusNoHit = new TH1F("radiusNoHit","radiusNoHit",100,0,tubeMax);
  TH1* posLengthHit = new TH1F("posLengthHit","posLengthHit",100,0,tubeLength);
  TH1* posLengthNoHit = new TH1F("posLengthNoHit","posLengthNoHit",100,0,tubeLength);
  for( int i=0;i<eventMax; ++i ){
    
    // generate drift radius + position along tube
    double radius = generateRadius();
    double positionAlongTube = generatePositionAlongTube();
    response.SetSegment(radius,positionAlongTube);
 
    // check whether tube fired as we also simulate tube inefficiency
    bool hasHit = response.GetSignal(&ranEngine);
    if( hasHit ){
      // get the drift time (TDC count) and ADC
      double driftTime = response.DriftTime();
      double adc       = response.AdcResponse();
      adcHist->Fill(adc);
      driftHist->Fill(driftTime);
      radiusHit->Fill(radius);
      posLengthHit->Fill(positionAlongTube);
    }else{
      radiusNoHit->Fill(radius);
      posLengthNoHit->Fill(positionAlongTube);
    }
    
  }
 
}

generateFlat()

double generateFlat

(

)

Definition at line 34 of file MDT_ResponseTest.cxx.

                      {
  return ranFlat.fire();
}

generatePositionAlongTube()

double generatePositionAlongTube

(

)

Definition at line 42 of file MDT_ResponseTest.cxx.

                                   {
  return tubeLength*generateFlat();
}

generateRadius()

double generateRadius

(

)

Definition at line 38 of file MDT_ResponseTest.cxx.

                        {
  return tubeMax*generateFlat();
}

main()

int main	(	int	argc,
		char *	argv[] )

Definition at line 79 of file MDT_ResponseTest.cxx.

                                 {
 
  int eventMax = 20000;
  if( argc > 1 ) {
    char *endptr;
    errno = 0;
    long convArg = strtol(argv[1], &endptr, 0);
    if(errno == ERANGE || *endptr != '\0' || argv[1] == endptr) {
      std::cout<<"Invalid parameter! Quit now!"<<std::endl;
      exit(1);
    }
    // Only needed if sizeof(int) < sizeof(long)
    if(convArg < INT_MIN || convArg > INT_MAX) {
      std::cout<<"Invalid parameter! Quit now!"<<std::endl;
      exit(1);
    }
    int tempMax = (int) convArg;
    if (tempMax > 0) eventMax = tempMax; 
  } 
  
  std::cout << " Starting simulation of MDT_response, events " << eventMax << std::endl;
 
  TFile* outputFile = new TFile("MDT_ResponseTest.root","RECREATE");
 
  // simulate long chamber
  tubeLength = 5000.;
  TDirectory* dir = outputFile->mkdir("LongTubes");
  dir->cd();
  generateEvents(eventMax);
 
  // simulate short chamber
  tubeLength = 500.;
  dir = outputFile->mkdir("ShortTubes");
  dir->cd();
  generateEvents(eventMax);
 
  // simulate short chamber, large incident angle
  double increasedPathLength = 0.25; // here one could use the angle
  response.SetClusterDensity(8.5*increasedPathLength);
  dir = outputFile->mkdir("ShortTubesMediumAngle");
  dir->cd();
  generateEvents(eventMax);
 
//   // simulate short chamber, large incident angle
//   increasedPathLength = 4; // here one could use the angle
//   response.SetClusterDensity(8.5*increasedPathLength);
//   dir = outputFile->mkdir("ShortTubesLargeAngle");
//   dir->cd();
//   generateEvents(eventMax);
 
  outputFile->Write();
  outputFile->Close();
 
  return 0;
}

ranFlat()

CLHEP::RandFlat ranFlat

(

ranEngine

)

Variable Documentation

ATLAS_NO_CHECK_FILE_THREAD_SAFETY

ATLAS_NO_CHECK_FILE_THREAD_SAFETY

Definition at line 21 of file MDT_ResponseTest.cxx.

ranEngine

CLHEP::RanluxEngine ranEngine

Definition at line 24 of file MDT_ResponseTest.cxx.

response

MDT_Response response

Definition at line 28 of file MDT_ResponseTest.cxx.

tubeLength

double tubeLength = 5000.

Definition at line 32 of file MDT_ResponseTest.cxx.

tubeMax

double tubeMax = 14.7

Definition at line 31 of file MDT_ResponseTest.cxx.