# last update on $Date: 2007/07/19 12:20:58 $ by $Author: flucke $
process Alignment = {
    
    # initialize  MessageLogger
    service = MessageLogger { 
#	untracked vstring debugModules = {"*"}
	untracked vstring destinations = { "alignment" }# {, "cout" }
	untracked vstring statistics = { "alignment" } #{, "cout" }
	untracked vstring categories = { "Alignment", "LogicError", "FwkReport", "TrackProducer"}
	#untracked PSet FwkReport     = { untracked string threshold = "WARNING" }
	#untracked PSet TrackProducer = { untracked string threshold = "WARNING" }
	
	untracked PSet cout    = { 
	    untracked string threshold = "DEBUG" # "ERROR" 
	    untracked PSet FwkReport = {
		untracked string threshold = "ERROR" 
	    }
	    untracked PSet TrackProducer = {
		untracked string threshold = "ERROR" 
	    }
	    #      	untracked bool noLineBreaks = true 
	}
	untracked PSet alignment  = { 
	    untracked string threshold = "DEBUG" 
	    untracked PSet INFO = { untracked int32 limit = 10 }
	    untracked PSet WARNING = { untracked int32 limit = 10 }
	    untracked PSet ERROR = { untracked int32 limit = -1 }
	    untracked PSet DEBUG = { untracked int32 limit = -1 }
	    untracked PSet Alignment = { untracked int32 limit = -1}
	    untracked PSet LogicError = { untracked int32 limit = -1}
	    # untracked bool noLineBreaks = true 
	}
    }
    
    # initialize magnetic field
    include "MagneticField/Engine/data/volumeBasedMagneticField.cfi"
    
    # ideal geometry and interface
    include "Geometry/CMSCommonData/data/cmsIdealGeometryXML.cfi"
    include "Geometry/TrackerNumberingBuilder/data/trackerNumberingGeometry.cfi" include "Geometry/CommonDetUnit/data/bareGlobalTrackingGeometry.cfi"
    include "Geometry/MuonNumbering/data/muonNumberingInitialization.cfi"
    
    # track selection for alignment
    include "Alignment/CommonAlignmentProducer/data/AlignmentTrackSelector.cfi"
    # replace AlignmentTrackSelector.src = ALCARECOTkAlZMuMu # ctfWithMaterialTracksBeamHaloMuon 
    replace AlignmentTrackSelector.src = ALCARECOTkAlZMuMu
    replace AlignmentTrackSelector.ptMin   = 10.
    # replace AlignmentTracks.ptMax   = 999.
    # replace AlignmentTracks.etaMin  = -5. # -2.4
    # replace AlignmentTracks.etaMax  =  5. # 2.4
    replace AlignmentTrackSelector.nHitMin = 8 # 8
    # replace AlignmentTracks.chi2nMax= 999999.
    # replace AlignmentTracks.applyNHighestPt = false
    replace AlignmentTrackSelector.nHighestPt = 2
   # replace AlignmentTrackSelector.minMultiplicity = 1 # 8

    include "Alignment/MillePedeAlignmentAlgorithm/data/PresigmaScenarios.cff"
    replace MillePedeAlignmentAlgorithm.pedeSteerer.Presigmas += TrackerORCAShortTermPresigmasDetBy10.Presigmas
    include "Alignment/ReferenceTrajectories/data/TrajectoryFactories.cff" 

    # Alignment producer
    #include "Alignment/CommonAlignmentProducer/data/AlignmentProducerCSA07t.cff"
    include "Alignment/CommonAlignmentProducer/data/AlignmentProducer.cff"

  replace AlignmentProducer.applyDbAlignment = true
  replace AlignmentProducer.ParameterBuilder.Selector = {
        vstring alignParams = {
  "PixelDets,111001,pixEnd", "BarrelDetsDS,111001", "BarrelDetsSS,101001", 
"TIDDets,111001,tecDS", "TIDDets,101001,tecSS", "TECDets,111001,tecDS", 
"TECDets,101001,tecSS"
 # "TOBHalfBarrels,111111"
 #, "PixelHalfBarrels,00rrrr"
 # "TIBHalfBarrels,111111",
 # "TECs,111111",
 # "TIDs,111111"
# "PixelHalfBarrels,00rrrr"
        }
 
   PSet pixEnd = {
            vdouble etaRanges = {-2.5,-0.0,0.0,2.5}  vdouble phiRanges = {} vdouble zRanges = {}
           vdouble rRanges   = {}
           vdouble xRanges = {}  vdouble yRanges = {}
        }

        PSet tecSS = {
            vdouble etaRanges = {-2.5,-0.1,0.1,2.5}  vdouble phiRanges = {} vdouble zRanges = {}
           vdouble rRanges   = {40., 60., 75., 999.}
           vdouble xRanges = {}  vdouble yRanges = {}
        }

        PSet tecDS = {
            vdouble etaRanges = {-2.5,-0.1,0.0,2.5}  vdouble phiRanges = {} vdouble zRanges = {}
           vdouble rRanges   = {0., 40., 60., 75.}
           vdouble xRanges = {}  vdouble yRanges = {}
        }

    }

    replace AlignmentProducer.algoConfig = {
	using MillePedeAlignmentAlgorithm
    }
    replace MillePedeAlignmentAlgorithm.mode = "mille" # "full" #"pede" #"mille"
    # replace MillePedeAlignmentAlgorithm.fileDir = "/afs/cern.ch/user/r/rmankel/scratch0/mpede-160b/jobData/jobm"
    replace MillePedeAlignmentAlgorithm.binaryFile = "milleBinaryISN.dat"
    replace MillePedeAlignmentAlgorithm.treeFile = "treeFile.root"
   replace MillePedeAlignmentAlgorithm.monitorFile = ""
    replace MillePedeAlignmentAlgorithm.beamspot = false
    replace MillePedeAlignmentAlgorithm.aliBeamspot = false
    replace MillePedeAlignmentAlgorithm.pedeSteerer.method = "spareGMRES 6 0.9"
    # replace MillePedeAlignmentAlgorithm.pedeSteerer.method = "inversion  9  0.8"
    #                                                                    <method>  n(iter)  Delta(F)
    replace MillePedeAlignmentAlgorithm.pedeSteerer.options = {
	"chisqcut  20.0  4.5","bandwidth 6"} #{,  "outlierdownweighting 3",   "dwfractioncut 0.1"}
   
    replace MillePedeAlignmentAlgorithm.pedeSteerer.pedeCommand = "/afs/cern.ch/user/m/mstoye/scratch0/MPS/pede_test"
    # replace MillePedeAlignmentAlgorithm.useTrackTsos = false
     replace MillePedeAlignmentAlgorithm.TrajectoryFactory = {
    #	using BzeroReferenceTrajectoryFactory
    	using TwoBodyDecayTrajectoryFactory
     }

#MILLEPEDEBLOCK
    
# refitting
## RM this is already included by the cff file
   include "RecoTracker/TrackProducer/data/RefitterWithMaterial.cff"
## these do not work any more, since already a replace in the cff
   replace TrackRefitter.src = "AlignmentTrackSelector"
   replace TrackRefitter.TrajectoryInEvent = true
    # needed for refit of hits:
    # include "CalibTracker/Configuration/data/SiStrip_FakeLorentzAngle.cff"
    # usually without refit: 
    replace TrackRefitter.TTRHBuilder = "WithoutRefit"# TransientTrackingRecHitBuilder: no refit of hits...
    include "RecoTracker/TransientTrackingRecHit/data/TransientTrackingRecHitBuilderWithoutRefit.cfi"
    # ... but matching for strip stereo should be redone: 
    replace ttrhbwor.Matcher = "StandardMatcher"
    # beam halo propagation needs larger phi changes going from one TEC to another
    # replace MaterialPropagator.MaxDPhi = 1000.

 source = PoolSource {
	untracked vstring fileNames = { 
"rfio:/castor/cern.ch/cms/store/mc/CSA08/Zmumu/ALCARECO/1PB_V2_RECO_TkAlZMuMu_v1/0042/22EC0C0D-4126-DD11-8EA5-001617E30D4A.root"

   }  
	untracked uint32 skipEvents = 0
    }

 #    untracked PSet maxEvents = {untracked int32 input = 100000}
 # add offline beamspot
   include "RecoVertex/BeamSpotProducer/data/BeamSpot.cff"
 #  include "RecoVertex/BeamSpotProducer/data/BeamSpotEarlyCollision.cff"
   include "Configuration/StandardSequences/data/FrontierConditions_GlobalTag.cff"
   replace GlobalTag.globaltag = "1PB_V2_RECO::All"
   path p = { AlignmentTrackSelector,offlineBeamSpot, TrackRefitter }

}