import glob, sys, os, string, io, ROOT, math, argparse, xlrd


def create_parfile(num, gasf, efie, gain_par, gain, clf, vd, wallc,clerrf):
    gainf=0
    if gain==70:
        gainf=76
    if gain==71:
        gainf=43
    if gain==72:
        gainf=44
    if gain==81:
        gainf=48
    if gain==82:
        gainf=50
    if gain==83:
        gainf=51
    if gain==84:
        gainf=52
    if gain==85:
        gainf=53
    if gain==86:
        gainf=54
    if gainf==0:
        print "No Gainfile found for gain", gain
        return 
    parfile=open("tpc/FOPI/par/tpc.run"+str(num)+".par", "w")
    parfile.write(
        "#################################################################\n"
        "# Digitisation and reconstruction parameters for TPC\n"
        "# Format: automatically ganerated\n"
        "# parameter_name:parameter_type(i-integer, d-double) parameter_value\n"
        "# TPCDigiPar\n"
        "# Description of parameters\n"
        "#################################################################)\n"
        "[TpcDigiPar]\n"
        "TpcGasFile:Int_t " + str(gasf) +"\n"
        "EField:Double_t " + str(efie) + "\n"
        "AttachFlag:Int_t 0\n"
        "LogitudinaldiffusionFlag:Int_t 1\n"
        "TransversediffusionFlag:Int_t 1\n"
        "DriftDistortionFlag:Int_t 0\n"
        "zMax:Double_t " + str(72.3) + "\n"
        "Gain:Double_t " + str(gain_par) + "\n"
        "Supression:Double_t 0.001\n"
        "Spread:Double_t 0.02\n"
        "MinSignalAmp:Double_t 6\n"  
        "rMin:Double_t 5.\n"
        "rMax:Double_t 15.5\n"
        "PadPlaneFile:Int_t 30\n"
        "PadShapeFile:Int_t 32\n"
        "PadShapeRange:Double_t 0.5\n"
        "PadShapeStep:Double_t 0.02\n"
        "PadShapeIntStep:Double_t 0.01\n"
        "ADCThreshold:Int_t 0\n"
        "ADCMax:Int_t 100000\n"
        "ADCBits:Int_t 11\n"
        "SamplingFreq:Double_t 15.55\n"
        "T0_wallclock:Double_t " + str(wallc) + "\n"
        "TimeBits:Int_t 9\n"
        "PSAThreshold:Int_t 1700\n"
        "Shaper_tDiff:Double_t 100\n"
        "Shaper_tInt:Double_t 300\n"
        "Shaper_tSig:Double_t 2\n"
        "TOTPSA:Int_t 2\n"
        "GAUSSIANNOISE:Int_t 0\n"
        "GAUSSIANNOISEAMP:Int_t 0\n"
        "Rate:Double_t 2e-2\n"
        "ManDriftVel: Double_t " + str(vd) + "\n"
        "AlignmentFile: Int_t 41\n"
        "SampleThresh: Double_t 0.\n"
        "GainCorrFile: Int_t " + str(gainf) + "\n"
        "ClusterCorrFile: Int_t " + str(clf) + "\n"
        "PadResponseFile: Int_t 80\n"
        "ClusterErrorFile: Int_t "+str(clerrf)+"\n"
        "#################################################################\n"
        "#TpcGasFile:     33: Ar/CO2 90/10 B=0.3; 34: Ne/CO2 90/10 B=0.6; 35: Ar/CO2 90/10 B=0.6;\n"
        "#T0_wallclock : 1000-20000\n")
    
    parfile.close()